Showing 1–50 of 63 results for author: Cheung, M C

Resolving moving heliospheric structures using interplanetary scintillation observations with the Murchison Widefield Array

Authors: A. Waszewski, J. S. Morgan, R. Chhetri, R. Ekers, M. C. M. Cheung, N. D. R Bhat, M. Johnston-Hollitt

Abstract: We have conducted a blind search in 49 consecutive days of interplanetary scintillation observations made by the Murchison Widefield Array from mid-2019, with overlapping daily observations approximately East and South-East of the Sun at an elongation of $\sim$30 degrees and a field of view of 30 degrees. These observations detect an unprecedented density of sources. In spite of these observations… ▽ More We have conducted a blind search in 49 consecutive days of interplanetary scintillation observations made by the Murchison Widefield Array from mid-2019, with overlapping daily observations approximately East and South-East of the Sun at an elongation of $\sim$30 degrees and a field of view of 30 degrees. These observations detect an unprecedented density of sources. In spite of these observations being taken at sunspot minimum, this search has revealed several interesting transitory features characterised by elevated scintillation levels. One solar wind enhancement is captured in two observations several hours apart, allowing its radial movement away from the Sun to be measured. We present here a methodology for measuring the plane-of-sky velocity for the moving heliospheric structure. The plane-of-sky velocity was inferred as $0.66\pm0.147\,^{\text{o}}\text{hr}^{-1}$, or $480\pm106\,\text{km}\,\text{s}^{-1}$ assuming a distance of 1AUえーゆー. After cross-referencing our observed structure with multiple catalogues of heliospheric events, we propose that the likely source of our observed structure is a stream-interaction region originating from a low-latitude coronal hole. This work demonstrates the power of widefield interplanetary scintillation observations to capture detailed features in the heliosphere which are otherwise unresolvable and go undetected. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: 14 pages, 5 figures

Machine learning in solar physics

Authors: A. Asensio Ramos, M. C. M. Cheung, I. Chifu, R. Gafeira

Abstract: The application of machine learning in solar physics has the potential to greatly enhance our understanding of the complex processes that take place in the atmosphere of the Sun. By using techniques such as deep learning, we are now in the position to analyze large amounts of data from solar observations and identify patterns and trends that may not have been apparent using traditional methods. Th… ▽ More The application of machine learning in solar physics has the potential to greatly enhance our understanding of the complex processes that take place in the atmosphere of the Sun. By using techniques such as deep learning, we are now in the position to analyze large amounts of data from solar observations and identify patterns and trends that may not have been apparent using traditional methods. This can help us improve our understanding of explosive events like solar flares, which can have a strong effect on the Earth environment. Predicting hazardous events on Earth becomes crucial for our technological society. Machine learning can also improve our understanding of the inner workings of the sun itself by allowing us to go deeper into the data and to propose more complex models to explain them. Additionally, the use of machine learning can help to automate the analysis of solar data, reducing the need for manual labor and increasing the efficiency of research in this field. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 100 pages, 13 figures, 286 references, accepted for publication as a Living Review in Solar Physics (LRSP)

Comprehensive radiative MHD simulations of eruptive flares above collisional polarity inversion lines

Authors: Matthias Rempel, Georgios Chintzoglou, Mark C. M. Cheung, Yuhong Fan, Lucia Kleint

Abstract: We present a new simulation setup using the MURaM radiative Magnetohydrodynamic (MHD) code that allows to study the formation of collisional polarity inversion lines (cPILs) in the photosphere and the coronal response including flares. In this scheme we start with a bipolar sunspot configuration and set the spots on collision course by imposing the appropriate velocity field at the footpoints in t… ▽ More We present a new simulation setup using the MURaM radiative Magnetohydrodynamic (MHD) code that allows to study the formation of collisional polarity inversion lines (cPILs) in the photosphere and the coronal response including flares. In this scheme we start with a bipolar sunspot configuration and set the spots on collision course by imposing the appropriate velocity field at the footpoints in the subphotospheric boundary. We produce different setups with the same initial spot separation by varying physical parameters such as the collision speed and minimum collision distance. While all setups lead to the formation of an EUV and X-ray sigmoid structure, only the cases with a close passing of the spots cause flares and mass eruptions. The energy release is in the $1-2 \times 10^{31}$ erg range, putting the simulated flares into the upper C-class to lower M-class range of GOES X-ray 1-8 A flux. While the setup with the more distant passing of the spots does not lead to a flare, the corona is nonetheless substantially heated, suggesting non-eruptive energy release mechanisms. We focus our discussion on two particular setups that differ in spot coherence and resulting cPIL length persistence. We find different timings in the transition from a sheared magnetic arcade (SMA) to magnetic flux rope (MFR); the setup with a large length but shorter duration cPIL produces a MFR during the eruption, while the MFR is pre-existing in the setup with a large length and longer duration cPIL. While both result in flares of comparable strength and the eruption of a CME, the setup with pre-existing MFR (and embedded filament) leads to an MFR eruption with a larger mass content. △ Less

Submitted 25 August, 2023; v1 submitted 9 March, 2023; originally announced March 2023.

Comments: 19 pages, 11 figures, 1 table, ApJ in press (replaced previous version with the accepted article)

Surface Flux Transport on the Sun

Authors: Anthony R. Yeates, Mark C. M. Cheung, Jie Jiang, Kristof Petrovay, Yi-Ming Wang

Abstract: We review the surface flux transport model for the evolution of magnetic flux patterns on the Sun's surface. Our underlying motivation is to understand the model's prediction of the polar field (or axial dipole) strength at the end of the solar cycle. The main focus is on the "classical" model: namely, steady axisymmetric profiles for differential rotation and meridional flow, and uniform supergra… ▽ More We review the surface flux transport model for the evolution of magnetic flux patterns on the Sun's surface. Our underlying motivation is to understand the model's prediction of the polar field (or axial dipole) strength at the end of the solar cycle. The main focus is on the "classical" model: namely, steady axisymmetric profiles for differential rotation and meridional flow, and uniform supergranular diffusion. Nevertheless, the review concentrates on recent advances, notably in understanding the roles of transport parameters and - in particular - the source term. We also discuss the physical justification for the surface flux transport model, along with efforts to incorporate radial diffusion, and conclude by summarizing the main directions where researchers have moved beyond the classical model. △ Less

Submitted 11 May, 2023; v1 submitted 2 March, 2023; originally announced March 2023.

Comments: 35 pages, 11 figures, accepted for publication in Space Science Reviews

Implications of Different Solar Photospheric Flux-Transport Models for Global Coronal and Heliospheric Modeling

Authors: Graham Barnes, Marc L. DeRosa, Shaela I. Jones, Charles N. Arge, Carl J. Henney, Mark C. M. Cheung

Abstract: The concept of surface-flux transport (SFT) is commonly used in evolving models of the large-scale solar surface magnetic field. These photospheric models are used to determine the large-scale structure of the overlying coronal magnetic field, as well as to make predictions about the fields and flows that structure the solar wind. We compare predictions from two SFT models for the solar wind, open… ▽ More The concept of surface-flux transport (SFT) is commonly used in evolving models of the large-scale solar surface magnetic field. These photospheric models are used to determine the large-scale structure of the overlying coronal magnetic field, as well as to make predictions about the fields and flows that structure the solar wind. We compare predictions from two SFT models for the solar wind, open magnetic field footpoints, and the presence of coronal magnetic null points throughout various phases of a solar activity cycle, focusing on the months of April in even-numbered years between 2012 and 2020, inclusive. We find that there is a solar cycle dependence to each of the metrics considered, but there is not a single phase of the cycle in which all the metrics indicate good agreement between the models. The metrics also reveal large, transient differences between the models when a new active region is rotating into the assimilation window. The evolution of the surface flux is governed by a combination of large scale flows and comparatively small scale motions associated with convection. Because the latter flows evolve rapidly, there are intervals during which their impact on the surface flux can only be characterized in a statistical sense, thus their impact is modeled by introducing a random evolution that reproduces the typical surface flux evolution. We find that the differences between the predicted properties are dominated by differences in the model assumptions and implementation, rather than selection of a particular realization of the random evolution. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: Accepted for publication in The Astrophysical Journal

Data-Driven Radiative Magnetohydrodynamics Simulations with the MURaM code

Authors: Feng Chen, Mark C. M. Cheung, Matthias Rempel, Georgios Chintzoglou

