Context: Hot coronal jets are a basic observed feature of the solar atmosphere whose physical origin is still being actively debated. Aims: We study six recurrent jets occurring in the active region NOAA 12644 on April 04, 2017. They are observed in all the hot filters of AIA as well as cool surges in IRIS slit-jaw high spatial and temporal resolution images. Methods: The AIA filters allow us to study the temperature and the emission measure of the jets using the filter ratio method. We study the pre-jet phases by analyzing the intensity oscillations at the base of the jets with the wavelet technique. Results: A fine co-alignment of the AIA and IRIS data shows that the jets are initiated at the top of a canopy-like, double-chambered structure with cool emission on one side and hot emission in the other. The hot jets are collimated in the hot temperature filters, have high velocities (around 250 km/s) and accompanied by the cool surges and ejected kernels both moving at about 45 km/s. In the pre-phase of the jets, at their base we find quasi-periodic intensity oscillations in phase with small ejections; they have a period between 2 and 6 minutes and are reminiscent of acoustic or MHD waves. Conclusions: This series of jets and surges provides a good case-study to test the 2D and 3D magnetohydrodynamic (MHD) models that result from magnetic flux emergence. The double-chambered structure found in the observations corresponds to the cold and hot loop regions found in the models beneath the current sheet that contains the reconnection site. The cool surge with kernels is comparable with the cool ejection and plasmoids that naturally appear in the models.

Ambipolar diffusion is a physical mechanism related to the drift between charged and neutral particles in a partially ionized plasma that is key in many different astrophysical systems. However, understanding its effects is challenging due to basic uncertainties concerning relevant microphysical aspects and the strong constraints it imposes on the numerical modeling. Our aim is to introduce a numerical tool that allows us to address complex problems involving ambipolar diffusion in which, additionally, departures from ionization equilibrium are important or high resolution is needed. The primary application of this tool is for solar atmosphere calculations, but the methods and results presented here may also have a potential impact on other astrophysical systems. We have developed a new module for the stellar atmosphere Bifrost code that improves its computational capabilities of the ambipolar diffusion term in the Generalized Ohm's Law. This module includes, among other things, collision terms adequate to processes in the coolest regions in the solar chromosphere. As a key feature of the module, we have implemented the Super Time-Stepping (STS) technique, that allows an important acceleration of the calculations. We have also introduced hyperdiffusion terms to guarantee the stability of the code. We show that to have an accurate value for the ambipolar diffusion coefficient in the solar atmosphere it is necessary to include as atomic elements in the equation of state not only hydrogen and helium but also the main electron donors like sodium, silicon and potassium. In addition, we establish a range of criteria to set up an automatic selection of the free parameters of the STS method that guarantees the best performance, optimizing the stability and speed for the ambipolar diffusion calculations. We validate the STS implementation by comparison with a self-similar analytical solution.

Magnetic flux emergence has been shown to be a key mechanism for unleashing a wide variety of solar phenomena. However, there are still open questions concerning the rise of the magnetized plasma through the atmosphere, mainly in the chromosphere, where the plasma departs from local thermodynamic equilibrium (LTE) and is partially ionized. We aim to investigate the impact of the nonequilibrium (NEQ) ionization and recombination and molecule formation of hydrogen, as well as ambipolar diffusion, on the dynamics and thermodynamics of the flux emergence process. Using the Bifrost code, we performed 2.5D numerical experiments of magnetic flux emergence from the convection zone up to the corona. The experiments include the NEQ ionization and recombination of atomic hydrogen, the NEQ formation and dissociation of H2 molecules, and the ambipolar diffusion term of the Generalized Ohm's Law. Our experiments show that the LTE assumption substantially underestimates the ionization fraction in most of the emerged region, leading to an artificial increase in the ambipolar diffusion and, therefore, in the heating and temperatures as compared to those found when taking the NEQ effects on the hydrogen ion population into account. We see that LTE also overestimates the number density of H2 molecules within the emerged region, thus mistakenly magnifying the exothermic contribution of the H2 molecule formation to the thermal energy during the flux emergence process. We find that the ambipolar diffusion does not significantly affect the amount of total unsigned emerged magnetic flux, but it is important in the shocks that cross the emerged region, heating the plasma on characteristic times ranging from 0.1 to 100 s. We also briefly discuss the importance of including elements heavier than hydrogen in the equation of state so as not to overestimate the role of ambipolar diffusion in the atmosphere.

