1 Introduction

It goes without saying that it is generally much easier to study the nearby Sun than it is to study much more distant solar-like stars. Nevertheless, stellar research can address questions about the Sun that observations of the Sun alone cannot answer. The Sun only provides one example of a cool main sequence star, so it cannot tell us by itself how its various properties relate to each other. By observing other solar-like stars, we can see how properties such as stellar activity, rotation, and age are correlated. This can teach us a lot about why the Sun has the properties it does today. It can also tell us what the Sun was like in the past and what it will be like in the future.

Many stellar analogs of solar phenomena are available for study: photospheres, chromospheres, coronae, starspots, magnetic fields, rotation, asteroseismology, etc. (see Skumanich, 1972; Linsky, 1980; Vogt et al., 1987; Gustafsson and Jørgensen, 1994; Johns-Krull and Valenti, 1996; Christensen-Dalsgaard, 2003; Favata and Micela, 2003; Güdel, 2004). Comparing solar properties with those observed for other stars provides a useful context for the solar measurements, improving our understanding of the Sun as well as for stars in general. However, one major solar phenomenon that has proven to be very difficult to study for other stars is the solar wind.

Some types of stellar winds are very easy to detect and study spectroscopically. The massive, radiation-pressure driven winds of hot stars and the cool, massive winds of red giants and super-giants both produce P Cygni emission line profiles that allow the measurement of wind properties with reasonable precision (Harper et al., 1995; Mullan et al., 1998; Kudritzki and Puls, 2000). However, these stars are not solar-like and the winds of these stars are not analogous to the much weaker wind that we see emanating from the Sun. The weak and fully ionized solar wind provides no spectral diagnostics analogous to those used to study more massive stellar winds. Directly detecting a truly solar-like wind around another solar-like star has therefore proven to be a formidable problem.

The first clear detections of winds around other solar-like stars have come from UV spectra of nearby stars from the Hubble Space Telescope (HST). Stellar H I Lyα lines at 1216 Å are always contaminated by very broad, saturated H I absorption. For a long time, this absorption was assumed to be entirely from interstellar H I. However, for some of the nearest stars, the interstellar medium (ISM) cannot account for all of the observed absorption. With the assistance of complex hydrodynamic models of the solar wind/ISM interaction, the excess Lyα absorption has been convincingly identified as being partly due to heated H I gas within our own heliosphere and partly due to analogous H I gas within the “astrospheres” surrounding the observed stars. Note that the word “astrosphere” is used here as the stellar analog for “heliosphere”, although “asterosphere” has also been used in the past (see Schrijver et al., 2003). “Astrosphere” has a longer history, with published usage in the literature dating at least back to 1978 (Fahr, 1978). The term “heliosphere” itself only dates back to the 1960s (Dessler, 1967).

The detection of astrospheric Lyα absorption represents an indirect detection of solar-like stellar winds, since astrospheres do not exist in the absence of a stellar wind. Furthermore, the amount of astrospheric absorption is dependent on the strength of the wind, so the astrospheric absorption has provided the first estimates of mass loss rates for solar-like stars. This article reviews the study of solar-like winds around other stars, especially results using the astrospheric Lyα absorption technique.

  • In Section 2, some background material is provided about the solar wind, ISM, and heliosphere.

  • In Section 3, techniques used to try to directly detect solar-like winds are reviewed.

  • In Section 4, the astrospheric Lyα absorption diagnostic is described in detail.

  • Section 5 provides a review of what the astrospheric analyses have taught us about the solar wind, and discusses some implications of these results both within and outside the realm of solar/stellar physics.

  • Finally, the article ends in Section 6 with some concluding comments about the future of this subject.

2 Background Material

2.1 The solar wind

Before describing how solar-like winds are detected around other stars, it is worthwhile to briefly review what is known about the solar wind and how we study its properties. The solar wind was first detected through its role in the formation of aurorae and the creation of comet tails. As far back as 1896, Kristian Birkeland proposed that aurorae were due to particles emanating from the Sun (see review by Stern, 1989), while Biermann (1951) first described how “corpuscular radiation” from the Sun was responsible for the plasma tails of comets. Today almost everything we know about the solar wind comes from in situ measurements of its properties from satellites. These measurements date back to the Soviet Luna missions in 1959 (Gringauz et al., 1962) and NASA’s Mariner 2 mission in 1962 (Neugebauer and Snyder, 1962). Numerous other spacecraft have participated in studying the solar wind since then, but of particular note are the venerable Voyager 1 and Voyager 2 satellites, which have returned data on the solar wind from 1977 through the present day (see Lazarus and McNutt Jr, 1990). More recently, the Ulysses spacecraft, launched in 1990 and still operating, has provided a first look at the solar wind outside of the ecliptic plane (McComas et al., 2000).

Within the ecliptic plane the solar wind is dominated by low speed streams with typical velocities, proton densities, and temperatures (at 1 AU) of V = 400 km s−1, n(H+) = 5 cm−3, and T = 105 K, respectively, although high speed streams with lower densities and V ≈ 800 km s−1 are not uncommon (see Feldman et al., 1977). At the maximum of the 11-year solar activity cycle, similar solar wind behavior is seen at almost all latitudes. However, Figure 1 shows that at solar minimum the wind above 30° ecliptic latitude is uniformly high speed, low density wind with V ≈ 800 km s−1 (McComas et al., 2000, 2002). This high speed wind originates from coronal holes, which are particularly prominent on the Sun during solar minimum conditions. These solar wind data imply a total mass loss rate for the Sun of ≈ 2 × 10−14 M yr−1. Although thermal temperatures for the wind at 1 AU are of order T = 105 K, densities are low enough that the wind cannot equilibrate to this temperature, and the ionization state of the wind is actually frozen in at coronal temperatures closer to T = 106 K. Hydrogen is fully ionized, meaning that protons are the dominant constituents of the wind by mass.

Figure 1:
figure 1

The solar wind velocity (red/blue line) and density (green line) observed by Ulysses as a function of ecliptic latitude (McComas et al., 2000). During solar minimum conditions, high latitudes are dominated by high speed, low density wind, while low latitudes see mostly lower speed wind with higher densities.

Perhaps the most fundamental question to ask about the solar wind is why it exists. Addressing this question is also necessary to assess whether similar winds should exist around other stars. Even before satellites proved the existence of a more-or-less steady solar wind, Parker (1958) predicted that such a wind should be present due to the existence of the 106 K corona surrounding the Sun. In Parker’s model, the solar wind exists because of thermal expansion from the hot corona. The predictions of this simple model agree remarkably well with the observed properties of the low speed wind that dominates in the ecliptic plane, although additional wind acceleration mechanisms invoking MHD waves have been proposed to explain the high speed streams (see MacGregor and Charbonneau, 1994; Cranmer, 2002). Thus, any star that has a hot corona analogous to that of the Sun should also have a wind analogous to that of the Sun. Observations with X-ray satellites such as Einstein and ROSAT clearly demonstrate that coronae are a ubiquitous phenomenon among cool main sequence stars (see Schmitt, 1997; Hünsch et al., 1999), so solar-like winds should be present around all solar-like stars. However, that does not mean that they are easy to detect (see Section 1).

2.2 The local interstellar medium

The solar wind does not expand forever. Eventually it runs into the local interstellar medium (LISM). The interaction region between the solar wind and LISM is the subject of Section 2.3, and it is through analogous interaction regions around other stars that solar-like stellar winds can be detected (see Section 4), but before the wind/ISM interaction regions are discussed, it is necessary to review the properties of the undisturbed LISM.

The principle method by which one studies the ISM is by observing absorption lines that interstellar material produces in spectra of distant stars. These studies reveal that ISM column densities remain rather low within about 100 pc of the Sun in most directions and then increase dramatically (Sfeir et al., 1999). This low density region is called the Local Bubble. Figure 2 shows a map of the Local Bubble in the Galactic plane (Lallement et al., 2003). The hot plasma within the Bubble has been detected directly from observations of the soft X-ray background (Snowden et al., 1995). Most locations within the Local Bubble are very hot (T ∼ 106 K) and rarified (ne ∼ 10−3 cm−3), and therefore completely ionized.

