Data-driven modeling of solar coronal magnetic field evolution and eruptions.

Magnetic fields play a fundamental role in the structure and dynamics of the solar corona. As they are driven by their footpoint motions on the solar surface, which transport energy from the interior of the Sun into its atmosphere, the coronal magnetic fields are stressed continuously with buildup of magnetic nonpotentiality in the form of topology complexity (magnetic helicity) and local electric currents (magnetic free energy). The accumulated nonpotentiality is often released explosively by solar eruptions, manifested as solar flares and coronal mass ejections, during which magnetic energy is converted into mainly kinetic, thermal, and nonthermal energy of the plasma, which can cause adverse space weather. To reveal the physical mechanisms underlying solar eruptions, it is vital to know the three-dimensional (3D) structure and evolution of the coronal magnetic fields. Because of a lack of direct measurements, the 3D coronal magnetic fields are commonly studied using numerical modeling, whereas traditional models mostly aim for a static extrapolation of the coronal field from the observable photospheric magnetic field data. Over the last decade, dynamic models that are driven directly by observation magnetograms have been developed and applied successfully to study solar coronal magnetic field evolution as well as its eruption, which offers a novel avenue for understanding their underlying magnetic topology and mechanism. In this paper, we review the basic methodology of the data-driven coronal models, state-of-the-art developments, their typical applications, and new physics that have been derived using these models. Finally, we provide an outlook for future developments and applications of the data-driven models.

On a limitation of Zeeman polarimetry and imperfect instrumentation in representing solar magnetic fields with weaker polarization signal.

Full disk vector magnetic fields are used widely for developing better understanding of large-scale structure, morphology, and patterns of the solar magnetic field. The data are also important for modeling various solar phenomena. However, observations of vector magnetic fields have one important limitation that may affect the determination of the true magnetic field orientation. This limitation stems from our ability to interpret the differing character of the Zeeman polarization signals which arise from the photospheric line-of-sight vs. the transverse components of the solar vector magnetic field, and is likely exacerbated by unresolved structure (non-unity fill fraction) as well as the disambiguation of the 180° degeneracy in the transverse-field azimuth. Here we provide a description of this phenomenon, and discuss issues, which require additional investigation.

NuSTAR observations of a repeatedly microflaring active region.

We investigate the spatial, temporal, and spectral properties of 10 microflares from AR12721 on 2018 September 9 and 10 observed in X-rays using the Nuclear Spectroscopic Telescope ARray and the Solar Dynamic Observatory's Atmospheric Imaging Assembly and Helioseismic and Magnetic Imager. We find GOES sub-A class equivalent microflare energies of 1026-1028 erg reaching temperatures up to 10 MK with consistent quiescent or hot active region (AR) core plasma temperatures of 3-4 MK. One microflare (SOL2018-09-09T10:33), with an equivalent GOES class of A0.1, has non-thermal hard X-ray emission during its impulsive phase (of non-thermal power ~7 × 1024 erg s-1) making it one of the faintest X-ray microflares to have direct evidence for accelerated electrons. In 4 of the 10 microflares, we find that the X-ray time profile matches fainter and more transient sources in the extreme-ultraviolet, highlighting the need for observations sensitive to only the hottest material that reaches temperatures higher than those of the AR core (>5 MK). Evidence for corresponding photospheric magnetic flux cancellation/emergence present at the footpoints of eight microflares is also observed.

Magnetoacoustic wave energy dissipation in the atmosphere of solar pores.

