Search

The heating of the chromosphere in internetwork regions remains one of the foremost open questions in solar physics. In the present study, we tackle this old problem by using a very-high-spatial-resolution simulation of quiet-Sun conditions performed with radiative MHD numerical models and interface region imaging spectrograph (IRIS) observations. We have expanded a previously existing 3D radiative MHD numerical model of the solar atmosphere, which included self-consistently locally driven magnetic amplification in the chromosphere, by adding ambipolar diffusion and time-dependent nonequilibrium hydrogen ionization to the model. The energy of the magnetic field is dissipated in the upper chromosphere, providing a large temperature increase due to ambipolar diffusion and nonequilibrium ionization (NEQI). At the same time, we find that adding the ambipolar diffusion and NEQI in the simulation has a minor impact on the local growth of the magnetic field in the lower chromosphere and its dynamics. Our comparison between synthesized Mg ii profiles from these high-spatial-resolution models, with and without ambipolar diffusion and NEQI, and quiet-Sun and coronal hole observations from IRIS now reveal a slightly better correspondence. The intensity of profiles is increased, and the line cores are slightly broader when ambipolar diffusion and NEQI effects are included. Therefore, the Mg ii profiles are closer to those observed than in previous models, though some differences still remain.

The ponderomotive force has been suggested to be the main mechanism to produce the so-called first ionization potential (FIP) effect—the enrichment of low-FIP elements observed in the outer solar atmosphere, in the solar wind, and in solar energetic events. It is well known that the ionization of these elements occurs within the chromosphere. Therefore, this phenomenon is intimately tied to the plasma state in the chromosphere and the corona. For this study, we combine IRIS observations, a single-fluid 2.5D radiative magnetohydrodynamics (MHD) model of the solar atmosphere, including ion–neutral interaction effects and nonequilibrium (NEQ) ionization effects, and a novel multifluid multispecies numerical model (based on the Ebysus code). Nonthermal velocities of Si iv measured from IRIS spectra can provide an upper limit for the strength of any high-frequency Alfvén waves. With the single-fluid model, we investigate the possible impact of NEQ ionization within the region where the FIP may occur, as well as the plasma properties in those regions. These models suggest that regions with strongly enhanced network and type II spicules are possible sites of large ponderomotive forces. We use the plasma properties of the single-fluid MHD model and the IRIS observations to initialize our multifluid models to investigate the multifluid effects on the ponderomotive force associated with Alfvén waves. Our multifluid analysis reveals that collisions and NEQ ionization effects dramatically impact the behavior of the ponderomotive force in the chromosphere, and existing theories may need to be revisited.

The Mg ii h and k lines are among the best diagnostic tools of the upper solar chromosphere. This region of the atmosphere is of particular interest, as it is the lowest region of the Sun’s atmosphere where the magnetic field is dominant in the energetics and dynamics, defining its structure. While highly successful in the photosphere and lower to mid-chromosphere, numerical models have produced synthetic Mg ii lines that do not match the observations well. We present a number of large-scale models with magnetic field topologies representative of the quiet Sun, ephemeral flux regions and plage, and also models where the numerical resolution is high and where we go beyond the MHD paradigm. The results of this study show models with a much improved correspondence with IRIS observations in terms of both intensities and widths, especially underscoring the importance of chromospheric mass loading and of capturing the magnetic field topology and evolution in simulations. This comes in addition to the importance of capturing the generation of small-scale velocity fields and including nonequilibrium ionization and ion−neutral interaction effects. However, it should be noted that difficulties in achieving a good correspondence remain, especially when considering the width of Mg ii h and k lines in plage. Understanding and modeling all these effects and their relative importance is necessary in order to reproduce observed spectral features and in isolating the missing pieces necessary to fully comprehend Mg ii formation.

