Source author record

Matthias Rempel

Matthias Rempel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

18 published item(s)

preprint2026arXiv

Measuring the Coronal Magnetic Field with 2D Coronal Seismology: A Forward-Modeling Validation

In recent years, a two-dimensional (2D) coronal seismology technique applied to spectral-imaging data from the Coronal Multi-channel Polarimeter (CoMP) and UCoMP has enabled routine measurement of the global coronal magnetic field. The technique combines coronal transverse wave phase speed from Doppler measurements with electron densities from the Fe \sc{xiii}\rm{} 10798/10747 Å intensity ratio to infer the magnetic field strength, while the wave propagation directions from Doppler measurements trace the magnetic field direction. To validate the accuracy and robustness of this method, we use forward modeling of a MURaM simulation that produces open and closed magnetic structures with excited waves. From the synthetic Doppler velocity, Fe \sc{xiii}\rm{} infrared line intensities, and linear polarization signals, we apply the 2D coronal seismology technique to estimate the magnetic field strength and direction. A comparison with the simulation ground truth shows close agreement, indicating that the technique can recover the line-of-sight emissivity-weighted magnetic field direction and strength with high accuracy. We also perform a parameter-space analysis to quantify sensitivities of the method to parameter choice. These findings provide practical guidance for CoMP/UCoMP-like analysis and demonstrate that 2D coronal seismology can deliver reliable, LOS emissivity-weighted measurements of the coronal magnetic field from coronal wave observations.

preprint2025arXiv

Multi-height probing of horizontal flows in the solar photosphere

We tested whether simultaneous spectropolarimetric imaging in two magnetically sensitive optical spectral lines, which probe two different layers of the solar atmosphere (the photosphere and the temperature minimum), can help constrain the depth variation of horizontal flows. We first tested the feasibility of our method using Fourier local correlation tracking (FLCT) to track physical quantities at different optical depths ($\logτ_{500}={-1,-2,-3,-4}$) in an atmosphere simulated with the MURaM code. We then inferred the horizontal distribution of the LOS magnetic field component from synthetic spectropolarimetric observations of Fe I 525.0 nm and Mg I b2 spectral lines, applied FLCT to the time sequence of these synthetic magnetograms, and compared our findings with the original height-dependent horizontal velocities. Tracking the LOS magnetic field component (which coincides with the vertical component at the disk center) yields horizontal velocities that, after appropriate temporal and spatial averaging, agree excellently with the horizontal component of the simulated velocities, both calculated at constant $τ_{500}$ surfaces, up to the temperature minimum ($\logτ_{500}=-3$). When tracking the temperature at constant $τ_{500}$ surfaces, this agreement already breaks down completely at the mid photosphere ($\logτ_{500}=-2$). Tracking the vertical component of the magnetic field inferred from synthetic observations of the Fe I 525.0 nm and the Mg I b2 spectral lines yields a satisfactory inference of the horizontal velocities in the mid-photosphere ($\logτ_{500}\approx-1$) and the temperature minimum ($\logτ_{500}\approx-3$), respectively. Our results indicate that high-spatial-resolution spectropolarimetric imaging in solar spectral lines can provide meaningful information about the horizontal plasma velocities over a range of heights.

preprint2021arXiv

Convolutional Neural Networks and Stokes Response Functions

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.

preprint2021arXiv

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

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2016arXiv

A low upper limit on the subsurface rise speed of solar active regions

Magnetic field emerges at the surface of the Sun as sunspots and active regions. This process generates a poloidal magnetic field from a rising toroidal flux tube, it is a crucial but poorly understood aspect of the solar dynamo. The emergence of magnetic field is also important because it is a key driver of solar activity. We show that measurements of horizontal flows at the solar surface around emerging active regions, in combination with numerical simulations of solar magnetoconvection, can constrain the subsurface rise speed of emerging magnetic flux. The observed flows imply that the rise speed of the magnetic field is no larger than 150 m/s at a depth of 20 Mm, that is, well below the prediction of the (standard) thin flux tube model but in the range expected for convective velocities at this depth. We conclude that convective flows control the dynamics of rising flux tubes in the upper layers of the Sun and cannot be neglected in models of flux emergence.

