Solar Wind Density Research Articles

Constraints to the drag-based reverse modeling

One of the most widely used space weather forecast models to simulate the propagation of coronal mass ejections (CMEs) is the analytical drag-based model (DBM). It predicts the CME arrival time and speed at Earth or at a specific target (planets, spacecraft) in the Solar System. The corresponding drag-based ensemble model (DBEM) additionally takes into account the uncertainty of the input parameters by making n ensembles and provides the most probable arrival time and speed as well as their uncertainty intervals. An important input parameter for DBM and DBEM is the drag parameter γ, which depends on the CME cross-section and mass, as well as the solar-wind density. The reverse-modeling technique applied to the DBM allows us to derive γ values that minimize transit time (TT) and arrival-speed (v_tar) errors. The present study highlights the limitations and constraints of such a procedure. We searched for optimal γ values that would yield the perfect TT within one hour of the actual observed CME transit time as well as perfect v_tar within ±75 km,s^-1. This optimal window for v_tar was found by increasing v_tar from ±10 to ±100 km,s^-1, where the ±75 km,s^-1 window gave the perfect TT and v_tar in the case of 87% of CMEs compared to the ±10 km,s^-1 window, which was used in some previous reverse-modeling studies and gave optimal results for only 45% of the events from our CME list. For our analysis, a 31 CME-ICME pair sample is used from the period spanning 1997 to 2018. The reverse-modeling method is applied using the DBEMv3 tool for different γ ranges from 0.01 to 10 times km^-1. We tested whether and how the obtained optimal γ depends on the chosen γ range. By increasing the γ range, we find that the optimal γ converges to a certain value for two thirds of the analyzed events. The highly constrained γ ranges resulted in shifted and skewed γ distributions. By using the largest γ range (0.01--10 times 10^ km^-1), the medians of the optimal γ distributions are obtained for two thirds of the events in the common operational DBEMv3 range of 0.01--0.5 times 10^ km^-1. We also found that the important quantity in determining the range of γ distribution and ability to find an optimal γ is the difference between the CME launch speed and the solar-wind speed (v_0-w), which together with γ define the drag acceleration in the DBM. For small v_0-w differences (e.g., < 200 km,s^-1), the reverse modeling may not be the appropriate method to find the optimal γ due to large divergence of γ values found, which may additionally be caused by larger input uncertainties and physical model limitations in turn leading to inappropriate γ values.

Mars Nightside Ionospheric Response During the Disappearing Solar Wind Event: First Results

AbstractWe investigated, for the first time, the impact of the disappearing solar wind (DSW) event [26–28 December 2022] on the deep nightside ionospheric species using MAVEN data sets. An enhanced plasma density has been observed in the Martian nightside ionosphere during extreme low solar wind density and pressure periods. At a given altitude, the electron density surged by ∼2.5 times, while for ions (NO+, O2+, CO2+, C+, N+, O+, and OH+), it enhanced by > 10 times, respectively, compared to their typical average quiet‐time periods. This investigation suggests that an upward ionospheric expansion likely took place in a direct consequence to the contrasting low dynamic/magnetic pressure and relatively higher nightside ionospheric pressure (by 1–2 orders) causing an increased ionospheric density. Moreover, the day‐to‐night plasma transport may also be a contributing factor to the increased plasma density. Thus, this study offers a new insight about planetary atmosphere/ionosphere during extreme quiescent solar wind periods.

The impact on astrometry by solar-wind effect in pulsar timing