Figure 2:
figure 2

Map of the Local Bubble in the Galactic plane, where the contours indicate 20 mÅ and 50 mÅ equivalent widths for the Na I D2 line (Lallement et al., 2003). The distance scale is in parsecs.

Although most of the volume of the Local Bubble consists of this hot material, the absorption line studies clearly demonstrate that there are cooler, partially neutral clouds embedded within the Local Bubble. Furthermore, since even the shortest lines of sight show absorption from H I and other low temperature species (Linsky, 1998; Linsky et al., 2000), the Sun must be located within one of these clouds. The cloud immediately surrounding the Sun has been called the Local Interstellar Cloud (LIC), which is roughly 5–7 pc across and has a total mass of about 0.32M (Redfield and Linsky, 2000). There are similar clouds that are apparently adjacent to the LIC (e.g. the “G” cloud and “Hyades” clouds, see Lallement and Bertin, 1992; Redfield and Linsky, 2001), although it is debatable whether the LIC is truly distinct from these clouds. Velocity gradients within a single cloud could in principle create the appearance of multiple clouds in absorption line studies.

The first evidence that the LIC is not entirely ionized came not from absorption line studies but from observations of solar Lyα emission scattering off interstellar H I gas flowing into the heliosphere (Bertaux and Blamont, 1971; Quémerais et al., 1999, 2000). Interstellar atoms have also been observed directly with particle detectors on board spacecraft such as Ulysses (Witte et al., 1993, 1996). Both measurements of LIC material flowing through the heliosphere and LISM absorption line studies have been used to estimate the direction and magnitude of the LIC vector, and the resulting vectors are in good agreement. The heliocentric vector derived from absorption lines has a magnitude of 25.7 km s−1 directed towards Galactic coordinates l = 186.1° and b = −16.4° (Lallement and Bertin, 1992; Lallement et al., 1995).

Other properties of the undisturbed LISM just beyond the heliosphere are less precisely known. Absorption line studies are hampered by probable variations of densities, temperatures, and ionization states within the LIC (see Cheng and Bruhweiler, 1990; Slavin and Frisch, 2002; Wood et al., 2003b), meaning that line-of-sight averages of these properties towards even the nearest stars are potentially different from the actual circumsolar LISM properties. Studies of LISM particles streaming through the heliosphere are hampered by the fact that the properties of these particles are often altered in the outer heliosphere, thereby requiring the assistance of models to extrapolate back to undisturbed LISM conditions (see Izmodenov et al., 2004). In any case, typical temperatures measured for the LIC are T = 6000–8000 K, typical hydrogen densities are n(H I) = 0.1–0.2 cm−3, and typical proton and electron densities are n(H+) ≈ ne = 0.04–0.2 cm−3 (Witte et al., 1993, 1996; Wood and Linsky, 1997, 1998; Izmodenov et al., 1999a; Redfield and Linsky, 2000; Frisch and Slavin, 2003).

2.3 The structure of the heliosphere

Wind/ISM interactions provide the means by which solar-like stellar winds can be detected (see Section 4). Our understanding of these interactions relies heavily on a long history of efforts to model the solar wind/ISM interaction. This heliospheric modeling is summarized briefly here, but for more comprehensive reviews see Holzer (1989); Baranov (1990); Suess (1990) and Zank (1999).

Modeling the large scale structure of the heliosphere began not long after the solar wind’s discovery (Parker, 1961, 1963). The basic structure of the heliosphere, which is shown schematically in Figure 3, is dominated by three prominent boundaries: the termination shock (TS), heliopause (HP), and bow shock (BS). The solar wind is highly supersonic, and the oval-shaped termination shock is where the radial wind is shocked to subsonic speeds. The 26 km s−1 laminar ISM flow is also generally believed to be supersonic, although in principle it could be subsonic if the poorly known ISM magnetic field is strong enough (Zank et al., 1996). Nevertheless, most heliospheric models assume the existence of a bow shock, where the ISM flow is shocked to subsonic speeds (see Figure 3). In between the TS and BS is the heliopause, which is a contact surface separating the plasma flows of the solar and interstellar winds.

Figure 3:
figure 3

Schematic picture of the heliospheric interface from Izmodenov et al. (2002), which can be divided into the 4 regions shown in the figure, with significantly different plasma properties. Region 1: supersonic solar wind; Region 2: subsonic solar wind; Region 3: disturbed interstellar gas and plasma; and Region 4: undisturbed interstellar medium.

The heliospheric structure shown in Figure 3 is inferred almost entirely from hydrodynamic models. However, in 2004, Voyager 1 crossed the TS at a distance of 94 AU from the Sun in roughly the upwind direction of the ISM flow (Stone et al., 2005). Precursors of this crossing were seen years prior to the accepted crossing date (Krimigis et al., 2003; Burlaga et al., 2003; McDonald et al., 2003). At the time of this writing, Voyager 1’s sister satellite Voyager 2 has begun to see these precursors but has not yet officially crossed the TS (Opher et al., 2006). The 94 AU TS distance measured by Voyager 1 is consistent with model predictions (Izmodenov et al., 2003). As for the HP and BS, recent models generally predict upwind distances of ∼ 140 AU and ∼ 240 AU, respectively. The Voyager satellites may not survive long enough to get out this far.

The plasma component of the LISM is diverted around the heliopause due to the strong plasma interactions, but neutrals in the LISM can penetrate into the solar system through the HP and TS. These neutrals were first detected through Lyα backscatter emission (Bertaux and Blamont, 1971). However, even after this discovery most hydrodynamic models of the heliosphere continued to ignore the neutrals since the collisional interactions involving neutrals are much weaker than those involving charged particles. Essentially, the assumption was made that the neutrals would pass through the heliosphere unimpeded, having little or no effect on the heliospheric structure.

It was recognized in the 1970s that the LISM neutrals could in fact play an important role in the solar wind/ISM collision through charge exchange interactions (Holzer, 1972; Wallis, 1975). However, trying to model this is very difficult, because the charge exchange sends the neutral H wildly out of thermal and ionization equilibrium. This means that simple fluid approximations break down and one has to resort to complex multi-fluid codes or ideally fully kinetic codes. It was not until much later that the first codes that treat the plasma and neutrals in a self-consistent manner were first developed (Baranov and Malama, 1993, 1995; Baranov and Zaitsev, 1995; Zank et al., 1996; Izmodenov et al., 1999a; Müller et al., 2000; Izmodenov et al., 2001). These models demonstrate that the heliospheric structure is indeed influenced significantly by the neutrals in many different ways.

For purposes of this article, the importance of the development of heliospheric codes that treat neutrals properly is that it is only because of the existence of neutral hydrogen in the LISM that heliospheric and astrospheric Lyα absorption is detectable, and it is only because of the existence of the self-consistent codes developed to model neutrals in the heliosphere that we can model the astrospheric absorption and extract stellar mass loss rates from the data. Figure 4 shows a heliospheric model that uses a hybrid kinetic code, where the protons are modeled as a fluid but a kinetic code is used for the neutrals (Lipatov et al., 1998; Müller et al., 2000; Wood et al., 2000b). The strong plasma interactions heat and compress LISM protons in between the HP and BS, and thanks to charge exchange processes these high temperatures and densities are transmitted to the neutral H. As a consequence, the heliosphere and astrospheres are permeated by a population of hot hydrogen, which produces a substantial amount of Lyα absorption in HST observations of nearby stars. Most of this absorption comes from the “hydrogen wall” region in between the HP and BS, where densities of the hot H I are particularly high (see Figure 4d). The heliospheric and astrospheric Lyα diagnostic is described in detail in Section 4.

Figure 4:
figure 4

(a) Proton temperature, (b) proton density, (c) neutral hydrogen temperature, and (d) neutral hydrogen density distributions for a heliospheric model from Wood et al. (2000b). The positions of the termination shock (TS), heliopause (HP), and bow shock (BS) are indicated in (a), and streamlines indicating the plasma flow direction are shown in (b). The distance scale is in AU.