preprint2016arXiv

Extension of the MURaM radiative MHD code for coronal simulations

We present a new version of the MURaM radiative MHD code that allows for simulations spanning from the upper convection zone into the solar corona. We implemented the relevant coronal physics in terms of optically thin radiative loss, field aligned heat conduction and an equilibrium ionization equation of state. We artificially limit the coronal Alfv{é}n and heat conduction speeds to computationally manageable values using an approximation to semi-relativistic MHD with an artificially reduced speed of light (Boris correction). We present example solutions ranging from quiet to active Sun in order to verify the validity of our approach. We quantify the role of numerical diffusivity for the effective coronal heating. We find that the (numerical) magnetic Prandtl number determines the ratio of resistive to viscous heating and that owing to the very large magnetic Prandtl number of the solar corona, heating is expected to happen predominantly through viscous dissipation. We find that reasonable solutions can be obtained with values of the reduced speed of light just marginally larger than the maximum sound speed. Overall this leads to a fully explicit code that can compute the time evolution of the solar corona in response to photospheric driving using numerical time steps not much smaller than $0.1$ seconds. Numerical simulations of the coronal response to flux emergence covering a time span of a few days are well within reach using this approach.

preprint2014arXiv

Photon mean free paths, scattering, and ever-increasing telescope resolution

We revisit an old question: what are the effects of observing stratified atmospheres on scales below a photon mean free path? The mean free path of photons emerging from the solar photosphere and chromosphere is near 100 km. Using current 1m-class telescopes, the mean free path is on the order of the angular resolution. But the Daniel K. Inoue Solar Telescope will have a diffraction limit of 0.020'' near the atmospheric cutoff at 310nm, corresponding to 14 km at the solar surface. Even a small amount of scattering in the source function leads to physical smearing due to this solar 'fog', with effects similar to a degradation of the telescope PSF. We discuss a unified picture that depends simply on the nature and amount of scattering in the source function. Scalings are derived from which the scattering in the solar atmosphere can be transcribed into an effective Strehl ratio, a quantity useful to observers. Observations in both permitted (e.g., Fe I 630.2 nm) and forbidden (Fe I 525.0 nm) lines will shed light on both instrumental performance as well as on small scale structures in the solar atmosphere.

preprint2012arXiv

Numerical sunspot models: Robustness of photospheric velocity and magnetic field structure

MHD simulations of sunspots have successfully reproduced many aspects of sunspot fine structure as consequence of magneto convection in inclined magnetic field. We study how global sunspot properties and penumbral fine structure depend on the magnetic top boundary condition as well as on grid spacing. The overall radial extent of the penumbra is subject to the magnetic top boundary condition. All other aspects of sunspot structure and penumbral fine structure are resolved at an acceptable level starting from a grid resolution of 48 [24] km (horizontal [vertical]). We find that the amount of inverse polarity flux and the overall amount of overturning convective motions in the penumbra are robust with regard to both, resolution and boundary conditions. At photospheric levels Evershed flow channels are strongly magnetized. We discuss in detail the relation between velocity and magnetic field structure in the photosphere and point out observational consequences.

preprint2012arXiv

On the Amplitude of Convective Velocities in the Deep Solar Interior