The suitability of solar pores as magnetic wave guides has been a key topic of discussion in recent years. Here, we present observational evidence of propagating magnetohydrodynamic wave activity in a group of five photospheric solar pores. Employing data obtained by the Facility Infrared Spectropolarimeter at the Dunn Solar Telescope, oscillations with periods of the order of 5 min were detected at varying atmospheric heights by examining Si ɪ 10827 Å line bisector velocities. Spectropolarimetric inversions, coupled with the spatially resolved root mean square bisector velocities, allowed the wave energy fluxes to be estimated as a function of atmospheric height for each pore. We find propagating magnetoacoustic sausage mode waves with energy fluxes on the order of 30 kW m-2 at an atmospheric height of 100 km, dropping to approximately 2 kW m-2 at an atmospheric height of around 500 km. The cross-sectional structuring of the energy fluxes reveals the presence of both body- and surface-mode sausage waves. Examination of the energy flux decay with atmospheric height provides an estimate of the damping length, found to have an average value across all five pores of Ld ≈ 268 km, similar to the photospheric density scale height. We find the damping lengths are longer for body mode waves, suggesting that surface mode sausage oscillations are able to more readily dissipate their embedded wave energies. This work verifies the suitability of solar pores to act as efficient conduits when guiding magnetoacoustic wave energy upwards into the outer solar atmosphere. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

High-resolution wave dynamics in the lower solar atmosphere.

The magnetic and convective nature of the Sun's photosphere provides a unique platform from which generated waves can be modelled, observed and interpreted across a wide breadth of spatial and temporal scales. As oscillations are generated in-situ or emerge through the photospheric layers, the interplay between the rapidly evolving densities, temperatures and magnetic field strengths provides dynamic evolution of the embedded wave modes as they propagate into the tenuous solar chromosphere. A focused science team was assembled to discuss the current challenges faced in wave studies in the lower solar atmosphere, including those related to spectropolarimetry and radiative transfer in the optically thick regions. Following the Theo Murphy international scientific meeting held at Chicheley Hall during February 2020, the scientific team worked collaboratively to produce 15 independent publications for the current Special Issue, which are introduced here. Implications from the current research efforts are discussed in terms of upcoming next-generation observing and high-performance computing facilities. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

Placing limits on long-term variations in quiet-Sun irradiance and their contribution to total solar irradiance and solar radiative forcing of climate.

Recent reconstructions of total solar irradiance (TSI) postulate that quiet-Sun variations could give significant changes to the solar power input to Earth's climate (radiative climate forcings of 0.7-1.1 W m-2 over 1700-2019) arising from changes in quiet-Sun magnetic fields that have not, as yet, been observed. Reconstructions without such changes yield solar forcings that are smaller by a factor of more than 10. We study the quiet-Sun TSI since 1995 for three reasons: (i) this interval shows rapid decay in average solar activity following the grand solar maximum in 1985 (such that activity in 2019 was broadly equivalent to that in 1900); (ii) there is improved consensus between TSI observations; and (iii) it contains the first modelling of TSI that is independent of the observations. Our analysis shows that the most likely upward drift in quiet-Sun radiative forcing since 1700 is between +0.07 and -0.13 W m-2. Hence, we cannot yet discriminate between the quiet-Sun TSI being enhanced or reduced during the Maunder and Dalton sunspot minima, although there is a growing consensus from the combinations of models and observations that it was slightly enhanced. We present reconstructions that add quiet-Sun TSI and its uncertainty to models that reconstruct the effects of sunspots and faculae.

Investigation of the Middle Corona with SWAP and a Data-Driven Non-Potential Coronal Magnetic Field Model.

The large field-of-view of the Sun Watcher using Active Pixel System detector and Image Processing (SWAP) instrument onboard the PRoject for Onboard Autonomy 2 (PROBA2) spacecraft provides a unique opportunity to study extended coronal structures observed in the EUV in conjunction with global coronal magnetic field simulations. A global non-potential magnetic field model is used to simulate the evolution of the global corona from 1 September 2014 to 31 March 2015, driven by newly emerging bipolar active regions determined from Helioseismic and Magnetic Imager (HMI) magnetograms. We compare the large-scale structure of the simulated magnetic field with structures seen off-limb in SWAP EUV observations. In particular, we investigate how successful the model is in reproducing regions of closed and open structures, the scale of structures, and compare the evolution of a coronal fan observed over several rotations. The model is found to accurately reproduce observed large-scale, off-limb structures. When discrepancies do arise they mainly occur off the east solar limb due to active regions emerging on the far side of the Sun, which cannot be incorporated into the model until they are observed on the Earth-facing side. When such "late" active region emergences are incorporated into the model, we find that the simulated corona self-corrects within a few days, so that simulated structures off the west limb more closely match what is observed. Where the model is less successful, we consider how this may be addressed, through model developments or additional observational products. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s11207-020-01668-2) contains supplementary material, which is available to authorized users.

