Researcher profile

M. Rempel

M. Rempel contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2022arXiv

Chromospheric Extension of the MURaM Code

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.

preprint2022arXiv

Heating of the solar chromosphere through current dissipation

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.

preprint2022arXiv

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

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.

preprint2021arXiv

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

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.

preprint2021arXiv

The Coronal Veil

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.

preprint2021arXiv

The Dimmest State of the Sun

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.