Source author record

J. Leenaarts

J. Leenaarts appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2025arXiv

The making of robust and highly performing imaging spectropolarimeters for large solar telescopes

We discuss the requirements, concepts, simulations, implementation, and calibration of two dual Fabry-Perot based imaging spectropolarimeters, CRISP and CHROMIS, at the Swedish 1-meter Solar Telescope, and CRISP2 that is under construction. These instruments use a combination of a high-resolution and a low-resolution etalon together with an order-sorting prefilter to define the bandpass. The overall design is made robust and stable by tailoring the low-resolution etalon reflectivity to accommodate expected cavity errors from both etalons, and by using a compact optical design that eliminates the need for folding mirrors. By using a telecentric design based on lenses rather than mirrors, image degradation by the FPI system is negligible, as shown in a previous publication, and the throughput of the system is maximised. Initial alignment, and maintaining that alignment over time, is greatly simplified. The telecentric design allows full calibration and/or modelling of essential system parameters to be carried out without interfering with the optical setup. We also discuss briefly the polarimeters developed for CRISP and CHROMIS. The high performance of CRISP and CHROMIS has been demonstrated in an earlier publication through measurements of the granulation contrast and comparisons with similar measurements simultaneously made through broadband continuum filters. Here, we focus on the aspects of the design that are central to enabling high performance and robustness, but also discuss the calibration and processing of the data, and use a few examples of processed data to demonstrate the achievable image and data quality. We put forward a proposal for a similar conceptual design for the European Solar Telescope and conclude by discussing potential problems of the proposed approach to designs of this type. Some aspects of these FPI systems may be of interest also outside the solar community.

preprint2022arXiv

Chromospheric Extension of the MURaM Code

Detailed numerical models of chromosphere and corona are required to understand the heating of the solar atmosphere. An accurate treatment of the solar chromosphere is complicated by the effects arising from Non Local Thermodynamic Equilibrium (NLTE) radiative transfer. A small number of strong, highly scattering lines dominate the cooling and heating in the chromosphere. Additionally, the recombination times of ionised hydrogen are longer than the dynamical timescales, requiring a non-equilibrium (NE) treatment of hydrogen ionisation. The MURaM code is extended to include the physical process required for accurate simulation of the solar chromosphere, as implemented in the Bifrost code. This includes a time-dependent treatment of hydrogen ionisation, a scattering multi-group radiation transfer scheme and approximations for NLTE radiative cooling. The inclusion of NE and NLTE physics has a large impact on the structure of the chromosphere; the NE treatment of hydrogen ionisation leads to a higher ionisation fraction and enhanced populations in the first excited state throughout cold inter-shock regions of the chromosphere. Additionally this prevents hydrogen ioniation from buffering energy fluctuations, leading to hotter shocks and cooler inter-shock regions. The hydrogen populations in the ground and first excited state are enhanced by $10^2-10^3$ in the upper chromosphere and up to $10^9$ near the transition region. Including the necessary NLTE physics leads to significant differences in chromospheric structure and dynamics. The thermodynamics and hydrogen populations calculated using the extended version of the MURaM code are consistent with previous non-equilibrium simulations. The electron number and temperature calculated using the non-equilibrium treatment of the chromosphere are required to accurately synthesise chromospheric spectral lines.

preprint2022arXiv

Heating of the solar chromosphere through current dissipation

The solar chromosphere is heated to temperatures higher than predicted by radiative equilibrium. This excess heating is greater in active regions where the magnetic field is stronger. We aim to investigate the magnetic topology associated with an area of enhanced millimeter (mm) brightness temperatures in a solar active region mapped by the Atacama Large Millimeter/submillimeter Array (ALMA) using spectropolarimetric co-observations with the 1-m Swedish Solar Telescope (SST). We used Milne-Eddington inversions, nonlocal thermodynamic equilibrium (non-LTE) inversions, and a magnetohydrostatic extrapolation to obtain constraints on the three-dimensional stratification of temperature, magnetic field, and radiative energy losses. We compared the observations to a snapshot of a magnetohydrodynamics simulation and investigate the formation of the thermal continuum at 3 mm using contribution functions. We find enhanced heating rates in the upper chromosphere of up to $\sim 5\rm\,kW\,m^{-2}$, where small-scale emerging loops interact with the overlying magnetic canopy leading to current sheets as shown by the magnetic field extrapolation. Our estimates are about a factor of two higher than canonical values, but they are limited by the ALMA spatial resolution ($\sim 1.2^{\prime\prime}$). Band 3 brightness temperatures reach about $\sim10^{4}\,$K in the region, and the transverse magnetic field strength inferred from the non-LTE inversions is on the order of $\sim 500\,$G in the chromosphere. We are able to quantitatively reproduce many of the observed features, including the integrated radiative losses in our numerical simulation. We conclude that the heating is caused by dissipation in current sheets. However, the simulation shows a complex stratification in the flux emergence region where distinct layers may contribute significantly to the emission in the mm continuum.

