Showing 1–50 of 92 results for author: Rempel, M

Observations of Locally Excited Waves in the Low Solar Atmosphere Using the Daniel K. Inouye Solar Telescope (DKIST)

Authors: Shah Mohammad Bahauddin, Catherine E. Fischer, Mark P. Rast, Ivan Milic, Friedrich Woeger, Matthias Rempel, Peter H. Keys, Thomas R. Rimmele

Abstract: We present an interpretation of the recent Daniel K. Inouye Solar Telescope (DKIST) observations of propagating wavefronts in the lower solar atmosphere. Using MPS/University of Chicago MHD (MURaM) radiative magnetohydrodynamic simulations spanning solar photosphere, overshoot region, and lower chromosphere, we identify three acoustic-wave source mechanisms, each occurring at a different atmospher… ▽ More We present an interpretation of the recent Daniel K. Inouye Solar Telescope (DKIST) observations of propagating wavefronts in the lower solar atmosphere. Using MPS/University of Chicago MHD (MURaM) radiative magnetohydrodynamic simulations spanning solar photosphere, overshoot region, and lower chromosphere, we identify three acoustic-wave source mechanisms, each occurring at a different atmospheric height. We synthesize the DKIST Visible Broadband Imager (VBI) G-band, blue-continuum, and CaIIK signatures of these waves at high spatial and temporal resolution, and conclude that the wavefronts observed by DKIST likely originate from acoustic sources at the top of the solar photosphere overshoot region and in the chromosphere proper. The overall importance of these local sources to the atmospheric energy and momentum budget of the solar atmosphere is unknown, but one of the excitation mechanism identified (upward propagating shock interaction with down-welling chromospheric plasma resulting in acoustic radiation) appears to be an important shock dissipation mechanism. Additionally, the observed wavefronts may prove useful for ultra-local helioseismological inversions and promise to play an important diagnostic role at multiple atmospheric heights. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: For associated video (mp4) files see https://drive.google.com/file/d/1vWrXVsaazyk--0WSaFXp_shNrZug0cu7/view?usp=drive_link, https://drive.google.com/file/d/12VL2O2O5QztxUVxyB0HqNeGdnKH-ti2Y/view?usp=drive_link, and https://drive.google.com/file/d/1ARmq-m7UTHofPAuQf-p4_9yooRFaJUe9/view?usp=drive_link

Magnetic flux in the Sun emerges unaffected by supergranular-scale surface flows

Authors: Prasad Mani, Chris S. Hanson, Siddarth Dhanpal, Shravan Hanasoge, Srijan Bharati Das, Matthias Rempel

Abstract: Magnetic flux emergence from the convection zone into the photosphere and beyond is a critical component of the behaviour of large-scale solar magnetism. Flux rarely emerges amid field-free areas at the surface, but when it does, the interaction between magnetism and plasma flows can be reliably explored. Prior ensemble studies identified weak flows forming near emergence locations, but the low si… ▽ More Magnetic flux emergence from the convection zone into the photosphere and beyond is a critical component of the behaviour of large-scale solar magnetism. Flux rarely emerges amid field-free areas at the surface, but when it does, the interaction between magnetism and plasma flows can be reliably explored. Prior ensemble studies identified weak flows forming near emergence locations, but the low signal-to-noise required averaging over the entire dataset, erasing information about variation across the sample. Here, we apply deep learning to achieve improved signal-to-noise, enabling a case-by-case study. We find that these associated flows are dissimilar across instances of emergence and also occur frequently in the quiet convective background. Our analysis suggests diminished influence of supergranular-scale convective flows and magnetic buoyancy on flux rise. Consistent with numerical evidence, we speculate that small-scale surface turbulence and / or deep-convective processes play an outsize role in driving flux emergence. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 20 pages, 13 figures, Accepted in ApJ

Spatial resolution effects on the solar open flux estimates

Authors: Ivan Milic, Rebecca Centeno, Xudong Sun, Matthias Rempel, Jaime de la Cruz Rodriguez

Abstract: Spectropolarimetric observations used to infer the solar magnetic fields are obtained with a limited spatial resolution. The effects of this limited resolution on the inference of the open flux over the observed region have not been extensively studied. We aim to characterize the biases that arise in the inference of the mean flux density by performing an end-to-end study that involves the generat… ▽ More Spectropolarimetric observations used to infer the solar magnetic fields are obtained with a limited spatial resolution. The effects of this limited resolution on the inference of the open flux over the observed region have not been extensively studied. We aim to characterize the biases that arise in the inference of the mean flux density by performing an end-to-end study that involves the generation of synthetic data, its interpretation (inversion), and a comparison of the results with the original model. We synthesized polarized spectra of the two magnetically sensitive lines of neutral iron around 630\,nm from a state-of-the-art numerical simulation of the solar photosphere. We then performed data degradation to simulate the effect of the telescope with a limited angular resolution and interpreted (inverted) the data using a Milne-Eddington spectropolarimetric inversion code. We then studied the dependence of the inferred parameters on the telescope resolution. The results show a significant decrease in the mean magnetic flux density -- related to the open flux observed at the disk center -- with decreasing telescope resolution. The original net magnetic field flux is fully resolved by a 1m telescope, but a 20\,cm aperture telescope yields a 30\% smaller value. Even in the fully resolved case, the result is still biased due to the corrugation of the photospheric surface. Even the spatially averaged quantities, such as the open magnetic flux in the observed region, are underestimated when the magnetic structures are unresolved. The reason for this is the presence of nonlinearities in the magnetic field inference process. This effect might have implications for the modeling of large-scale solar magnetic fields; for example, those corresponding to the coronal holes, or the polar magnetic fields, which are relevant to our understanding of the solar cycle. △ Less

Submitted 4 February, 2024; originally announced February 2024.

Comments: To be published in A&A

Advancing solar magnetic field extrapolations through multi-height magnetic field measurements

Authors: Robert Jarolim, Benoit Tremblay, Matthias Rempel, Momchil Molnar, Astrid M. Veronig, Julia K. Thalmann, Tatiana Podladchikova

Abstract: Non-linear force-free extrapolations are a common approach to estimate the 3D topology of coronal magnetic fields based on photospheric vector magnetograms. The force-free assumption is a valid approximation at coronal heights, but for the dense plasma conditions in the lower atmosphere, this assumption is not satisfied. In this study, we utilize multi-height magnetic field measurements in combina… ▽ More Non-linear force-free extrapolations are a common approach to estimate the 3D topology of coronal magnetic fields based on photospheric vector magnetograms. The force-free assumption is a valid approximation at coronal heights, but for the dense plasma conditions in the lower atmosphere, this assumption is not satisfied. In this study, we utilize multi-height magnetic field measurements in combination with physics-informed neural networks to advance solar magnetic field extrapolations. We include a flexible height-mapping, which allows us to account for the different formation heights of the observed magnetic field measurements. The comparison to analytical and simulated magnetic fields demonstrates that including chromospheric magnetic field measurements leads to a significant improvement of our magnetic field extrapolations. We also apply our method to chromospheric line-of-sight magnetograms, from the Vector Spectromagnetograph (VSM) on the Synoptic Optical Long-term Investigations of the Sun (SOLIS) observatory, in combination with photospheric vector magnetograms, from the Helioseismic Magnetic Imager (HMI) onboard the Solar Dynamic Observatory (SDO). The comparison to observations in extreme ultraviolet wavelengths shows that the additional chromospheric information leads to a better agreement with the observed coronal structures. In addition, our method intrinsically provides an estimate of the corrugation of the observed magnetograms. With this new approach, we make efficient use of multi-height magnetic field measurements and advance the realism of coronal magnetic field simulations. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Quantifying Poynting flux in the Quiet Sun Photosphere

Authors: Dennis Tilipman, Maria Kazachenko, Benoit Tremblay, Ivan Milic, Valentin Martinez Pillet, Matthias Rempel

Abstract: Poynting flux is the flux of magnetic energy, which is responsible for chromospheric and coronal heating in the solar atmosphere. It is defined as a cross product of electric and magnetic fields, and in ideal MHD conditions it can be expressed in terms of magnetic field and plasma velocity. Poynting flux has been computed for active regions and plages, but estimating it in the quiet Sun (QS) remai… ▽ More Poynting flux is the flux of magnetic energy, which is responsible for chromospheric and coronal heating in the solar atmosphere. It is defined as a cross product of electric and magnetic fields, and in ideal MHD conditions it can be expressed in terms of magnetic field and plasma velocity. Poynting flux has been computed for active regions and plages, but estimating it in the quiet Sun (QS) remains challenging due to resolution effects and polarimetric noise. However, with upcoming DKIST capabilities, these estimates will become more feasible than ever before. Here, we study QS Poynting flux in Sunrise/IMaX observations and MURaM simulations. We explore two methods for inferring transverse velocities from observations - FLCT and a neural network based method DeepVel - and show DeepVel to be the more suitable method in the context of small-scale QS flows. We investigate the effect of azimuthal ambiguity on Poynting flux estimates, and we describe a new method for azimuth disambiguation. Finally, we use two methods for obtaining the electric field. The first method relies on idealized Ohm's law, whereas the second is a state-of-the-art inductive electric field inversion method PDFI SS. We compare the resulting Poynting flux values with theoretical estimates for chromospheric and coronal energy losses and find that some of Poynting flux estimates are sufficient to match the losses. Using MURaM simulations, we show that photospheric Poynting fluxes vary significantly with optical depth, and that there is an observational bias that results in underestimated Poynting fluxes due to unaccounted shear term contribution. △ Less

Submitted 5 July, 2023; originally announced July 2023.

A study of the capabilities for inferring atmospheric information from high-spatial-resolution simulations

Authors: C. Quintero Noda, E. Khomenko, M. Collados, B. Ruiz Cobo, R. Gafeira, N. Vitas, M. Rempel, R. J. Campbell, A. Pastor Yabar, H. Uitenbroek, D. Orozco Suárez

Abstract: In this work, we study the accuracy that can be achieved when inferring the atmospheric information from realistic numerical magneto-hydrodynamic simulations that reproduce the spatial resolution we will obtain with future observations made by the 4m class telescopes DKIST and EST. We first study multiple inversion configurations using the SIR code and the Fe I transitions at 630 nm until we obtai… ▽ More In this work, we study the accuracy that can be achieved when inferring the atmospheric information from realistic numerical magneto-hydrodynamic simulations that reproduce the spatial resolution we will obtain with future observations made by the 4m class telescopes DKIST and EST. We first study multiple inversion configurations using the SIR code and the Fe I transitions at 630 nm until we obtain minor differences between the input and the inferred atmosphere in a wide range of heights. Also, we examine how the inversion accuracy depends on the noise level of the Stokes profiles. The results indicate that when the majority of the inverted pixels come from strongly magnetised areas, there are almost no restrictions in terms of the noise, obtaining good results for noise amplitudes up to 1$\times10^{-3}$ of $I_c$. At the same time, the situation is different for observations where the dominant magnetic structures are weak, and noise restraints are more demanding. Moreover, we find that the accuracy of the fits is almost the same as that obtained without noise when the noise levels are on the order of 1$\times10^{-4}$of $I_c$. We, therefore, advise aiming for noise values on the order of or lower than 5$\times10^{-4}$ of $I_c$ if observers seek reliable interpretations of the results for the magnetic field vector reliably. We expect those noise levels to be achievable by next-generation 4m class telescopes thanks to an optimised polarisation calibration and the large collecting area of the primary mirror. △ Less

Submitted 2 June, 2023; originally announced June 2023.

Comments: 14 pages, 13 figures

The Puzzling Structure of Solar Convection: Window into the Dynamo