Abstract: We present a method of conducting data-driven simulations of solar active regions and flux emergence with the MURaM radiative magnetohydrodynamics (MHD) code. The horizontal electric field derived from the full velocity and magnetic vectors, is implemented at the photospheric (bottom) boundary to drive the induction equation. The energy equation accounts for thermal conduction along magnetic field… ▽ More We present a method of conducting data-driven simulations of solar active regions and flux emergence with the MURaM radiative magnetohydrodynamics (MHD) code. The horizontal electric field derived from the full velocity and magnetic vectors, is implemented at the photospheric (bottom) boundary to drive the induction equation. The energy equation accounts for thermal conduction along magnetic fields, optically-thin radiative loss, and heating of coronal plasma by viscous and resistive dissipation, which allows for a realistic presentation of the thermodynamic properties of coronal plasma that are key to predicting the observational features of solar active regions and eruptions. To validate the method, the photospheric data from a comprehensive radiative MHD simulation of solar eruption (the ground truth) are used to drive a series of numerical experiments. The data-driven simulation reproduces the accumulation of free magnetic energy over the course of flux emergence in the ground truth with an error of 3\%. The onset time is approximately 8\,min delayed compared to the ground truth. However, a precursor-like signature can be identified at the correct onset time. The data-driven simulation captures key eruption-related emission features and plasma dynamics of the ground truth flare over a wide temperature span from $\log_{10}T{=}4.5$ to $\log_{10}T{>}8$. The evolution of the flare and coronal mass ejection as seen in synthetic extreme ultraviolet images is also reproduced with high fidelity. The method helps to understand the evolution of magnetic field in a more realistic coronal environment and to link the magnetic structures to observable diagnostics. △ Less

Submitted 29 March, 2023; v1 submitted 18 January, 2023; originally announced January 2023.

Comments: 35 pages, 24 figures, accepted for publication in ApJ

Science Platforms for Heliophysics Data Analysis

Authors: Monica G. Bobra, Will T. Barnes, Thomas Y. Chen, Mark C. M. Cheung, Laura A. Hayes, Jack Ireland, Miho Janvier, Michael S. F. Kirk, James P. Mason, Stuart J. Mumford, Paul J. Wright

Abstract: 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. 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. △ Less

Submitted 2 January, 2023; originally announced January 2023.

Comments: Heliophysics 2050 White Paper

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

Authors: M. Rempel, Y. Fan, M. Dikpati, A. Malanushenko, M. D. Kazachenko, M. C. M. Cheung, G. Chintzoglou, X. Sun, G. H. Fisher, T. Y. Chen

Abstract: 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 orientatio… ▽ More 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. △ Less

Submitted 29 December, 2022; originally announced December 2022.

Comments: Heliophysics 2050 White Paper

Signatures of dynamic fibrils at the coronal base: Observations from Solar Orbiter/EUI

Authors: Sudip Mandal, Hardi Peter, Lakshmi Pradeep Chitta, Regina A. Cuadrado, Udo Schühle, Luca Teriaca, Sami K. Solanki, Louise Harra, David Berghmans, Frédéric Auchère, Susanna Parenti, Andrei N. Zhukov, Éric Buchlin, Cis Verbeeck, Emil Kraaikamp, Luciano Rodriguez, David M. Long, Conrad Schwanitz, Krzysztof Barczynski, Gabriel Pelouze, Philip J. Smith, Wei Liu, Mark C. Cheung

Abstract: The solar chromosphere hosts a wide variety of transients, including dynamic fibrils (DFs) that are characterised as elongated, jet-like features seen in active regions, often through H$αあるふぁ$ diagnostics. So far, these features have been difficult to identify in coronal images primarily due to their small size and the lower spatial resolution of the current EUV imagers. Here we present the first unam… ▽ More The solar chromosphere hosts a wide variety of transients, including dynamic fibrils (DFs) that are characterised as elongated, jet-like features seen in active regions, often through H$αあるふぁ$ diagnostics. So far, these features have been difficult to identify in coronal images primarily due to their small size and the lower spatial resolution of the current EUV imagers. Here we present the first unambiguous signatures of DFs in coronal EUV data using high-resolution images from the Extreme Ultraviolet Imager (EUI) on board Solar Orbiter. Using the data acquired with the 174~Å High Resolution Imager (HRI$_{EUV}$) of EUI, we find many bright dot-like features (of size 0.3-0.5 Mm) that move up and down (often repeatedly) in the core of an active region. In a space-time map, these features produce parabolic tracks akin to the chromospheric observations of DFs. Properties such as their speeds (14 km~s$^{-1}$), lifetime (332~s), deceleration (82 m~s$^{-2}$) and lengths (1293~km) are also reminiscent of the chromospheric DFs. The EUI data strongly suggest that these EUV bright dots are basically the hot tips (of the cooler chromospheric DFs) that could not be identified unambiguously before because of a lack of spatial resolution. △ Less

Submitted 9 December, 2022; originally announced December 2022.

Comments: Accepted for publication in A&A Letters. Event movie can be downloaded from https://drive.google.com/file/d/1o_4jHA5JbyQtrpUBtB3ItE_s3HjF6ncc/view?usp=sharing

Journal ref: A&A 670, L3 (2023)

Solar Orbiter and SDO Observations, and Bifrost MHD Simulations of Small-scale Coronal Jets

Authors: Navdeep K. Panesar, Viggo H. Hansteen, Sanjiv K. Tiwari, Mark C. M. Cheung, David Berghmans, Daniel Müller

Abstract: We report high-resolution, high-cadence observations of five small-scale coronal jets in an on-disk quiet Sun region observed with Solar Orbiter's EUI/\hri\ in 174 Å. We combine the \hri\ images with the EUV images of SDO/AIA and investigate magnetic setting of the jets using co-aligned line-of-sight magnetograms from SDO/HMI. The \hri\ jets are miniature versions of typical coronal jets as they s… ▽ More We report high-resolution, high-cadence observations of five small-scale coronal jets in an on-disk quiet Sun region observed with Solar Orbiter's EUI/\hri\ in 174 Å. We combine the \hri\ images with the EUV images of SDO/AIA and investigate magnetic setting of the jets using co-aligned line-of-sight magnetograms from SDO/HMI. The \hri\ jets are miniature versions of typical coronal jets as they show narrow collimated spires with a base brightening. Three out of five jets result from a detectable minifilament eruption following flux cancelation at the neutral line under the minifilament, analogous to coronal jets. To better understand the physics of jets, we also analyze five small-scale jets from a high-resolution Bifrost MHD simulation in synthetic \FeIX/\FeX\ emissions. The jets in the simulation reside above neutral lines and four out of five jets are triggered by magnetic flux cancelation. The temperature maps show the evidence of cool gas in the same four jets. Our simulation also shows the signatures of opposite Doppler shifts (of the order of $\pm$10s of \kms) in the jet spire, which is evidence of untwisting motion of the magnetic field in the jet spire. The average jet duration, spire length, base width, and speed in our observations (and in synthetic \FeIX/\FeX\ images) are 6.5$\pm$4.0 min (9.0$\pm$4.0 min), 6050$\pm$2900 km (6500$\pm$6500 km), 2200$\pm$850 km, (3900$\pm$2100 km), and 60$\pm$8 \kms\ (42$\pm$20 \kms), respectively. Our observation and simulation results provide a unified picture of small-scale solar coronal jets driven by magnetic reconnection accompanying flux cancelation. This picture also aligns well with the most recent reports of the formation and eruption mechanisms of larger coronal jets. △ Less

Submitted 11 November, 2022; originally announced November 2022.

Comments: 15 pages, 9 Figures, 1 Table, Accepted for publication in ApJ

Solar coronal heating from small-scale magnetic braids

Authors: L. P. Chitta, H. Peter, S. Parenti, D. Berghmans, F. Auchère, S. K. Solanki, R. Aznar Cuadrado, U. Schühle, L. Teriaca, S. Mandal, K. Barczynski, É. Buchlin, L. Harra, E. Kraaikamp, D. M. Long, L. Rodriguez, C. Schwanitz, P. J. Smith, C. Verbeeck, A. N. Zhukov, W. Liu, M. C. M. Cheung

Abstract: Relaxation of braided coronal magnetic fields through reconnection is thought to be a source of energy to heat plasma in active region coronal loops. However, observations of active region coronal heating associated with an untangling of magnetic braids remain sparse. One reason for this paucity could be the lack of coronal observations with a sufficiently high spatial and temporal resolution to c… ▽ More Relaxation of braided coronal magnetic fields through reconnection is thought to be a source of energy to heat plasma in active region coronal loops. However, observations of active region coronal heating associated with an untangling of magnetic braids remain sparse. One reason for this paucity could be the lack of coronal observations with a sufficiently high spatial and temporal resolution to capture this process in action. Using new observations with high spatial resolution (250-270 km on the Sun) and high cadence (3-10 s) from the Extreme Ultraviolet Imager (EUI) on board Solar Orbiter, we observed the untangling of small-scale coronal braids in different active regions. The untangling is associated with impulsive heating of the gas in these braided loops. We assess that coronal magnetic braids overlying cooler chromospheric filamentary structures are perhaps more common. Furthermore, our observations show signatures of spatially coherent and intermittent coronal heating during the relaxation of the magnetic braids. Our study reveals the operation of gentle and impulsive modes of magnetic reconnection in the solar corona. △ Less