preprint2022arXiv

Quiet Sun Center to Limb Variation of the Linear Polarization Observed by CLASP2 Across the Mg II h & k Lines

The CLASP2 (Chromospheric LAyer SpectroPolarimeter 2) sounding rocket mission was launched on 2019 April 11. CLASP2 measured the four Stokes parameters of the Mg II h & k spectral region around 2800 Angstroms along a 200 arcsecond slit at three locations on the solar disk, achieving the first spatially and spectrally resolved observations of the solar polarization in this near ultraviolet region. The focus of the work presented here is the center-to-limb variation of the linear polarization across these resonance lines, which is produced by the scattering of anisotropic radiation in the solar atmosphere. The linear polarization signals of the Mg II h & k lines are sensitive to the magnetic field from the low to the upper chromosphere through the Hanle and magneto-optical effects. We compare the observations to theoretical predictions from radiative transfer calculations in unmagnetized semi-empirical models, arguing that magnetic fields and horizontal inhomogeneities are needed to explain the observed polarization signals and spatial variations. This comparison is an important step in both validating and refining our understanding of the physical origin of these polarization signatures, and also in paving the way toward future space telescopes for probing the magnetic fields of the solar upper atmosphere via ultraviolet spectropolarimetry.

preprint2022arXiv

Spatio-temporal analysis of chromospheric heating in a plage region

Our knowledge of the heating mechanisms that are at work in the chromosphere of plage regions remains highly unconstrained from observational studies. The purpose of our study is to estimate the chromospheric heating terms from a plage dataset, characterize their spatio-temporal distribution and set constraints on the heating processes that are at work. We make use of NLTE inversions to infer a model of the photosphere and chromosphere of a plage dataset acquired with the Swedish 1-m Solar Telescope. We use this model atmosphere to calculate the chromospheric radiative losses from H i, Ca ii and Mg ii atoms. We approximate the chromospheric heating terms by the net radiative losses predicted by the inverted model. In order to make the analysis of time-series over a large field-of-view computationally tractable, we make use of a neural network. In the lower chromosphere, the contribution from the Ca ii lines is dominant and located in the surroundings of the photospheric footpoints. In the upper chromosphere, the H i contribution is dominant. Radiative losses in the upper chromosphere form an homogeneous patch that covers the plage region. The net radiative losses can be split in a periodic component with an average amplitude of ampQ = 7.6 kW m^{-2} and a static (or very slowly evolving) component with a mean value of -26.1 kW m^{-2}. Our interpretation is that in the lower chromosphere, the radiative losses are tracing the sharp lower edge of the hot magnetic canopy, where the electric current is expected to be large. In the upper chromosphere, both the magnetic field and the distribution of net radiative losses are room-filling, whereas the amplitude of the periodic component is largest. Our results suggest that acoustic wave heating may be responsible for one third of the energy deposition in the upper chromosphere, whereas other heating mechanisms must be responsible for the rest.

preprint2022arXiv

Subarcsecond imaging of a solar active region filament with ALMA and IRIS

