- Predictive Science Inc., San Diego, CA, United States
Understanding how coronal structure propagates and evolves from the Sun and into the heliosphere has been thoroughly explored using sophisticated MHD models. From these, we have a reasonably good working understanding of the dynamical processes that shape the formation and evolution of stream interaction regions and rarefactions, including their locations, orientations, and structure. However, given the technical expertize required to produce, maintain, and run global MHD models, their use has been relatively restricted. In this study, we refine a simple Heliospheric eXtrapolation Technique (HUX) to include not only forward mapping from the Sun to 1 AU (or elsewhere), but backward mapping toward the Sun. We demonstrate that this technique can provide substantially more accurate mappings than the standard, and often applied “ballistic” approximation. We also use machine learning (ML) methods to explore whether the HUX approximation to the momentum equation can be refined without loss of simplicity, finding that it likely provides the optimum balance. We suggest that HUX can be used, in conjunction with coronal models (PFSS or MHD) to more accurately connect measurements made at 1 AU, Stereo-A, Parker Solar Probe, and Solar Orbiter with their solar sources. In particular, the HUX technique: 1) provides a substantial improvement over the “ballistic” approximation for connecting to the source longitude of streams; 2) is almost as accurate, but considerably easier to implement than MHD models; and 3) can be applied as a general tool to magnetically connect different regions of the inner heliosphere together, as well as providing a simple 3-D reconstruction.
1 Introduction
Plasma is heated in the corona and accelerates away from it to form the solar wind. It is convenient (although probably an oversimplification) to separate what we believe to be intrinsically spatial variations from temporal variations [1]. The former is the repeating large-scale structure we observe in the interplanetary medium from one rotation to the next, while the latter manifests as sporadic eruptions of transient phenomena, including coronal mass ejections (CMEs) as well a rich set of other effects, such as jets, plumes, small-scale flux ropes, etc., We recognize, however, that even the structure we label as “spatial”, is undoubtedly driven, at least to some extent, by time-dependent phenomena. Nevertheless, it remains useful to treat the large-scale, quasi-corotating structure as a predominantly time-stationary process. This is supported by many observations over the last 50+ years, which show that, particularly during the declining phase of the solar activity cycle, and the interval surrounding solar minimum, in particular, global models relying on steady-state solutions can capture the essential features of solar wind streams (e.g., [2–4]).
Being able to connect in situ measurements back to their solar sources, or explore the evolution of plasma as it propagates away from the Sun is both scientifically and operationally important. For example, the source of the slow solar wind remains the topic of considerable debate; thus, being able to estimate from which regions on the Sun specific plasma came from could shed light on its origin (e.g., [5]). Operationally, both ambient solar wind and empirically based coronal mass ejections (CMEs) can be mapped from the Sun to the Earth to improve space weather forecasts (e.g., [6]). Over the years, a number of techniques have been developed to map solar wind plasma from one point in the solar wind to another. The simplest is the so-called “ballistic” approximation [7], which assumes that parcels of plasma travel radially, and maintain a constant speed as they propagate from one location to another (either inwards or outwards). The technique is simple to implement, but can lead to substantial errors, particularly at compression regions, where the interaction between plasma of different speeds is dynamic. Sophisticated global MHD models (e.g., [8], address these limitations by providing three-dimensional dynamic evolution of the plasma. However, they are difficult to develop and require significant resources to run them.
[9]; herein referred to as Paper 1, introduced the concept of Heliospheric Upwinding eXtrapolation (HUX). This simplified method bridged the gap between kinematic (i.e., constant speed) mapping (e.g., [7]) and full MHD calculations (e.g., [2]), providing a surprisingly accurate extrapolation of solar wind speeds close to the Sun to 1 AU or beyond. By neglecting the pressure gradient and gravity terms, they were able to show that the momentum equation reduced to the inviscid Burger’s equation, which could be easily converted to a difference algorithm and used to evolve speed profiles near the Sun radially outwards. Additionally, to account for the general acceleration of the solar wind near the Sun, a simple expression was derived and applied. Together, these two components allowed MHD solutions to be reproduced with high accuracy (Pearson correlation coefficient,
Several studies have leveraged the HUX technique to improve space weather forecasts of solar wind streams. [6]; for example, developed ensemble solutions of solar wind conditions at Earth by sampling a range of latitudes about the sub-Earth point, resulting in more accurate forecasts and variances that accurately captured the model uncertainty [10]. used HUX to propagate several competing models for defining the speed of the solar wind in the high corona out to 1 AU, allowing them to assess forecasting performance among them. They also found that the ability to incorporate multiple realizations (which would not have been computationally feasible with global MHD model solutions), they were able to improve forecasting performance of all the models considered. More recently [11], refined HUX by generalizing it to a form similar to that of the viscous Burger’s equation (Tunable HUX, or THUX), essentially adding an additional term to mimic the effects of viscosity. It is worth noting that this also introduced a new free parameter, η, that can be tuned to make the model comparisons better, and, the changes in the forecasted speeds were found to be modest, and not obviously improvements (see their Figure 1) [12]. combined the HUX approach with a machine learning (ML) algorithm (specifically, Gradient Boosting Regressors), to produce a fast and accurate model for forecasting ambient solar wind conditions at Earth.
The HUX technique is also being incorporated into space weather programs. [13]; for example, have integrated HUX into a pilot program for developing space weather capabilities for the Indian Space Research Organization (ISRO). Additionally [14], have used HUX within an ensemble-based CME arrival prediction model (ELEvoHI), finding that the combination of the WSA model for prescribing the speed at the inner boundary and HUX resulted in a reduction in the Mean Absolute Error (MAE) by almost 2 h, in comparison to other techniques using the observed solar wind speed at L1.
The HUX technique has also been generalized to study time-dependent phenomena [15]. interpreted the
In this study, we describe a simple refinement to the HUX allowing the user to map solar wind streams from 1 AU (or elsewhere) back to the Sun, which can have a substantial impact on the inferred source longitudes of the observed plasma. Additionally, using a data-driven sparse regression method, we explore whether there are any simple improvements that can be applied to the HUX approach. Finally, we explore the possible impacts of including differential rotation, as well as any constraints imposed by resolution.
2 Models
Here, we introduce the three main modeling techniques that are applied in our study. Specifically: 1) ballistic mapping; 2) MHD modeling; and 3) the HUX technique.
2.1 Ballistic Mapping
The simplest (but least useful) approximation we could make for evolving the speed of the solar wind as it propagates away from the Sun is to assume that it does not change speed in response to dynamical interactions between adjacent parcels of plasma. This has a number of obvious problems, perhaps largest of which is that in mapping the speeds out, it becomes possible for faster parcels to outrun slower ones. By virtue of the fact that the solar wind can be reasonably approximated as a fluid, this is patently nonphysical. When the procedure is instead applied in the reverse direction, this can lead to the well known problems of “dwells”, where parcels of plasma observed at a later times in the solar wind map back to earlier launch times at the Sun, again by apparently “crossing paths” [16].
In spite of its simplicity, for completeness, we formally define the ballistic approximation as:
which simply states that the speed at
where
2.2 Magnetohydrodynamic Model
For our numerical experiments, we use the Magnetohydrodynamic Algorithm outside a Sphere (MAS) code, which solves the time-dependent resistive magnetohydrodynamic (MHD) equations (e.g., [8]). Although the model can be run with a range of increasingly complex energy transport processes, for the purposes of our study, we rely on either polytropic or thermodynamic solutions [4]. The inner boundary of the calculation is set to
2.3 Heliospheric Upwinding eXtrapolation Technique
In an earlier study [9], herein referred to as paper 1, we developed a simple mapping technique that approaches the simplicity of the ballistic approach but retains most of the accuracy of the MHD approach, by reducing the momentum equation to the much simpler Burger equation.
Briefly, the solar wind motion can be described as the fluid momentum equation in a corotating frame of reference:
where ρ is the proton mass density, V is the velocity, p is the thermal pressure, G is the gravitational constant, Ms is the mass of the Sun, and γ is the polytropic index. This contrasts the more common form of the momentum equation in that the time derivative,
Using the forward upwind difference algorithm, we can represent this as:
We denote this form as HUX-backward, or HUX-b.
As discussed in more detail in paper 1, we can apply an acceleration term to the forward mapped speeds of the form:
where
Similarly, by inspection, we can define a deceleration term for the HUX-b mapping as follows:
3 Results
3.1 A “Ground Truth” From the Magnetohydrodynamic Results and Comparison With the Ballistic Approximation
As in paper 1, we use a global MHD solution as the “ground truth” for validating the more approximate HUX results. We focus again on CR 2068, but also add a number of other Carrington rotations to the analysis, to ensure that we do not over-generalize the results from one solution. In Figure 2A we summarize the equatorial magnetic field lines from CR 2068. In contrast to the usual way that field lines are drawn by tracing along the magnetic field vectors, in this case, we use the velocity field directly to evolve the field lines since they are convected out along the plasma flow, and there are no time-dependent effects to complicate matters. Thus, strictly speaking, these “field lines” are really streamlines. They are arbitrarily colored, but, importantly, each source longitude is drawn with the same color in subsequent plots, allowing us to directly compare the relative evolution of the “field lines” between approaches.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/679497/fphy-09-679497-HTML-r2/image_m/fphy-09-679497-g001.jpg)
FIGURE 1. (A) Field lines traced out from the inner boundary of the MHD solution using the MHD velocity solution. (B) Field lines traced out from the inner equatorial boundary assuming that the speed for that field line remains the same (i.e., the ballistic approximation).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/679497/fphy-09-679497-HTML-r2/image_m/fphy-09-679497-g002.jpg)
FIGURE 2. Probability density (histogram) of the difference between the ballistically mapped longitudes and the MHD-mapped longitudes of the field lines from the inner boundary at 30Rs to 1 AU. The blue blocks show the actual frequency while the orange curve provides a Gaussian fit to the data. The Gaussian fit is used only as a practical way to estimate the spread in the values, and not to suggest that the errors are normally distributed.
In Figure 2B we have drawn field lines from the same sources under the assumption that the speed does not evolve as the field is dragged radially out with the flow. This is the so-called ballistic approximation. We note that while field lines generally terminate near the correct location, there is considerable variability in the accuracy of the mapping. More importantly, we see that the “field lines” overlap one another; a physically forbidden result.
To more quantitatively assess the difference between the MHD mapping and the ballistic mapping, in Figure 1, we show the difference between the ballistically mapped longitudes and the MHD-mapped longitudes of the field lines. The orange curve is a Gaussian fit to the data to enable an estimate of the statistical parameters for the distribution of errors. We find that the mean/median error in the ballistic mapping is
3.2 Mapping Solar Wind Streams Out From the Sun With HUX-f
Figure 3 (top panel) shows the radial speed at the inner boundary of the MHD calculation (
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/679497/fphy-09-679497-HTML-r2/image_m/fphy-09-679497-g003.jpg)
FIGURE 3. (A) Solar wind radial speed as a function of Carrington longitude and sin (latitude) at the inner boundary of the MHD simulation. (B) solar wind speed at 1 AU as computed by the MHD solution. (C) Solar wind speed at 1 AU as estimated using the HUX-f technique.
Before doing this, however, it is useful to assess how much evolution has taken place from the inner radial boundary to 1 AU in the MHD model. In Figure 4A, we compare the equatorial solar wind at
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/679497/fphy-09-679497-HTML-r2/image_m/fphy-09-679497-g004.jpg)
FIGURE 4. (A) Comparison of solar wind speed at the equator at 30 solar radii (blue) and 1 AU (red) from MHD solution. (B) Comparison of solar speed at the equator at 1 AU from the MHD solution (blue) and HUX-f technique (green).
We explore the errors more quantitatively in Figure 5, which shows the difference between the HUX-f and MHD speeds, normalized to the MHD speeds at 1 AU. On average, the errors associated with the HUX-f approach are 0.086, or approximately 9%. Comparison with Figure 3 also suggests that the sense of the errors are systematic, with positive errors associated with compression regions, and negative errors associated with rarefactions. This is reinforced qualitatively in Figure 4, which, as already noted, showed that the HUX-f technique underestimates the peak of the high-speed streams, but also tends to overestimate the troughs. Nevertheless, the key point is that the average errors are substantially smaller than, say, a ballistic approximation would produce.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/679497/fphy-09-679497-HTML-r2/image_m/fphy-09-679497-g005.jpg)
FIGURE 5. Estimate of the error in the HUX-f mapping as compared to the MHD result at 1 AU. Positive errors are shown in blue, negative errors are red, and an exact match appears as white.
As a final consideration of the errors associated with the mapped speeds, in Figure 6A we show the point-by-point relationship between all grid points (i.e., all latitudes and longitudes) at 1 AU for the HUX-f and MHD solutions. The computed Pearson correlation coefficient for this was found to be 0.997. However, it should be noted that this is likely heavily biased by the large fraction of the solar wind volume that is not undergoing any significant evolution, that is, the fastest and slowest wind pinning the extremes of the correlation. A more appropriate estimate is to limit the comparison to the equator. This is accomplished in Figure 6B, where the equatorial values as well as one slice on either side have been plotted. The correlation coefficient for this subset, while still remarkably high (0.96) is lower than for the full comparison.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/679497/fphy-09-679497-HTML-r2/image_m/fphy-09-679497-g006.jpg)
FIGURE 6. Point by point scatter plot between MHD and HUX-f results. (A) For all latitudes, the Pearson correlation coefficient (PCC) was 0.994. (B) For a narrow band about the heliographic equator, the PCC was 0.958.
3.3 Outward Extrapolation of Field Lines Using HUX-f
We can also map out field lines using the HUX-f technique and compare them with the MHD results. In Figure 7B we have drawn field lines from the same sources as in Figure 1. At least qualitatively, comparison between panels (a) and (b) suggests a very close match between the MHD and HUX-f approaches.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/679497/fphy-09-679497-HTML-r2/image_m/fphy-09-679497-g007.jpg)
FIGURE 7. (A) Field lines traced out from the inner boundary of the MHD solution using the MHD velocity solution. (B) Field lines traced out from the inner equatorial boundary assuming that the speed for that field line remains the same.
Again, to more quantitatively assess the difference between the MHD mapping and the HUX-f mapping, in Figure 8, we show the difference between the HUX-f-mapped longitudes and the MHD-mapped longitudes of the field lines. This time, we find that the mean/median error in the ballistic mapping is
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/679497/fphy-09-679497-HTML-r2/image_m/fphy-09-679497-g008.jpg)
FIGURE 8. Probability density (histogram) of difference between the ballistically mapped longitudes and the MHD-mapped longitudes of the field lines. The blue blocks show the actual frequency while the orange curve provides a Gaussian fit to these data.
3.4 Inward Extrapolation of Field Lines Using Heliospheric Upwinding eXtrapolation-B
Finally, we can repeat the analysis using the HUX-b technique. In this case, we start the mapping at 1 AU and draw the field lines back to the Sun. In Figure 9B we have drawn field lines from the same MHD sources at 1 AU using the HUX-b formulation. Again, at least qualitatively, comparison between panels (a) and (b) suggest a close match between the MHD and HUX-b approaches. However, it can be noted that in the high-speed streams (where the field lines are more radial and show more space between them), the HUX-b field lines are not separated to the same degree. Additionally, where the field lines are closer together, the effect is reduced in the HUX-b solution. We can interpret this in terms of compression and rarefaction regions, suggesting that the HUX-b solutions would produce weaker compression regions and less tenuous rarefaction regions that would, in reality, be present. Nevertheless, a novel application of such a visualization is that it allows one to infer the location of compression/rarefaction regions with only a knowledge of the velocity field.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/679497/fphy-09-679497-HTML-r2/image_m/fphy-09-679497-g009.jpg)
FIGURE 9. (A) Field lines traced out from the inner boundary of the of the MHD solution using the MHD velocity solution. (B) Field lines traced out from the inner equatorial boundary assuming that the speed for that field line remains the same.
Considering the distribution in the errors, in Figure 10, we show the difference between the HUX-b-mapped longitudes and the MHD-b-mapped longitudes of the field lines. This time, we find that the mean/median error in the ballistic mapping is
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/679497/fphy-09-679497-HTML-r2/image_m/fphy-09-679497-g010.jpg)
FIGURE 10. Probability density (histogram) of difference between the ballistically mapped longitudes and the MHD-mapped longitudes of the field lines. The blue blocks show the actual frequency while the orange curve provides a Gaussian fit to these data.
3.5 Potential Improvements to the Heliospheric Upwinding eXtrapolation Approach
The HUX approach is simple. In fact, it is probably the simplest physically based improvement to a ballistic approximation that could be made. As noted earlier, several refinements have been made, principally to better address the mixing of temporal and spatial variations at the source. However, here we would like to ask specific questions. First, does differential rotation affect the tuning of the free parameters in the HUX model? Second, does the inclusion of more terms from the momentum equation would improve the approach? And third, is the HUX mapping sensitive to grid resolution?
3.5.1 Effects of Differential Rotation
In paper 1, we estimated the two free parameters of the HUX technique (
where θ is colatitude. Thus, at the equator, this produces a rotation rate corresponding to a period of 25.38 days, while at the pole, the period is 31.5 days.
To estimate the error associated with these solutions, we define the residuals between the MHD model (the “data”, or “ground truth” in this case) and the HUX model to be:
and use the mean square error (MSE) as a measure of the goodness-of-fit of the model:
For CR 2068 (Figure 3), using the values of α and
3.5.2 Incorporation of the Pressure Term Into the Heliospheric Upwinding eXtrapolation Model
As shown in paper 1, we could reduce the momentum equation to the inviscid Burgers’ equation under the assumption that the magnetic field, pressure gradient, and gravity can be neglected. But, to what extent is this a reasonable approximation, or, stated another way, what errors are likely introduced? As a secondary question, we can ask whether there is a more speculative formalism that might produce better mappings of the streams from the Sun to 1 AU?
To investigate this, we consider the following generalized expression for
where:(14) Here,
To fully explore the parameter space defined by Eq. 14, we employed the PDE-FIND package, developed by [20]; which searches through a library of terms to fit the optional subset for time-series-like data (or model results).
The details of the analysis are provided in the supplemental information (SI, GitHub), but briefly, we used ridge sparse regression to find the optimal number of terms:
The results (SI) showed that the terms
Finally, it is worth noting that PDE-FIND is a physics-agnostic approach, thus, our inclusion of various combinations of first and second-order partial derivatives of was a useful test for validating the approach since it correctly and independently identified the term
3.5.3 Numerical Convergence of the Heliospheric Upwinding eXtrapolation Mapping
Although the HUX approach is extremely simple, it is still derived from a partial differential equation, and subject to potential convergence issues, such as the Courant–Friedrichs–Lewy (CFL) condition. For the HUX technique, this requires that:
In Figure 11 we compare speed traces as a function of longitude for a range of radial resolutions, from a maximum of
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/679497/fphy-09-679497-HTML-r2/image_m/fphy-09-679497-g011.jpg)
FIGURE 11. Comparison of MHD model speeds with HUX-f solutions of different radial resolution. In this case
4 Conclusions and Discussion
In this study, we have further developed the Heliospheric Upwinding eXtrapolation (HUX) technique. Specifically, we have: 1) Demonstrated how the approach can be used to map solar wind streams back to the Sun; and 2) Shown that the current formalism is probably as complicated as it needs to be to produce useful results. In addressing these questions, we were able to show how much the HUX technique would improve the accuracy of ballistic mapping studies, which are typically used to identify the source regions of solar wind streams or energetic particles, for example.
This study is not without limitations. In particular, we used MHD model results as the “ground truth”. Although this is reasonable in the sense that we have no other global dataset of solar wind speeds, it assumes that the MHD formalism is accurate enough to evolve solar wind streams through the heliosphere. This, of course breaks down at sufficiently high frequencies. Thus we can only claim that the HUX approach is reasonable on the largest spacial and temporal scales (i.e., macroscopic structure). Additionally, any artifacts introduced by the MHD model, such as numerical diffusion, which are also mimicked by the HUX technique would also contribute to a higher accuracy of the results, but which might not exist in practice. However, given the many studies that have validated the MHD approach for studying solar wind evolution (e.g., [2, 18, 21, 22]), it seems reasonable to conclude that the MHD solutions provide a sufficiently good “ground truth” for such tests.
It is worth noting that our application of PDE-FIND included the viscous term (
In closing, we reiterate that in this study, we have focused on the procedure that could be applied to various in situ datasets, allowing them to be mapped back toward the Sun, or outward from one spacecraft location to another. With the successful launch and commissioning of PSP and Solar Orbiter, as well as the presence of 1 AU spacecraft, including ACE and DSCOVR at Earth, and STEREO-A offset from the Earth Sun line to varying degrees, the technique we have described should be a useful tool for exploring the evolution of streams from one location to another. The approach can also be applied to a broad set of historical datasets including Helios, Ulysses, Pioneer, and Voyagers 1 and 2, just to mention a few. In a future study, we will describe the application of HUX-b and HUX-f to a variety of heliospheric in situ datasets.
Data Availability Statement
The model results used in this study can be found in the HUX repository, located at https://github.com/predsci/HUX.
Author Contributions
PR developed the theoretical formalism for the ideas presented here, with contributions from OI. OI wrote the Python code to test the concepts with minor contributions from PR. PR wrote the paper with input from OI.
Funding
The authors gratefully acknowledge support from NASA (80NSSC18K0100, NNX16AG86G, 80NSSC18K1129, 80NSSC18K0101, 80NSSC20K1285, 80NSSC18K1201, and NNN06AA01C)
Conflict of Interest
Authors PR and OI were employed by Predictive Science Inc.
Acknowledgments
Additionally, the authors would like to express their gratitude to both reviewers for constructive comments that improved the clarity of the manuscript.
References
1. Gosling JT. Corotating and Transient Solar Wind Flows in Three Dimensions. Annu Rev Astron Astrophys (1996) 34:35–73. doi:10.1146/annurev.astro.34.1.35
2. Riley P, Linker JA, Mikić Z. An Empirically-Driven Global MHD Model of the Solar Corona and Inner Heliosphere. J Geophys Res (2001) 106:15889–901. doi:10.1029/2000JA000121
3. Riley P, Lionello R, Linker JA, Mikic Z, Luhmann J, Wijaya J. Global MHD Modeling of the Solar Corona and Inner Heliosphere for the Whole Heliosphere Interval. Sol Phys (2011b) 274:361–77. doi:10.1007/s11207-010-9698-x
4. Riley P, Downs C, Linker JA, Mikic Z, Lionello R, Caplan RM. Predicting the Structure of the Solar Corona and Inner Heliosphere during Parker Solar Probe 's First Perihelion Pass. Astrophysical J Lett (2019) 874:L15. doi:10.3847/2041-8213/ab0ec3
5. Riley P, Linker JA, Arge CN. On the Role Played by Magnetic Expansion Factor in the Prediction of Solar Wind Speed. Space Weather (2015) 13:154–69. doi:10.1002/2014sw001144
6. Owens MJ, Riley P. Probabilistic Solar Wind Forecasting Using Large Ensembles of Near-Sun Conditions with a Simple One-Dimensional "Upwind" Scheme. Space Weather (2017) 15:1461–74. doi:10.1002/2017sw001679
7. Snyder CW, Neugebauer M. The Relation of Mariner-2 Plasma Data to Solar Phenomena. In: RJ Mackin, and M Neugebauer, editors. The Solar Wind. Oxford: Permanon Press (1966). p. 25.
8. Riley P, Linker JA, Lionello R, Mikic Z. Corotating Interaction Regions during the Recent Solar Minimum: The Power and Limitations of Global MHD Modeling. J Atmos Solar-Terrestrial Phys (2012a) 83:1–10. doi:10.1016/j.jastp.2011.12.013
9. Riley P, Lionello R. Mapping Solar Wind Streams from the Sun to 1 AU: A Comparison of Techniques. Sol Phys (2011) 270:575–92. doi:10.1007/s11207-011-9766-x
10. Reiss MA, MacNeice PJ, Mays LM, Arge CN, Möstl C, Nikolic L, et al. Forecasting the Ambient Solar Wind with Numerical Models. I. On the Implementation of an Operational Framework. Astrophysical J Suppl Ser (2019) 240:35. doi:10.3847/1538-4365/aaf8b3
11. Reiss MA, MacNeice PJ, Muglach K, Arge CN, Möstl C, Riley P, et al. Forecasting the Ambient Solar Wind with Numerical Models. Ii. An Adaptive Prediction System for Specifying Solar Wind Speed Near the Sun. Astrophysical J (2020) 891:165. doi:10.3847/1538-4357/ab78a0
12. Bailey RL, Reiss MA, Arge CN, Möstl C, Henney CJ, Owens MJ, et al. Using gradient boosting regression to improve ambient solar wind model predictions. Space Weather (2021). doi:10.1029/2020SW002673
13. Kumar S, Paul A, Vaidya B. A Comparison Study of Extrapolation Models and Empirical Relations in Forecasting Solar Wind. Front Astron Space Sci (2020) 7:92. doi:10.3389/fspas.2020.572084
14. Amerstorfer T, Hinterreiter J, Reiss MA, Möstl C, Davies JA, Bailey RL, et al. Evaluation of Cme Arrival Prediction Using Ensemble Modeling Based on Heliospheric Imaging Observations. Space Weather (2021) 19:e2020SW002553. doi:10.1029/2020sw002553
15. Owens M, Lang M, Barnard L, Riley P, Ben-Nun M, Scott CJ, et al. A Computationally Efficient, Time-dependent Model of the Solar Wind for Use as a Surrogate to Three-Dimensional Numerical Magnetohydrodynamic Simulations. Solar Phys (2020) 295:1–17. doi:10.1007/s11207-020-01605-3
16. Nolte JT, Roelof EC. Large-scale Structure of the Interplanetary Medium. Sol Phys (1973) 33:241–57. doi:10.1007/bf00152395
17. Riley P, Lionello R, Mikić Z, Linker J. Using Global Simulations to Relate the Three‐Part Structure of Coronal Mass Ejections to In Situ Signatures. Astrophysical J (2008) 672:1221–7. doi:10.1086/523893
18. Riley P, Stevens M, Linker JA, Lionello R, Mikic Z, Luhmann JG. Modeling the Global Structure of the Heliosphere during the Recent Solar Minimum: Model Improvements and Unipolar Streamers. AIP Conf Proc (2012) 1436:337–43. doi:10.1063/1.4723628
19. Riley P. CME Dynamics in a Structured Solar Wind. In: SR Habbal, R Esser, V Hollweg, and PA Isenberg, editors. Solar Wind Nine, Am. Inst. Phys. Conf. Proc. (1999). p. 131
20. Rudy SH, Brunton SL, Proctor JL, Kutz JN. Data-driven Discovery of Partial Differential Equations. Sci Adv (2017) 3:e1602614. doi:10.1126/sciadv.1602614
21. Riley P, Luhmann J, Opitz A, Linker JA, Mikic Z. Interpretation of the Cross-Correlation Function of ACE and STEREO Solar Wind Velocities Using a Global MHD Model. J Geophys Res (2010) 115:11104–+. doi:10.1029/2010JA015717
Keywords: heliosphere (711), solar wind streams, coronal mass ejection, magnetohydrodynamics, space weather, in situ measurements
Citation: Riley P and Issan O (2021) Using a Heliospheric Upwinding eXtrapolation Technique to Magnetically Connect Different Regions of the Heliosphere. Front. Phys. 9:679497. doi: 10.3389/fphy.2021.679497
Received: 11 March 2021; Accepted: 28 April 2021;
Published: 19 May 2021.
Edited by:
Emiliya Yordanova, Institute for Space Physics (Uppsala), SwedenReviewed by:
Chaowei Jiang, Harbin Institute of Technology, Shenzhen, ChinaJulia E. Stawarz, Imperial College London, United Kingdom
Copyright © 2021 Riley and Issan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Pete Riley, pete@predsci.com