We obtain lower limits on the amplitude of convective velocities in the deep solar convection zone based only on the observed properties of the differential rotation and meridional circulation together with simple and robust dynamical balances obtained from the fundamental MHD equations. The linchpin of the approach is the concept of gyroscopic pumping whereby the meridional circulation across isosurfaces of specific angular momentum is linked to the angular momentum transport by the convective Reynolds stress. We find that the amplitude of the convective velocity must be at least 30 m s$^{-1}$ in the upper CZ ($r \sim 0.95 R$) and at least 8 m s$^{-1}$ in the lower CZ ($r \sim 0.75 R$) in order to be consistent with the observed mean flows. Using the base of the near-surface shear layer as a probe of the rotational influence, we are further able to show that the characteristic length scale of deep convective motions must be no smaller than 5.5--30 Mm. These results are compatible with convection models but suggest that the efficiency of the turbulent transport assumed in advection-dominated flux-transport dynamo models is generally not consistent with the mean flows they employ.

preprint2011arXiv

A more realistic representation of overshoot at the base of the solar convective envelope as seen by helioseismology

The stratification near the base of the Sun's convective envelope is governed by processes of convective overshooting and element diffusion, and the region is widely believed to play a key role in the solar dynamo. The stratification in that region gives rise to a characteristic signal in the frequencies of solar p modes, which has been used to determine the depth of the solar convection zone and to investigate the extent of convective overshoot. Previous helioseismic investigations have shown that the Sun's spherically symmetric stratification in this region is smoother than that in a standard solar model without overshooting, and have ruled out simple models incorporating overshooting, which extend the region of adiabatic stratification and have a more-or-less abrupt transition to subadiabatic stratification at the edge of the overshoot region. In this paper we consider physically motivated models which have a smooth transition in stratification bridging the region from the lower convection zone to the radiative interior beneath. We find that such a model is in better agreement with the helioseismic data than a standard solar model.

preprint2011arXiv

Can overturning motions in penumbral filaments be detected?

Numerical simulations indicate that the filamentation of sunspot penumbrae and the associated systematic outflow (the Evershed effect) are due to convectively driven fluid motions constrained by the inclined magnetic field. We investigate whether these motions, in particular the upflows in the bright filaments and the downflows at their edges can be reliably observed with existing instrumentation. We use a snapshot from a sunspot simulation to calculate 2D maps of synthetic line profiles for the spectral lines Fe\sci 7090.4 Å~ and C\sci 5380.34 Å. The maps are spatially and spectrally degraded according to typical instrument properties. Line-of-sight velocities are determined from line bisector shifts. We find that the detectability of the convective flows is strongly affected by spatial smearing, particularly so for the downflows. Furthermore, the line-of-sight velocities are dominated by the Evershed flow unless the observation is made very near to disk center. These problems may have compromised recent attempts to detect overturning penumbral convection. Lines with a low formation height are best suited to detect the convective flows.

preprint2011arXiv

Penumbral fine structure and driving mechanisms of large-scale flows in simulated sunspots

We analyze in detail the penumbral structure found in a recent radiative MHD simulation. Near tau=1, the simulation produces penumbral fine structure consistent with the observationally inferred interlocking comb structure. Fast outflows exceeding 8 km/s are present along almost horizontal stretches of the magnetic field; in the outer half of the penumbra, we see opposite polarity flux indicating flux returning beneath the surface. The bulk of the penumbral brightness is maintained by small-scale motions turning over on scales shorter than the length of a typical penumbral filament. The resulting vertical rms velocity at tau=1 is about half of that found in the quiet Sun. Radial outflows in the sunspot penumbra have two components. In the uppermost few 100 km, fast outflows are driven primarily through the horizontal component of the Lorentz force, which is confined to narrow boundary layers beneath tau=1, while the contribution from horizontal pressure gradients is reduced in comparison to granulation as a consequence of anisotropy. The resulting Evershed flow reaches its peak velocity near tau=1 and falls off rapidly with height. Outflows present in deeper layers result primarily from a preferred ring-like alignment of convection cells surrounding the sunspot. These flows reach amplitudes of about 50% of the convective rms velocity rather independent of depth. A preference for the outflow results from a combination of Lorentz force and pressure driving. While the Evershed flow dominates by velocity amplitude, most of the mass flux is present in deeper layers and likely related to a large-scale moat flow.

preprint2011arXiv