We report high-resolution, high-cadence observations of five small-scale coronal jets in an on-disk quiet Sun region observed with Solar Orbiter’s EUI/HRI _EUV in 174 Å. We combine the HRI _EUV images with the EUV images of SDO/AIA and investigate the magnetic setting of the jets using coaligned line-of-sight magnetograms from SDO/HMI. The HRI _EUV jets are miniature versions of typical coronal jets as they show narrow collimated spires with a base brightening. Three out of five jets result from a detectable minifilament eruption following flux cancelation at the neutral line under the minifilament, analogous to coronal jets. To better understand the physics of jets, we also analyze five small-scale jets from a high-resolution Bifrost MHD simulation in synthetic Fe ix /Fe x emissions. The jets in the simulation reside above neutral lines and four out of five jets are triggered by magnetic flux cancelation. The temperature maps show evidence of cool gas in the same four jets. Our simulation also shows the signatures of opposite Doppler shifts (of the order of ±10 s of km s ^−1 ) in the jet spire, which is evidence of untwisting motion of the magnetic field in the jet spire. The average jet duration, spire length, base width, and speed in our observations (and in synthetic Fe ix /Fe x images) are 6.5 ± 4.0 min (9.0 ± 4.0 minutes), 6050 ± 2900 km (6500 ± 6500 km), 2200 ± 850 km, (3900 ± 2100 km), and 60 ± 8 km s ^−1 (42 ± 20 km s ^−1 ), respectively. Our observation and simulation results provide a unified picture of small-scale solar coronal jets driven by magnetic reconnection accompanying flux cancelation. This picture also aligns well with the most recent reports of the formation and eruption mechanisms of larger coronal jets.

The thermal structure of the chromosphere is regulated through a complex interaction of various heating processes, radiative cooling, and the ionization degree of the plasma. Here, we study the impact on the thermal properties of the chromosphere when including the combined action of nonequilibrium ionization (NEI) of hydrogen and helium and ion–neutral interaction effects. We have performed a 2.5D radiative magnetohydrodynamic simulation using the Bifrost code. This model includes ion–neutral interaction effects by solving the generalized Ohm’ s law (GOL) as well as NEI for hydrogen and helium. The GOL equation includes ambipolar diffusion and the Hall term. We compare this simulation with another simulation that computes the ionization in local thermodynamic equilibrium (LTE) including ion–neutral interaction effects. Our numerical models reveal substantial thermal differences in magneto-acoustic shocks, the wake behind the shocks, spicules, low-lying magnetic loops, and the transition region. In particular, we find that heating through ambipolar diffusion in shock wakes is substantially less efficient, while in the shock fronts themselves it is more efficient, under NEI conditions than when assuming LTE. [ABSTRACT FROM AUTHOR]

The role of type II spicules in the corona has been a much debated topic in recent years. This paper aims to shed light on the impact of type II spicules in the corona using novel 2.5D radiative MHD simulations, including ion–neutral interaction effects with the Bifrost code. We find that the formation of simulated type II spicules, driven by the release of magnetic tension, impacts the corona in various manners. Associated with the formation of spicules, the corona exhibits (1) magneto-acoustic shocks and flows, which supply mass to coronal loops, and (2) transversal magnetic waves and electric currents that propagate at Alfvén speeds. The transversal waves and electric currents, generated by the spicule’s driver and lasting for many minutes, are dissipated and heat the associated loop. These complex interactions in the corona can be connected with blueshifted secondary components in coronal spectral lines (red–blue asymmetries) observed with Hinode/EIS and SOHO/SUMER, as well as the EUV counterpart of type II spicules and propagating coronal disturbances observed with the 171 Å and 193 Å SDO/AIA channels. [ABSTRACT FROM AUTHOR]

We investigate the effects of interactions between ions and neutrals on the chromosphere and overlying corona using 2.5D radiative MHD simulations with the Bifrost code. We have extended the code capabilities implementing ion–neutral interaction effects using the generalized Ohm’s law, i.e., we include the Hall term and the ambipolar diffusion (Pedersen dissipation) in the induction equation. Our models span from the upper convection zone to the corona, with the photosphere, chromosphere, and transition region partially ionized. Our simulations reveal that the interactions between ionized particles and neutral particles have important consequences for the magnetothermodynamics of these modeled layers: (1) ambipolar diffusion increases the temperature in the chromosphere; (2) sporadically the horizontal magnetic field in the photosphere is diffused into the chromosphere, due to the large ambipolar diffusion; (3) ambipolar diffusion concentrates electrical currents, leading to more violent jets and reconnection processes, resulting in (3a) the formation of longer and faster spicules, (3b) heating of plasma during the spicule evolution, and (3c) decoupling of the plasma and magnetic field in spicules. Our results indicate that ambipolar diffusion is a critical ingredient for understanding the magnetothermodynamic properties in the chromosphere and transition region. The numerical simulations have been made publicly available, similar to previous Bifrost simulations. This will allow the community to study realistic numerical simulations with a wider range of magnetic field configurations and physics modules than previously possible. [ABSTRACT FROM AUTHOR]

