Flow-tube dynamics and coronal holes

We reexamine the well-known polytropic flow-tube model of the expanding solar corona, and find that as the divergence of the flow tube increases the expansion speed increases throughout the flow, over a stated parameter range. Corresponding to a specified flow-tube geometry the terminal speed of the fluid may be far in excess of the value corresponding to purely spherically symmetric flow. The implications of the results for the modelling of high-speed streams emanating from coronal holes are discussed.     

A solar wind model with the power spectrum of Alfvénic fluctuations

A new solar wind model has been developed by including in the model the Alfvénic fluctuation power spectrum equation proposed by Tu et al. (1984). The basic assumptions of the model are as follows: (1) for heliocentric distances r > 10 R _⊙, the radial variation of the power spectrum of Alfvénic fluctuations is controlled by the spectrum equation (1), (2) for heliocentric distances r < 10 R _⊙, the radial variation of the fluctuation amplitude is determined by the Alfvén wave WKB solution, (3) no energy cascades from the low-frequency boundary of the Alfvénic fluctuation power spectrum into the fluctuation frequency range, and the energy which cascades from the high-energy boundary of the spectrum into the higher frequency range is transported to heat of the solar wind flow. Some solutions of this model which, on one hand, describe the major properties of the Alfvénic fluctuations and the high-speed flow observed by Helios in the space range between 0.3–1 AU and, on the other hand, are consistent with the observational constraints at the coronal base have been obtained under the following conditions: (1) the spectrum index of the fluctuations is near to -1 for almost the whole frequency range at 10 R _⊙, (2) the particle flux density at 1 AU is not greater than 3 × 10⁸ cm^?2 s^?1, (3) the solution is for spherically-symmetric flow geometry or the solution passes through the outermost of the three critical points of the rapidly diverging flow geometry with f _max = 7. Some solutions passing through the innermost critical point of the rapidly diverging flow geometry with f _max = 7 have been found, however, with too low pressure at the coronal base to compare with the observational constraints. Heat addition or other kind of momentum addition for r < 10 R _⊙ is required to modify this model to yield better agreement with observations. A cascade energy flux function which leads to Kolmogorov power law in the high-frequency range of Alfvénic fluctuations is presented in Appendix A. More detailed discussions about the characteristics, the boundary conditions and the solution of the spectrum equation (1) are given in Appendix B.     

Parameters of the Magnetic Flux inside Coronal Holes

The parameters of the magnetic flux distribution inside low-latitude coronal holes (CHs) were analyzed. A statistical study of 44 CHs based on Solar and Heliospheric Observatory (SOHO)/MDI full disk magnetograms and SOHO/EIT 284?Å images showed that the density of the net magnetic flux, B _net, does not correlate with the associated solar wind speeds, V _x . Both the area and net flux of CHs correlate with the solar wind speed and the corresponding spatial Pearson correlation coefficients are 0.75 and 0.71, respectively. A possible explanation for the low correlation between B _net and V _x is proposed. The observed non-correlation might be rooted in the structural complexity of the magnetic field. As a measure of the complexity of the magnetic field, the filling factor, f(r), was calculated as a function of spatial scales. In CHs, f(r) was found to be nearly constant at scales above 2 Mm, which indicates a monofractal structural organization and smooth temporal evolution. The magnitude of the filling factor is 0.04 from the Hinode SOT/SP data and 0.07 from the MDI/HR data. The Hinode data show that at scales smaller than 2 Mm, the filling factor decreases rapidly, which means a multifractal structure and highly intermittent, burst-like energy release regime. The absence of the necessary complexity in CH magnetic fields at scales above 2 Mm seems to be the most plausible reason why the net magnetic flux density does not seem to be related to the solar wind speed: the energy release dynamics, needed for solar wind acceleration, appears to occur at small scales below 1 Mm.     

Energy balance of the corona and the origin of quasi-stationary high-speed solar wind streams