Average motion of emerging solar active region polarities: II. Joy's law.

CONTEXT: The tilt of solar active regions described by Joy's law is essential for converting a toroidal field to a poloidal field in Babcock-Leighton dynamo models. In thin flux tube models the Coriolis force causes what we observe as Joy's law, acting on east-west flows as they rise towards the surface. AIMS: Our goal is to measure the evolution of the average tilt angle of hundreds of active regions as they emerge, so that we can constrain the origins of Joy's law. METHODS: We measured the tilt angle of the primary bipoles in 153 emerging active regions (EARs) in the Solar Dynamics Observatory Helioseismic Emerging Active Region survey. We used line-of-sight magnetic field measurements averaged over 6 h to define the polarities and measure the tilt angle up to four days after emergence. RESULTS: We find that at the time of emergence the polarities are on average aligned east-west, and that neither the separation nor the tilt depends on latitude. We do find, however, that EARs at higher latitudes have a faster north-south separation speed than those closer to the equator at the emergence time. After emergence, the tilt angle increases and Joy's law is evident about two days later. The scatter in the tilt angle is independent of flux until about one day after emergence, when we find that higher-flux regions have a smaller scatter in tilt angle than lower-flux regions. CONCLUSIONS: Our finding that active regions emerge with an east-west alignment is consistent with earlier observations, but is still surprising since thin flux tube models predict that tilt angles of rising flux tubes are generated below the surface. Previously reported tilt angle relaxation of deeply anchored flux tubes can be largely explained by the change in east-west separation. We conclude that Joy's law is caused by an inherent north-south separation speed present when the flux first reaches the surface, and that the scatter in the tilt angle is consistent with buffeting of the polarities by supergranulation.

Average surface flows before the formation of solar active regions and their relationship to the supergranulation pattern.

CONTEXT: The emergence of solar active regions is an important but poorly understood aspect of the solar dynamo. AIMS: Knowledge of the flows associated with the rise of active-region-forming magnetic concentrations through the near-surface layers will help determine the mechanisms of active region formation. METHODS: We used helioseismic holography and granulation tracking to measure the horizontal flows at the surface that precede the emergence of active regions. We then averaged these flows over about sixty emerging active regions to reduce the noise, selecting active regions that emerge into relatively quiet Sun. To help interpret the results, we constructed a simple model flow field by generating synthetic "emergence locations" that are probabilistically related to the locations of supergranulation-scale convergence regions in the quiet Sun. RESULTS: The flow maps obtained from helioseismology and granulation tracking are very similar (correlation coefficients for single maps around 0.96). We find that active region emergence is, on average, preceded by converging horizontal flows of amplitude about 40 ms-1. The convergence region extends over about 40 Mm in the east-west direction and about 20 Mm in the north-south direction and is centered in the retrograde direction relative to the emergence location. This flow pattern is largely reproduced by a model in which active region emergence occurs preferentially in the prograde direction relative to supergranulation inflows. CONCLUSIONS: Averaging over many active regions reveals a statistically significant pattern of near-surface flows prior to emergence. The qualitative success of our simple model suggests that rising flux concentrations and supergranule-scale flows interact during the emergence process.

Flows around Averaged Solar Active Regions.

