The physics of twisted magnetic tubes rising in a stratified medium: two dimensional results

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.Comment: To be published in ApJ, Vol. 492, Jan 10th, 1998; 25 pages, 16 figures. NEW VERSION: THE PREVIOUS ONE DIDN'T PRINT CORRECTLY. The style file overrulehere.sty is include

Plasma Jets and Eruptions in Solar Coronal Holes: a 3D flux emergence experiment

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.Comment: Accepted for publication in The Astrophysical Journal (vol 770, June 2013

Are tornado-like magnetic structures able to support solar prominence plasma?

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.Comment: Accepted for publication in ApJ

The power spectrum of solar convection flows from high-resolution observations and 3D simulations

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.Comment: 12 pages, 9 figures, accepted in A&

Turbulent Erosion of Magnetic Flux Tubes

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.Comment: 17 pages, 11 figure

Solar Fe abundance and magnetic fields - Towards a consistent reference metallicity

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.Comment: 12 pages, 7 figures, A&A in pres

Numerical Experiments on the Two-step Emergence of Twisted Magnetic Flux Tubes in the Sun

We present the new results of the two-dimensional numerical experiments on the cross-sectional evolution of a twisted magnetic flux tube rising from the deeper solar convection zone (-20,000 km) to the corona through the surface. The initial depth is ten times deeper than most of previous calculations focusing on the flux emergence from the uppermost convection zone. We find that the evolution is illustrated by the two-step process described below: the initial tube rises due to its buoyancy, subject to aerodynamic drag due to the external flow. Because of the azimuthal component of the magnetic field, the tube maintains its coherency and does not deform to become a vortex roll pair. When the flux tube approaches the photosphere and expands sufficiently, the plasma on the rising tube accumulates to suppress the tube's emergence. Therefore, the flux decelerates and extends horizontally beneath the surface. This new finding owes to our large scale simulation calculating simultaneously the dynamics within the interior as well as above the surface. As the magnetic pressure gradient increases around the surface, magnetic buoyancy instability is triggered locally and, as a result, the flux rises further into the solar corona. We also find that the deceleration occurs at a higher altitude than in our previous experiment using magnetic flux sheets (Toriumi and Yokoyama). By conducting parametric studies, we investigate the conditions for the two-step emergence of the rising flux tube: field strength > 1.5x10^4 G and the twist > 5.0x10^-4 km^-1 at -20,000 km depth.Comment: 42 pages, 13 figures, 2 tables, accepted for publication in ApJ. High-resolution figures will appear in the published versio