Authors: Nicholas A. Featherstone, Kyle C. Augustson, Jonathan M. Aurnou, Catherine Blume, Benjamin P. Brown, Nicholas Brummell, Keaton J. Burns, Michael A. Calkins, Maria Camisassa, Mausumi Dikpati, Yuhong Fan, J. R. Fuentes, Gustavo Guerrero, Bradley W. Hindman, Keith Julien, Irina N. Kitiashvili, Lydia Korre, Daniel Lecoanet, Bhishek Manek, Loren Matilsky, Mark Miesch, Nicholas J. Nelson, Jeffrey S. Oishi, Whitney T. Powers, Matthias Rempel , et al. (3 additional authors not shown)

Abstract: The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective… ▽ More The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective flows. Those flows are in-turn thought to sustain both the Sun's turbulent EMF and the large-scale flows of differential rotation and meridional circulation suspected of influencing the dynamo's organization and timing. Continued progress in this area is complicated by (i) inconsistencies between helioseismic measurements of convective and meridional flow made with different techniques and instruments, and (ii) a lack of high-latitude data for convection, differential rotation, and meridional flow. We suggest that the path forward to resolving these difficulties is twofold. First, the acquisition of long-term helioseismic and emissivity measurements obtained from a polar vantage point is vital to complete our picture of the Sun's outer convection zone. Second, sustained and expanded investment in theory-oriented and combined theory/observational research initiatives will be crucial to fully exploit these new observations and to resolve inconsistencies between existing measurements. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Comments: White paper submitted to the Decadal Survey for Solar and Space Physics 2024-2033. 7 pages; 4 figures

Small-scale dynamos: From idealized models to solar and stellar applications

Authors: Matthias Rempel, Tanayveer Bhatia, Luis Bellot Rubio, Maarit J. Korpi-Lagg

Abstract: In this article we review small-scale dynamo processes that are responsible for magnetic field generation on scales comparable to and smaller than the energy carrying scales of turbulence. We provide a review of critical observation of quiet Sun magnetism, which have provided strong support for the operation of a small-scale dynamo in the solar photosphere and convection zone. After a review of ba… ▽ More In this article we review small-scale dynamo processes that are responsible for magnetic field generation on scales comparable to and smaller than the energy carrying scales of turbulence. We provide a review of critical observation of quiet Sun magnetism, which have provided strong support for the operation of a small-scale dynamo in the solar photosphere and convection zone. After a review of basic concepts we focus on numerical studies of kinematic growth and non-linear saturation in idealized setups, with special emphasis on the role of the magnetic Prandtl number for dynamo onset and saturation. Moving towards astrophysical applications we review convective dynamo setups that focus on the deep convection zone and the photospheres of solar-like stars. We review the critical ingredients for stellar convection setups and discuss their application to the Sun and solar-like stars including comparison against available observations. △ Less

Submitted 4 May, 2023; originally announced May 2023.

Comments: 54 pages, 10 figures, submitted to Space Science Reviews

Report number: NORDITA 2023-024

Limitations and biases in the retrieval of the polar magnetic field I: the role of the magnetic filling factor in Milne-Eddington inversions of simulated Hinode/SP data

Authors: Rebecca Centeno, Ivan Milić, Matthias Rempel, Nariaki V. Nitta, Xudong Sun

Abstract: We study the extent to which Milne-Eddington inversions are able to retrieve and characterize the magnetic landscape of the solar poles from observations by the spectropolarimeter onboard Hinode. In particular, we evaluate whether a variable magnetic filling factor is an adequate modeling technique for retrieving the intrinsic magnetic properties from every pixel in the polar field of view. We fir… ▽ More We study the extent to which Milne-Eddington inversions are able to retrieve and characterize the magnetic landscape of the solar poles from observations by the spectropolarimeter onboard Hinode. In particular, we evaluate whether a variable magnetic filling factor is an adequate modeling technique for retrieving the intrinsic magnetic properties from every pixel in the polar field of view. We first generate synthetic spectra emerging from a numerical simulation of a "plage" region at an inclined line of sight of 65$^{\circ}$, and degrade the data to emulate real observations. Then, we invert the synthetic spectra with two Milne-Eddington inversion codes that feature different treatments of the magnetic filling factor, and relate the retrieved magnetic quantities back to their original values in the simulation cube. We find that while the apparent retrieved magnetic properties map well the spatially-degraded simulation, the intrinsic magnetic quantities bear little relation to the magnetic field at the native resolution of the simulation. We discuss the systematic biases caused by line-of-sight foreshortening, spatial degradation, photon noise and modeling assumptions embedded in the inversion algorithm. △ Less

Submitted 1 May, 2023; originally announced May 2023.

Comments: Accepted for publication in ApJ

Eruption of a Magnetic Flux Rope in a Comprehensive Radiative Magnetohydrodynamic Simulation of flare-productive active regions

Authors: Feng Chen, Matthias Rempel, Yuhong Fan

Abstract: Radiative magnetohydrodynamic simulation includes sufficiently realistic physics to allow for the synthesis of remote sensing observables that can be quantitatively compared with observations. We analyze the largest flare in a simulation of the emergence of large flare-productive active regions described by Chen et al. The flare releases $4.5\times10^{31}$ erg of magnetic energy and is accompanied… ▽ More Radiative magnetohydrodynamic simulation includes sufficiently realistic physics to allow for the synthesis of remote sensing observables that can be quantitatively compared with observations. We analyze the largest flare in a simulation of the emergence of large flare-productive active regions described by Chen et al. The flare releases $4.5\times10^{31}$ erg of magnetic energy and is accompanied by a spectacular coronal mass ejection. Synthetic soft X-ray flux of this flare reaches M2 class. The eruption reproduces many key features of observed solar eruptions. A pre-existing magnetic flux rope is formed along the highly sheared polarity inversion line between a sunspot pair and is covered by an overlying multi-pole magnetic field. During the eruption, the progenitor flux rope actively reconnects with the canopy field and evolves to the large-scale multi-thermal flux rope that is observed in the corona. Meanwhile, the magnetic energy released via reconnection is channeled down to the lower atmosphere and gives rise to bright soft X-ray post-flare loops and flare ribbons that reproduce the morphology and dynamic evolution of observed flares. The model helps to shed light on questions of where and when the a flux rope may form and how the magnetic structures in an eruption are related to observable emission properties. △ Less

Submitted 31 May, 2023; v1 submitted 9 March, 2023; originally announced March 2023.

Comments: 13 Pages, 6 figures, accepted for publication in ApJL

Comprehensive radiative MHD simulations of eruptive flares above collisional polarity inversion lines

Authors: Matthias Rempel, Georgios Chintzoglou, Mark C. M. Cheung, Yuhong Fan, Lucia Kleint

Abstract: We present a new simulation setup using the MURaM radiative Magnetohydrodynamic (MHD) code that allows to study the formation of collisional polarity inversion lines (cPILs) in the photosphere and the coronal response including flares. In this scheme we start with a bipolar sunspot configuration and set the spots on collision course by imposing the appropriate velocity field at the footpoints in t… ▽ More We present a new simulation setup using the MURaM radiative Magnetohydrodynamic (MHD) code that allows to study the formation of collisional polarity inversion lines (cPILs) in the photosphere and the coronal response including flares. In this scheme we start with a bipolar sunspot configuration and set the spots on collision course by imposing the appropriate velocity field at the footpoints in the subphotospheric boundary. We produce different setups with the same initial spot separation by varying physical parameters such as the collision speed and minimum collision distance. While all setups lead to the formation of an EUV and X-ray sigmoid structure, only the cases with a close passing of the spots cause flares and mass eruptions. The energy release is in the $1-2 \times 10^{31}$ erg range, putting the simulated flares into the upper C-class to lower M-class range of GOES X-ray 1-8 A flux. While the setup with the more distant passing of the spots does not lead to a flare, the corona is nonetheless substantially heated, suggesting non-eruptive energy release mechanisms. We focus our discussion on two particular setups that differ in spot coherence and resulting cPIL length persistence. We find different timings in the transition from a sheared magnetic arcade (SMA) to magnetic flux rope (MFR); the setup with a large length but shorter duration cPIL produces a MFR during the eruption, while the MFR is pre-existing in the setup with a large length and longer duration cPIL. While both result in flares of comparable strength and the eruption of a CME, the setup with pre-existing MFR (and embedded filament) leads to an MFR eruption with a larger mass content. △ Less

Submitted 25 August, 2023; v1 submitted 9 March, 2023; originally announced March 2023.

Comments: 19 pages, 11 figures, 1 table, ApJ in press (replaced previous version with the accepted article)

Data-Driven Radiative Magnetohydrodynamics Simulations with the MURaM code

Authors: Feng Chen, Mark C. M. Cheung, Matthias Rempel, Georgios Chintzoglou

Abstract: We present a method of conducting data-driven simulations of solar active regions and flux emergence with the MURaM radiative magnetohydrodynamics (MHD) code. The horizontal electric field derived from the full velocity and magnetic vectors, is implemented at the photospheric (bottom) boundary to drive the induction equation. The energy equation accounts for thermal conduction along magnetic field… ▽ More We present a method of conducting data-driven simulations of solar active regions and flux emergence with the MURaM radiative magnetohydrodynamics (MHD) code. The horizontal electric field derived from the full velocity and magnetic vectors, is implemented at the photospheric (bottom) boundary to drive the induction equation. The energy equation accounts for thermal conduction along magnetic fields, optically-thin radiative loss, and heating of coronal plasma by viscous and resistive dissipation, which allows for a realistic presentation of the thermodynamic properties of coronal plasma that are key to predicting the observational features of solar active regions and eruptions. To validate the method, the photospheric data from a comprehensive radiative MHD simulation of solar eruption (the ground truth) are used to drive a series of numerical experiments. The data-driven simulation reproduces the accumulation of free magnetic energy over the course of flux emergence in the ground truth with an error of 3\%. The onset time is approximately 8\,min delayed compared to the ground truth. However, a precursor-like signature can be identified at the correct onset time. The data-driven simulation captures key eruption-related emission features and plasma dynamics of the ground truth flare over a wide temperature span from $\log_{10}T{=}4.5$ to $\log_{10}T{>}8$. The evolution of the flare and coronal mass ejection as seen in synthetic extreme ultraviolet images is also reproduced with high fidelity. The method helps to understand the evolution of magnetic field in a more realistic coronal environment and to link the magnetic structures to observable diagnostics. △ Less

Submitted 29 March, 2023; v1 submitted 18 January, 2023; originally announced January 2023.

Comments: 35 pages, 24 figures, accepted for publication in ApJ

Towards data-driven modeling and real-time prediction of solar flares and coronal mass ejections

Authors: M. Rempel, Y. Fan, M. Dikpati, A. Malanushenko, M. D. Kazachenko, M. C. M. Cheung, G. Chintzoglou, X. Sun, G. H. Fisher, T. Y. Chen

