Solar Wind and its Evolution
By using our previous results of magnetohydrodynamical simulations for the solar wind from open flux tubes, I discuss how the solar wind in the past is different from the current solar wind. The simulations are performed in fixed one-dimensional super-radially open magnetic flux tubes by inputing various types of fluctuations from the photosphere, which automatically determines solar wind properties in a forward manner. The three important parameters which determine physical properties of the solar wind are surface fluctuation, magnetic field strengths, and the configuration of magnetic flux tubes. Adjusting these parameters to the sun at earlier times in a qualitative sense, I infer that the quasi-steady-state component of the solar wind in the past was denser and slightly slower if the effect of the magneto-centrifugal force is not significant. I also discuss effects of magneto-centrifugal force and roles of coronal mass ejections.
100 \publishyear2011 \frompageSolar Wind and its Evolution \topageReferences \copyrighttext
Young solar-type stars are generally very active: X-ray flux is up to times larger than the present-day sun (Güdel et al., 1997; Güdel, 2004), and the X-ray temperature is also higher (Ribas et al., 2005; Telleschi et al., 2005). Observations of young main sequence stars show very strong magnetic field strengths with an order of kG or even larger (Donati & Collier Cameron, 1997; Saar & Brandenburg, 1999; Saar, 2001, see also Donati & Landstreet 2009 for recent review), which are much stronger than the average strength of 1-10 G of the current sun. Young stars are generally fast rotators and strong magnetic fields are generated by strong differential rotation through dynamo activities of magnetoconvection (e.g., Brun et al., 2004). As time goes on, a star loses its angular momentum through magnetic stellar wind (Weber & Davis, 1967) and the magnetic activities become weak (e.g., Ayres, 1997).
Solar wind is hot plasma emanating from the Sun, and the mass loss
rate amounts to g s
(yr) at present.
As inferred from high activity of young solar-type stars,
I expect that the solar wind was stronger at earlier times.
Wood et al. (2002, 2005) observed asterospheres of low-mass stars by the
Hubble space telescope and determined the mass loss rates.
The estimated mass loss rates show a decreasing trend with time,
except for very young stars, thought there is still a large scatter in the
data, partly because the observed stars range from very low-mass to near
The solar wind roughly consists of two components. The first component is called fast solar wind with speed of 700 – 800 km s at the earth orbit. Fast solar wind is more steady and streams out from coronal holes which generally corresponds to open magnetic flux tube regions. The other component is slow solar wind with speed km s. Slow wind is more complex and transient, and in most cases comes from lower latitude regions. Recent HINODE observations show that some portions of the slow wind appear to be originating from open flux tubes near active regions (Sakao et al., 2007; Imada et al., 2007; Harra et al., 2008), while there are still debates of other sources for slow wind.
As other types of magnetized plasma from the sun, coronal mass ejections (CMEs) also are supposed to affect atmospheres of planets especially at early times. In this paper, however, I mainly focus on solar wind from open magnetic flux tubes, and briefly mention effects of CMEs afterward. Our group has carried out forward numerical modeling of the solar wind by using direct magnetohydrodynamical (MHD) simulations, in which we can directly test the response of the solar atmosphere to the surface fluctuations and properties of magnetic fields. By using the results of the simulations with parameters suitable for the younger sun, I discuss the solar wind evolution.
2 Simulation Model
In this section, I briefly describe the simulations of the solar wind. For detail, please refer to Suzuki & Inutsuka (2006, SI06, hereafter). In order to cover the region with huge density contrast from g cm at the photosphere to g cm at the outer simulation boundary located at astronomical unit (AU), simple one dimensional (1-D) open flux tubes are adopted. The effects of super-radial expansion of flux tubes are incorporated by taking into account an expansion factor in the conservation of magnetic flux of radial () component, :
where is a super-radial expansion factor. Most of the solar surface is covered by closed magnetic loop structure. Then, open flux tubes rapidly open above these loops; is introduced to consider this effect. The same function as in (Kopp & Holzer, 1976) is adopted for (see SI06 for detail).
In SI 06 we input the transverse fluctuations of the field line by the granulations at the photosphere. Amplitude, , at the photosphere is chosen to be compatible with the observed photospheric velocity amplitude km s (Holweger et al., 1978). SI06 tested various types of spectra; in this paper I discuss the results of the power spectra in proportion to , where is frequency. The surface fluctuations generate upgoing Alfvén waves which contribute to the acceleration of the solar wind in upper regions. At the outer boundaries, a non-reflecting condition is imposed for all the MHD waves, which enables us to carry out simulations for a long time until quasi-steady state solutions are obtained without unphysical wave reflection.
SI06 dynamically treat the propagation and dissipation of the waves and the heating and acceleration of the plasma by solving ideal MHD equations with the relevant physical processes, including the sun’s gravity, radiative cooling, and thermal conduction (Suzuki & Inutsuka, 2005, SI06). I do not take into account stellar rotation, which may be important for very young stars. For the initial condition, we assume a static atmosphere with temperature K in order to see whether the atmosphere is heated up to coronal temperature and accelerated to accomplish the transonic flow. From , the transverse fluctuations are injected from the photosphere and continue the simulations until the quasi-steady states are achieved.
3 Solar Wind Response to Surface Condition
By injecting fluctuations from the photospheric surface, the initially cool and static atmosphere is effectively heated and accelerated by the dissipation of the generated upgoing Alfvén waves. The sharp transition region which divides the cool chromosphere with K and the hot corona with K is formed owing to a thermally unstable region around K in the radiative cooling function (Landini & Monsignori-Fossi, 1990). The hot corona streams out as the transonic solar wind (see Suzuki & Inutsuka, 2005, & SI06 for detail).
The heating and acceleration of the solar wind plasma in inner heliosphere results from the dissipation of Alfvén waves. In the simple 1-D treatment, Alfvén waves mainly dissipate via nonlinear mode conversion; slow MHD waves are nonlinearly generated from outgoing Alfvén waves and the slow MHD waves are damped by shock formation as a result of steepening of wave shape. The shocks also heat up surrounding plasma, which play a central role in the heating of the solar wind plasma. Magnetic pressure associated with the Alfvén waves decreases with height as a result of the successive dissipation process, which directly pushes the plasma outward (momentum input) in addition to gas pressure.
The young sun is more active than the present-day sun, and therefore, I expect that surface fluctuations are stronger. The properties of the magnetic fields are also supposed to be very different. I discuss how properties of the solar wind are affected by fluctuation amplitudes and magnetic fields at the solar surface.
3.1 Dependence on Surface Fluctuation Amplitude
Figure 1 shows the response of the solar wind plasma to the surface fluctuations with different amplitudes, km s, 0.7km s and 0.4km s. The intermediate case (km s) is a reference case which explains the present-day solar wind.
The maximum temperature of km s case is K (Table 1), which is cooler than the usual corona. The density is lower than the fiducial case by 1-2 orders of magnitude because the sufficient mass cannot be supplied to the corona by the chromospheric evaporation owing to the low temperature; the evaporation is drastically suppressed as decreases, since the conductive flux () sensitively depends on . As a result, the mass flux () becomes much lower, by a factor of 50, than that of the present solar wind.
On the other hand, the larger km s case gives larger coronal density, although the coronal temperature is only slightly higher than the temperature of the reference case. As a result of the dense corona, the mass flux of the solar wind is 10 times larger (Table 1). This case shows complicated structure of temperature and density near the solar surface. The temperature starts to increase from a deeper location around than the other cases. Thanks to this, the decrease of the density is slower (larger pressure scale height) so that the density around is two orders of magnitude larger than that of the reference case with km s. However, the temperature decreases slightly instead of showing a monotonical increase; it cannot go over the peak of the radiative cooling function at K (Landini & Monsignori-Fossi, 1990) because the radiative loss is efficient, owing to the large density. The temperature again increases from and above there the corona forms. The coronal density is larger than the that of the reference case as explained above.
A striking feature is that the mass flux of solar winds sensitively depends on the input wave amplitude at the solar photosphere. Since the injected energy is proportional to the square of , the difference of the input energy between the small amplitude case (km s) and the large amplitude case (km s) is 12. However, the output mass flux of the large amplitude case is 400 times larger than the mass flux of the small amplitude case.
The sensitive behavior on the surface amplitude can be understood by reflection and nonlinear dissipation of Alfvén waves (SI06). Because the heating is smaller, the temperature is lower in the smaller cases. Then, the scale height becomes smaller and the density decreases rapidly. The Alfvén speed, , changes more rapidly due to the rapid decrease of density and the wave shape is largely deformed, which enhances the reflection. When the input wave energy decreases, a positive feedback operates; a smaller fraction of the energy can reach the coronal region. As a result, the density and mass flux of the solar wind becomes much smaller than the decreasing factor of the input energy.
In addition to the effect of the wave reflection, the nonlinear dissipation of the Alfvén waves also plays a role in the sensitive dependence on the input wave energy (SI06). When the input wave energy becomes smaller, the density becomes smaller as explained above. Then, the nonlinearity of Alfvén wave, , decreases, not only because the amplitude, , is smaller but also because the Alfvén speed, , is larger. Therefore, the Alfvén waves do not dissipate and the heating is reduced, which further decreases the density; this is another type of positive feedback, which also results in the sensitive dependence of the density on the input wave energy.
3.2 Magnetic Fields
At earlier epochs the sun is expected to have possessed stronger magnetic field than today, as inferred from observations of young stars (Donati & Collier Cameron, 1997; Saar & Brandenburg, 1999) The configuration of magnetic fields as well as field strength plays an essential role in determining the physical solar wind conditions. In the 1-D treatment as adopted in this paper, a super-radial expansion factor, , determines the geometry of flux tubes. Wang & Sheeley (1990, 1991) showed that the solar wind speed at 1AU is anti-correlated with from their long-term observations. This can be naively understood by energetics consideration; since directly determines an adiabatic loss of plasma in flux tubes, solar wind plasma in largely expanding (larger ) flux tubes undergoes larger energy loss.
By the comparison of the outflow speed obtained by interplanetary scintillation measurements with observed photospheric field strength Kojima et al. (2005) have found that the solar wind velocity is better correlated with surface magnetic field strength, divided by the expansion factor, : They claim that is a better index for the speed of solar wind than individual or . The correlation between and solar wind speed is well explained by Alfvén waves in expanding flux tubes (Suzuki, 2006). Because the energy flux of input Alfvén waves is proportional to , the positive correlation of solar wind speed with is quite natural. The negative correlation with reflects the adiabatic loss of Alfvén waves in expanding tubes as explained above.
One can understood the correlation more specifically from nonlinear dissipation of Alfvén waves as well. The nonlinearity of the Alfvén waves, is controlled by in the outer region where the flux tube is already super-radially open. Wave energy does not effectively dissipate in the larger case in the subsonic region because of relatively small nonlinearity, and more energy remains in the supersonic region. In general, energy and momentum inputs in the supersonic region result in higher wind speed, while those in the subsonic region raise the mass flux () of the wind by increasing density (Lamers & Cassinelli, 1999). This indicates that the solar wind speed is positively correlated with .
4 Evolution of Solar Wind
I have described based on our numerical simulations that the solar wind from open magnetic flux tubes is mainly controlled by the three parameters: surface fluctuation amplitude, , radial magnetic field strength at the photosphere, , and super-radial expansion factor, , of flux tubes. Thus, if I can determine these parameters for the young sun, I can estimate the properties of solar wind at that time.
4.1 Surface Fluctuation Amplitude
I can speculate that was larger at earlier times from circumstantial evidences. First, the rotation of the sun should be faster at earlier times because of loss of angular momentum with time. Hence, in addition to surface convection, interior differential rotation will lead to stronger surface fluctuations. From observational facts, younger stars shows higher X-ray activities as discussed in the introduction section, which also anticipates larger . An important issue is that a small increase of leads to a large increase of mass flux of solar wind.
4.2 Magnetic Field ––
is also supposed to be larger than today as implied from observations of young stars (Donati & Collier Cameron, 1997; Saar & Brandenburg, 1999). Rossby number, , which is defined as the ratio of stellar rotation period to convective turn-over time, is an important index which describes stellar magnetic activities (Noyes et al., 1984). Surface magnetic field strength, , multiplied by a magnetic areal filling factor, , is well-correlated with . Stars with small , which correspond to fast rotating young stars, have 1-10 kG (Saar, 2001). In some stars ’s are obtained and the values are 50-70%, which shows that is an order of 10 kG in these stars.
If all the magnetic flux is open to the interplanetary space, I can set
; in reality, however, closed loop structure is more dominant,
and then, it is expected that .
I speculate that super-radial expansion, , of the young sun was
larger as discussed from now, although it is not simple to pin down
the specific value.
Recently, the configuration of magnetic field has been observed
in a number of low-mass stars (e.g. Donati et al., 2008).
Fast rotating solar mass stars possess non-axisymmetric poloidal components
with substantial toroidal fields (see Fig.3 of Donati & Landstreet, 2009).
This implies that the field configuration is quite complicated, and is very
different from ordered dipole structure.
In comparison with 11-year periodicities of present-day solar activity
During the solar minimum periods, the solar wind consists of fast wind which is from mid- to high-latitude regions, and slow wind from lower-latitude regions. On the other hand, during the solar maximum periods, most of the regions are occupied by transient slow wind (McComas et al., 2008). At the solar maximum magnetic fields are dominated by a dipole component, while at the solar minimum higher multi-pole moments become more important showing more complex field configuration (Hakamada et al., 2005). Because the surface is mostly covered with closed magnetic fields during the solar maximum, the field strength in the outer regions is not as strong as that at the solar minimum although the field strength at the photosphere is stronger. This fact shows that is larger during the solar maximum. The speed of the solar wind during the solar maximum is slower, which is consistent with the observational trend (Kojima et al., 2005), because of larger adiabatic loss in open magnetic flux tubes.
Based on these considerations, namely (i) the magnetic field structure of the young sun more resembles the current solar maximum than the solar minimum, and (ii) open flux tubes show larger expansion during the solar maximum because larger surface area is occupied by closed structure, I infer that open magnetic flux tubes of the young sun have much larger . I speculate that, even though the magnetic field strength at the footpoints is stronger, , is smaller at earlier time.
4.3 Solar Wind in the Past
I would like to discuss how the properties of solar wind change in the past based on the simulation results. I have stated that was larger and was smaller at earlier times. These conditions imply considerably denser but slightly slower solar wind in the past, although it is quite difficult to quantitatively determine the physical condition.
In this paper I neglect the effect of magneto-centrifugal force in driving winds (Weber & Davis, 1967; Mestel, 1968). Under the typical solar condition, if the rotation period is 4 days (6-7 times faster than the present sun) or less, magneto-centrifugal force plays an important role and the wind speed becomes significantly higher (Newkirk, 1980).
Recently, Holzwarth & Jardine (2007) calculated the evolution of solar wind by considering megneto-centrifugal force. They concluded that the solar wind in the past was faster but not so much denser than today, which are different from my conclusion. Holzwarth & Jardine (2007) adopted the density at a ’reference point’, 1.1 times of the stellar radius, in a parameterized way; they assume power-law dependence on stellar rotation rate. Since the basal density strongly constrains the mass loss rate (see, e.g. section 3 of Lamers & Cassinelli, 1999) and the adopted dependence in their standard case is rather weak, the derived mass loss rate shows weak dependence with time. In a strict sense, the density at the coronal base, accordingly the density at the reference level, is determined by the energy balance among heating, thermal conduction, and radiative cooling (SI06; see also Suzuki 2007 for different types of stars). The base densities of young active stars are larger than those of quiescent stars, and the mass loss rates are supposed to be significantly larger. In this case the wind speed becomes slower because more mass needs to be pushed away. The same tendency on the coronal base density is also reported by Stereborg et al. (2011) who performed 3D MHD simulations for a young sun.
It is not still possible to discuss the wind speed in a quantitative sense at the moment. My conclusion of the slightly slower wind in the past might be corrected because magneto-centrifugal force could be significant. On the other hand, I suppose that the wind speed is not so high as estimated by Holzwarth & Jardine (2007) bacause the coronal base density should larger at that time than their assumption. Accordingly, the density of the solar wind could be larger than those calculated in Holzwarth & Jardine (2007).
In this paper I focus on quasi-steady-state solar wind from open magnetic flux regions, and do not consider CMEs. While the quasi-steady-state component dominantly contributes to the total mass loss from the present-day sun than CMEs, at early times CMEs might play a more important role because of stronger magnetic activities. From the fossil record on the lunar surface, it is inferred that the speed of solar wind was faster in the past (Ray & Heymann, 1980; Newkirk, 1980). I infer that the fossil record comprises effects of both quasi-steady state component and CMEs. If strong CMEs dominantly contribute to the mass loss, the bulk speed of plasma flow from the sun would be recognized as higher. I speculate that the fossil inprint on the moon is connected with high CME activities at early epochs.
Based on our previous results of MHD simulations, I discuss the evolution of the quasi-steady-state component of solar wind from open flux tubes. The properties of the solar wind are determined by the three parameters, surface fluctuation amplitude, magnetic field strength, and expansion of open flux tube. Referring to observational data of young low-mass stars, I suppose that the surface amplitude was larger and the ratio of field strength to flux tube expansion was smaller. Following this speculation, the solar wind in the past was dense and slightly slower than today, while if the effect of magneto-centrifugal force is taken into accout, the speed of the solar wind might be higher at very early epoch than my conclusion. CMEs are also supposed play a more important role in the past than today because of higher magnetic activities. In order to understand the effects of solar magnetized plasma on the atmospheres of planets, we need to understand the evolution of CMEs as well as the evoltuion of quasi-steady-state solar wind.
Acknowledgements.acknowledgment This work was supported in part by Grants-in-Aid for Scientific Research from the MEXT of Japan, 19015004, 20740100, 22864006, and Inamori Foundation. The author thanks two anonymous reviewers for many constructive suggestions to improve the paper. Guest editor M. Yamauchi thanks two anonymous reviewers in evaluating this paper. (Received February 3, 2011; Revised April 6, 2011; Accepted April 11, 2011).
- Although the observed stars include not only G-type (solar-analogue) stars but also lower mass K- and M-type stars, it is worth discussing trends because the mechanism of mass loss of these stars are supposed to be similar; the stellar wind is accelerated by wave and turbulent motions associated with magnetic field driven by fluctuation in a surface convective zone.
- I would like to note that the 11-year duration is modulated with time depending activity even in last several hundred years (Miyahara et al., 2008)
- Ayres, T. R., Evolution of the solar ionizing flux, J. Geophys. Res., 102, 1641 - 1652, 1997
- Brun, A., S., Miesch, M. S., Toomre, J., Global-Scale Turbulent Convection and Magnetic Dynamo Action in the Solar Envelope, Astrophys. J., 614, 1073 - 1098, 2004
- Donati, J.-F. & Collier Cameron, A., Differential rotation and magnetic polarity patterns on AB Doradus, Mon. Not. Roy. Astron. Soc., 291, 1 - 19, 1997
- Donati, J.-F. et al., Large-scale magnetic topologies of early M dwarfs, Mon. Not. Roy. Astron. Soc., 390, 545 - 560
- Donati, J.-F. & Landstreet, J. D., Magnetic Field of Nondegenerate Stars, Ann. Rev. Astron. Astrophys., 47, 333 - 370
- Güdel, M. X-ray astronomy of stellar coronae, Astron. Astrophys. Rev., 12, 71 - 237, 2004
- Güdel, M., Guinan, E. F., Skinner, S. L., The X-Ray Sun in Time: A Study of the Long-Term Evolution of Coronae of Solar-Type Stars, Astrophys. J., 483, 947 - 960, 1997
- Hakamada, K., Kojima, M., Ohmi, T., Tokumaru, M, Fujiki, K., Correlation between Expansion Rate of the Coronal Magnetic Field and Solar Wind Speed in a Solar Activity Cycle Sol. Phys., 227, 387 - 399, 2005
- Harra, L. K. et al. Outflows at the Edges of Active Regions: Contribution to Solar Wind Formation?, Astrophys. J. Lett., 676, L147 - L150, 2008
- Holweger, H., Gehlsen, M., & Ruland, F.: Spatially-averaged properties of the photospheric velocity field, Astron. Astrophys., 70, 537 - 542, 1978
- Holzwarth, V. & Jardine, M., Theoretical Mass Loss Rates of Cool Main-sequence Stars, Astron. Astrophys., 463, 11 - 21, 2007
- Imada, S. et al., Discovery of a Temperature-Dependent Upflow in the Plage Region During a Gradual Phase of the X-Class Flare., Pub. Astron. Soc. Japan, 59, S793 - 799, 2007
- Johns-Jrull, C. M., The Magnetic Fields of Classical T Tauri Stars, Astrphys. J., 664, 975 - 985, 2007
- Kojima, M., K. Fujiki, M. Hirano, M. Tokumaru, T. Ohmi, and K. Hakamada, Solar Wind properties from IPS observations, The Sun and the heliosphere as an Integrated System, Giannina Poletto and Steven T. Suess, Eds. Kluwer Academic Publishers, 147 - 181, 2005
- Kopp, R. A. & Holzer, T. E., Dynamics of coronal hole regions. I – Steady polytropic flows with multiple critical points, Sol. Phys., 49, 43 - 56, 1976
- Lamers, H. J. G. L. M. & Cassinelli, J. P., Introduction to Stellar Wind, Cambridge, 1999
- Landini, M. & Monsignori-Fossi, B. C.: The X-UV spectrum of thin plasmas, Astron. Astrophys. Supp., 82, 229 - 260, 1990
- McComas, Ebert, R. W., Elliot, H. A., Golestein, B. E., Cosling, J. T., Schwadron, N. A., Skoug, R. M., Weaker solar wind from the polar coronal holes and the whole Sun., Geophys. Res. Lett., 35, L18103
- Mestel, L., Magnetic Braking by a Stellar Wind - I., Mon. Not. Roy. Astron. Soc., 138, 359 - 391
- Miyahara, H., Yokoyama, Y., Masuda, K., Possible link between multi-decadal climate cycles and periodic reversals of solar magnetic field polarity, Earth. Planet. Sci. Lett., 272, 290 - 295
- Newkirk, Gordon Jr., Solar Variability on Time Scale of 10 years to 10 years, Proc. Conf. Ancient Sun, 293 - 320, 1980
- Noyes, R. W., Hartmann, L. W., Baliunas, S. L., Duncan, D. K., & Vaughan, A. H., Rotation, Convection, and Magnetic Activity in Main-sequence Stars, Astrophys. J., 279, 763 - 777, 1984
- Ray, J. & Heymann, D., A Model for Nitrogen Isotopic Variations in the Lunar Regolith: Possible Solar System Contributions from a Nearby Planetary Nebula, Proc. Conf. Ancient Sun, 491 - 512, 1980
- Ribas, I., Guinan, E. F., Güdel, M., Audard, M., Evolution of the Solar Activity over Time and Effects on Planetary Atmospheres. I. High-Energy Irradiances (1-1700øA)., Astrophys. J., 622, 680 - 694, 2005
- Saar, S. H., Recent Measurements of (and Inferences about) Magnetic Fields on K and M stars, the 11th Cool Stars, Stellar Systems and the Sun, ASP Conf. Series, 223, 292 - 299, 2001
- Saar, S. H. & Brandenburg, A., The Evolution of the Magnetic Activity Cycle Period. II. Results for an Expanded Stellar Sample, Astrophys. J., 524, 295 - 310, 1999
- Sakao, T. et al., Continuous Plasma Outflows from the Edge of a Solar Active Region as a Possible Source of Solar Wind., Science, 318, 1585 - 1588, 2007
- Stereborg, M. G., Cohen, O., Drake, J. J., Gombosi, T. I., Modeling the Young Sun’s Solar Wind and Its Interaction with Earth’s Paleomagnetosphere, J. Geophys. Res., 116, A01217, 2011
- Suzuki, T. K., Forcasting Solar Wind Speeds, Astrophys. J. Lett., 640, L75 - L78, 2006
- Suzuki, T. K., Structured Red Giant Winds with Magnetized Hot Bubbles and the Corona/Cool Wind Dividing Line, Astrophys. J., 659, 1592 - 1610
- Suzuki, T. K. & Inutsuka, S.: Making the corona and the fast solar wind: a self-consistent simulation for the low-frequency Alfvén waves from photosphere to 0.3AU, Astrophys. J. Lett., 632, L49 - L52, 2005
- Suzuki,T. K. & Inutsuka, S.:Solar Winds Driven by Nonlinear Low-Frequency Alfvén Waves from the Photosphere : Parametric Study for Fast/Slow Winds and Disappearance of Solar Winds, J. Geophys. Res., 111, A6, A06101, 2006 (SI06)
- Telleschi, A. Güdel, M., Briggs, K., Audard, M., Ness, J.-U., Skinner, S. L., Coronal Evolution of the Sun in Time: High-Resolution X-Ray Spectroscopy of Solar Analogs with Different Ages, Astrophys. J., 622, 653 - 679, 2005
- Wang, Y.-M. & Sheeley, Jr, N. R., Solar wind speed and coronal flux-tube expansion, Astrophys. J., 355, 726 - 732, 1990
- Wang, Y.-M. & Sheeley, Jr, N. R., Why fast solar wind originates from slowly expanding coronal flux tubes, Astrophys. J. Lett., 372, L45 - L48, 1991
- Weber, E. J., Davis, L.,Jr. The Angular Momentum of the Solar Wind, Astorphys. J., 148, 217 - 227, 1967
- Wood, B. E., Müller, H.-R., Zank, G. P., Linsky, J. L., Measured Mass-Loss Rates of Solar-like Stars as a Function of Age and Activity, Astorphys. J., 574, 412 - 425, 2002
- Wood, B. E., Müller, H.-R., Zank, G. P., Linsky, J. L., Redfield, S., New Mass-Loss Measurements from Astrospheric Ly Absorption, Astrophys. J. Lett., 628, L143 - L146, 2005 (Author), (article’s title), (Journal), (Volume), (page number)–(page number), (year).