We investigate the temperature enhancements and formation heights of impulsive heating phenomena in solar active regions such as Ellerman bombs (EBs), ultraviolet bursts (UVBs), and flaring active-region fibrils (FAFs) using interferometric observations in the millimeter (mm) continuum provided by the Atacama Large Millimeter/submillimeter Array (ALMA). We examined 3 mm signatures of heating events identified in Solar Dynamics Observatory (SDO) observations of an active region and compare the results with synthetic spectra from a 3D radiative magnetohydrodynamic simulation. We estimated the contribution from the corona to the mm brightness using differential emission measure analysis. We report the null detection of EBs in the 3 mm continuum at $\sim1.2$" spatial resolution, which is evidence that they are sub-canopy events that do not significantly contribute to heating the upper chromosphere. In contrast, we find the active region to be populated with multiple compact, bright, flickering mm bursts -- reminiscent of UVBs. The high brightness temperatures of up to $\sim14200$ K in some events may have a significant contribution (up to $\sim$7%) from the corona. We also detect FAF-like events in the 3 mm continuum that show rapid motions of $>10000$ K plasma launched with high plane-of-sky velocities ($37-340\rm\,km\,s^{-1}$) from nanoflare kernels. The mm FAFs are the brightest class of warm canopy fibrils that connect magnetic regions of opposite polarity. The simulation confirms that ALMA should be able to detect the mm counterparts of UVBs and small flares and thus provide a complementary diagnostic for impulsive heating in the solar chromosphere.

NASA's Interface Region Imaging Spectrograph (IRIS) provides high resolution observations of the solar atmosphere through UV spectroscopy and imaging. Since the launch of IRIS in June 2013, we have conducted systematic observation campaigns in coordination with the Swedish 1-m Solar Telescope (SST) on La Palma. The SST provides complementary high-resolution observations of the photosphere and chromosphere. The SST observations include spectro-polarimetric imaging in photospheric Fe I lines and spectrally-resolved imaging in the chromospheric Ca II 8542 A, H-alpha, and Ca II K lines. We present a database of co-aligned IRIS and SST datasets that is open for analysis to the scientific community. The database covers a variety of targets including active regions, sunspots, plage, quiet Sun, and coronal holes.

We present observations of the solar chromosphere obtained simultaneously with the Atacama Large Millimeter/submillimeter Array (\emph{ALMA}) and the Interface Region Imaging Spectrograph (\emph{IRIS}). The observatories targeted a chromospheric plage region of which the spatial distribution (split between strongly and weakly magnetized regions) allowed the study of linear-like structures in isolation, free of contamination from background emission. Using these observations in conjunction with a radiative magnetohydrodynamic 2.5D model covering the upper convection zone all the way to the corona that considers non-equilibrium ionization effects, we report the detection of an on-disk chromospheric spicule with \emph{ALMA} and confirm its multithermal nature. In addition, we discuss the strikingly high degree of similarity between chromospheric plage features observed in \emph{ALMA}/Band6 and \emph{IRIS}/\ion{Si}{4} (also reproduced in our model) suggesting that \emph{ALMA}/Band6 does not observe in the low chromosphere as previously thought but rather observes the upper chromospheric parts of structures such as spicules and other bright structures above plage at geometric heights near transition region temperatures. We also show that \emph{IRIS}/\ion{Mg}{2} is not as well correlated with \emph{ALMA}/Band6 as was previously thought. For these comparisons, we propose and employ a novel empirical method for the determination of plage regions, which seems to better isolate plage from its surrounding regions as compared to other methods commonly used. We caution that isolating plage from its immediate surroundings must be done with care to mitigate stastistical bias in quantitative comparisons between different chromospheric observables. Lastly, we report indications for chromospheric heating due to traveling shocks supported by the \emph{ALMA}/Band6 observations.

Coronal rain consists of cool and dense plasma condensations formed in coronal loops as a result of thermal instability. Previous numerical simulations of thermal instability and coronal rain formation have relied on artificially adding a coronal heating term to the energy equation. To reproduce large-scale characteristics of the corona, using more realistic coronal heating prescription is necessary. We analyse coronal rain formation and evolution in a 3-dimensional radiative magnetohydrodynamic simulation spanning from convection zone to corona which is self-consistently heated by magnetic field braiding as a result of convective motions. We investigate the spatial and temporal evolution of energy dissipation along coronal loops which become thermally unstable. Ohmic dissipation in the model leads to the heating events capable of inducing sufficient chromospheric evaporation into the loop to trigger thermal instability and condensation formation. The cooling of the thermally unstable plasma occurs on timescales comparable to the duration of the individual impulsive heating events. The impulsive heating has sufficient duration to trigger thermal instability in the loop but does not last long enough to lead to coronal rain limit cycles. We show that condensations can either survive and fall into the chromosphere or be destroyed by strong bursts of Joule heating associated with a magnetic reconnection events. In addition, we find that condensations can also form along open magnetic field lines.

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 non-equilibrium ionization (NEI) of hydrogen and helium and ion-neutral interaction effects. We have performed a 2.5D radiative magnetohydrodynamic simulation including ion-neutral interaction effects by solving the generalized Ohm's law (GOL) as well as NEI for hydrogen and helium using the Bifrost code. 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.

This textbook gives a perspective of heliophysics in a way that emphasizes universal processes from a perspective that draws attention to what provides Earth (and similar (exo-)planets) with a relatively stable setting in which life as we know it can thrive. The book is intended for students in physical sciences in later years of their university training and for beginning graduate students in fields of solar, stellar, (exo-)planetary, and planetary-system sciences.

The emergence of magnetic flux into the solar atmosphere produces, among other phenomena, Ellerman bombs (EBs), which are observed in H alpha and are due to magnetic reconnection in the photosphere below the chromospheric canopy. Signs of reconnection are also observed in other spectral lines, typical of the chromosphere or the transition region. An example are the UV bursts observed in the transition region lines of Si IV and the upper chromospheric lines of Mg II. In this work we analyze high cadence, high resolution coordinated observations between the Swedish 1-m Solar Telescope (SST) and the Interface Region Imaging Spectrograph (IRIS) spacecraft. H alpha images from the SST provide us with the positions, timings and trajectories of EBs in an emerging flux region. Simultaneous, co-aligned IRIS slit-jaw images at 133 (C II, transition region), 140 (Si IV, transition region) and 279.6 (Mg II k, core, upper chromosphere) nm, as well as spectroscopy in the far and near ultraviolet from the fast spectrograph raster, allow us to study the possible chromospheric/transition region counterparts of those EBs. Our main goal is to study the possible temporal and spatial relationship between several reconnection events at different layers in the atmosphere (namely EBs and UV bursts), the timing history between them, and the connection of these dynamical phenomena to the ejection of surges in the chromosphere. We also investigate the properties of an extended UV burst and their variations across the burst domain. Our results suggest a scenario where simultaneous and co-spatial EBs and UV bursts are part of the same reconnection system occurring sequentially along a vertical or nearly vertical current sheet. Heating and bidirectional jets trace the location where reconnection takes place. These results support and expand those obtained from recent numerical simulations of magnetic flux emergence.

The Multi-slit Solar Explorer (MUSE) is a proposed mission aimed at understanding the physical mechanisms driving the heating of the solar corona and the eruptions that are at the foundation of space weather. MUSE contains two instruments, a multi-slit EUV spectrograph and a context imager. It will simultaneously obtain EUV spectra (along 37 slits) and context images with the highest resolution in space (0.33-0.4 arcsec) and time (1-4 s) ever achieved for the transition region and corona. The MUSE science investigation will exploit major advances in numerical modeling, and observe at the spatial and temporal scales on which competing models make testable and distinguishable predictions, thereby leading to a breakthrough in our understanding of coronal heating and the drivers of space weather. By obtaining spectra in 4 bright EUV lines (Fe IX 171A, Fe XV 284A, Fe XIX-XXI 108A) covering a wide range of transition region and coronal temperatures along 37 slits simultaneously, MUSE will be able to "freeze" the evolution of the dynamic coronal plasma. We describe MUSE's multi-slit approach and show that the optimization of the design minimizes the impact of spectral lines from neighboring slits, generally allowing line parameters to be accurately determined. We also describe a Spectral Disambiguation Code to resolve multi-slit ambiguity in locations where secondary lines are bright. We use simulations of the corona and eruptions to perform validation tests and show that the multi-slit disambiguation approach allows accurate determination of MUSE observables in locations where significant multi-slit contamination occurs.

The Spectral Imaging of the Coronal Environment (SPICE) instrument is a high-resolution imaging spectrometer operating at extreme ultraviolet (EUV) wavelengths. In this paper, we present the concept, design, and pre-launch performance of this facility instrument on the ESA/NASA Solar Orbiter mission. The goal of this paper is to give prospective users a better understanding of the possible types of observations, the data acquisition, and the sources that contribute to the instrument's signal. The paper discusses the science objectives, with a focus on the SPICE-specific aspects, before presenting the instrument's design, including optical, mechanical, thermal, and electronics aspects. This is followed by a characterisation and calibration of the instrument's performance. The paper concludes with descriptions of the operations concept and data processing. The performance measurements of the various instrument parameters meet the requirements derived from the mission's science objectives. The SPICE instrument is ready to perform measurements that will provide vital contributions to the scientific success of the Solar Orbiter mission.