Quiescent filaments appear as absorption features on the solar disk when observed in chromospheric lines and at continuum wavelengths in the millimeter (mm) range. Active region (AR) filaments are their small-scale, low-altitude analogues, but they could not be resolved in previous mm observations. This spectral diagnostic can provide insight into the details of the formation and physical properties of their fine threads, which are still not fully understood. Here, we shed light on the thermal structure of an AR filament using high-resolution brightness temperature ($T_{\rm b}$) maps taken with ALMA Band 6 complemented by simultaneous IRIS near-UV spectra, Hinode/SOT photospheric magnetograms, and SDO/AIA extreme-UV images. Some of the dark threads visible in the AIA 304 Å passband and in the core of Mg II resonance lines have dark ($T_{\rm b}<5000$K) counterparts in the 1.25 mm maps, but their visibility significantly varies across the filament spine and in time. These opacity changes are possibly related to variations in temperature and electron density in filament fine structures. The coolest $T_{\rm b}$ values ($<$5000 K) coincide with regions of low integrated intensity in the Mg II h and k lines. ALMA Band 3 maps taken after the Band 6 ones do not clearly show the filament structure, contrary to the expectation that the contrast should increase at longer wavelengths based on previous observations of quiescent filaments. The ALMA maps are not consistent with isothermal conditions, but the temporal evolution of the filament may partly account for this.

preprint2021arXiv

Active region chromospheric magnetic fields

Context. A proper estimate of the chromospheric magnetic fields is believed to improve modelling of both active region and coronal mass ejection evolution. Aims. We investigate the similarity between the chromospheric magnetic field inferred from observations and the field obtained from a magnetohydrostatic (MHS) extrapolation. Methods. Based Fe i 6173 Å and Ca ii 8542 Å observations of NOAA active region 12723, we employed the spatially-regularised weak-field approximation (WFA) to derive the vector magnetic field in the chromosphere from Ca ii, as well as non-LTE inversions of Fe i and Ca ii to infer a model atmosphere for selected regions. Milne-Eddington inversions of Fe i serve as photospheric boundary for the MHS model that delivers the three-dimensional field, gas pressure and density. Results. For the line-of-sight component, the MHS chromospheric field generally agrees with the non-LTE inversions and WFA, but tends to be weaker than those when larger in magnitude than 300 G. The observationally inferred transverse component is stronger, especially in magnetically weaker regions, yet the qualitative distribution with height is similar to the MHS results. For either field component the MHS chromospheric field lacks the fine structure derived from the inversions. Furthermore, the MHS model does not recover the magnetic imprint from a set of high fibrils connecting the main polarities. Conclusions. The MHS extrapolation and WFA provide a qualitatively similar chromospheric field, where the azimuth of the former is better aligned with Ca ii 8542 Å fibrils than that of the WFA, especially outside strong-field concentrations. The amount of structure as well as the transverse field strengths are underestimated by the MHS extrapolation. This underscores the importance of considering a chromospheric magnetic field constraint in data-driven modelling of active regions.

preprint2021arXiv

The formation and heating of chromospheric fibrils in a radiation-MHD simulation

Aims: We examine the movements of mass elements within dense fibrils using passive tracer particles (corks) in order to understand fibril creation and destruction processes. Methods: Simulated fibrils were selected at times when they were visible in an H$α$ image proxy. The corks were selected within fibril H$α$ formation regions. From this set, a cork was selected, and the field line passing through it was constructed. Other fibrilar corks close to this fieldline were also selected. Pathlines were constructed, revealing the locations of the mass elements forward and backward in time. The forces acting on these mass elements were analysed. Results: The main process of fibrilar loading in the simulation is different to the mass loading scenario in which waves steepen into shocks and push material upwards along the fieldlines from near their footpoints. Twisted low lying fieldlines were destabilised and then they untwisted, lifting the material trapped above their apexes via the Lorentz force. Subsequently, the majority of the mass drained down the fieldlines towards one or both footpoints under gravity. Material with large horizontal velocities could also elevated in rising fieldlines, creating somewhat parabolic motions, but material was not generally moving upward along a stationary magnetic fieldline during loading. Conclusions: The processes observed in the simulation are plausible additional scenarios. Criteria for observing such events are described. It is desirable that our simulations can also form more densely-packed fibrils from material fed from the base of field footpoints. Experimental parameters required to achieve this are discussed.

preprint2020arXiv

Non-LTE inversions of a confined X2.2 flare: I. Vector magnetic field in the photosphere and chromosphere