Submitted 26 November, 2022; v1 submitted 25 September, 2022; originally announced September 2022.

Comments: Published in Astronomy & Astrophysics

Journal ref: A&A 667, A166 (2022)

Exploring the Limits of Synthetic Creation of Solar EUV Images via Image-to-Image Translation

Authors: Valentina Salvatelli, Luiz F. G. dos Santos, Souvik Bose, Brad Neuberg, Mark C. M. Cheung, Miho Janvier, Meng Jin, Yarin Gal, Atilim Gunes Baydin

Abstract: The Solar Dynamics Observatory (SDO), a NASA multi-spectral decade-long mission that has been daily producing terabytes of observational data from the Sun, has been recently used as a use-case to demonstrate the potential of machine learning methodologies and to pave the way for future deep-space mission planning. In particular, the idea of using image-to-image translation to virtually produce ext… ▽ More The Solar Dynamics Observatory (SDO), a NASA multi-spectral decade-long mission that has been daily producing terabytes of observational data from the Sun, has been recently used as a use-case to demonstrate the potential of machine learning methodologies and to pave the way for future deep-space mission planning. In particular, the idea of using image-to-image translation to virtually produce extreme ultra-violet channels has been proposed in several recent studies, as a way to both enhance missions with less available channels and to alleviate the challenges due to the low downlink rate in deep space. This paper investigates the potential and the limitations of such a deep learning approach by focusing on the permutation of four channels and an encoder--decoder based architecture, with particular attention to how morphological traits and brightness of the solar surface affect the neural network predictions. In this work we want to answer the question: can synthetic images of the solar corona produced via image-to-image translation be used for scientific studies of the Sun? The analysis highlights that the neural network produces high-quality images over three orders of magnitude in count rate (pixel intensity) and can generally reproduce the covariance across channels within a 1% error. However the model performance drastically diminishes in correspondence of extremely high energetic events like flares, and we argue that the reason is related to the rareness of such events posing a challenge to model training. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Comments: 16 pages, 8 figures. To be published on ApJ (submitted on Feb 21st, accepted on July 28th)

Journal ref: ApJ 937 (2022) 100

Ultraviolet Spectropolarimetry With Polstar: Using Polstar to test Magnetospheric Mass-loss Quenching

Authors: M. E. Shultz, R. Casini, M. C. M. Cheung, A. David-Uraz, T. del Pino Alemán, C. Erba, C. P. Folsom, K. Gayley, R. Ignace, Z. Keszthelyi, O. Kochukhov, Y. Nazé, C. Neiner, M. Oksala, V. Petit, P. A. Scowen, N. Sudnik, A. ud-Doula, J. S. Vink, G. A. Wade

Abstract: Polstar is a proposed NASA MIDEX space telescope that will provide high-resolution, simultaneous full-Stokes spectropolarimetry in the far ultraviolet, together with low-resolution linear polarimetry in the near ultraviolet. This observatory offers unprecedented capabilities to obtain unique information on the magnetic and plasma properties of the magnetospheres of hot stars. We describe an observ… ▽ More Polstar is a proposed NASA MIDEX space telescope that will provide high-resolution, simultaneous full-Stokes spectropolarimetry in the far ultraviolet, together with low-resolution linear polarimetry in the near ultraviolet. This observatory offers unprecedented capabilities to obtain unique information on the magnetic and plasma properties of the magnetospheres of hot stars. We describe an observing program making use of the known population of magnetic hot stars to test the fundamental hypothesis that magnetospheres should act to rapidly drain angular momentum, thereby spinning the star down, whilst simultaneously reducing the net mass-loss rate. Both effects are expected to lead to dramatic differences in the evolution of magnetic vs. non-magnetic stars. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: 20 pages, 10 figures, accepted for publication in ApSS. arXiv admin note: substantial text overlap with arXiv:2111.06434

Ultraviolet Spectropolarimetric Diagnostics of Hot Star Magnetospheres

Authors: Asif ud-Doula, M. C. M. Cheung, A. David-Uraz, C. Erba, C. P. Folsom, K. Gayley, Y. Naze, C. Neiner, V. Petit, R. Prinja, M. E. Shultz, N. Sudnik, J. S. Vink, G. A. Wade

Abstract: Several space missions and instruments for UV spectropolarimetry are in preparation, such as the proposed NASA MIDEX Polstar project, the proposed ESA M mission Arago, and the Pollux instrument on the future LUVOIR-like NASA flagship mission. In the frame of Polstar, we have studied the capabilities these observatories would offer to gain information on the magnetic and plasma properties of the ma… ▽ More Several space missions and instruments for UV spectropolarimetry are in preparation, such as the proposed NASA MIDEX Polstar project, the proposed ESA M mission Arago, and the Pollux instrument on the future LUVOIR-like NASA flagship mission. In the frame of Polstar, we have studied the capabilities these observatories would offer to gain information on the magnetic and plasma properties of the magnetospheres of hot stars, helping us test the fundamental hypothesis that magnetospheres should act to rapidly drain angular momentum, thereby spinning the star down, whilst simultaneously reducing the net mass-loss rate. Both effects are expected to lead to dramatic differences in the evolution of magnetic vs. non-magnetic stars. △ Less

Submitted 26 June, 2022; originally announced June 2022.

Comments: Accepted for publication in Astrophysics and Space Science. arXiv admin note: substantial text overlap with arXiv:2111.06434

Report number: ASTR-D-22-00095R3

Global geomagnetic perturbation forecasting using Deep Learning

Authors: Vishal Upendran, Panagiotis Tigas, Banafsheh Ferdousi, Teo Bloch, Mark C. M. Cheung, Siddha Ganju, Asti Bhatt, Ryan M. McGranaghan, Yarin Gal

Abstract: 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 increasi… ▽ More 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. △ Less

Submitted 12 May, 2022; originally announced May 2022.

Comments: 23 pages, 8 figures, 5 tables; accepted for publication in AGU: Spaceweather

Coronal Mass Ejections and Dimmings: A Comparative Study using MHD Simulations and SDO Observations

Authors: Meng Jin, Mark C. M. Cheung, Marc L. DeRosa, Nariaki V. Nitta, Carolus J. Schrijver

Abstract: 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… ▽ More 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. △ Less

Submitted 25 February, 2022; originally announced February 2022.

Comments: 16 pages, 9 figures, accepted for publication in ApJ

Ultraviolet Spectropolarimetry With Polstar: Hot Star Magnetospheres

Authors: M. E. Shultz, R. Casini, M. C. M. Cheung, A. David-Uraz, T. del Pino Alemán, C. Erba, C. P. Folsom, K. Gayley, R. Ignace, Z. Keszthelyi, O. Kochukhov, Y. Nazé, C. Neiner, M. Oksala, V. Petit, P. A. Scowen, N. Sudnik, A. ud-Doula, J. S. Vink, G. A. Wade

Abstract: Polstar is a proposed NASA MIDEX space telescope that will provide high-resolution, simultaneous full-Stokes spectropolarimetry in the far ultraviolet, together with low-resolution linear polarimetry in the near ultraviolet. In this white paper, we describe the unprecedented capabilities this observatory would offer in order to obtain unique information on the magnetic and plasma properties of the… ▽ More Polstar is a proposed NASA MIDEX space telescope that will provide high-resolution, simultaneous full-Stokes spectropolarimetry in the far ultraviolet, together with low-resolution linear polarimetry in the near ultraviolet. In this white paper, we describe the unprecedented capabilities this observatory would offer in order to obtain unique information on the magnetic and plasma properties of the magnetospheres of hot stars. This would enable a test of the fundamental hypothesis that magnetospheres should act to rapidly drain angular momentum, thereby spinning the star down, whilst simultaneously reducing the net mass-loss rate. Both effects are expected to lead to dramatic differences in the evolution of magnetic vs. non-magnetic stars. △ Less

Submitted 9 December, 2021; v1 submitted 11 November, 2021; originally announced November 2021.

Comments: White paper, 40 pages

The solar internetwork. III. Unipolar versus bipolar flux appearance

Authors: Milan Gošić, Luis R. Bellot Rubio, Mark C. M. Cheung, David Orozco Suárez, Yukio Katsukawa, Jose Carlos Del Toro Iniesta

Abstract: 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 i… ▽ More 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. △ Less

Submitted 4 November, 2021; originally announced November 2021.

Comments: 14 pages, 11 figures. Accepted for publication in ApJ. Animations are available at https://www.lmsal.com/~mgosic/download/animations_apj_2021b.tgz

The Magnetic Origin of Solar Campfires

Authors: Navdeep K. Panesar, Sanjiv K. Tiwari, David Berghmans, Mark C. M. Cheung, Daniel Muller, Frederic Auchere, Andrei Zhukov