3 Direct Wind Detection Techniques

The astrospheric Lyα absorption diagnostic described in detail in Section 4 represents only an indirect detection of stellar winds, since the H I that produces the absorption is essentially LISM rather than wind material. The H I is nevertheless heated by its interaction with the stellar wind, and the astrospheric absorption has therefore proven to be very useful for successfully detecting and measuring solar-like winds around other stars. However, there have been attempts to detect these winds more directly, and these mostly unsuccessful efforts are briefly summarized here.

One can try to look for free-free radio emission from solar-like winds, since they are fully ionized and should therefore produce emission at some level. However, current radio telescopes can only detect winds if they are much stronger than that of the Sun. There have been some claims of very high mass loss rates for a few very active stars using observations at millimeter wavelengths (Mullan et al., 1992), but these interpretations of the data are highly controversial (Lim and White, 1996; van den Oord and Doyle, 1997). The problem is that the coronae of these active stars are also sources of radio emission, which makes problematic the identification of a wind as the source of the emission. Furthermore, it has been argued that massive winds around active stars should absorb the flaring coronal emission that is often observed from these stars, suggesting that massive winds cannot be present (Lim and White, 1996). Nondetections of radio emission have been used to derive upper limits to the mass loss rates of various stars, but these upper limits are typically 2–3 orders of magnitude higher than the solar mass loss rate, so these are not very stringent constraints (Brown et al., 1990; Drake et al., 1993; Lim et al., 1996b; Gaidos et al., 2000).

Variable ultraviolet absorption features observed from the close, eclipsing binary V471 Tau (K2 V+DA) have been interpreted as being due to a wind from the K2 V star (Mullan et al., 1989; Bond et al., 2001). Even if this interpretation is correct, it is questionable whether the wind produced by this star can be considered to be truly “solar-like” given the close presence of the white dwarf companion. In addition, instead of a spherically symmetric wind it has been proposed that the UV absorption could instead be indicative of coronal material being funneled directly from the K2 V star to the white dwarf through the magnetospheric interaction of the two stars (Lim et al., 1996a).

One final wind detection technique that has been proposed is to look for X-ray emission surrounding nearby stars, caused by charge exchange between highly ionized heavy atoms in the stellar wind and inflowing LISM neutrals. This is very analogous to the process by which comets produce X-rays (see Lisse et al., 2001; Cravens, 2002). In the heliosphere, this charge exchange X-ray emission may be responsible for a significant fraction of the observed soft X-ray background (Cravens, 2000). Wargelin and Drake (2002) searched for circumstellar X-ray emission in Chandra observations of the nearest star, Proxima Cen, but they failed to detect any. Based on this nondetection, they quote an upper limit for Proxima Cen’s mass loss rate of < 14. This can be compared with the upper limit of < 350 derived from a nondetection of radio emission from Proxima Cen (Lim et al., 1996b), and the upper limit of < 0.2 derived from the nondetection of astrospheric Lyα absorption (Wood et al., 2001). The astrospheric Lyα absorption diagnostic (see Section 4) is roughly two orders of magnitude more sensitive than the X-ray diagnostic and roughly three orders of magnitude more sensitive than the radio measurement.

4 Detecting Winds Through Astrospheric Absorption

4.1 Analyzing H I Lyman-alpha lines

Spectroscopic analyses of stellar H I Lyman-α lines have proven to be the best way so far to clearly detect and measure weak solar-like winds, but analysis of this line is complex. The history of Lyα absorption observations and analyses is summarized here, with emphasis on how these studies eventually led to the detection of solar-like stellar winds.

The Lyα line at 1216 Å is the most fundamental transition of the most abundant atom in the universe. As such, the line is a valuable diagnostic for many purposes. For cool stars, the Lyα line is a strong radiative coolant for stellar chromospheres, and is therefore an important chromospheric diagnostic. Cool stars produce very little continuum emission at 1216 Å, so Lyα lines observed from cool stars are strong, isolated emission lines. However, these lines are always heavily contaminated by interstellar absorption.

Figure 5 shows the Lyα lines observed from the two members of the a Cen binary, which is the nearest star system to us at a distance of only 1.3 pc. The centers of the stellar emission lines are obliterated by broad, saturated H I absorption. Narrower absorption from deuterium (D I) is visible −0. 33 Å from the H I absorption. The D I absorption is entirely from interstellar material in between α Cen and the Sun, and interstellar absorption also accounts for much of the H I absorption. The H I and D I absorption lines are therefore valuable diagnostics for the LISM, and stellar Lyα data like that in Figure 5 have proven to be useful for measuring the properties of warm, neutral ISM gas, and also for mapping out the distribution of this gas in the vicinity of the Sun (see Linsky et al., 2000; Redfield and Linsky, 2000).

Figure 5:
figure 5

HST/GHRS spectra of the Lyα lines of α Cen A and B, showing broad absorption from interstellar H I and narrow absorption from D I (Linsky and Wood, 1996).

The importance of these data is magnified further by the measurements the data provide of D/H ratios. The LISM D/H ratio has important applications for both cosmology and our understanding of Galactic chemical evolution (see McCullough, 1992; Linsky, 1998; Moos et al., 2002; Wood et al., 2004). The light element abundances in the universe are powerful diagnostics for Big Bang nucleosynthesis calculations, with the exact abundances depending on the cosmic baryon density, Ωb. The deuterium abundance is particularly sensitive to Ωb, so measuring the primordial D/H ratio and thereby constraining cosmological models has been a major goal in astronomy (see Boesgaard and Steigman, 1985; Burles et al., 2001). The most accurate meaningful D/H measurements come from analyses of LISM absorption like that in Figure 5. Unfortunately, the LISM D/H ratio only provides a lower limit to the primordial D/H ratio. Since deuterium is destroyed in stellar interiors, the D/H ratio is expected to have decreased with time, so Galactic chemical evolution models are required to extrapolate back to a primordial value (see Prantzos, 1996; Tosi et al., 1998; Chiappini et al., 2002). Primordial D/H values can also be measured more directly by measuring D/H in more pristine intergalactic material (see Kirkman et al., 2003).

The desire to improve our understanding of the LISM and measure D/H has provided the primary impetus behind stellar Lyα analyses. This work dates back to Copernicus (York and Rogerson, 1976; Dupree et al., 1977), which was the first ultraviolet astronomical satellite to provide high quality stellar Lyα spectra. Copernicus was followed by the long-lived International Ultraviolet Explorer (IUE) satellite, which also provided Lyα spectra that could be analyzed for LISM and D/H purposes (see Murthy et al., 1987; Diplas and Savage, 1994). However, it was the Goddard High Resolution Spectrograph (GHRS) instrument aboard HST that was the first UV spectrometer capable of fully resolving the D I and H I absorption line profiles. The GHRS spectrometer was replaced by the Space Telescope Imaging Spectrometer (STIS) in 1997.

The first Lyα analyses from HST data were for the lines of sight to Capella and Procyon (Linsky et al., 1993, 1995). However, the third analysis, which was of the α Cen data shown in Figure 5, presented a dilemma. The observed H I absorption is simply inconsistent with D I and other ISM absorption lines (Mg II, Fe II, etc.). The D I absorption and the other non-H I lines are centered at a heliocentric velocity of v = −18.0 ± 0.2 km s−1, and the widths of these lines suggest an interstellar temperature of T = 5400 ± 500 K. However, the H I absorption implies v = −15.8 ± 0.2 km s−1 and T = 8350K. In other words, the H I absorption is broader than it should be, and it is also redshifted by 2.2 km s−1 from where it should be. (Linsky and Wood, 1996). Interestingly enough, the H I redshift was also discerned earlier in much lower quality IUE spectra (Landsman et al., 1984).

