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.

In the past few years, the ALMA radio telescope has become available for solar observations. ALMA diagnostics of the solar atmosphere are of high interest because of the theoretically expected linear relationship between the brightness temperature at mm wavelengths and the local gas temperature in the solar atmosphere. Key for the interpretation of solar ALMA observations is understanding where in the solar atmosphere the ALMA emission originates. Recent theoretical studies have suggested that ALMA bands at 1.2 (band 6) and 3 mm (band 3) form in the middle and upper chromosphere at significantly different heights. We study the formation of ALMA diagnostics using a 2.5D radiative MHD model that includes the effects of ion-neutral interactions (ambipolar diffusion) and non-equilibrium ionization of hydrogen and helium. Our results suggest that in active regions and network regions, observations at both wavelengths most often originate from similar heights in the upper chromosphere, contrary to previous results. Non-equilibrium ionization increases the opacity in the chromosphere so that ALMA mostly observe spicules and fibrils along the canopy fields. We combine these modeling results with observations from IRIS, SDO and ALMA to suggest a new interpretation for the recently reported "dark chromospheric holes", regions of very low temperatures in the chromosphere.

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.

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.

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).

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.

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.

Spicules are the most ubuiquitous type of jets in the solar atmosphere. The advent of high-resolution imaging and spectroscopy from the Interface Region Imaging Spectrograph (IRIS) and ground-based observatories has revealed the presence of very high apparent motions of order 100-300 km/s in spicules, as measured in the plane of the sky. However, line-of-sight measurements of such high speeds have been difficult to obtain, with values deduced from Doppler shifts in spectral lines typically of order 30-70 km/s. In this work we resolve this long-standing discrepancy using recent 2.5D radiative MHD simulations. This simulation has revealed a novel driving mechanism for spicules in which ambipolar diffusion resulting from ion-neutral interactions plays a key role. In our simulation we often see that the upward propagation of magnetic waves and electrical currents from the low chromosphere into already existing spicules can lead to rapid heating when the currents are rapidly dissipated by ambipolar diffusion. The combination of rapid heating and the propagation of these currents at Alfv\'enic speeds in excess of 100 km/s leads to the very rapid apparent motions, and often wholesale appearance, of spicules at chromospheric and transition region temperatures. In our simulation, the observed fast apparent motions in such jets are actually a signature of a heating front, and much higher than the mass flows, which are of order 30-70 km/s. Our results can explain the behavior of transition region "network jets" and the very high apparent speeds reported for some chromospheric spicules.

Spectral observations of the solar transition region (TR) and corona show broadening of spectral lines beyond what is expected from thermal and instrumental broadening. The remaining non-thermal broadening is significant (5-30 km/s) and correlated with the intensity. Here we study spectra of the TR Si IV 1403 Angstrom line obtained at high resolution with the Interface Region Imaging Spectrograph (IRIS). We find that the large improvement in spatial resolution (0.33 arcsec) of IRIS compared to previous spectrographs (2 arcsec) does not resolve the non-thermal line broadening which in most regions remains at pre-IRIS levels of about 20 km/s. This invariance to spatial resolution indicates that the processes behind the broadening occur along the line-of-sight (LOS) and/or on spatial scales (perpendicular to the LOS) smaller than 250 km. Both effects appear to play a role. Comparison with IRIS chromospheric observations shows that, in regions where the LOS is more parallel to the field, magneto-acoustic shocks driven from below impact the TR and can lead to significant non-thermal line broadening. This scenario is supported by MHD simulations. While these do not show enough non-thermal line broadening, they do reproduce the long-known puzzling correlation between non-thermal line broadening and intensity. This correlation is caused by the shocks, but only if non-equilibrium ionization is taken into account. In regions where the LOS is more perpendicular to the field, the prevalence of small-scale twist is likely to play a significant role in explaining the invariance and correlation with intensity.