Abstract: Solar campfires are fine-scale heating events, recently observed by Extreme Ultraviolet Imager (EUI), onboard Solar Orbiter. Here we use EUI 174Å images, together with EUV images from SDO/AIA, and line-of-sight magnetograms from SDO/HMI to investigate the magnetic origin of 52 randomly selected campfires in the quiet solar corona. We find that (i) the campfires are rooted at the edges of photosphe… ▽ More Solar campfires are fine-scale heating events, recently observed by Extreme Ultraviolet Imager (EUI), onboard Solar Orbiter. Here we use EUI 174Å images, together with EUV images from SDO/AIA, and line-of-sight magnetograms from SDO/HMI to investigate the magnetic origin of 52 randomly selected campfires in the quiet solar corona. We find that (i) the campfires are rooted at the edges of photospheric magnetic network lanes; (ii) most of the campfires reside above the neutral line between majority-polarity magnetic flux patch and a merging minority-polarity flux patch, with a flux cancelation rate of $\sim$10$^{18}$Mx hr$^{-1}$; (iii) some of the campfires occur repeatedly from the same neutral line; (iv) in the large majority of instances, campfires are preceded by a cool-plasma structure, analogous to minifilaments in coronal jets; and (v) although many campfires have `complex' structure, most campfires resemble small-scale jets, dots, or loops. Thus, `campfire' is a general term that includes different types of small-scale solar dynamic features. They contain sufficient magnetic energy ($\sim$10$^{26}$-10$^{27}$ erg) to heat the solar atmosphere locally to 0.5--2.5MK. Their lifetimes range from about a minute to over an hour, with most of the campfires having a lifetime of $<$10 minutes. The average lengths and widths of the campfires are 5400$\pm$2500km and 1600$\pm$640km, respectively. Our observations suggest that (a) the presence of magnetic flux ropes may be ubiquitous in the solar atmosphere and not limited to coronal jets and larger-scale eruptions that make CMEs, and (b) magnetic flux cancelation is the fundamental process for the formation and triggering of most campfires. △ Less

Submitted 13 October, 2021; originally announced October 2021.

Comments: Accepted for publication in ApJ Letters, 20 Pages, 1 Table, 12 Figures

SynthIA: A Synthetic Inversion Approximation for the Stokes Vector Fusing SDO and Hinode into a Virtual Observatory

Authors: Richard E. L. Higgins, David F. Fouhey, Spiro K. Antiochos, Graham Barnes, Mark C. M. Cheung, J. Todd Hoeksema, KD Leka, Yang Liu, Peter W. Schuck, Tamas I. Gombosi

Abstract: 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 an… ▽ More 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. △ Less

Submitted 27 August, 2021; originally announced August 2021.

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

Authors: Mark C. M. Cheung, Juan Martínez-Sykora, Paola Testa, Bart De Pontieu, Georgios Chintzoglou, Matthias Rempel, Vanessa Polito, Graham S. Kerr, Katharine K. Reeves, Lyndsay Fletcher, Meng Jin, Daniel Nóbrega-Siverio, Sanja Danilovic, Patrick Antolin, Joel Allred, Viggo Hansteen, Ignacio Ugarte-Urra, Edward DeLuca, Dana Longcope, Shinsuke Takasao, Marc DeRosa, Paul Boerner, Sarah Jaeggli, Nariaki Nitta, Adrian Daw , et al. (3 additional authors not shown)

Abstract: 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 at… ▽ More 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. △ Less

Submitted 30 June, 2021; v1 submitted 29 June, 2021; originally announced June 2021.

Comments: 42 pages, 22 figures, submitted to ApJ

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

Authors: Bart De Pontieu, Paola Testa, Juan Martinez-Sykora, Patrick Antolin, Konstantinos Karampelas, Viggo Hansteen, Matthias Rempel, Mark C. M. Cheung, Fabio Reale, Sanja Danilovic, Paolo Pagano, Vanessa Polito, Ineke De Moortel, Daniel Nobrega-Siverio, Tom Van Doorsselaere, Antonino Petralia, Mahboubeh Asgari-Targhi, Paul Boerner, Mats Carlsson, Georgios Chintzoglou, Adrian Daw, Ed DeLuca, Leon Golub, Takuma Matsumoto, Ignacio Ugarte-Urra , et al. (2 additional authors not shown)

Abstract: 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 t… ▽ More 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. △ Less

Submitted 29 June, 2021; originally announced June 2021.

Comments: 46 pages, 22 figures, submitted to ApJ

The Coronal Veil

Authors: A. Malanushenko, M. C. M. Cheung, C. E. DeForest, J. A. Klimchuk, M. Rempel

Abstract: 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 p… ▽ More 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. △ Less

Submitted 29 November, 2021; v1 submitted 28 June, 2021; originally announced June 2021.

Comments: ApJ (in press)

Plasma Heating Induced by Tadpole-Like Downflows in the Flaring Solar Corona

Authors: Tanmoy Samanta, Hui Tian, Bin Chen, Katharine K. Reeves, Mark C. M. Cheung, Angelos Vourlidas, Dipankar Banerjee

Abstract: As one of the most spectacular energy release events in the solar system, solar flares are generally powered by magnetic reconnection in the solar corona. As a result of the re-arrangement of magnetic field topology after the reconnection process, a series of new loop-like magnetic structures are often formed and are known as flare loops. A hot diffuse region, consisting of around 5-10 MK plasma,… ▽ More As one of the most spectacular energy release events in the solar system, solar flares are generally powered by magnetic reconnection in the solar corona. As a result of the re-arrangement of magnetic field topology after the reconnection process, a series of new loop-like magnetic structures are often formed and are known as flare loops. A hot diffuse region, consisting of around 5-10 MK plasma, is also observed above the loops and is called a supra-arcade fan. Often, dark, tadpole-like structures are seen to descend through the bright supra-arcade fans. It remains unclear what role these so-called supra-arcade downflows (SADs) play in heating the flaring coronal plasma. Here we show a unique flare observation, where many SADs collide with the flare loops and strongly heat the loops to a temperature of 10-20 MK. Several of these interactions generate clear signatures of quasi-periodic enhancement in the full-Sun-integrated soft X-ray emission, providing an alternative interpretation for quasi-periodic pulsations that are commonly observed during solar and stellar flares. △ Less

Submitted 26 March, 2021; originally announced March 2021.

Comments: Published in The Innovation; Animations are available in the Journal's Website

Journal ref: The Innovation 2(1), 100083 (2021)

Global Earth Magnetic Field Modeling and Forecasting with Spherical Harmonics Decomposition

Authors: Panagiotis Tigas, Téo Bloch, Vishal Upendran, Banafsheh Ferdoushi, Mark C. M. Cheung, Siddha Ganju, Ryan M. McGranaghan, Yarin Gal, Asti Bhatt

Abstract: 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 spac… ▽ More 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. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Comments: Third Workshop on Machine Learning and the Physical Sciences (NeurIPS 2020), Vancouver, Canada

Multi-Channel Auto-Calibration for the Atmospheric Imaging Assembly using Machine Learning

Authors: Luiz F. G. dos Santos, Souvik Bose, Valentina Salvatelli, Brad Neuberg, Mark C. M. Cheung, Miho Janvier, Meng Jin, Yarin Gal, Paul Boerner, Atılım Güneş Baydin

Abstract: Solar activity plays a quintessential role in influencing the interplanetary medium and space-weather around the Earth. Remote sensing instruments onboard heliophysics space missions provide a pool of information about the Sun's activity via the measurement of its magnetic field and the emission of light from the multi-layered, multi-thermal, and dynamic solar atmosphere. Extreme UV (EUV) waveleng… ▽ More Solar activity plays a quintessential role in influencing the interplanetary medium and space-weather around the Earth. Remote sensing instruments onboard heliophysics space missions provide a pool of information about the Sun's activity via the measurement of its magnetic field and the emission of light from the multi-layered, multi-thermal, and dynamic solar atmosphere. Extreme UV (EUV) wavelength observations from space help in understanding the subtleties of the outer layers of the Sun, namely the chromosphere and the corona. Unfortunately, such instruments, like the Atmospheric Imaging Assembly (AIA) onboard NASA's Solar Dynamics Observatory (SDO), suffer from time-dependent degradation, reducing their sensitivity. Current state-of-the-art calibration techniques rely on periodic sounding rockets, which can be infrequent and rather unfeasible for deep-space missions. We present an alternative calibration approach based on convolutional neural networks (CNNs). We use SDO-AIA data for our analysis. Our results show that CNN-based models could comprehensively reproduce the sounding rocket experiments' outcomes within a reasonable degree of accuracy, indicating that it performs equally well compared with the current techniques. Furthermore, a comparison with a standard "astronomer's technique" baseline model reveals that the CNN approach significantly outperforms this baseline. Our approach establishes the framework for a novel technique to calibrate EUV instruments and advance our understanding of the cross-channel relation between different EUV channels. △ Less