Obtaining the magnetic field vector accurately in the solar atmosphere is essential for studying changes in field topology during flares and to reliably model space weather. We tackle this problem by applying various inversion methods to a confined X2.2 flare in NOAA AR 12673 on September 6, 2017, comparing the photospheric and chromospheric magnetic field vector with those from two numerical models of this event. We obtain the photospheric field from Milne-Eddington (ME) and (non-)local thermal equilibrium (non-LTE) inversions of Hinode SOT/SP Fe I 6301.5Å and 6302.5Å. The chromospheric field is obtained from a spatially-regularised weak field approximation (WFA) and non-LTE inversions of Ca II 8542Å observed with CRISP at the Swedish 1-m Solar Telescope. The LTE- and non-LTE-inferred photospheric field components are strongly correlated throughout the atmosphere, with stronger field and higher temperatures in the non-LTE inversions. For the chromospheric field, the non-LTE inversions correlate well with the spatially-regularised WFA. We find strong-field patches of over 4.5 kG in the photosphere, co-located with similar concentrations exceeding 3 kG in the chromosphere. The obtained field strengths are up to 2-3 times higher than in the numerical models, with more concentrated and structured photosphere-to-chromosphere shear close to the polarity inversion line. The LTE and non-LTE Fe I inversions yield essentially the same photospheric field, while ME inversions fail to reproduce the field vector orientation where Fe I is in emission. Our inversions confirm the locations of flux rope footpoints that are predicted by numerical models. However, pre-processing and lower spatial resolution lead to weaker and smoother field in the models than what the data indicate. This emphasises the need for higher spatial resolution in the models to better constrain pre-eruptive flux ropes.

preprint2020arXiv

Observational constraints on the origin of the elements II. 3D non-LTE formation of Ba ii lines in the solar atmosphere

Context. The pursuit of more realistic spectroscopic modelling and consistent abundances has led us to begin a new series of papers designed to improve current solar and stellar abundances of various atomic species. To achieve this, we have began updating the three-dimensional (3D) non-local thermodynamic equilibrium (non-LTE) radiative transfer code, Multi3D, and the equivalent one-dimensional (1D) non-LTE radiative transfer code, MULTI. Aims. We examine our improvements to these codes by redetermining the solar barium abundance. Barium was chosen for this test as it is an important diagnostic element of the s-process in the context of galactic chemical evolution. New Ba II + H collisional data for excitation and charge exchange reactions computed from first principles had recently become available and were included in the model atom. The atom also includes the effects of isotopic line shifts and hyperfine splitting. Method. A grid of 1D LTE barium lines were constructed with MULTI and fit to the four Ba II lines available to us in the optical region of the solar spectrum. Abundance corrections were then determined in 1D non-LTE, 3D LTE, and 3D non-LTE. A new 3D non-LTE solar barium abundance was computed from these corrections. Results. We present for the first time the full 3D non-LTE barium abundance of $A({\rm Ba})=2.27\pm0.02\pm0.01$, which was derived from four individual fully consistent barium lines. Errors here represent the systematic and random errors, respectively.

preprint2019arXiv

The multi-thermal chromosphere: inversions of ALMA and IRIS data

Numerical simulations of the solar chromosphere predict a diverse thermal structure with both hot and cool regions. Observations of plage regions in particular feature broader and brighter chromospheric lines, which suggest that they are formed in hotter and denser conditions than in the quiet-Sun, but also implies a non-thermal component whose source is unclear. We revisit the problem of the stratification of temperature and microturbulence in plage now adding millimeter continuum observations provided by ALMA to inversions of near-ultraviolet IRIS spectra as a powerful new diagnostic to disentangle the two parameters. We fit cool chromospheric holes and track the fast evolution of compact mm brightenings in the plage region. We use the STiC non-LTE inversion code to simultaneously fit real ultraviolet and millimeter spectra in order to infer the thermodynamic parameters of the plasma. We confirm the anticipated constraining potential of ALMA in non-LTE inversions of the solar chromosphere. We find significant differences between the inversion results of IRIS data alone compared to the results of a combination with the mm data: the IRIS+ALMA inversions have increased contrast and temperature range, and tend to favor lower values of microturbulence in the chromosphere of plage. The average brightness temperature of the plage region at 1.25 mm is 8500 K, but the ALMA maps also show much cooler ($\sim3000$ K) and hotter ($\sim11\,000$ K) evolving features partially seen in other diagnostics. To explain the former, the inversions require the existence of localized, low temperature regions in the chromosphere where molecules such as CO could form. The hot features could sustain such high temperatures due to non-equilibrium hydrogen ionization effects in a shocked chromosphere - a scenario that is supported by low-frequency shock wave patterns found in the MgII lines probed by IRIS.