Rotating magnetic structures are common in astrophysics, from vortex tubes and tornados in the Sun all the way to jets in different astrophysical systems. The physics of these objects often combine inertial, magnetic, gas pressure and gravitational terms. Also, they often show approximate symmetries that help simplify the otherwise rather intractable equations governing their morphology and evolution. Here we propose a general formulation of the equations assuming axisymmetry and a self-similar form for all variables: in spherical coordinates $(r,\theta,\phi)$, the magnetic field and plasma velocity are taken to be of the form: ${\bf B}={\bf f}(\theta)/r^n$ and ${\bf v}={\bf g}(\theta)/r^m$, with corresponding expressions for the scalar variables like pressure and density. Solutions are obtained for potential, force-free, and non-force-free magnetic configurations. Potential-field solutions can be found for all values of~$n$. Non-potential force-free solutions possess an azimuthal component $B_\phi$ and exist only for $n\ge2$; the resulting structures are twisted and have closed field lines but are not collimated around the system axis. In the non-force free case, including gas pressure, the magnetic field lines acquire an additional curvature to compensate for an outward pointing pressure gradient force. We have also considered a pure rotation situation with no gravity, in the zero-$\beta$ limit: the solution has cylindrical geometry and twisted magnetic field lines. The latter solutions can be helpful in producing a collimated magnetic field structure; but they exist only when $n<0$ and $m<0$: for applications they must be matched to an external system at a finite distance from the origin.

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

Surges are ubiquitous cool ejections in the solar atmosphere that often appear associated with transient phenomena like UV bursts or coronal jets. Recent observations from the Interface Region Imaging Spectrograph (IRIS) show that surges, although traditionally related to chromospheric lines, can exhibit enhanced emission in Si IV with brighter spectral profiles than for the average transition region (TR). In this paper, we explain why surges are natural sites to show enhanced emissivity in TR lines. We performed 2.5D radiative-MHD numerical experiments using the Bifrost code including the nonequilibrium ionization of silicon and oxygen. A surge is obtained as a by-product of magnetic flux emergence; the TR enveloping the emerged domain is strongly affected by nonequilibrium effects: assuming statistical equilibrium would produce an absence of Si IV and O IV ions in most of the region. Studying the properties of the surge plasma emitting in the Si IV 1402.77 {\AA} and O IV 1401.16 {\AA} lines, we find that a) the timescales for the optically-thin losses and heat conduction are very short, leading to departures from statistical equilibrium, and b) the surge emits in Si IV more and has an emissivity ratio of Si IV to O IV larger than a standard TR. Using synthetic spectra, we conclude the importance of line-of-sight effects: given the involved geometry of the surge, the line-of-sight can cut the emitting layer at small angles and/or cross it multiple times, causing prominent, spatially intermittent brightenings both in Si IV and O IV.

Surges often appear as a result of the emergence of magnetized plasma from the solar interior. Traditionally, they are observed in chromospheric lines such as H$\alpha$ 6563 \AA and Ca II 8542 \AA. However, whether there is a response to the surge appearance and evolution in the Si IV lines or, in fact, in many other transition region lines has not been studied. In this paper we analyze a simultaneous episode of an H$\alpha$ surge and a Si IV burst that occurred on 2016 September 03 in active region AR12585. To that end, we use coordinated observations from the Interface Region Imaging Spectrograph (IRIS) and the Swedish 1-m Solar Telescope (SST). For the first time, we report emission of Si IV within the surge, finding profiles that are brighter and broader than the average. Furthermore, the brightest Si IV patches within the domain of the surge are located mainly near its footpoints. To understand the relation between the surges and the emission in transition region lines like Si IV, we have carried out 2.5D radiative MHD (RMHD) experiments of magnetic flux emergence episodes using the Bifrost code and including the non-equilibrium ionization of silicon. Through spectral synthesis we explain several features of the observations. We show that the presence of Si IV emission patches within the surge, their location near the surge footpoints and various observed spectral features are a natural consequence of the emergence of magnetized plasma from the interior to the atmosphere and the ensuing reconnection processes.