We study, for the first time, the ascent of granular-sized magnetic bubbles from the solar photosphere through the chromosphere into the transition region and above. Such events occurred in a flux emerging region in NOAA 11850 on 2013 September 25. During that time, the first co-observing campaign between the Swedish 1-m Solar Telescope (SST) and the Interface Region Imaging Spectrograph (IRIS) spacecraft was carried out. Simultaneous observations of the chromospheric Hα 656.28 nm and Ca ii 854.2 nm lines, plus the photospheric Fe i 630.25 nm line, were made with the CRISP spectropolarimeter at the Spitzer Space Telescope (SST) reaching a spatial resolution of 0.″14. At the same time, IRIS was performing a four-step dense raster of the emerging flux region, taking slit jaw images at 133 (C ii, transition region), 140 (Si iv, transition region), 279.6 (Mg ii k, core, upper chromosphere), and 283.2 nm (Mg ii k, wing, photosphere). Spectroscopy of several lines was performed by the IRIS spectrograph in the far- and near-ultraviolet, of which we have used the Si iv 140.3 and the Mg ii k 279.6 nm lines. Coronal images from the Atmospheric Imaging Assembly of the Solar Dynamics Observatory were used to investigate the possible coronal signatures of the flux emergence events. The photospheric and chromospheric properties of small-scale emerging magnetic bubbles have been described in detail in Ortiz et al. Here we are able to follow such structures up to the transition region. We describe the properties, including temporal delays, of the observed flux emergence in all layers. We believe this may be an important mechanism of transporting energy and magnetic flux from subsurface layers to the transition region and corona. [ABSTRACT FROM AUTHOR]

The properties of nonstatistical equilibrium ionization of silicon and oxygen ions are analyzed in this work. We focus on five solar targets (quiet Sun; coronal hole; plage; quiescent active region, AR; and flaring AR) as observed with the Interface Region Imaging Spectrograph (IRIS). IRIS is best suited for this work owing to the high cadence (up to 0.5 s), high spatial resolution (up to 0.″32), and high signal-to-noise ratios for O iv λ1401 and Si iv λ1402. We find that the observed intensity ratio between lines of three times ionized silicon and oxygen ions depends on their total intensity and that this correlation varies depending on the region observed (quiet Sun, coronal holes, plage, or active regions) and on the specific observational objects present (spicules, dynamic loops, jets, microflares, or umbra). In order to interpret the observations, we compare them with synthetic profiles taken from 2D self-consistent radiative MHD simulations of the solar atmosphere, where the statistical equilibrium or nonequilibrium treatment of silicon and oxygen is applied. These synthetic observations show vaguely similar correlations to those in the observations, i.e., between the intensity ratios and their intensities, but only in the nonequilibrium case do we find that (some of) the observations can be reproduced. We conclude that these lines are formed out of statistical equilibrium. We use our time-dependent nonequilibrium ionization simulations to describe the physical mechanisms behind these observed properties. [ABSTRACT FROM AUTHOR]

The Interface Region Imaging Spectrograph (IRIS) reveals small-scale rapid brightenings in the form of bright grains all over coronal holes and the quiet Sun. These bright grains are seen with the IRIS 1330, 1400, and 2796 Å slit-jaw filters. We combine coordinated observations with IRIS and from the ground with the Swedish 1 m Solar Telescope (SST) which allows us to have chromospheric (Ca ii 8542 Å, Ca ii H 3968 Å, Hα, and Mg ii k 2796 Å) and transition region (C ii 1334 Å, Si iv 1403 Å) spectral imaging, and single-wavelength Stokes maps in Fe i 6302 Å at high spatial (), temporal, and spectral resolution. We conclude that the IRIS slit-jaw grains are the counterpart of so-called acoustic grains, i.e., resulting from chromospheric acoustic waves in a non-magnetic environment. We compare slit-jaw images (SJIs) with spectra from the IRIS spectrograph. We conclude that the grain intensity in the 2796 Å slit-jaw filter comes from both the Mg ii k core and wings. The signal in the C ii and Si iv lines is too weak to explain the presence of grains in the 1300 and 1400 Å SJIs and we conclude that the grain signal in these passbands comes mostly from the continuum. Although weak, the characteristic shock signatures of acoustic grains can often be detected in IRIS C ii spectra. For some grains, a spectral signature can be found in IRIS Si iv. This suggests that upward propagating acoustic waves sometimes reach all the way up to the transition region. [ABSTRACT FROM AUTHOR]