Subsurface magnetic field and flow structure of simulated sunspots

We present a series of numerical sunspot models addressing the subsurface field and flow structure in up to 16 Mm deep domains covering up to 2 days of temporal evolution. Changes in the photospheric appearance of the sunspots are driven by subsurface flows in several Mm depth. Most of magnetic field is pushed into a downflow vertex of the subsurface convection pattern, while some fraction of the flux separates from the main trunk of the spot. Flux separation in deeper layers is accompanied in the photosphere with light bridge formation in the early stages and formation of pores separating from the spot at later stages. Over a time scale of less than a day we see the development of a large scale flow pattern surrounding the sunspots, which is dominated by a radial outflow reaching about 50% of the convective rms velocity in amplitude. Several components of the large scale flow are found to be independent from the presence of a penumbra and the associated Evershed flow. While the simulated sunspots lead to blockage of heat flux in the near surface layers, we do not see compelling evidence for a brightness enhancement in their periphery. We further demonstrate that the influence of the bottom boundary condition on the stability and long-term evolution of the sunspot is significantly reduced in a 16 Mm deep domain compared to the shallower domains considered previously.

preprint2010arXiv

3D numerical MHD modeling of sunspots with radiation transport

Sunspot fine structure has been modeled in the past by a combination of idealized magneto-convection simulations and simplified models that prescribe the magnetic field and flow structure to a large degree. Advancement in numerical methods and computing power has enabled recently 3D radiative MHD simulations of entire sunspots with sufficient resolution to address details of umbral dots and penumbral filaments. After a brief review of recent developments we focus on the magneto-convective processes responsible for the complicated magnetic structure of the penumbra and the mechanisms leading to the driving of strong horizontal outflows in the penumbra (Evershed effect). The bulk of energy and mass is transported on scales smaller than the radial extent of the penumbra. Strong horizontal outflows in the sunspot penumbra result from a redistribution of kinetic energy preferring flows along the filaments. This redistribution is facilitated primarily through the Lorentz force, while horizontal pressure gradients play only a minor role. The Evershed flow is strongly magnetized: While we see a strong reduction of the vertical field, the horizontal field component is enhanced within filaments.

preprint2006arXiv

Influence of random fluctuations in the Lambda-effect on meridional flow and differential rotation

We present a mean field model based on the approach taken by Rempel (astro-ph/0604451) in order to investigate the influence of stochastic fluctuations in the Reynolds stresses on meridional flow and differential rotation. The stochastic fluctuations found in the meridional flow pattern directly resemble the stochastic fluctuations of the Reynolds stresses, while the stochastic fluctuations in the differential rotation are smaller by almost two orders of magnitude. It is further found that the correlation length and time scale of the stochastic fluctuations have only a weak influence on meridional flow, but a significant influence on the magnitude of variations in the differential rotation. We analyze the energy fluxes within the model to estimate time scales for the replenishment of differential rotation and meridional flow. We find that the time scale for the replenishment of differential rotation (~10 years) is nearly four orders of magnitude longer than the time scale for the replenishment of meridional flow, which explains the differences in the response to stochastic fluctuations of the Reynolds stress found for both flow fields.

preprint2006arXiv

Transport of toroidal magnetic field by the meridional flow at the base of the solar convection zone

In this paper we discuss the transport of toroidal magnetic field by a weak meridional flow at the base of the convection zone. We utilize the differential rotation and meridional flow model developed by Rempel and incorporate feedback of a purely toroidal magnetic field in two ways: directly through the Lorentz force (magnetic tension) and indirectly through quenching of the turbulent viscosity, which affects the parametrized turbulent angular momentum transport in the model. In the case of direct Lorentz force feedback we find that a meridional flow with an amplitude of around 2 m/s can transport a magnetic field with a strength of 20 to 30 kG. Quenching of turbulent viscosity leads to deflection of the meridional flow from the magnetized region and a significant reduction of the transport velocity if the magnetic field is above equipartition strength.