{We report on the first of a series of studies on coronal bright points investigating the physical mechanism that generates these phenomena.} {The aim of this paper is to understand the magnetic-field structure that hosts the bright points.} {We use longitudinal magnetograms taken by the Solar Optical Telescope with the Narrowband Filter Imager. For a single case, magnetograms from the Helioseismic and Magnetic Imager were added to the analysis. The longitudinal magnetic field component is used to derive the potential magnetic fields of the large regions around the bright points. A magneto-static field extrapolation method is tested to verify the accuracy of the potential field modelling. The three dimensional magnetic fields are investigated for the presence of magnetic null points and their influence on the local magnetic domain.} {In 9 out of 10 cases the bright point resides in areas where the coronal magnetic field contains an opposite polarity intrusion defining a magnetic null point above it. It is found that X-ray bright points reside, in these 9 cases, in a limited part of the projected fan dome area, either fully inside the dome or expanding over a limited area below which typically a dominant flux concentration resides. The 10th bright point is located in a bipolar loop system without an overlying null point.} {All bright points in coronal holes and two out of tree bright points in quiet Sun regions are seen to reside in regions containing a magnetic null point. An yet unidentified process(es) generates the BPs in specific regions of the fan-dome structure. }

Coronal jets represent important manifestations of ubiquitous solar transients, which may be the source of significant mass and energy input to the upper solar atmosphere and the solar wind. While the energy involved in a jet-like event is smaller than that of "nominal" solar flares and coronal mass ejections (CMEs), jets share many common properties with these phenomena, in particular, the explosive magnetically driven dynamics. Studies of jets could, therefore, provide critical insight for understanding the larger, more complex drivers of the solar activity. On the other side of the size-spectrum, the study of jets could also supply important clues on the physics of transients close or at the limit of the current spatial resolution such as spicules. Furthermore, jet phenomena may hint to basic process for heating the corona and accelerating the solar wind; consequently their study gives us the opportunity to attack a broad range of solar-heliospheric problems.

Cool and dense ejections, typically H$_{\alpha}$ surges, often appear alongside EUV or X-Ray coronal jets as a result of the emergence of magnetized plasma from the solar interior. Idealized numerical experiments explain those ejections as being indirectly associated with the magnetic reconnection taking place between the emerging and preexisting systems. However, those experiments miss basic elements that can importantly affect the surge phenomenon. In this paper we study the cool surges using a realistic treatment of the radiation transfer and material plasma properties. To that end, the Bifrost code is used, which has advanced modules for the equation of state of the plasma, photospheric and chromospheric radiation transfer, heat conduction and optically thin radiative cooling. We carry out a 2.5D experiment of the emergence of magnetized plasma through (meso)granular convection cells and the low atmosphere to the corona. Through detailed Lagrange tracing, we study the formation and evolution of the cool ejection and, in particular, the role of the entropy sources: this allows us to discern families of evolutionary patterns for the plasma elements. In the launch phase many elements suffer accelerations well in excess of gravity; when nearing the apex of their individual trajectories, instead, the plasma elements follow quasi-parabolic trajectories with acceleration close to g$_{\odot^{}}$. We show how the formation of the cool ejection is mediated by a wedge-like structure composed of two shocks, one of which leads to the detachment of the surge from the original emerged plasma dome.

