An official website of the United States government
Abstract We describe, test, and apply a technique to incorporate full-Sun, surface flux evolution into an MHD model of the global solar corona. Requiring only maps of the evolving surface flux, our method is similar to that of Lionello et al., but we introduce two ways to correct the electric field at the lower boundary to mitigate spurious currents. We verify the accuracy of our procedures by comparing to a reference simulation, driven with known flows and electric fields. We then present a thermodynamic MHD calculation lasting one solar rotation driven by maps from the magnetic flux evolution model of Schrijver & DeRosa. The dynamic, time-dependent nature of the model corona is illustrated by examining the evolution of the open flux boundaries and forward-modeled EUV emission, which evolve in response to surface flows and the emergence and cancellation flux. Although our main goal is to present the method, we briefly investigate the relevance of this evolution to properties of the slow solar wind, examining the mapping of dipped field lines to the topological signatures of the “S-Web” and comparing charge state ratios computed in the time-dependently driven run to a steady-state equivalent. Interestingly, we find that driving on its own does not significantly improve the charge state ratios, at least in this modest resolution run that injects minimal helicity. Still, many aspects of the time-dependently driven model cannot be captured with traditional steady-state methods, and such a technique may be particularly relevant for the next generation of solar wind and coronal mass ejection models.
more »
« less
Time-dependent Dynamics of the Corona
Abstract We present in this Letter the first global comparison between traditional line-tied steady-state magnetohydrodynamic models and a new, fully time-dependent thermodynamic magnetohydrodynamic simulation of the global corona. To approximate surface magnetic field distributions and magnitudes around solar minimum, we use the Lockheed Evolving Surface-Flux Assimilation Model to obtain input maps that incorporate flux emergence and surface flows over a full solar rotation, including differential rotation and meridional flows. Each time step evolves the previous state of the plasma with a new magnetic field input boundary condition, mimicking photospheric driving on the Sun. We find that this method produces a qualitatively different corona compared to steady-state models. The magnetic energy levels are higher in the time-dependent model, and coronal holes evolve more along the following edge than they do in steady-state models. Coronal changes, as illustrated with forward-modeled emission maps, evolve on longer timescales with time-dependent driving. We discuss implications for active and quiet Sun scenarios, solar wind formation, and widely used steady-state assumptions like potential field source surface calculations.
more »
« less
- Award ID(s):
- 1854790
- PAR ID:
- 10478597
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal Letters
- Volume:
- 959
- Issue:
- 1
- ISSN:
- 2041-8205
- Format(s):
- Medium: X Size: Article No. L4
- Size(s):
- Article No. L4
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract One systematic limitation of solar coronal hole (CH) detection at extreme ultraviolet (EUV) wavelengths is the obscuration of dark regions of the corona by brighter structures along the line of sight. Another problem arises when using CHs to compute the Sun’s open magnetic flux, where surface measurements of the radial magnetic field, , are situated slightly below the effective height of coronal EUV emission. In this paper, we explore these two limitations utilizing a thermodynamic magnetohydrodynamic (MHD) model of the corona for Carrington rotation (CR) 2101, where we generate CH detections from EUV 193 Å images of the corona forward-modeled from the MHD solution, and where the modeled open field is known. We demonstrate a method to combine EUV images into a full Sun map that helps alleviate CH obscuration called theminimum intensity diskmerge(MIDM). We also show the variation in measured open flux and CH area that is due to the effective height differences between EUV and measurements. We then apply the MIDM method to SDO/AIA 193 Å observations from CR 2101, and conduct an analogous analysis. In this case, the MIDM method uses time-varying images, the effects of which are discussed. We show that overall, the MIDM method and an appreciation of the effective height mismatch provide a useful new way to extract a broader view of CHs, especially near the poles. In turn, they enable improved estimates of the open magnetic flux, and help facilitate comparisons between models and observations.more » « less
-
The NASA-NSF sponsored Space Weather with Quantified Uncertainty (SWQU) project's main objective is to develop a data-driven, time-dependent, open source model of the solar corona and heliosphere. One key component of the SWQU effort is using a data-assimilation flux transport model to generate an ensemble of synchronic radial magnetic field maps as boundary conditions for the coronal field model. To accomplish this goal, we are developing a new Open-source Flux Transport (OFT) software suite. While there are a number of established flux transport models in the community, OFT is distinguished from many of these efforts in 3 key attributes: (1) It is based on modern computing techniques that will allow many realizations to be rapidly computed on multi-core systems and/or GPUs, (2) it is designed to be easily extensible, and (3) OFT will be released as an open source project. OFT consists of three software packages: 1) OFTpy: a python package for data acquisition, database organization, and Carrington map processing, 2) ConFlow: a Fortran code that generates super granular convective flows, and 3) High-Performance Flux Transport (HipFT): a modular, GPU-accelerated Fortran code for modeling surface flux transport with data assimilation. Here, we present the current state of the OFT project, key features and methods of OFTpy, ConFlow, and HipFt, and real-world examples of data-assimilation and flux transport with HipFT. Validation and performance tests are shown, including generating an ensemble of OFT maps.more » « less
-
Context.Coronal mass ejections (CMEs) are eruptions of plasma from the Sun that travel through interplanetary space and may encounter Earth. CMEs often enclose a magnetic flux rope (MFR), the orientation of which largely determines the CMEs’ geoeffectiveness. Current operational CME models do not model MFRs, but a number of research ones do, including the Open Solar Physics Rapid Ensemble Information (OSPREI) model. Aims.We report the sensitivity of OSPREI to a range of user-selected photospheric and coronal conditions. Methods.We modeled four separate CMEs observed in situ by Parker Solar Probe (PSP). We varied the input photospheric conditions using four input magnetograms (HMI Synchronic, HMI Synoptic, GONG Synoptic Zero-Point Corrected, and GONG ADAPT). To vary the coronal field reconstruction, we employed the Potential Field Source Surface (PFSS) model and varied its source-surface height in the range 1.5–3.0R⊙with 0.1R⊙increments. Results.We find that both the input magnetogram and PFSS source surface often affect the evolution of the CME as it propagates through the Sun’s corona into interplanetary space, and therefore the accuracy of the MFR prediction compared to in situ data at PSP. There is no obvious best combination of input magnetogram and PFSS source surface height. Conclusions.The OSPREI model is moderately sensitive to the input photospheric and coronal conditions. Based on where the source region of the CME is located on the Sun, there may be best practices when selecting an input magnetogram to use.more » « less
-
Abstract A plethora of coronal models, from empirical to more complex magnetohydrodynamic (MHD) ones, are being used for reconstructing the coronal magnetic field topology and estimating the open magnetic flux. However, no individual solution fully agrees with coronal hole observations and in situ measurements of open flux at 1 au, as there is a strong deficit between the model and observations contributing to the known problem of the missing open flux. In this paper, we investigate the possible origin of the discrepancy between modeled and observed magnetic field topology by assessing the effect on the simulation output by the choice of the input boundary conditions and the simulation setup, including the choice of numerical schemes and the parameter initialization. In the frame of this work, we considered four potential field source surface-based models and one fully MHD model, different types of global magnetic field maps, and model initiation parameters. After assessing the model outputs using a variety of metrics, we conclude that they are highly comparable regardless of the differences set at initiation. When comparing all models to coronal hole boundaries extracted by extreme-ultraviolet filtergrams, we find that they do not compare well. This mismatch between observed and modeled regions of the open field is a candidate contributing to the open flux problem.more » « less