preprint2015arXiv

Non-LTE oxygen line formation in 3D hydrodynamic model stellar atmospheres

The OI 777 nm lines are among the most commonly used diagnostics for the oxygen abundances in the atmospheres of FGK-type stars. However, they form in conditions that are far from local thermodynamic equilibrium (LTE). We explore the departures from LTE of atomic oxygen, and their impact on OI lines, across the Stagger-grid of three-dimensional hydrodynamic model atmospheres. For the OI 777 nm triplet we find significant departures from LTE. These departures are larger in stars with larger effective temperatures, smaller surface gravities, and larger oxygen abundances. We present grids of predicted 3D non-LTE based equivalent widths for the OI616nm, [OI] 630 nm, [OI] 636 nm, and OI 777 nm lines, as well as abundance corrections to 1D LTE based results.

preprint2015arXiv

The Galactic chemical evolution of oxygen inferred from 3D non-LTE spectral line formation calculations

We revisit the Galactic chemical evolution of oxygen, addressing the systematic errors inherent in classical determinations of the oxygen abundance that arise from the use of one dimensional hydrostatic (1D) model atmospheres and from the assumption of local thermodynamic equilibrium (LTE). We perform detailed 3D non-LTE radiative transfer calculations for atomic oxygen lines across a grid of 3D hydrodynamic stag- ger model atmospheres for dwarfs and subgiants. We apply our grid of predicted line strengths of the [OI] 630 nm and OI 777 nm lines using accurate stellar parameters from the literature. We infer a steep decay in [O/Fe] for [Fe/H] $\gtrsim$ -1.0, a plateau [O/Fe] $\approx$ 0.5 down to [Fe/H] $\approx$ -2.5 and an increasing trend for [Fe/H] $\lesssim$ -2.5. Our 3D non-LTE calculations yield overall concordant results from the two oxygen abundance diagnostics.

preprint2014arXiv

The Interface Region Imaging Spectrograph (IRIS)

The Interface Region Imaging Spectrograph (IRIS) small explorer spacecraft provides simultaneous spectra and images of the photosphere, chromosphere, transition region, and corona with 0.33-0.4 arcsec spatial resolution, 2 s temporal resolution and 1 km/s velocity resolution over a field-of-view of up to 175 arcsec x 175 arcsec. IRIS was launched into a Sun-synchronous orbit on 27 June 2013 using a Pegasus-XL rocket and consists of a 19-cm UV telescope that feeds a slit-based dual-bandpass imaging spectrograph. IRIS obtains spectra in passbands from 1332-1358, 1389-1407 and 2783-2834 Angstrom including bright spectral lines formed in the chromosphere (Mg II h 2803 Angstrom and Mg II k 2796 Angstrom) and transition region (C II 1334/1335 Angstrom and Si IV 1394/1403 Angstrom). Slit-jaw images in four different passbands (C II 1330, Si IV 1400, Mg II k 2796 and Mg II wing 2830 Angstrom) can be taken simultaneously with spectral rasters that sample regions up to 130 arcsec x 175 arcsec at a variety of spatial samplings (from 0.33 arcsec and up). IRIS is sensitive to emission from plasma at temperatures between 5000 K and 10 MK and will advance our understanding of the flow of mass and energy through an interface region, formed by the chromosphere and transition region, between the photosphere and corona. This highly structured and dynamic region not only acts as the conduit of all mass and energy feeding into the corona and solar wind, it also requires an order of magnitude more energy to heat than the corona and solar wind combined. The IRIS investigation includes a strong numerical modeling component based on advanced radiative-MHD codes to facilitate interpretation of observations of this complex region. Approximately eight Gbytes of data (after compression) are acquired by IRIS each day and made available for unrestricted use within a few days of the observation.

