BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

M. Rempel

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

ResearcherUnclaimed source record
astro-ph.SRphysics.space-phastro-ph.IMphysics.plasm-phastro-ph.EPastro-phastro-ph.HEComputationphysics.data-anphysics.flu-dyn

Catalog footprint

What is connected

30works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

30 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Observed and Simulated Power Spectra of Kinetic and Magnetic Energy retrieved with 2D inversions

We try to retrieve the power spectra with certainty to the highest spatial frequencies allowed by current instrumentation. For this, we use 2D inversion code that were able to recover information up to the instrumental diffraction limit. The retrieved power spectra have shallow slopes extending further down to much smaller scales than found before. They seem not to show any power law. The observed slopes at subgranular scales agree with those obtained from recent local dynamo simulations. Small differences are found for vertical component of kinetic energy that suggest that observations suffer from an instrumental effect that is not taken into account.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Efficient small-scale dynamo in solar convection zone

We investigate small-scale dynamo action in the solar convection zone through a series of high resolution MHD simulations in a local Cartesian domain with 1$R_\odot$ (solar radius) of horizontal extent and a radial extent from 0.715 to 0.96$R_\odot$. The dependence of the solution on resolution and diffusivity is studied. For a grid spacing of less than 350 km, the root mean square magnetic field strength near the base of the convection zone reaches 95% of the equipartition field strength (i.e. magnetic and kinetic energy are comparable). For these solutions the Lorentz force feedback on the convection velocity is found to be significant. The velocity near the base of the convection zone is reduced to 50% of the hydrodynamic one. In spite of a significant decrease of the convection velocity, the reduction in the enthalpy flux is relatively small, since the magnetic field also suppresses the horizontal mixing of the entropy between up- and downflow regions. This effect increases the amplitude of the entropy perturbation and makes convective energy transport more efficient. We discuss potential implications of these results for solar global convection and dynamo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical simulations of sunspot decay: On the penumbra -- Evershed flow -- moat flow connection

We present a series of high-resolution sunspot simulations that cover a time span of up to 100 hours. The simulation domain extends about 18 Mm in depth beneath the photosphere and 98 Mm horizontally. We use open boundary conditions that do not maintain the initial field structure against decay driven by convective motions. We consider two setups: A sunspot simulation with penumbra, and a "naked-spot" simulation in which we removed the penumbra after 20 hours through a change in the magnetic top boundary condition. While the sunspot has an Evershed outflow of 3-4 km/s, the naked spot is surrounded by an inflow of 1-2 km/s in close proximity. However, both spots are surrounded by an outflow on larger scales with a few 100 m/s flow speed in the photosphere. While the sunspot has almost constant magnetic flux content for the simulated time span of 3-4 days, the naked spot decays steadily at a rate of $10^{21}$ Mx/day. A region with reduced downflow filling factor, which is more extended for the sunspot, surrounds both spots. The absence of downflows perturbs the upflow/downflow massflux balance and leads to a large-scale radially overturning flow system, the photospheric component of this flow is to the observable moat flow. The reduction of the downflow filling factor also inhibits submergence of magnetic field in the proximity of the spots, which stabilizes them against decay. While this effect is present for both spots, it is more pronounced for the sunspot and explains the almost stationary magnetic flux content.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Comparison of inversion codes for polarized line formation in MHD simulations. I. Milne-Eddington codes

Milne-Eddington (M-E) inversion codes for the radiative transfer equation are the most widely used tools to infer the magnetic field from observations of the polarization signals in photospheric and chromospheric spectral lines. Unfortunately, a comprehensive comparison between the different M-E codes available to the solar physics community is still missing, and so is a physical interpretation of their inferences. In this contribution we offer a comparison between three of those codes (VFISV, ASP/HAO, and HeLIx$^+$). These codes are used to invert synthetic Stokes profiles that were previously obtained from realistic non-grey three-dimensional magnetohydrodynamical (3D MHD) simulations. The results of the inversion are compared with each other and with those from the MHD simulations. In the first case, the M-E codes retrieve values for the magnetic field strength, inclination and line-of-sight velocity that agree with each other within $σ_B \leq 35$ (Gauss), $σ_γ\leq 1.2°$, and $σ_{\rm v} \leq 10$ ms$^{-1}$, respectively. Additionally, M-E inversion codes agree with the numerical simulations, when compared at a fixed optical depth, within $σ_B \leq 130$ (Gauss), $σ_γ\leq 5°$, and $σ_{\rm v} \leq 320$ ms$^{-1}$. Finally, we show that employing generalized response functions to determine the height at which M-E codes measure physical parameters is more meaningful than comparing at a fixed geometrical height or optical depth. In this case the differences between M-E inferences and the 3D MHD simulations decrease to $σ_B \leq 90$ (Gauss), $σ_γ\leq 3°$, and $σ_{\rm v} \leq 90$ ms$^{-1}$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-resolution calculation of the solar global convection with the reduced speed of sound technique: II. Near surface shear layer with the rotation