Abstract Astrometry of pulsars, particularly their distances, serves as a critical input for various astrophysical experiments using pulsars. Pulsar timing is a primary approach for determining a pulsar’s position, parallax, and distance. In this paper, we explore the influence of the solar wind on astrometric measurements obtained through pulsar timing, focusing on its potential to affect the accuracy of these parameters. Using both theoretical calculation and mock-data simulations, we demonstrate a significant correlation between the pulsar position, annual parallax and the solar-wind density parameters. This correlation strongly depends on the pulsar’s ecliptic latitude. We show that fixing solar-wind density to an arbitrary value in the timing analysis can introduce significant bias in the estimated pulsar position and parallax, and its significance is highly dependent on the ecliptic latitude of the pulsar and the timing precision of the data. For pulsars with favourable ecliptic latitude and timing precision, the astrometric and solar-wind parameters can be measured jointly with other timing parameters using single-frequency data. The parameter correlation can be mitigated by using multi-frequency data, which also significantly improves the measurement precision of these parameters; this is particularly important for pulsars at a medium or high ecliptic latitude. Additionally, for a selection of pulsars we reprocess their EPTA Data Release 2 data to include modelling of solar-wind effect in the timing analysis. This delivers significant measurements of both parallax and solar-wind density, the latter of which are consistent with those obtained at low-frequency band. In the future, combining pulsar timing data at gigahertz and lower frequencies will probably deliver the most robust and precise measurements of astrometry and solar wind properties in pulsar timing.

Global Simulation of the Solar Wind: A Comparison with Parker Solar Probe Observations during 2018–2022

Global magnetohydrodynamic (MHD) models play an important role in the infrastructure of space weather forecasting. Validating such models commonly utilizes in situ solar wind measurements made near the Earth’s orbit. The purpose of this study is to test the performance of G3DMHD (a data driven, time-dependent, 3D MHD model of the solar wind) with Parker Solar Probe (PSP) measurements. Since its launch in 2018 August, PSP has traversed the inner heliosphere at different radial distances sunward of the Earth (the closest approach ∼13.3 R ⊙), thus providing a good opportunity to study evolution of the solar wind and to validate heliospheric models of the solar wind. The G3DMHD model simulation is driven by a sequence of maps of the photospheric field extrapolated to the assumed source surface (2.5 R ⊙) using the potential field model from 2018 to 2022, which covers the first 15 PSP orbits. The Pearson correlation coefficient (cc) and the mean absolute scaled error (MASE) are used as the metrics to evaluate the model performance. It is found that the model performs better for both magnetic intensity (cc = 0.75; MASE = 0.60) and the solar wind density (cc = 0.73; MASE = 0.50) than for the solar wind speed (cc = 0.15; MASE = 1.29) and temperature (cc = 0.28; MASE = 1.14). This is due primarily to lack of accurate boundary conditions. The well-known underestimate of the magnetic field in solar minimum years is also present. Assuming that the radial magnetic field becomes uniformly distributed with latitude at or below 18 R ⊙ (the inner boundary of the computation domain), the agreement in the magnetic intensity significantly improves (cc = 0.83; MASE = 0.49).

Heavy Ion Escape at Mars during the Disappearing Solar Wind Event in 2022 December

A unique event known as the disappearing solar wind (DSW), characterized by an extremely low-density solar wind stream, occurred at Mars on 2022 December 26–27. As this stream flowed past Mars, several properties of the Mars–solar wind interaction changed in response to the density drop. We utilize in situ plasma measurements from the Mars Atmosphere and Volatile EvolutioN spacecraft to examine heavy ion escape during this event. We find that escaping ions with energies above and below 30 eV responded differently to the reduction in solar wind density. High-energy ions experienced a decrease in flux, whereas low-energy ions experienced an increase, regardless of whether they were escaping through the plume or tailward channel. Furthermore, we observe a net reduction in the flux of plume escaping ions during the DSW period, primarily due to a considerable reduction in the high-energy component, which typically dominates plume escape. In contrast, the overall flux of tailward escaping ions increased during this period. These variations in heavy ion escape are mainly attributed to substantial reductions in solar wind dynamic pressure and momentum density. This paper provides new insights into the dynamics of heavy ion escape under an exceptional state of the Mars–solar wind interaction.

Magnetosonic waves in the Martian ionosphere driven by upstream proton cyclotron waves: Two-point observations by MAVEN and Mars Express