The energy balance of open-field regions of the corona and solar wind and the influence of the flow geometry in the corona upon the density and temperature, are analyzed. It is found that the energy flux arriving at the corona is constant for the corona's open regions with different flow geometries. For the waves heating the corona and solar wind, the dependence of the absorption coefficient on the corona's plasma density is found to be within the range of distances r=1.05–1.5R . It is shown that the wave absorption is more dependent on electron density than the coronal emission. It is this difference that causes lower-density coronal holes to be colder than quiet regions. It is found that the additional energy flux necessary for providing energy balance of the corona and for producing solar wind is a flux of Alfvén waves, which can provide the energy needed for producing quasi-stationary high-speed solar wind streams. Theoretical models of coronal holes and the question of why the high-speed solar wind streams are precisely flowing out of coronal holes, are discussed.     

Non-radial solar wind flows and IMF Bz during 1973-2003

The characteristics of latitudinal angles of solar wind flow (θ_v) observed near earth have been studied during the period 1973-2003. The average magnitude of θ_v shows distinct enhancements during the declining and maximum phases of the sunspot cycles. A close association of B_z component of IMF in the GSE system and the orientation of meridional flows in the solar wind is found which depends on the IMF sector polarity. This effect has been studied in typical geomagnetic storm periods. The occurrence of non-radial flows is also found to exhibit heliolatitudinal dependence during the years 1975 and 1985 as a characteristic feature of non-radial solar wind expansion from polar coronal holes.     

Terrestrial low-frequency bursts: Escape paths of radio waves through the bow shock

From the analysis of 119 low-frequency (LF) burst spectra observed onboard the Wind spacecraft, we propose an interpretation of the frequency-time characteristics including the low frequency cutoff of the LF burst spectra, and we use these characteristics to sound the bow shock structure at large tailward distances from Earth. When observed from within the solar wind, LF bursts appear to be made of two spectral components. The high frequency one is bursty and observed above twice the solar wind plasma frequency f_psw. The low frequency one is diffuse (ITKR) and its spectrum extends from about 2f_psw to a cutoff frequency f_c not much higher than f_psw; its onset time δt(f) increases as the frequency f decreases. For each of the 119 events observed from near the Lagrange point L1, the solar wind density variations were measured and the variations of the density jump across the shock calculated from plasma data all along a shock model over more than 2000R_E. But, except for a few events, neither the solar wind overdensities nor the shock density barrier can prevent waves with frequencies below f_c from reaching the spacecraft. Scattering on plasma density inhomogeneities was then introduced to account for the propagation of the LF burst waves in the magnetosheath, from near Earth to their escape point through the bow shock at a frequency-dependent distance |X_esc(f)| (GSE), and then in the solar wind to the spacecraft. In such media, at frequencies between 2f_psw and f_psw, the bulk speed of the scattered waves decreases rapidly as f decreases, and this accounts for the observed variations of the onset time δt(f). Angular scattering can also account for the observed cutoff at f_c if the distance |X_esc(f)| increases exponentially when f/f_psw decreases. As the shock model we used meets that requirement, we consider that this model is valid, which implies that the bow shock still exists beyond 1000R_E from the Earth. The observed decrease of the average spectral intensity of the LF burst between about 1.5f_psw and 2f_psw can also be explained by the scattering in the solar wind if we take into account the angular distribution of the rays when they leave the bow shock.     

Transient response of the solar wind to changes in flow geometry

The transient response of the solar wind to changes in geometry is examined. An initial stationary flow in a configuration that diverges as r ² is assumed. This state corresponds to the usual solar wind solution. The effect on the flow through a tube whose area A(r, t) diverges faster than r ², with the degree of divergence increasing in time, is considered. The asymptotic form of A(r, t) is chosen to mimic the form inferred in coronal holes. A detailed parameter study relating the form of A(r, t) to the pattern of flow in the tube is presented. It is observed that in the limit of large time (large compared to , the time constant for change in geometry of a flow tube) the solutions obtained from a time-dependent analysis can depend upon . For sufficiently large , the asymptotic solution is the same as the steady state solution obeying the correct boundary conditions and possessing a smooth sonic transition. However, if the geometry changes rapidly enough, solutions exhibiting shock-like discontinuities can also exist. This is essentially a new feature that emerges from the present investigation. Finally, it is suggested that this study may be useful in describing flows in evolving coronal holes.     

Small Solar Wind Transients and Their Connection to the Large-Scale Coronal Structure