We present a high-resolution, highly stratified numerical simulation of rotating thermal convection in a spherical shell. Our aim is to study in detail the processes that can maintain a near surface shear layer (NSSL) as inferred from helioseismology. Using the reduced speed of sound technique we can extend our global convection simulation to $0.99\,R_{\odot}$ and include near the top of our domain small-scale convection with short time scales that is only weakly influenced by rotation. We find the formation of a NSSL preferentially in high latitudes in the depth range $r=0.95-0.975R_\odot$. The maintenance mechanisms are summarized as follows. Convection under weak influence of rotation leads to Reynolds stresses that transport angular momentum radially inward in all latitudes. This leads to the formation of a strong poleward directed meridional flow and a NSSL, which is balanced in the meridional plane by forces resulting from the $\langle v'_r v'_θ\rangle$ correlation of turbulent velocities. The origin of the required correlations depends to some degree on latitude. In high latitudes a positive correlation $\langle v'_rv'_θ\rangle$ is induced in the NSSL by the poleward meridional flow whose amplitude increases with the radius, while a negative correlation is generated by the Coriolis force in bulk of the convection zone. In low latitudes a positive correlation $\langle v'_rv'_θ\rangle$ results from rotationally aligned convection cells ("banana cells"). The force caused by these Reynolds stresses is in balance with the Coriolis force in the NSSL.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-resolution calculations of the solar global convection with the reduced speed of sound technique: I. The structure of the convection and the magnetic field without the rotation

We carry out non-rotating high-resolution calculations of the solar global convection, which resolve convective scales of less than 10 Mm. To cope with the low Mach number conditions in the lower convection zone, we use the reduced speed of sound technique (RSST), which is simple to implement and requires only local communication in the parallel computation. In addition, the RSST allows us to expand the computational domain upward to about $0.99 R_{\odot}$ as it can also handle compressible flows. Using this approach, we study the solar convection zone on the global scale, including small-scale near-surface convection. In particular, we investigate the influence of the top boundary condition on the convective structure throughout the convection zone as well as on small-scale dynamo action. Our main conclusions are: 1. The small-scale downflows generated in the near-surface layer penetrate into deeper layers to some extent and excite small-scale turbulence in the region $>0.9R_\odot$, where $R_\odot$ is the solar radius. 2. In the deeper convection zone ($<0.9R_\odot$), the convection is not influenced by the location of the upper boundary. 3. Using an LES approach, we can achieve small-scale dynamo action and maintain a field of about $0.15-0.25 B_\mathrm{eq}$ throughout the convection zone, where $B_\mathrm{eq}$ is the equipartition magnetic field to the kinetic energy. 4. The overall dynamo efficiency varies significantly in the convection zone as a consequence of the downward directed Poynting flux and the depth variation of the intrinsic convective scales.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical simulations of active region scale flux emergence: From spot formation to decay

We present numerical simulations of active region scale flux emergence covering a time span of up to 6 days. Flux emergence is driven by a bottom boundary condition that advects a semi-torus of magnetic field with $1.7 \times 10^{22}$ Mx flux into the computational domain. The simulations show that, even in the absense of twist, the magnetic flux is able the rise through the upper $15.5$ Mm of the convection zone and emerge into the photosphere to form spots. We find that spot formation is sensitive to the persistence of upflows at the bottom boundary footpoints, i.e. a continuing upflow would prevent spot formation. In addition, the presence of a torus-aligned flow (such flow into the retrograde direction is expected from angular momentum conservation during the rise of flux ropes through the convection zone) leads to a significant asymmetry between the pair of spots, with the spot corresponding to the leading spot on the Sun being more axisymmetric and coherent, but also forming with a delay relative to the following spot. The spot formation phase transitions directly into a decay phase. Subsurface flows fragment the magnetic field and lead to intrusions of almost field free plasma underneath the photosphere. When such intrusions reach photospheric layers, the spot fragments. The time scale for spot decay is comparable to the longest convective time scales present in the simulation domain. We find that the dispersal of flux from a simulated spot in the first two days of the decay phase is consistent with self-similar decay by turbulent diffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical simulations of quiet Sun magnetism: On the contribution from a small-scale dynamo