preprint2013arXiv

The formation of IRIS diagnostics I. A quintessential model atom of Mg II and general formation properties of the Mg II h&k lines

NASA's Interface Region Imaging Spectrograph (IRIS) space mission will study how the solar atmosphere is energized. IRIS contains an imaging spectrograph that covers the Mg II h&k lines as well as a slit-jaw imager centered at Mg II k. Understanding the observations will require forward modeling of Mg II h&k line formation from 3D radiation-MHD models. This paper is the first in a series where we undertake this forward modeling. We discuss the atomic physics pertinent to h&k line formation, present a quintessential model atom that can be used in radiative transfer computations and discuss the effect of partial redistribution (PRD) and 3D radiative transfer on the emergent line profiles. We conclude that Mg II h&k can be modeled accurately with a 4-level plus continuum Mg II model atom. Ideally radiative transfer computations should be done in 3D including PRD effects. In practice this is currently not possible. A reasonable compromise is to use 1D PRD computations to model the line profile up to and including the central emission peaks, and use 3D transfer assuming complete redistribution to model the central depression.

preprint2013arXiv

The formation of IRIS diagnostics II. The formation of the Mg II h&k lines in the solar atmosphere

NASA's Interface Region Imaging Spectrograph (IRIS) small explorer mission will study how the solar atmosphere is energized. IRIS contains an imaging spectrograph that covers the Mg II h&k lines as well as a slit-jaw imager centered at Mg II k. Understanding the observations requires forward modeling of Mg II h&k line formation from 3D radiation-MHD models. We compute the vertically emergent h&k intensity from a snapshot of a dynamic 3D radiation-MHD model of the solar atmosphere, and investigate which diagnostic information about the atmosphere is contained in the synthetic line profiles. We find that the Doppler shift of the central line depression correlates strongly with the vertical velocity at optical depth unity, which is typically located less than 200 km below the transition region (TR). By combining the Doppler shifts of the h and the k line we can retrieve the sign of the velocity gradient just below the TR. The intensity in the central line depression is anticorrelated with the formation height, especially in subfields of a few square Mm. This intensity could thus be used to measure the spatial variation of the height of the transition region. The intensity in the line-core emission peaks correlates with the temperature at its formation height, especially for strong emission peaks. The peaks can thus be exploited as a temperature diagnostic. The wavelength difference between the blue and red peaks provides a diagnostic of the velocity gradients in the upper chromosphere. The intensity ratio of the blue and red peaks correlates strongly with the average velocity in the upper chromosphere. We conclude that the Mg II h&k lines are excellent probes of the very upper chromosphere just below the transition region, a height regime that is impossible to probe with other spectral lines.

preprint2012arXiv

Chromospheric magnetic fields. Observations, simulations and their interpretation

The magnetic field of the quiet-Sun chromosphere remains a mystery for solar physicists. The reduced number of chromospheric lines are intrinsically hard to model and only a few of them are magnetically sensitive. In this work, we use a 3D numerical simulation of the outer layers of the solar atmosphere, to asses the reliability of non-LTE inversions, in this case applied to the Ca II 8542 Å line. We show that NLTE inversions provide realistic estimates of physical quantities from synthetic observations.

preprint2011arXiv

Quiet-Sun imaging asymmetries in NaI D1 compared with other strong Fraunhofer lines

Imaging spectroscopy of the solar atmosphere using the NaI D1 line yields marked asymmetry between the blue and red line wings: sampling a quiet-Sun area in the blue wing displays reversed granulation, whereas sampling in the red wing displays normal granulation. The MgI b2 line of comparable strength does not show this asymmetry, nor does the stronger CaII 8542 line. We demonstrate the phenomenon with near-simultaneous spectral images in NaI D1, MgI b2, and CaII 8542 from the Swedish 1-m Solar Telescope. We then explain it with line-formation insights from classical 1D modeling and with a 3D magnetohydrodynamical simulation combined with NLTE spectral line synthesis that permits detailed comparison with the observations in a common format. The cause of the imaging asymmetry is the combination of correlations between intensity and Dopplershift modulation in granular overshoot and the sensitivity to these of the steep profile flanks of the NaI D1 line. The MgI b2 line has similar core formation but much wider wings due to larger opacity buildup and damping in the photosphere. Both lines obtain marked core asymmetry from photospheric shocks in or near strong magnetic concentrations, less from higher-up internetwork shocks that produce similar asymmetry in the spatially averaged CaII 8542 profile.

