Simulations Of Magnetic Flux Emergence Research Articles

Turbulent convection as a significant hidden provider of magnetic helicity in solar eruptions

Solar flares and coronal mass ejections, the primary space weather disturbances affecting the entire heliosphere and near-Earth environment, mainly emanate from sunspot regions harbouring high degrees of magnetic twist. However, it is not clear how magnetic helicity, the quantity for measuring the magnetic twist, is supplied to the upper solar atmosphere via the emergence of magnetic flux from the turbulent convection zone. Here, we report state-of-the-art numerical simulations of magnetic flux emergence from the deep convection zone. By controlling the twist of emerging flux, we find that with the support of convective upflow, the untwisted emerging flux can reach the solar surface without collapsing, in contrast to previous theoretical predictions, and eventually create sunspots. Because of the turbulent twisting of magnetic flux, the produced sunspots exhibit rotation and inject magnetic helicity into the upper atmosphere, amounting to a substantial fraction of injected helicity in the twisted cases that is sufficient to produce flare eruptions. This result indicates that the turbulent convection is responsible for supplying a non-negligible amount of magnetic helicity and potentially contributes to solar flares.

Forward modelling of brightness variations in Sun-like stars

Context. The amplitude and morphology of light curves of Sun-like stars change substantially with increasing rotation rate: brightness variations are amplified and become more regular. This has not been explained so far. Aims. We develop a modelling approach for calculating brightness variations of stars with various rotation rates and use it to explain the observed trends in stellar photometric variability. Methods. We combined numerical simulations of magnetic flux emergence and transport with a model for stellar brightness variability to calculate synthetic light curves of stars as observed by the Kepler telescope. We computed the distribution of the magnetic flux on the stellar surface for various rotation rates and degrees of active-region nesting (i.e. the tendency of active regions to emerge in the vicinity of recently emerged regions). Using the resulting maps of the magnetic flux, we computed the rotational variability of our simulated stellar light curves as a function of rotation rate and nesting of magnetic features and compared our calculations to Kepler observations. Results. We show that both the rotation rate and the degree of nesting have a strong impact on the amplitude and morphology of stellar light curves. In order to explain the variability of most of the Kepler targets with known rotation rates, we need to increase the degree of nesting to values that are much higher than the values on the Sun. Conclusions. The suggested increase in nesting with the rotation rate can provide clues about the flux emergence process for high levels of stellar activity.

Association Between Magnetic Pressure Difference and the Movement of Solar Pores

Solar pores are closely related to the concentration, dissipation, and transportation of solar magnetic flux. Their observable characteristics can provide constraints on models and simulations of magnetic flux emergence and formation. The specific property investigated in this study is their horizontal movement. The aim is to investigate whether the movement is correlated with any observable quantities. Our statistical analysis of 61 compact pores identified from the Spaceweather HMI Active Region Patches (SHARP) from 2011 to 2018 indicates that the direction of movement is often either parallel or antiparallel to the direction of maximum magnetic pressure difference (dP mag) at the opposite sides of the edges of the pores. The correlation coefficients for both the parallel and antiparallel cases are higher than 0.74. Despite the high correlation, our analysis using the transfer entropy indicates no significant causal relationship between the direction of motion and the direction of maximum dP mag.

Ellerman bombs and UV bursts: reconnection at different atmospheric layers

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

Spontaneous Generation of δ-sunspots in Convective Magnetohydrodynamic Simulation of Magnetic Flux Emergence

Abstract Observations reveal that strong solar flares and coronal mass ejections tend to occur in complex active regions characterized by δ-sunspots, spot rotation, sheared polarity inversion lines (PILs), and magnetic flux ropes. Here we report on the first modeling of spontaneous δ-spot generation as a result of flux emergence from the turbulent convection zone. Utilizing state-of-the-art radiative magnetohydrodynamics code R2D2, we simulate the emergence of a force-free flux tube in the convection zone that stretches down to −140 Mm. Elevated by large-scale convective upflows, the tube appears on the photosphere as two emerging bipoles. The opposite polarities collide against each other due to the subsurface connectivity, and they develop into a pair of closely packed δ-spots. The Lorentz force drives the spot rotation and a strong counter-streaming flow of 10 km s−1 at the PIL in δ-spots, which, in tandem with local convection, strengthens the horizontal field to 4 kG and builds up a highly sheared PIL. In the atmosphere above the PIL, a flux rope structure is created. All these processes follow the multi-buoyant segment theory of the δ-spot formation, and they occur as a natural consequence of interaction between magnetic flux and turbulent convection, suggesting that the generation of δ-spots and the resultant flare eruptions may be a stochastically determined process.