We present a series of radiative MHD simulations addressing the origin and distribution of mixed polarity magnetic field in the solar photosphere. To this end we consider numerical simulations that cover the uppermost 2-6 Mm of the solar convection zone and we explore scales ranging from 2 km to 25 Mm. We study how the strength and distribution of magnetic field in the photosphere and subsurface layers depend on resolution, domain size and boundary conditions. We find that 50% of the magnetic energy at the τ=1 level comes from field with the less than 500 G strength and that 50% of the energy resides on scales smaller than about 100 km. While probability distribution functions are essentially independent of resolution, properly describing the spectral energy distribution requires grid spacings of 8 km or smaller. The formation of flux concentrations in the photosphere exceeding 1 kG requires a mean vertical field strength greater than 30-40 G at τ=1. The filling factor of kG flux concentrations increases with overall domain size as magnetic field becomes organized by larger, longer lived flow structures. A solution with a mean vertical field strength of around 85 G at τ=1 requires a subsurface RMS field strength increasing with depth at the same rate as the equipartition field strength. We consider this an upper limit for the quiet Sun field strength, which implies that most of the convection zone is magnetized close to equipartition. We discuss these findings in view of recent high-resolution spectropolarimetric observations of quiet Sun magnetism.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Role of Subsurface Flows in Solar Surface Convection: Modeling the Spectrum of Supergranular and Larger Scale Flows

We model the solar horizontal velocity power spectrum at scales larger than granulation using a two-component approximation to the mass continuity equation. The model takes four times the density scale height as the integral (driving) scale of the vertical motions at each depth. Scales larger than this decay with height from the deeper layers. Those smaller are assumed to follow a Kolomogorov turbulent cascade, with the total power in the vertical convective motions matching that required to transport the solar luminosity in a mixing length formulation. These model components are validated using large scale radiative hydrodynamic simulations. We reach two primary conclusions: 1. The model predicts significantly more power at low wavenumbers than is observed in the solar photospheric horizontal velocity spectrum. 2. Ionization plays a minor role in shaping the observed solar velocity spectrum by reducing convective amplitudes in the regions of partial helium ionization. The excess low wavenumber power is also seen in the fully nonlinear three-dimensional radiative hydrodynamic simulations employing a realistic equation of state. This adds to other recent evidence suggesting that the amplitudes of large scale convective motions in the Sun are significantly lower than expected. Employing the same feature tracking algorithm used with observational data on the simulation output, we show that the observed low wavenumber power can be reproduced in hydrodynamic models if the amplitudes of large scale modes in the deep layers are artificially reduced. Since the large scale modes have reduced amplitudes, modes on the scale of supergranulation and smaller remain important to convective heat flux even in the deep layers, suggesting that small scale convective correlations are maintained through the bulk of the solar convection zone.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Time-Distance Helioseismology of Two Realistic Sunspot Simulations

Linear time-distance helioseismic inversions are carried out using several filtering schemes to determine vector flow velocities within two $\sim100^2\,{\rm Mm^2}\times 20\,{\rm Mm}$ realistic magnetohydrodynamic sunspot simulations of 25~hr. One simulation domain contains a model of a full sunspot (i.e. one with both an umbra and penumbra), while the other contains a pore (i.e. a spot without a penumbra). The goal is to test current helioseismic methods using these state-of-the-art simulations of magnetic structures. We find that horizontal flow correlations between inversion and simulation flow maps are reasonably high ($\sim0.5$--0.8) in the upper 3~Mm at distances exceeding 25--30~Mm from spot center, but are substantially lower at smaller distances and larger depths. Inversions of forward-modeled travel times consistently outperform those of our measured travel times in terms of horizontal flow correlations, suggesting that our inability to recover flow structure near these active regions is largely due to the fact that we are unable to accurately measure travel times near strong magnetic features. In many cases the velocity amplitudes from the inversions underestimate those of the simulations by up to 50\%, possibly indicating nonlinearity of the forward problem. In every case, we find that our inversions are unable to recover the vertical flow structure of the simulations at any depth.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Validating Time-Distance Helioseismology With Realistic Quiet Sun Simulations