Submitted 1 February, 2021; v1 submitted 27 December, 2020; originally announced December 2020.

Comments: 12 pages, 7 figures, 8 tables. This is a pre-print of an article submitted and accepted by A&A Journal

Journal ref: A&A 648, A53 (2021)

Non-Neutralized Electric Current of Active Regions Explained as a Projection Effect

Authors: Xudong Sun, Mark C. M. Cheung

Abstract: 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… ▽ More 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. △ Less

Submitted 2 December, 2020; v1 submitted 23 November, 2020; originally announced November 2020.

Comments: 19 pages, 7 figures. Accepted for publication in Solar Physics. Minor updates from v1

Sun-as-a-star Spectral Irradiance Observations of Transiting Active Regions

Authors: Shin Toriumi, Vladimir S. Airapetian, Hugh S. Hudson, Carolus J. Schrijver, Mark C. M. Cheung, Marc L. DeRosa

Abstract: Major solar flares are prone to occur in active region atmospheres associated with large, complex, dynamically-evolving sunspots. This points to the importance of monitoring the evolution of starspots, not only in visible but also in ultra violet (UV) and X-rays, in understanding the origin and occurrence of stellar flares. To this end, we perform spectral irradiance analysis on different types of… ▽ More Major solar flares are prone to occur in active region atmospheres associated with large, complex, dynamically-evolving sunspots. This points to the importance of monitoring the evolution of starspots, not only in visible but also in ultra violet (UV) and X-rays, in understanding the origin and occurrence of stellar flares. To this end, we perform spectral irradiance analysis on different types of transiting solar active regions by using a variety of full-disk synoptic observations. The target events are an isolated sunspot, spotless plage, and emerging flux in prolonged quiet-Sun conditions selected from the past decade. We find that the visible continuum and total solar irradiance become darkened when the spot is at the central meridian, whereas it is bright near the solar limb; UV bands sensitive to the chromosphere correlate well with the variation of total unsigned magnetic flux in the photosphere; amplitudes of EUV and soft X-ray increase with the characteristic temperature, whose light curves are flat-topped due to their sensitivity to the optically thin corona; the transiting spotless plage does not show the darkening in the visible irradiance, while the emerging flux produces an asymmetry in all light curves about the central meridian. The multi-wavelength sun-as-a-star study described here indicates that such time lags between the coronal and photospheric light curves have the potential to probe the extent of coronal magnetic fields above the starspots. In addition, EUV wavelengths that are sensitive to the transition-region temperature sometimes show anti-phased variations, which may be used for diagnosing plasmas around starspots. △ Less

Submitted 10 August, 2020; originally announced August 2020.

Comments: 20 pages, 11 figures, 3 tables, accepted for publication in ApJ

The Coronal Global Evolutionary Model: Using HMI Vector Magnetogram and Doppler Data to Determine Coronal Magnetic Field Evolution

Authors: J. Todd Hoeksema, William P. Abbett, David J. Bercik, Mark C. M. Cheung, Marc L. DeRosa, George H. Fisher, Keiji Hayashi, Maria D. Kazachenko, Yang Liu, Erkka Lumme, Benjamin J. Lynch, Xudong Sun, Brian T. Welsch

Abstract: The Coronal Global Evolutionary Model (CGEM) provides data-driven simulations of the magnetic field in the solar corona to better understand the build-up of magnetic energy that leads to eruptive events. The CGEM project has developed six capabilities. CGEM modules (1) prepare time series of full-disk vector magnetic field observations to (2) derive the changing electric field in the solar photosp… ▽ More The Coronal Global Evolutionary Model (CGEM) provides data-driven simulations of the magnetic field in the solar corona to better understand the build-up of magnetic energy that leads to eruptive events. The CGEM project has developed six capabilities. CGEM modules (1) prepare time series of full-disk vector magnetic field observations to (2) derive the changing electric field in the solar photosphere over active-region scales. This local electric field is (3) incorporated into a surface flux transport model that reconstructs a global electric field that evolves magnetic flux in a consistent way. These electric fields drive a (4) 3D spherical magneto-frictional (SMF) model, either at high-resolution over a restricted range of solid angle or at lower resolution over a global domain, to determine the magnetic field and current density in the low corona. An SMF-generated initial field above an active region and the evolving electric field at the photosphere are used to drive (5) detailed magneto-hydrodynamic (MHD) simulations of active regions in the low corona. SMF or MHD solutions are then used to compute emissivity proxies that can be compared with coronal observations. Finally, a lower-resolution SMF magnetic field is used to initialize (6) a global MHD model that is driven by an SMF electric-field time series to simulate the outer corona and heliosphere, ultimately connecting Sun to Earth. As a demonstration, this report features results of CGEM applied to observations of the evolution of NOAA Active Region 11158 in February 2011. △ Less

Submitted 2 October, 2020; v1 submitted 25 June, 2020; originally announced June 2020.

Comments: 19 pages, 8 figures

Journal ref: Astrophysical Journal Supplement Series 250:28 (15pp), 2020 October

Solar wind prediction using deep learning

Authors: Vishal Upendran, M. C. M Cheung, Shravan Hanasoge, Ganapathi Krishnamurthi

Abstract: Emanating from the base of the Sun's corona, the solar wind fills the interplanetary medium with a magnetized stream of charged particles whose interaction with the Earth's magnetosphere has space-weather consequences such as geomagnetic storms. Accurately predicting the solar wind through measurements of the spatio-temporally evolving conditions in the solar atmosphere is important but remains an… ▽ More Emanating from the base of the Sun's corona, the solar wind fills the interplanetary medium with a magnetized stream of charged particles whose interaction with the Earth's magnetosphere has space-weather consequences such as geomagnetic storms. Accurately predicting the solar wind through measurements of the spatio-temporally evolving conditions in the solar atmosphere is important but remains an unsolved problem in heliophysics and space-weather research. In this work, we use deep learning for prediction of solar wind (SW) properties. We use Extreme Ultraviolet images of the solar corona from space based observations to predict the SW speed from the NASA OMNIWEB dataset, measured at Lagragian point 1. We evaluate our model against autoregressive and naive models, and find that our model outperforms the benchmark models, obtaining a best-fit correlation of 0.55 $\pm$ 0.03 with the observed data. Upon visualization and investigation of how the model uses data to make predictions, we find higher activation at the coronal holes for fast wind prediction ($\approx$ 3 to 4 days prior to prediction), and at the active regions for slow wind prediction. These trends bear an uncanny similarity to the influence of regions potentially being the sources of fast and slow wind, as reported in literature. This suggests that our model was able to learn some of the salient associations between coronal and solar wind structure without built-in physics knowledge. Such an approach may help us discover hitherto unknown relationships in heliophysics datasets. △ Less

Submitted 10 June, 2020; originally announced June 2020.

Comments: 29 pages, 12 figures, accepted for publication in AGU: Spaceweather

Coronal Dimming as a Proxy for Stellar Coronal Mass Ejections

Authors: Meng Jin, Mark C. M. Cheung, Marc L. DeRosa, Nariaki V. Nitta, Carolus J. Schrijver, Kevin France, Adam Kowalski, James P. Mason, Rachel Osten

Abstract: Solar coronal dimmings have been observed extensively in the past two decades and are believed to have close association with coronal mass ejections (CMEs). Recent study found that coronal dimming is the only signature that could differentiate powerful ares that have CMEs from those that do not. Therefore, dimming might be one of the best candidates to observe the stellar CMEs on distant Sun-like… ▽ More Solar coronal dimmings have been observed extensively in the past two decades and are believed to have close association with coronal mass ejections (CMEs). Recent study found that coronal dimming is the only signature that could differentiate powerful ares that have CMEs from those that do not. Therefore, dimming might be one of the best candidates to observe the stellar CMEs on distant Sun-like stars. In this study, we investigate the possibility of using coronal dimming as a proxy to diagnose stellar CMEs. By simulating a realistic solar CME event and corresponding coronal dimming using a global magnetohydrodynamics model (AWSoM: Alfven-wave Solar Model), we first demonstrate the capability of the model to reproduce solar observations. We then extend the model for simulating stellar CMEs by modifying the input magnetic flux density as well as the initial magnetic energy of the CME flux rope. Our result suggests that with improved instrument sensitivity, it is possible to detect the coronal dimming signals induced by the stellar CMEs. △ Less

Submitted 14 February, 2020; originally announced February 2020.

Comments: 8 pages, 4 figures, to appear in the Proceedings of IAU Symposium No. 354 - Solar and Stellar Magnetic Fields: Origins and Manifestations

Comparative Study of Data-driven Solar Coronal Field Models Using a Flux Emergence Simulation as a Ground-truth Data Set