Forward modelling of brightness variations in Sun-like stars

Context. The latitudinal distribution of starspots deviates from the solar pattern with increasing rotation rate. Numerical simulations of magnetic flux emergence and transport can help model the observed stellar activity patterns and the associated brightness variations. Aims. We set up a composite model for the processes of flux emergence and transport on Sun-like stars to simulate stellar brightness variations for various levels of magnetic activity and rotation rates. Methods. Assuming that the distribution of magnetic flux at the base of the convection zone follows solar scaling relations, we calculate the emergence latitudes and tilt angles of bipolar regions at the surface for various rotation rates, using thin-flux-tube simulations. Taking these two quantities as input to a surface flux transport (SFT) model, we simulate the diffusive-advective evolution of the radial field at the stellar surface, including effects of active region nesting. Results. As the rotation rate increases, (1) magnetic flux emerges at higher latitudes and an inactive gap opens around the equator, reaching a half-width of 20° for 8 Ω⊙; and (2) the tilt angles of freshly emerged bipolar regions show stronger variations with latitude. Polar spots can form at 8 Ω⊙ by accumulation of follower-polarity flux from decaying bipolar regions. From 4 Ω⊙ to 8 Ω⊙, the maximum spot coverage changes from 3 to 20%, respectively, compared to 0.4% in the solar model. Nesting of activity can lead to strongly non-axisymmetric spot distributions. Conclusions. On Sun-like stars rotating at 8 Ω⊙ (Prot ≃ 3 days), polar spots can form, owing to higher levels of flux emergence rate and tilt angles. Defining spots by a threshold field strength yields global spot coverages that are roughly consistent with stellar observations.

Photospheric Velocity Structures during the Emergence of Small Active Regions on the Sun

Abstract We study the plasma flows in the solar photosphere during the emergence of two small active regions, NOAA 9021 and 10768. Using Solar and Heliospheric Observatory/Michelson Doppler Imager data, we find that the strong plasma upflows appear at the initial stage of active region formation, with maximum upflow velocities of −1650 and −1320 m s−1. The structures with enhanced upflows have size ∼8 Mm in diameter, and they exist for 1–2 hr. The parameters of the enhanced upflows are consistent with those of the large active region NOAA 10488, which may suggest the possibility that the elementary emerging magnetic loops that appear at the earliest phase of active region formation have similar properties, irrespective of scales of active regions. Comparison between the observations and a numerical simulation of magnetic flux emergence shows a striking consistency. We find that the driving force of the plasma upflow is at first the gas pressure gradient and later the magnetic pressure gradient.

Bombs and Flares at the Surface and Lower Atmosphere of the Sun

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

THE ORIGIN OF NET ELECTRIC CURRENTS IN SOLAR ACTIVE REGIONS

There is a recurring question in solar physics about whether or not electric currents are neutralized in active regions (ARs). This question was recently revisited using three-dimensional (3D) magnetohydrodynamic (MHD) numerical simulations of magnetic flux emergence into the solar atmosphere. Such simulations showed that flux emergence can generate a substantial net current in ARs. Another source of AR currents are photospheric horizontal flows. Our aim is to determine the conditions for the occurrence of net vs. neutralized currents with this second mechanism. Using 3D MHD simulations, we systematically impose line-tied, quasi-static, photospheric twisting and shearing motions to a bipolar potential magnetic field. We find that such flows: (1) produce both {\it direct} and {\it return} currents, (2) induce very weak compression currents - not observed in 2.5D - in the ambient field present in the close vicinity of the current-carrying field, and (3) can generate force-free magnetic fields with a net current. We demonstrate that neutralized currents are in general produced only in the absence of magnetic shear at the photospheric polarity inversion line - a special condition rarely observed. We conclude that, as magnetic flux emergence, photospheric flows can build up net currents in the solar atmosphere, in agreement with recent observations. These results thus provide support for eruption models based on pre-eruption magnetic fields possessing a net coronal current.