Linear time-distance helioseismic inversions are carried out for vector flow velocities using travel times measured from two $\sim 100^2\,{\rm Mm^2}\times 20\,{\rm Mm}$ realistic magnetohydrodynamic quiet-Sun simulations of about 20 hr. The goal is to test current seismic methods on these state-of-the-art simulations. Using recent three-dimensional inversion schemes, we find that inverted horizontal flow maps correlate well with the simulations in the upper $\sim 3$ Mm of the domains for several filtering schemes, including phase-speed, ridge, and combined phase-speed and ridge measurements. In several cases, however, the velocity amplitudes from the inversions severely underestimate those of the simulations, possibly indicating nonlinearity of the forward problem. We also find that, while near-surface inversions of the vertical velocites are best using phase-speed filters, in almost all other example cases these flows are irretrievable due to noise, suggesting a need for statistical averaging to obtain better inferences.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Magnetic field intensification by three-dimensional explosion process

We investigate an intensification mechanism for the magnetic field near the base of the solar convection zone that does not rely on differential rotation. Such mechanism in addition to differential rotation has been suggested by studies of flux emergence, which typically require field strength in excess of those provided by differential rotation alone. We study here a process in which potential energy of the superadiabatically stratified convection zone is converted into magnetic energy. This mechanism, know as explosion of magnetic flux tubes, has been previously studied in the thin flux tube approximation as well as two-dimensional MHD simulations, we expand the investigation to three-dimensional MHD simulations. Our main result is that enough intensification can be achieved in a three-dimensional magnetic flux sheet as long as the spatial scale of the imposed perturbation normal to the magnetic field is sufficiently large. When this spatial scale is small, the flux sheet tends to rise toward the surface, resulting in a significant decrease of the magnetic field amplification.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical Calculation of Convection with Reduced Speed of Sound Technique

Context. The anelastic approximation is often adopted in numerical calculation with low Mach number, such as stellar internal convection. This approximation requires frequent global communication, because of an elliptic partial differential equation. Frequent global communication is negative factor for the parallel computing with a large number of CPUs. Aims. The main purpose of this paper is to test the validity of a method that artificially reduces the speed of sound for the compressible fluid equations in the context of stellar internal convection. The reduction of speed of sound allows for larger time steps in spite of low Mach number, while the numerical scheme remains fully explicit and the mathematical system is hyperbolic and thus does not require frequent global communication. Methods. Two and three dimensional compressible hydrodynamic equations are solved numerically. Some statistical quantities of solutions computed with different effective Mach numbers (due to reduction of speed of sound) are compared to test the validity of our approach. Results. Numerical simulations with artificially reduced speed of sound are a valid approach as long as the effective Mach number (based on the reduced speed of sound) remains less than 0.7.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Waves as the source of apparent twisting motions in sunspot penumbrae

The motion of dark striations across bright filaments in a sunspot penumbra has become an important new diagnostic of convective gas flows in penumbral filaments. The nature of these striations has, however, remained unclear. Here we present an analysis of small scale motions in penumbral filaments in both simulations and observations. The simulations, when viewed from above, show fine structure with dark lanes running outwards from the dark core of the penumbral filaments. The dark lanes either occur preferentially on one side or alternate between both sides of the filament. We identify this fine structure with transverse (kink) oscillations of the filament, corresponding to a sideways swaying of the filament. These oscillations have periods in the range of 5-7 min and propagate outward and downward along the filament. Similar features are found in observed G-band intensity time series of penumbral filaments in a sunspot located near disk center obtained by the Broadband Filter Imager (BFI) on board {\it Hinode}. We also find that some filaments show dark striations moving to both sides of the filaments. Based on the agreement between simulations and observations we conclude that the motions of these striations are caused by transverse oscillations of the underlying bright filaments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Properties of Umbral Dots as Measured from the New Solar Telescope Data and MHD Simulations

We studied bright umbral dots (UDs) detected in a moderate size sunspot and compared their statistical properties to recent MHD models. The study is based on high resolution data recorded by the New Solar Telescope at the Big Bear Solar Observatory and 3D MHD simulations of sunspots. Observed UDs, living longer than 150 s, were detected and tracked in a 46 min long data set, using an automatic detection code. Total 1553 (620) UDs were detected in the photospheric (low chromospheric) data. Our main findings are: i) none of the analyzed UDs is precisely circular, ii) the diameter-intensity relationship only holds in bright umbral areas, and iii) UD velocities are inversely related to their lifetime. While nearly all photospheric UDs can be identified in the low chromospheric images, some small closely spaced UDs appear in the low chromosphere as a single cluster. Slow moving and long living UDs seem to exist in both the low chromosphere and photosphere, while fast moving and short living UDs are mainly detected in the photospheric images. Comparison to the 3D MHD simulations showed that both types of UDs display, on average, very similar statistical characteristics. However, i) the average number of observed UDs per unit area is smaller than that of the model UDs, and ii) on average, the diameter of model UDs is slightly larger than that of observed ones.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Discovery of a 1.6-year Magnetic Activity Cycle in the Exoplanet Host Star iota Horologii