Linsky and Wood (1996) interpreted the problem as being due to the presence of an extra H I absorption component that contributes no absorption to any of the weaker ISM lines. This was also one of the hypotheses suggested by Landsman et al. (1984) based on the IUE data. Two-component fits to the HST/GHRS data from Linsky and Wood (1996) suggest that the extra H I absorption can be explained by the existence of a redshifted H I absorption component with a temperature of T ≈ 30000K and a column density (in cm−2) of log N(H I) ≈ 15.0. This temperature is much hotter than typical LISM material, and the column density is almost three orders of magnitude lower than the ISM H I column density towards α Cen. The low column density would explain why the absorption component is only seen in the H I line and not in any of the other ISM lines from atomic species with much lower abundances. However, the interpretation for this absorption component was initially a mystery.

Fortuitously, the α Cen Lyα analysis was being performed at about the same time as the first heliospheric models including neutrals in a self-consistent manner were being developed (see Section 2.3). One session of the 1995 IUGG (International Union of Geodesy and Geophysics) General Assembly brought together interstellar and heliospheric experts, and it was realized during that meeting that the heated heliospheric hydrogen predicted by the new heliospheric models had precisely the right properties to account for the extra a Cen absorption component. It was quickly realized that if hot hydrogen existed around the Sun, then it should exist around other solar-like stars as well, so the initial a Cen analysis suggested that the excess H I absorption could be partly due to astrospheric as well as heliospheric absorption (Linsky and Wood, 1996). Thus was born a new way to detect and study stellar winds.

More work was required to verify this interpretation of the Lyα data. Gayley et al. (1997) made the first direct comparison between the α Cen Lyα data and the predictions of various heliospheric models. This work demonstrated that heliospheric Lyα absorption can only account for the excess absorption seen on the red side of the Lyα line, and that astrospheric absorption is required to explain the excess absorption seen on the blue side of the line. A more refined interpretation for the a Cen Lyα absorption was therefore developed, and is shown schematically in Figure 6. The four middle panels in Figure 6 show what happens to the α Cen B Lyα line profile as it journeys from the star towards the Sun. The Lyα profile first has to traverse the hot hydrogen in the star’s astrosphere, which erases the central part of the line. The Lyα emission then makes the long interstellar journey, resulting in additional absorption, including some from D I. Finally, the profile has to travel through the heliosphere, resulting in additional absorption on the red side of the line. Most of the astrospheric and heliospheric absorption is from material in the “hydrogen wall” region mentioned at the end of Section 2.3 (see Figure 4).

Figure 6:
figure 6

Schematic diagram showing how a stellar Lyα profile changes from its initial appearance at the star and then through various regions that absorb parts of the profile before it reaches an observer at Earth: the stellar astrosphere, the LISM, and finally the heliosphere (Wood et al., 2003b). The lower panel shows the actual observed Lyα profile of a Cen B. The upper solid line is the assumed stellar emission profile and the dashed line is the ISM absorption alone. The excess absorption is due to heliospheric H I (green shading) and astrospheric H I (red shading).

Why is the heliospheric absorption redshifted relative to the interstellar absorption? For the α Cen line of sight, which is roughly in the upwind direction relative to the LISM flow seen by the Sun, the heliospheric absorption is redshifted primarily because of the deceleration and deflection of interstellar material as it crosses the bow shock. Heliospheric models predict that heliospheric Lyα absorption should always be redshifted relative to the ISM absorption, even in downwind directions, although the physical explanation for the redshift in the downwind direction is more complicated (Izmodenov et al., 1999b; Wood et al., 2000b). Conversely, astrospheric absorption will always be blueshifted relative to the ISM absorption, since we are viewing that absorption from outside the astrosphere rather than inside. It is very fortunate that heliospheric and astrospheric material produce excess absorption on opposite sides of the Lyα line, as this makes it possible to identify the source of the absorption.

The bottom panel of Figure 6 shows the observed Lyα profile of a Cen B (Linsky and Wood, 1996). As mentioned above, the non-H I ISM lines observed towards α Cen suggest υ = −18.0 ± 0.2 km s−1 and T = 5400 ± 500 K for the ISM material in this direction. The dashed line in the bottom panel of Figure 6 shows what the H I absorption looks like when forced to be consistent with these results. No matter what is assumed for the stellar line profile, and no matter what is assumed for the ISM H I column density, there is always excess absorption on both sides of the H I absorption feature that cannot be explained by ISM absorption. As suggested above, the red side excess is best interpreted as heliospheric absorption and the blue side excess is best interpreted as astrospheric absorption.

This example illustrates how heliospheric and astrospheric absorption is detected. The ISM H I absorption is estimated by forcing the H I fit parameters to be consistent with D I and other ISM lines. In many cases, this still leads to excellent fits to the data, but in some cases there is evidence for excess H I absorption on one or both sides of the line, indicating the presence of heliospheric and/or astrospheric absorption. The Lyα analyses can be simplified even further if one assumes that D/H = 1.5 × 10−5, in addition to forcing υ and T to be consistent for D I and H I. This assumption should be valid for nearby stars, since recent work suggests that D/H ≈ 1.5 × 10−5 throughout the Local Bubble, with no evidence for variation (Linsky, 1998; Moos et al., 2002; Wood et al., 2004).

The heliospheric/astrospheric interpretation of the excess Lyα absorption has strong theoretical support, but additional evidence for the validity of this interpretation is still valuable. The best purely empirical demonstration that it is correct comes by comparing the Lyα absorption observed towards a Cen with that observed towards a distant M dwarf companion of the α Cen system called Proxima Cen (Wood et al., 2001). This comparison is made in Figure 7. The Lyα absorption profiles agree well on the red side of the line where the heliospheric absorption is located. However, the blue-side excess absorption seen towards α Cen is not seen towards Proxima Cen. This means that the blue-side excess absorption seen towards α Cen has to be from circumstellar material surrounding α Cen that does not extend as far as the distant companion Proxima Cen (∼ 12 000 AU away), consistent with the astrospheric interpretation. Apparently, Proxima Cen must have a weaker wind than the α Cen binary, which results in a much smaller astrosphere and much less astrospheric Lyα absorption. (The two α Cen stars are close enough that they will share the same astrosphere, and the astrospheric absorption will therefore be characteristic of the combined winds of both stars.) This example suggests how the astrospheric absorption might be used as a diagnostic for the mass loss rates of solar-like stars, which is a subject that is discussed in detail in Section 4.3.

Figure 7:
figure 7

Comparison between the Lyα spectra of α Cen B (green histogram) and Proxima Cen (red histogram) from Wood et al. (2001). The inferred ISM absorption is shown as a green dashed line. The Alpha/Proxima Cen data agree well on the red side of the H I absorption, but on the blue side the Proxima Cen data do not show the excess Lyα absorption seen toward α Cen (i.e. the astrospheric absorption).

4.2 Comparing heliospheric absorption with model predictions

In this article, we are more interested in the astrospheric Lyα absorption than the heliospheric absorption, since we are concerned with using the astrospheric detections to measure the properties of solar-like stellar winds. However, hydrodynamic models of the astrospheres are required to do this (see Section 4.3). In order to believe the results of these models, it is crucial to demonstrate that they are properly extrapolated from heliospheric models that can reproduce observed heliospheric absorption. Thus, the heliospheric absorption and efforts to reproduce it using models are reviewed in this section.

Gayley et al. (1997) were the first to clearly demonstrate that heliospheric models are capable of reproducing the observed excess Lyα absorption, at least on the red side of the line (see Section 4.1). However, the exact amount of absorption predicted depends on exactly what properties are assumed for the surrounding LISM. Since many of these LISM properties are not precisely known (see Section 2.2), there is the hope that the heliospheric absorption can be used as a diagnostic for the LISM properties, such that the heliospheric absorption is only reproduced when the correct LISM parameters are assumed. In practice, this has proven to be very difficult.