Recent high-resolution and high-cadence observations have surprisingly suggested that prominence barbs exhibit apparent rotating motions suggestive of a tornado-like structure. Additional evidence has been provided by Doppler measurements. The observations reveal opposite velocities for both hot and cool plasma on the two sides of a prominence barb. This motion is persistent for several hours and has been interpreted in terms of rotational motion of prominence feet. Several authors suggest that such barb motions are rotating helical structures around a vertical axis similar to tornadoes on Earth. One of the difficulties of such a proposal is how to support cool prominence plasma in almost-vertical structures against gravity. In this work we model analytically a tornado-like structure and try to determine possible mechanisms to support the prominence plasma. We have found that the Lorentz force can indeed support the barb plasma provided the magnetic structure is sufficiently twisted and/or significant poloidal flows are present.

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 importance of magnetic fields in three-dimensional magnetoconvection models of the Sun's photosphere is investigated in terms of their influence on the continuum intensity at different viewing inclination angles, and on the intensity profile of two [O I] spectral lines. We use the RH numerical radiative transfer code to perform a posteriori spectral synthesis on the same time-series of magnetoconvection models used in our publications on the effect of magnetic fields on abundance determination. We obtain a good match of the synthetic disc-centre continuum intensity to the absolute continuum values from the FTS observational spectrum; the match of the centre-to-limb variation (CLV) synthetic data to observations is also good, thanks, in part, to the 3D radiation transfer capabilities of the RH code. The different levels of magnetic flux in the numerical time-series do not modify the quality of the match. Concerning the targetted [O I] spectral lines, we find, instead, that magnetic fields lead to non-negligible changes in the synthetic spectrum, with larger average magnetic flux causing the line to become noticeably weaker. The photospheric oxygen abundance that one would derive if instead using non-magnetic numerical models would thus be lower by a few to several centidexes. The inclusion of magnetic fields is confirmed to be important for improving the current modelling of the Sun, here in particular in terms of spectral line formation and of deriving consistent chemical abundances. These results may shed further light on the still controversial issue regarding the precise value of the solar oxygen abundance.

We compare Fourier spectra of photospheric velocity fields from very high resolution IMaX observations to those from recent 3D numerical magnetoconvection models. We carry out a proper comparison by synthesizing spectral lines from the numerical models and then applying to them the adequate residual instrumental degradation that affects the observational data. Also, the validity of the usual observational proxies is tested by obtaining synthetic observations from the numerical boxes and comparing the velocity proxies to the actual velocity values from the numerical grid. For the observations, data from the SUNRISE/IMaX instrument with about 120 km spatial resolution are used, thus allowing the calculation of observational Fourier spectra well into the subgranular range. For the simulations, we use four series of runs obtained with the STAGGER code and synthesize the IMaX spectral line (FeI 5250.2 A) from them. Proxies for the velocity field are obtained via Dopplergrams (vertical component) and local correlation tracking (horizontal component). A very good match between observational and simulated Fourier power spectra is obtained for the vertical velocity data for scales between 200 km and 6 Mm. Instead, a clear vertical shift is obtained when the synthetic observations are not degraded. The match for the horizontal velocity data is much less impressive because of the inaccuracies of the LCT procedure. Concerning the internal comparison of the direct velocity values of the numerical boxes with those from the synthetic observations, a high correlation (0.96) is obtained for the vertical component when using the velocity values on the log($tau_{500}$) = -1 surface in the box. The corresponding Fourier spectra are near each other. A lower maximum correlation (0.5) is reached (at $tau_{500}$ = 1) for the horizontal velocities as a result of the coarseness of the LCT procedure.

Numerical 3D radiative (M)HD simulations of solar convection are used to understand the physical properties of the solar photosphere. To validate this approach, it is important to check that no excessive thermodynamic fluctuations arise as a consequence of the partially incomplete treatment of radiative transfer. We investigate the realism of 3D convection simulations carried out with the Stagger code. We compared the characteristic properties of several spectral lines in solar disc centre observations with spectra synthesized from the simulations. We degraded the synthetic spectra to the spatial resolution of the observations using the continuum intensity distribution. We estimated the necessary spectral degradation by comparing atlas spectra with averaged observed spectra. In addition to deriving a set of line parameters directly, we used the SIR code to invert the spectra. Most of the line parameters from the observational data are matched well by the degraded simulation spectra. The inversions predict a macroturbulent velocity below 10 m/s for the simulation at full spatial resolution, whereas they yield ~< 1000 m/s at a spatial resolution of 0.3". The temperature fluctuations in the inversion of the degraded simulation do not exceed those from the observational data (of the order of 100-200 K rms for -2<log tau<-0.5). The comparison of line parameters in spatially averaged profiles with the averaged values of line parameters in spatially resolved profiles indicates a significant change of (average) line properties at a spatial scale between 0.13" and 0.3". Up to a spatial resolution of 0.3", we find no indications of the presence of excessive thermodynamic fluctuations in the 3D HD simulation. To definitely confirm that simulations without spatial degradation contain fully realistic thermodynamic fluctuations requires observations at even better spatial resolution.

