skip to main content


Title: MSWIM2D: Two-dimensional Outer Heliosphere Solar Wind Modeling
Abstract

The vast size of the Sun’s heliosphere, combined with sparse spacecraft measurements over that large domain, makes numerical modeling a critical tool to predict solar wind conditions where there are no measurements. This study models the solar wind propagation in 2D using the BATSRUS MHD solver to form the MSWIM2D data set of solar wind in the outer heliosphere. Representing the solar wind from 1 to 75 au in the ecliptic plane, a continuous model run from 1995–present has been performed. The results are available for free athttp://csem.engin.umich.edu/mswim2d/. The web interface extracts output at desired locations and times. In addition to solar wind ions, the model includes neutrals coming from the interstellar medium to reproduce the slowing of the solar wind in the outer heliosphere and to extend the utility of the model to larger radial distances. The inclusion of neutral hydrogen is critical to recreating the solar wind accurately outside of ∼4 au. The inner boundary is filled by interpolating and time-shifting in situ observations from L1 and STEREO spacecraft when available. Using multiple spacecraft provides a more accurate boundary condition than a single spacecraft with time shifting alone. Validations of MSWIM2D are performed using MAVEN and New Horizons observations. The results demonstrate the efficacy of this model to propagate the solar wind to large distances and obtain practical, useful solar wind predictions. For example, the rms error of solar wind speed prediction at Mars is only 66 km s−1and at Pluto is a mere 25 km s−1.

 
more » « less
Award ID(s):
1663800 2031019
NSF-PAR ID:
10368011
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal Supplement Series
Volume:
260
Issue:
2
ISSN:
0067-0049
Format(s):
Medium: X Size: Article No. 43
Size(s):
["Article No. 43"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Small-scale interplanetary magnetic flux ropes (SMFRs) are similar to ICMEs in magnetic structure, but are smaller and do not exhibit coronal mass ejection plasma signatures. We present a computationally efficient and GPU-powered version of the single-spacecraft automated SMFR detection algorithm based on the Grad–Shafranov (GS) technique. Our algorithm can process higher resolution data, eliminates selection bias caused by a fixed 〈B〉 threshold, has improved detection criteria demonstrated to have better results on an MHD simulation, and recovers full 2.5D cross sections using GS reconstruction. We used it to detect 512,152 SMFRs from 27 yr (1996–2022) of 3 s cadence Wind measurements. Our novel findings are the following: (1) the SMFR filling factor (∼ 35%) is independent of solar activity, distance to the heliospheric current sheet, and solar wind plasma type, although the minority of SMFRs with diameters greater than ∼0.01 au have a strong solar activity dependence; (2) SMFR diameters follow a log-normal distribution that peaks below the resolved range (≳104km), although the filling factor is dominated by SMFRs between 105and 106km; (3) most SMFRs at 1 au have strong field-aligned flows like those from Parker Solar Probe measurements; (4) the radial density (generally ∼1 detected per 106km) and axial magnetic flux density of SMFRs are higher in faster solar wind types, suggesting that they are more compressed. Implications for the origin of SMFRs and switchbacks are briefly discussed. The new algorithm and SMFR dataset are made freely available.

     
    more » « less
  2. Abstract

    Interstellar neutral atoms enter the heliosphere at a relatively slow speed corresponding to the motion of the Sun through the local interstellar medium, which is approximately 25 km s−1. Neutral hydrogen atoms enter from the approximate location of the Voyager spacecraft and are eventually ionized primarily by collision with thermal solar wind ions. An earlier analysis by Hollick et al. examined low-frequency magnetic waves observed by the Voyager spacecraft from launch through 1990 that are thought to arise from the scattering of newborn interstellar pickup H+and He+. We report an analysis of Voyager 1 observations in 1991, which is the last year of high-resolution magnetic field data that are publicly available, and find 70 examples of low-frequency waves with the characteristics that suggest excitation by pickup H+and 10 examples of waves consistent with excitation by pickup He+. We find a particularly dense cluster of observations at the tail end of what is thought to be a Merged Interaction Region (MIR) that was previously studied by Burlaga & Ness using Voyager 2 observations. This is not unexpected if the MIR is followed by a large rarefaction region, as they tend to be regions of reduced turbulence levels that permit the growth of the waves over the long time periods that are generally required of this instability.

     
    more » « less
  3. Abstract

    We present the results of processing the effects of the powerful gamma-ray burst GRB221009A captured by the charged particle detectors (electrostatic analyzers and solid-state detectors) on board spacecraft at different points in the heliosphere on 2022 October 9. To follow the GRB221009A propagation through the heliosphere, we used the electron and proton flux measurements from solar missions Solar Orbiter and STEREO-A; Earth’s magnetosphere and solar wind missions THEMIS and Wind; meteorological satellites POES15, POES19, and MetOp3; and MAVEN—a NASA mission orbiting Mars. GRB221009A had a structure of four bursts: the less intense Pulse 1—the triggering impulse—was detected by gamma-ray observatories atT0= 13:16:59 UT (near the Earth); the most intense Pulses 2 and 3 were detected on board all the spacecraft from the list; and Pulse 4 was detected in more than 500 s after Pulse 1. Due to their different scientific objectives, the spacecraft, whose data were used in this study, were separated by more than 1 au (Solar Orbiter and MAVEN). This enabled the tracking of GRB221009A as it was propagating across the heliosphere. STEREO-A was the first to register Pulse 2 and 3 of the GRB, almost 100 s before their detection by spacecraft in the vicinity of Earth. MAVEN detected GRB221009A Pulses 2, 3, and 4 at the orbit of Mars about 237 s after their detection near Earth. By processing the observed time delays, we show that the source location of the GRB221009A was at R.A. 288.°5, decl. 18.°5 ± 2° (J2000).

     
    more » « less
  4. Abstract Interstellar pickup ions are an ubiquitous and thermodynamically important component of the solar wind plasma in the heliosphere. These PUIs are born from the ionization of the interstellar neutral gas, consisting of hydrogen, helium, and trace amounts of heavier elements, in the solar wind as the heliosphere moves through the local interstellar medium. As cold interstellar neutral atoms become ionized, they form an energetic ring beam distribution comoving with the solar wind. Subsequent scattering in pitch angle by intrinsic and self-generated turbulence and their advection with the radially expanding solar wind leads to the formation of a filled-shell PUI distribution, whose density and pressure relative to the thermal solar wind ions grows with distance from the Sun. This paper reviews the history of in situ measurements of interstellar PUIs in the heliosphere. Starting with the first detection in the 1980s, interstellar PUIs were identified by their highly nonthermal distribution with a cutoff at twice the solar wind speed. Measurements of the PUI distribution shell cutoff and the He focusing cone, a downwind region of increased density formed by the solar gravity, have helped characterize the properties of the interstellar gas from near-Earth vantage points. The preferential heating of interstellar PUIs compared to the core solar wind has become evident in the existence of suprathermal PUI tails, the nonadiabatic cooling index of the PUI distribution, and PUIs’ mediation of interplanetary shocks. Unlike the Voyager and Pioneer spacecraft, New Horizon’s Solar Wind Around Pluto (SWAP) instrument is taking the only direct measurements of interstellar PUIs in the outer heliosphere, currently out to $\sim47~\text{au}$ ∼ 47 au from the Sun or halfway to the heliospheric termination shock. 
    more » « less
  5. ABSTRACT

    In this work, we present numerical simulations of Stream Interaction Regions (SIRs) and Corotating Interaction Regions (CIRs) using the sunrunner3d tool that employs as a coronal model the boundary conditions obtained by corhel/mas with the pluto code that describes the global 3D structure of the solar wind using the magnetohydrodynamics (MHD) approach in the inner heliosphere. Specifically, we selected a set of SIRs and CIRs observed by the Parker Solar Probe (PSP) and STEREO-A (STA) missions during the Carrington rotations (CRs) 2207 to 2210 and CRs from 2020 to 2022. In order to describe the dynamics of the plasma that constitutes the solar wind background conditions for the selected CRs, we solve the ideal MHD equations in an inertial frame of reference, managing the solar rotation by rotating the boundary values in ϕ (longitude) at a rate corresponding to the sidereal rotation rate of the solar equator. We show that our results using sunrunner3d can globally reproduce the plasma parameters, such as radial velocity, number proton density, and radial magnetic field strength of these large-scale structures, observed by PSP and STA at distances near the Sun and around 1 au, respectively. These results allow exploring the global evolution of SIRs/CIRs in the inner heliosphere using sunrunner3d.

     
    more » « less