It has been realized for some time that the slow solar wind with its embedded heliospheric current sheet often exhibits complex features suggesting at least partially transient origin. In this paper we investigate the structure of the slow solar wind using the observations by the Wind and STEREO spacecraft during two Carrington rotations (2054 and 2055). These occur at the time of minimum solar activity when the interplanetary medium is dominated by recurrent high-speed streams and large-scale interplanetary coronal mass ejections (ICMEs) are rare. However, the signatures of transients with small scale-sizes and/or low magnetic field strength (comparable with the typical solar wind value, ～?5 nT) are frequently found in the slow solar wind at these times. These events do not exhibit significant speed gradients across the structure, but instead appear to move with the surrounding flow. Source mapping using models based on GONG magnetograms suggests that these transients come from the vicinity of coronal source surface sector boundaries. In situ they are correspondingly observed in the vicinity of high density structures where the dominant electron heat flux reverses its flow polarity. These weak transients might be indications of dynamical changes at the coronal hole boundaries or at the edges of the helmet streamer belt previously reported in coronagraph observations. Our analysis supports the idea that even at solar minimum, a considerable fraction of the slow solar wind is transient in nature.     

A Two-Dimensional Alfvén-Wave-Driven Solar Wind Model

This paper presents a two-dimensional, Alfvén-wave-driven solar wind model, in which the wave energy is assumed to cascade from the low-frequency Alfvén waves to high-frequency ion cyclotron waves and to be transferred to the solar wind protons by cyclotron resonance at the Kolmogorov rate. A typical structure in the meridional plane consisting of a coronal streamer near the Sun, a fast wind in high latitudes, and a slow wind across the heliospheric current sheet, is found. The fast wind obtained in the polar region is essentially similar to that derived by previous one-dimensional flow-tube models, and its density profile in the vicinity of the Sun roughly matches relevant observations. The proton conditions at 1 AU are also consistent with observations for both the fast and slow winds. The Alfvén waves appear in the fast- and slow-wind regions simultaneously and have comparable amplitudes, which agrees with Helios observations. The acceleration and heating of the solar wind by the Alfvén waves are found to occur mainly in the near-Sun region. It is demonstrated in terms of one-dimensional calculations that the distinct properties of the fast and slow winds are mainly attributed to different geometries of the flow tubes associated with the two sorts of winds. In addition, the 2-D and 1-D simulations give essentially the same results for both the fast and the slow winds.     

First order latitude effects in the solar wind

The Weber-Davis model of the solar wind is generalized to include the effects of latitude. The principal assumptions of perfect electrical conductivity, rotational symmetry, a polytropic relation between pressure and density, and a flow aligned magnetic field in a system rotating with the Sun, are retained. A flow aligned magnetic field in the rotating system may be expressed in terms of the flow velocity and density. Rotational symmetry fixes the longitudinal flow velocity V_φ in terms of the flow in the r?θ plane. Thus, the original three dimensional magnetohydrodynamic flow problem is reduced to a two dimensional hydrodynamic flow problem in the r?θ plane.There are three critical surfaces associated with the equations which supply conditions to determine three of six required boundary conditions. The specified boundary conditions at the base of the corona are the temperature, density, and magnitude of the magnetic field. The equations are then expanded about the radial, nonrotating Parker solution and an analytic solution is obtained for the resulting first order equations. The results show that for constant coronal boundary conditions there is a latitudinal flow toward the solar poles, as a result of magnetic stresses, which persists out to large distances for the Sun. Associated with this flow is a latitudinal component of the magnetic field. The radial flow parameters are, to within small first order differences, in agreement with those of the Parker and the Weber-Davis models of the solar wind.The equations are further generalized to permit first order latitudinal variations in the specified coronal boundary conditions. Results at 1 a.u. are presented for 5 per cent latitudinal differences between the equatorial and polar values. These results show that the solution at 1 a.u. is most sensitive to a latitudinal dependence in the boundary temperature and least sensitive to a latitudinal dependence in the magnetic field magnitude.A solution is then obtained for an approximate dipolar variation in the coronal magnetic field magnitude. This solution predicts that the latitudinal flow is initially toward the Equator due to magnetic channeling; however, this effect is rapidly overcome and the latitudinal flow at 1 a.u. is toward the pole and not significantly different from the solution for constant boundary conditions.     