Authors: Shin Toriumi, Shinsuke Takasao, Mark C. M. Cheung, Chaowei Jiang, Yang Guo, Keiji Hayashi, Satoshi Inoue

Abstract: 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 ma… ▽ More 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. △ Less

Submitted 11 January, 2020; originally announced January 2020.

Comments: 17 pages, 11 figures, 1 table, accepted for publication in ApJ. The ground-truth dataset (photospheric slices of the flux emergence simulation) is available online at https://doi.org/10.5281/zenodo.3591984

The PDFI_SS Electric Field Inversion Software

Authors: George H. Fisher, Maria D. Kazachenko, Brian T. Welsch, Xudong Sun, Erkka Lumme, David J. Bercik, Marc L. DeRosa, Mark C. M. Cheung

Abstract: 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 Kazache… ▽ More 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. △ Less

Submitted 24 April, 2020; v1 submitted 17 December, 2019; originally announced December 2019.

Comments: 68 pages, 28 figures, 3 tables. This version (arXiv) will be updated when major changes to the software occur

Journal ref: Astrophysical Journal Supplement Series, 248:2 (52pp) 2020 May

The Multi-slit Approach to Coronal Spectroscopy with the Multi-slit Solar Explorer (MUSE)

Authors: Bart De Pontieu, Juan Martinez-Sykora, Paola Testa, Amy Winebarger, Adrian Daw, Viggo Hansteen, Mark C. M. Cheung, Patrick Antolin, the MUSE team

Abstract: 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… ▽ More 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. △ Less

Submitted 19 September, 2019; originally announced September 2019.

Comments: 51 pages, 17 figures, submitted to Astrophysical Journal

Three-dimensional modeling of chromospheric spectral lines in a simulated active region

Authors: Johan P. Bjørgen, Jorrit Leenaarts, Matthias Rempel, Mark C. M. Cheung, Sanja Danilovic, Jaime de la Cruz Rodríguez, Andrii V. Sukhorukov

Abstract: Because of the complex physics that governs the formation of chromospheric lines, interpretation of solar chromospheric observations is difficult. The origin and characteristics of many chromospheric features are, because of this, unresolved. We focus here on studying two prominent features: long fibrils and flare ribbons. To model them, we use a 3D MHD simulation of an active region which self-co… ▽ More Because of the complex physics that governs the formation of chromospheric lines, interpretation of solar chromospheric observations is difficult. The origin and characteristics of many chromospheric features are, because of this, unresolved. We focus here on studying two prominent features: long fibrils and flare ribbons. To model them, we use a 3D MHD simulation of an active region which self-consistently reproduces both of them. We model the H$αあるふぁ$, Mg II k, Ca II K, and Ca II 8542 Å lines using the 3D non-LTE radiative transfer code Multi3D. This simulation reproduces long fibrils that span between the opposite-polarity sunspots and go up to 4 Mm in height. They can be traced in all lines due to density corrugation. Opposite to previous studies, H$αあるふぁ$, Mg II h&k, and Ca II H&K, are formed at similar height in this model. Magnetic field lines are aligned with the H$αあるふぁ$ fibrils, but the latter holds to a lesser extent for the Ca II 8542 Å line. The simulation shows structures in the H$αあるふぁ$ line core that look like flare ribbons. The emission in the ribbons is caused by a dense chromosphere and a transition region at high column mass. The ribbons are visible in all chromospheric lines, but least prominent in Ca II 8542 Å line. In some pixels, broad asymmetric profiles with a single emission peak are produced, similar to the profiles observed in flare ribbons. They are caused by a deep onset of the chromospheric temperature rise and large velocity gradients. The simulation produces long fibrils similar to what is seen in observations. It also produces structures similar to flare ribbons despite the lack of non-thermal electrons in the simulation. The latter suggests that thermal conduction might be a significant agent in transporting flare energy to the chromosphere in addition to non-thermal electrons. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: 17 pages, 15 figures, accepted for publication by A&A

Journal ref: A&A 631, A33 (2019)

Machine learning reveals systematic accumulation of electric current in lead-up to solar flares

Authors: Dattaraj B. Dhuri, Shravan M. Hanasoge, Mark C. M. Cheung

Abstract: Solar flares - bursts of high-energy radiation responsible for severe space-weather effects - are a consequence of the occasional destabilization of magnetic fields rooted in active regions (ARs). The complexity of AR evolution is a barrier to a comprehensive understanding of flaring processes and accurate prediction. Though machine learning (ML) has been used to improve flare predictions, the pot… ▽ More Solar flares - bursts of high-energy radiation responsible for severe space-weather effects - are a consequence of the occasional destabilization of magnetic fields rooted in active regions (ARs). The complexity of AR evolution is a barrier to a comprehensive understanding of flaring processes and accurate prediction. Though machine learning (ML) has been used to improve flare predictions, the potential for revealing precursors and associated physics has been underexploited. Here, we train ML algorithms to classify between vector-magnetic-field observations from flaring ARs, producing at least one M-/X-class flare, and non-flaring ARs. Analysis of magnetic-field observations accurately classified by the machine presents statistical evidence for (1) ARs persisting in flare-productive states --- characterized by AR area --- for days, before and after M- and X-class flare events, (2) systematic pre-flare build-up of free energy in the form of electric currents, suggesting that associated subsurface magnetic field is twisted, (3) intensification of Maxwell stresses in the corona above newly emerging ARs, days before first flares. These results provide new insights into flare physics and improving flare forecasting. △ Less

Submitted 24 May, 2019; originally announced May 2019.

Comments: Related code is available at https://github.com/dbdhuri/ML-SolarFlare-Analysis

Journal ref: Published in Proceedings of the National Academy of Sciences of the United States of America (PNAS) May, 2019

A Machine Learning Dataset Prepared From the NASA Solar Dynamics Observatory Mission

Authors: Richard Galvez, David F. Fouhey, Meng Jin, Alexandre Szenicer, Andrés Muñoz-Jaramillo, Mark C. M. Cheung, Paul J. Wright, Monica G. Bobra, Yang Liu, James Mason, Rajat Thomas

Abstract: In this paper we present a curated dataset from the NASA Solar Dynamics Observatory (SDO) mission in a format suitable for machine learning research. Beginning from level 1 scientific products we have processed various instrumental corrections, downsampled to manageable spatial and temporal resolutions, and synchronized observations spatially and temporally. We illustrate the use of this dataset w… ▽ More In this paper we present a curated dataset from the NASA Solar Dynamics Observatory (SDO) mission in a format suitable for machine learning research. Beginning from level 1 scientific products we have processed various instrumental corrections, downsampled to manageable spatial and temporal resolutions, and synchronized observations spatially and temporally. We illustrate the use of this dataset with two example applications: forecasting future EVE irradiance from present EVE irradiance and translating HMI observations into AIA observations. For each application we provide metrics and baselines for future model comparison. We anticipate this curated dataset will facilitate machine learning research in heliophysics and the physical sciences generally, increasing the scientific return of the SDO mission. This work is a direct result of the 2018 NASA Frontier Development Laboratory Program. Please see the appendix for access to the dataset. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Accepted to The Astrophysical Journal Supplement Series; 11 pages, 8 figures

Multi-component Decomposition of Astronomical Spectra by Compressed Sensing

Authors: Mark C. M. Cheung, Bart De Pontieu, Juan Martınez-Sykora, Paola Testa, Amy R. Winebarger, Adrian Daw, Viggo Hansteen, Patrick Antolin, Theodore D. Tarbell, Jean-Pierre Wuelser, Peter Young, the MUSE team

Abstract: The signal measured by an astronomical spectrometer may be due to radiation from a multi-component mixture of plasmas with a range of physical properties (e.g. temperature, Doppler velocity). Confusion between multiple components may be exacerbated if the spectrometer sensor is illuminated by overlapping spectra dispersed from different slits, with each slit being exposed to radiation from a diffe… ▽ More The signal measured by an astronomical spectrometer may be due to radiation from a multi-component mixture of plasmas with a range of physical properties (e.g. temperature, Doppler velocity). Confusion between multiple components may be exacerbated if the spectrometer sensor is illuminated by overlapping spectra dispersed from different slits, with each slit being exposed to radiation from a different portion of an extended astrophysical object. We use a compressed sensing method to robustly retrieve the different components. This method can be adopted for a variety of spectrometer configurations, including single-slit, multi-slit (e.g., the proposed MUlti-slit Solar Explorer mission; MUSE) and slot spectrometers (which produce overlappograms). △ Less

Submitted 11 February, 2019; originally announced February 2019.

Comments: submitted to ApJ

The Origin of Major Solar Activity - Collisional Shearing Between Nonconjugated Polarities of Multiple Bipoles Emerging Within Active Regions

Authors: Georgios Chintzoglou, Jie Zhang, Mark C. M. Cheung, Maria Kazachenko