One problem has been finding enough detections of heliospheric Lyα absorption to provide proper constraints for the models, although the situation has been gradually improving. The first heliospheric absorption detection was for the α Cen line of sight described in detail in Section 4.1. The second detection was for the downwind line of sight towards Sirius (Izmodenov et al., 1999b), though the analysis of Hébrard et al. (1999) suggests that this detection is not very secure. The upwind line of sight towards 36 Oph (Wood et al., 2000a) provided the third detection, and there is a marginal detection for the crosswind line of sight towards HZ 43 (Kruk et al., An HST archival Lyα survey by Wood et al. (2005b) resulted in four more detections (70 Oph, ξ Boo, 61 Vir, and HD 16. Most recently, a detailed analysis of all reconstructed stellar Lyα lines based on HST data has found evidence for very broad, weak heliospheric absorption for three lines of sight (χ1 Ori, HD 28205, HD 28568) observed within 20° of the downwind direction (Wood et al., 2007b). This brings the grand total of absorption detections to 11. In addition, Lemoine et al. (2002) and Vidal-Madjar and Ferlet (2002) have claimed to find evidence for weak heliospheric absorption towards the similar Capella and G191-B2B lines of sight, but these claims rely on subtle statistical arguments rather than clearly visible excess absorption.

The heliospheric absorption detections can be supplemented with other lines of sight observed by HST that at least provide useful upper limits for the amount of heliospheric absorption. Figure 8 shows the Lyα absorption profiles observed for three of the lines of sight with detected heliospheric absorption (36 Oph, α Cen, and Sirius) and three additional lines of sight with nondetections (31 Com, β Cas, and ϵ Eri). The figure zooms in on the red sides of the profiles where the heliospheric absorption should be located. The θ angles in the figure indicate the angle of the line of sight with respect to the upwind direction of the interstellar flow. The six sight lines sample diverse orientiation angles, ranging from the nearly upwind direction towards 36 Oph (θ = 12°) to the nearly downwind line of sight to ϵ Eri (θ = 148°). The dotted lines show the ISM absorption alone, constrained by forcing consistency with D I (see Section 4.1) — only 36 Oph, α Cen, and Sirius show excesses that reveal the presence of heliospheric absorption. A successful heliospheric model should accurately reproduce the amount of heliospheric absorption detected towards these three stars, while predicting no detectable absorption towards the other three.

Figure 8:
figure 8

Comparison of the H I absorption predicted by a four-fluid heliospheric model (dashed lines) and the observations, where the model heliospheric absorption is shown after having been added to the ISM absorption (dotted lines). Reasonably good agreement is observed, although there is a slight underprediction of absorption towards 36 Oph and Sirius, and a slight overprediction towards ϵ Eri (Wood et al., 2000b).

Figure 8 shows a model that agrees reasonably well with the data, despite a slight overprediction of absorption towards ϵ Eri and slight underpredictions for 36 Oph and Sirius (Wood et al., 2000b). This model assumes T = 8000K, n(HI) = 0.14 cm−3, and n(H+) = 0.10 cm−3 for the ambient LISM, parameters well within the range of values inferred by other means (see Section 2.2). However, not all heliospheric models that assume these parameters find agreement with the data.

This brings us to the second problem with trying to infer ambient LISM parameters from the heliospheric Lyα data: Results currently seem to be very model dependent. It is mentioned in Section 2.3 just how difficult it is to properly consider neutrals in heliospheric models due to charge exchange processes driving the neutral H out of thermal equilibrium. The model used in Figure 8 is a “four-fluid” model of the type developed by Zank et al. (1996), where one fluid represents the protons, and three distinct fluids are used to represent the neutral hydrogen, one fluid for each distinct region where charge exchange occurs (inside the TS, between the TS and HP, and between the HP and BS). However, there are other approaches, such as the hybrid kinetic code of Müller et al. (2000) and the Monte Carlo kinetic code of Baranov and Malama (1993, 1995). The heliospheric absorption predicted by these kinetic models is not identical to that predicted by the four-fluid models (Wood et al., 2000b; Izmodenov et al., 2002). Currently the kinetic models seem to have more difficulty reproducing the observed heliospheric absorption than the four-fluid models, especially in downwind directions where they tend to predict too much absorption. However, the kinetic models should in principle yield more accurate velocity distributions for the neutral H than codes with multi-fluid approximations. A complex multi-component treatment of the protons in the heliosphere seems to improve the kinetic models’ ability to fit the data (Malama et al., 2006; Wood et al., 2007b). Clearly more work is required to attain some sort of convergence in the models before LISM parameters can be unambiguously derived from the data.

However, a third difficulty with using the heliospheric absorption to infer ambient LISM properties is that the absorption may not be as sensitive to these properties as one might wish. Izmodenov et al. (2002) experiment with different LISM proton and neutral hydrogen densities and find surprisingly little change in the predicted Lyα absorption, at least in upwind and sidewind directions. This may be bad news for the diagnostic power of the heliospheric absorption, but it is actually good news for the astrospheric analyses that are described in Section 4.3. In using astrospheric models to help extract stellar mass loss rates from the astrospheric absorption, one has to assume that the LISM does not vary much from one location to another. The results of Izmodenov et al. (2002) suggest that the models are not very sensitive to the modest variations in LISM properties that one might expect to be present in the solar neighborhood.

Finally, there are some aspects of heliospheric physics that are only beginning to be considered in the models. The models mentioned previously do not consider either the heliospheric magnetic field carried outwards by the solar wind (see Nerney et al., 1991), or the poorly known interstellar magnetic field. Not only is a proper MHD treatment of the heliosphere difficult, but the problem is inherently three dimensional, whereas the models mentioned previously assume a 2D axisymmetric geometry. Using a 2D approach, Florinski et al. (2004) find that a strong ISM field oriented parallel to the LISM flow does not yield significantly different predictions for heliospheric Lyα absorption than models without magnetic fields. However, 3D models are required to include the heliospheric field, and 3D models are also required to consider ISM field orientations other than parallel to the flow vector.

Initial 3D models developed to model these effects (see Linde et al., 1998) did not include neutrals in a self-consistent manner. Dealing with both neutrals and magnetic fields properly in a 3D model is a very formidable problem. Nevertheless, the 3D models without neutrals do suggest that MHD effects could in principle lead to changes in the heliospheric structure that could affect the Lyα absorption. Examples include the unstable jet sheet and north-south asymmetries predicted by Opher et al. (2003, 2004, 2006). In addition, Ratkiewicz et al. (1998) find that if the LISM magnetic field is skewed with respect to the ISM flow, the effective nose of the heliosphere could be significantly shifted from the upwind direction. Even in the absence of magnetic fields, latitudinal variations in solar wind properties could also cause asymmetries in the heliosphere (Pauls and Zank, 1997). It is possible that all these asymmetries suggested by 3D models could be detectable in Lyα absorption. However, neutrals must be included properly in the models to make clearer predictions. Only very recently has this been done (Izmodenov et al., 2005; Pogorelov et al., 2006). Wood et al. (2007a) have made the first comparison between such models and the Lyα data, finding that the absorption predicted by the models is modestly affected by the assumed LISM field strength and orientation, allowing some constraints on these quantities to be inferred from the data. However, with the absorption being modestly dependent on so many uncertain LISM properties, including particle densities, it is probably unreasonable to expect analysis of the absorption by itself to yield a single set of acceptable LISM parameters.

4.3 Measuring stellar mass loss rates

The α Cen line of sight provided the first detection of astrospheric Lyα absorption (Linsky and Wood, 1996), but the use of models by Gayley et al. (1997) was necessary to clearly demonstrate that heliospheric absorption could not explain the blue-side excess as well as the stronger redside excess (see Figure 6). By that time, there were already two other lines of sight, ϵ Ind and λ And, with excess Lyα absorption found only on the blue side of the absorption line, which clearly could not be heliospheric in origin (Wood et al., 1996). This excess absorption was immediately interpreted as being solely astrospheric. Thus, a case could be made for either α Cen AB or ϵ Ind/λ And being the first stars with detected solar-like coronal winds.

There are now a total of 13 published detections of astrospheric absorption. These detections are listed in Table 1, in addition to Proxima Cen since the nondetection for Proxima Cen is used to derive an upper limit for its mass loss rate (see below).

Table 1: Published Astrospheric Detections and Mass Loss Measurements

Some of these stars are binaries where the individual components are clearly close enough that both stars will reside within the same astrosphere, meaning that the observed astrospheric absorption is indicative of the combined mass loss of both stars. For these binaries (α Cen, 36 Oph, and λ And), the spectral types of both stars are listed in Table 1 and the listed stellar surface areas are the combined areas of both stars. In contrast, 61 Cyg A’s companion is far enough away that 61 Cyg A should have an astrosphere all to itself, so it is listed alone in Table 1.

Three detections are flagged as uncertain in Table 1 for various reasons that will not be discussed here (see Wood et al., 2002, 2005b). Results regarding these three lines of sight should be considered with caution. There was originally another uncertain astrospheric detection, 40 Eri A (Wood and Linsky, 1998; Wood et al., 2002), which has now been dropped entirely from the list of astrospheric detections. Since astrospheric hydrogen scatters stellar Lyα photons, solar-like stars should in principle be surrounded by faint nebulae of astrospheric Lyα emission. There was an unsuccessful attempt to detect this emission surrounding 40 Eri A, and the lack of success was used to argue that the tentative detection of astrospheric absorption for 40 Eri A could not be correct (Wood et al., 2003a).

The astrospheric absorption detections represent an indirect detection of solar-like winds from the observed stars, and the amount of absorption observed has diagnostic power. Frisch (1993) envisioned using astrospheres as probes for the interstellar medium, but so far work has focused on the use of the astrospheric absorption for estimating stellar mass loss rates. The higher the stellar mass loss rate, the larger the astrosphere will be, and the more Lyα absorption it will produce. However, extracting quantitative mass loss measurements from the data requires the assistance of astrospheric models analogous to the heliospheric models discussed in Section 2.3.

The first step in modeling an astrosphere is to determine what the interstellar wind is like in the rest frame of the star. The proper motions and radial velocities of the nearby stars in Table 1 are known very accurately, as are their distances (see Perryman et al., 1997). The ISM flow vector is generally assumed to be the same as the LIC flow vector seen by the Sun (Lallement et al., 1995). Although multiple ISM velocity components are often seen towards even very nearby stars, the components are never separated by more than 5–10 km s−1, meaning that the LIC vector should be a reasonable approximation for these other clouds. An example is the nearby “G” cloud. Since a Cen, 70 Oph, and 36 Oph are known to lie within this cloud rather than the LIC, the “G” cloud vector from Lallement and Bertin (1992) is used instead of the LIC vector, but this does not change things very much. In any case, Table 1 lists the ISM wind velocity seen by each star (VISM). The θ value in Table 1 indicates the angle between our line of sight to the star and the upwind direction of the ISM flow seen by the star.

Astrospheric models are extrapolated from a heliospheric model that fits the observed heliospheric Lyα absorption. The principle heliospheric model used in the past for this purpose is the four-fluid model described in Section 4.2, which is the source of the predicted heliospheric absorption in Figure 8. In order to convert this to an astrospheric model, the model is recomputed using the same wind and ISM input parameters, except for the ISM wind speed which is changed to the VISM value appropriate for the star (see Table 1). The stellar wind proton density is varied to experiment with mass loss rates different from that of the Sun. Figure 9 shows four models of the α Cen astrosphere computed assuming different mass loss rates in the range of = 0.2–2 (Wood et al., 2001). The model astrospheres naturally become larger as the mass loss is increased. The figure shows the H I density distribution, but the models also provide temperature distributions and flow patterns. From these models astrospheric absorption can be computed for the θ = 79° line of sight to the star, for comparison with the data. Figure 10 shows this comparison for α Cen, along with similar comparisons for five other stars in Table 1 (Wood et al., 2002). The model with M = 2 agrees best with the data, so this is the estimated mass loss rate of α Cen. Figure 11 shows similar data-model comparisons for six additional stars from Table 1 (Wood et al., 2005a).

Figure 9:
figure 9

Distribution of H I density predicted by hydrodynamic models of the Alpha/Proxima Cen astrospheres, assuming stellar mass loss rates of (from top to bottom) 0.2, 0.5, 1.0, and 2.0 (Wood et al., 2001). The distance scale is in AU. Streamlines show the H I flow pattern.

Figure 10:
figure 10

Closeups of the blue side of the H I Lyα absorption lines for six stars with detected astrospheric absorption, plotted on a heliocentric velocity scale. Narrow D I ISM absorption is visible in all the spectra just blueward of the saturated H I absorption. Green dashed lines indicate the interstellar absorption alone, and blue lines in each panel show the additional astrospheric absorption predicted by hydrodynamic models of the astrospheres assuming various mass loss rates (Wood et al., 2002).

Figure 11:
figure 11

A figure analogous to Figure 10 , but for six other lines of sight (Wood et al., 2005a).

Table 1 lists all mass loss rates measured in this manner, and Figures 12 and 13 show the astrospheric models that lead to the best fits to the data. The astrospheres vary greatly in size, both due to different mass loss rates and to different ISM wind speeds. It is worth noting that the largest of these astrospheres (ϵ Eri, 70 Oph) would be comparable in size to the full Moon in the night sky, if we could see them.

Figure 12:
figure 12

Maps of H I density from hydrodynamic models of stellar astrospheres (Wood et al., 2002). The models shown are the ones that lead to the best fits to the data in Figure 10 . The distance scale is in AU. The star is at coordinate (0,0) and the ISM wind is from the right. The dashed lines indicate the Sun-star line of sight.

Figure 13:
figure 13

Maps of H I density from hydrodynamic models of stellar astrospheres (Wood et al., 2005a), analogous to Figure 12 . The models shown are the ones that lead to the best fits to the data in Figure 11 .

Given the model dependence and other difficulties described in Section 4.2, one might wonder if mass loss rates derived using the astrospheric models are at all reliable. However, the crucial point is that the heliospheric model from which the astrospheric models are extrapolated successfully reproduces the heliospheric absorption. For the astrospheric modeling purposes described here, this is more important than whether the input parameters of that model are actually correct. In essence, the observed heliospheric absorption is calibrating the models before they are being applied to modeling astrospheres, and in this way the mass loss measurement technique is actually semi-empirical. Nevertheless, the mass loss rates measured using the astrospheric Lyα absorption technique should not be considered to be terribly precise. Uncertainties in the mass loss rates are probably of order a factor of 2 (Wood et al., 2002).

Other assumptions are implicitly made in this mass loss measurement procedure. One is that the LISM does not vary greatly from one star to the next. However, for these very nearby stars LISM variations should be modest, and as discussed in Section 4.2 modest variations probably do not greatly affect the predicted astrospheric absorption. Another assumption is that the stellar wind speeds of the stars in Table 1 are similar to that of the Sun. An argument in favor of this assumption is that stellar wind velocities are generally not very different from the surface escape speeds based on past experience, regardless of the type of wind or star one is considering (see Section 2.1). This includes the Sun, which has a surface escape speed of 619 km s−1, very similar to observed solar wind velocities. All the main sequence stars in Table 1 have similar escape speeds, so one might expect them to have similar wind velocities. However, the magneto-centrifugal wind acceleration models of Holzwarth and Jardine (2007) suggest that rapidly rotating stars might have wind speeds much faster than their surface escape speeds. Thus, the appropriateness of the constant wind speed assumption is still a matter for debate. The one star in Table 1 that has an escape velocity significantly different from that of the Sun is λ And. The G8 IV-III primary that surely dominates the wind from this binary is not a main sequence star, and it has a surface escape speed about 3 times lower than that of the Sun. Perhaps in the future the mass loss rate listed in Table 1 should be remeasured assuming a lower wind velocity.

In addition to the mass loss measurements for stars with detected astrospheric absorption, Table 1 also lists an upper limit derived from a nondetection for Proxima Cen. Upper limits cannot be computed for most stars with nondetections, because a likely interpretation for most nondetections is that the star is surrounded by hot, fully ionized ISM material, rather than partially neutral gas like that which surrounds the Sun. The Local Bubble in which the Sun is located is mostly filled with this hot material (see Section 2.2). The Sun just happens to be located within one of the cooler, partially neutral clouds that lie within the Bubble. For Proxima Cen, the hot ISM explanation for the astrospheric nondetection can be discarded, because Proxima Cen’s distant companion α Cen has detected astrospheric absorption, proving that Proxima Cen is not located within the hot ISM. Thus, for Proxima Cen a meaningful upper limit to the stellar mass loss rate can be derived.

No mass loss measurement is reported for HD 128987 in Table 1, since no astrospheric model seems to be able to reproduce the apparent astrospheric absorption, bringing into question the astrospheric interpretation of the absorption for this star (Wood et al., 2005a). The fundamental problem is the very low LISM wind speed seen by this star (VISM = 8 km s−1; see Table 1). This low speed is not sufficient to result in much deceleration and heating of neutral H within the astrosphere, thereby yielding little predicted absorption. Thus, the apparent astrospheric detection for HD 128987 remains a mystery at this time.

5 Implications for the Sun and Solar System

5.1 Inferring the mass loss history of the Sun

In addition to stellar mass loss rates, Table 1 lists coronal X-ray luminosities from the ROSAT PSPC instrument (see Hünsch et al., 1999; Schmitt and Liefke, 2004). Solar-like winds have their origins in stellar coronae (see Section 2.1), so one might expect the winds to be correlated with coronal properties such as X-ray emission. Thus, in Figure 14 the mass loss rates measured from the astrospheric Lyα absorption (per unit surface area) are plotted versus X-ray surface fluxes (Wood et al., 2005a). For the main sequence stars mass loss increases with coronal activity. A power law is fitted to these GK stars in Figure 14. Quantitatively, this relation is

$$\dot M \propto F_{\rm{X}}^{1.34 \pm 0.18}.$$
(1)

The saturation line in Figure 14 indicates the maximum X-ray flux observed from solar-like stars (Güdel et al., 1997).

Figure 14:
figure 14

Measured mass loss rates (per unit surface area) plotted versus X-ray surface flux (Wood et al., 2005a). The filled and open circles are main sequence and evolved stars, respectively. For the main sequence stars with log FX < 8 × 105 ergs cm−2 s−1, mass loss appears to increase with coronal activity, so a power law has been fitted to these stars, and the shaded region is the estimated uncertainty in the fit. The saturation line represents the maximum FX value observed from solar-like stars.

It is interesting to note that during the solar cycle, the Sun’s wind strength is actually anticorrelated with its X-ray flux. The solar wind is weaker at solar maximum than at solar minimum despite coronal X-ray fluxes being much higher (Lazarus and McNutt Jr, 1990). This is presumably due to the fact that winds are more associated with the large scale dipole component of the solar magnetic field instead of the small scale active regions responsible for most of the Sun’s X-ray emission. The dipole field actually weakens at solar maximum along with the wind. However, the interior magnetic dynamo is ultimately responsible for both the small scale and large scale fields, so as a whole both field components should increase with increasing dynamo activity, consistent with the mass loss/activity correlation in Figure 14 (Schrijver et al., 2003).

The evolved stars are clearly inconsistent with the main sequence stars in Figure 14. The very active coronae of λ And and DK UMa produce surprisingly weak winds, though it should be noted that both of these astrospheric detections are flagged as being questionable in Table 1. There are three main sequence stars with log FX > 8 × 105 ergs cm−2 s−1, which have low mass loss measurements that are not consistent with the wind-activity correlation that seems to exist for the low activity main sequence stars. Two of these stars (Proxima Cen and EV Lac) are tiny M dwarf stars. If these were the only discrepant data points one could perhaps argue that the discrepancy is due to these M dwarfs being significantly less solar-like than the G and K dwarfs that make up the rest of the main sequence sample of stars. However, this interpretation is invalidated by the third discrepant measurement, that of ξ Boo. Being a binary with two rather solar-like stars (G8 V+K4 V), there is no easy way to dismiss the ξ Boo measurement, which implies that the power law relation does not extend to high activity levels for any type of star. More mass loss measurements of active stars would clearly be helpful to better define the characteristics of solar-like winds at high coronal activity levels.

Based on the available data, the mass-loss/activity relation appears to change its character at log FX ≈ 8 × 105 ergs cm−2 s−1. One possible explanation for this concerns the existence of polar spots for very active stars. Low activity stars presumably have starspot patterns like that of the Sun, where spots are confined to low latitudes. However, for very active stars not only are spots detected at high latitudes, but a majority of these stars show evidence for large polar spots (Strassmeier, 2002). The existence of high latitude and polar spots represents a fundamental change in the stellar magnetic geometry (Schrijver and Title, 2001), and it is possible that this dramatic change in magnetic field structure could affect the winds emanating from these stars. Perhaps stars with polar spots might have a magnetic field with a strong dipolar component that could envelope the entire star and inhibit stellar outflows, thereby explaining why active stars have weaker winds than the mass-loss/activity relation of less active main sequence stars would predict. For ξ Boo A, high latitude spots of some sort have been detected (Toner and Gray, 1988). Petit et al. (2005) have detected a strong global dipole field component for ξ Boo A, consistent with the picture presented above. They also detected a large-scale toroidal field component, which would have no solar analog whatsoever, consistent with the idea that very active solar-like stars have significantly different magnetic field structures from those of the Sun and other low-activity stars.

Figure 14 illustrates how mass loss varies with coronal activity. But what about age? There is a known connection between activity and age, for the following reasons. The gravitational contraction of interstellar clouds that results in star formation leads to rapid rotation for young, newly born stars. This rapid rotation leads to vigorous dynamo activity and therefore high surface magnetic activity and high coronal X-ray emission. However, the magnetic fields of these young, rapidly rotating stars drag against their winds, and this magnetic braking gradually slows the stellar rotation. This in turn leads to lower activity levels and X-ray fluxes. An enormous amount of effort has been expended in the past few decades to observationally establish exactly how rotation relates to stellar age (see Skumanich, 1972; Soderblom et al., 1993) and how rotation relates to stellar activity, which is most easily measured through X-ray emission (see Pallavicini et al., 1981; Walter, 1982, 1983; Caillault and Helfand, 1985; Micela et al., 1985; Fleming et al., 1989; Stauffer et al., 1994). For solar-like stars, Ayres (1997) finds

$${V_{{\rm{rot}}}} \propto {t^{- 0.6 \pm 0.1}}$$
(2)

for the rotation/age relation, while X-ray flux and rotation are related by

$${F_{\rm{X}}} \propto V_{{\rm{rot}}}^{2.9 \pm 0.3}.$$
(3)

Equations (1), (2), and (3) can be combined to obtain the following relation between mass loss and age for solar-like stars:

$$\dot M \propto {t^{- 2.33 \pm 0.55}}.$$
(4)

This is the first empirically determined mass loss evolution law for solar-like stars, and Figure 15 shows what this relation implies for the mass loss history of the Sun in particular (Wood et al., 2005a).

Figure 15:
figure 15

The mass loss history of the Sun suggested by the power law relation from Figure 14 (Wood et al., 2005a). The low mass-loss rate measurement for ξ Boo implies that the wind weakens at t ≈ 0.7 Gyr as one goes back in time.

The truncation of the power law relation in Figure 14 leads to the truncation of the mass-loss/age relation in Figure 15 at about t = 0.7 Gyr. The location of ξ Boo is shown in order to infer what the solar wind might have been like at earlier times. Despite the high activity cutoff, the mass loss measurements obtained so far clearly suggest that winds are generally stronger for young solar-like stars, and as a consequence the solar wind was presumably much stronger early in the Sun’s lifetime. This has many important implications, some of which are discussed in Sections 5.2, 5.3, and 5.4.

5.2 Magnetic braking

The magnetic braking process by which stars like the Sun shed angular momentum is the reason why stellar activity decreases with time and why the solar mass loss/age relation in Figure 15 can be inferred from the mass loss/activity relation in Figure 14 (see Section 5.1). Winds play an important role in this process, because it is the wind that the stellar field drags against in slowing down the star’s rotation. The efficiency of this braking is clearly related to the density of the wind, and therefore to the mass loss rate. Thus, the wind evolution law in Equation (4) has consequences for how the effectiveness of the magnetic braking changes with time. Models for magnetic braking suggest relations of the form

$${{\dot \Omega} \over \Omega} \propto {{\dot M} \over M}{\left({{{{R_{\rm{A}}}} \over R}} \right)^m},$$
(5)

where Ω is the angular rotation rate and RA is the Alfvén radius (Weber and Davis Jr, 1967; Stepień, 1988; Gaidos et al., 2000). The exponent m is a number between 0 and 2, where m = 2 corresponds to a purely radial magnetic field. Mestel (1984) claims that more reasonable magnetic geometries suggest m = 0–1. The Alfvén radius is

$${R_{\rm{A}}} = \sqrt {{{{V_{\rm{w}}}\dot M} \over {B_{\rm{r}}^2}}},$$
(6)

where Vw is the stellar wind speed and Br is the disk-averaged radial magnetic field.

For a star like the Sun, the star’s mass and radius are relatively invariant. If Vw does not vary with time, which was also an assumption used in the derivation of mass loss rates from the astrospheric absorption (see Section 4.3), then the time dependence of all quantities in Equations (5) and (6) are known except for that of Br. If Br is expressed as a power law, Brtα, Equations (4), (5), and (6) combined suggest

$$\alpha = 1/m - (1.17 \pm 0.28)(m + 2)/m.$$
(7)

Assuming that m is in the physically allowable range of m = 0–2 yields the upper limit α < −1.3, while the more likely range of m = 0–1 suggested by Mestel (1984) implies α < −1.7. In any case, the empirical mass loss evolution law in Equation (4) is consistent with theoretical descriptions of magnetic braking only if disk-averaged stellar magnetic fields decline at least as fast as t−1.3.

5.3 The Faint Young Sun problem

Evolutionary models of the solar interior lead to the prediction that the Sun should have been about 25% fainter 3.8 Gyr ago (Gough, 1981; Bahcall et al., 2001). This is a fairly robust result that has been known for a long time. However, this creates significant difficulties for those who study the climate history of planets in our solar system. With the young Sun this faint, the planets should have been significantly colder than they are now. On Earth and Mars, surface water should have been frozen, but this contradicts geologic evidence that abundant water existed and flowed on the surfaces of both of these planets. This is the so-called “Faint Young Sun” paradox (Sagan and Mullen, 1972; Kasting, 1991). Most attempts to solve this problem have appealed to increased amounts of greenhouse gases in the early atmospheres of Earth and Mars, which allowed the surface temperatures to remain similar to present-day temperatures despite a fainter Sun (see Walker, 1985; Kasting and Ackerman, 1986; Kasting, 1991).

However, a couple theories have been advanced that involve the solar wind. One of these is simply that the young solar wind was strong enough to significantly decrease the mass of the Sun, meaning that the young Sun was more massive and thus more luminous than standard solar evolutionary models predict (Guzik et al., 1987; Sackmann and Boothroyd, 2003). If the Sun was only about 2% more massive 3.8 Gyr ago, then the Sun would still have been bright enough to maintain sufficiently warm temperatures on Earth and Mars. The mass loss evolution law derived from stellar wind measurements in Equation (4) does indeed suggest that the solar wind was stronger in the past. Unfortunately, it was not strong enough. If Equation (4) is correct, the Sun could not have been more than 0.2% more massive 3.8 Gyr ago. Thus, this theory ultimately is not supported by the stellar wind measurements (Minton and Malhotra, 2007).

However, there is a less direct method that has been proposed by which a stronger young solar wind could contribute to a solution of the Faint Young Sun problem. There have been several claims that cosmic rays have important effects on the Earth’s climate by stimulating cloud formation, which is believed to generally cool the Earth’s atmosphere (see Svensmark, 1988). However, cosmic rays are modulated by the solar wind, and a stronger solar wind would reduce the flux of cosmic rays into the Earth’s atmosphere. Thus, the idea is that the stronger wind of the young Sun led to lower cosmic ray fluxes, helping the Earth (and perhaps Mars) to maintain warm temperatures (Shaviv, 2003). However, the whole idea of cosmic rays influencing the terrestrial climate to this degree remains very controversial (see Carslaw et al., 2002).

5.4 Erosion of planetary atmospheres

Even if the solar wind plays no role in the solution to the Faint Young Sun paradox described in Section 5.3, the solar wind may still have had a dramatic effect on planets in the past through the eroding effects of the wind on planetary atmospheres. Solar wind sputtering processes have been proposed as having had important effects on the atmospheres of both Venus and Titan (Chassefière, 1997; Lammer et al., 2000), but it is the Martian atmosphere that may have been the most dramatically affected by the solar wind. The connection between the Martian atmosphere and the question of whether life once existed on the planet makes the Mars example even more interesting.

Mars apparently had running water on its surface in the distant past, but it is now very dry. Isotopic evidence strongly suggests that the Martian atmosphere was once much thicker and more conducive to the existence of surface water (see Carr, 1996; Jakosky and Phillips, 2001). What caused the loss of the early Martian atmosphere, and presumably the loss of surface water and habitability as well? Solar wind sputtering is a leading candidate for the cause of these changes (Luhmann et al., 1992; Perez de Tejada, 1992; Jakosky et al., 1994; Kass and Yung, 1995; Lundin, 2001). Unlike Earth, the Martian atmosphere is not currently protected from the solar wind by a strong magnetosphere. There is evidence that Mars once had a magnetic field, but it disappeared at least 3.9 Gyr ago (Acuña et al., 1999). At that point, Mars would have been exposed to a solar wind about 80 times stronger than the current wind according to Figure 15. The stronger solar wind makes wind erosion an even more likely culprit behind the disappearance of most of the Martian atmosphere.

Studies of Martian atmosphere evolution that consider the wind evolution law inferred from the astrospheric measurements, along with modern estimates of the history of solar X-ray and UV fluxes, have already begun (Guinan et al., 2003; Lammer et al., 2003; Ribas et al., 2005). Investigations into how the solar wind may have affected other planets will undoubtedly also be an active area of future research, stimulated in part by the discovery of planets around other stars. Many of these extrasolar planets orbit very close to their stars, which means that they will be exposed to wind fluxes orders of magnitude higher than the Earth or Mars have ever seen. Knowledge of the evolution of solar-like winds is crucial for understanding the atmospheric evolution of these extrasolar planets, and efforts are already underway to try to model wind erosion effects on such planets (Grießmeier et al., 2004; Preusse et al., 2005).

6 Conclusions

The study of solar-like winds around other stars is not an easy one, due to the difficulty in detecting these winds. However, such measurements are the only way to empirically infer the history of the solar wind and to assess its potential effects on atmospheres of planets in our solar system. At this time, the indirect astrospheric Lyα absorption technique is the only way to detect weak solar-like stellar winds. Radio and X-ray observations can place limits on wind fluxes (see Section 3), but there is no reason to believe that these diagnostics will approach the sensitivity level of the Lyα diagnostic in the near future.

Currently there are only about a dozen mass loss measurements for solar-like stars, a rather small data sample. Additional astrospheric detections are required to better constrain relations between stellar winds, activity, and age. Unfortunately, the STIS instrument on HST failed in August 2004, and at the time of this writing is still unavailable. Furthermore, there is no future mission even in the planning stages that would be capable of the high resolution UV spectrometry necessary to detect absorption signatures of stellar astrospheres. Thus, the future of this subject area is very uncertain, at least on the data side. On the theory side, there is hope that improvements in our ability to numerically model the heliosphere and astrospheres will allow more precise analyses of existing data, and in this way improve our understanding of solar-like stellar winds and astrospheres.