Ellerman bombs (EBs) and UV bursts are both brightenings related to flux emergence regions and specifically to magnetic flux of opposite polarity that meet in the photosphere. These two reconnection-related phenomena, nominally formed far apart, occasionally occur in the same location and at the same time, thus challenging our understanding of reconnection and heating of the lower solar atmosphere. We consider the formation of an active region, including long fibrils and hot and dense coronal plasma. The emergence of a untwisted magnetic flux sheet, injected $2.5$~Mm below the photosphere, is studied as it pierces the photosphere and interacts with the preexisting ambient field. Specifically, we aim to study whether EBs and UV bursts are generated as a result of such flux emergence and examine their physical relationship. The Bifrost radiative magnetohydrodynamics code was used to model flux emerging into a model atmosphere that contained a fairly strong ambient field, constraining the emerging field to a limited volume wherein multiple reconnection events occur as the field breaks through the photosphere and expands into the outer atmosphere. Synthetic spectra of the different reconnection events were computed using the $1.5$D RH code and the fully 3D MULTI3D code. The formation of UV bursts and EBs at intensities and with line profiles that are highly reminiscent of observed spectra are understood to be a result of the reconnection of emerging flux with itself in a long-lasting current sheet that extends over several scale heights through the chromosphere. Synthetic diagnostics suggest that there are no compelling reasons to assume that UV bursts occur in the photosphere. Instead, EBs and UV bursts are occasionally formed at opposite ends of a long current sheet that resides in an extended bubble of cool gas.

The presence of magnetic field is crucial in the transport of energy through the solar atmosphere. Recent ground-based and space-borne observations of the quiet Sun have revealed that magnetic field accumulates at photospheric heights, via a local dynamo or from small-scale flux emergence events. However, most of this small-scale magnetic field may not expand into the chromosphere due to the entropy drop with height at the photosphere. Here we present a study that uses a high resolution 3D radiative MHD simulation of the solar atmosphere with non-grey and non-LTE radiative transfer and thermal conduction along the magnetic field to reveal that: 1) the net magnetic flux from the simulated quiet photosphere is not sufficient to maintain a chromospheric magnetic field (on average), 2) processes in the lower chromosphere, in the region dominated by magneto-acoustic shocks, are able to convert kinetic energy into magnetic energy, 3) the magnetic energy in the chromosphere increases linearly in time until the r.m.s. of the magnetic field strength saturates at roughly 4 to 30 G (horizontal average) due to conversion from kinetic energy, 4) and that the magnetic features formed in the chromosphere are localized to this region.

The signal measured by an astronomical spectrometer may be due to radiation from a multi-component mixture of plasmas with a range of physical properties (e.g. temperature, Doppler velocity). Confusion between multiple components may be exacerbated if the spectrometer sensor is illuminated by overlapping spectra dispersed from different slits, with each slit being exposed to radiation from a different portion of an extended astrophysical object. We use a compressed sensing method to robustly retrieve the different components. This method can be adopted for a variety of spectrometer configurations, including single-slit, multi-slit (e.g., the proposed MUlti-slit Solar Explorer mission; MUSE) and slot spectrometers (which produce overlappograms).

Small bipolar magnetic features are observed to appear in the interior of individual granules in the quiet Sun, signaling the emergence of tiny magnetic loops from the solar interior. We study the origin of those features as part of the magnetoconvection process in the top layers of the convection zone. Two quiet-Sun magnetoconvection models, calculated with the radiation-magnetohydrodynamic (MHD) Bifrost code and with domain stretching from the top layers of the convection zone to the corona, are analyzed. Using 3D visualization as well as a posteriori spectral synthesis of Stokes parameters, we detect the repeated emergence of small magnetic elements in the interior of granules, as in the observations. Additionally, we identify the formation of organized horizontal magnetic sheets covering whole granules. Our approach is twofold, calculating statistical properties of the system, like joint probability density functions (JPDFs), and pursuing individual events via visualization tools. We conclude that the small magnetic loops surfacing within individual granules in the observations may originate from sites at or near the downflows in the granular and mesogranular levels, probably in the first 1 or 1.5 Mm below the surface. We also document the creation of granule-covering magnetic sheet-like structures through the sideways expansion of a small subphotospheric magnetic concentration picked up, and pulled out of the interior, by a nascent granule. The sheet-like structures we found in the models may match the recent observations of Centeno et al. (2017).

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\'en 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 blue shifted 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 (PCDs) observed with the 171~\AA\ and 193~\AA\ SDO/AIA channels.

The measured average velocities in solar and stellar spectral lines formed at transition region temperatures have been difficult to interpret. However, realistic three-dimensional radiation magnetohydrodynamics (3D rMHD) models of the solar atmosphere are able to reproduce the observed dominant line shifts and may thus hold the key to resolve these issues. Our new 3D rMHD simulations aim to shed light on how mass flows between the chromosphere and corona and on how the coronal mass is maintained. Passive tracer particles, so-called corks, allow the tracking of parcels of plasma over time and thus the study of changes in plasma temperature and velocity not only locally, but also in a co-moving frame. By following the trajectories of the corks, we can investigate mass and energy flows and understand the composition of the observed velocities. Our findings show that most of the transition region mass is cooling. The preponderance of transition region redshifts in the model can be explained by the higher percentage of downflowing mass in the lower and middle transition region. The average upflows in the upper transition region can be explained by a combination of both stronger upflows than downflows and a higher percentage of upflowing mass. The most common combination at lower and middle transition region temperatures are corks that are cooling and traveling downward. For these corks, a strong correlation between the pressure gradient along the magnetic field line and the velocity along the magnetic field line has been observed, indicating a formation mechanism that is related to downward propagating pressure disturbances. Corks at upper transition region temperatures are subject to a rather slow and highly variable but continuous heating process.

Low-lying loops have been discovered at the solar limb in transition region temperatures by the Interface Region Imaging Spectrograph (IRIS). They do not appear to reach coronal temperatures, and it has been suggested that they are the long-predicted unresolved fine structures (UFS). These loops are dynamic and believed to be visible during both heating and cooling phases. Making use of coordinated observations between IRIS and the Swedish 1-m Solar Telescope, we study how these loops impact the solar chromosphere. We show for the first time that there is indeed a chromospheric signal of these loops, seen mostly in the form of strong Doppler shifts and a conspicuous lack of chromospheric heating. In addition, we find that several instances have a inverse Y-shaped jet just above the loop, suggesting that magnetic reconnection is driving these events. Our observations add several puzzling details to the current knowledge of these newly discovered structures; this new information must be considered in theoretical models.

CHROMIS, a new imaging spectrometer at the Swedish 1-m Solar Telescope (SST), can observe the chromosphere in the H and K lines of Ca II at high spatial and spectral resolution. Accurate modeling as well as an understanding of the formation of these lines are needed to interpret the SST/CHROMIS observations. Such modeling is computationally challenging because these lines are influenced by strong departures from local thermodynamic equilibrium, three-dimensional radiative transfer, and partially coherent resonance scattering of photons. We aim to model the CaII H&K lines in 3D model atmospheres to understand their formation and to investigate their diagnostic potential for probing the chromosphere. We model the synthetic spectrum of Ca II using the radiative transfer code Multi3D in three different radiation-magnetohydrodynamic model atmospheres computed with the Bifrost code. We classify synthetic intensity profiles according to their shapes and study how their features are related to the physical properties in the model atmospheres. We investigate whether the synthetic data reproduce the observed spatially-averaged line shapes, center-to-limb variation and compare with SST/CHROMIS images. The spatially-averaged synthetic line profiles show too low central emission peaks, and too small separation between the peaks. The trends of the observed center-to-limb variation of the profiles properties are reproduced by the models. The Ca II H&K line profiles provide a temperature diagnostic of the temperature minimum and the temperature at the formation height of the emission peaks. The Doppler shift of the central depression is an excellent probe of the velocity in the upper chromosphere.

Magnetic reconnection is thought to drive a wide variety of dynamic phenomena in the solar atmosphere. Yet the detailed physical mechanisms driving reconnection are difficult to discern in the remote sensing observations that are used to study the solar atmosphere. In this paper we exploit the high-resolution instruments Interface Region Imaging Spectrograph (IRIS) and the new CHROMIS Fabry-Perot instrument at the Swedish 1-m Solar Telescope (SST) to identify the intermittency of magnetic reconnection and its association with the formation of plasmoids in so-called UV bursts in the low solar atmosphere. The Si IV 1403A UV burst spectra from the transition region show evidence of highly broadened line profiles with often non-Gaussian and triangular shapes, in addition to signatures of bidirectional flows. Such profiles had previously been linked, in idealized numerical simulations, to magnetic reconnection driven by the plasmoid instability. Simultaneous CHROMIS images in the chromospheric Ca II K 3934A line now provide compelling evidence for the presence of plasmoids, by revealing highly dynamic and rapidly moving brightenings that are smaller than 0.2 arcsec and that evolve on timescales of order seconds. Our interpretation of the observations is supported by detailed comparisons with synthetic observables from advanced numerical simulations of magnetic reconnection and associated plasmoids in the chromosphere. Our results highlight how subarcsecond imaging spectroscopy sensitive to a wide range of temperatures combined with advanced numerical simulations that are realistic enough to compare with observations can directly reveal the small-scale physical processes that drive the wide range of phenomena in the solar atmosphere.

In the lower solar atmosphere, the chromosphere is permeated by jets known as spicules, in which plasma is propelled at speeds of 50 to 150 kilometers per second into the corona. The origin of the spicules is poorly understood, although they are expected to play a role in heating the million-degree corona and are associated with Alfv\'enic waves that help drive the solar wind. We compare magnetohydrodynamic simulations of spicules with observations from the Interface Region Imaging Spectrograph and the Swedish 1-m Solar Telescope. Spicules are shown to occur when magnetic tension is amplified and transported upward through interactions between ions and neutrals or ambipolar diffusion. The tension is impulsively released to drive flows, heat plasma (through ambipolar diffusion), and generate Alfv\'enic waves.