A broad jet was observed in a weak magnetic field area at the edge of active region NOAA 11106. The peculiar shape and magnetic environment of the broad jet raised the question of whether it was created by the same physical processes of previously studied jets with reconnection occurring high in the corona. We carried out a multi-wavelength analysis using the EUV images from the Atmospheric Imaging Assembly (AIA) and magnetic fields from the Helioseismic and Magnetic Imager (HMI) both on-board the SDO satellite. The jet consisted of many different threads that expanded in around 10 minutes to about 100 Mm in length, with the bright features in later threads moving faster than in the early ones, reaching a maximum speed of about 200 km s^{-1}. Time-slice analysis revealed a striped pattern of dark and bright strands propagating along the jet, along with apparent damped oscillations across the jet. This is suggestive of a (un)twisting motion in the jet, possibly an Alfven wave. A topological analysis of an extrapolated field was performed. Bald patches in field lines, low-altitude flux ropes, diverging flow patterns, and a null point were identified at the basis of the jet. Unlike classical lambda or Eiffel-tower shaped jets that appear to be caused by reconnection in current sheets containing null points, reconnection in regions containing bald patches seems to be crucial in triggering the present jet. There is no observational evidence that the flux ropes detected in the topological analysis were actually being ejected themselves, as occurs in the violent phase of blowout jets; instead, the jet itself may have gained the twist of the flux rope(s) through reconnection. This event may represent a class of jets different from the classical quiescent or blowout jets, but to reach that conclusion, more observational and theoretical work is necessary.

A three-dimensional numerical experiment of the launching of a hot and fast coronal jet followed by several violent eruptions is analyzed in detail. These events are initiated through the emergence of a magnetic flux rope from the solar interior into a coronal hole. We explore the evolution of the emerging magnetically-dominated plasma dome surmounted by a current sheet and the ensuing pattern of reconnection. A hot and fast coronal jet with inverted-Y shape is produced that shows properties comparable to those frequently observed with EUV and X-Ray detectors. We analyze its 3D shape, its inhomogeneous internal structure, and its rise and decay phases, lasting for some 15-20 min each. Particular attention is devoted to the field-line connectivities and the reconnection pattern. We also study the cool and high-density volume that appears encircling the emerged dome. The decay of the jet is followed by a violent phase with a total of five eruptions. The first of them seems to follow the general pattern of tether-cutting reconnection in a sheared arcade, although modified by the field topology created by the preceding reconnection evolution. The two following eruptions take place near and above the strong field-concentrations at the surface. They show a twisted, \Omega-loop like rope expanding in height, with twist being turned into writhe, thus hinting at a kink instability (perhaps combined with a torus-instability) as the cause of the eruption. The succession of a main jet ejection and a number of violent eruptions that resemble mini-CME's and their physical properties suggest that this experiment may provide a model for the blowout jets recently proposed in the literature.

We investigate the impact on Fe abundance determination of including magnetic flux in series of 3D radiation-MHD simulations of solar convection which we used to synthesize spectral intensity profiles corresponding to disc centre. A differential approach is used to quantify the changes in theoretical equivalent width of a set of 28 iron spectral lines spanning a wide range in lambda, excitation potential, oscillator strength, Land\'e factor, and formation height. The lines were computed in LTE using the spectral synthesis code LILIA. We used input magnetoconvection snapshots covering 50 minutes of solar evolution and belonging to series having an average vertical magnetic flux density of 0, 50, 100 and 200 G. For the relevant calculations we used the Copenhagen Stagger code. The presence of magnetic fields causes both a direct (Zeeman-broadening) effect on spectral lines with non-zero Land\'e factor and an indirect effect on temperature-sensitive lines via a change in the photospheric T-tau stratification. The corresponding correction in the estimated atomic abundance ranges from a few hundredths of a dex up to |Delta log(Fe)| ~ 0.15 dex, depending on the spectral line and on the amount of average magnetic flux within the range of values we considered. The Zeeman-broadening effect gains relatively more importance in the IR. The largest modification to previous solar abundance determinations based on visible spectral lines is instead due to the indirect effect, i.e., the line-weakening caused by a warmer stratification on an optical depth scale. Our results indicate that the average solar iron abundance obtained when using magnetoconvection models can be 0.03-0.11 dex higher than when using the simpler HD convection approach. We demonstrate that accounting for magnetic flux is important in state-of-the-art solar photospheric abundance determinations based on 3D simulations.