Recent observations from the Mars Atmosphere and Volatile EvolutioN (MAVEN) and Mars Express (MEX) spacecraft have suggested that pressure pulses originating from upstream proton cyclotron waves (PCWs) can “ring” the Martian magnetopause at the same frequency and drive magnetosonic waves in the upper ionosphere of Mars, thereby transporting energy from the solar wind into the ionosphere. However, the limitation of single-spacecraft measurements prevents simultaneous observations of the driver and response of this “ringing” process of the Martian magnetosphere. Here we utilize two-point measurements from MAVEN and MEX to characterize the ringing probability at which upstream PCWs drive compressional fluctuations in the ionospheric magnetic field. We develop an algorithm to identify PCW-driven magnetosonic waves in the upper ionosphere of Mars from the two-point magnetic field data. The derived ringing probability is higher on the dayside, outside strong crustal magnetic fields, and under high solar wind density conditions. We also show that the median power of dayside ionospheric magnetic field fluctuations is enhanced by a factor of ∼2 at corresponding frequencies in the presence of upstream PCWs compared to the median power in the absence of upstream PCWs. These results demonstrate the prevalence of energy deposits into the dayside Martian ionosphere from the solar wind mediated by the PCW-driven ringing of the magnetosphere. Future studies, possibly with new multi-point observations, should address the detailed processes of wave propagation and energy transport through the system and the long-term impact of this chain of processes on the planetary ion heating in the ionosphere and atmospheric loss from Mars.

Occurrence of Kelvin–Helmholtz Instability at Lunar Distance Magnetopause: ARTEMIS Observation

Kelvin–Helmholtz waves can be observed frequently at the near-Earth magnetopause and play an important role in the transport of particles, momentum, and energy from the solar wind to the magnetosphere. This work analyzes the occurrence of Kelvin–Helmholtz instability (KHI) at lunar distance magnetopause, which has not been thoroughly studied currently based on Acceleration, Reconnection, Turbulence, and Electrodynamics of the Moon's Interaction with the Sun satellite observations, and it also investigates the effect of the upstream solar wind and interplanetary magnetic field (IMF). Statistical results show that (1) the occurrence rate is about 15% of the time at lunar distance, lower than at the flank magnetopause, and (2) the occurrence rate decreases with the magnetoacoustic Mach number, Alfvén Mach number, solar wind velocity, and dynamic pressure but only shows a slightly positive correlation with solar wind density. Unlike at the dayside magnetopause, the occurrence rate of KHI diminishes as the solar wind velocity increases at the lunar distance magnetopause, and (3) the occurrence rate decreases with IMF amplitude and is influenced by IMF orientation. As a function of the IMF clock angle, the occurrence rate reaches its maximum at ∼24% when the clock angle is zero. The statistical results are basically consistent with the currently accepted linear theory of KHI, except for a lower rate for higher-speed solar wind. This work contributes to understanding the excitation and evolution of KHI along the magnetopause and plasma transport process in the tail magnetopause.

Magnetohydrodynamic Modeling of Background Solar Wind near Mars: Comparison with MAVEN and Tianwen-1

Combined with data assimilation methods, a three-dimensional magnetohydrodynamic (MHD) numerical model is an effective tool to explore the mechanism of space weather. As a driver of space weather, the dynamic development of stream interaction regions (SIRs) near the orbit of Mars is an area of active research. In this study, we use the interplanetary total variation diminishing (TVD) MHD model to simulate solar wind parameters and model SIRs near Mars from 2021 November 15 to 2021 December 31. In this model, the MHD equations are solved by the conservation TVD Lax–Friedrichs scheme in a rotating spherical coordinate system with six component meshes used on the spherical shell. Solar wind velocity, density, temperature, and magnetic field strength are given at the inner boundary due to the characteristic waves propagating outward. We compared modeled results with observations from Mars Atmospheric Volatile EvolutioN (MAVEN) and Tianwen-1 (China’s first Mars exploration mission). Statistical analysis shows that the simulated results can capture SIRs and are in good agreement with observations; moreover, the assimilated results based on the Kalman filter improve the accuracy of numerical prediction compared with simulated results. This paper is the first attempt to simulate SIR events combined with MAVEN and Tianwen-1 in situ observations. Our work demonstrates that using the MHD model with the Kalman filter to reconstruct solar wind parameters can help us study the characteristics of SIRs near Mars, improve the capabilities of space weather forecasting, and understand the background solar wind environment.