Observations of the upper chromosphere shows an enormous amount of intricate fine structure. Much of this comes in the form of linear features which are most often assumed to be well aligned with the direction of the magnetic field in the low plasma beta regime thought to dominate the upper chromosphere. We use advanced radiative MHD simulations including the effects of ion-neutral interactions (using the generalized Ohm's law) in the partially ionized chromosphere to show that the magnetic field is often not well aligned with chromospheric features. This occurs where the ambipolar diffusion is large, i.e., ions and neutral populations decouple as the ion-neutral collision frequency drops allowing the field to slip through the neutral population, currents perpendicular to the field are strong, and thermodynamic timescales are longer than or similar to the those of ambipolar diffusion. We find this often happens in dynamic spicule or fibril-like features at the top of the chromosphere. This has important consequences for field extrapolation methods which increasingly use such upper chromospheric features to help constrain the chromospheric magnetic field: our results invalidate the underlying assumption that these features are aligned with the field. In addition, our results cast doubt on results from 1D hydrodynamic models, which assume that plasma remains on the same field lines. Finally, our simulations show that ambipolar diffusion significantly alters the amount of free energy available in the coronal part of our simulated volume, which is likely to have consequences for studies of flare initiation.

It is typically assumed that emission in the passbands of the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) is dominated by single or several strong lines from ions that under equilibrium conditions are formed in a narrow range of temperatures. However, most SDO/AIA channels also contain contributions from lines of ions that have formation temperatures that are significantly different from the "dominant" ion(s). We investigate the importance of these lines by forward modeling the emission in the SDO/AIA channels with 3D radiative MHD simulations of a model that spans the upper layer of the convection zone to the low corona. The model is highly dynamic. In addition, we pump a steadily increasing magnetic flux into the corona, in order to increase the coronal temperature through the dissipation of magnetic stresses. As a consequence, the model covers different ranges of coronal temperatures as time progresses. The model covers coronal temperatures that are representative of plasma conditions in coronal holes and quiet sun. The 131, 171, and 304 \AA{} AIA passbands are found to be least influenced by the so-called "non-dominant" ions, and the emission observed in these channels comes mostly from plasma at temperatures near the formation temperature of the dominant ion(s). On the other hand, the other channels are strongly influenced by the non-dominant ions, and therefore significant emission in these channels comes from plasma at temperatures that are different from the "canonical" values. We have also studied the influence of non-dominant ions on the AIA passbands when different element abundances are assumed (photospheric and coronal), and when the effects of the electron density on the contribution function are taken into account.

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 magneto-thermodynamics 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 magneto-thermo-dynamic 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.

A spectacular manifestation of solar activity is the appearance of transient brightenings in the far wings of the H$\alpha$ line, known as Ellerman bombs (EBs). Recent observations obtained by the Interface Region Imaging Spectrograph (IRIS) have revealed another type of plasma "bombs" (UV bursts) with high temperatures of perhaps up to $8\times 10^4$ K within the cooler lower solar atmosphere. Realistic numerical modeling showing such events is needed to explain their nature. Here, we report on 3D radiative magneto-hydrodynamic simulations of magnetic flux emergence in the solar atmosphere. We find that ubiquitous reconnection between emerging bipolar magnetic fields can trigger EBs in the photosphere, UV bursts in the mid/low chromosphere and small (nano-/micro-) flares ($10^6$ K) in the upper chromosphere. These results provide new insights on the emergence and build up of the coronal magnetic field and the dynamics and heating of the solar surface and lower atmosphere.

Using radiation magnetohydrodynamic simulations of the solar atmospheric layers from the upper convection zone to the lower corona, we investigate the self-consistent excitation of slow magneto-acoustic body waves (slow modes) in a magnetic flux concentration. We find that the convective downdrafts in the close surroundings of a two-dimensional flux slab "pump" the plasma inside it in the downward direction. This action produces a downflow inside the flux slab, which encompasses ever higher layers, causing an upwardly propagating rarefaction wave. The slow mode, excited by the adiabatic compression of the downflow near the optical surface, travels along the magnetic field in the upward direction at the tube speed. It develops into a shock wave at chromospheric heights, where it dissipates, lifts the transition region, and produces an offspring in the form of a compressive wave that propagates further into the corona. In the wake of downflows and propagating shock waves, the atmosphere inside the flux slab in the chromosphere and higher tends to oscillate with a period of $\nu\approx 4$~mHz. We conclude that this process of "magnetic pumping" is a most plausible mechanism for the direct generation of longitudinal chromospheric and coronal compressive waves within magnetic flux concentrations, and it may provide an important heat source in the chromosphere. It may also be responsible for certain types of dynamic fibrils.

