A correlation between solar wind speed at Earth and the amount of magnetic field line expansion in the corona was verified in 1989 using 22 years of solar and interplanetary observations. We trace the evolution of this relationship from its birth 15 years earlier in the Skylab era to its current use as a space weather forecasting technique. This paper is the transcript of an invited talk at the joint session of the Historical Astronomy Division and the Solar Physics Division of the American Astronomical Society during its 224th meeting in Boston, MA, on 3 June 2014.
This is the 25th anniversary of the origin of the Wang–Sheeley–Arge (WSA)
model. It occurred in 1989 when there was no internet as we now know it, no
cell phones, and no Google, and Apple stock was selling for
What is the WSA model? It is a procedure for deriving solar wind speed from maps of the photospheric magnetic field. It is based on the idea that, far from the Sun, the speed of the solar wind depends on the path that the wind took as it passed through the lower corona, as illustrated in Fig. 1. If the flux tube that guided the flow flares out, then the distant speed is slow. But if the flux tube remains focused, then the distant speed is fast. The amount of expansion is obtained from a current-free extension of the photospheric field.
A cartoon, illustrating that fast wind comes from magnetic flux tubes with small coronal expansion and that slow wind comes from flux tubes with large coronal expansion.
This road through 1989 had four principal milestones. The first milestone occurred on 1 June (25 years prior to the start of this meeting) when I was defending the Office of Naval Research (ONR) Contract Research Program in Solar Physics before an external review committee. I linked the ONR work to our own Naval Research Laboratory (NRL) research by saying that we used the contractor magnetograms to derive the locations of open field regions and thus coronal holes. And we used our new flux-transport program, created at NRL by Jay Boris and Rick DeVore, to evolve the magnetic maps forward in time by 1 or 2 months so that we could calculate the locations of future coronal holes. Because coronal holes were sources of high-speed solar wind streams and associated geomagnetic disturbances, I argued that this procedure ought to be useful for forecasting geomagnetic activity.
The committee liked the leveraging to the NRL program, but it wanted the procedure to be available at the US Air Force (USAF) and National Oceanic and Atmospheric Administration (NOAA) forecasting centers. Consequently, an action item was created, asking me to devise a plan to make our coronal-hole forecasting procedure operational at those centers.
The second, even more important milestone occurred during August as
Yi-Ming Wang was trying to find a magnetic-field quantity that was
correlated with the solar wind speed. Previously, we had some success
comparing coronal holes and wind speed in Bartels displays
The situation changed when Yi-Ming read an old paper by
The original Bartels comparison of simulated and observed wind speed
(left columns) and interplanetary magnetic field (right columns) from the paper by
The
High Altitude Observatory image of the 12 November 1966 solar eclipse.
As Levine et al. knew and Yi-Ming recognized, the important quantity
was not the expansion of the entire coronal hole, but the expansion
of the individual flux tubes within the hole. As shown in Fig. 4,
the individual flux tubes of an axisymmetric dipole field have
relatively small expansions, except toward the edge of the
hole where the tubes flare out considerably. Because this equatorial
streamer-belt/current-sheet region was known to be associated with slow
wind, it followed that the solar wind speed was inversely related to
the amount of expansion that occurs in individual flux tubes of the
coronal hole. However, the question of how much slow wind comes from
rapidly expanding flux tubes at the edge of a coronal hole and how much
comes from the tops of the stretched coronal loops in the streamer belt
remains an active subject of study
Current-free field lines of a photospheric dipole field, constrained to be radial at 2.5 solar radii to simulate the stretching effect of the coronal expansion. These field lines show the increased flux-tube expansion toward the edge of the polar coronal holes.
Yi-Ming's discovery was the break that we needed. When I told my ONR contact, Ken Davis, about the new result, he asked me to brief Bruce Robinson, the ONR Director of Science, in a meeting on 24 August. As you can see from my briefing notes in Fig. 5, the words “new result last week” pinpoint the discovery to a date around 17 August 1989. In addition, the notes contain three points that summarized our status at that time: “seem to have found a way to derive wind speed from magnetograms”, “a missing link in our forecasting capability”, and “relevant to the request to make our model operational”. This briefing led to the third milestone 4 days later on 28 August 1989, which was to give ONR the requested plan to make our forecasting model operational at the USAF and NOAA forecasting centers. The expansion factor–wind speed relation was a central part of that plan.
The fourth milestone was to publish the result and make it familiar to the
scientific community. As shown in Fig. 6, our paper,
Notes prepared for my 24 August 1989 meeting with the ONR Science Director. They pinpoint Yi-Ming's discovery around 17 August 1989.
The title page of our paper relating solar wind speed to the coronal expansion factor as published in the 1 June 1990 issue of The Astrophysical Journal.
During the next year, we pursued the implications of the inverse correlation
between wind speed and coronal expansion factor
We published this result on 1 May 1991
September 1991 photo, showing (left to right) Yi-Ming Wang, the author, Ken Dere, Ana Nash, and Ken Widing in Goslar, Germany, during Solar Wind Seven. The photo was taken by Don Michels, who was coordinating international efforts for the Large Angle Spectrometric Coronagraph (LASCO). LASCO was launched as part of the Solar and Heliospheric Observatory (SOHO) 4 years later on 2 December 1995, and is still operating successfully.
The web site created by Nick Arge, showing a running 3-day forecast of solar wind speed and interplanetary magnetic field (IMF) polarity derived from the Wang–Sheeley–Arge model during a Carrington rotation in April–May 2014.
The title page of the announcement by
The National Centers for Environmental Prediction web site (
The original plan of making our forecasting procedure operational was eventually carried out. During 1993–1995, we sent monthly updates of our calculations to Capt. Paul Gehred at the USAF/AWS in Boulder, CO, who compared geomagnetic forecasts based on those calculations with the recurrence-based forecasts that were being used at the operational centers at that time. I lost contact with Paul in mid-November 1995, when he moved to the Air Force Research Lab in Dayton, OH, but his reports indicated that he had obtained slightly better results with the model-based forecasts than with the recurrence-based ones.
The process entered a new phase around the time of the AGU meeting in December 1995, probably in an encounter between Bob McCoy (who had been managing the ONR program since I left it in 1994), Ernie Hildner (the NOAA/Space Environment Laboratory (SEL) Director), and Todd Hoeksema (a member of the Wilcox Solar Observatory (WSO) contractor team). I think that Bob came up with the idea to reduce the funding of each contractor team by a little bit and use the savings to fund a post doc at NOAA/SEL to make the Wang–Sheeley (WS) model operational. It was really a brilliant idea. The plan would not only provide additional work on the model, but it would also save the dwindling ONR solar program for several more years. For Ernie, it would bring additional funding to NOAA/SEL, and for Todd and the other contractors, it would mean the continuation of most of their ONR support. So when NOAA/SEL scientist Vic Pizzo faxed me a copy of a 30 January 1996 NOAA proposal to fund a post doc for upgrading the Wang–Sheeley solar wind prediction model, I strongly recommended it to Bob McCoy.
Nick Arge became the new post doc and joined NOAA/SEL in early October 1996.
When he visited NRL in late November, we taught him to use our codes and gave
him the software to take back to Boulder. Sometime in 1997, he created the WS
web site with the ONR and NOAA logos shown in Fig. 8. During 1997–2003, Nick
continued to make essential contributions to the development, implementation,
and improvement of the model
In 2004, Nick suggested that the coronal portion of the WSA model be coupled with the 3-D MHD solar wind model that Dusan Odstrcil had developed for calculating the outward evolution of the solar wind. The combined effort could be done as part of the Center for Integrated Space Weather Modeling (CISM) effort led by Boston University and supported by the National Science Foundation's Science and Technology Center. The idea was to use the WSA model to derive the wind speed and magnetic field in the vicinity of 20–30 solar radii, where the solar wind is both super-sonic and super-Alfvenic, and then to use this information as the inner boundary condition for the MHD model. Dusan and the CISM leadership team endorsed Nick's idea and supported the effort to combine the two models and develop them further. The MHD model was called ENLIL after the Sumerian god of winds and storms (from EN meaning “lord” and LIL meaning “wind”), so the combined model became known as the WSA-ENLIL model.
Finally, in 2011, Vic Pizzo and his colleagues succeeded in making the WSA-ENLIL
model operational at the NOAA/SEL, by that time called the NOAA Space Weather
Prediction Center (NOAA/SWPC). This transition to operations is described in a paper
in Space Weather
According to Ernie Hildner, the WSA–ENLIL model is the first operational space weather model. Its output is an official National Weather Service (NWS) product and has all the aspects of other NWS products, including reliability, legal status, and resistance to hacker's corruptions, and it is run on the same NWS computers that run the meteorological models. Nick Arge is proud that the Sun now appears on the same page as snow, floods, tornados, and hurricanes, as you can see on the NWS National Centers for Environmental Prediction (NCEP) web site shown in Fig. 10. Bob McCoy, now Director of the Geophysical Institute in Fairbanks, told me that he regards the WSA–ENLIL model as one of our biggest success stories.
The WSA-ENLIL forecasting model as it appeared on the NOAA/NWS web site
in 2014, providing running forecasts of solar wind speed and ion density at Earth
and at the STEREO A and B spacecraft. The link now forwards to a new address given
by
If you click on “Space Weather Prediction Center”, you will find a list of “Products and Data”. In that list, under “Models”, click on “WSA–Enlil Solar Wind Prediction”. You will obtain a running 7-day display of predicted plasma density and wind speed at Earth, at the STEREO A and STEREO B spacecraft, and throughout the inner heliosphere, as shown by the movie in Fig. 11.
(Harry Warren downloaded
the relevant frames from the NWS web site and created the movie used in the
talk. Readers of this article can obtain the corresponding movie by going
to the NWS web address,
In summary, during the past 25 years, we have progressed from a position defending long-term synoptic observations of the Sun's magnetic fields to a forecasting position between floods and tornados on the NWS/NCEP web site. But in this glow of success, let us not forget that there were no web sites in 1989 when the model was created, and that the expansion factor–wind speed relation resulted because researchers like Yi-Ming Wang (and Randy Levine and his colleagues before that) were curious about those solar observations. In this audience, I see people who are using similar observations to study the sunspot cycle and solar activity. It is easy to imagine that their research will ultimately lead to a prominent position on a forecasting medium that has not even been invented yet.
WSA-ENLIL movies mentioned in reference to Fig. 11 can be obtained at
The solar magnetic field data used in the creation of Fig. 2 were obtained from the
Mount Wilson Observatory and may be found at
The author declares that he has no conflict of interest.
I gleaned most of this historical material from my own notes and records, and from my email messages from the 1990s. I am especially grateful to Nick Arge, Bob McCoy, Ernie Hildner, Todd Hoeksema, Vic Pizzo, and Yi-Ming Wang for recent conversations and email messages that helped to fill in some gaps in this historical material. I hope that I have expressed their recollections and opinions accurately. Also, I am grateful to the referees for providing additional help and suggestions. One of the referees, Nick Arge, helped me to describe the formation of the WSA-ENLIL model and the transition to operations during 2004–2011 – activities with which he was directly involved. Financial support was obtained from NASA and CNR. Edited by: T. Kikuchi Reviewed by: N. Arge and one anonymous referee