We explore the general properties of near-surface flows around solar active regions. Helioseismic holography is applied to HMI Dopplergrams yielding nearly 5000 flow measurements of 336 unique active regions observed by the Solar Dynamics Observatory between 2010 and 2014. Ensemble averages of the flows, over subsets of regions sorted on the basis of magnetic flux, are performed. These averages show that converging flows, with speeds of about 10 m s-1 and extending up to 10° from the active-region centers, are prevalent and have similar properties for all regions with magnetic flux above 1021 Mx. Retrograde flows are also detected, with amplitudes around 10 m s-1, which predominantly, but not exclusively, flank the polar side of the active regions. We estimate the expected contribution of these active-region flows to longitudinal averages of zonal and meridional flows and demonstrate the plausibility that they are responsible for at least some component of the time-varying global-scale flows. The reliability of our flow determination is tested using publicly available MHD simulations of both quiet-Sun convection and of a sunspot. While validating the overall methodology in general, the sunspot simulation demonstrates the presence of artifacts that may compromise quantitative flow inferences from some helioseismic measurements.

The Role of Twist in Kinked Flux Rope Emergence and Delta-spot Formation.

It has been observationally well established that the magnetic configurations most favorable for producing energetic flaring events reside in δ-spots, a class of sunspots defined as having opposite-polarity umbrae sharing a common penumbra. They are frequently characterized by extreme compactness, strong rotation, and anti-Hale orientation. Numerous studies have shown that nearly all of the largest solar flares originate in δ-spots, making the understanding of these structures a fundamental step in predicting space weather. Despite their important influence on the space environment, surprisingly little is understood about the origin and behavior of δ-spots. In this paper, we perform a systematic study of the behavior of emerging flux ropes to test a theoretical model for the formation of δ-spots: the kink instability of emerging flux ropes. We simulated the emergence of highly twisted, kink-unstable flux ropes from the convection zone into the corona, and we compared their photospheric properties to those of emerged weakly twisted, kink-stable flux ropes. We show that the photospheric manifestations of the emergence of highly twisted flux ropes closely match the observed properties of δ-spots, and we discuss the resulting implications for observations. Our results strongly support and extend previous theoretical work that suggested that the kink instability of emerging flux ropes is a promising candidate to explain δ-spot formation, as it reproduces their key characteristics very well.

Reflection of Fast Magnetosonic Waves near a Magnetic Reconnection Region.

Magnetic reconnection in the solar corona is thought to be unstable with the formation of multiple interacting plasmoids, and previous studies have shown that plasmoid dynamics can trigger MHD waves of different modes propagating outward from the reconnection site. However, variations in plasma parameters and magnetic field strength in the vicinity of a coronal reconnection site may lead to wave reflection and mode conversion. In this paper we investigate the reflection and refraction of fast magnetoacoustic waves near a reconnection site. Under a justified assumption of an analytically specified Alfvén speed profile, we derive and solve analytically the full wave equation governing the propagation of fast-mode waves in a non-uniform background plasma without recourse to the small wavelength approximation. We show that the waves undergo reflection near the reconnection current sheet due to the Alfvén speed gradient and that the reflection efficiency depends on the plasma-ß parameter, as well as on the wave frequency. In particular, we find that waves are reflected more efficiently near reconnection sites in a low-ß plasma, which is typical under solar coronal conditions. Also, the reflection is larger for lower-frequency waves while high-frequency waves propagate outward from the reconnection region almost without the reflection. We discuss the implications of efficient wave reflection near magnetic reconnection sites in strongly magnetized coronal plasma for particle acceleration, and also the effect this might have on first ionization potential (FIP) fractionation by the ponderomotive force of these waves in the chromosphere.

Reconnection-Driven Magnetohydrodynamic Turbulence in a Simulated Coronal-Hole Jet.