We use UV spectral observations of active regions with the Interface Region Imaging Spectrograph (IRIS) to investigate the properties of the coronal FeXII 1349.4A emission at unprecedented high spatial resolution (~0.33"). We find that by using appropriate observational strategies (i.e., long exposures, lossless compression), FeXII emission can be studied with IRIS at high spatial and spectral resolution, at least for high density plasma (e.g., post-flare loops, and active region moss). We find that upper transition region (moss) FeXII emission shows very small average Doppler redshifts (v_Dop ~3 km/s), as well as modest non-thermal velocities (with an average ~24 km/s, and the peak of the distribution at ~15 km/s). The observed distribution of Doppler shifts appears to be compatible with advanced 3D radiative MHD simulations in which impulsive heating is concentrated at the transition region footpoints of a hot corona. While the non-thermal broadening of FeXII 1349.4A peaks at similar values as lower resolution simultaneous Hinode/EIS measurements of FeXII 195A, IRIS observations show a previously undetected tail of increased non-thermal broadening that might be suggestive of the presence of subarcsecond heating events. We find that IRIS and EIS non-thermal line broadening measurements are affected by instrumental effects that can only be removed through careful analysis. Our results also reveal an unexplained discrepancy between observed 195.1/1349.4A FeXII intensity ratios and those predicted by the CHIANTI atomic database.

We study the ascent of granular-sized magnetic bubbles from the solar photosphere through the chromosphere into the transition region and above, for the first time. Such events occurred in a flux emerging region in NOAA 11850 on September 25, 2013. During that time, the first co-observing campaign between the Swedish 1-m Solar Telescope and the IRIS spacecraft was carried out. Simultaneous observations of the chromospheric H$\alpha$ 656.28 nm and \ion{Ca}{2} 854.2 nm lines, plus the photospheric \ion{Fe}{1} 630.25 nm line, were made with the CRISP spectropolarimeter at the SST reaching a spatial resolution of 0."14. At the same time, IRIS was performing a four-step dense raster of the said emerging flux region, taking slit-jaw images at 133 (C~{\sc ii}, transition region), 140 (\ion{Si}{4}, transition region), 279.6 (\ion{Mg}{2} k, core, upper chromosphere), and 283.2 nm (\ion{Mg}{2} 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 \ion{Si}{4} 140.3 and the \ion{Mg}{2} 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. (2014). 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.

Active regions (AR) appearing on the surface of the Sun are classified into $\alpha$, $\beta$, $\gamma$, and $\delta$ by the rules of the Mount Wilson Observatory, California on the basis of their topological complexity. Amongst these, the $\delta$-sunspots are known to be super-active and produce the most X-ray flares. Here, we present results from a simulation of the Sun by mimicking the upper layers and the corona, but starting at a more primitive stage than any earlier treatment. We find that this initial state consisting of only a thin sub-photospheric magnetic sheet breaks into multiple flux-tubes which evolve into a colliding-merging system of spots of opposite polarity upon surface emergence, similar to those often seen on the Sun. The simulation goes on to produce many exotic $\delta$-sunspot associated phenomena: repeated flaring in the range of typical solar flare energy release and ejective helical flux ropes with embedded cool-dense plasma filaments resembling solar coronal mass ejections.

The properties of non-statistical equilibrium ionization of silicon and oxygen ions are analyzed in this work. We focus on four 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 due to the high cadence (up to 0.5s), high spatial resolution (up to 0.32"), and high signal to noise ratios for O IV and Si IV. 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, micro-flares 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 non-equilibrium treatment silicon and oxygen is applied. These synthetic observations show vaguely similar correlations as in the observations, i.e. between the intensity ratios and their intensities, but only in the non-equilibrium 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 non-equilibrium ionization simulations to describe the physical mechanisms behind these observed properties.

Context. The solar chromosphere is the interface between the solar surface and the solar corona. Modelling of this region is difficult because it represents the transition from optically thick to thin radiation escape, from gas-pressure domination to magnetic-pressure domination, from a neutral to an ionised state, from MHD to plasma physics, and from near-equilibrium (LTE) to non-equilibrium conditions. Aims. Our aim is to provide the community with realistic simulations of the magnetic solar outer atmosphere. This will enable detailed comparison of existing and upcoming observations with synthetic observables from the simulations, thereby elucidating the complex interactions of magnetic fields and plasma that are crucial for our understanding of the dynamic outer atmosphere. Methods. We used the radiation magnetohydrodynamics code Bifrost to perform simulations of a computational volume with a magnetic field topology similar to an enhanced network area on the Sun. Results. The full simulation cubes are made available online. The general properties of the simulation are discussed, and limitations are discussed.

Advanced 3D radiative MHD simulations now reproduce many properties of the outer solar atmosphere. When including a domain from the convection zone into the corona, a hot chromosphere and corona are self-consistently maintained. Here we study two realistic models, with different simulated area, magnetic field strength and topology, and numerical resolution. These are compared in order to characterize the heating in the 3D-MHD simulations which self-consistently maintains the structure of the atmosphere. We analyze the heating at both large and small scales and find that heating is episodic and highly structured in space, but occurs along loop shaped structures, and moves along with the magnetic field. On large scales we find that the heating per particle is maximal near the transition region and that widely distributed opposite-polarity field in the photosphere leads to a greater heating scale height in the corona. On smaller scales, heating is concentrated in current sheets, the thicknesses of which are set by the numerical resolution. Some current sheets fragment in time, this process occurring more readily in the higher-resolution model leading to spatially highly intermittent heating. The large scale heating structures are found to fade in less than about five minutes, while the smaller, local, heating shows time scales of the order of 2 minutes in one model and 1 minutes in the other, higher-resolution, model.

A triplet of subordinate lines of Mg II exists in the region around the h&k lines. In solar spectra these lines are seen mostly in absorption, but in some cases can become emission lines. The aim of this work is to study the formation of this triplet, and investigate any diagnostic value they can bring. Using 3D radiative magnetohydrodynamic simulations of quiet Sun and flaring flux emergence, we synthesize spectra and investigate how spectral features respond to the underlying atmosphere. We find that emission in the lines is rare and is typically caused by a steep temperature increase in the lower chromosphere (above 1500 K, with electron densities above 10$^{18}$ m$^{-3}$). In both simulations the lines are sensitive to temperature increases taking place at column masses >= 5e-4 g cm$^{-2}$. Additional information can also be inferred from the peak-to-wing ratio and shape of the line profiles. Using observations from NASA's Interface Region Imaging Spectrograph we find both absorption and emission line profiles with similar shapes to the synthetic spectra, which suggests that these lines represent a useful diagnostic that complements the MgII h&k lines.

Magnetic flux emergence into the outer layers of the Sun is a fundamental mechanism for releasing energy into the chromosphere and the corona. In this paper, we study the emergence of granular-sized flux concentrations and the structuring of the corresponding physical parameters and atmospheric diagnostics in the upper photo- sphere and in the chromosphere. We make use of a realistic 3D MHD simulation of the outer layers of the Sun to study the formation of the Ca II 8542 line. We also derive semi-empirical 3D models from non-LTE inversions of our observations. These models contain depth-dependent information of the temperature and line-of-sight stratification. Our analysis explains the peculiar Ca II 8542 profiles observed in the flux-emerging region. In addition, we derive detailed temperature and velocity maps describing the ascent of magnetic bubbles from the photosphere to the chromosphere. The inversions suggest that, in active regions, granular-sized bubbles emerge up to the lower chromosphere where the existing large-scale field hinders their ascent. We report hints of heating when the field reaches the chromosphere.

The energy for the coronal heating must be provided from the convection zone. The amount and the method by which this energy is transferred into the corona depends on the properties of the lower atmosphere and the corona itself. We review: 1) how the energy could be built in the lower solar atmosphere; 2) how this energy is transferred through the solar atmosphere; and 3) how the energy is finally dissipated in the chromosphere and/or corona. Any mechanism of energy transport has to deal with the various physical processes in the lower atmosphere. We will focus on a physical process that seems to be highly important in the chromosphere and not deeply studied until recently: the ion-neutral interaction effects (INIE) in the chromosphere. We review the relevance and the role of the partial ionization in the chromosphere and show that this process actually impacts considerably the outer solar atmosphere. We include analysis of our 2.5D radiative MHD simulations with the Bifrost code (Gudiksen et al. 2011) including the partial ionization effects on the chromosphere and corona and thermal conduction along magnetic field lines. The photosphere, chromosphere and transition region are partially ionized and the interaction between ionized particles and neutral particles has important consequences on the magneto-thermodynamics of these layers. The INIE are treated using generalized Ohm's law, i.e., we consider the Hall term and the ambipolar diffusion in the induction equation. The interaction between the different species affects the modeled atmosphere as follows: 1) the ambipolar diffusion dissipates magnetic energy and increases the minimum temperature in the chromosphere; 2) the upper chromosphere may get heated and expanded over a greater range of heights. These processes reveal appreciable differences between the modeled atmospheres of simulations with and without INIE.

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 \AA, 1400 \AA\ and 2796 \AA\ 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 \AA, Ca II H 3968 \AA, H\alpha, and Mg II k 2796 \AA), and transition region (C II 1334 \AA, Si IV 1402) spectral imaging, and single-wavelength Stokes maps in Fe I 6302 \AA at high spatial (0.33"), 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 with spectra from the IRIS spectrograph. We conclude that the grain intensity in the 2796 \AA\ 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 \AA\ slit-jaw images and we conclude that the grain signal in these passbands comes mostly from the continuum. Even though weak, the characteristic shock signatures of acoustic grains can often be detected in IRIS C II spectra. For some grains, 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.

We report on observations of recurrent jets by instruments onboard the Interface Region Imaging Spectrograph (IRIS), Solar Dynamics Observatory (SDO) and Hinode spacecrafts. Over a 4-hour period on July 21st 2013, recurrent coronal jets were observed to emanate from NOAA Active Region 11793. FUV spectra probing plasma at transition region temperatures show evidence of oppositely directed flows with components reaching Doppler velocities of +/- 100 km/s. Raster Doppler maps using a Si IV transition region line show all four jets to have helical motion of the same sense. Simultaneous observations of the region by SDO and Hinode show that the jets emanate from a source region comprising a pore embedded in the interior of a supergranule. The parasitic pore has opposite polarity flux compared to the surrounding network field. This leads to a spine-fan magnetic topology in the coronal field that is amenable to jet formation. Time-dependent data-driven simulations are used to investigate the underlying drivers for the jets. These numerical experiments show that the emergence of current-carrying magnetic field in the vicinity of the pore supplies the magnetic twist needed for recurrent helical jet formation.

As the interface between the Sun's photosphere and corona, the chromosphere and transition region play a key role in the formation and acceleration of the solar wind. Observations from the Interface Region Imaging Spectrograph reveal the prevalence of intermittent small-scale jets with speeds of 80-250 km/s from the narrow bright network lanes of this interface region. These jets have lifetimes of 20-80 seconds and widths of 300 km or less. They originate from small-scale bright regions, often preceded by footpoint brightenings and accompanied by transverse waves with ~20 km/s amplitudes. Many jets reach temperatures of at least ~100000 K and constitute an important element of the transition region structures. They are likely an intermittent but persistent source of mass and energy for the solar wind.

We use coordinated observations with the Interface Region Imaging Spectrograph (IRIS) and the Swedish 1-m Solar Telescope (SST) to identify the disk counterpart of type II spicules in upper-chromospheric and transition region (TR) diagnostics. These disk counterparts were earlier identified through short-lived asymmetries in chromospheric spectral lines: rapid blue- or red-shifted excursions (RBEs or RREs). We find clear signatures of RBEs and RREs in Mg II h & k, often with excursions of the central h3 and k3 absorption features in concert with asymmetries in co-temporal and co-spatial H-alpha spectral profiles. We find spectral signatures for RBEs and RREs in C II 1335 and 1336 A and Si IV 1394 and 1403 A spectral lines and interpret this as a sign that type II spicules are heated to at least TR temperatures, supporting other recent work. These C II and Si IV spectral signals are weaker for a smaller network region than for more extended network regions in our data. A number of bright features around extended network regions observed in IRIS slit-jaw imagery SJI 1330 and 1400, recently identified as network jets, can be clearly connected to H-alpha RBEs and/or RREs in our coordinated data. We speculate that at least part of the diffuse halo around network regions in the IRIS SJI 1330 and 1400 images can be attributed to type II spicules with insufficient opacity in the C II and Si IV lines to stand out as single features in these passbands.

The heating of the outer solar atmospheric layers, i.e., the transition region and corona, to high temperatures is a long standing problem in solar (and stellar) physics. Solutions have been hampered by an incomplete understanding of the magnetically controlled structure of these regions. The high spatial and temporal resolution observations with the Interface Region Imaging Spectrograph (IRIS) at the solar limb reveal a plethora of short, low lying loops or loop segments at transition-region temperatures that vary rapidly, on the timescales of minutes. We argue that the existence of these loops solves a long standing observational mystery. At the same time, based on comparison with numerical models, this detection sheds light on a critical piece of the coronal heating puzzle.

The solar chromosphere and transition region (TR) form an interface between the Sun's surface and its hot outer atmosphere. Here most of the non-thermal energy that powers the solar atmosphere is transformed into heat, although the detailed mechanism remains elusive. High-resolution (0.33-arcsec) observations with NASA's Interface Region Imaging Spectrograph (IRIS) reveal a chromosphere and TR that are replete with twist or torsional motions on sub-arcsecond scales, occurring in active regions, quiet Sun regions, and coronal holes alike. We coordinated observations with the Swedish 1-m Solar Telescope (SST) to quantify these twisting motions and their association with rapid heating to at least TR temperatures. This view of the interface region provides insight into what heats the low solar atmosphere.

The physical processes causing energy exchange between the Sun's hot corona and its cool lower atmosphere remain poorly understood. The chromosphere and transition region (TR) form an interface region between the surface and the corona that is highly sensitive to the coronal heating mechanism. High resolution observations with the Interface Region Imaging Spectrograph (IRIS) reveal rapid variability (about 20 to 60 seconds) of intensity and velocity on small spatial scales at the footpoints of hot dynamic coronal loops. The observations are consistent with numerical simulations of heating by beams of non-thermal electrons, which are generated in small impulsive heating events called "coronal nanoflares". The accelerated electrons deposit a sizable fraction of their energy in the chromosphere and TR. Our analysis provides tight constraints on the properties of such electron beams and new diagnostics for their presence in the nonflaring corona.

The solar atmosphere was traditionally represented with a simple one-dimensional model. Over the past few decades, this paradigm shifted for the chromosphere and corona that constitute the outer atmosphere, which is now considered a dynamic structured envelope. Recent observations by IRIS (Interface Region Imaging Spectrograph) reveal that it is difficult to determine what is up and down even in the cool 6000-K photosphere just above the solar surface: this region hosts pockets of hot plasma transiently heated to almost 100,000 K. The energy to heat and accelerate the plasma requires a considerable fraction of the energy from flares, the largest solar disruptions. These IRIS observations not only confirm that the photosphere is more complex than conventionally thought, but also provide insight into the energy conversion in the process of magnetic reconnection.

Solar spicules have eluded modelers and observers for decades. Since the discovery of the more energetic type II, spicules have become a heated topic but their contribution to the energy balance of the low solar atmosphere remains unknown. Here we give a first glimpse of what quiet Sun spicules look like when observed with NASA's recently launched Interface Region Imaging Spectrograph (IRIS). Using IRIS spectra and filtergrams that sample the chromosphere and transition region we compare the properties and evolution of spicules as observed in a coordinated campaign with Hinode and the Atmospheric Imaging Assembly. Our IRIS observations allow us to follow the thermal evolution of type II spicules and finally confirm that the fading of Ca II H spicules appears to be caused by rapid heating to higher temperatures. The IRIS spicules do not fade but continue evolving, reaching higher and falling back down after 500-800 s. Ca II H type II spicules are thus the initial stages of violent and hotter events that mostly remain invisible in Ca II H filtergrams. These events have very different properties from type I spicules, which show lower velocities and no fading from chromospheric passbands. The IRIS spectra of spicules show the same signature as their proposed disk counterparts, reinforcing earlier work. Spectroheliograms from spectral rasters also confirm that quiet Sun spicules originate in bushes from the magnetic network. Our results suggest that type II spicules are indeed the site of vigorous heating (to at least transition region temperatures) along extensive parts of the upward moving spicular plasma.

IRIS data allow us to study the solar transition region (TR) with an unprecedented spatial resolution of 0.33 arcsec. On 2013 August 30, we observed bursts of high Doppler shifts suggesting strong supersonic downflows of up to 200 km/s and weaker, slightly slower upflows in the spectral lines Mg II h and k, C II 1336 \AA, Si IV 1394 \AA, and 1403 \AA, that are correlated with brightenings in the slitjaw images (SJIs). The bursty behavior lasts throughout the 2 hr observation, with average burst durations of about 20 s. The locations of these short-lived events appear to be the umbral and penumbral footpoints of EUV loops. Fast apparent downflows are observed along these loops in the SJIs and in AIA, suggesting that the loops are thermally unstable. We interpret the observations as cool material falling from coronal heights, and especially coronal rain produced along the thermally unstable loops, which leads to an increase of intensity at the loop footpoints, probably indicating an increase of density and temperature in the TR. The rain speeds are on the higher end of previously reported speeds for this phenomenon, and possibly higher than the free-fall velocity along the loops. On other observing days, similar bright dots are sometimes aligned into ribbons, resembling small flare ribbons. These observations provide a first insight into small-scale heating events in sunspots in the TR.

We present the first results of sunspot oscillations from observations by the Interface Region Imaging Spectrograph. The strongly nonlinear oscillation is identified in both the slit-jaw images and the spectra of several emission lines formed in the transition region and chromosphere. We first apply a single Gaussian fit to the profiles of the Mgii 2796.35 {\AA}, Cii 1335.71 {\AA}, and Si iv 1393.76 {\AA} lines in the sunspot. The intensity change is about 30%. The Doppler shift oscillation reveals a sawtooth pattern with an amplitude of about 10 km/s in Si iv. In the umbra the Si iv oscillation lags those of Cii and Mgii by about 3 and 12 s, respectively. The line width suddenly increases as the Doppler shift changes from redshift to blueshift. However, we demonstrate that this increase is caused by the superposition of two emission components. We then perform detailed analysis of the line profiles at a few selected locations on the slit. The temporal evolution of the line core is dominated by the following behavior: a rapid excursion to the blue side, accompanied by an intensity increase, followed by a linear decrease of the velocity to the red side. The maximum intensity slightly lags the maximum blueshift in Si iv, whereas the intensity enhancement slightly precedes the maximum blueshift in Mgii. We find a positive correlation between the maximum velocity and deceleration, a result that is consistent with numerical simulations of upward propagating magnetoacoustic shock waves.

We report on the formation of small solar flares produced by patchy magnetic reconnection between interacting magnetic loops. A three-dimensional (3D) magnetohydrodynamic (MHD) numerical experiment was performed, where a uniform magnetic flux sheet was injected into a fully developed convective layer. The gradual emergence of the field into the solar atmosphere results in a network of magnetic loops, which interact dynamically forming current layers at their interfaces. The formation and ejection of plasmoids out of the current layers leads to patchy reconnection and the spontaneous formation of several small (size ? 1-2Mm) flares. We find that these flares are short-lived (30 s - 3 min) bursts of energy in the range O(10^25 - 10^27) ergs, which is basically the nanoflare-microflare range. Their persistent formation and co-operative action and evolution leads to recurrent emission of fast EUV/X-ray jets and considerable plasma heating in the active corona.

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.

We study a granular-sized magnetic flux emergence event that occurred in NOAA 11024 in July 2009. The observations were made with the CRISP spectropolarimeter at the Swedish 1 m Solar Telescope achieving a spatial resolution of 0.14". Simultaneous full Stokes observations of the two photospheric Fe I lines at 630.2 nm and the chromospheric Ca II 854.2 nm line allow us to describe in detail the emergence process across the solar atmosphere. We report here on 3D semi-spherical bubble events, where instead of simple magnetic footpoints, we observe complex semi-circular feet straddling a few granules. The most characteristic signature in these events is the observation of a dark bubble in filtergrams taken in the wings of the Ca II 854.2 nm line. We can infer how the bubble rises through the solar atmosphere as we see it progressing from the wings to the core of Ca II 854.2 nm. In the photosphere, the magnetic bubble shows mean upward Doppler velocities of 2 km/s. In about 3.5 minutes it travels some 1100 km to reach the mid chromosphere, implying an average ascent speed of 5.2 km/s. To aid the interpretation of the observations, we carry out 3D numerical simulations of the evolution of a horizontal, untwisted magnetic flux sheet injected in the convection zone, using the Bifrost code. The computational domain spans from the upper convection zone to the lower corona. In the modeled chromosphere the rising flux sheet produces a large, cool, magnetized bubble. We compare this bubble with the observed ones and find excellent agreement, including similar field strengths and velocity signals in the photosphere and chromosphere, temperature deficits, ascent speeds, expansion velocities, and lifetimes.

We performed a 3D radiative MHD simulation of the solar atmosphere. This simulation shows a jet-like feature that shows similarities to the type II spicules observed for the first time with Hinode. Rapid Blueshifted Events (RBEs) on the solar disk are associated with these spicules. Observational results suggest they may contribute significantly in supplying the corona with hot plasma. We perform a detailed comparison of the properties of the simulated jet with those of type II spicules (observed with Hinode) and RBEs (with ground-based instruments). We analyze variety of synthetic emission and absorption lines from the simulations including chromospheric Ca II and Ha to TR and coronal temperatures (10E4 to several 10E6K). We compare their synthetic intensities, line profiles, Doppler shifts, line widths and asymmetries with observations from Hinode/SOT and EIS, SOHO/SUMER, SST and SDO/AIA. Many properties of the synthetic observables resemble the observations, and we describe in detail the physical processes that lead to these observables. Detailed analysis of the synthetic observables provides insight into how observations should be analyzed to derive information about physical variables in such a dynamic event. For example, we find that LOS superposition in the optically thin atmosphere requires the combination of Doppler shifts and spectral line asymmetry to determine the velocity in the jet. Other properties differ from the observations, especially in the chromospheric lines. The mass density of the part of the spicule with a chromospheric temperature is too low to produce significant opacity in chromospheric lines. These and other discrepancies are described in detail, and we discuss which mechanisms and physical processes may need to be included in the MHD simulations to mimic the thermodynamic processes of the chromosphere and corona, in particular to reproduce type II spicules.

The High-resolution Coronal Imager (Hi-C) has provided Fe XII 193A images of the upper transition region moss at an unprecedented spatial (~0.3-0.4 arcsec) and temporal (5.5s) resolution. The Hi-C observations show in some moss regions variability on timescales down to ~15s, significantly shorter than the minute scale variability typically found in previous observations of moss, therefore challenging the conclusion of moss being heated in a mostly steady manner. These rapid variability moss regions are located at the footpoints of bright hot coronal loops observed by SDO/AIA in the 94A channel, and by Hinode/XRT. The configuration of these loops is highly dynamic, and suggestive of slipping reconnection. We interpret these events as signatures of heating events associated with reconnection occurring in the overlying hot coronal loops, i.e., coronal nanoflares. We estimate the order of magnitude of the energy in these events to be of at least a few $10^{23}rg, also supporting the nanoflare scenario. These Hi-C observations suggest that future observations at comparable high spatial and temporal resolution, with more extensive temperature coverage are required to determine the exact characteristics of the heating mechanism(s).

Determining the temperature distribution of coronal plasmas can provide stringent constraints on coronal heating. Current observations with the Extreme ultraviolet Imaging Spectrograph onboard Hinode and the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory provide diagnostics of the emission measure distribution (EMD) of the coronal plasma. Here we test the reliability of temperature diagnostics using 3D radiative MHD simulations. We produce synthetic observables from the models, and apply the Monte Carlo Markov chain EMD diagnostic. By comparing the derived EMDs with the "true" distributions from the model we assess the limitations of the diagnostics, as a function of the plasma parameters and of the signal-to-noise of the data. We find that EMDs derived from EIS synthetic data reproduce some general characteristics of the true distributions, but usually show differences from the true EMDs that are much larger than the estimated uncertainties suggest, especially when structures with significantly different density overlap along the line-of-sight. When using AIA synthetic data the derived EMDs reproduce the true EMDs much less accurately, especially for broad EMDs. The differences between the two instruments are due to the: (1) smaller number of constraints provided by AIA data, (2) broad temperature response function of the AIA channels which provide looser constraints to the temperature distribution. Our results suggest that EMDs derived from current observatories may often show significant discrepancies from the true EMDs, rendering their interpretation fraught with uncertainty. These inherent limitations to the method should be carefully considered when using these distributions to constrain coronal heating.

Spicules are long, thin, highly dynamic features that jut out ubiquitously from the solar limb. They dominate the interface between the chromosphere and corona and may provide significant mass and energy to the corona. We use high-quality observations with the Swedish 1-m Solar Telescope to establish that so-called type II spicules are characterized by the simultaneous action of three different types of motion: 1. field-aligned flows of order 50-100 km/s, 2. swaying motions of order 15-20 km/s, and 3. torsional motions of order 25-30 km/s. The first two modes have been studied in detail before, but not the torsional motions. Our analysis of many near-limb and off-limb spectra and narrow-band images yields strong evidence that most, if not all, type-II spicules undergo large torsional modulation and that these motions, like spicule swaying, represent Alfvenic waves propagating outward at several hundred km/s. The combined action of the different motions explains the similar morphology of spicule bushes in the outer red and blue wings of chromospheric lines, and needs to be taken into account when interpreting Doppler motions to derive estimates for field-aligned flows in spicules and determining the Alfvenic wave energy in the solar atmosphere. Our results also suggest that large torsional motion is an ingredient in the production of type-II spicules and that spicules play an important role in the transport of helicity through the solar atmosphere.

The solar chromosphere is weakly ionized and interactions between ionized particles and neutral particles likely have significant consequences for the thermodynamics of the plasma. We investigate the importance of introducing neutral particles using numerical 2.5D radiative MHD simulations obtained with the Bifrost code. The models span from the upper layers of the convection zone to the low corona, and solve the full MHD equations with non-grey and NLTE radiative transfer, and thermal conduction. The effects of partial ionization are implemented using the generalized Ohm's law. The approximations required in going from three fluids to the generalized Ohm's law are tested in our simulations. The Ohmic diffusion, the Hall term, and ambipolar diffusion show strong variations in the chromosphere. These strong variations of the various magnetic diffusivities are absent or significantly underestimated when, as has been common for these types of studies, using the VAL-C model as a basis for estimates. In addition, we find that differences in estimating the magnitude of ambipolar diffusion arise depending on which method is used to calculate the ion-neutral collision frequency. These differences cause uncertainties in the different magnetic diffusivity terms. In the chromosphere, we find that the ambipolar diffusion is of the same order of magnitude or even larger than the numerical diffusion used to stabilize our code. As a consequence, ambipolar diffusion produces a strong impact on the modeled atmosphere. Perhaps more importantly, it suggests that at least in the chromospheric domain, self-consistent simulations of the solar atmosphere driven by magneto-convection can accurately describe the impact of the dominant form of resistivity (ambipolar diffusion). This suggests that such simulations may be more realistic in their approach to the lower solar atmosphere than previously assumed.

We present the results of numerical simulations of wave propagation and jet formation in solar atmosphere models with different magnetic field configurations. The presence in the chromosphere of waves with periods longer than the acoustic cutoff period has been ascribed to either strong inclined magnetic fields, or changes in the radiative relaxation time. Our simulations include a sophisticated treatment of radiative losses, as well as fields with different strengths and inclinations. Using Fourier and wavelet analysis techniques, we investigate the periodicity of the waves that travel through the chromosphere. We find that the velocity signal is dominated by waves with periods around 5 minutes in regions of strong, inclined field, including at the edges of strong flux tubes where the field expands, whereas 3-minute waves dominate in regions of weak or vertically oriented fields. Our results show that the field inclination is very important for long-period wave propagation, whereas variations in the radiative relaxation time have little effect. Furthermore, we find that atmospheric conditions can vary significantly on timescales of a few minutes, meaning that a Fourier analysis of wave propagation can be misleading. Wavelet techniques take variations with time into account and are more suitable analysis tools. Finally, we investigate the properties of jets formed by the propagating waves once they reach the transition region, and find systematic differences between the jets in inclined field regions and those in vertical field regions, in agreement with observations of dynamic fibrils.

Context: Numerical simulations of stellar convection and photospheres have been developed to the point where detailed shapes of observed spectral lines can be explained. Stellar atmospheres are very complex, and very different physical regimes are present in the convection zone, photosphere, chromosphere, transition region and corona. To understand the details of the atmosphere it is necessary to simulate the whole atmosphere since the different layers interact strongly. These physical regimes are very diverse and it takes a highly efficient massively parallel numerical code to solve the associated equations. Aims: The design, implementation and validation of the massively parallel numerical code Bifrost for simulating stellar atmospheres from the convection zone to the corona. Methods: The code is subjected to a number of validation tests, among them the Sod shock tube test, the Orzag-Tang colliding shock test, boundary condition tests and tests of how the code treats magnetic field advection, chromospheric radiation, radiative transfer in an isothermal scattering atmosphere, hydrogen ionization and thermal conduction. Results: Bifrost completes the tests with good results and shows near linear efficiency scaling to thousands of computing cores.

have forced the definition of a new type of spicule, "type II's", that are characterized by rising rapidly, having short lives, and by fading away at the end of their lifetimes. Here, we report on features found in realistic 3D simulations of the outer solar atmosphere that resemble the observed type II spicules. These features evolve naturally from the simulations as a consequence of the magnetohydrodynamical evolution of the model atmosphere. The simulations span from the upper layer of the convection zone to the lower corona and include the emergence of horizontal magnetic flux. The state-of-art Oslo Staggered Code (OSC) is used to solve the full MHD equations with non-grey and non-LTE radiative transfer and thermal conduction along the magnetic field lines. We describe in detail the physics involved in a process which we consider a possible candidate as a driver mechanism to produce type II spicules. The modeled spicule is composed of material rapidly ejected from the chromosphere that rises into the corona while being heated. Its source lies in a region with large field gradients and intense electric currents, which lead to a strong Lorentz force that squeezes the chromospheric material, resulting in a vertical pressure gradient that propels the spicule along the magnetic field, as well as Joule heating, which heats the the jet material, forcing it to fade.

aims: We aim to provide an estimate of the minimum temperature of the quiet solar chromosphere. methods: We perform a 2D radiation-MHD simulation spanning the upper convection zone to the lower corona. The simulation includes non-LTE radiative transfer and {an equation-of-state that includes non-equilibrium ionization of hydrogen and non-equilibrium H_2 molecule formation}. We analyze the reliability of the various assumptions made in our model in order to assess the realism of the simulation. results: Our simulation contains pockets of cool gas with down to 1660 K from 1 Mm up to 3.2 Mm height. It overestimates the radiative heating, and contains non-physical heating below 1660 K. Therefore we conclude that cool pockets in the quiet solar chromosphere might have even lower temperatures than in the simulation, provided that there exist areas in the chromosphere without significant magnetic heating. We suggest off-limb molecular spectroscopy to look for such cool pockets and 3D simulations including a local dynamo and a magnetic carpet to investigate Joule heating in the quiet chromosphere.

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.

We present the implementation of a radiative transfer solver with coherent scattering in the new BIFROST code for radiative magneto-hydrodynamical (MHD) simulations of stellar surface convection. The code is fully parallelized using MPI domain decomposition, which allows for large grid sizes and improved resolution of hydrodynamical structures. We apply the code to simulate the surface granulation in a solar-type star, ignoring magnetic fields, and investigate the importance of coherent scattering for the atmospheric structure. A scattering term is added to the radiative transfer equation, requiring an iterative computation of the radiation field. We use a short-characteristics-based Gauss-Seidel acceleration scheme to compute radiative flux divergences for the energy equation. The effects of coherent scattering are tested by comparing the temperature stratification of three 3D time-dependent hydrodynamical atmosphere models of a solar-type star: without scattering, with continuum scattering only, and with both continuum and line scattering. We show that continuum scattering does not have a significant impact on the photospheric temperature structure for a star like the Sun. Including scattering in line-blanketing, however, leads to a decrease of temperatures by about 350\,K below log tau < -4. The effect is opposite to that of 1D hydrostatic models in radiative equilibrium, where scattering reduces the cooling effect of strong LTE lines in the higher layers of the photosphere. Coherent line scattering also changes the temperature distribution in the high atmosphere, where we observe stronger fluctuations compared to a treatment of lines as true absorbers.

An empirical model has been developed to reproduce the drift of the spectrum recorded by EIS on board Hinode using instrumental temperatures and relative motion of the spacecraft. The EIS spectrum shows an artificial drift in wavelength dimension in sync with the revolution of the spacecraft, which is caused by temperature variations inside the spectrometer. The drift amounts to 70 km s$^{-1}$ in Doppler velocity and introduces difficulties in velocity measurements. An artificial neural network is incorporated to establish a relationship between the instrumental temperatures and the spectral drift. This empirical model reproduces observed spectrum shift with an rms error of 4.4 km s$^{-1}$. This procedure is robust and applicable to any spectrum obtained with EIS, regardless of of the observing field. In addition, spectral curvatures and spatial offset in the North - South direction are determined to compensate for instrumental effects.

We present evidence that transition region red-shifts are naturally produced in episodically heated models where the average volumetric heating scale height lies between that of the chromospheric pressure scale height of 200 km and the coronal scale height of 50 Mm. In order to do so we present results from 3d MHD models spanning the upper convection zone up to the corona, 15 Mm above the photosphere. Transition region and coronal heating in these models is due both the stressing of the magnetic field by photospheric and convection `zone dynamics, but also in some models by the injection of emerging magnetic flux.

We have performed 3D radiation MHD simulations extending from the convection zone to the corona covering a box 16 Mm$^3$ at 32 km spatial resolution. The simulations show very fine structure in the chromosphere with acoustic shocks interacting with the magnetic field. Magnetic flux concentrations have a temperature lower than the surroundings in the photosphere but higher in the low chromosphere. The heating is there mostly through ohmic dissipation preferentially at the edges of the flux concentrations. The magnetic field is often wound up around the flux concentrations. When acoustic waves travel up along the field this topology leads to swirling motions seen in chromospheric diagnostic lines such as the calcium infrared triplet.

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.

Many models for coronal loops have difficulty explaining the observed EUV brightness of the transition region, which is often significantly less than theoretical models predict. This discrepancy has been addressed by a variety of approaches including filling factors and time-dependent heating. Here we focus on an effect that has been ignored so far: the absorption of EUV light with wavelengths below 912 {\AA} by the resonance continua of neutral hydrogen and helium. Such absorption is expected to occur in the low-lying transition region of hot, active region loops, that is co-located with cool chromospheric features and called ``moss'' as a result of the reticulated appearance resulting from the absorption. We use co-temporal and co-spatial spectroheliograms obtained with SOHO/SUMER and Hinode/EIS of Fe XII 1242 {\AA}, 195 {\AA} and 186.88 {\AA}, and compare the density determination from the 186/195 {\AA} line ratio to that resulting from the 195/1242 {\AA} line ratio. We find significant absorption of 195 {\AA} emission caused by the chromospheric inclusions in the moss. We find that the amount of absorption is generally of order a factor of 2. We compare to numerical models and show that the observed effect is well reproduced by 3D radiative MHD models of the transition region and corona. We use STEREO A/B data of the same active region and find that increased angles between line-of-sight and local vertical cause additional absorption. Our determination of the amount of chromospheric absorption of TR emission can be used to better constrain coronal heating models.

3D simulations of magnetic flux emergence are carried out in a computational domain spanning the upper layers of the convection zone to the lower corona. We use the Oslo Staggered Code to solve the full MHD equations with non-grey and NLTE radiative transfer and thermal conduction along the magnetic field lines. In this paper we concentrate on the later stages of the simulations and study the evolution of the structure of the rising flux in the upper chromosphere and corona, the interaction between the emerging flux and the weak coronal magnetic field initially present, and the associated dynamics. The flux tube injected at the bottom boundary rises to the photosphere where it largely remains. However, some parts of the flux tube become unstable and expand in patches into the upper chromosphere. The flux rapidly expands towards the corona, pushing the coronal and transition region material aside, lifting and maintaining the transition region at heights greater than 5 Mm above the photosphere for extensive periods of time. The pre-existing magnetic field in the corona and transition region is perturbed by the incoming flux and reoriented by a series of high Joule heating events. Low density structures form in the corona while at later times a high density filamentary structure appears in the lower part of the expanding flux. The dynamics of these and other structures is discussed. While Joule heating due to the expanding flux is episodic, it increases in relative strength as fresh magnetic field rises and becomes energetically important in the upper chromosphere and corona at later times. Chromospheric, transition region and coronal lines are computed and their response to the perturbation caused by the expanding emerging flux is discussed.

The Sun's corona is millions of degrees hotter than its 5,000 K photosphere. This heating enigma is typically addressed by invoking the deposition at coronal heights of non-thermal energy generated by the interplay between convection and magnetic field near the photosphere. However, it remains unclear how and where coronal heating occurs and how the corona is filled with hot plasma. We show that energy deposition at coronal heights cannot be the only source of coronal heating, by revealing a significant coronal mass supply mechanism that is driven from below, in the chromosphere. We quantify the asymmetry of spectral lines observed with Hinode and SOHO and identify faint but ubiquitous upflows with velocities that are similar (50-100 km/s) across a wide range of magnetic field configurations and for temperatures from 100,000 to several million degrees. These upflows are spatio-temporally correlated with and have similar upward velocities as recently discovered, cool (10,000 K) chromospheric jets or (type II) spicules. We find these upflows to be pervasive and universal. Order of magnitude estimates constrained by conservation of mass and observed emission measures indicate that the mass supplied by these spicules can play a significant role in supplying the corona with hot plasma. The properties of these events are incompatible with coronal loop models that only include nanoflares at coronal heights. Our results suggest that a significant part of the heating and energizing of the corona occurs at chromospheric heights, in association with chromospheric jets.

We analyze features that resemble type i spicules in two different 3D numerical simulations in which we include horizontal magnetic flux emergence in a computational domain spanning the upper layers of the convection zone to the lower corona. The two simulations differ mainly in the preexisting ambient magnetic field strength and in the properties of the inserted flux tube. We use the Oslo Staggered Code (OSC) to solve the full MHD equations with non-grey and non-LTE radiative transfer and thermal conduction along the magnetic field lines. We find a multitude of features that show a spatiotemporal evolution that is similar to that observed in type i spicules, which are characterized by parabolic height vs. time profiles, and are dominated by rapid upward motion at speeds of 10-30 km/s, followed by downward motion at similar velocities. We measured the parameters of the parabolic profile of the spicules and find similar correlations between the parameters as those found in observations. The values for height (or length) and duration of the spicules found in the simulations are more limited in range than those in the observations. The spicules found in the simulation with higher preexisting ambient field have shorter length and smaller velocities. From the simulations, it appears that these kinds of spicules can, in principle, be driven by a variety of mechanisms that include p-modes, collapsing granules, magnetic energy release in the photosphere and lower chromosphere and convective buffeting of flux concentrations.

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.

We present the results of numerical simulations of wave-induced magnetic reconnection in a model of the solar atmosphere. In the magnetic field geometry we study in this article, the waves, driven by a monochromatic piston and a driver taken from Hinode observations, induce periodic reconnection of the magnetic field, and this reconnection appears to help drive long-period chromospheric jets. By synthesizing observations for a variety of wavelengths that are sensitive to a wide range of temperatures, we shed light on the often confusing relationship between the plethora of jet-like phenomena in the solar atmosphere, e.g., explosive events, spicules, blinkers, and other phenomena thought to be caused by reconnection.

Recently, De Pontieu et al. (2007b) discovered a class of spicules that evolves more rapidly than previously known spicules, with rapid apparent motions of 50--150 km s${}^{-1}$, thickness of a few 100 km, and lifetimes of order 10--60 seconds. These so-called type II spicules have been difficult to study because of limited spatio-temporal and thermal resolution. Here we use the IBIS instrument to search for the high velocities in the disk counterpart of type II spicules. We have detected rapidly evolving events, with lifetimes that are less than a minute and often equal to the cadence of the instrument (19 secs). These events are characterized by a Doppler shift that only appears in the blue wing of the Ca II IR line. Furthermore the spatial extent, lifetime, and location near network, all suggest a link to type II spicules. However, the magnitude of the measured Doppler velocity is significantly lower than the apparent motions seen at the limb. We use Monte Carlo simulations to show that this discrepancy can be explained by a forward model in which the visibility on the disk of the high-velocity flows in these events is limited by a combination of line-of-sight projection and reduced opacity in upward propelled plasma, especially in reconnection driven jets that are powered by a roughly constant energy supply.

3D numerical simulations of a horizontal magnetic flux tube emergence with different twist are carried out in a computational domain spanning the upper layers of the convection zone to the lower corona. We use the Oslo Staggered Code to solve the full MHD equations with non-grey and non-LTE radiative transfer and thermal conduction along the magnetic field lines. The emergence of the magnetic flux tube input at the bottom boundary into a weakly magnetized atmosphere is presented. The photospheric and chromospheric response is described with magnetograms, synthetic images and velocity field distributions. The emergence of a magnetic flux tube into such an atmosphere results in varied atmospheric responses. In the photosphere the granular size increases when the flux tube approaches from below. In the convective overshoot region some 200km above the photosphere adiabatic expansion produces cooling, darker regions with the structure of granulation cells. We also find collapsed granulation in the boundaries of the rising flux tube. Once the flux tube has crossed the photosphere, bright points related with concentrated magnetic field, vorticity, high vertical velocities and heating by compressed material are found at heights up to 500km above the photosphere. At greater heights in the magnetized chromosphere, the rising flux tube produces a cool, magnetized bubble that tends to expel the usual chromospheric oscillations. In addition the rising flux tube dramatically increases the chromospheric scale height, pushing the transition region and corona aside such that the chromosphere extends up to 6Mm above the photosphere. The emergence of magnetic flux tubes through the photosphere to the lower corona is a relatively slow process, taking of order 1 hour.

We use coordinated Hinode SOT/EIS observations that include high-resolution magnetograms, chromospheric and TR imaging and TR/coronal spectra in a first test to study how the dynamics of the TR are driven by the highly dynamic photospheric magnetic fields and the ubiquitous chromospheric waves. Initial analysis shows that these connections are quite subtle and require a combination of techniques including magnetic field extrapolations, frequency-filtered time-series and comparisons with synthetic chromospheric and TR images from advanced 3D numerical simulations. As a first result, we find signatures of magnetic flux emergence as well as 3 and 5 mHz wave power above regions of enhanced photospheric magnetic field in both chromospheric, transition region and coronal emission.

We use high-resolution observations of the Sun in Ca II H 3968 A from the Solar Optical Telescope on Hinode to show that there are at least two types of spicules that dominate the structure of the magnetic solar chromosphere. Both types are tied to the relentless magnetoconvective driving in the photosphere, but have very different dynamic properties. ``Type-I'' spicules are driven by shock waves that form when global oscillations and convective flows leak into the upper atmosphere along magnetic field lines on 3-7 minute timescales. ``Type-II'' spicules are much more dynamic: they form rapidly (in ~10s), are very thin (<200km wide), have lifetimes of 10-150s (at any one height) and seem to be rapidly heated to (at least) transition region temperatures, sending material through the chromosphere at speeds of order 50-150 km/s. The properties of Type II spicules suggest a formation process that is a consequence of magnetic reconnection, typically in the vicinity of magnetic flux concentrations in plage and network. Both types of spicules are observed to carry Alfven waves with significant amplitudes of order 20 km/s.

We present high spatial resolution spectroscopic measurements of dynamic fibrils (DFs) in the Ca {\small{II}} 8662 {\AA} line. These data show clear Doppler shifts in the identified DFs, which demonstrates that at least a subset of DFs are actual mass motions in the chromosphere. A statistical analysis of 26 DFs reveals a strong and statistically significant correlation between the maximal velocity and the deceleration. The range of the velocities and the decelerations are substantially lower, about a factor two, in our spectroscopic observations compared to the earlier results based on proper motion in narrow band images. There are fundamental differences in the different observational methods; when DFs are observed spectroscopically the measured Doppler shifts are a result of the atmospheric velocity, weighted with the response function to velocity over an extended height. When the proper motion of DFs is observed in narrow band images, the movement of the top of the DF is observed. This point is sharply defined because of the high contrast between the DF and the surroundings. The observational differences between the two methods are examined by several numerical experiments using both numerical simulations and a time series of narrow band H$\alpha$ images. With basis in the simulations we conclude that the lower maximal velocity is explained by the low formation height of the Ca IR line. We conclude that the present observations support the earlier result that DFs are driven by magneto-acoustic shocks exited by convective flows and p-modes.

The ionization of hydrogen in the solar chromosphere and transition region does not obey LTE or instantaneous statistical equilibrium because the timescale is long compared with important hydrodynamical timescales, especially of magneto-acoustic shocks. We implement an algorithm to compute non-equilibrium hydrogen ionization and its coupling into the MHD equations within an existing radiation MHD code, and perform a two-dimensional simulation of the solar atmosphere from the convection zone to the corona. Analysis of the simulation results and comparison to a companion simulation assuming LTE shows that: a) Non-equilibrium computation delivers much smaller variations of the chromospheric hydrogen ionization than for LTE. The ionization is smaller within shocks but subsequently remains high in the cool intershock phases. As a result, the chromospheric temperature variations are much larger than for LTE because in non-equilibrium, hydrogen ionization is a less effective internal energy buffer. The actual shock temperatures are therefore higher and the intershock temperatures lower. b) The chromospheric populations of the hydrogen n = 2 level, which governs the opacity of Halpha, are coupled to the ion populations. They are set by the high temperature in shocks and subsequently remain high in the cool intershock phases. c) The temperature structure and the hydrogen level populations differ much between the chromosphere above photospheric magnetic elements and above quiet internetwork. d) The hydrogen n = 2 population and column density are persistently high in dynamic fibrils, suggesting that these obtain their visibility from being optically thick in Halpha also at low temperature.

We use Hinode/SOT Ca II H-line and blue continuum broadband observations to study the presence and power of high frequency acoustic waves at high spatial resolution. We find that there is no dominant power at small spatial scales; the integrated power using the full resolution of Hinode (0.05'' pixels, 0.16'' resolution) is larger than the power in the data degraded to 0.5'' pixels (TRACE pixel size) by only a factor of 1.2. At 20 mHz the ratio is 1.6. Combining this result with the estimates of the acoustic flux based on TRACE data of Fossum & Carlsson (2006), we conclude that the total energy flux in acoustic waves of frequency 5-40 mHz entering the internetwork chromosphere of the quiet Sun is less than 800 W m$^{-2}$, inadequate to balance the radiative losses in a static chromosphere by a factor of five.

A major goal in solar physics has during the last five decades been to find how energy flux generated in the solar convection zone is transported and dissipated in the outer solar layers. Progress in this field has been slow and painstaking. However, advances in computer hardware and numerical methods, vastly increased observational capabilities and growing physical insight seem finally to be leading towards understanding. Here we present exploratory numerical MHD models that span the entire solar atmosphere from the upper convection zone to the lower corona. These models include non-grey, non-LTE radiative transport in the photosphere and chromosphere, optically thin radiative losses as well as magnetic field-aligned heat conduction in the transition region and corona.

We present high spatial and high temporal resolution observations of the quiet Sun in H-alpha obtained with the Swedish 1-m Solar Telescope on La Palma. We observe that many mottles, jet-like features in the quiet Sun, display clear up- and downward motions along their main axis. In addition, many mottles show vigorous transverse displacements. Unique identification of the mottles throughout their lifetime is much harder than for their active region counterpart, dynamic fibrils. This is because many seem to lack a sharply defined edge at their top, and significant fading often occurs throughout their lifetime. For those mottles that can be reliably tracked, we find that the mottle tops often undergo parabolic paths. We find a linear correlation between the deceleration these mottles undergo and the maximum velocity they reach, similar to what was found earlier for dynamic fibrils. Combined with an analysis of oscillatory properties, we conclude that at least part of the quiet Sun mottles are driven by magnetoacoustic shocks. In addition, the mixed polarity environment and vigorous dynamics suggest that reconnection may play a significant role in the formation of some quiet Sun jets.

We present the results of numerical simulations of shock wave-driven jets in the solar atmosphere. The dependence of observable quantities like maximum velocity and deceleration on parameters such as the period and amplitude of initial disturbances and the inclination of the magnetic field is investigated. Our simulations show excellent agreement with observations, and shed new light on the correlation between velocity and deceleration and on the regional differences found in observations.

With the recent advent of the Swedish 1-m Solar Telescope (SST), advanced image processing techniques, as well as numerical simulations that provide a more realistic view of the chromosphere, a comprehensive understanding of chromospheric jets such as spicules, mottles and fibrils is now within reach. In this paper, we briefly summarize results from a recent analysis of dynamic fibrils, short-lived jet-like features that dominate the chromosphere (as imaged in H-alpha) above and about active region plage. Using extremely high-resolution observations obtained at the SST, and advanced numerical 2D radiative MHD simulations, we show that fibrils are most likely formed by chromospheric shock waves that occur when convective flows and global oscillations leak into the chromosphere along the field lines of magnetic flux concentrations. In addition, we present some preliminary observations of quiet Sun jets or mottles. We find that the mechanism that produces fibrils in active regions is most likely also at work in quiet Sun regions, although it is modified by the weaker magnetic field and the presence of more mixed-polarity. A comparison with numerical simulations suggests that the weaker magnetic field in quiet Sun allows for significantly stronger (than in active regions) transverse motions that are superposed on the field-aligned, shock-driven motions. This leads to a more dynamic, and much more complex environment than in active region plage. In addition, our observations of the mixed polarity environment in quiet Sun regions suggest that other mechanisms, such as reconnection, may well play a significant role in the formation of some quiet Sun jets.

We present unprecedented high resolution H-alpha observations, obtained with the Swedish 1-m Solar Telescope, that, for the first time, spatially and temporally resolve dynamic fibrils in active regions on the Sun. These jet-like features are similar to mottles or spicules in quiet Sun. We find that most of these fibrils follow almost perfect parabolic paths in their ascent and descent. We measure the properties of the parabolic paths taken by 257 fibrils, and present an overview of the deceleration, maximum velocity, maximum length and duration, as well as their widths and the thickness of a bright ring that often occurs above dynamic fibrils. We find that the observed deceleration of the projected path is typically only a fraction of solar gravity, and incompatible with a ballistic path at solar gravity. We report on significant differences of fibril properties between those occurring above a dense plage region, and those above a less dense plage region where the magnetic field seems more inclined from the vertical. We compare these findings to advanced numerical 2D radiative MHD simulations, and find that fibrils are most likely formed by chromospheric shock waves that occur when convective flows and global oscillations leak into the chromosphere along the field lines of magnetic flux concentrations. Detailed comparison of observed and simulated fibril properties shows striking similarities of the values for deceleration, maximum velocity, maximum length and duration. We compare our results with observations of mottles and find that a similar mechanism is most likely at work in the quiet Sun.

The formation of jets such as dynamic fibrils, mottles, and spicules in the solar chromosphere is one of the most important, but also most poorly understood, phenomena of the Sun's magnetized outer atmosphere. We use extremely high-resolution observations from the Swedish 1-m Solar Telescope combined with advanced numerical modeling to show that in active regions these jets are a natural consequence of upwardly propagating slow mode magnetoacoustic shocks. These shocks form when waves generated by convective flows and global p-mode oscillations in the lower lying photosphere leak upward into the magnetized chromosphere. We find excellent agreement between observed and simulated jet velocities, decelerations, lifetimes and lengths. Our findings suggest that previous observations of quiet sun spicules and mottles may also be interpreted in light of a shock driven mechanism.

This work addresses the problem of plasma condensation and ``catastrophic cooling'' in solar coronal loops. We have carried out numerical calculations of coronal loops and find several classes of time-dependent solutions (static, periodic, irregular), depending on the spatial distribution of a temporally constant energy deposition in the loop. Dynamic loops exhibit recurrent plasma condensations, accompanied by high-speed downflows and transient brightenings of transition region lines, in good agreement with features observed with TRACE. Furthermore, these results also offer an explanation for the recent EIT observations of De Groof et al. (2004) of moving bright blobs in large coronal loops. In contrast to earlier models, we suggest that the process of catastrophic cooling is not initiated by a drastic decrease of the total loop heating but rather results from a loss of equilibrium at the loop apex as a natural consequence of heating concentrated at the footpoints of the loop, but constant in time.