Search for: All records

Abstract The spectra of coronal mass ejections (CMEs) in the low corona play a crucial role in understanding their origins and physical mechanisms and enhancing space weather forecasting. However, capturing these spectra faces significant challenges. This paper introduces a scheme of a multislit spectrometer design with five slits, acquiring the global spectra of the solar corona simultaneously with a focus on the spectra of CMEs in the low corona. The chosen wavelength range of the spectrometer (170–180 Å) includes four extreme ultraviolet emission lines (Fex174.53 Å, Feix171.07 Å, Fex175.26 Å, Fex177.24 Å), which provides information on the plasma velocity, density, and temperature. Utilizing a numerical simulation of the global corona for both the on-disk and the off-limb scenarios, we focus on resolving the ambiguity associated with various Doppler velocity components of CMEs, particularly for a fast CME in the low corona. A new application of our decomposition technique is adopted, enabling the successful identification of multiple discrete CME velocity components. Our findings demonstrate a strong correlation between the synthetic model spectra and the inverted results, indicating the robustness of our decomposition method and its significant potential for global monitoring of the solar corona, including CMEs. 

Abstract Global solar photospheric magnetic maps play a critical role in solar and heliospheric physics research. Routine magnetograph measurements of the field occur only along the Sun–Earth line, leaving the far side of the Sun unobserved. Surface flux transport (SFT) models attempt to mitigate this by modeling the surface evolution of the field. While such models have long been established in the community (with several releasing public full-Sun maps), none are open source. The Open-source Flux Transport (OFT) model seeks to fill this gap by providing an open and user-extensible SFT model that also builds on the knowledge of previous models with updated numerical and data acquisition/assimilation methods along with additional user-defined features. In this first of a series of papers on OFT, we introduce its computational core: the High-performance Flux Transport (HipFT) code (https://github.com/predsci/hipft). HipFT implements advection, diffusion, and data assimilation in a modular design that supports a variety of flow models and options. It can compute multiple realizations in a single run across model parameters to create ensembles of maps for uncertainty quantification and is high-performance through the use of multi-CPU and multi-GPU parallelism. HipFT is designed to enable users to write extensions easily, enhancing its flexibility and adaptability. We describe HipFT’s model features, validations of its numerical methods, performance of its parallel and GPU-accelerated code implementation, analysis/postprocessing options, and example use cases. 

Abstract We generalize a magnetogram-matching Biot–Savart law (BSl) from planar to spherical geometry. For a given coronal current densityJ, this law determines the magnetic field $\tilde{\mathit{B}}$whose radial component vanishes at the surface. The superposition of $\tilde{\mathit{B}}$with a potential field defined by a given surface radial field,B_r, provides the entire configuration whereB_rremains unchanged by the currents. Using this approach, we (1) upgrade our regularized BSls for constructing coronal magnetic flux ropes (MFRs) and (2) propose a new method for decomposing a measured photospheric magnetic field as $\mathit{B}={\mathit{B}}_{\mathrm{pot}}+{\mathit{B}}_T+{\mathit{B}}_{\tilde{S}}$, where the potential,B_pot, toroidal,B_T, and poloidal, ${\mathit{B}}_{\tilde{S}}$, fields are determined byB_r,J_r, and the surface divergence ofB–B_pot, respectively, all derived from magnetic data. OurB_Tis identical to the one in the alternative Gaussian decomposition by P. W. Schuck et al., whileB_potand ${\mathit{B}}_{\tilde{S}}$are different from their poloidal fields ${\mathit{B}}_{\mathrm{P}}^{<}$and ${\mathit{B}}_{\mathrm{P}}^{>}$, which arepotentialin the infinitesimal proximity to the upper and lower side of the surface, respectively. In contrast, our ${\mathit{B}}_{\tilde{S}}$has no such constraints and, asB_potandB_T, refers to thesameupper side of the surface. In spite of these differences, for a continuousJdistribution across the surface,B_potand ${\mathit{B}}_{\tilde{S}}$are linear combinations of ${\mathit{B}}_{\mathrm{P}}^{<}$and ${\mathit{B}}_{\mathrm{P}}^{>}$. We demonstrate that, similar to the Gaussian method, our decomposition allows one to identify the footprints and projected surface-location of MFRs in the solar corona, as well as the direction and connectivity of their currents. 

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. 

Abstract On 3 February 2022, SpaceX launched 49 Starlink satellites, 38 of which unexpectedly de‐orbited. Although this event was attributed to space weather, definitive causality remained elusive because space weather conditions were not extreme. In this study, we identify solar sources of the interplanetary coronal mass ejections that were responsible for the geomagnetic storms around the time of launch of the Starlink satellites and for the first time, investigate their impact on Earth's magnetosphere using magnetohydrodynamic modeling. The model results demonstrate that the satellites were launched into an already disturbed space environment that persisted over several days. However, on performing comparative satellite orbital decay analyses, we find that space weather alone was not responsible but conspired together with a low‐altitude insertion and low satellite mass‐to‐area ratio to precipitate this unusual loss. Our work bridges space weather causality across the Sun–Earth system—with relevance for space‐based human technologies. 

Abstract Full-disk spectroscopic observations of the solar corona are highly desired to forecast solar eruptions and their impact on planets and to uncover the origin of solar wind. In this paper, we introduce a new multislit design (five slits) to obtain extreme-ultraviolet (EUV) spectra simultaneously. The selected spectrometer wavelength range (184–197 Å) contains several bright EUV lines that can be used for spectral diagnostics. The multislit approach offers an unprecedented way to efficiently obtain the global spectral data but the ambiguity from different slits should be resolved. Using a numerical simulation of the global corona, we primarily concentrate on the optimization of the disambiguation process, with the objective of extracting decomposed spectral information of six primary lines. This subsequently facilitates a comprehensive series of plasma diagnostics, including density (Fexii195.12/186.89 Å), Doppler velocity (Fexii193.51 Å), line width (Fexii193.51 Å), and temperature diagnostics (Feviii185.21 Å, Fex184.54 Å, Fexi188.22 Å, and Fexii193.51 Å). We find a good agreement between the forward modeling parameters and the inverted results at the initial eruption stage of a coronal mass ejection, indicating the robustness of the decomposition method and its immense potential for global monitoring of the solar corona. 

Abstract Coronal Mass Ejections (CMEs) are immense eruptions of plasma and magnetic fields that are propelled outward from the Sun, sometimes with velocities greater than 2000 km/s. They are responsible for some of the most severe space weather at Earth, including geomagnetic storms and solar energetic particle (SEP) events. We have developed CORHEL-CME, an interactive tool that allows non-expert users to routinely model multiple CMEs in a realistic coronal and heliospheric environment. The tool features a web-based user interface that allows the user to select a time period of interest, and employs Regularized Biot-Savart Law (RBSL) flux ropes to create stable and unstable pre-eruptive configurations within a background global magnetic field. The properties of these configurations can first be explored in a zero-beta magnetohydrodynamic (MHD) model, followed by complete CME simulations in thermodynamic MHD, with propagation out to 1 AU. We describe design features of the interface and computations, including the innovations required to efficiently compute results on practical timescales with moderate computational resources. CORHEL-CME is now implemented at NASA's Community Coordinated Modeling Center (CCMC) using NASA Amazon Web Services (AWS). It will be available to the public in early 2024. 

Abstract Unconditionally stable time stepping schemes are useful and often practically necessary for advancing parabolic operators in multi-scale systems. However, serious accuracy problems may emerge when taking time steps that far exceed the explicit stability limits. In our previous work, we compared the accuracy and performance of advancing parabolic operators in a thermodynamic MHD model using an implicit method and an explicit super time-stepping (STS) method. We found that while the STS method outperformed the implicit one with overall good results, it was not able to damp oscillatory behavior in the solution efficiently, hindering its practical use. In this follow-up work, we evaluate an easy-to-implement method for selecting a practical time step limit (PTL) for unconditionally stable schemes. This time step is used to ‘cycle’ the operator-split thermal conduction and viscosity parabolic operators. We test the new time step with both an implicit and STS scheme for accuracy, performance, and scaling. We find that, for our test cases here, the PTL dramatically improves the STS solution, matching or improving the solution of the original implicit scheme, while retaining most of its performance and scaling advantages. The PTL shows promise to allow more accurate use of unconditionally stable schemes for parabolic operators and reliable use of STS methods. 

Abstract Observations have shown a clear association of filament/prominence eruptions with the emergence of magnetic flux in or near filament channels. Magnetohydrodynamic (MHD) simulations have been employed to systematically study the conditions under which such eruptions occur. These simulations to date have modeled filament channels as 2D flux ropes or 3D uniformly sheared arcades. Here we present MHD simulations of flux emergence into a more realistic configuration consisting of a bipolar active region containing a line-tied 3D flux rope. We use the coronal flux-rope model of Titov et al. as the initial condition and drive our simulations by imposing boundary conditions extracted from a flux emergence simulation by Leake et al. We identify three mechanisms that determine the evolution of the system: (i) reconnection displacing footpoints of field lines overlying the coronal flux rope, (ii) changes of the ambient field due to the intrusion of new flux at the boundary, and (iii) interaction of the (axial) electric currents in the preexisting and newly emerging flux systems. The relative contributions and effects of these mechanisms depend on the properties of the preexisting and emerging flux systems. Here we focus on the location and orientation of the emerging flux relative to the coronal flux rope. Varying these parameters, we investigate under which conditions an eruption of the latter is triggered. 

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.