The Mount Wilson Ca HK survey revealed magnetic activity variations in a large sample of solar-type stars with timescales ranging from 2.5 to 25 years. This broad range of cycle periods is thought to reflect differences in the rotational properties and the depths of the surface convection zones for stars with various masses and ages. In 2007 we initiated a long-term monitoring campaign of Ca II H and K emission for a sample of 57 southern solar-type stars to measure their magnetic activity cycles and their rotational properties when possible. We report the discovery of a 1.6-year magnetic activity cycle in the exoplanet host star iota Horologii, and we obtain an estimate of the rotation period that is consistent with Hyades membership. This is the shortest activity cycle so far measured for a solar-type star, and may be related to the short-timescale magnetic variations recently identified in the Sun and HD49933 from helio- and asteroseismic measurements. Future asteroseismic observations can be compared to those obtained near the magnetic minimum in 2006 to search for cycle-induced shifts in the oscillation frequencies. If such short activity cycles are common in F stars, then NASA's Kepler mission should observe their effects in many of its long-term asteroseismic targets.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling the Subsurface Structure of Sunspots

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Simulation of the Formation of a Solar Active Region

We present a radiative magnetohydrodynamics simulation of the formation of an Active Region on the solar surface. The simulation models the rise of a buoyant magnetic flux bundle from a depth of 7.5 Mm in the convection zone up into the solar photosphere. The rise of the magnetic plasma in the convection zone is accompanied by predominantly horizontal expansion. Such an expansion leads to a scaling relation between the plasma density and the magnetic field strength such that $B\propto\varrho^{1/2}$. The emergence of magnetic flux into the photosphere appears as a complex magnetic pattern, which results from the interaction of the rising magnetic field with the turbulent convective flows. Small-scale magnetic elements at the surface first appear, followed by their gradual coalescence into larger magnetic concentrations, which eventually results in the formation of a pair of opposite polarity spots. Although the mean flow pattern in the vicinity of the developing spots is directed radially outward, correlations between the magnetic field and velocity field fluctuations allow the spots to accumulate flux. Such correlations result from the Lorentz-force driven, counter-streaming motion of opposite-polarity fragments. The formation of the simulated Active Region is accompanied by transient light bridges between umbrae and umbral dots. Together with recent sunspot modeling, this work highlights the common magnetoconvective origin of umbral dots, light bridges and penumbral filaments.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Solar convection zone dynamics

A comprehensive understanding of the solar magnetic cycle requires detailed modeling of the solar interior including the maintenance and variation of large scale flows (differential rotation and meridional flow), the solar dynamo and the flux emergence process connecting the magnetic field in the solar convection zone with magnetic field in the photosphere and above. Due to the vast range of time and length scales encountered, a single model of the entire convection zone is still out of reach. However, a variety of aspects can be modeled through a combined approach of 3D MHD models and simplified descriptions. We will briefly review our current theoretical understanding of these processes based on numerical models of the solar interior.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Radiative MHD simulation of sunspot structure

Results of a 3D MHD simulation of a sunspot with a photospheric size of about 20 Mm are presented. The simulation has been carried out with the MURaM code, which includes a realistic equation of state with partial ionization and radiative transfer along many ray directions. The largely relaxed state of the sunspot shows a division in a central dark umbral region with bright dots and a penumbra showing bright filaments of about 2 to 3 Mm length with central dark lanes. By a process similar to the formation of umbral dots, the penumbral filaments result from magneto-convection in the form of upflow plumes, which become elongated by the presence of an inclined magnetic field: the upflow is deflected in the outward direction while the magnetic field is weakened and becomes almost horizontal in the upper part of the plume near the level of optical depth unity. A dark lane forms owing to the piling up of matter near the cusp-shaped top of the rising plume that leads to an upward bulging of the surfaces of constant optical depth. The simulated penumbral structure corresponds well to the observationally inferred interlocking-comb structure of the magnetic field with Evershed outflows along dark-laned filaments with nearly horizontal magnetic field and overturning perpendicular (`twisting') motion, which are embedded in a background of stronger and less inclined field. Photospheric spectral lines are formed at the very top and somewhat above the upflow plumes, so that they do not fully sense the strong flow as well as the large field inclination and significant field strength reduction in the upper part of the plume structures.

0 citations0 reviewsNo code linkedOpen access