δ-SUNSPOT FORMATION IN SIMULATION OF ACTIVE-REGION-SCALE FLUX EMERGENCE

$\delta$-sunspots, with highly complex magnetic structures, are very productive in energetic eruptive events, such as X-class flares and homologous eruptions. We here study the formation of such complex magnetic structures by numerical simulations of magnetic flux emergence from the convection zone into the corona in an active-region-scale domain. In our simulation, two pairs of bipolar sunspots form on the surface, originating from two buoyant segments of a single subsurface twisted flux rope, following the approach of Toriumi et al. (2014). Expansion and rotation of the emerging fields in the two bipoles drive the two opposite polarities into each other with apparent rotating motion, producing a compact $\delta$-sunspot with a sharp polarity inversion line. The formation of the $\delta$-sunspot in such a realistic-scale domain produces emerging patterns similar to those formed in observations, e.g. the inverted polarity against Hale's law, the curvilinear motion of the spot, strong transverse field with highly sheared magnetic and velocity fields at the PIL. Strong current builds up at the PIL, giving rise to reconnection, which produces a complex coronal magnetic connectivity with non-potential fields in the Delta-spot overlaid by more relaxed fields connecting the two polarities at the two ends.

Observations and 3D MHD simulations of a solar active region jet

Aims. We study an active region jet originating from NOAA 8531 on May 15 1999. We perform 3D MHD numerical simulations of magnetic flux emergence and its subsequent reconnection with preexisting magnetic flux. Then, we compare the physical properties of the observed jet with the reconnecting outflow produced in the numerical model. Methods. We report observations of this jet using a series of TRACE 171 A filtergrams, simultaneous observations from SUMER in Ne viii 770A and C iv 1548 A as well as MDI magnetograms. In the numerical simulation, the full compressible and resistive MHD equations are solved, including viscous and Ohmic heating. Results. A high-velocity upflow ($\simeq$100 km s -1 ) is observed after the emergence of new magnetic flux at the edge of the active region. The jet is recorded over a range of temperatures between 10 5 K and $1.5 \times 10^6$ K. In our numerical experiments, we find that the jet is the result of magnetic reconnection between newly emerging flux and the preexisting magnetic field of the active region. Conclusions. The hot and high-velocity bidirectional flows occur as a result of the interaction between oppositely directed magnetic fields. Observations and numerical results are strongly suggestive of effective reconnection process being responsible for producing jets when emerging flux appears in solar active regions.

TWISTED FLUX TUBE EMERGENCE FROM THE CONVECTION ZONE TO THE CORONA. II. LATER STATES

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.

A Three-dimensional Study of Reconnection, Current Sheets, and Jets Resulting from Magnetic Flux Emergence in the Sun

We present the results of a set of three-dimensional numerical simulations of magnetic flux emergence from below the photosphere and into the corona. The corona includes a uniform and horizontal magnetic field as a model for a preexisting 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 is formed at the interface that marks the positions of maximum jump in the field vector between the two systems. Relative angles above 90° yield abundant magnetic reconnection and plasma heating. The reconnection is seen to be intrinsically three-dimensional in nature and to be accompanied by marked local heating. It generates collimated high-speed outflows only a short distance from the reconnection site, and these propagate along the ambient magnetic field lines as jets. As a result of the reconnection, magnetic field lines from the magnetized plasma below the surface end up connecting to coronal field lines, thus causing a profound change in the connectivity of the magnetic regions in the corona. The experiments presented here yield a number of features repeatedly observed with the TRACE and Yohkoh satellites, such as the establishment of connectivity between emergent and preexisting active regions, local heating, and high-velocity outflows.