The 3D structure of an active region (AR) filament is studied using nonlinear force-free field (NLFFF) extrapolations based on simultaneous observations at a photospheric and a chromospheric height. To that end, we used the Si I 10827 \AA\ line and the He I 10830 \AA\ triplet obtained with the Tenerife Infrared Polarimeter (TIP) at the VTT (Tenerife). The two extrapolations have been carried out independently from each other and their respective spatial domains overlap in a considerable height range. This opens up new possibilities for diagnostics in addition to the usual ones obtained through a single extrapolation from, typically, a photospheric layer. Among those possibilities, this method allows the determination of an average formation height of the He I 10830 \AA\ signal of \approx 2 Mm above the surface of the sun. It allows, as well, to cross-check the obtained 3D magnetic structures in view of verifying a possible deviation from the force- free condition especially at the photosphere. The extrapolations yield a filament formed by a twisted flux rope whose axis is located at about 1.4 Mm above the solar surface. The twisted field lines make slightly more than one turn along the filament within our box, which results in 0.055 turns/Mm. The convex part of the field lines (as seen from the solar surface) constitute dips where the plasma can naturally be supported. The obtained 3D magnetic structure of the filament depends on the choice of the observed horizontal magnetic field as determined from the 180\circ solution of the azimuth. We derive a method to check for the correctness of the selected 180\circ ambiguity solution.

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.

The relation of the solar surface magnetic field with mesogranular cells is studied using high spatial (~ 100 km) and temporal (~ 30 sec) resolution data obtained with the IMaX instrument aboard SUNRISE. First, mesogranular cells are identified using Lagrange tracers (corks) based on horizontal velocity fields obtained through Local Correlation Tracking. After ~ 20 min of integration, the tracers delineate a sharp mesogranular network with lanes of width below about 280 km. The preferential location of magnetic elements in mesogranular cells is tested quantitatively. Roughly 85% of pixels with magnetic field higher than 100 G are located in the near neighborhood of mesogranular lanes. Magnetic flux is therefore concentrated in mesogranular lanes rather than intergranular ones. Secondly, magnetic field extrapolations are performed to obtain field lines anchored in the observed flux elements. This analysis, therefore, is independent of the horizontal flows determined in the first part. A probability density function (PDF) is calculated for the distribution of distances between the footpoints of individual magnetic field lines. The PDF has an exponential shape at scales between 1 and 10 Mm, with a constant characteristic decay distance, indicating the absence of preferred convection scales in the mesogranular range. Our results support the view that mesogranulation is not an intrinsic convective scale (in the sense that it is not a primary energy-injection scale of solar convection), but also give quantitative confirmation that, nevertheless, the magnetic elements are preferentially found along mesogranular lanes.

We explore the effect of the magnetic field when using realistic three-dimensional convection experiments to determine solar element abundances. By carrying out magnetoconvection simulations with a radiation-hydro code (the Copenhagen stagger code) and through a-posteriori spectral synthesis of three Fe I lines, we obtain evidence that moderate amounts of mean magnetic flux cause a noticeable change in the derived equivalent widths compared with those for a non-magnetic case. The corresponding Fe abundance correction for a mean flux density of 200 G reaches up to ~0.1 dex in magnitude. These results are based on space- and time-averaged line profiles over a time span of 2.5 solar hours in the statistically stationary regime of the convection. The main factors causing the change in equivalent widths, namely the Zeeman broadening and the modification of the temperature stratification, act in different amounts and, for the iron lines considered here, in opposite directions; yet, the resulting $|\Delta\log\epsilon_{\odot}(Fe)|$ coincides within a factor two in all of them, even though the sign of the total abundance correction is different for the visible and infrared lines. We conclude that magnetic effects should be taken into account when discussing precise values of the solar and stellar abundances and that an extended study is warranted.