Extreme-ultraviolet and X-ray jets occur frequently in magnetically open coronal holes on the Sun, especially at high solar latitudes. Some of these jets are observed by white-light coronagraphs as they propagate through the outer corona toward the inner heliosphere, and it has been proposed that they give rise to microstreams and torsional Alfvén waves detected in situ in the solar wind. To predict and understand the signatures of coronal-hole jets, we have performed a detailed statistical analysis of such a jet simulated with an adaptively refined magnetohydrodynamics model. The results confirm the generation and persistence of three-dimensional, reconnection-driven magnetic turbulence in the simulation. We calculate the spatial correlations of magnetic fluctuations within the jet and find that they agree best with the Müller-Biskamp scaling model including intermittent current sheets of various sizes coupled via hydrodynamic turbulent cascade. The anisotropy of the magnetic fluctuations and the spatial orientation of the current sheets are consistent with an ensemble of nonlinear Alfvén waves. These properties also reflect the overall collimated jet structure imposed by the geometry of the reconnecting magnetic field. A comparison with Ulysses observations shows that turbulence in the jet wake is in quantitative agreement with that in the fast solar wind.

Magnetism, dynamo action and the solar-stellar connection.

The Sun and other stars are magnetic: magnetism pervades their interiors and affects their evolution in a variety of ways. In the Sun, both the fields themselves and their influence on other phenomena can be uncovered in exquisite detail, but these observations sample only a moment in a single star's life. By turning to observations of other stars, and to theory and simulation, we may infer other aspects of the magnetism-e.g., its dependence on stellar age, mass, or rotation rate-that would be invisible from close study of the Sun alone. Here, we review observations and theory of magnetism in the Sun and other stars, with a partial focus on the "Solar-stellar connection": i.e., ways in which studies of other stars have influenced our understanding of the Sun and vice versa. We briefly review techniques by which magnetic fields can be measured (or their presence otherwise inferred) in stars, and then highlight some key observational findings uncovered by such measurements, focusing (in many cases) on those that offer particularly direct constraints on theories of how the fields are built and maintained. We turn then to a discussion of how the fields arise in different objects: first, we summarize some essential elements of convection and dynamo theory, including a very brief discussion of mean-field theory and related concepts. Next we turn to simulations of convection and magnetism in stellar interiors, highlighting both some peculiarities of field generation in different types of stars and some unifying physical processes that likely influence dynamo action in general. We conclude with a brief summary of what we have learned, and a sampling of issues that remain uncertain or unsolved.

Magnetic Nulls and Super-radial Expansion in the Solar Corona.

Magnetic fields in the Sun's outer atmosphere-the corona-control both solar-wind acceleration and the dynamics of solar eruptions. We present the first clear observational evidence of coronal magnetic nulls in off-limb linearly polarized observations of pseudostreamers, taken by the Coronal Multichannel Polarimeter (CoMP) telescope. These nulls represent regions where magnetic reconnection is likely to act as a catalyst for solar activity. CoMP linear-polarization observations also provide an independent, coronal proxy for magnetic expansion into the solar wind, a quantity often used to parameterize and predict the solar wind speed at Earth. We introduce a new method for explicitly calculating expansion factors from CoMP coronal linear-polarization observations, which does not require photospheric extrapolations. We conclude that linearly polarized light is a powerful new diagnostic of critical coronal magnetic topologies and the expanding magnetic flux tubes that channel the solar wind.

SEARCH FOR A SIGNATURE OF TWIST-REMOVAL IN THE MAGNETIC FIELD OF SUNSPOTS IN RELATION WITH MAJOR FLARES.

We investigate the restructuring of the magnetic field in sunspots associated with two flares: the X6.5 flare on 6 December 2006 and the X2.2 flare on 15 February 2011. The observed changes were evaluated with respect to the so-called twist-removal model, in which helicity (twist) is removed from the corona as the result of an eruption. Since no vector magnetograms were available for the X6.5 flare, we applied the azimuthal symmetry approach to line-of-sight magnetograms to reconstruct the pseudo-vector magnetic field and investigate the changes in average twist and inclination of magnetic field in the sunspot around the time of the flare. For the X2.2 flare, results from the full vector magnetograms were compared with the pseudo-vector field data. For both flares, the data show changes consistent with the twist-removal scenario. We also evaluate the validity of the azimuthal symmetry approach on simple isolated round sunspots. In general, the derivations based on the azimuthal symmetry approach agree with true-vector field data though we find that even for symmetric sunspots the distribution of the magnetic field may deviate from an axially symmetric distribution.