Spicules have been proposed as significant contributors to the mass and energy balance of the corona. While previous observations have provided a glimpse of short-lived transient brightenings in the corona that are associated with spicules, these observations have been contested and are the subject of a vigorous debate both on the modeling and the observational side. Therefore, it remains unclear whether plasma is heated to coronal temperatures in association with spicules. We use high-resolution observations of the chromosphere and transition region with the Interface Region Imaging Spectrograph (IRIS) and of the corona with the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) to show evidence of the formation of coronal structures associated with spicular mass ejections and heating of plasma to transition region and coronal temperatures. Our observations suggest that a significant fraction of the highly dynamic loop fan environment associated with plage regions may be the result of the formation of such new coronal strands, a process that previously had been interpreted as the propagation of transient propagating coronal disturbances (PCD)s. Our observations are supported by 2.5D radiative MHD simulations that show heating to coronal temperatures in association with spicules. Our results suggest that heating and strong flows play an important role in maintaining the substructure of loop fans, in addition to the waves that permeate this low coronal environment.

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.

We study the buoyant rise of magnetic flux tubes embedded in an adiabatic stratification using two-and three-dimensional, MHD simulations. We analyze the dependence of the tube evolution on the field line twist and on the curvature of the tube axis in different diffusion regimes. To be able to achieve a comparatively high spatial resolution we use the FLASH code, which has a built-in Adaptive Mesh Refinement (AMR) capability. Our 3D experiments reach Reynolds numbers that permit a reasonable comparison of the results with those of previous 2D simulations. When the experiments are run without AMR, hence with a comparatively large diffusivity, the amount of longitudinal magnetic flux retained inside the tube increases with the curvature of the tube axis. However, when a low-diffusion regime is reached by using the AMR algorithms, the magnetic twist is able to prevent the splitting of the magnetic loop into vortex tubes and the loop curvature does not play any significant role. We detect the generation of vorticity in the main body of the tube of opposite sign on the opposite sides of the apex. This is a consequence of the inhomogeneity of the azimuthal component of the field on the flux surfaces. The lift force associated with this global vorticity makes the flanks of the tube move away from their initial vertical plane in an antisymmetric fashion. The trajectories have an oscillatory motion superimposed, due to the shedding of vortex rolls to the wake, which creates a Von Karman street.

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.

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.

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 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 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.

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.

The existence of asymmetries in the circular polarization (Stokes V) profiles emerging from the solar photosphere is known since the 1970s. These profiles require the presence of a velocity gradient along the line of sight, possibly associated with gradients of magnetic field strength, inclination and/or azimuth. We have focused our study on the Stokes V profiles showing extreme asymmetry in the from of only one lobe. Using Hinode spectropolarimetric measurements we have performed a statistical study of the properties of these profiles in the quiet sun. We show their spatial distribution, their main physical properties, how they are related with several physical observables and their behavior with respect to their position on the solar disk. The single lobed Stokes V profiles occupy roughly 2% of the solar surface. For the first time, we have observed their temporal evolution and have retrieved the physical conditions of the atmospheres from which they emerged using an inversion code implementing discontinuities of the atmospheric parameters along the line of sight. In addition, we use synthetic Stokes profiles from 3D magnetoconvection simulations to complement the results of the inversion. The main features of the synthetic single-lobed profiles are in general agreement with the observed ones, lending support to the magnetic and dynamic topologies inferred from the inversion. The combination of all these different analysis suggests that most of the single-lobed Stokes V profiles are signals coming from magnetic flux emergence and/or submergence processes taking place in small patches in the photospheric of the quiet sun.

Recent spectroscopic observations have revealed the ubiquitous presence of blueward asymmetries of emission lines formed in the solar corona and transition region. These asymmetries are most prominent in loop footpoint regions, where a clear correlation of the asymmetry with the Doppler shift and line width determined from the single Gaussian fit is found. Such asymmetries suggest at least two emission components: a primary component accounting for the background emission and a secondary component associated with high-speed upflows. The latter has been proposed to play a vital role in the coronal heating process and there is no agreement on its properties. Here we slightly modify the initially developed technique of Red-Blue (RB) asymmetry analysis and apply it to both artificial spectra and spectra observed by the EUV Imaging Spectrometer onboard Hinode, and demonstrate that the secondary component usually contributes a few percent of the total emission, has a velocity ranging from 50 to 150 km s-1 and a Gaussian width comparable to that of the primary one in loop footpoint regions. The results of the RB asymmetry analysis are then used to guide a double Gaussian fit and we find that the obtained properties of the secondary component are generally consistent with those obtained from the RB asymmetry analysis. Through a comparison of the location, relative intensity, and velocity distribution of the blueward secondary component with the properties of the upward propagating disturbances revealed in simultaneous images from the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory, we find a clear association of the secondary component with the propagating disturbances.

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.

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.