preprint2010arXiv

Numerical simulations of the quiet chromosphere

Numerical simulations of the solar chromosphere have become increasingly realistic over the past 5 years. However, many observed chromospheric structures and behavior are not reproduced. Current models do not show fibrils in Ca II 8542, and neither reproduce the Ca II 8542 bisector. The emergent H-alpha line core intensity computed from the models show granulation instead of chromospheric shocks or fibrils. I discuss these deficiencies and speculate about what physics should be included to alleviate these shortcomings.

preprint2009arXiv

On-disk counterparts of type II spicules in the Ca II 854.2 nm and Halpha lines

Recently a second type of spicules was discovered at the solar limb with the Solar Optical Telescope (SOT) onboard the Japanese Hinode spacecraft. These previously unrecognized type II spicules are thin chromospheric jets that are shorter-lived (10-60 s) and that show much higher apparent upward velocities (of order 50-100 km/s) than the classical spicules. Since they have been implicated in providing hot plasma to coronal loops, their formation, evolution and properties are important ingredients for a better understanding of the mass and energy balance of the low solar atmosphere. Here we report on the discovery of the disk counterparts of type II spicules using spectral imaging data in the Ca II 854.2 nm and Halpha lines with the CRisp Imaging SpectroPolarimeter (CRISP) at the Swedish Solar Telescope (SST) in La Palma. We find rapid blueward excursions in the line profiles of both chromospheric lines that correspond to thin, jet-like features that show apparent velocities of order 50 km/s. These blueward excursions seem to form a separate absorbing component with Doppler shifts of order 20 and 50 km/s for the Ca II 854.2 nm and Halpha line respectively. We show that the appearance, lifetimes, longitudinal and transverse velocities and occurrence rate of these rapid blue excursions on the disk are very similar to those of the type II spicules at the limb. A detailed study of the spectral line profiles in these events suggests that plasma is accelerated along the jet, and plasma is being heated throughout the short lifetime of the event.

preprint2009arXiv

The quiet solar atmosphere observed and simulated in NaD1

The NaD1 line in the solar spectrum is sometimes attributed to the solar chromosphere. We study its formation in quiet-Sun network and internetwork. We first present high-resolution profile-resolved images taken in this line with the imaging spectrometer IBIS at the Dunn Solar Telescope and compare these to simultaneous chromospheric images taken in Ca 8542 and Halpha. We then model NaD1 formation by performing 3D NLTE profile synthesis for a snapshot from a 3D radiation-magnetohydrodynamics simulation. We find that most NaD1 brightness is not chromospheric but samples the magnetic concentrations that make up the quiet-Sun network in the photosphere, well below the height where they merge into chromospheric canopies, with aureoles from 3D resonance scattering. The line core is sensitive to magneto-acoustic shocks in and near magnetic concentrations, where shocks occur deeper than elsewhere, and may provide evidence of heating sited deep within magnetic concentrations.

preprint2009arXiv

Three-dimensional non-LTE radiative transfer computation of the Ca 8542 infrared line from a radiation-MHD simulation

Interpretation of imagery of the solar chromosphere in the widely used \CaIIIR infrared line is hampered by its complex, three-dimensional and non-LTE formation. Forward modelling is required to aid understanding. We use a 3D non-LTE radiative transfer code to compute synthetic \CaIIIR images from a radiation-MHD simulation of the solar atmosphere spanning from the convection zone to the corona. We compare the simulation with observations obtained with the CRISP filter at the Swedish 1--m Solar Telescope. We find that the simulation reproduces dark patches in the blue line wing caused by Doppler shifts, brightenings in the line core caused by upward-propagating shocks and thin dark elongated structures in the line core that form the interface between upward and downward gas motion in the chromosphere. The synthetic line core is narrower than the observed one, indicating that the sun exhibits both more vigorous large-scale dynamics as well as small scale motions that are not resolved within the simulation, presumably owing to a lack of spatial resolution.