Recent observations of coronal hole areas with the XRT and EIS instruments onboard the Hinode satellite have shown with unprecedented detail the launching of fast, hot jets away from the solar surface. In some cases these events coincide with episodes of flux emergence from beneath the photosphere. In this letter we show results of a 3D numerical experiment of flux emergence from the solar interior into a coronal hole and compare them with simultaneous XRT and EIS observations of a jet-launching event that accompanied the appearance of a bipolar region in MDI magnetograms. The magnetic skeleton and topology that result in the experiment bear a strong resemblance to linear force-fee extrapolations of the SOHO/MDI magnetograms. A thin current sheet is formed at the boundary of the emerging plasma. A jet is launched upward along the open reconnected field lines with values of temperature, density and velocity in agreement with the XRT and EIS observations. Below the jet, a split-vault structure results with two chambers: a shrinking one containing the emerged field loops and a growing one with loops produced by the reconnection. The ongoing reconnection leads to a horizontal drift of the vault-and-jet structure. The timescales, velocities, and other plasma properties in the experiment are consistent with recent statistical studies of this type of events made with Hinode data.

The emergence of magnetic flux from the convection zone into the corona is an important process for the dynamical evolution of the coronal magnetic field. In this paper we extend our previous numerical investigations, by looking at the process of flux interaction as an initially twisted flux tube emerges into a plane parallel, coronal magnetic field. Significant differences are found in the dynamical appearance and evolution of the emergence process depending on the relative orientation between the rising flux system and any preexisting coronal field. When the flux systems are nearly anti-parallel, the experiments show substantial reconnection and demonstrate clear signatures of a high temperature plasma located in the high velocity outflow regions extending from the reconnection region. However, the cases that have a more parallel orientation of the flux systems show very limited reconnection and none of the associated features. Despite the very different amount of reconnection between the two flux systems, it is found that the emerging flux that is still connected to the original tube, reaches the same height as a function of time. As a compensation for the loss of tube flux, a clear difference is found in the extent of the emerging loop in the direction perpendicular to the main axis of the initial flux tube. Increasing amounts of magnetic reconnection decrease the volume, which confines the remaining tube flux.

We study the emergence of magnetic flux from the near-surface layers of the solar convection zone into the photosphere. To model magnetic flux emergence, we carried out a set of numerical radiative magnetohydrodynamics simulations. Our simulations take into account the effects of compressibility, energy exchange via radiative transfer, and partial ionization in the equation of state. All these physical ingredients are essential for a proper treatment of the problem. Furthermore, the inclusion of radiative transfer allows us to directly compare the simulation results with actual observations of emerging flux. We find that the interaction between the magnetic flux tube and the external flow field has an important influence on the emergent morphology of the magnetic field. Depending on the initial properties of the flux tube (e.g. field strength, twist, entropy etc.), the emergence process can also modify the local granulation pattern. The emergence of magnetic flux tubes with a flux of $10^{19}$ Mx disturbs the granulation and leads to the transient appearance of a dark lane, which is coincident with upflowing material. These results are consistent with observed properties of emerging magnetic flux.

We study the structure and reveal the physical nature of the reversed granulation pattern in the solar photosphere by means of 3-dimensional radiative hydrodynamics simulations. We used the MURaM code to obtain a realistic model of the near-surface layers of the convection zone and the photosphere. The pattern of horizontal temperature fluctuations at the base of the photosphere consists of relatively hot granular cells bounded by the cooler intergranular downflow network. With increasing height in the photosphere, the amplitude of the temperature fluctuations diminishes. At a height of z=130-140 km in the photosphere, the pattern of horizontal temperature fluctuations reverses so that granular regions become relatively cool compared to the intergranular network. Detailed analysis of the trajectories of fluid elements through the photosphere reveal that the motion of the fluid is non-adiabatic, owing to strong radiative cooling when approaching the surface of optical depth unity followed by reheating by the radiation field from below. The temperature structure of the photosphere results from the competition between expansion of rising fluid elements and radiative heating. The former acts to lower the temperature of the fluid whereas the latter acts to increase it towards the radiative equilibrium temperature with a net entropy gain. After the fluid overturns and descends towards the convection zone, radiative energy loss again decreases the entropy of the fluid. Radiative heating and cooling of fluid elements that penetrate into the photosphere and overturn do not occur in equal amounts. The imbalance in the cumulative heating and cooling of these fluid elements is responsible for the reversal of temperature fluctuations with respect to height in the photosphere.