Azimuthal Size Scales of Solar Wind Periodic Density Structures

Periodic density structures (PDSs) are quasiperiodic variations of solar wind density ranging from a few minutes to a few hours. PDSs advect with the solar wind and have radial length scales (L x ) of tens to several thousand megameters, thus belonging to the class of “mesoscale structures.” Current interplanetary multispacecraft observations are not at spatial separations capable of directly measuring the 3D size scale of PDSs or other mesoscale structures. Instead, previous investigations estimated characteristic spatial scales in solar wind parameters using cross-correlation and/or coherence analysis applied to multispacecraft observations. For the solar wind density and interplanetary magnetic field (IMF) intensity, the reported size scales perpendicular to the Sun–Earth line (L y ) ranged between ≈30 and ≈200 Earth Radii (R E). Here, we implemented a similar approach for the same parameters, but focused on high-density, slow-solar-wind intervals with PDSs observed by the Wind and ARTEMIS-P1 spacecraft. Additionally, this is the first statistical study of the IMF intensity periodicities in relation to PDSs. We identified intervals in which the two spacecraft observed the same periodicity, obtaining two PDS groups based on their radial length scale: L x1 ≈ 86R E and L x2 ≈ 35R E. Then, we classified the events based on the periodic variations’ coherence level. Reproducing the results with simulations of the PDSs’ transit, we inferred the L y order of magnitudes for the two PDS groups: L y1 ≈ 340R E and L y2 ≈ 187R E. Knowing the PDSs’ size scales is fundamental for constraining models aimed at reproducing these structures and is critical for better understanding the PDS–magnetosphere coupling.

Magnetic Storm‐Time Red Aurora as Seen From Hokkaido, Japan on 1 December 2023 Associated With High‐Density Solar Wind

AbstractWe report a citizen science‐motivated study on the cause of an unusually bright red aurora as witnessed from Hokkaido, Japan during a magnetic storm on 1 December 2023. The auroral brightness of 5 kR is unusual for the Dst index peak of only −107 nT. In spite of the moderate storm amplitude, the extremely high solar wind density of >50/cc and dynamic pressure of >25 nPa caused the aurora oval extension to 53 magnetic latitudes (L = 2.8). We discuss that the drift loss of the ring current particles across the small‐size magnetopause is important, and Hokkaido was at the right position to see the direct effect of the large particle injection of the storm‐time substorm.

A Pressure Pulse‐Driven Transient Magnetospheric Event

AbstractBursty reconnection models predict that flux transfer events (FTEs) moving along the magnetopause launch fast mode compressional waves into the magnetosheath that push the bow shock outward. By contrast, increases in the solar wind density striking the bow shock should push that boundary inward and launch fast mode compressional waves that propagate across the magnetosheath, drive waves on the magnetopause, and generate transient events in the outer magnetosphere. Multipoint ACE, Wind, THEMIS, and GOES‐11/12 solar wind, bow shock, and magnetospheric observations on 14 October 2008 provide direct evidence for solar wind pressure pulses producing a large amplitude indentation with crater FTE‐like properties on the magnetopause.

Orientation of IMF Discontinuity Normals Across the Solar Cycles

