Science Platforms for Heliophysics Data Analysis
We recommend that NASA maintain and fund science platforms that enable interactive and scalable data analysis in order to maximize the scientific return of data collected from space-based instruments.
Source author record
Mark C. M. Cheung
Mark C. M. Cheung 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
21works
10topics
4close collaborators
Research graph
See the researcher in context
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
21 published item(s)
preprint2022arXiv
Coronal Mass Ejections and Dimmings: A Comparative Study using MHD Simulations and SDO Observations
Solar coronal dimmings have been observed extensively in the past two decades. Due to their close association with coronal mass ejections (CMEs), there is a critical need to improve our understanding of the physical processes that cause dimmings as well as their relationship with CMEs. In this study, we investigate coronal dimmings by combining simulation and observational efforts. By utilizing a data-constrained global magnetohydrodynamics model (AWSoM: Alfven-wave Solar Model), we simulate coronal dimmings resulting from different CME energetics and flux rope configurations. We synthesize the emissions of different EUV spectral bands/lines and compare with SDO/AIA and EVE observations. A detailed analysis of the simulation and observation data suggests that the transient dimming / brightening are related to plasma heating processes, while the long-lasting core and remote dimmings are caused by mass loss process induced by the CME. Moreover, the interaction between the erupting flux rope with different orientations and the global solar corona could significantly influence the coronal dimming patterns. Using metrics such as dimming depth and dimming slope, we investigate the relationship between dimmings and CME properties (e.g., CME mass, CME speed) in the simulation. Our result suggests that coronal dimmings encode important information about the associated CMEs, which provides a physical basis for detecting stellar CMEs from distant solar-like stars.
0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv
Global geomagnetic perturbation forecasting using Deep Learning
Geomagnetically Induced Currents (GICs) arise from spatio-temporal changes to Earth's magnetic field which arise from the interaction of the solar wind with Earth's magnetosphere, and drive catastrophic destruction to our technologically dependent society. Hence, computational models to forecast GICs globally with large forecast horizon, high spatial resolution and temporal cadence are of increasing importance to perform prompt necessary mitigation. Since GIC data is proprietary, the time variability of horizontal component of the magnetic field perturbation (dB/dt) is used as a proxy for GICs. In this work, we develop a fast, global dB/dt forecasting model, which forecasts 30 minutes into the future using only solar wind measurements as input. The model summarizes 2 hours of solar wind measurement using a Gated Recurrent Unit, and generates forecasts of coefficients which are folded with a spherical harmonic basis to enable global forecasts. When deployed, our model produces results in under a second, and generates global forecasts for horizontal magnetic perturbation components at 1-minute cadence. We evaluate our model across models in literature for two specific storms of 5 August 2011 and 17 March 2015, while having a self-consistent benchmark model set. Our model outperforms, or has consistent performance with state-of-the-practice high time cadence local and low time cadence global models, while also outperforming/having comparable performance with the benchmark models. Such quick inferences at high temporal cadence and arbitrary spatial resolutions may ultimately enable accurate forewarning of dB/dt for any place on Earth, resulting in precautionary measures to be taken in an informed manner.
0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv
Global Earth Magnetic Field Modeling and Forecasting with Spherical Harmonics Decomposition
Modeling and forecasting the solar wind-driven global magnetic field perturbations is an open challenge. Current approaches depend on simulations of computationally demanding models like the Magnetohydrodynamics (MHD) model or sampling spatially and temporally through sparse ground-based stations (SuperMAG). In this paper, we develop a Deep Learning model that forecasts in Spherical Harmonics space 2, replacing reliance on MHD models and providing global coverage at one minute cadence, improving over the current state-of-the-art which relies on feature engineering. We evaluate the performance in SuperMAG dataset (improved by 14.53%) and MHD simulations (improved by 24.35%). Additionally, we evaluate the extrapolation performance of the spherical harmonics reconstruction based on sparse ground-based stations (SuperMAG), showing that spherical harmonics can reliably reconstruct the global magnetic field as evaluated on MHD simulation.
0 citations0 reviewsNo code linkedOpen access
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.
0 citations0 reviewsNo code linkedOpen access
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.
0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv
SynthIA: A Synthetic Inversion Approximation for the Stokes Vector Fusing SDO and Hinode into a Virtual Observatory
Both NASA's Solar Dynamics Observatory (SDO) and the JAXA/NASA Hinode mission include spectropolarimetric instruments designed to measure the photospheric magnetic field. SDO's Helioseismic and Magnetic Imager (HMI) emphasizes full-disk high-cadence and good spatial resolution data acquisition while Hinode's Solar Optical Telescope Spectro-Polarimeter (SOT-SP) focuses on high spatial resolution and spectral sampling at the cost of a limited field of view and slower temporal cadence. This work introduces a deep-learning system named SynthIA (Synthetic Inversion Approximation), that can enhance both missions by capturing the best of each instrument's characteristics. We use SynthIA to produce a new magnetogram data product, SynodeP (Synthetic Hinode Pipeline), that mimics magnetograms from the higher spectral resolution Hinode/SOT-SP pipeline, but is derived from full-disk, high-cadence, and lower spectral-resolution SDO/HMI Stokes observations. Results on held-out data show that SynodeP has good agreement with the Hinode/SOT-SP pipeline inversions, including magnetic fill fraction, which is not provided by the current SDO/HMI pipeline. SynodeP further shows a reduction in the magnitude of the 24-hour oscillations present in the SDO/HMI data. To demonstrate SynthIA's generality, we show the use of SDO/AIA data and subsets of the HMI data as inputs, which enables trade-offs between fidelity to the Hinode/SOT-SP inversions, number of observations used, and temporal artifacts. We discuss possible generalizations of SynthIA and its implications for space weather modeling. This work is part of the NASA Heliophysics DRIVE Science Center (SOLSTICE) at the University of Michigan under grant NASA 80NSSC20K0600E, and will be open-sourced.
0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv
The solar internetwork. III. Unipolar versus bipolar flux appearance
Small-scale internetwork (IN) magnetic fields are considered to be the main building blocks of the quiet Sun magnetism. For this reason, it is crucial to understand how they appear on the solar surface. Here, we employ a high-resolution, high-sensitivity, long-duration Hinode/NFI magnetogram sequence to analyze the appearance modes and spatio-temporal evolution of individual IN magnetic elements inside a supergranular cell at the disk center. From identification of flux patches and magnetofrictional simulations, we show that there are two distinct populations of IN flux concentrations: unipolar and bipolar features. Bipolar features tend to be bigger and stronger than unipolar features. They also live longer and carry more flux per feature. Both types of flux concentrations appear uniformly over the solar surface. However, we argue that bipolar features truly represent the emergence of new flux on the solar surface, while unipolar features seem to be formed by coalescence of background flux. Magnetic bipoles appear at a faster rate than unipolar features (68 as opposed to 55 Mx cm$^{-2}$ day$^{-1}$), and provide about 70% of the total instantaneous IN flux detected in the interior of the supergranule.
0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv
Comparative Study of Data-driven Solar Coronal Field Models Using a Flux Emergence Simulation as a Ground-truth Data Set
For a better understanding of magnetic field in the solar corona and dynamic activities such as flares and coronal mass ejections, it is crucial to measure the time-evolving coronal field and accurately estimate the magnetic energy. Recently, a new modeling technique called the data-driven coronal field model, in which the time evolution of magnetic field is driven by a sequence of photospheric magnetic and velocity field maps, has been developed and revealed the dynamics of flare-productive active regions. Here we report on the first qualitative and quantitative assessment of different data-driven models using a magnetic flux emergence simulation as a ground-truth (GT) data set. We compare the GT field with those reconstructed from the GT photospheric field by four data-driven algorithms. It is found that, at least, the flux rope structure is reproduced in all coronal field models. Quantitatively, however, the results show a certain degree of model dependence. In most cases, the magnetic energies and relative magnetic helicity are comparable to or at most twice of the GT values. The reproduced flux ropes have a sigmoidal shape (consistent with GT) of various sizes, a vertically-standing magnetic torus, or a packed structure. The observed discrepancies can be attributed to the highly non-force-free input photospheric field, from which the coronal field is reconstructed, and to the modeling constraints such as the treatment of background atmosphere, the bottom boundary setting, and the spatial resolution.
0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv
Non-Neutralized Electric Current of Active Regions Explained as a Projection Effect
Active regions (ARs) often possess an observed net electric current in a single magnetic polarity. We show that such "non-neutralized" currents can arise from a geometric projection effect when a twisted flux tube obliquely intersects the photosphere. To this end, we emulate surface maps of an emerging AR by sampling horizontal slices of a semi-torus flux tube at various heights. Although the tube has no net toroidal current, its poloidal current, when projected along the vertical direction, amounts to a significant non-neutralized component on the surface. If the tube emerges only partially as in realistic settings, the non-neutralized current will 1) develop as double ribbons near the sheared polarity inversion line, (2) positively correlate with the twist, and 3) reach its maximum before the magnetic flux. The projection effect may be important to the photospheric current distribution, in particular during the early stages of flux emergence.
0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv
The PDFI_SS Electric Field Inversion Software
We describe the PDFI_SS software library, which is designed to find the electric field at the Sun's photosphere from a sequence of vector magnetogram and Doppler velocity measurements, and estimates of horizontal velocities obtained from local correlation tracking using the recently upgraded FLCT code. The library, a collection of Fortran subroutines, uses the "PDFI" technique described by Kazachenko et al. (2014), but modified for use in spherical, Plate-Carrée geometry on a staggered grid. The domain over which solutions are found is a subset of the global spherical surface, defined by user-specified limits of colatitude and longitude. Our staggered-grid approach, based on that of Yee (1966), is more conservative and self-consistent compared to the centered, Cartesian grid used by Kazachenko et al. (2014). The library can be used to compute an end-to-end solution for electric fields from data taken by the HMI instrument aboard NASA's SDO Mission. This capability has been incorporated into the HMI pipeline processing system operating at SDO's JSOC. The library is written in a general and modular way so that the calculations can be customized to modify or delete electric field contributions, or used with other data sets. Other applications include "nudging" numerical models of the solar atmosphere to facilitate assimilative simulations. The library includes an ability to compute "global" (whole-Sun) electric field solutions. The library also includes an ability to compute Potential Magnetic Field solutions in spherical coordinates. This distribution includes a number of test programs which allow the user to test the software.
0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv
The Multi-slit Approach to Coronal Spectroscopy with the Multi-slit Solar Explorer (MUSE)
The Multi-slit Solar Explorer (MUSE) is a proposed mission aimed at understanding the physical mechanisms driving the heating of the solar corona and the eruptions that are at the foundation of space weather. MUSE contains two instruments, a multi-slit EUV spectrograph and a context imager. It will simultaneously obtain EUV spectra (along 37 slits) and context images with the highest resolution in space (0.33-0.4 arcsec) and time (1-4 s) ever achieved for the transition region and corona. The MUSE science investigation will exploit major advances in numerical modeling, and observe at the spatial and temporal scales on which competing models make testable and distinguishable predictions, thereby leading to a breakthrough in our understanding of coronal heating and the drivers of space weather. By obtaining spectra in 4 bright EUV lines (Fe IX 171A, Fe XV 284A, Fe XIX-XXI 108A) covering a wide range of transition region and coronal temperatures along 37 slits simultaneously, MUSE will be able to "freeze" the evolution of the dynamic coronal plasma. We describe MUSE's multi-slit approach and show that the optimization of the design minimizes the impact of spectral lines from neighboring slits, generally allowing line parameters to be accurately determined. We also describe a Spectral Disambiguation Code to resolve multi-slit ambiguity in locations where secondary lines are bright. We use simulations of the corona and eruptions to perform validation tests and show that the multi-slit disambiguation approach allows accurate determination of MUSE observables in locations where significant multi-slit contamination occurs.
0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv
Homologous Helical Jets: Observations by IRIS, SDO and Hinode and Magnetic Modeling with Data-Driven Simulations
We report on observations of recurrent jets by instruments onboard the Interface Region Imaging Spectrograph (IRIS), Solar Dynamics Observatory (SDO) and Hinode spacecrafts. Over a 4-hour period on July 21st 2013, recurrent coronal jets were observed to emanate from NOAA Active Region 11793. FUV spectra probing plasma at transition region temperatures show evidence of oppositely directed flows with components reaching Doppler velocities of +/- 100 km/s. Raster Doppler maps using a Si IV transition region line show all four jets to have helical motion of the same sense. Simultaneous observations of the region by SDO and Hinode show that the jets emanate from a source region comprising a pore embedded in the interior of a supergranule. The parasitic pore has opposite polarity flux compared to the surrounding network field. This leads to a spine-fan magnetic topology in the coronal field that is amenable to jet formation. Time-dependent data-driven simulations are used to investigate the underlying drivers for the jets. These numerical experiments show that the emergence of current-carrying magnetic field in the vicinity of the pore supplies the magnetic twist needed for recurrent helical jet formation.
0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv
Light Bridge in a Developing Active Region. I. Observation of Light Bridge and its Dynamic Activity Phenomena
Light bridges, the bright structures that divide the umbra of sunspots and pores into smaller pieces, are known to produce wide variety of activity events in solar active regions (ARs). It is also known that the light bridges appear in the assembling process of nascent sunspots. The ultimate goal of this series of papers is to reveal the nature of light bridges in developing ARs and the occurrence of activity events associated with the light bridge structures from both observational and numerical approaches. In this first paper, exploiting the observational data obtained by Hinode, IRIS, and Solar Dynamics Observatory (SDO), we investigate the detailed structure of the light bridge in NOAA AR 11974 and its dynamic activity phenomena. As a result, we find that the light bridge has a weak, horizontal magnetic field, which is transported from the interior by large-scale convective upflow and is surrounded by strong, vertical fields of adjacent pores. In the chromosphere above the bridge, a transient brightening occurs repeatedly and intermittently, followed by a recurrent dark surge ejection into higher altitudes. Our analysis indicates that the brightening is the plasma heating due to magnetic reconnection at lower altitudes, while the dark surge is the cool, dense plasma ejected from the reconnection region. From the observational results, we conclude that the dynamic activity observed in a light bridge structure such as chromospheric brightenings and dark surge ejections are driven by magnetoconvective evolution within the light bridge and its interaction with surrounding magnetic fields.
0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv
Light Bridge in a Developing Active Region. II. Numerical Simulation of Flux Emergence and Light Bridge Formation
Light bridges, the bright structure dividing umbrae in sunspot regions, show various activity events. In Paper I, we reported on analysis of multi-wavelength observations of a light bridge in a developing active region (AR) and concluded that the activity events are caused by magnetic reconnection driven by magnetconvective evolution. The aim of this second paper is to investigate the detailed magnetic and velocity structures and the formation mechanism of light bridges. For this purpose, we analyze numerical simulation data from a radiative magnetohydrodynamics model of an emerging AR. We find that a weakly-magnetized plasma upflow in the near-surface layers of the convection zone is entrained between the emerging magnetic bundles that appear as pores at the solar surface. This convective upflow continuously transports horizontal fields to the surface layer and creates a light bridge structure. Due to the magnetic shear between the horizontal fields of the bridge and the vertical fields of the ambient pores, an elongated cusp-shaped current layer is formed above the bridge, which may be favorable for magnetic reconnection. The striking correspondence between the observational results of Paper I and the numerical results of this paper provides a consistent physical picture of light bridges. The dynamic activity phenomena occur as a natural result of the bridge formation and its convective nature, which has much in common with those of umbral dots and penumbral filaments.
0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv
Multi-Parametric Study of Rising 3D Buoyant Flux Tubes in an Adiabatic Stratification Using AMR
We study the buoyant rise of magnetic flux tubes embedded in an adiabatic stratification using two-and three-dimensional, MHD simulations. We analyze the dependence of the tube evolution on the field line twist and on the curvature of the tube axis in different diffusion regimes. To be able to achieve a comparatively high spatial resolution we use the FLASH code, which has a built-in Adaptive Mesh Refinement (AMR) capability. Our 3D experiments reach Reynolds numbers that permit a reasonable comparison of the results with those of previous 2D simulations. When the experiments are run without AMR, hence with a comparatively large diffusivity, the amount of longitudinal magnetic flux retained inside the tube increases with the curvature of the tube axis. However, when a low-diffusion regime is reached by using the AMR algorithms, the magnetic twist is able to prevent the splitting of the magnetic loop into vortex tubes and the loop curvature does not play any significant role. We detect the generation of vorticity in the main body of the tube of opposite sign on the opposite sides of the apex. This is a consequence of the inhomogeneity of the azimuthal component of the field on the flux surfaces. The lift force associated with this global vorticity makes the flanks of the tube move away from their initial vertical plane in an antisymmetric fashion. The trajectories have an oscillatory motion superimposed, due to the shedding of vortex rolls to the wake, which creates a Von Karman street.
0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv
The Coronal Global Evolutionary Model: Using HMI Vector Magnetogram and Doppler Data to Model the Buildup of Free Magnetic Energy in the Solar Corona
The most violent space weather events (eruptive solar flares and coronal mass ejections) are driven by the release of free magnetic energy stored in the solar corona. Energy can build up on timescales of hours to days, and then may be suddenly released in the form of a magnetic eruption, which then propagates through interplanetary space, possibly impacting the Earth's space environment. Can we use the observed evolution of the magnetic and velocity fields in the solar photosphere to model the evolution of the overlying solar coronal field, including the storage and release of magnetic energy in such eruptions? The objective of CGEM, the Coronal Global Evolutionary Model, funded by the NASA/NSF Space Weather Modeling program, is to develop and evaluate such a model for the evolution of the coronal magnetic field. The evolving coronal magnetic field can then be used as a starting point for magnetohydrodynamic (MHD) models of the corona, which can then be used to drive models of heliospheric evolution and predictions of magnetic field and plasma density conditions at 1AU.
0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv
Thermal Diagnostics with the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory: A Validated Method for Differential Emission Measure Inversions
We present a new method for performing differential emission measure (DEM) inversions on narrow-band EUV images from the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO). The method yields positive definite DEM solutions by solving a linear program. This method has been validated against a diverse set of thermal models of varying complexity and realism. These include (1) idealized gaussian DEM distributions, (2) 3D models of NOAA Active Region 11158 comprising quasi-steady loop atmospheres in a non-linear force-free field, and (3) thermodynamic models from a fully-compressible, 3D MHD simulation of AR corona formation following magnetic flux emergence. We then present results from the application of the method to AIA observations of Active Region 11158, comparing the region's thermal structure on two successive solar rotations. Additionally, we show how the DEM inversion method can be adapted to simultaneously invert AIA and XRT data, and how supplementing AIA data with the latter improves the inversion result. The speed of the method allows for routine production of DEM maps, thus facilitating science studies that require tracking of the thermal structure of the solar corona in time and space.
0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv
A Method for Data-Driven Simulations of Evolving Solar Active Regions
We present a method for performing data-driven simulations of solar active region formation and evolution. The approach is based on magnetofriction, which evolves the induction equation assuming the plasma velocity is proportional to the Lorentz force. The simulations of active region coronal field are driven by temporal sequences of photospheric magnetograms from the Helioseismic Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO). Under certain conditions, the data-driven simulations produce flux ropes that are ejected from the modeled active region due to loss of equilibrium. Following the ejection of flux ropes, we find an enhancement of the photospheric horizontal field near the polarity inversion line. We also present a method for the synthesis of mock coronal images based on a proxy emissivity calculated from the current density distribution in the model. This method yields mock coronal images that are somewhat reminiscent of images of active regions taken by instruments such as SDO's Atmospheric Imaging Assembly (AIA) at extreme ultraviolet wavelengths.
0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv
Magnetohydrodynamics of the Weakly Ionized Solar Photosphere
We investigate the importance of ambipolar diffusion and Hall currents for high-resolution comprehensive ('realistic') photospheric simulations. To do so we extended the radiative magnetohydrodynamics code \emph{MURaM} to use the generalized Ohm's law under the assumption of local thermodynamic equilibrium. We present test cases comparing analytical solutions with numerical simulations for validation of the code. Furthermore, we carried out a number of numerical experiments to investigate the impact of these neutral-ion effects in the photosphere. We find that, at the spatial resolutions currently used (5-20 km per grid point), the Hall currents and ambipolar diffusion begin to become significant -- with flows of 100 m/s in sunspot light bridges, and changes of a few percent in the thermodynamic structure of quiet-Sun magnetic features. The magnitude of the effects is expected to increase rapidly as smaller-scale variations are resolved by the simulations.
0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv
A Critical Assessment of Nonlinear Force-Free Field Modeling of the Solar Corona for Active Region 10953
Nonlinear force-free field (NLFFF) models are thought to be viable tools for investigating the structure, dynamics and evolution of the coronae of solar active regions. In a series of NLFFF modeling studies, we have found that NLFFF models are successful in application to analytic test cases, and relatively successful when applied to numerically constructed Sun-like test cases, but they are less successful in application to real solar data. Different NLFFF models have been found to have markedly different field line configurations and to provide widely varying estimates of the magnetic free energy in the coronal volume, when applied to solar data. NLFFF models require consistent, force-free vector magnetic boundary data. However, vector magnetogram observations sampling the photosphere, which is dynamic and contains significant Lorentz and buoyancy forces, do not satisfy this requirement, thus creating several major problems for force-free coronal modeling efforts. In this article, we discuss NLFFF modeling of NOAA Active Region 10953 using Hinode/SOT-SP, Hinode/XRT, STEREO/SECCHI-EUVI, and SOHO/MDI observations, and in the process illustrate the three such issues we judge to be critical to the success of NLFFF modeling: (1) vector magnetic field data covering larger areas are needed so that more electric currents associated with the full active regions of interest are measured, (2) the modeling algorithms need a way to accommodate the various uncertainties in the boundary data, and (3) a more realistic physical model is needed to approximate the photosphere-to-corona interface in order to better transform the forced photospheric magnetograms into adequate approximations of nearly force-free fields at the base of the corona. We make recommendations for future modeling efforts to overcome these as yet unsolved problems.
0 citations0 reviewsNo code linkedOpen access