Helioseismic Holography of Simulated Sunspots: dependence of the travel time on magnetic field strength and Wilson depression.

Improving methods for determining the subsurface structure of sunspots from their seismic signature requires a better understanding of the interaction of waves with magnetic field concentrations. We aim to quantify the impact of changes in the internal structure of sunspots on local helioseismic signals. We have numerically simulated the propagation of a stochastic wave field through sunspot models with different properties, accounting for changes in the Wilson depression between 250 and 550 km and in the photospheric umbral magnetic field between 1500 and 3500 G. The results show that travel-time shifts at frequencies above approximately 3.50 mHz (depending on the phase-speed filter) are insensitive to the magnetic field strength. The travel time of these waves is determined exclusively by the Wilson depression and sound-speed perturbation. The travel time of waves with lower frequencies is affected by the direct effect of the magnetic field, although photospheric field strengths below 1500 G do not leave a significant trace on the travel-time measurements. These results could potentially be used to develop simplified travel-time inversion methods.

Three-Dimensional Simulations of Tearing and Intermittency in Coronal Jets.

Observations of coronal jets increasingly suggest that local fragmentation and intermittency play an important role in the dynamics of these events. In this work we investigate this fragmentation in high-resolution simulations of jets in the closed-field corona. We study two realizations of the embedded-bipole model, whereby impulsive helical outflows are driven by reconnection between twisted and untwisted field across the domed fan plane of a magnetic null. We find that the reconnection region fragments following the onset of a tearing-like instability, producing multiple magnetic null points and flux-rope structures within the current layer. The flux ropes formed within the weak-field region in the center of the current layer are associated with "blobs" of density enhancement that become filamentary threads as the flux ropes are ejected from the layer, whereupon new flux ropes form behind them. This repeated formation and ejection of flux ropes provides a natural explanation for the intermittent outflows, bright blobs of emission, and filamentary structure observed in some jets. Additional observational signatures of this process are discussed. Essentially all jet models invoke reconnection between regions of locally closed and locally open field as the jet-generation mechanism. Therefore, we suggest that this repeated tearing process should occur at the separatrix surface between the two flux systems in all jets. A schematic picture of tearing-mediated jet reconnection in three dimensions is outlined.

Solar Coronal Jets: Observations, Theory, and Modeling.

Chromospheric and 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 major 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.

HARD X-RAY EMISSION DURING FLARES AND PHOTOSPHERIC FIELD CHANGES.

We study the correlation between abrupt permanent changes of magnetic field during X-class flares observed by the Global Oscillation Network Group and Helioseismic and Magnetic Imager instruments, and the hard X-ray (HXR) emission observed by RHESSI, to relate the photospheric field changes to the coronal restructuring and investigate the origin of the field changes. We find that spatially the early RHESSI emission corresponds well to locations of the strong field changes. The field changes occur predominantly in the regions of strong magnetic field near the polarity inversion line (PIL). The later RHESSI emission does not correspond to significant field changes as the flare footpoints are moving away from the PIL. Most of the field changes start before or around the start time of the detectable HXR signal, and they end at about the same time or later than the detectable HXR flare emission. Some of the field changes propagate with speed close to that of the HXR footpoint at a later phase of the flare. The propagation of the field changes often takes place after the strongest peak in the HXR signal when the footpoints start moving away from the PIL, i.e., the field changes follow the same trajectory as the HXR footpoint, but at an earlier time. Thus, the field changes and HXR emission are spatio-temporally related but not co-spatial nor simultaneous. We also find that in the strongest X-class flares the amplitudes of the field changes peak a few minutes earlier than the peak of the HXR signal. We briefly discuss this observed time delay in terms of the formation of current sheets during eruptions.