Interplanetary Signatures of Unipolar Streamers and the Origin of the Slow Solar Wind

Unipolar streamers (also known as pseudo-streamers) are coronal structures that, at least in coronagraph images, and when viewed at the correct orientation, are often indistinguishable from dipolar (or “standard”) streamers. When interpreted with the aid of a coronal magnetic field model, however, they are shown to consist of a pair of loop arcades. Whereas dipolar streamers separate coronal holes of the opposite polarity and whose cusp is the origin of the heliospheric current sheet, unipolar streamers separate coronal holes of the same polarity and are therefore not associated with a current sheet. In this study, we investigate the interplanetary signatures of unipolar streamers. Using a global MHD model of the solar corona driven by the observed photospheric magnetic field for Carrington rotation 2060, we map the ACE trajectory back to the Sun. The results suggest that ACE fortuitously traversed through a large and well-defined unipolar streamer. We also compare heliospheric model results at 1 AU with ACE in-situ measurements for Carrington rotation 2060. The results strongly suggest that the solar wind associated with unipolar streamers is slow. We also compare predictions using the original Wang–Sheeley (WS) empirically determined inverse relationship between solar wind speed and expansion factor. Because of the very low expansion factors associated with unipolar streamers, the WS model predicts high speeds, in disagreement with the observations. We discuss the implications of these results in terms of theories for the origin of the slow solar wind. Specifically, premises relying on the expansion factor of coronal flux tubes to modulate the properties of the plasma (and speed, in particular) must address the issue that while the coronal expansion factors are significantly different at dipolar and unipolar streamers, the properties of the measured solar wind are, at least qualitatively, very similar.     

Solution of one-fluid model equations with short range retarding magnetic forces for the quiet solar wind

The discrepancy of the low predicted versus the observed coronal particle densities is investigated by considering radial magnetic forces acting at the base of the corona in the one fluid model equations with anomalous thermal conductivity for the quiet solar wind. If the short range retarding magnetic force is taken to fall asr ^?5,r being the heliocentric distance, then in order to obtain satisfactory agreement between the predicted and observed (about 3×10⁸ cm^?3 at 1R _⊙) coronal densities, the strength of the retarding magnetic force at 1R _⊙ should be 1.2 times that of the gravitational force.     

Solar wind interaction with interstellar helium

Solar wind interaction with neutral interstellar helium focused by the Sun's gravity in the downwind solar cavity is discussed in a hydrodynamical approach. Upon ionization the helium atoms “picked up” by the (single fluid) solar wind plasma cause a slight decrease in the wind speed and a corresponding marked temperature increase. For neutral helium density outside the cavity n_He^∞ = 0.01 atoms cm^?3 and for interstellar kinetic temperature T_He= 10,000 K, the reduction is speed of the solar wind on the downwind axis at 10 AU from the Sun amounts to about 2kms^?1; the solar wind temperature excess attains 7000 K. The resulting pressure excess leads to a non-radial flow of the order of 0.25 km s^?1. The possibility of experimental confirmation is discussed.     

Magnetohydrodynamic tube waves and high speed solar wind streams

It has been widely conjectured that magnetohydrodynamic (MHD) waves may provide the extra momentum or energy required to explain the high speed solar wind streams that originate in coronal holes. Although the magnetic structuring inherent in this problem has been incorporated into models of the bulk flow, this is not generally true of the associated treatments of wave propagation. In particular, as pointed out by Davila (1985), we might generally expect the magnetic geometry to substantially modify those waves whose wavelength is comparable to the hole width. Using both a geometrical optics and an eigenmode approach, Davila addressed the question of wave propagation in a simple uniform width flux slab model of a coronal hole and concluded

1. the hole may act as a ‘leaky wave guide’, i.e., waves travelling along it may leak into the surrounding corona, but
2. the group velocity of waves with periods in a physically relevant range (around 100 s) is downward, indicating that such waves cannot carry energy into the solar wind and therefore cannot be driving it.

We agree with (i) but argue that (ii) results from a mistaken interpretation of a dispersion relation, and is incorrect. Furthermore, we apply the cylindrical tube leaky wave approach of Cally (1986) to a simple coronal hole model, and find two wavetypes with substantial upward energy fluxes. However, of these, we argue that the so-called ‘trig modes’ (geometry modified fast waves) leak so profusely that they are unable to transport energy over the distance required; the non-axisymmetric ‘thin tube’ modes, though, do not suffer from this disability.     