AbstractTransient events like hot flow anomalies and foreshock bubbles are common in the Earth's foreshock. These foreshock transients may play an important role in the solar wind‐magnetosphere interaction. They typically occur when backstreaming ions accumulate at the intersection of interplanetary magnetic field (IMF) discontinuities with the bow shock. Discontinuity orientations play a key role in determining when and where transients form in the foreshock and strike the magnetopause. They also play a role in determining the amplitude, and significance, of individual transients. We investigate properties of the IMF discontinuity normals across two solar cycles. We compare Advanced Composition Explorer (ACE) IMF discontinuity observations for 2001, 2002, 2012, and 2014 (solar maximum) and 2008, 2009, 2019, and 2020 (solar minimum). We employ both the Minimum variance analysis (MVA) and Cross‐product methods to determine phase front normals. Most discontinuity normals point earthward and dawnward at longitudes ranging from 180° to 250°. A greater range of normal can be seen during solar minimum than during solar maximum. Normals lie further from the Sun‐Earth line when solar wind densities are low and velocities are high. This provides a natural explanation for the observed tendency of transients to occur for low solar wind densities and high velocities in terms of longer discontinuity interaction times with the bow shock.

WawHelioIonMP: A Semiempirical Tool for the Determination of Latitudinal Variation in the Ionization Rate of Interstellar Hydrogen and the Solar Wind

The latitudinal structure of the solar wind varies during the cycle of solar activity. Analysis of this variation is important for understanding the solar activity and interpretation of observations of heliospheric energetic neutral hydrogen atoms and interstellar neutral (ISN) atoms inside the heliosphere, which yield information on the heliosphere and its interaction with the interstellar medium. Existing methods of retrieving this information from indirect remote-sensing measurements of phenomena, including the heliospheric backscatter glow and interplanetary scintillations of remote radio sources, are challenging to apply in real time. Here, we propose a method WawHelioIonMP of approximate retrieval of latitudinal profiles of the ionization rates of ISN H using a machine-learning-based interpretation of the helioglow. Assuming that we know their history during two past solar cycles and have observations of the helioglow for close-to-circumsolar circles with a radius close to 90°, we derive statistically an algebraic relation between the ionization profiles and lightcurves. With the relation reversed, we are then able to derive the ionization rate profiles based on observed light curves, such as those planned for the GLObal solar Wind Structure (GLOWS) experiment on the forthcoming NASA mission Interstellar Mapping and Acceleration Probe (IMAP). The application of this method is straightforward and rapid because complex simulations are no longer needed. We present the method of retrieval of the profiles of the ionization rates, leaving the discussion of details of the decomposition of the retrieved ionization rate profiles into profiles of the solar wind speed and density to a future paper.

Sliding-window cross-correlation and mutual information methods in the analysis of solar wind measurements

Context. When describing the relationships between two data sets, four crucial aspects must be considered, namely: timescales, intrinsic lags, linear relationships, and non-linear relationships. We present a tool that combines these four aspects and visualizes the underlying structure where two data sets are highly related. The basic mathematical methods used here are cross-correlation and mutual information (MI) analyses. As an example, we applied these methods to a set of two-month’s worth of solar wind density and total magnetic field strength data. Aims. Two neighboring solar wind parcels may have undergone different heating and acceleration processes and may even originate from different source regions. However, they may share very similar properties, which would effectively “hide” their different origins. When this hidden information is mixed with noise, describing the relationships between two solar wind parameters becomes challenging. Time lag effects and non-linear relationships between solar wind parameters are often overlooked, while simple time-lag-free linear relationships are sometimes insufficient to describe the complex processes in space physics. Thus, we propose this tool to analyze the monotonic (or linear) and non-monotonic (or non-linear) relationships between a pair of solar wind parameters within a certain time period, taking into consideration the effects of different timescales and possible time lags. Methods. Our tool consists of two parts: the sliding-window cross-correlation (SWCC) method and sliding-window mutual information (SWMI) method. As their names suggest, both parts involve a set of sliding windows. By independently sliding these windows along the time axis of the two time series, this technique can assess the correlation coefficient (and mutual information) between any two windowed data sets with any time lags. Visualizing the obtained results enables us to identify structures where two time series are highly correlated, while providing information on the relevant timescales and time lags. Results. We applied our proposed tool to solar wind density and total magnetic field strength data. Structures with distinct timescales were identified. Our tool also detected the presence of short-term anti-correlations coexisting with long-term positive correlations between solar wind density and magnetic field strength. Some non-monotonic relationships were also found. Conclusions. The visual products of our tool (the SWCC+SWMI maps) represent an innovative extension of traditional numerical methods, offering users a more intuitive perspective on the data. The SWCC and SWMI methods can be used to identify time periods where one parameter has a strong influence on the other. Of course, they can also be applied to other data, such as multi-wavelength photometric and spectroscopic time series, thus providing a new tool for solar physics analyses.