Abstract: Active Regions (ARs) that exhibit compact Polarity Inversion Lines (PILs) are known to be very flare-productive. However, the physical mechanisms behind this statistical inference have not been demonstrated conclusively. We show that such PILs can occur due to the collision between two emerging flux tubes nested within the same AR. In such multipolar ARs, the flux tubes may emerge simultaneously o… ▽ More Active Regions (ARs) that exhibit compact Polarity Inversion Lines (PILs) are known to be very flare-productive. However, the physical mechanisms behind this statistical inference have not been demonstrated conclusively. We show that such PILs can occur due to the collision between two emerging flux tubes nested within the same AR. In such multipolar ARs, the flux tubes may emerge simultaneously or sequentially, each initially producing a bipolar magnetic region (BMR) at the surface. During each flux tube's emergence phase, the magnetic polarities can migrate such that opposite polarities belonging to different BMRs collide, resulting in shearing and cancellation of magnetic flux. We name this process 'collisional shearing' to emphasize that the shearing and flux cancellation develops due to the collision. Collisional shearing is a process different from the known concept of flux cancellation occurring between polarities of a single bipole, a process that has been commonly used in many numerical models. High spatial and temporal resolution observations from the Solar Dynamics Observatory for two emerging ARs, AR11158 and AR12017, show the continuous cancellation of up to 40% of the unsigned magnetic flux of the smallest BMR, which occurs at the collisional PIL for as long as the collision persists. The flux cancellation is accompanied by a succession of solar flares and CMEs, products of magnetic reconnection along the collisional PIL. Our results suggest that the quantification of magnetic cancellation driven by collisional shearing needs to be taken into consideration in order to improve the prediction of solar energetic events and space weather. △ Less

Submitted 9 January, 2019; v1 submitted 6 November, 2018; originally announced November 2018.

Comments: 68 pages, 15 figures, and 7 online movies. Accepted for publication in the Astrophysical Journal

A Truly Global Extreme Ultraviolet Wave from the SOL2017-09-10 X8.2+ Solar Flare-Coronal Mass Ejection

Authors: Wei Liu, Meng Jin, Cooper Downs, Leon Ofman, Mark C. M. Cheung, Nariaki V. Nitta

Abstract: We report SDO/AIA observations of an extraordinary global extreme ultraviolet (EUV) wave triggered by the X8.2+ flare-CME eruption on 2017 September 10. This was one of the best EUV waves ever observed with modern instruments, yet it was likely the last one of such magnitudes of Solar Cycle 24 as the Sun heads toward the minimum. Its remarkable characteristics include the following. (1) The wave w… ▽ More We report SDO/AIA observations of an extraordinary global extreme ultraviolet (EUV) wave triggered by the X8.2+ flare-CME eruption on 2017 September 10. This was one of the best EUV waves ever observed with modern instruments, yet it was likely the last one of such magnitudes of Solar Cycle 24 as the Sun heads toward the minimum. Its remarkable characteristics include the following. (1) The wave was observed, for the first time, to traverse the full-Sun corona over the entire visible solar disk and off-limb circumference, manifesting a truly global nature, owing to its exceptionally large amplitude, e.g., with EUV enhancements by up to 300% at 1.1 Rsun from the eruption. (2) This leads to strong transmissions (in addition to commonly observed reflections) in and out of both polar coronal holes, which are usually devoid of EUV waves. It has elevated wave speeds >2000 km/s within them, consistent with the expected higher fast-mode magnetosonic wave speeds. The coronal holes essentially serve as new "radiation centers" for the waves being refracted out of them, which then travel toward the equator and collide head-on, causing additional EUV enhancements. (3) The wave produces significant compressional heating to local plasma upon its impact, indicated by long-lasting EUV intensity changes and differential emission measure increases at higher temperatures (e.g., log T=6.2) accompanied by decreases at lower temperatures (e.g., log T=6.0). These characteristics signify the potential of such EUV waves for novel magnetic and thermal diagnostics of the solar corona {\it on global scales}. △ Less

Submitted 13 August, 2018; v1 submitted 25 July, 2018; originally announced July 2018.

Comments: Accepted by ApJ Letters (as of July 24, 2018), 8 pages, 5 figures; Update: minor edits made on August 14, 2018, after submitting ApJL proof corrections

Solar ultraviolet bursts

Authors: Peter R. Young, Hui Tian, Hardi Peter, Robert J. Rutten, Chris J. Nelson, Zhenghua Huang, Brigitte Schmieder, Gregal J. M. Vissers, Shin Toriumi, Luc H. M. Rouppe van der Voort, Maria S. Madjarska, Sanja Danilovic, Arkadiusz Berlicki, L. P. Chitta, Mark C. M. Cheung, Chad Madsen, Kevin P. Reardon, Yukio Katsukawa, Petr Heinzel

Abstract: The term "ultraviolet (UV) burst" is introduced to describe small, intense, transient brightenings in ultraviolet images of solar active regions. We inventorize their properties and provide a definition based on image sequences in transition-region lines. Coronal signatures are rare, and most bursts are associated with small-scale, canceling opposite-polarity fields in the photosphere that occur i… ▽ More The term "ultraviolet (UV) burst" is introduced to describe small, intense, transient brightenings in ultraviolet images of solar active regions. We inventorize their properties and provide a definition based on image sequences in transition-region lines. Coronal signatures are rare, and most bursts are associated with small-scale, canceling opposite-polarity fields in the photosphere that occur in emerging flux regions, moving magnetic features in sunspot moats, and sunspot light bridges. We also compare UV bursts with similar transition-region phenomena found previously in solar ultraviolet spectrometry and with similar phenomena at optical wavelengths, in particular Ellerman bombs. Akin to the latter, UV bursts are probably small-scale magnetic reconnection events occurring in the low atmosphere, at photospheric and/or chromospheric heights. Their intense emission in lines with optically thin formation gives unique diagnostic opportunities for studying the physics of magnetic reconnection in the low solar atmosphere. This paper is a review report from an International Space Science Institute team that met in 2016-2017. △ Less

Submitted 3 October, 2018; v1 submitted 15 May, 2018; originally announced May 2018.

Comments: Review article accepted for publication in Space Science Reviews

An Observationally-Constrained Model of a Flux Rope that Formed in the Solar Corona

Authors: Alexander W. James, Gherardo Valori, Lucie M. Green, Yang Liu, Mark C. M. Cheung, Yang Guo, Lidia van Driel-Gesztelyi

Abstract: Coronal mass ejections (CMEs) are large-scale eruptions of plasma from the coronae of stars. Understanding the plasma processes involved in CME initiation has applications to space weather forecasting and laboratory plasma experiments. James et al. (Sol. Phys. 292, 71, 2017) used EUV observations to conclude that a magnetic flux rope formed in the solar corona above NOAA Active Region 11504 before… ▽ More Coronal mass ejections (CMEs) are large-scale eruptions of plasma from the coronae of stars. Understanding the plasma processes involved in CME initiation has applications to space weather forecasting and laboratory plasma experiments. James et al. (Sol. Phys. 292, 71, 2017) used EUV observations to conclude that a magnetic flux rope formed in the solar corona above NOAA Active Region 11504 before it erupted on 14 June 2012 (SOL2012-06-14). In this work, we use data from the Solar Dynamics Observatory to model the coronal magnetic field of the active region one hour prior to eruption using a nonlinear force-free field extrapolation, and find a flux rope reaching a maximum height of 150 Mm above the photosphere. Estimations of the average twist of the strongly asymmetric extrapolated flux rope are between 1.35 and 1.88 turns, depending on the choice of axis, although the erupting structure was not observed to kink. The decay index near the apex of the axis of the extrapolated flux rope is comparable to typical critical values required for the onset of the torus instability, so we suggest that the torus instability drove the eruption. △ Less

Submitted 22 February, 2018; originally announced February 2018.

Comments: 8 pages, 5 figures

Witnessing a Large-scale Slipping Magnetic Reconnection along a Dimming Channel during a Solar Flare

Authors: Ju Jing, Rui Liu, Mark C. M. Cheung, Jeongwoo Lee, Yan Xu, Chang Liu, Chunming Zhu, Haimin Wang