On the dispersal of artificially-injected gases in the night-time atmosphere

We have studied the extent to which various transport processes affect the dispersal of a gas artificially injected into the night-time atmosphere at F-region altitudes. In addition to diffusion, we have found that nonlinear acceleration, viscous stress, and thermospheric winds affect the dispersal of the injected gas. The magnitude of the effect depends on the atmospheric density, which is a function of solar activity. For an injected H₂ gas, non-linear acceleration and viscous stress rapidly become more important than diffusion above about 300 km for low solar activity (T_∞ = 750K), 340 km for moderate solar activity (T_∞ = 1000K), and 400 km for high solar activity (T_∞ = 1500K). For an injected H₂O gas, the corresponding altitudes are 350, 400, and 470 km for low, moderate and high solar activity, respectively. The effect of nonlinear acceleration and viscous stress is to retard the expansion of the injected gas. Thermospheric winds of 150–400 m s^?1 are important at altitudes near and below the F-region peak electron density. These winds act to transport the injected gas in the wind direction and this affects the shape and temporal development of the subsequent ionospheric hole. Because the H₂O diffusion coefficient is smaller than the H₂ diffusion coefficient, winds are more important for H₂O than for H₂.     

Coronal Temperature Profiles Obtained from Kinetic Models and from Coronal Brightness Measurements Obtained During Solar Eclipses

Coronal density, temperature, and heat-flux distributions for the equatorial and polar corona have been deduced from Saito’s model of averaged coronal white-light (WL) brightness and polarization observations. These distributions are compared with those determined from a kinetic collisionless/exospheric model of the solar corona. This comparison indicates similar distributions at large radial distances (>?7 R_⊙) in the collisionless region. However, rather important differences are found close to the Sun in the acceleration region of the solar wind. The exospheric heat flux is directed away from the Sun, while that inferred from all WL coronal observations is in the opposite direction, i.e. conducting heat from the inner corona toward the chromosphere. This could indicate that the source of coronal heating extends up into the inner corona, where it maximizes at r>1.5 R_⊙, well above the transition region.     

Expansion of the solar wind from a two-hole corona

A one-fluid model is employed to study the global expansion of the solar wind from a two-hole corona, under the assumptions that the holes are confined to polar caps within 30° of heliographic colatitude, the flow is steady and axisymmetric, and the geometry of streamlines is prescribed. The boundary conditions are adjusted in such a way that the calculated solar wind properties at 1 AU are in a reasonable agreement with observational results. A series of numerical solutions are obtained, the series produces a maximum terminal speed of 829 km s^?1 at the pole. The calculated solar wind speeds are strongly latitude dependent and are positively correlated with local divergence factor of a stream tube. The solutions imply that most plasma properties are highly inhomogeneous at the polar caps. The flow velocity, the temperature, the proton number flux and the conduction heat flux all increase towards the hole center.     

Signatures of Slow Solar Wind Streams from Active Regions in the Inner Corona

The identification of solar-wind sources is an important question in solar physics. The existing solar-wind models (e.g., the Wang–Sheeley–Arge model) provide the approximate locations of the solar wind sources based on magnetic field extrapolations. It has been suggested recently that plasma outflows observed at the edges of active regions may be a source of the slow solar wind. To explore this we analyze an isolated active region (AR) adjacent to small coronal hole (CH) in July/August 2009. On 1 August, Hinode/EUV Imaging Spectrometer observations showed two compact outflow regions in the corona. Coronal rays were observed above the active-region coronal hole (ARCH) region on the eastern limb on 31 July by STEREO-A/EUVI and at the western limb on 7 August by CORONAS-Photon/TESIS telescopes. In both cases the coronal rays were co-aligned with open magnetic-field lines given by the potential field source surface model, which expanded into the streamer. The solar-wind parameters measured by STEREO-B, ACE, Wind, and STEREO-A confirmed the identification of the ARCH as a source region of the slow solar wind. The results of the study support the suggestion that coronal rays can represent signatures of outflows from ARs propagating in the inner corona along open field lines into the heliosphere.