Solar Wind Driven from GONG Magnetograms in the Last Solar Cycle

In a previous study, Huang et al. used the Alfvén Wave Solar atmosphere Model, one of the widely used solar wind models in the community, driven by ADAPT-GONG magnetograms to simulate the solar wind in the last solar cycle and found that the optimal Poynting flux parameter can be estimated from either the open field area or the average unsigned radial component of the magnetic field in the open field regions. It was also found that the average energy deposition rate (Poynting flux) in the open field regions is approximately constant. In the current study, we expand the previous work by using GONG magnetograms to simulate the solar wind for the same Carrington rotations and determine if the results are similar to the ones obtained with ADAPT-GONG magnetograms. Our results indicate that similar correlations can be obtained from the GONG maps. Moreover, we report that ADAPT-GONG magnetograms can consistently provide better comparisons with 1 au solar wind observations than GONG magnetograms, based on the best simulations selected by the minimum of the average curve distance for the solar wind speed and density.

Large TEC Variations Between Mars and Earth: Simulation and Observation Comparisons

AbstractLarge total electron content (TEC) variations along the line from Mars to Earth have been demonstrated by analyzing the ground received Tianwen‐1 differential one‐way range (DOR) signals when corona mass ejections (CMEs) and co‐rotating interaction regions (CIRs) pass through the signal path. Here, the TEC variations along the line from Mars to Earth are calculated from the Wang‐Sheeley‐Arge (WSA)‐Enlil space weather forecast model. The results are compared with the Tianwen‐1 line of sight (LOS) observations as well as the STEREO‐A in situ solar wind density measurements. It is found that the TEC variations calculated from the WSA‐Enlil model are usually much smaller in amplitude, especially in the cases of CIRs passing through the signal path. In addition, the timings of CMEs and CIRs passing through the DOR signal path in the WSA‐Enlil model can be half a day earlier or later than the observed ones. In comparison with in situ solar wind measurements, the ground received Tianwen‐1 DOR signals carry solar wind density information over a large spatial region along the signal path. Such remote sensing measurement can be valuable for constraining and improving global solar wind forecast models.

Solar Wind Density and Core Temperature Derived from the PSP Quasi-thermal Noise Measurements

Quasi-thermal noise (QTN) spectroscopy is a valuable method to deduce important parameters in space plasma, such as plasma density and temperature, especially when direct particle measurements are not available. The present study develops a new fitting method to fit the QTN spectra observed by the Parker Solar Probe (PSP) with a comprehensive theoretical QTN spectral model. By combining the steepest descent and Levenberg–Marquardt algorithms, the new method is more flexible with initial guess values but still yields reliable solar wind electron density and temperature values. The new method is applied to derive the solar wind density and core temperature from the QTN measurements during 10 encounters of PSP. The electron density and temperature values obtained vary with the radial distance from the Sun as n e ∝ r −2.12 and T e ∝ r −0.71, both of which are consistent with existing models and previous results.

The Martian Ionospheric Response to the Co‐Rotating Interaction Region That Caused the Disappearing Solar Wind Event at Mars