Abstract: Modeling of transient events in the solar atmosphere requires the confluence of 3 critical elements: (1) model sophistication, (2) data availability, and (3) data assimilation. This white paper describes required advances that will enable statistical flare and CME forecasting (e.g. eruption probability and timing, estimation of strength, and CME details, such as speed and magnetic field orientatio… ▽ More Modeling of transient events in the solar atmosphere requires the confluence of 3 critical elements: (1) model sophistication, (2) data availability, and (3) data assimilation. This white paper describes required advances that will enable statistical flare and CME forecasting (e.g. eruption probability and timing, estimation of strength, and CME details, such as speed and magnetic field orientation) similar to weather prediction on Earth. △ Less

Submitted 29 December, 2022; originally announced December 2022.

Comments: Heliophysics 2050 White Paper

Small-scale loops heated to transition region temperatures and their chromospheric signatures in the simulated solar atmosphere

Authors: M. Skan, S. Danilovic, J. Leenaarts, F. Calvo, M. Rempel

Abstract: Recent observations revealed loop-like structures at very small scales visible in observables that sample transition region (TR) and coronal temperatures. Their formation remains unclear. We study an example of a bipolar system in realistic magnetohydrodynamic simulations and forward synthesis of spectral lines to investigate how these features occur. Computations are done using the MURaM code… ▽ More Recent observations revealed loop-like structures at very small scales visible in observables that sample transition region (TR) and coronal temperatures. Their formation remains unclear. We study an example of a bipolar system in realistic magnetohydrodynamic simulations and forward synthesis of spectral lines to investigate how these features occur. Computations are done using the MURaM code to generate model atmospheres. The synthetic H$α$ and Si IV spectra are calculated at two angles ($μ= 1$, $μ= 0.66$) using the Multi3D code. Magnetic field lines are traced in the model and the evolution of the underlying field topology is examined. The synthetic H$α$ dopplergrams reveal loops that evolve dramatically within a few minutes. The synthetic H$α$ line profiles show observed asymmetries and doppler shifts in the line core. They, however, also show strong emission peaks in the line wings, even at the slated view. The synthetic Si IV emission features partly coincide with structures visible in H$α$ dopplergrams and partly follow separate magnetic field threads. Some are even visible in the emission measure maps for the lg$(T /K)= [5.0, 5.5]$ temperature interval. The emission areas trace out the magnetic field lines rooted in opposite polarities in a bipolar region. We find that our results largely reproduce the observed features and their characteristics. A bipolar system with footpoints undergoing rapid movement and shuffling can produce many small-scale recurrent events heated to high temperatures. The morphology and evolution of the resulting observable features can vary depending on the viewing angle. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: Online material: https://dubshen.astro.su.se/~mosk1498/mov/fig16_4p.mp4

Journal ref: A&A 672, A47 (2023)

Characterization of magneto-convection in sunspots -- The Gough & Tayler stability criterion in MURaM sunspot simulations

Authors: M. Schmassmann, M. Rempel, N. Bello González, R. Schlichenmaier, J. Jurčák

Abstract: We analyse a sunspot simulation in an effort to understand the origin of the convective instabilities giving rise to the penumbral and umbral distinct regimes. We applied the criterion from Gough \& Tayler (1966), accounting for the stabilising effect of the vertical magnetic field to investigate the convective instabilities in a MURaM sunspot simulation. We find: (1) a highly unstable shallow lay… ▽ More We analyse a sunspot simulation in an effort to understand the origin of the convective instabilities giving rise to the penumbral and umbral distinct regimes. We applied the criterion from Gough \& Tayler (1966), accounting for the stabilising effect of the vertical magnetic field to investigate the convective instabilities in a MURaM sunspot simulation. We find: (1) a highly unstable shallow layer right beneath the surface extending all over the simulation box in which convection is triggered by radiative cooling in the photosphere; (2) a deep umbral core (beneath -5 Mm) stabilised against overturning convection that underlies a region with stable background values permeated by slender instabilities coupled to umbral dots; (3) filamentary instabilities below the penumbra nearly parallel to the surface and undulating instabilities coupled to the penumbra which originate in the deep layers. These deep-rooted instabilities result in the vigorous magneto-convection regime characteristic of the penumbra; (4) convective downdrafts in the granulation, penumbra, and umbra develop at about 2 km/s, 1 km/s, and 0.1 km/s, respectively, indicating that the granular regime of convection is more vigorous than the penumbra convection regime which, in turn, is more vigorous than the close-to-steady umbra; (5) the GT criterion outlines both the sunspot magnetopause and peripatopause, highlighting the tripartite nature of the sub-photospheric layers of magnetohydrodynamic (MHD) sunspot models; and, finally, (6) the Jurčák criterion is the photospheric counterpart of the GT criterion in deep layers. The GT criterion as a diagnostic tool reveals the tripartite nature of sunspot structure with distinct regimes of magneto-convection in the umbra, penumbra, and granulation operating in realistic MHD simulations. △ Less

Submitted 31 October, 2022; originally announced October 2022.

Comments: Published in Astronomy and Astrophysics, Volume 656, id.A92, 8 pp

Report number: Volume 656, id.A92, 8 pp

Journal ref: 2021A&A...656A..92S

Rapid Blue- and Red-shifted Excursions in H$α$ line profiles synthesized from realistic 3D MHD simulations

Authors: S. Danilovic, J. P. Bjørgen, J. Leenaarts, M. Rempel

Abstract: Rapid blue- and red-shifted events (RBEs/RREs) may have an important role in mass-loading and heating the solar corona, but their nature and origin are still debatable. We aim to model these features to learn more about their properties, formation and origin. A realistic three-dimensional (3D) magneto-hydrodynamic (MHD) model of a solar plage region is created. Synthetic H$α$ spectra are generated… ▽ More Rapid blue- and red-shifted events (RBEs/RREs) may have an important role in mass-loading and heating the solar corona, but their nature and origin are still debatable. We aim to model these features to learn more about their properties, formation and origin. A realistic three-dimensional (3D) magneto-hydrodynamic (MHD) model of a solar plage region is created. Synthetic H$α$ spectra are generated and the spectral signatures of these features are identified. The magnetic field lines associated with these events are traced and the underlying dynamic is studied. The model reproduces well many properties of RBEs and RREs, such as spatial distribution, lateral movement, length and lifetimes. Synthetic H$α$ line profiles, similarly to observed ones, show strong blue- or red-shift and asymmetries. These line profiles are caused by the vertical component of velocity with magnitudes larger than $30-40$ km/s that appear mostly in the height range of $2-4$ Mm. By tracing magnetic field lines, we show that the vertical velocity that causes the appearance of RBE/RREs to appear is always associated with the component of velocity perpendicular to the magnetic field line. The study confirms the hypothesis that RBEs and RREs are signs of Alfv{é}nic waves with, in some cases, a significant contribution from slow magneto-acoustic mode. △ Less

Submitted 29 August, 2022; originally announced August 2022.

Comments: submitted to A&A, movies in https://dubshen.astro.su.se/~sdani/rbes/

Journal ref: A&A 670, A50 (2023)

Effects of spectral resolution on simple magnetic field diagnostics of the Mg II h & k lines

Authors: Rebecca Centeno, Matthias Rempel, Roberto Casini, Tanausu del Pino Aleman

Abstract: We study the effects of finite spectral resolution on the magnetic field values retrieved through the weak field approximation (WFA) from the cores of the Mg II h & k lines. The retrieval of the line-of-sight (LOS) component of the magnetic field, $B_{\rm LOS}$, from synthetic spectra generated in a uniformly magnetized FAL-C atmosphere are accurate when restricted to the inner lobes of Stokes V.… ▽ More We study the effects of finite spectral resolution on the magnetic field values retrieved through the weak field approximation (WFA) from the cores of the Mg II h & k lines. The retrieval of the line-of-sight (LOS) component of the magnetic field, $B_{\rm LOS}$, from synthetic spectra generated in a uniformly magnetized FAL-C atmosphere are accurate when restricted to the inner lobes of Stokes V. As we degrade the spectral resolution, partial redistribution (PRD) effects, that more prominently affect the outer lobes of Stokes V, are brought into the line core through spectral smearing, degrading the accuracy of the WFA and resulting in an inference bias, which is more pronounced the poorer the resolution. When applied to a diverse set of spectra emerging from a sunspot simulation, we find a good accuracy in the retrieved $B_{\rm LOS}$ when comparing it to the model value at the height where the optical depth in the line core is unity. The accuracy is preserved up to field strengths of B~1500 G. Limited spectral resolution results in a small bias toward weaker retrieved fields. The WFA for the transverse component of the magnetic field is also evaluated. Reduced spectral resolution degrades the accuracy of the inferences because spectral mixing results in the line effectively probing deeper layers of the atmosphere. △ Less

Submitted 15 August, 2022; originally announced August 2022.

Comments: Accepted for publication in ApJ

Chromospheric Extension of the MURaM Code

Authors: D. Przybylski, R. Cameron, S. K. Solanki, M. Rempel, J. Leenaarts, L. S. Anusha, V. Witzke, A. I. Shapiro

Abstract: Detailed numerical models of chromosphere and corona are required to understand the heating of the solar atmosphere. An accurate treatment of the solar chromosphere is complicated by the effects arising from Non Local Thermodynamic Equilibrium (NLTE) radiative transfer. A small number of strong, highly scattering lines dominate the cooling and heating in the chromosphere. Additionally, the recombi… ▽ More Detailed numerical models of chromosphere and corona are required to understand the heating of the solar atmosphere. An accurate treatment of the solar chromosphere is complicated by the effects arising from Non Local Thermodynamic Equilibrium (NLTE) radiative transfer. A small number of strong, highly scattering lines dominate the cooling and heating in the chromosphere. Additionally, the recombination times of ionised hydrogen are longer than the dynamical timescales, requiring a non-equilibrium (NE) treatment of hydrogen ionisation. The MURaM code is extended to include the physical process required for accurate simulation of the solar chromosphere, as implemented in the Bifrost code. This includes a time-dependent treatment of hydrogen ionisation, a scattering multi-group radiation transfer scheme and approximations for NLTE radiative cooling. The inclusion of NE and NLTE physics has a large impact on the structure of the chromosphere; the NE treatment of hydrogen ionisation leads to a higher ionisation fraction and enhanced populations in the first excited state throughout cold inter-shock regions of the chromosphere. Additionally this prevents hydrogen ioniation from buffering energy fluctuations, leading to hotter shocks and cooler inter-shock regions. The hydrogen populations in the ground and first excited state are enhanced by $10^2-10^3$ in the upper chromosphere and up to $10^9$ near the transition region. Including the necessary NLTE physics leads to significant differences in chromospheric structure and dynamics. The thermodynamics and hydrogen populations calculated using the extended version of the MURaM code are consistent with previous non-equilibrium simulations. The electron number and temperature calculated using the non-equilibrium treatment of the chromosphere are required to accurately synthesise chromospheric spectral lines. △ Less

Submitted 6 April, 2022; originally announced April 2022.

Comments: 24 pages, 12 figures, to be published in A&A

Journal ref: A&A 664, A91 (2022)

Heating of the solar chromosphere through current dissipation

Authors: J. M. da Silva Santos, S. Danilovic, J. Leenaarts, J. de la Cruz Rodríguez, X. Zhu, S. M. White, G. J. M. Vissers, M. Rempel

Abstract: The solar chromosphere is heated to temperatures higher than predicted by radiative equilibrium. This excess heating is greater in active regions where the magnetic field is stronger. We aim to investigate the magnetic topology associated with an area of enhanced millimeter (mm) brightness temperatures in a solar active region mapped by the Atacama Large Millimeter/submillimeter Array (ALMA) using… ▽ More The solar chromosphere is heated to temperatures higher than predicted by radiative equilibrium. This excess heating is greater in active regions where the magnetic field is stronger. We aim to investigate the magnetic topology associated with an area of enhanced millimeter (mm) brightness temperatures in a solar active region mapped by the Atacama Large Millimeter/submillimeter Array (ALMA) using spectropolarimetric co-observations with the 1-m Swedish Solar Telescope (SST). We used Milne-Eddington inversions, nonlocal thermodynamic equilibrium (non-LTE) inversions, and a magnetohydrostatic extrapolation to obtain constraints on the three-dimensional stratification of temperature, magnetic field, and radiative energy losses. We compared the observations to a snapshot of a magnetohydrodynamics simulation and investigate the formation of the thermal continuum at 3 mm using contribution functions. We find enhanced heating rates in the upper chromosphere of up to $\sim 5\rm\,kW\,m^{-2}$, where small-scale emerging loops interact with the overlying magnetic canopy leading to current sheets as shown by the magnetic field extrapolation. Our estimates are about a factor of two higher than canonical values, but they are limited by the ALMA spatial resolution ($\sim 1.2^{\prime\prime}$). Band 3 brightness temperatures reach about $\sim10^{4}\,$K in the region, and the transverse magnetic field strength inferred from the non-LTE inversions is on the order of $\sim 500\,$G in the chromosphere. We are able to quantitatively reproduce many of the observed features, including the integrated radiative losses in our numerical simulation. We conclude that the heating is caused by dissipation in current sheets. However, the simulation shows a complex stratification in the flux emergence region where distinct layers may contribute significantly to the emission in the mm continuum. △ Less

Submitted 19 February, 2022; v1 submitted 8 February, 2022; originally announced February 2022.

Comments: 18 pages, 13 figures. Accepted for publication in A&A. Typos and bibtex corrected

Journal ref: A&A 661, A59 (2022)

A solar coronal loop in a box: Energy generation and heating

Authors: C. Breu, H. Peter, R. Cameron, S. K. Solanki, D. Przybylski, M. Rempel, L. P. Chitta

Abstract: Coronal loops are the basic building block of the upper solar atmosphere. Comprehending how these are energized, structured, and evolve is key to understanding stellar coronae. Here we investigate how the energy to heat the loop is generated by photospheric magneto-convection, transported into the upper atmosphere, and how the internal structure of a coronal loop forms. In a 3D magnetohydrodyn… ▽ More Coronal loops are the basic building block of the upper solar atmosphere. Comprehending how these are energized, structured, and evolve is key to understanding stellar coronae. Here we investigate how the energy to heat the loop is generated by photospheric magneto-convection, transported into the upper atmosphere, and how the internal structure of a coronal loop forms. In a 3D magnetohydrodynamics (MHD) model, we study an isolated coronal loop rooted with both footpoints in a shallow layer within the convection zone using the MURaM code. To resolve its internal structure, we limited the computational domain to a rectangular box containing a single coronal loop as a straightened magnetic flux tube. Field-aligned heat conduction, gray radiative transfer in the photosphere and chromosphere, and optically thin radiative losses in the corona were taken into account. The footpoints were allowed to interact self-consistently with the granulation surrounding them. The loop is heated by a Poynting flux that is self-consistently generated through small-scale motions within individual magnetic concentrations in the photosphere. Turbulence develops in the upper layers of the atmosphere as a response to the footpoint motions. We see little sign of heating by large-scale braiding of magnetic flux tubes from different photospheric concentrations at a given footpoint. The synthesized emission, as it would be observed by the Atmospheric Imaging Assembly or the X-ray Telescope, reveals transient bright strands that form in response to the heating events. Overall, our model roughly reproduces the properties and evolution of the plasma as observed within coronal loops. With this model we can build a coherent picture of how the energy flux to heat the upper atmosphere is generated near the solar surface and how this process drives and governs the heating and dynamics of a coronal loop. △ Less

Submitted 21 December, 2021; originally announced December 2021.

Comments: 19 pages, 16 figures

Convolutional Neural Networks and Stokes Response Functions

Authors: Rebecca Centeno, Natasha Flyer, Lipi Mukherjee, Ricky Egeland, Roberto Casini, Tanausu del Pino Aleman, Matthias Rempel

Abstract: In this work, we study the information content learned by a convolutional neural network (CNN) when trained to carry out the inverse mapping between a database of synthetic Ca II intensity spectra and the vertical stratification of the temperature of the atmospheres used to generate such spectra. In particular, we evaluate the ability of the neural network to extract information about the sensitiv… ▽ More In this work, we study the information content learned by a convolutional neural network (CNN) when trained to carry out the inverse mapping between a database of synthetic Ca II intensity spectra and the vertical stratification of the temperature of the atmospheres used to generate such spectra. In particular, we evaluate the ability of the neural network to extract information about the sensitivity of the spectral line to temperature as a function of height. By training the CNN on sufficiently narrow wavelength intervals across the Ca II spectral profiles, we find that the error in the temperature prediction shows an inverse relationship to the response function of the spectral line to temperature, this is, different regions of the spectrum yield a better temperature prediction at their expected regions of formation. This work shows that the function that the CNN learns during the training process contains a physically-meaningful mapping between wavelength and atmospheric height. △ Less

Submitted 7 December, 2021; originally announced December 2021.

Comments: Accepted for publication in ApJ

Efficient numerical treatment of ambipolar and Hall drift as hyperbolic system

Authors: M. Rempel, D. Przybylski

Abstract: Partially ionized plasmas, such as the solar chromosphere, require a generalized Ohm's law including the effects of ambipolar and Hall drift. While both describe transport processes that arise from the multifluid equations and are therefore of hyperbolic nature, they are often incorporated in models as a diffusive, i.e. parabolic process. While the formulation as such is easy to include in standar… ▽ More Partially ionized plasmas, such as the solar chromosphere, require a generalized Ohm's law including the effects of ambipolar and Hall drift. While both describe transport processes that arise from the multifluid equations and are therefore of hyperbolic nature, they are often incorporated in models as a diffusive, i.e. parabolic process. While the formulation as such is easy to include in standard MHD models, the resulting diffusive time-step constraints do require often a computationally more expensive implicit treatment or super-time-stepping approaches. In this paper we discuss an implementation that retains the hyperbolic nature and allows for an explicit integration with small computational overhead. In the case of ambipolar drift, this formulation arises naturally by simply retaining a time derivative of the drift velocity that is typically omitted. This alone leads to time-step constraints that are comparable to the native MHD time-step constraint for a solar setup including the region from photosphere to lower solar corona. We discuss an accelerated treatment that can further reduce time-step constraints if necessary. In the case of Hall drift we propose a hyperbolic formulation that is numerically similar to that for the ambipolar drift and we show that the combination of both can be applied to simulations of the solar chromosphere at minimal computational expense. △ Less

Submitted 26 November, 2021; originally announced November 2021.

Comments: 21 pages, 13 figures, accepted for publication in ApJ

Solar Atmosphere Radiative Transfer Model Comparison based on 3D MHD simulations

Authors: M. Haberreiter, S. Criscuoli, M. Rempel, T. M. D. Pereira

Abstract: The reconstruction of the solar spectral irradiance (SSI) on various time scales is essential for the understanding of the Earth's climate response to the SSI variability. The driver of the SSI variability is understood to be the intensity contrast of magnetic features present on the Sun with respect to the largely non-magnetic quiet Sun. However, different spectral synthesis codes lead to divergi… ▽ More The reconstruction of the solar spectral irradiance (SSI) on various time scales is essential for the understanding of the Earth's climate response to the SSI variability. The driver of the SSI variability is understood to be the intensity contrast of magnetic features present on the Sun with respect to the largely non-magnetic quiet Sun. However, different spectral synthesis codes lead to diverging projections of SSI variability. In this study we compare three different radiative transfer codes and carry out a detailed analysis of their performance. We perform the spectral synthesis at the continuum wavelength of 665 nm with the Code for Solar Irradiance (COSI), and the Rybicki-Hummer (RH), and Max Planck University of Chicago Radiative MHD (MURaM) codes for three 3D MHD simulations snapshots, a non-magnetic case, and MHD simulations with 100 G, and 200 G magnetic field strength. We determine the intensity distributions, the intensity differences and ratios for the spectral synthesis codes. We identify that the largest discrepancies originate in the intergranular lanes where the most field concentration occurs. Overall, the applied radiative transfer codes give consistent intensity distributions. Also, the intensity variation as a function of magnetic field strength for the particular 100 G and 200 G snapshots agree within the 2-3% range. △ Less

Submitted 6 September, 2021; originally announced September 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 653, A161 (2021)

Probing the Physics of the Solar Atmosphere with the Multi-slit Solar Explorer (MUSE): II. Flares and Eruptions

Authors: Mark C. M. Cheung, Juan Martínez-Sykora, Paola Testa, Bart De Pontieu, Georgios Chintzoglou, Matthias Rempel, Vanessa Polito, Graham S. Kerr, Katharine K. Reeves, Lyndsay Fletcher, Meng Jin, Daniel Nóbrega-Siverio, Sanja Danilovic, Patrick Antolin, Joel Allred, Viggo Hansteen, Ignacio Ugarte-Urra, Edward DeLuca, Dana Longcope, Shinsuke Takasao, Marc DeRosa, Paul Boerner, Sarah Jaeggli, Nariaki Nitta, Adrian Daw , et al. (3 additional authors not shown)

Abstract: Current state-of-the-art spectrographs cannot resolve the fundamental spatial (sub-arcseconds) and temporal scales (less than a few tens of seconds) of the coronal dynamics of solar flares and eruptive phenomena. The highest resolution coronal data to date are based on imaging, which is blind to many of the processes that drive coronal energetics and dynamics. As shown by IRIS for the low solar at… ▽ More Current state-of-the-art spectrographs cannot resolve the fundamental spatial (sub-arcseconds) and temporal scales (less than a few tens of seconds) of the coronal dynamics of solar flares and eruptive phenomena. The highest resolution coronal data to date are based on imaging, which is blind to many of the processes that drive coronal energetics and dynamics. As shown by IRIS for the low solar atmosphere, we need high-resolution spectroscopic measurements with simultaneous imaging to understand the dominant processes. In this paper: (1) we introduce the Multi-slit Solar Explorer (MUSE), a spaceborne observatory to fill this observational gap by providing high-cadence (<20 s), sub-arcsecond resolution spectroscopic rasters over an active region size of the solar transition region and corona; (2) using advanced numerical models, we demonstrate the unique diagnostic capabilities of MUSE for exploring solar coronal dynamics, and for constraining and discriminating models of solar flares and eruptions; (3) we discuss the key contributions MUSE would make in addressing the science objectives of the Next Generation Solar Physics Mission (NGSPM), and how MUSE, the high-throughput EUV Solar Telescope (EUVST) and the Daniel K Inouye Solar Telescope (and other ground-based observatories) can operate as a distributed implementation of the NGSPM. This is a companion paper to De Pontieu et al. (2021; arXiv:2106.15584), which focuses on investigating coronal heating with MUSE. △ Less

Submitted 30 June, 2021; v1 submitted 29 June, 2021; originally announced June 2021.

Comments: 42 pages, 22 figures, submitted to ApJ

Probing the physics of the solar atmosphere with the Multi-slit Solar Explorer (MUSE): I. Coronal Heating

Authors: Bart De Pontieu, Paola Testa, Juan Martinez-Sykora, Patrick Antolin, Konstantinos Karampelas, Viggo Hansteen, Matthias Rempel, Mark C. M. Cheung, Fabio Reale, Sanja Danilovic, Paolo Pagano, Vanessa Polito, Ineke De Moortel, Daniel Nobrega-Siverio, Tom Van Doorsselaere, Antonino Petralia, Mahboubeh Asgari-Targhi, Paul Boerner, Mats Carlsson, Georgios Chintzoglou, Adrian Daw, Ed DeLuca, Leon Golub, Takuma Matsumoto, Ignacio Ugarte-Urra , et al. (2 additional authors not shown)

Abstract: The Multi-slit Solar Explorer (MUSE) is a proposed NASA MIDEX mission, currently in Phase A, composed of a multi-slit EUV spectrograph (in three narrow spectral bands centered around 171A, 284A, and 108A) and an EUV context imager (in two narrow passbands around 195A and 304A). MUSE will provide unprecedented spectral and imaging diagnostics of the solar corona at high spatial (<0.5 arcsec), and t… ▽ More The Multi-slit Solar Explorer (MUSE) is a proposed NASA MIDEX mission, currently in Phase A, composed of a multi-slit EUV spectrograph (in three narrow spectral bands centered around 171A, 284A, and 108A) and an EUV context imager (in two narrow passbands around 195A and 304A). MUSE will provide unprecedented spectral and imaging diagnostics of the solar corona at high spatial (<0.5 arcsec), and temporal resolution (down to ~0.5s) thanks to its innovative multi-slit design. By obtaining spectra in 4 bright EUV lines (Fe IX 171A , Fe XV 284A, Fe XIX-Fe XXI 108A) covering a wide range of transition region and coronal temperatures along 37 slits simultaneously, MUSE will for the first time be able to "freeze" (at a cadence as short as 10 seconds) with a spectroscopic raster the evolution of the dynamic coronal plasma over a wide range of scales: from the spatial scales on which energy is released (~0.5 arcsec) to the large-scale often active-region size (170 arcsec x 170 arcsec) atmospheric response. We use advanced numerical modeling to showcase how MUSE will constrain the properties of the solar atmosphere on the spatio-temporal scales (~0.5 arcsec, ~20 seconds) and large field-of-view on which various state-of-the-art models of the physical processes that drive coronal heating, solar flares and coronal mass ejections (CMEs) make distinguishing and testable predictions. We describe how the synergy between MUSE, the single-slit, high-resolution Solar-C EUVST spectrograph, and ground-based observatories (DKIST and others) can address how the solar atmosphere is energized, and the critical role MUSE plays because of the multi-scale nature of the physical processes involved. In this first paper, we focus on how comparisons between MUSE observations and theoretical models will significantly further our understanding of coronal heating mechanisms. △ Less

Submitted 29 June, 2021; originally announced June 2021.

Comments: 46 pages, 22 figures, submitted to ApJ

The Coronal Veil

Authors: A. Malanushenko, M. C. M. Cheung, C. E. DeForest, J. A. Klimchuk, M. Rempel

Abstract: Coronal loops, seen in solar coronal images, are believed to represent emission from magnetic flux tubes with compact cross-sections. We examine the 3D structure of plasma above an active region in a radiative magnetohydrodynamic simulation to locate volume counterparts for coronal loops. In many cases, a loop cannot be linked to an individual thin strand in the volume. While many thin loops are p… ▽ More Coronal loops, seen in solar coronal images, are believed to represent emission from magnetic flux tubes with compact cross-sections. We examine the 3D structure of plasma above an active region in a radiative magnetohydrodynamic simulation to locate volume counterparts for coronal loops. In many cases, a loop cannot be linked to an individual thin strand in the volume. While many thin loops are present in the synthetic images, the bright structures in the volume are fewer, and of complex shape. We demonstrate that this complexity can form impressions of thin bright loops, even in the absence of thin bright plasma strands. We demonstrate the difficulty of discerning from observations whether a particular loop corresponds to a strand in the volume, or a projection artifact. We demonstrate how apparently isolated loops could deceive observers, even when observations from multiple viewing angles are available. While we base our analysis on a simulation, the main findings are independent from a particular simulation setup and illustrate the intrinsic complexity involved in interpreting observations resulting from line-of-sight integration in an optically thin plasma. We propose alternative interpretation for strands seen in EUV images of the corona. The "coronal veil" hypothesis is mathematically more generic, and naturally explains properties of loops that are difficult to address otherwise -- such as their constant cross section and anomalously high density scale height. We challenge the paradigm of coronal loops as thin magnetic flux tubes, offering new understanding of solar corona and, by extension, of other magnetically confined bright, hot plasmas. △ Less

Submitted 29 November, 2021; v1 submitted 28 June, 2021; originally announced June 2021.

Comments: ApJ (in press)

A Comprehensive Radiative Magnetohydrodynamics Simulation of Active Region Scale Flux Emergence from the Convection Zone to the Corona

Authors: Feng Chen, Matthias Rempel, Yuhong Fan

Abstract: We present a comprehensive radiative magnetohydrodynamic simulation of the quiet Sun and large solar active regions. The 197 Mm wide simulation domain spans from 18 (10) Mm beneath the photosphere to 113 Mm in the solar corona. Radiative transfer assuming local thermal equilibrium, optically-thin radiative losses, and anisotropic conduction transport provide the necessary realism for synthesizing… ▽ More We present a comprehensive radiative magnetohydrodynamic simulation of the quiet Sun and large solar active regions. The 197 Mm wide simulation domain spans from 18 (10) Mm beneath the photosphere to 113 Mm in the solar corona. Radiative transfer assuming local thermal equilibrium, optically-thin radiative losses, and anisotropic conduction transport provide the necessary realism for synthesizing observables to compare with remote sensing observations of the photosphere and corona. This model self-consistently reproduces observed features of the quiet Sun, emerging and developed active regions, and solar flares up to M class. Here, we report an overview of the first results. The surface magnetoconvection yields an upward Poynting flux that is dissipated in the corona and heats the plasma to over one million K. The quiescent corona also presents ubiquitous propagating waves, jets, and bright points with sizes down to 2 Mm. Magnetic flux bundles emerge into the photosphere and give rise to strong and complex active regions with over $10^{23}$ Mx magnetic flux. The coronal free magnetic energy, which is approximately 18\% of the total magnetic energy, accumulates to approximately $10^{33}$ erg. The coronal magnetic field is clearly non-force-free, as the Lorentz force needs to balance the pressure force and viscous stress as well as drive magnetic field evolution. The emission measure from $\log_{10}T{=}4.5$ to $\log_{10}T{>}7$ provides a comprehensive view of the active region corona, such as coronal loops of various lengths and temperatures, mass circulation by evaporation and condensation, and eruptions from jets to large-scale mass ejections. △ Less

Submitted 5 September, 2022; v1 submitted 26 June, 2021; originally announced June 2021.

Comments: 42 Pages, 19 figures, accepted for publication in ApJ

Measuring the magnetic origins of solar flares, CMEs and Space Weather

Authors: Philip Judge, Matthias Rempel, Rana Ezzeddine, Lucia Kleint, Ricky Egeland, Svetlana Berdyugina, Thomas Berger, Joan Burkepile, Rebecca Centeno, Giuliana de Toma, Mausumi Dikpati, Yuhong Fan, Holly Gilbert, Daniela Lacatus

Abstract: We take a broad look at the problem of identifying the magnetic solar causes of space weather. With the lackluster performance of extrapolations based upon magnetic field measurements in the photosphere, we identify a region in the near UV part of the spectrum as optimal for studying the development of magnetic free energy over active regions. Using data from SORCE, Hubble Space Telescope, and SKY… ▽ More We take a broad look at the problem of identifying the magnetic solar causes of space weather. With the lackluster performance of extrapolations based upon magnetic field measurements in the photosphere, we identify a region in the near UV part of the spectrum as optimal for studying the development of magnetic free energy over active regions. Using data from SORCE, Hubble Space Telescope, and SKYLAB, along with 1D computations of the near-UV (NUV) spectrum and numerical experiments based on the MURaM radiation-MHD and HanleRT radiative transfer codes, we address multiple challenges. These challenges are best met through a combination of near UV lines of bright \ion{Mg}{2}, and lines of \ion{Fe}{2} and \ion{Fe}{1} (mostly within the $4s-4p$ transition array) which form in the chromosphere up to $2\times10^4$ K. Both Hanle and Zeeman effects can in principle be used to derive vector magnetic fields. However, for any given spectral line the $τ=1$ surfaces are generally geometrically corrugated owing to fine structure such as fibrils and spicules. By using multiple spectral lines spanning different optical depths, magnetic fields across nearly-horizontal surfaces can be inferred in regions of low plasma $β$, from which free energies, magnetic topology and other quantities can be derived. Based upon the recently-reported successful suborbital space measurements of magnetic fields with the CLASP2 instrument, we argue that a modest space-borne telescope will be able to make significant advances in the attempts to predict solar eruptions. Difficulties associated with blended lines are shown to be minor in an Appendix. △ Less

Submitted 16 June, 2021; v1 submitted 14 June, 2021; originally announced June 2021.

Comments: Accepted in ApJ, V2 fixed author spelling error

On the (in)stability of sunspots

Authors: Hanna Strecker, Wolfgang Schmidt, Rolf Schlichenmaier, Matthias Rempel

Abstract: The stability of sunspots is one of the long-standing unsolved puzzles in the field of solar magnetism. We study the effects that destabilise and stabilise the flux tube of a simulated sunspot in the upper convection zone. The depth-varying effects of fluting instability, buoyancy forces, and timescales on the flux tube are analysed. The simulation was calculated with the MURaM code. The domain ha… ▽ More The stability of sunspots is one of the long-standing unsolved puzzles in the field of solar magnetism. We study the effects that destabilise and stabilise the flux tube of a simulated sunspot in the upper convection zone. The depth-varying effects of fluting instability, buoyancy forces, and timescales on the flux tube are analysed. The simulation was calculated with the MURaM code. The domain has a lateral extension of 98 Mm x 98 Mm and extends almost 18 Mm below the solar surface. The analysed data set of 30 hours shows a stable sunspot at the solar surface. We studied the evolution of the flux tube at horizontal layers by means of the relative change in perimeter and area with a linear stability analysis. We find a corrugation along the perimeter of the flux tube that proceeds fastest at a depth of about 8 Mm below the surface. Towards the surface and towards deeper layers, the decrease in compactness is damped. From the stability analysis, we find that above a depth of 2 Mm, the sunspot is stabilised by buoyancy forces. The spot is least stable at a depth of about 3 Mm because of fluting instability. The stability of the sunspot at the surface affects the behaviour of the field lines in deeper layers by magnetic tension. Therefore the fluting instability is damped at depths of about 3 Mm, and the decrease in compactness is strongest at a depth of about 8 Mm. The more vertical orientation of the magnetic field and the longer convective timescale slow down the corrugation process in layers deeper than 10 Mm. The formation of large intrusions of field-free plasma below the surface destabilises the sunspot, and eventually lead to the disruption and decay of the sunspot. This process is not visible at the surface, where the sunspot is stabilised by buoyancy forces. The onset of sunspot decay occurs in deeper layers, while the sunspot still appears stable in the photosphere. △ Less

Submitted 21 March, 2021; originally announced March 2021.

Comments: 12 pages, 14 figures; appendix: 3 pages, 2 figures, 1 table; accepted for publication in A&A

Journal ref: A&A 649, A123 (2021)

The Dimmest State of the Sun

Authors: K. L. Yeo, S. K. Solanki, N. A. Krivova, M. Rempel, L. S. Anusha, A. I. Shapiro, R. V. Tagirov, V. Witzke

Abstract: How the solar electromagnetic energy entering the Earth's atmosphere varied since pre-industrial times is an important consideration in the climate change debate. Detrimental to this debate, estimates of the change in total solar irradiance (TSI) since the Maunder minimum, an extended period of weak solar activity preceding the industrial revolution, differ markedly, ranging from a drop of 0.75 Wm… ▽ More How the solar electromagnetic energy entering the Earth's atmosphere varied since pre-industrial times is an important consideration in the climate change debate. Detrimental to this debate, estimates of the change in total solar irradiance (TSI) since the Maunder minimum, an extended period of weak solar activity preceding the industrial revolution, differ markedly, ranging from a drop of 0.75 Wm-2 to a rise of 6.3 Wm-2. Consequently, the exact contribution by solar forcing to the rise in global temperatures over the past centuries remains inconclusive. Adopting a novel approach based on state-of-the-art solar imagery and numerical simulations, we establish the TSI level of the Sun when it is in its least-active state to be 2.0 +/- 0.7 Wm-2 below the 2019 level. This means TSI could not have risen since the Maunder minimum by more than this amount, thus restricting the possible role of solar forcing in global warming. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: Supporting information: https://agupubs.onlinelibrary.wiley.com/action/downloadSupplement?doi=10.1029%2F2020GL090243&file=grl61316-sup-0001-Text_SI-S01.pdf

Journal ref: Geophysical Research Letters, 47, e2020GL090243 (2020)

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

Authors: H. Ji, J. Karpen, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, A. Bhattacharjee, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, B. Chen, L. -J. Chen, Y. Chen, A. Chien, L. Comisso, D. Craig, J. Dahlin, W. Daughton, E. DeLuca , et al. (83 additional authors not shown)

Abstract: Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagen… ▽ More Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 4 pages (including a title page), white paper submitted to Helio2050 workshop at https://www.hou.usra.edu/meetings/helio2050/. arXiv admin note: substantial text overlap with arXiv:2004.00079

Critical Science Plan for the Daniel K. Inouye Solar Telescope (DKIST)

Authors: Mark P. Rast, Nazaret Bello González, Luis Bellot Rubio, Wenda Cao, Gianna Cauzzi, Edward DeLuca, Bart De Pontieu, Lyndsay Fletcher, Sarah E. Gibson, Philip G. Judge, Yukio Katsukawa, Maria D. Kazachenko, Elena Khomenko, Enrico Landi, Valentin Martínez Pillet, Gordon J. D. Petrie, Jiong Qiu, Laurel A. Rachmeler, Matthias Rempel, Wolfgang Schmidt, Eamon Scullion, Xudong Sun, Brian T. Welsch, Vincenzo Andretta, Patrick Antolin , et al. (62 additional authors not shown)

Abstract: The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With… ▽ More The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With this Critical Science Plan (CSP) we attempt to anticipate some of what those capabilities will enable, providing a snapshot of some of the scientific pursuits that the Daniel K. Inouye Solar Telescope hopes to engage as start-of-operations nears. The work builds on the combined contributions of the DKIST Science Working Group (SWG) and CSP Community members, who generously shared their experiences, plans, knowledge and dreams. Discussion is primarily focused on those issues to which DKIST will uniquely contribute. △ Less

Submitted 20 August, 2020; v1 submitted 18 August, 2020; originally announced August 2020.

Acoustic-gravity wave propagation characteristics in 3D radiation hydrodynamic simulations of the solar atmosphere

Authors: B. Fleck, M. Carlsson, E. Khomenko, M. Rempel, O. Steiner, G. Vigeesh

Abstract: There has been tremendous progress in the degree of realism of three-dimensional radiation magneto-hydrodynamic simulations of the solar atmosphere in the past decades. Four of the most frequently used numerical codes are Bifrost, CO5BOLD, MANCHA3D, and MURaM. Here we test and compare the wave propagation characteristics in model runs from these four codes by measuring the dispersion relation of a… ▽ More There has been tremendous progress in the degree of realism of three-dimensional radiation magneto-hydrodynamic simulations of the solar atmosphere in the past decades. Four of the most frequently used numerical codes are Bifrost, CO5BOLD, MANCHA3D, and MURaM. Here we test and compare the wave propagation characteristics in model runs from these four codes by measuring the dispersion relation of acoustic-gravity waves at various heights. We find considerable differences between the various models. The height dependence of wave power, in particular of high-frequency waves, varies by up to two orders of magnitude between the models, and the phase difference spectra of several models show unexpected features, including $\pm180^\circ$ phase jumps. △ Less

Submitted 11 July, 2020; originally announced July 2020.

Comments: 19 pages, 15 figures

A distinct magnetic property of the inner penumbral boundary III. Analysis of simulated sunspots

Authors: Jan Jurcak, Markus Schmassmann, Matthias Rempel, Nazaret Bello Gonzalez, Rolf Schlichenmaier

Abstract: The analyses of sunspot observations revealed a fundamental magnetic property of the umbral boundary, the invariance of the vertical component of the magnetic field. We aim to analyse the magnetic properties of the umbra-penumbra boundary in simulated sunspots and thus assess their similarity to observed sunspots. Also, we aim to investigate the role of plasma $β$ and the ratio of kinetic to magne… ▽ More The analyses of sunspot observations revealed a fundamental magnetic property of the umbral boundary, the invariance of the vertical component of the magnetic field. We aim to analyse the magnetic properties of the umbra-penumbra boundary in simulated sunspots and thus assess their similarity to observed sunspots. Also, we aim to investigate the role of plasma $β$ and the ratio of kinetic to magnetic energy in simulated sunspots on the convective motions. We use a set of non-grey simulation runs of sunspots with the MURaM code. These data are used to synthesise the Stokes profiles that are then degraded to the Hinode spectropolarimeter-like observations. Then, the data are treated like real Hinode observations of a sunspot and magnetic properties at the umbral boundaries are determined. Simulations with potential field extrapolation produce a realistic magnetic field configuration on their umbral boundaries. Two simulations with potential field upper boundary, but different subsurface magnetic field structures, differ significantly in the extent of their penumbrae. Increasing the penumbra width by forcing more horizontal magnetic fields at the upper boundary results in magnetic properties that are not consistent with observations. This implies that the size of the penumbra is given by the subsurface structure of the magnetic field. None of the sunspot simulations is consistent with observed properties of the magnetic field and direction of the Evershed flow at the same time. Strong outward directed Evershed flows are only found in setups with artificially enhanced horizontal component of the magnetic field at the top boundary that are not consistent with the observed magnetic field properties at the UP boundary. We want to stress out that the `photospheric' boundary of simulated sunspots is defined by a magnetic field strength of equipartition field value. △ Less

Submitted 8 April, 2020; originally announced April 2020.

Journal ref: A&A 638, A28 (2020)

On the contribution of quiet Sun magnetism to solar irradiance variations: Constraints on quiet Sun variability and grand minimum scenarios

Authors: M. Rempel

Abstract: While the quiet Sun magnetic field shows only little variation with the solar cycle, long-term variations cannot be completely ruled out from first principles. We investigate the potential effect of quiet Sun magnetism on spectral solar irradiance through a series of small-scale dynamo simulations with zero vertical flux imbalance ($\langle B_z\rangle=0$) and varying levels of small-scale magnetic… ▽ More While the quiet Sun magnetic field shows only little variation with the solar cycle, long-term variations cannot be completely ruled out from first principles. We investigate the potential effect of quiet Sun magnetism on spectral solar irradiance through a series of small-scale dynamo simulations with zero vertical flux imbalance ($\langle B_z\rangle=0$) and varying levels of small-scale magnetic field strength, and one weak network case with an additional flux imbalance corresponding to a flux density of $\langle B_z\rangle=100$ G. From these setups we compute the dependence of the outgoing radiative energy flux on the mean vertical magnetic field strength in the photosphere at continuum optical depth $τ=1$ ($\langle \vert B_z\vert\rangle_{τ=1}$). We find that a quiet Sun setup with a mean vertical field strength of $\langle \vert B_z\vert\rangle_{τ=1}=69$ G is about $0.6~\%$ brighter than a non-magnetic reference case. We find a linear dependence of the outgoing radiative energy flux on the mean field strength $\langle \vert B_z\vert\rangle_{τ=1}$ with a relative slope of $1.4\cdot 10^{-4}$ G$^{-1}$. With this sensitivity, only a moderate change of the quiet Sun field strength by $10\%$ would lead to a total solar irradiance variation comparable to the observed solar cycle variation. While this does provide strong indirect constraints on possible quiet Sun variations during a regular solar cycle, it also emphasizes that potential variability over longer time scales could make a significant contribution to longer-term solar irradiance variations. △ Less

Submitted 3 April, 2020; originally announced April 2020.

Comments: 15 pages, 10 figures, accepted for publication in ApJ

Journal ref: ApJ, Volume 894, Issue 2, id.140, May 2020

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena throughout the Universe

Authors: H. Ji, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, L. -J. Chen, Y. Chen, A. Chien, D. Craig, J. Dahlin, W. Daughton, E. DeLuca, C. F. Dong, S. Dorfman, J. Drake, F. Ebrahimi , et al. (75 additional authors not shown)

Abstract: This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process. This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process. △ Less

Submitted 31 March, 2020; originally announced April 2020.

Comments: 6 pages, 1 figure, white paper submitted to both Plasma 2020 and Astro 2020 Decadal Surveys

Combining magneto-hydrostatic constraints with Stokes profiles inversions

Authors: J. M. Borrero, A. Pastor Yabar, M. Rempel, B. Ruiz Cobo

Abstract: Inversion codes for the polarized radiative transfer equation can be used to infer the temperature $T$, line-of-sight velocity $v_{\rm los}$, and magnetic field $\rm{\bf B}$ as a function of the continuum optical-depth $τ_{\rm c}$. However, they do not directly provide the gas pressure $P_{\rm g}$ or density $ρ$. In order to obtain these latter parameters, inversion codes rely instead on the assum… ▽ More Inversion codes for the polarized radiative transfer equation can be used to infer the temperature $T$, line-of-sight velocity $v_{\rm los}$, and magnetic field $\rm{\bf B}$ as a function of the continuum optical-depth $τ_{\rm c}$. However, they do not directly provide the gas pressure $P_{\rm g}$ or density $ρ$. In order to obtain these latter parameters, inversion codes rely instead on the assumption of hydrostatic equilibrium (HE) in addition to the equation of state (EOS). Unfortunately, the assumption of HE is rather unrealistic across magnetic field lines. This is because the role of the Lorentz force, among other factors, is neglected. This translates into an inaccurate conversion from optical depth $τ_{\rm c}$ to geometrical height $z$. We aim at improving this conversion via the application of magneto-hydrostatic (MHS) equilibrium instead of HE. We develop a method to solve the momentum equation under MHS equilibrium (i.e., taking the Lorentz force into account) in three dimensions. The method is based on the solution of a Poisson-like equation. Considering the gas pressure $P_{\rm g}$ and density $ρ$ from three-dimensional magneto-hydrodynamic (MHD) simulations of sunspots as a benchmark, we compare the results from the application of HE and MHS equilibrium. We find that HE retrieves the gas pressure and density within an order of magnitude of the MHD values in only about 47 \% of the domain. This translates into an error of about $160-200$ km in the determination of the $z-τ_{\rm c}$ conversion. On the other hand, the application of MHS equilibrium allows determination of $P_{\rm g}$ and $ρ$ within an order of magnitude in 84 \% of the domain. In this latter case, the $z-τ_{\rm c}$ conversion is obtained with an accuracy of $30-70$ km. △ Less

Submitted 30 October, 2019; originally announced October 2019.

Comments: Accepted for publication in Astronomy \& Astrophysics. 12 pages, 10 figures

Journal ref: A&A 632, A111 (2019)

Principles Of Heliophysics: a textbook on the universal processes behind planetary habitability

Authors: Karel Schrijver, Fran Bagenal, Tim Bastian, Juerg Beer, Mario Bisi, Tom Bogdan, Steve Bougher, David Boteler, Dave Brain, Guy Brasseur, Don Brownlee, Paul Charbonneau, Ofer Cohen, Uli Christensen, Tom Crowley, Debrah Fischer, Terry Forbes, Tim Fuller-Rowell, Marina Galand, Joe Giacalone, George Gloeckler, Jack Gosling, Janet Green, Nick Gross, Steve Guetersloh , et al. (37 additional authors not shown)

Abstract: Heliophysics is the system science of the physical connections between the Sun and the solar system. As the physics of the local cosmos, it embraces space weather and planetary habitability. The wider view of comparative heliophysics forms a template for conditions in exoplanetary systems and provides a view over time of the aging Sun and its magnetic activity, of the heliosphere in different sett… ▽ More Heliophysics is the system science of the physical connections between the Sun and the solar system. As the physics of the local cosmos, it embraces space weather and planetary habitability. The wider view of comparative heliophysics forms a template for conditions in exoplanetary systems and provides a view over time of the aging Sun and its magnetic activity, of the heliosphere in different settings of the interstellar medium and subject to stellar impacts, of the space physics over evolving planetary dynamos, and of the long-term influence on planetary atmospheres by stellar radiation and wind. Based on a series of NASA-funded summer schools for early-career researchers, this textbook is intended for students in physical sciences in later years of their university training and for beginning graduate students in fields of solar, stellar, (exo-)planetary, and planetary-system sciences. The book emphasizes universal processes from a perspective that draws attention to what provides Earth (and similar (exo-)planets) with a relatively stable setting in which life as we know it could thrive. The text includes 200 "Activities" in the form of exercises, explorations, literature readings, "what if" challenges, and group discussion topics; many of the Activities provide additional information complementing the main text. Solutions and discussions are included in an Appendix for a selection of the exercises. △ Less

Submitted 27 March, 2024; v1 submitted 30 October, 2019; originally announced October 2019.

Comments: 429 pages, 121 figures, and 200 "activities" in the form of problems, exercises, literature readings, and "what if" challenges; V2.0: activities modified; added 'solutions' to selected activities; clarified text and corrected typos; added subject index; V2.1 added Activity hyperlinks and corrected typos (see 'Version history' in the manuscript)

Superstrong photospheric magnetic fields in sunspot penumbrae

Authors: A. Siu-Tapia, A. Lagg, M. van Noort, M. Rempel, S. K. Solanki

Abstract: Recently, there have been some reports of unusually strong photospheric magnetic fields (which can reach values of over 7 kG) inferred from Hinode SOT/SP sunspot observations within penumbral regions. These superstrong penumbral fields are even larger than the strongest umbral fields on record and appear to be associated with supersonic downflows. The finding of such fields has been controversial… ▽ More Recently, there have been some reports of unusually strong photospheric magnetic fields (which can reach values of over 7 kG) inferred from Hinode SOT/SP sunspot observations within penumbral regions. These superstrong penumbral fields are even larger than the strongest umbral fields on record and appear to be associated with supersonic downflows. The finding of such fields has been controversial since they seem to show up only when spatially coupled inversions are performed. Here, we investigate and discuss the reliability of those findings by studying in detail observed spectra associated with particularly strong magnetic fields at the inner edge of the penumbra of active region 10930. We apply classical diagnostic methods and various inversions with different model atmospheres to the observed Stokes profiles in two selected pixels with superstrong magnetic fields, and compare the results with a magnetohydrodynamic simulation of a sunspot whose penumbra contains localized regions with strong fields (nearly 5 kG at $τ=1$) associated with supersonic downflows... △ Less

Submitted 30 September, 2019; originally announced September 2019.

Comments: 19 pages, 14 figures, 3 tables. Now accepted for publication in Section 9. The Sun of Astronomy and Astrophysics. Date of acceptance is 23/09/2019

Journal ref: A&A 631, A99 (2019)

Opposite polarity magnetic field and convective downflows in a simulated sunspot penumbra

Authors: Lokesh Bharti, Matthias Rempel

Abstract: Recent numerical simulations and observations of sunspots show a significant amount of opposite polarity magnetic field within the sunspot penumbra. Most of the opposite polarity field is associated with convective downflows. We present an analysis of 3D MHD simulations through forward modeling of synthetic Stokes profiles of the Fe\sci 6301.5 Å~ and Fe\sci 6302.5 Å~ lines). The synthetic Stokes p… ▽ More Recent numerical simulations and observations of sunspots show a significant amount of opposite polarity magnetic field within the sunspot penumbra. Most of the opposite polarity field is associated with convective downflows. We present an analysis of 3D MHD simulations through forward modeling of synthetic Stokes profiles of the Fe\sci 6301.5 Å~ and Fe\sci 6302.5 Å~ lines). The synthetic Stokes profiles are spatially and spectrally degraded considering typical instrument properties. Line bisector shifts of the Fe\sci 6301.5 Å~ line are used to determine line-of-sight velocities. Far wing magnetograms are constructed from the Stokes V profiles of the Fe\sci 6302.5 Å~ line. While we find an overall good agreement between observations and simulations, the fraction of opposite polarity magnetic field, the downflow filling factor and the opposite polarity-downflow association are strongly affected by spatial smearing and presence of strong gradients in the line-of-sight magnetic field and velocity. A significant fraction of opposite polarity magnetic field and downflows are hidden in the observations due to typical instrumental noise. Comparing simulations that differ by more than a factor of two in grid spacing we find that these quantities are robust within the simulations. △ Less

Submitted 18 August, 2019; originally announced August 2019.

Comments: 10 pages, 4 figures, Accepted for publication in Astrophysical Journal

Three-dimensional modeling of chromospheric spectral lines in a simulated active region

Authors: Johan P. Bjørgen, Jorrit Leenaarts, Matthias Rempel, Mark C. M. Cheung, Sanja Danilovic, Jaime de la Cruz Rodríguez, Andrii V. Sukhorukov

Abstract: Because of the complex physics that governs the formation of chromospheric lines, interpretation of solar chromospheric observations is difficult. The origin and characteristics of many chromospheric features are, because of this, unresolved. We focus here on studying two prominent features: long fibrils and flare ribbons. To model them, we use a 3D MHD simulation of an active region which self-co… ▽ More Because of the complex physics that governs the formation of chromospheric lines, interpretation of solar chromospheric observations is difficult. The origin and characteristics of many chromospheric features are, because of this, unresolved. We focus here on studying two prominent features: long fibrils and flare ribbons. To model them, we use a 3D MHD simulation of an active region which self-consistently reproduces both of them. We model the H$α$, Mg II k, Ca II K, and Ca II 8542 Å lines using the 3D non-LTE radiative transfer code Multi3D. This simulation reproduces long fibrils that span between the opposite-polarity sunspots and go up to 4 Mm in height. They can be traced in all lines due to density corrugation. Opposite to previous studies, H$α$, Mg II h&k, and Ca II H&K, are formed at similar height in this model. Magnetic field lines are aligned with the H$α$ fibrils, but the latter holds to a lesser extent for the Ca II 8542 Å line. The simulation shows structures in the H$α$ line core that look like flare ribbons. The emission in the ribbons is caused by a dense chromosphere and a transition region at high column mass. The ribbons are visible in all chromospheric lines, but least prominent in Ca II 8542 Å line. In some pixels, broad asymmetric profiles with a single emission peak are produced, similar to the profiles observed in flare ribbons. They are caused by a deep onset of the chromospheric temperature rise and large velocity gradients. The simulation produces long fibrils similar to what is seen in observations. It also produces structures similar to flare ribbons despite the lack of non-thermal electrons in the simulation. The latter suggests that thermal conduction might be a significant agent in transporting flare energy to the chromosphere in addition to non-thermal electrons. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: 17 pages, 15 figures, accepted for publication by A&A

Journal ref: A&A 631, A33 (2019)

Nonlinear dynamo models using quasi-biennial oscillations constrained by sunspot area data

Authors: Fadil Inceoglu, Rosaria Simoniello, Rainer Arlt, Matthias Rempel

Abstract: Contex: Solar magnetic activity exhibits variations with periods between 1.5--4 years, the so-called quasi-biennial oscillations (QBOs), in addition to the well-known 11-year Schwabe cycles. Solar dynamo is thought to be the responsible mechanism for generation of the QBOs. Aims: In this work, we analyse sunspot areas to investigate the spatial and temporal behaviour of the QBO signal and study… ▽ More Contex: Solar magnetic activity exhibits variations with periods between 1.5--4 years, the so-called quasi-biennial oscillations (QBOs), in addition to the well-known 11-year Schwabe cycles. Solar dynamo is thought to be the responsible mechanism for generation of the QBOs. Aims: In this work, we analyse sunspot areas to investigate the spatial and temporal behaviour of the QBO signal and study the responsible physical mechanisms using simulations from fully nonlinear mean-field flux-transport dynamos. Methods: We investigated the behaviour of the QBOs in the sunspot area data in full disk, and northern and southern hemispheres, using wavelet and Fourier analyses. We also ran solar dynamos with two different approaches to generating a poloidal field from an existing toroidal field, Babcock-Leighton and turbulent $α$ mechanisms. We then studied the simulated magnetic field strengths as well as meridional circulation and differential rotation rates using the same methods. Results: The results from the sunspot areas show that the QBOs are present in the full disk and hemispheric sunspot areas and they show slightly different spatial and temporal behaviours, indicating a slightly decoupled solar hemispheres. The QBO signal is generally intermittent and in-phase with the sunspot area data, surfacing when the solar activity is in maximum. The results from the BL-dynamos showed that they are neither capable of generating the slightly decoupled behaviour of solar hemispheres nor can they generate QBO-like signals. The turbulent $α$-dynamos, on the other hand, generated decoupled hemispheres and some QBO-like shorter cycles. Conclusions: In conclusion, our simulations show that the turbulent $α$-dynamos with the Lorentz force seems more efficient in generating the observed temporal and spatial behaviour of the QBO signal compared with those from the BL-dynamos. △ Less

Submitted 7 April, 2019; originally announced April 2019.

Reversed dynamo at small scales and large magnetic Prandtl number

Authors: Axel Brandenburg, Matthias Rempel

Abstract: We show that at large magnetic Prandtl numbers, the Lorentz force does work on the flow at small scales and drives fluid motions, whose energy is dissipated viscously. This situation is opposite to that in a normal dynamo, where the flow does work against the Lorentz force. We compute the spectral conversion rates between kinetic and magnetic energies for several magnetic Prandtl numbers and show… ▽ More We show that at large magnetic Prandtl numbers, the Lorentz force does work on the flow at small scales and drives fluid motions, whose energy is dissipated viscously. This situation is opposite to that in a normal dynamo, where the flow does work against the Lorentz force. We compute the spectral conversion rates between kinetic and magnetic energies for several magnetic Prandtl numbers and show that normal (forward) dynamo action occurs on large scales over a progressively narrower range of wavenumbers as the magnetic Prandtl number is increased. At higher wavenumbers, reversed dynamo action occurs, i.e., magnetic energy is converted back into kinetic energy at small scales. We demonstrate this in both direct numerical simulations forced by volume stirring and in large eddy simulations of solar convectively driven small-scale dynamos. Low density plasmas such as stellar coronae tend to have large magnetic Prandtl numbers, i.e., the viscosity is large compared with the magnetic diffusivity. The regime in which viscous dissipation dominates over resistive dissipation for large magnetic Prandtl numbers was also previously found in large eddy simulations of the solar corona, i.e., our findings are a more fundamental property of MHD that is not just restricted to dynamos. Viscous energy dissipation is a consequence of positive Lorentz force work, which may partly correspond to particle acceleration in close-to-collisionless plasmas. This is, however, not modeled in the MHD approximation employed. By contrast, resistive energy dissipation on current sheets is expected to be unimportant in stellar coronae. △ Less

Submitted 23 May, 2019; v1 submitted 28 March, 2019; originally announced March 2019.

Comments: 9 pages, 6 figures, 2 tables, Astrophys. J., in press

Report number: NORDITA-2019-029

Journal ref: Astrophys. J. 879, 57 (2019)

What the sudden death of solar cycles can tell us about the nature of the solar interior

Authors: Scott W. McIntosh, Robert J. Leamon, Ricky Egeland, Mausumi Dikpati, Yuhong Fan, Matthias Rempel

Abstract: We observe the abrupt end of solar activity cycles at the Sun's equator by combining almost 140 years of observations from ground and space. These "terminator" events appear to be very closely related to the onset of magnetic activity belonging to the next sunspot cycle at mid-latitudes and the polar-reversal process at high-latitudes. Using multi-scale tracers of solar activity we examine the tim… ▽ More We observe the abrupt end of solar activity cycles at the Sun's equator by combining almost 140 years of observations from ground and space. These "terminator" events appear to be very closely related to the onset of magnetic activity belonging to the next sunspot cycle at mid-latitudes and the polar-reversal process at high-latitudes. Using multi-scale tracers of solar activity we examine the timing of these events in relation to the excitation of new activity and find that the time taken for the solar plasma to communicate this transition is of the order of one solar rotation, but could be shorter. Utilizing uniquely comprehensive solar observations from the Solar Terrestrial Relations Observatory (STEREO), and Solar Dynamics Observatory (SDO) we see that this transitional event is strongly longitudinal in nature. Combined, these characteristics imply that magnetic information is communicated through the solar interior rapidly. A range of possibilities exist to explain such behavior: the presence of magnetic reconnection in the deep interior, internal gravity waves on the solar tachocline, or that the magnetic fields present in the Sun's convection zone could be very large, with a poloidal field strengths reaching 50k - considerably larger than conventional explorations of solar and stellar dynamos estimate. Regardless of mechanism responsible, the rapid timescales demonstrated by the Sun's global magnetic field reconfiguration present strong constraints on first-principles numerical simulations of the solar interior and, by extension, other stars. △ Less

Submitted 3 June, 2019; v1 submitted 25 January, 2019; originally announced January 2019.

Comments: Submitted [10/2017] to, and rejected [01/2019] by, Nature Astronomy. Accepted Solar Physics [6/1/2019]

Small-scale dynamo simulations: Magnetic field amplification in exploding granules and the role of deep and shallow recirculation

Authors: Matthias Rempel

Abstract: We analyze recent high resolution photospheric small-scale dynamo simulations that were computed with the MURaM radiative MHD code. We focus the analysis on newly forming downflow lanes in exploding granules since they show how weakly magnetized regions in the photosphere (center of granules) evolve into strongly magnetized regions (downflow lanes). We find that newly formed downflow lanes exhibit… ▽ More We analyze recent high resolution photospheric small-scale dynamo simulations that were computed with the MURaM radiative MHD code. We focus the analysis on newly forming downflow lanes in exploding granules since they show how weakly magnetized regions in the photosphere (center of granules) evolve into strongly magnetized regions (downflow lanes). We find that newly formed downflow lanes exhibit initially mostly a laminar converging flow that amplifies the vertical magnetic field embedded in the granule from a few 10 G to field strengths exceeding 800 G. This results in extended magnetic sheets that have a length comparable to granular scales. Field amplification by turbulent shear happens first a few 100 km beneath the visible layers of the photosphere. Shallow recirculation transports the resulting turbulent field into the photosphere within minutes, after which the newly formed downflow lane shows a mix of strong magnetic sheets and turbulent field components. We stress in particular the role of shallow and deep recirculation for the organization and strength of magnetic field in the photosphere and discuss the photospheric and sub-photospheric energy conversion associated with the small-scale dynamo process. While the energy conversion through the Lorentz force depends only weakly on the saturation field strength (and therefore deep or shallow recirculation), it is strongly dependent on the magnetic Prandtl number. We discuss the potential of these findings for further constraining small-scale dynamo models through high resolution observations. △ Less

Submitted 22 May, 2018; originally announced May 2018.

Comments: 17 pages, 11 figures, ApJ in press

Evershed and counter-Evershed flows in sunspot MHD simulations

Authors: A. L. Siu-Tapia, M. Rempel, A. Lagg, S. K. Solanki

Abstract: There have been a few reports in the literature of counter-Evershed flows observed in well developed sunspot penumbrae, i.e. flows directed towards the umbra along penumbral filaments. Here we investigate the driving forces of such counter-Evershed flows in a radiative magnetohydrodynamic simulation of a sunspot and compare them with the forces acting on the normal Evershed flow. The simulation co… ▽ More There have been a few reports in the literature of counter-Evershed flows observed in well developed sunspot penumbrae, i.e. flows directed towards the umbra along penumbral filaments. Here we investigate the driving forces of such counter-Evershed flows in a radiative magnetohydrodynamic simulation of a sunspot and compare them with the forces acting on the normal Evershed flow. The simulation covers a timespan of 100 solar hours and generates an Evershed outflow exceeding 8 km s$^{-1}$ in the penumbra along radially aligned filaments where the magnetic field is almost horizontal. Additionally, the simulation produces a fast counter-Evershed flow (i.e., an inflow near $τ= 1$) in some regions within the penumbra, reaching peak flow speeds of $\sim$12 km s$^{-1}$. The counter-Evershed flows are transient and typically last a few hours before they turn into outflows again. By using the kinetic energy equation and evaluating its various terms in the simulation box, we found that the Evershed flow occurs due to overturning convection in a strongly inclined magnetic field while the counter-Evershed flows can be well described as siphon flows. △ Less

Submitted 4 December, 2017; originally announced December 2017.

Comments: Paper accepted for publication in The Astrophysical Journal

Transport of internetwork magnetic flux elements in the solar photosphere

Authors: Piyush Agrawal, Mark P. Rast, Milan Gosic, Luis R. Bellot Rubio, Matthias Rempel

Abstract: The motions of small-scale magnetic flux elements in the solar photosphere can provide some measure of the Lagrangian properties of the convective flow. Measurements of these motions have been critical in estimating the turbulent diffusion coefficient in flux-transport dynamo models and in determining the Alfven wave excitation spectrum for coronal heating models. We examine the motions of interne… ▽ More The motions of small-scale magnetic flux elements in the solar photosphere can provide some measure of the Lagrangian properties of the convective flow. Measurements of these motions have been critical in estimating the turbulent diffusion coefficient in flux-transport dynamo models and in determining the Alfven wave excitation spectrum for coronal heating models. We examine the motions of internetwork flux elements in a 24 hour long Hinode/NFI magnetogram sequence with 90 second cadence, and study both the scaling of their mean squared displacement and the shape of their displacement probability distribution as a function of time. We find that the mean squared displacement scales super-diffusively with a slope of about 1.48. Super-diffusive scaling has been observed in other studies for temporal increments as small as 5 seconds, increments over which ballistic scaling would be expected. Using high-cadence MURaM simulations, we show that the observed super-diffusive scaling at short temporal increments is a consequence of random changes in the barycenter positions due to flux evolution. We also find that for long temporal increments, beyond granular lifetimes, the observed displacement distribution deviates from that expected for a diffusive process, evolving from Rayleigh to Gaussian. This change in the distribution can be modeled analytically by accounting for supergranular advection along with motions due to granulation. These results complicate the interpretation of magnetic element motions as strictly advective or diffusive on short and long timescales and suggest that measurements of magnetic element motions must be used with caution in turbulent diffusion or wave excitation models. We propose that passive trace motions in measured photospheric flows may yield more robust transport statistics. △ Less

Submitted 12 December, 2017; v1 submitted 3 November, 2017; originally announced November 2017.

Comments: submitted to ApJ

Nature of grand minima and maxima from fully non-linear Flux-Transport Dynamos

Authors: F. Inceoglu, R. Arlt, M. Rempel

Abstract: We aim to investigate the nature and occurrence characteristics of grand solar minimum and maximum periods, which are observed in the solar proxy records such as 10Be and 14C, using a fully non-linear Babcock-Leighton type flux-transport dynamo including momentum and entropy equations. The differential rotation and meridional circulation are generated from the effect of turbulent Reynolds stress a… ▽ More We aim to investigate the nature and occurrence characteristics of grand solar minimum and maximum periods, which are observed in the solar proxy records such as 10Be and 14C, using a fully non-linear Babcock-Leighton type flux-transport dynamo including momentum and entropy equations. The differential rotation and meridional circulation are generated from the effect of turbulent Reynolds stress and are subjected to back-reaction from the magnetic field. To generate grand minimum and maximum-like periods in our simulations, we used random fluctuations in the angular momentum transport process, namely the Lambda-mechanism, and in the Babcock-Leighton mechanism. To characterise the nature and occurrences of the identified grand minima and maxima in our simulations, we used the waiting time distribution analyses, which reflects whether the underlying distribution arises from a random or a memory-bearing process. The results show that, in majority of the cases, the distributions of grand minima and maxima reveal that the nature of these events originates from memoryless processes. We also found that in our simulations the meridional circulation speed tends to be smaller during grand maximum, while it is faster during grand minimum periods. The radial differential rotation tend to be larger during grand maxima, while it is smaller during grand minima. The latitudinal differential rotation on the other hand is found to be larger during grand minima. △ Less

Submitted 24 October, 2017; originally announced October 2017.

Comments: accepted for publication in Astrophysical Journal

Journal ref: Astrophys. J. 848, 93 (2017)

Emergence of magnetic flux generated in a solar convective dynamo. I: Formation of Sunspots and Active regions, and Origin of Their Asymmetries

Authors: Feng Chen, Matthias Rempel, Yuhong Fan

Abstract: We present a realistic numerical model of sunspot and active region formation based on the emergence of flux bundles generated in a solar convective dynamo. To this end we use the magnetic and velocity fields in a horizontal layer near the top boundary of the solar convective dynamo simulation to drive realistic radiative-magnetohydrodynamic simulations of the upper most layers of the convection z… ▽ More We present a realistic numerical model of sunspot and active region formation based on the emergence of flux bundles generated in a solar convective dynamo. To this end we use the magnetic and velocity fields in a horizontal layer near the top boundary of the solar convective dynamo simulation to drive realistic radiative-magnetohydrodynamic simulations of the upper most layers of the convection zone. The main results are: (1) The emerging flux bundles rise with the mean speed of convective upflows, and fragment into small-scale magnetic elements that further rise to the photosphere, where bipolar sunspot pairs are formed through the coalescence of the small-scale magnetic elements. (2) Filamentary penumbral structures form when the sunspot is still growing through ongoing flux emergence. In contrast to the classical Evershed effect, the inflow seems to prevail over the outflow in a large part of the penumbra. (3) A well formed sunspot is a mostly monolithic magnetic structures that is anchored in a persistent deep-seated downdraft lane. The flow field outside the spot shows a giant vortex ring that comprises of an inflow below 15 Mm depth and an outflow above 15 Mm depth. (4) The sunspots successfully reproduce the fundamental properties of the observed solar active regions, including the more coherent leading spots with a stronger field strength, and the correct tilts of bipolar sunspot pairs. These asymmetries can be linked to the intrinsic asymmetries in the magnetic and flow fields adapted from the convective dynamo simulation. △ Less

Submitted 16 August, 2017; v1 submitted 19 April, 2017; originally announced April 2017.

Comments: 30 pages, 16 figures, accepted for publication in ApJ

Lower solar atmosphere and magnetism at ultra-high spatial resolution

Authors: Remo Collet, Serena Criscuoli, Ilaria Ermolli, Damian Fabbian, Nuno Guerreiro, Margit Haberreiter, Courtney Peck, Tiago M. D. Pereira, Matthias Rempel, Sami K. Solanki, Sven Wedemeyer-Boehm

Abstract: We present the scientific case for a future space-based telescope aimed at very high spatial and temporal resolution imaging of the solar photosphere and chromosphere. Previous missions (e.g., HINODE, SUNRISE) have demonstrated the power of observing the solar photosphere and chromosphere at high spatial resolution without contamination from Earth's atmosphere. We argue here that increased spatial… ▽ More We present the scientific case for a future space-based telescope aimed at very high spatial and temporal resolution imaging of the solar photosphere and chromosphere. Previous missions (e.g., HINODE, SUNRISE) have demonstrated the power of observing the solar photosphere and chromosphere at high spatial resolution without contamination from Earth's atmosphere. We argue here that increased spatial resolution (from currently 70 km to 25 km in the future) and high temporal cadence of the observations will vastly improve our understanding of the physical processes controlling solar magnetism and its characteristic scales. This is particularly important as the Sun's magnetic field drives solar activity and can significantly influence the Sun-Earth system. At the same time a better knowledge of solar magnetism can greatly improve our understanding of other astrophysical objects. △ Less

Submitted 7 December, 2016; originally announced December 2016.