We present the results of a set of three-dimensional numerical simulations of magnetic flux emergence from below the photosphere into the corona that include a uniform and horizontal coronal magnetic field mimicking a pre-existing large-scale coronal magnetic system. Cases with different relative orientations of the upcoming and coronal fields are studied. Upon contact, a concentrated current sheet with the shape of an arch or bridge is formed at the interface which marks the positions of maximum jump in the field vector between the two systems. Relative angles above 90 degrees yield abundant magnetic reconnection and plasma heating. The reconnection is seen to be intrisincally three-dimensional in nature, except at singular positions along the current sheet. It drives collimated high-speed and high-temperature outflows only a short distance from the reconnection site that propagate along the ambient magnetic field lines as jets. Due to the low plasma density in the corona, these jets may propagate over large distances and, therefore help distribute high-density and high-temperature plasma along these newly reconnected field lines. The experiments permit to discern and visualize the three-dimensional shape and relative position of the upcoming plasma hill, high-speed jets and subphotospheric flux system. As a result of the reconnection, magnetic field lines from the magnetized plasma below the surface end up as coronal field lines, thus causing a profound change in the connectivity of the magnetic regions in the corona. The experiments presented here thus yield a number of features repeatedly observed with the TRACE satellite and the YOHKOH-SXT detector, like the establishment of connectivity between emergent and pre-existing active regions, local heating and high-velocity outflows.

A systematic study of the linear thermal stability of a medium subject to cooling, self-gravity and thermal conduction is carried out for the case when the unperturbed state is subject to global cooling and expansion. A general, recursive WKB solution for the perturbation problem is obtained which can be applied to a large variety of situations in which there is a separation of time-scales for the different physical processes. Solutions are explicitly given and discussed for the case when sound propagation and/or self-gravity are the fastest processes, with cooling, expansion and thermal conduction operating on slower time-scales. A brief discussion is also added for the solutions in the cases in which cooling or conduction operate on the fastest time-scale. The general WKB solution obtained in this paper permits solving the problem of the effect of thermal conduction and self-gravity on the thermal stability of a globally cooling and expanding medium. As a result of the analysis, the critical wavelength (often called Field length) above which cooling makes the perturbations unstable against the action of thermal conduction is generalized to the case of an unperturbed background with net cooling. As an astrophysical application, the generalized Field length is calculated for a hot (10^4 - 10^8 K), optically thin medium (as pertains, for instance, for the hot interstellar medium of SNRs or superbubbles) using a realistic cooling function and including a weak magnetic field. The stability domains are compared with the predictions made on the basis of models for which the background is in thermal equilibrium. The instability domain of the sound waves, in particular, is seen to be much larger in the case with net global cooling.

Results from a numerical and analytical investigation of the solution of a nonlinear axially symmetric diffusion equation for the magnetic field are presented for the case when the nonlinear dependence of the diffusivity nu(B) on the magnetic field satisfies basic physical requirements. We find that for sufficiently strong nonlinearity (i.e. for sufficiently strong reduction of nu inside the tube) a current sheet is spontaneously formed around the tube within one diffusion timescale. This sheet propagates inwards with a velocity inversely proportional to the ratio of the field strength just inside the current sheet to the equipartition field strength B0/Be, so the lifetime of a tube with constant internal flux density is increased approximately by a factor not exceeding B0/Be, even for infinitely effective inhibition of turbulence inside the tube. Among the applications of these results we point out that toroidal flux tubes in the solar convective zone are subject to significant flux loss owing to turbulent erosion on a timescale of about 1 month, and that turbulent erosion may be responsible for the formation of a current sheet around a sunspot. It is further proposed that, despite the simplifying assumptions involved, our solutions correctly reflect the essential features of the sunspot decay process.

The physics of a twisted magnetic flux tube rising in a stratified medium is studied using a numerical MHD code. The problem considered is fully compressible (no Boussinesq approximation), includes ohmic resistivity, and is two dimensional, i.e., there is no variation of the variables in the direction of the tube axis. We study a high plasma beta case with small ratio of radius to external pressure scaleheight. The results obtained can therefore be of relevance to understand the transport of magnetic flux across the solar convection zone.