Abstract: We report the intriguing large-scale dynamic phenomena associated with the M6.5 flare~(SOL2015-06-22T18:23) in NOAA active region 12371, observed by RHESSI, Fermi, and the Atmospheric Image Assembly (AIA) and Magnetic Imager (HMI) on the Solar Dynamic Observatory (SDO). The most interesting feature of this event is a third ribbon (R3) arising in the decay phase, propagating along a dimming channel… ▽ More We report the intriguing large-scale dynamic phenomena associated with the M6.5 flare~(SOL2015-06-22T18:23) in NOAA active region 12371, observed by RHESSI, Fermi, and the Atmospheric Image Assembly (AIA) and Magnetic Imager (HMI) on the Solar Dynamic Observatory (SDO). The most interesting feature of this event is a third ribbon (R3) arising in the decay phase, propagating along a dimming channel (seen in EUV passbands) towards a neighboring sunspot. The propagation of R3 occurs in the presence of hard X-ray footpoint emission, and is broadly visible at temperatures from 0.6~MK to over 10~MK through the Differential Emission Measure (DEM) analysis. The coronal loops then undergo an apparent slipping motion following the same path of R3, after a $\sim$ 80~min delay. To understand the underlying physics, we investigate the magnetic configuration and the thermal structure of the flaring region. Our results are in favor of a slipping-type reconnection followed by the thermodynamic evolution of coronal loops. In comparison with those previously reported slipping reconnection events, this one proceeds across a particularly long distance ($\sim$60 Mm) over a long period of time ($\sim$50 min), and shows two clearly distinguished phases: the propagation of the footpoint brightening driven by nonthermal particle injection and the apparent slippage of loops governed by plasma heating and subsequent cooling. △ Less

Submitted 5 June, 2017; originally announced June 2017.

Various Local Heating Events in the Earliest Phase of Flux Emergence

Authors: Shin Toriumi, Yukio Katsukawa, Mark C. M. Cheung

Abstract: Emerging flux regions (EFRs) are known to exhibit various sporadic local heating events in the lower atmosphere. To investigate the characteristics of these events, especially to link the photospheric magnetic fields and atmospheric dynamics, we analyze Hinode, IRIS, and SDO data of a new EFR in NOAA AR 12401. Out of 151 bright points (BPs) identified in Hinode/SOT Ca images, 29 are overlapped by… ▽ More Emerging flux regions (EFRs) are known to exhibit various sporadic local heating events in the lower atmosphere. To investigate the characteristics of these events, especially to link the photospheric magnetic fields and atmospheric dynamics, we analyze Hinode, IRIS, and SDO data of a new EFR in NOAA AR 12401. Out of 151 bright points (BPs) identified in Hinode/SOT Ca images, 29 are overlapped by an SOT/SP scan. Seven BPs in the EFR center possess mixed-polarity magnetic backgrounds in the photosphere. Their IRIS UV spectra (e.g., Si IV 1402.8 A) are strongly enhanced and red- or blue-shifted with tails reaching +/- 150 km/s, which is highly suggestive of bi-directional jets, and each brightening lasts for 10 - 15 minutes leaving flare-like light curves. Most of this group show bald patches, the U-shaped photospheric magnetic loops. Another 10 BPs are found in unipolar regions at the EFR edges. They are generally weaker in UV intensities and exhibit systematic redshifts with Doppler speeds up to 40 km/s, which could exceed the local sound speed in the transition region. Both types of BPs show signs of strong temperature increase in the low chromosphere. These observational results support the physical picture that heating events in the EFR center are due to magnetic reconnection within cancelling undular fields like Ellerman bombs, while the peripheral heating events are due to shocks or strong compressions caused by fast downflows along the overlying arch filament system. △ Less

Submitted 5 January, 2017; originally announced January 2017.

Comments: 17 pages, 10 figures, accepted for publication in ApJ. Corresponding movie for Figure 1 can be found at http://www-space.eps.s.u-tokyo.ac.jp/~toriumi/toriumi2017b.mov

A Numerical Study of Long-Range Magnetic Impacts during Coronal Mass Ejections

Authors: M. Jin, C. J. Schrijver, M. C. M. Cheung, M. L. DeRosa, N. V. Nitta, A. M. Title

Abstract: With the global view and high-cadence observations from SDO/AIA and STEREO, many spatially separated solar eruptive events appear to be coupled. However, the mechanisms for "sympathetic" events are still largely unknown. In this study, we investigate the impact of an erupting flux rope on surrounding solar structures through large-scale magnetic coupling. We build a realistic environment of the so… ▽ More With the global view and high-cadence observations from SDO/AIA and STEREO, many spatially separated solar eruptive events appear to be coupled. However, the mechanisms for "sympathetic" events are still largely unknown. In this study, we investigate the impact of an erupting flux rope on surrounding solar structures through large-scale magnetic coupling. We build a realistic environment of the solar corona on 2011 February 15 using a global magnetohydrodynamics (MHD) model and initiate coronal mass ejections (CMEs) in active region (AR) 11158 by inserting Gibson-Low analytical flux ropes. We show that a CME's impact on the surrounding structures depends not only on the magnetic strength of these structures and their distance to the source region, but also on the interaction between the CME with the large-scale magnetic field. Within the CME expansion domain where the flux rope field directly interacts with the solar structures, expansion-induced reconnection often modifies the overlying field, thereby increasing the decay index. This effect may provide a primary coupling mechanism underlying the sympathetic eruptions. The magnitude of the impact is found to depend on the orientation of the erupting flux rope, with the largest impacts occurring when the flux rope is favorably oriented for reconnecting with the surrounding regions. Outside the CME expansion domain, the influence of the CME is mainly through field line compression or post-eruption relaxation. Based on our numerical experiments, we discuss a way to quantify the eruption impact, which could be useful for forecasting purposes. △ Less

Submitted 15 March, 2016; originally announced March 2016.

Comments: 38 pages, 13 figures

Journal ref: The Astrophysical Journal, Volume 820, Issue 1, article id. 16, 2016

Numerical Study on Emergence of Kinked Flux Tube for Understanding of Possible Origin of Delta-spot Regions

Authors: Shinsuke Takasao, Yuhong Fan, Mark. C. M. Cheung, Kazunari Shibata

Abstract: We carried out a magnetohydrodynamics simulation where a subsurface twisted kink-unstable flux tube emerges from the solar interior to the corona. Unlike the previous expectations based on the bodily emergence of a knotted tube, we found that the kinked tube can spontaneously form a complex quadrupole structure at the photosphere. Due to the development of the kink instability before the emergence… ▽ More We carried out a magnetohydrodynamics simulation where a subsurface twisted kink-unstable flux tube emerges from the solar interior to the corona. Unlike the previous expectations based on the bodily emergence of a knotted tube, we found that the kinked tube can spontaneously form a complex quadrupole structure at the photosphere. Due to the development of the kink instability before the emergence, the magnetic twist at the kinked apex of the tube is greatly reduced, although the other parts of the tube is still strongly twisted. This leads to the formation of a complex quadrupole structure: a pair of the coherent, strongly twisted spots and a narrow complex bipolar pair between it. The quadrupole is formed by the submergence of a portion of emerged magnetic fields. This result is relevant for understanding of the origin of the complex multipolar $δでるた$-spot regions that have a strong magnetic shear and emerge with polarity orientations not following Hale-Nicholson and Joy Laws. △ Less

Submitted 9 November, 2015; originally announced November 2015.

Comments: 34 pages, 17 figures, published in ApJ

Journal ref: ApJ, 813, 112, 2015 November 10

Light Bridge in a Developing Active Region. II. Numerical Simulation of Flux Emergence and Light Bridge Formation

Authors: Shin Toriumi, Mark C. M. Cheung, Yukio Katsukawa

Abstract: 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 mag… ▽ More 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. △ Less

Submitted 1 September, 2015; originally announced September 2015.

Comments: 39 pages, 11 figures, 2 tables, accepted for publication in ApJ

Light Bridge in a Developing Active Region. I. Observation of Light Bridge and its Dynamic Activity Phenomena

Authors: Shin Toriumi, Yukio Katsukawa, Mark C. M. Cheung

Abstract: 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… ▽ More 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. △ Less

Submitted 1 September, 2015; originally announced September 2015.

Comments: 38 pages, 9 figures, 1 table, accepted for publication in ApJ

Multi-Parametric Study of Rising 3D Buoyant Flux Tubes in an Adiabatic Stratification Using AMR

Authors: Juan Martinez-Sykora, Fernando Moreno-Insertis, Mark C. M. Cheung

Abstract: 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… ▽ More 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. △ Less

Submitted 6 July, 2015; originally announced July 2015.

Comments: 20 pages, 13 Figures, Accepted to be published in ApJ

The Coronal Global Evolutionary Model: Using HMI Vector Magnetogram and Doppler Data to Model the Buildup of Free Magnetic Energy in the Solar Corona

Authors: George H. Fisher, William. P. Abbett, David J. Bercik, Maria D. Kazachenko, Benjamin J. Lynch, Brian T. Welsch, J. Todd Hoeksema, Keiji Hayashi, Yang Liu, Aimee A. Norton, Alberto Sainz Dalda, Xudong Sun, Marc L. DeRosa, Mark C. M. Cheung

Abstract: 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 us… ▽ More 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えーゆー. △ Less

Submitted 9 July, 2015; v1 submitted 22 May, 2015; originally announced May 2015.

Comments: 7 pages, 2 figures. Feature article in Space Weather 13