AbstractAn unusually low density solar wind event was observed in December 2022 moving past both Earth and Mars. The source was traced back to a coronal hole and active region on the Sun's surface. The resulting solar wind lead to the development of a co‐rotating interaction region (CIR) and trailing rarefaction region that lasted for multiple solar rotations. Within this structure, the solar wind conditions, including density, velocity, and magnetic field magnitude and orientation drastically changed. In this study we analyze the response of the Martian ionosphere using MAVEN data to these changing solar wind conditions. The low density solar wind region associated with the December event resulted in the expansion of the Martian ionospheric boundaries. We show that the ion composition boundary (ICB) is located at extreme altitudes that are beyond previously observed locations from the MAVEN mission between 2015 and 2018. Furthermore, the boundary between shocked solar wind and the Martian ionosphere identified using electron and ion data moved together on the dayside of the planet with the changing solar wind conditions. However, at the flank region these boundaries do not move together, and we show here that the decoupling of the two boundaries may be the result of a change in the interplanetary magnetic field azimuthal angle.

Simulating Compressive Stream Interaction Regions during Parker Solar Probe’s First Perihelion Using Stream-aligned Magnetohydrodynamics

We used the stream-aligned magnetohydrodynamics (SA-MHD) model to simulate Carrington rotation 2210, which contains Parker Solar Probe’s (PSP) first perihelion at 36.5 R ⊙ on 2018 November 6, to provide context to the in situ PSP observations by FIELDS and SWEAP. The SA-MHD model aligns the magnetic field with the velocity vector at each point, thereby allowing for clear connectivity between the spacecraft and the source regions on the Sun, without unphysical magnetic field structures. During this Carrington rotation, two stream interaction regions (SIRs) form, due to the deep solar minimum. We include the energy partitioning of the parallel and perpendicular ions and the isotropic electrons to investigate the temperature anisotropy through the compression regions to better understand the wave energy amplification and proton thermal energy partitioning in a global context. Overall, we found good agreement in all in situ plasma parameters between the SA-MHD results and the observations at PSP, STEREO-A, and Earth, including at PSP’s perihelion and through the compression region of the SIRs. In the typical solar wind, the parallel proton temperature is preferentially heated, except in the SIR, where there is an enhancement in the perpendicular proton temperature. This is further showcased in the ion cyclotron relaxation time, which shows a distinct decrease through the SIR compression regions. This work demonstrates the success of the Alfvén wave turbulence theory for predicting interplanetary magnetic turbulence levels, while self-consistently reproducing solar wind speeds, densities, and overall temperatures, including at small heliocentric distances and through SIR compression regions.

Constraints on Solar Wind Density and Velocity Based on Coronal Tomography and Parker Solar Probe Measurements

Previous work has established an empirical relationship between densities gained from coronal rotational tomography near the ecliptic plane with solar wind outflow speeds at heliocentric distance r 0 = 8R ⊙. This work aims to include solar wind acceleration, and thus velocity profiles out to 1 au. Inner boundary velocities are given as a function of normalized tomographic densities, ρ N , as V0=75∗e−5.2*ρN+108 , and typically range from 100 to 180 km s−1. The subsequent acceleration is defined as V(r)=V01+αIP1−e−r−r0/rH , with α IP ranging between 1.75 and 2.7, and r H between 50 and 35 R ⊙ dependent on V 0. These acceleration profiles approximate the distribution of in situ measurements by Parker Solar Probe (PSP) and other measurements at 1 au. Between 2018 November and 2021 September these constraints are applied using the HUXt model and give good agreement with in situ observations at PSP, with a ∼6% improvement compared with using a simpler constant acceleration model previously considered. Given the known tomographical densities at 8 R ⊙, we extrapolate density to 1 au using the model velocities and assuming mass flux conservation. Extrapolated densities agree well with OMNI measurements. Thus coronagraph-based estimates of densities define the ambient solar wind outflow speed, acceleration, and density from 8 R ⊙ to at least 1 au. This sets a constraint on more advanced models, and a framework for forecasting that provides a valid alternative to the use of velocities derived from magnetic field extrapolations.