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.
Researcher profile
Mark C. M. Cheung
Mark C. M. Cheung contributes to research discovery and scholarly infrastructure.
ResearcherAffiliation not importedOpen to collaborate
Trust snapshot
Quick read
Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
8topics
4close collaborators
Actions
Decide how to stay connected
Identity and collaboration
How to connect with this researcher
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Open a focused conversation when the fit is right
Claim this author entity first to unlock direct invitations.
Research graph
See the researcher in context
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
Building this graph slice
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
13 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
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