Issue 
A&A
Volume 592, August 2016



Article Number  A65  
Number of page(s)  11  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201628599  
Published online  26 July 2016 
Fluxtube geometry and solar wind speed during an activity cycle
^{1} Université de Toulouse, UPSOMP, IRAP, 31400 Toulouse, France
email: rui.pinto@irap.omp.eu
^{2} CNRS, IRAP, 9 Av. colonel Roche, BP 44346, 31028 Toulouse Cedex 4, France
^{3} Laboratoire AIM ParisSaclay, CEA/Irfu Université ParisDiderot CNRS/ INSU, 91191 GifsurYvette, France
Received: 28 March 2016
Accepted: 23 May 2016
Context. The solar wind speed at 1 AU shows cyclic variations in latitude and in time which reflect the evolution of the global background magnetic field during the activity cycle. It is commonly accepted that the terminal (asymptotic) wind speed in a given magnetic fluxtube is generally anticorrelated with its total expansion ratio, which motivated the definition of widely used semiempirical scaling laws relating one to the other. In practice, such scaling laws require ad hoc corrections (especially for the slow wind in the vicinities of streamer/coronal hole boundaries) and empirical fits to in situ spacecraft data. A predictive law based solely on physical principles is still missing.
Aims. We test whether the fluxtube expansion is the controlling factor of the wind speed at all phases of the cycle and at all latitudes (close to and far from streamer boundaries) using a very large sample of windcarrying open magnetic fluxtubes. We furthermore search for additional physical parameters based on the geometry of the coronal magnetic field which have an influence on the terminal wind flow speed.
Methods. We use numerical magnetohydrodynamical simulations of the corona and wind coupled to a dynamo model to determine the properties of the coronal magnetic field and of the wind velocity (as a function of time and latitude) during a whole 11yr activity cycle. These simulations provide a large statistical ensemble of open fluxtubes which we analyse conjointly in order to identify relations of dependence between the wind speed and geometrical parameters of the fluxtubes which are valid globally (for all latitudes and moments of the cycle).
Results. Our study confirms that the terminal (asymptotic) speed of the solar wind depends very strongly on the geometry of the open magnetic fluxtubes through which it flows. The total fluxtube expansion is more clearly anticorrelated with the wind speed for fast rather than for slow wind flows, and effectively controls the locations of these flows during solar minima. Overall, the actual asymptotic wind speeds attained – especially those of the slow wind – are also strongly dependent on fieldline inclination and magnetic field amplitude at the footpoints. We suggest ways of including these parameters in future predictive scaling laws for the solar wind speed.
Key words: Sun: corona / solar wind / Sun: magnetic fields
© ESO, 2016
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
The solar wind flow exhibits a largescale distribution of fast and slow wind streams which evolves markedly during the solar activity cycle. During solar minima, slow wind flows are essentially confined to a region between 20 deg and −20 deg around the ecliptic plane and emanate from the vicinity of the streamer/coronal hole boundaries (S/CH), while fast wind flows stream out from the polar coronal holes and fill all the polar and midlatitude regions faraway from the Sun. During solar maxima, these two wind components get mixed in latitude as a consequence of streamers (and pseudostreamers) appearing at high solar latitudes and of coronal holes making incursions into the lowlatitude regions (McComas et al. 2003). It is clear that the spatial (latitudinal) distribution of slow and fast wind flows follows closely the cyclic variations of the underlying coronal magnetic field structure, itself a consequence of the 11yr cycle of the solar dynamo (McComas et al. 2008; Smith 2011; Richardson & Kasper 2008). This, together with the notion that the properties of the surface motions (assumed as energy sources for the heating and acceleration of the wind) are much more uniform across the solar disk than the amplitude of the wind flows above, suggests that the coronal environment causes the segregation between fast and slow solar wind flows. In particular, the wind terminal speeds seem to be determined to a great extent by the geometrical properties of the magnetic fluxtubes through which the solar wind flows (Wang & Sheeley 1990). Theories that predict the solar wind speed often make use of simple parameters describing the variations of the crosssections of the fluxtubes as a function of height, namely the expansion factor (1)where A_{0} and A_{1} are the crosssection of a given elemental fluxtube respectively at the surface of the Sun (r = r_{0}) and at some point higher in the corona (r = r_{1}>r_{0}) above which the fluxtubes expand radially outwards (and not superradially). A radially expanding fluxtube has a total expansion ratio f = 1, while a very strongly expanding fluxtube has f ≫ 1. The use of potential field extrapolations with sourcesurface (PFSS) from magnetogram data lead to associating r_{1} with the radius of the sourcesurface, commonly placed at a fixed height of r_{SS} = 2.5 R_{⊙}, and to setting the expansion factor f ≡ f_{SS} with respect to this height (Wang & Sheeley 1990). This value for the sourcesurface radius was determined to be the one which produced the best match between the geometry of the extrapolated magnetic fields and the shapes of the coronal structures observed in whitelight during solar eclipses, especially the size of the streamers and coronal hole boundaries (Wang et al. 2010; Wang 2009). However, matching quantities such as the open magnetic flux requires defining r_{SS} as a function of the solar activity (Lee et al. 2011; Arden et al. 2014), or more generally as a function of the properties of the global coronal magnetic field (Réville et al. 2015).
Fig. 1
Structure of the coronal magnetic field at five representative instants of the modelled solar cycle. The times are, from left to right, t = 0, 2.75 ,3.30, 3.80, 4.40 yr corresponding, respectively, to the activity minimum, the early rising phase, the late rising phase, the maximum, and the declining phase of the cycle. The field lines are coloured according to the terminal speed of the wind (at the outer boundary of the domain, r = 15 R_{⊙}) flowing along each one (dark red and orange represent fast wind flows, and blue slow wind flows; cf. the colourbar in Fig. 2). The grey lines represent closed field lines, for completeness. 
Suzuki (2006) also suggested that the terminal wind speed would be better predicted by a combination of the expansion factor and the magnetic field amplitude at the footpoint of any given fluxtube or, equivalently, to the open magnetic flux. Other authors also invoke empirically derived parameters such as the angular distance from the footpoint of a given magnetic fluxtube to the nearest streamer/coronalhole (S/CH) boundary (parameter θ_{b}; Arge et al. 2003, 2004; McGregor et al. 2011). Recent studies by Li et al. (2011) and Lionello et al. (2014) also brought forward that fieldline curvature may have an impact on the wind speed based on solar wind simulations constrained by idealised magnetic geometries. Analysis of Ulysses data combined with PFSS extrapolations by Peleikis et al. (2016) also suggests that fieldline bending correlates to some extent with terminal wind speed. Here, we analyse in greater detail the distributions of terminal wind speed with respect to the geometrical properties of a very large set of individual fluxtubes. We use the numerical method described by Pinto et al. (2011) to couple a solar dynamo to the corona and the solar wind in order to have access to a large sample of solar wind profiles representative of all solar latitudes and moments of the activity cycle.
2. Methods
We use globalscale magnetohydrodynamical (MHD) simulations of the solar dynamo, corona and wind to investigate how the solar wind speed relates with fluxtube geometry at all latitudes and moments of the solar cycle. The model provides us with maps of the timevarying coronal magnetic field and of the wind velocity in the meridional plane from the solar surface up to 15 R_{⊙} during an 11 yr activity cycle. The simulation method, described in detail in Pinto et al. (2011), couples two 2.5D axisymmetric MHD codes. The first − STELEM (Jouve & Brun 2007) − computes the temporal evolution of the surface magnetic field using a meanfield kinematic dynamo (with a BabcockLeighton source term at the surface) driven by the meridional circulation and differential rotation in the convection zone. The second code − DIP (Grappin et al. 2000) − computes the temporal evolution of an MHD solar corona with a selfconsistent wind. The model assumes the corona and the wind to be isothermal, with a uniform coronal temperature T_{0} = 1.3 MK and a specific heat ratio γ = 1 (hence, we do not model the heating and cooling mechanisms in detail). Figure 1 shows some snapshots of the evolution of the coronal magnetic field during the solar cycle in our model. The open magnetic field lines are coloured according to the terminal speed of the wind flow which streams along each one; dark red corresponds to the fastest wind flows (~520 km s^{1} at 15 R_{⊙}) and dark blue to the slowest wind flows (~250 km s^{1} at 15 R_{⊙}). The grey lines are closed field lines with no bulk flow along them. The first instant represented in the figure (in the panel to the left) corresponds to the minimum of activity, and shows a large equatorial streamer surrounded by two wide polar coronal holes. A single equatorial current sheet extends from the top of the large streamer outwards. Fast wind flows well inside the polar coronal holes, hence originating at high latitudes at the surface of the Sun and spreading out to lower latitudes higher up in the corona. Slow wind flows are restricted to the regions nears the boundaries between the equatorial streamer and the polar coronal holes. This general picture is maintained during the rising phase of the cycle (second and third panels of Fig. 1, up to about 3.5 yr), although the slow wind spreads progressively into higher latitudes. The fourth instant represented in Fig. 1 corresponds roughly to the maximum of activity, which occurs just before the global polarity reversal. At this stage, the large equatorial streamer has given place to multiple smaller streamers lying at low latitudes, and new highlatitude streamers have appeared. Several current sheets are now present at different latitudes. Slow wind flows now occupy a much wider range of latitudes, and the overall contrast between the maximum and minimum wind speeds become smaller. A very moderate slow to fast wind gradient from the equator to the poles still persists, nevertheless. The solar wind now originates at multiple and noncontiguous locations over the surface, but each with a latitudinal extent which is much smaller than those of the polar coronal holes present at the minimum of activity. Other important topological features such as pseudostreamers (different from streamers in that they are formed in unipolar regions in coronal holes and are not associated with a heliospheric current sheet) also appear both during the cycle rise and decay phases, as shown respectively in the third and fifth panels of the figure (see also Pinto et al. 2011).
Fig. 2
Timelatitude diagram of the terminal wind speed covering 11 yr of the solar cycle (yaxis) and 180 deg of latitude (from north to south pole, xaxis). The colourscheme is the same as in the field lines plotted in Fig. 1. 
Figure 2 shows the temporal evolution of the latitudinal distribution of fast/slow wind distributions at 15 R_{⊙} during the cycle in a timelatitude diagram. The colour scheme is the same as in Fig. 1, with orange tones representing fast wind flows and blue tones representing slow wind flows. The diagram shows that the fast wind is more prevalent during activity minima, in opposition to the slow wind. The width of the boundaries between slow and fast wind flows varies during the cycle. The transition from slow to fast wind is much sharper during activity minima. This fast wind – slow wind pattern shows good qualitative agreement with the IPS radio maps of Manoharan (2009) and Tokumaru et al. (2010), and with the estimations by Wang & Sheeley (2006) using ULYSSES data and semiempirical methods. Overall, the contrast between the highest and the lowest wind speeds is highest at the minimum of activity and lowest at the maximum. The wind speeds we obtain for the fast wind flows are, however, lower than those measured in the solar wind. Adding additional sources of acceleration known to be efficient in the fast wind regime, such as the ponderomotive force resulting from the propagation and variation of amplitude of Alfvén waves (e.g. Oran et al. 2013; Gressl et al. 2013; van der Holst et al. 2014), would potentially solve this problem. However, we do not want to rely here on physical mechanisms whose inclusion depends on the outcome of the model itself, which is the distribution of fast and slow wind flows.
The goal of this study is to relate the wind speed distributions in Fig. 2 to the geometry of the magnetic features of the low corona shown in Fig. 1. We consider a large sample of magnetic field lines which probe all openfield zones at different instants of the cycle (hence sampling the full time and latitude intervals). The field lines are equally spaced in latitude (Δθ ≈ 1 deg) at the outer boundary of the numerical domain (r = 15 R_{⊙}), where they are rooted, and a new set of field lines is sampled each 0.5 yr. We then extract a wind velocity profile along each line, and perform an ensemble analysis in order to correlate the terminal wind speeds with parameters of the correspondent magnetic fluxtubes and to identify trends valid for all latitudes and moments of the cycle.
We do not make direct comparisons with spacecraft data in this manuscript, as our model does not represent real solar dynamo data. Our solutions are generically representative of a solar cycle but do not mimic the conditions of one specific cycle. However, our simulations let us study the evolution of the wind speed and fluxtube expansion ratios much more coherently than can be achieved by combinations of observational and extrapolation methods, which can be very sensitive to small variations in magnetic connectivity (especially in zones where the magnetic field lines are very strongly divergent), and systematically produce over and underestimations of the expansion ratios (Cohen 2015). Furthermore, in our model the geometry of the closedfield structures (streamers, pseudostreamers) result from the interaction of the wind flow with the coronal magnetic field, and so their heights are not limited to a preestablished sourcesurface height.
3. Results
3.1. Wind speed and expansion
Fig. 3
Terminal wind speed V_{wind} as a function of the total expansion factor f_{tot} for the whole cycle and for all latitudes. 
Fig. 4
Histograms of the terminal wind speed V_{wind} for three different intervals of total expansion factor (1 ≤ f_{tot}< 10 in black, 10 ≤ f_{tot}< 20 in blue, and f_{tot} ≥ 20 in red). The data used covers the whole activity cycle and all latitudes. 
Fig. 5
Total expansion factors (f_{tot}; see Eq. (2)) as a function of latitude at five instants of the cycle (from the top to the bottom panel: t = 0, 2.75 ,3.30, 3.85, 4.40 yr). The field lines and points are coloured as a function of the terminal wind speed (with the same colour table as in Fig. 1). The wind speed is anticorrelated with f_{tot} for a large fraction of the simulated activity cycle, but the are exceptions. The blue and orange peaks (corresponding to slow and fast wind) in the fourth panel (t = 3.85 yr) show that fluxtubes with equal values of f_{tot} can bear wind flows with different speeds. 
Figure 3 shows the dependence of the terminal wind speed V_{wind} on the total expansion factor f_{tot} for all points in the data set (all latitudes and times sampled). The total expansion factor is defined as (2)where B_{0} and B_{1} are evaluated respectively at the surface and at the outer boundary of the domain (15 R_{⊙}), with the same notation for A_{0} and A_{1}. There is, as expected, a general negative correlation between f_{tot} and V_{wind} (cf. Wang & Sheeley 1990), but the terminal wind speed V_{wind} does not seem to be a simple function of f_{tot} alone (see also Woolsey & Cranmer 2014). The spread in the scatterplot is large, especially for the slow wind part of the diagram where the data points are roughly regularly spaced over more than an order of magnitude in f_{tot}, but only by a factor ~2 in V_{wind}. Furthermore, there is a break in the diagram separating the low wind speed/highexpansion part of the diagram from the high wind speed/lowexpansion part.
Figure 4 shows three histograms of the terminal wind speed for three contiguous intervals of the total expansion factor. The black bars correspond to the interval 1 ≤ f_{tot}< 10, the blue bars to 10 ≤ f_{tot}< 20, and the red bars to f_{tot} ≥ 20. The figure shows that the asymptotic wind speeds are indeed inversely correlated with the expansion factor, but the spreads in the distributions are large and, furthermore, the different distributions overlap each other to a great extent. This implies that the parameter f_{tot} alone is not enough to determine the terminal wind speed attained on a given fluxtube.
The discrepancy is more evident when comparing smaller subsets of the simulation data. Figure 5 represents the total fluxtube expansion factors f_{tot} as a function of fieldline latitude (measured at the outer boundary, where the field lines are rooted) for five subsamples of open magnetic field lines, corresponding to the five instants represented in Fig. 1. The plot symbols are coloured as a function of terminal wind speed using the same colourtable (dark red for fast wind flows, dark blue for slow wind flows). During the minimum of activity (first panel), the slow wind flows streams uniquely at the periphery of the coronal holes, while fast wind flows within the wide coronal holes placed at high latitude. The transition from fast to slow wind streams corresponds to a sharp increase in f_{tot}. In this configuration, the usual semiempirical hypothesis relating the wind speed with the geometry of the flow applies (at least qualitatively), with the wind speed being strongly dependent on the fluxtube total expansion ratio and/or proximity to the nearest coronal hole boundary (the wind is indeed slower in the fluxtubes with higher expansion rates and/or closer to a S/CH boundary). However, the situation changes as the activity cycle progress, and becomes very different during the maximum of activity (fourth panel). The slow wind flows cover a much larger latitudinal extent at the solar maximum, there are many more smaller streamers and current sheets spread in latitude, and the transition between fast and slow winds becomes much smoother. Overall, the total expansion ratios are higher than at the minimum (f_{tot}> 20 almost everywhere), the angular distances to S/CH boundaries at the surface (parameter θ_{b}) are much lower (everywhere lower than 1 deg), but the minimum wind speed is actually higher (by a factor ~1.2). Overall, there is one sharp peak in expansion ratio (f_{tot} ≈ 120) related to a wind speed of about 250 km s^{1} at the minimum, and four expansion peaks reaching f_{tot} ≈ 400 corresponding to wind speeds in the range of 350−400 km s^{1} Our simulations hence suggest that the simple parameters f_{tot} (similar to the parameter f_{SS} found in the literature which makes use of PFSS methods, except that here expansion can occur for larger height ranges) and θ_{b} are not sufficient to predict terminal wind speeds accurately.
Fig. 6
Profiles of the expansion factors as a function of the distance to the surface corresponding to fluxtubes at different latitudes at five different instants (from the top to the bottom panel: t = 0, 2.75, 3.30, 3.85, 4.40 yr). The curves are coloured according to the asymptotic wind speed, as in the previous figures. The figures show that the terminal wind speed is anticorrelated with the total expansion factor (the value of at r ≫ 1 R_{⊙}), especially when close to the solar minimum, but does not so much on the details of at low altitudes. 
Let us now examine the expansion profiles of our set of fluxtubes, that is, the way the crosssection of each fluxtube varies with height. Figure 6 represents several expansion profiles (3)where l is the distance to the surface measured along the field line, such that f_{tot} = f(l ≫ R_{⊙}). The curves are colourcoded as those in the preceding figures, with orange representing flows with high terminal speed and blue representing slow wind flows. The figure shows that most of the fluxtubes display smooth and steadily increasing expansion profiles, especially during the minimum of the activity cycle. Conversely, fluxtubes that are in the vicinity of small streamers and pseudostreamers (more common during the rising phase and the maximum of the activity cycle) often display a sharp increase in followed by an abrupt inflexion at low coronal altitudes, which then evolve into a more conventional and smoother overexpansion (cf. Wang et al. 2012). The terminal speeds of the wind flows in these particular fluxtubes is in most cases related to the total expansion ( far from the surface) rather to the maximum values of occurring below the height of the streamers (clearly visible in the third and fifth panels of Fig. 6). In fact, fluxtubes showing such strong inflexions can produce both slow and fast wind flows, as shown in the third and fifth panels of Fig. 6. It would be physically sound to search for an additional parameter describing the radial expansion of the fluxtubes (as the height or height interval over which strong expansion occurs rather than just the total expansion), but our results do not reveal the effect any such parameter clearly and unambiguously (even though these features seem at first to be more often related to fast wind flows when they occur lower down in the corona, and to slower flows when they occur higher up in our simulations).
3.2. Wind speed, expansion, and magnetic field amplitude
Other predictive theories for the solar wind speed combine the magnetic field amplitude at the footpoint of the fluxtubes in addition to their expansion ratios (e.g. Kovalenko 1978, 1981; Wang & Sheeley 1991; Suzuki 2006). This is justified, in part, by considering that the energy input to the solar wind takes the form of a Poynting flux resulting from horizontal surface motions at the surface of the Sun acting on the magnetic field lines which cross it and extend into the corona. The corresponding energy flux density (the Poynting flux) can in that case be expressed as at the surface, where B_{0} is the amplitude of the magnetic field there, v_{⊥} is the amplitude of the transverse motions, and ρ is the mass density. The factor B_{0} is the most significant one in the expression, as long as it is assumed that has the same statistical properties across the whole solar surface. Suzuki (2006), in particular, suggests that the wind speed should scale as (keeping the coronal temperature and the wave power input constant). The amplitude of the magnetic field and the way it decays with height can also influence the rate and position at which this energy input is dissipated in the corona, which in turn affects the terminal wind speed (see Pinto et al. 2009; Cranmer 2002,and many others). Alternatively, one can interpret this relationship between wind speed and B_{0}/f_{tot} as resulting from a conjugation of geometrical properties of the coronal magnetic field. It is easy to see from Eq. (2) that this ratio is an expression of the magnetic field amplitude away from the Sun (or of the open magnetic flux), whose average value is known to vary cyclically, being smaller near the solar minima and larger near the solar maxima. Fujiki et al. (2015) have shown, in fact, that the degree of correlation between B_{0}/f and wind speed itself varies systematically over the solar cycle using a combination of PFSS extrapolations and IPS radio data. They attributed this to the fact that the magnetic field amplitude at the sourcesurface and the average wind speed both vary during the solar cycle, albeit with different relative variations (and possibly different causes).
Fig. 7
Scatterplots of the total fluxtube expansion rate as a function of the magnetic field strength at the surface. The top panel represent these quantities calculated from our simulations for an 11 yr period. The points are colourcoded according to the speed attained by the wind flow in each fluxtube (orange represents fast wind, blue represents slow wind). The bottom panel represents only the simulation data points relative to the minimum and to the maximum (first and third panels in Fig. 5). 
Fig. 8
Terminal wind speed V_{wind} as a function of the ratio B_{0}/f_{tot} for the full data set. The green dashed line corresponds to the curve (B_{0}/f_{tot})^{0.5} proposed by Suzuki (2006) normalised to an arbitrary coefficient for easier representation. The brown dashed and blue continuous lines represent the curves (B_{0}/f_{tot})^{0.1} and (B_{0}/f_{tot})^{2.2}, which are empirical fits to the upper and lower envelopes of the data points (correspondingly, fast and slow wind limits; the curves are also normalised to arbitrary coefficients). 
Figure 7 shows the total fluxtube expansion rate f_{tot} plotted against the footpoint magnetic field amplitudes B_{0} in our simulations. Each point in the scatterplot corresponds to one of the fluxtubes sampled. The points are coloured using the same colour scheme as in Figs. 5 and 6: red corresponds to high terminal wind speeds and blue to low terminal wind speeds. The top panel shows the whole set of fluxtubes, while the bottom panel shows only the data points for instants relative to the minimum and to the maximum (at the same two instants as in Fig. 5). The main features of the full f_{tot}−B_{0} diagram (top panel) are a general positive trend/slope (although with considerable scatter), a welldefined cutoff at the low end of the scatter plot, a sharp transition from fast to slow wind regimes (more clearly visible in the second panel), and a dependence of the wind speed on both quantities. The f_{tot}−B_{0} curves in the bottom panel show that at a given moment of the cycle these two quantities relate to each another and follow a broken power law of the kind , with different values of the index α for the slow and for the fast wind. The fast wind regime corresponds to a moderate α, while the slow wind regime corresponds to a much steeper (and almost undefined) index α. These α indexes are roughly invariant throughout the whole cycle. The strongest deviation occurs when the background magnetic field is at its most multipolar state, at about the maximum of activity. The spread in the f_{tot}−B_{0} diagram increases when all instants of the cycle are represented, but the properties described above are maintained throughout. The “lowend cutoff” of the f_{tot}−B_{0} diagram corresponds to the configuration of the coronal magnetic field at the solar minimum, when it reaches its simpler topological configuration (i.e. quasidipolar). This is the state at which the smallest expansion ratios are attained (f_{tot} ≈ 1, close to the poles). The minimal expansion ratios progressively increase as the cycle proceeds towards the maximum of activity. Equivalently, the latitudinal extent of the coronal holes at the surface decreases from minimum to maximum. The largest deviation with respect to the lowend cutoff line corresponds to the solar maximum. The overall red – blue gradient in Fig. 7 indicates that the terminal wind speed is anticorrelated with f_{tot}; the transition from fast to slow wind is almost horizontal (or orthogonal to the f_{tot} axis), but there also are clear variations of wind speed in the direction of the B_{0} axis (which are particularly visible in the slow wind part of the diagram, plotted in different shades of blue). The figure therefore indicates that these two parameters should indeed be combined in order to predict the solar wind speed.
Figure 8 relates the terminal wind speed V_{wind} with the ratio B_{0}/f_{tot} for all the data points in our simulations. The figure shows three curves which represent the power law V_{wind} ∝ (B_{0}/f_{tot})^{ν} with three distinct indices ν. The green dashed line corresponds to ν = 0.5, as proposed by Suzuki (2006), the brown dashed line to ν = 0.1, and the continuous blue line to ν = 2.2. The blue and orange lines are empirical fits to the upper (fast wind) and lower (slow wind) envelopes of our data points. All three curves are normalised and offset to arbitrary coefficients. It is clear from the figure that we cannot fit a single power law of the kind V_{wind} ∝ (B_{0}/f_{tot})^{ν} to our data set. We could perhaps fit a piecewise power law with at least two indices, one for the slow wind and the other for the fast wind. However, it would still be necessary to specify where the transition between the two should be placed as a function of the activity cycle and/or of the level of complexity of the background coronal field. A small subset of our wind solution (moderately fast winds) seems to follow a curve corresponding to the intermediate index ν = 0.5, even though it requires different normalisation and offset coefficients at each given moment of the cycle. This suggests that additional parameters should be taken into consideration.
Fig. 9
Profiles of the fieldline inclination angle (angle to the vertical direction, in degrees, from −90 to 90) as a function of the distance to the surface for fluxtubes at different latitudes and at five different instants of the cycles (from the top to the bottom panel: t = 0, 2.75 ,3.30, 3.85, 4.40 yr). The diamonds mark the position of the maximum inclination for each fluxtube. The times and colour scheme are the same as in Figs. 5 and 6. 
3.3. Influence of the fieldline curvature
We now investigate the influence of the fieldline curvature and inclination on the wind speed. Large deviations of the fluxtubes from the vertical direction are expected to have several effects on the wind flow. The wind flowing along inclined portions of a fluxtube will see a reduced effective gravity (which can in extreme cases be null or even flip sign on fieldline switchbacks). Fieldline curvature will also increase the flow path length between two given heights (e.g. between the surface of the Sun and 1 AU) and decrease the pressure gradient felt by the wind flow.
Figure 9 shows the inclination of a large sample of field lines as a functions of radius. The inclination angle is the angle between the vertical direction and the direction parallel to the magnetic field at any given point of a field line. The different panels represent different instants of the cycle and orange/blue lines represent fast/slow winds, as in Fig. 6. The diamonds indicate the positions of the maximum inclination of each field line represented. The fieldline inclination profiles vary considerably throughout the activity cycle. The lowlatitude field lines are significantly inclined with respect to the vertical direction from the surface up to more than 8 R_{⊙} during the minimum, but only up to 3−4 R_{⊙} during the maximum. This difference relates to the height of the streamers at those moments (see Fig. 1). Field lines passing close to the pseudostreamer boundaries display the most extreme inclinations (reaching π/ 2 in some cases), as seen in the third and fifth panels of Fig. 9 (the third and fifth panels of Fig. 1 show the geometry of these structures, respectively midlatitude and polar pseudostreamers). In all cases, field lines with large inclinations are systematically related to slower wind speeds than field lines with small deviations to the vertical direction. The only exceptions to this rule are field lines with strongly inclined parts at lower coronal heights, below ~1.5 R_{⊙}. In other words, the effect of fieldline inclination seems to be important at coronal heights at which the wind flow has already started being accelerated – reaching at least 25 to 50% of its terminal speed – rather than at heights at which the wind speed is still negligible. Fluxtubes which have longer inclined portions always have the slowest wind flows (the darker blue lines in Fig. 9).
Figure 10 shows three histograms of the terminal wind speed for three contiguous intervals of the maximum inclination angle α for each of the fluxtubes. The black bars correspond to the interval 0° ≤ α< 15°, the blue bars to 15° ≤ α ≥ 30°, and the red bars to α ≥ 30°. The figure shows that the slowest wind flows are always associated with fluxtubes with large deviations from the vertical direction, and that fieldline inclination is generally anticorrelated with terminal wind speed. However, as in Fig. 4, the spreads in the distributions are large and the different distributions overlap. It must be noted, however, that extreme inclinations (e.g., α ≥ 60°) can be associated with moderate and fast wind flows, as long as the field lines are bent only in the very low corona (see Fig. 9) where the wind speeds are still very small and effect of inclination on the wind speed is minimised.
Fig. 10
Histograms of the terminal wind speed V_{wind} for three different intervals of the maximum fieldline inclination (0° ≤ α< 15° in black, 15° ≤ α< 30° in blue, and α ≥ 30° in red). The data used covers the whole activity cycle and all latitudes. 
Fig. 11
Terminal wind speed V_{wind} as a function of , where L is the fieldline length and α is the maximum inclination of each field line with respect to the vertical direction. 
Figure 11 shows the terminal wind speed v_{wind} plotted against , where the ratio B_{0}/f_{tot} is the same as in Fig. 8, L is the total fieldline length from the surface up to r = 15 R_{⊙}, and α is the maximum inclination angle attained by the fluxtube. This empirically derived expression reduces the scatter along the xaxis in the diagram in Fig. 8, especially for the fast wind part of the diagram. The slow wind part still maintains a large dispersion.
4. Discussion and conclusions
This study confirms that the terminal (asymptotic) speed of the solar wind depends very strongly on the geometry of the open magnetic fluxtubes through which it flows. Our results indicate that the fluxtube inclination should be added to the wellknown and often invoked fluxtube expansion factor on future predictive laws for the solar wind speed. In our simulations, the wind speed is mildly anticorrelated with the total expansion factor (denoted f_{tot} here, to be distinguished from the f_{ss} used in observational studies based on PFSS extrapolations), but in a way which does not allow a general fit to a power law of the kind V_{wind} = a + b·f^{α} (see Fig. 3). Furthermore, the spread in this relation is large, and a given expansion factor can correspond to different values of wind velocity, especially close to solar maximum (see Figs. 4 and 5). While the variations in speed in the fast wind regime (especially during solar minimum) seem to be effectively controlled by the fluxtube expansion, those in the mid to slow wind regimes require the additional deceleration provided by the fieldline inclination in the low to mid corona.
To get an insight on the combined effect of these two geometrical parameters (fluxtube expansion and inclination), let us express the momentum equation of a wind flow along the direction of an arbitrary field line as (4)where l is the coordinate parallel to the magnetic field, ρ is the plasma density, v the wind speed, g = GM_{⊙}/r^{2} is the gravitational acceleration vector, α is the angle between the magnetic field and the vertical direction, P is the gas pressure, and F the combination of any other source terms involved. The inclination of the fluxtube intervenes directly on Eq. (4) by reducing the amplitude of the gravitational term, while the expansion acts by modulating all the fluxes at play (such as the mass flux involving the density ρ and the wind speed v). The relative influence of the expansion and inclination of the fluxtube on the wind acceleration and deceleration becomes more clear if Eq. (4) is rewritten as (5)where M is the wind Mach number, T = T_{0} is the plasma temperature, A is the fluxtube crosssection, and the other symbols have their usual meanings (cf. Wang 1994). The wind is assumed to be isothermal and composed of an ideal and fully ionised hydrogen gas, and other source terms (summarised by the term F in Eq. (4)) were discarded for simplicity. Below the sonic point (where M< 1), moderate inclination angles (0 ≤ α< 90 deg) contribute to decelerating the wind flow by reducing the amplitude of the gravitational term. More extreme situations with angles α> 90 deg (e.g. in fieldline switchbacks) would lead to a local acceleration of the flow (although the net contribution to the terminal wind speed should be negligible because of the deceleration occurring in the turning back of the field line into the upward direction). A growing fluxtube crosssection – with a positive – will also contribute to decelerating the wind flow. The wind flow will be less decelerated (or even accelerated) on a “reconverging” section of a fluxtube. Overall, the amplitude of all these effects will be minimised on sections of a fluxtube where the wind is either very slow (v ≈ 0) or very fast (M^{2} ≫ 1), and will be maximised at wind speeds close to the sound speed. The fluxtube inclination and expansion terms change sign in the supersonic regime, where M> 1, hence reversing their effects on the acceleration/deceleration of the wind flow there.
Fig. 12
The two terms on the righthand side of Eq. (5) as a function of the distance to the surface for the fluxtubes which carry the fastest and the slowest wind flow at three different instants of the cycles (from the top to the bottom panel: t = 0, 3.30, 3.85 yr), as in Fig. 13. The colour scheme for the wind speed is the same as in Figs. 6 and 9. The continuous lines represent the first term (gravitation and inclination), and the dashed lines represent the second term (expansion). The terms switch sign at the sonic point (cf. the denominator in Eq. (5)) at a height varying between 2 and 6 R_{⊙}. 
Fig. 13
Ratio of the absolute value of the two terms on the righthand side of Eq. (5) as a function of the distance to the surface for the fluxtubes which carry the fastest and the slowest wind flow at three different instants of the cycles (the same as in Fig. 12). The black dotted line shows the value of this ratio for a vertically aligned and radially expanding wind flow (a Parker wind) for comparison. 
Figure 12 shows the contributions of the two terms on the righthand side of Eq. (5) for the wind acceleration and deceleration ∂_{l}v, for the fluxtubes with maximum and minimum terminal wind speed, hereafter and A for simplicity. The continuous lines represent the term , and the dashed lines represent the term A. The colours represent the terminal wind speed, as before. As explained above, the effect of both terms is greatest in the vicinity of the sonic point (where M = 1, and both terms switch sign), and smallest close to the surface and at greater coronal heights. The dashed and continuous lines are nearly symmetrical at all times, indicating that one term cannot in general be neglected with respect to the other. Figure 13 represents the ratio of the absolute value of the two terms at the same instants as in Fig. 12. For comparison, this ratio is equal to ^{(}GM_{⊙}m_{p}^{)}/^{(}4k_{b}Tr^{)} for a verticallyaligned and radially expanding fluxtube, which would drive a Parker wind (represented by a black dotted line in the figure). The blue and orange lines (our solutions) approach this limiting case at high coronal altitudes, where fluxtube inclination and overexpansion become negligible. The two terms are indeed of the same order of magnitude at all times, with the being usually larger than A (but rarely by a factor greater than 2). The slowest wind streams are, in general, associated with rather irregular curves for the ; the dips in these curves are due to local increases in fieldline inclination (or decreases in ). Figures 5 and 6 indeed show that the wind speed is generally anticorrelated with the expansion ratios for the fast wind flows close to the minimum of activity, but not as much for the slow wind or during the maximum. Figure 7 shows that in a slow wind regime (plotted in shades of blue) the wind speed can actually vary more strongly with the magnetic field amplitude B_{0} than with f_{tot}. Furthermore, some of the peaks in Fig. 5 (especially in the fourth panel) show that exceptions to the f_{tot} – V_{wind} relation exist even in the fast wind regime. All these deviations, i.e. winds slower than expected for their total expansion ratios, are related to fluxtubes which have strong inclinations in a range of heights r = 2−4 R_{⊙} (cf. Fig. 9). Inclination does not have a significant impact on the wind flows when it occurs below these heights, as the wind flow is still very subsonic there (see Eq. (5)). The slowest of the wind flows found exist in fluxtubes for which the inclined portions are particularly long (dark blue lines in the first panel of Fig. 9).
These results motivate the search for predictive wind speed scaling laws which use a combination of the geometrical parameters discussed here: fluxtube expansion, magnetic field amplitude, and inclination. The plot in Fig. 11 illustrates a first step in this direction, and consists of a factor similar to that suggested by Suzuki (2006) corrected by , where α is the maximum inclination angle for a given fluxtube (marked in Fig. 9 with diamonds), and L its length (a rough measure of the height interval over which the field line is bent). The resulting expression provides a much better fit to our wind solutions, even though some scatter remains in the slow wind part. It is clear that the relative heights of the inclined and overexpanding fractions of the fluxtubes with respect to the height of the sonic point (or rather, the factor v/^{(}1−M^{2}^{)} in Eq. (5)) is a key parameter that regulates the asymptotic wind speeds; unfortunately, the positions of the sonic point are not known in advance (they depend on the wind speed profile) and cannot be used as a predictor.
The analysis carried out in this manuscript is based on results of a series of numerical simulations using an MHD model of an axisymmetric solar corona at fixed temperature T_{0}. The choice of the parameter T_{0} has a direct effect on the wind solutions: a higher coronal temperature will produce a globally faster and denser wind (Parker 1964; Leer & Holzer 1979; Hansteen & Leer 1995; Pinto et al. 2011). However, the model will still correctly reproduce the relative variations of the properties of the wind in latitude (as long as T_{0} remains fixed throughout and is set to a reasonable value, in the range 1−1.5 MK). Fast and slow solar wind streams are known to have different temperature profiles, with slow winds reaching higher temperatures in the low corona than fast winds, but lower temperatures far from the Sun. The effect of these temperature variations could perhaps be simulated in our model by setting a nonuniform coronal T_{0} (varying in latitude and in time), but – as discussed above – the distribution of fast and slow wind flows is the outcome of the model, on which its parameters cannot depend. An indepth investigation of the thermodynamics of the solar wind flows requires a different kind of model (as in Suzuki & Inutsuka 2006; Pinto et al. 2009; Schwadron & McComas 2003; Woolsey & Cranmer 2014), which will be the subject of future work. Nevertheless, the approach we adopt here strongly suggests that the absolute values and spatial distribution of the solar wind speed is mostly controlled by the geometry of the coronal magnetic field (independently of the heating scenario considered).
The coronal magnetic field is obtained via a kinematic meanfield dynamo model, and is not meant to reproduce a specific solar cycle. Instead, it produces many coronal features representative of the evolution of any solar cycle. More importantly, our simulations provide a very large statistical ensemble of open magnetic fluxtubes with wind flows covering a very wide range of fluxtube geometries.
5. Summary
We have investigated how the solar wind speed relates with fluxtube geometry at all latitudes and moments of the solar cycle by means of globalscale MHD simulations of the solar dynamo, corona, and wind. The model generates maps of the slowly varying coronal magnetic field and of the wind velocity in the meridional plane from the solar surface up to 15 R_{⊙} during an 11 yr activity cycle. We analysed a large sample of individual magnetic fluxtubes covering the full latitude and time intervals in order to derive correlations between their geometrical parameters and the resulting wind speeds that remain valid for all latitudes and moments of the cycle. In addition to the total expansion factors f_{tot} and absolute magnetic field amplitudes B_{0}, we found that the wind speed also depends strongly on fluxtube inclination. Future work will focus on improving the empirically fitted scaling law relating the wind speed to these three fluxtube parameters, on testing more sophisticated heating scenarios, and on extending the analysis to sets of magnetic fluxtubes derived from observations.
Acknowledgments
R.F.P. acknowledges funding by the FP7 project #606692 (HELCATS). The dynamo–wind coupling method was developed within the frame of the STARS2 project (www.stars2.eu). A.S.B. acknowledges funding by the ERC grant Solar Predict, INSU/PNST and the CNES Solar Orbiter grant. We are grateful to Y.M. Wang and R. Grappin for enlightening discussions and comments on the manuscript.
References
 Arden, W. M., Norton, A. A., & Sun, X. 2014, J. Geophys. Res. (Space Phys.), 119, 1476 [NASA ADS] [CrossRef] [Google Scholar]
 Arge, C. N., Odstrcil, D., Pizzo, V. J., & Mayer, L. R. 2003, in AIP Conf. Proc., 679, 190 [Google Scholar]
 Arge, C. N., Luhmann, J. G., Odstrcil, D., Schrijver, C. J., & Li, Y. 2004, J. Atm. SolarTerrestrial Phys., 66, 1295 [Google Scholar]
 Cohen, O. 2015, Sol. Phys., 290, 2245 [NASA ADS] [CrossRef] [Google Scholar]
 Cranmer, S. R. 2002, Space Sci. Rev., 101, 229 [Google Scholar]
 Fujiki, K., Tokumaru, M., Iju, T., Hakamada, K., & Kojima, M. 2015, Sol. Phys., 290, 2491 [NASA ADS] [CrossRef] [Google Scholar]
 Grappin, R., Léorat, J., & Buttighoffer, A. 2000, A&A, 362, 342 [NASA ADS] [Google Scholar]
 Gressl, C., Veronig, A. M., Temmer, M., et al. 2013, Sol. Phys. [Google Scholar]
 Hansteen, V. H., & Leer, E. 1995, J. Geophys. Res., 100, 21577 [NASA ADS] [CrossRef] [Google Scholar]
 Jouve, L., & Brun, A. S. 2007, A&A, 474, 239 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kovalenko, V. A. 1978, Geomagnetism and Aeronomy, 18, 769 [NASA ADS] [Google Scholar]
 Kovalenko, V. A. 1981, Sol. Phys., 73, 383 [NASA ADS] [CrossRef] [Google Scholar]
 Lee, C. O., Luhmann, J. G., Hoeksema, J. T., et al. 2011, Sol. Phys., 269, 367 [Google Scholar]
 Leer, E., & Holzer, T. E. 1979, Sol. Phys., 63, 143 [NASA ADS] [CrossRef] [Google Scholar]
 Li, B., Xia, L. D., & Chen, Y. 2011, A&A, 529, A148 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lionello, R., Velli, M., Downs, C., Linker, J. A., & Mikic, Z. 2014, ApJ, 796, 111 [Google Scholar]
 Manoharan, P. K. 2009, in Proc. International Astronomical Union, 5, Solar and Stellar Variability: Impact on Earth and Planets, 356 [Google Scholar]
 McComas, D. J., Elliott, H. A., Schwadron, N. A., et al. 2003, Geophys. Res. Lett., 30, 24 [Google Scholar]
 McComas, D. J., Ebert, R. W., Elliott, H. A., et al. 2008, Geophys. Res. Lett., 35, L18103 [NASA ADS] [CrossRef] [Google Scholar]
 McGregor, S. L., Hughes, W. J., Arge, C. N., Owens, M. J., & Odstrcil, D. 2011, J. Geophys. Res. (Space Phys.), 116, 3101 [Google Scholar]
 Oran, R., van der Holst, B., Landi, E., et al. 2013, ApJ, 778, 176 [Google Scholar]
 Parker, E. N. 1964, ApJ, 139, 72 [NASA ADS] [CrossRef] [Google Scholar]
 Peleikis, T., Kruse, M., Berger, L., Drews, C., & WimmerSchweingruber, R. F. 2016, in AIP Conf. Proc. 1720, 020003 [Google Scholar]
 Pinto, R., Grappin, R., Wang, Y.M., & Léorat, J. 2009, A&A, 497, 537 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pinto, R. F., Brun, A. S., Jouve, L., & Grappin, R. 2011, ApJ, 737, 72 [Google Scholar]
 Réville, V., Brun, A. S., Strugarek, A., et al. 2015, ApJ, 814, 99 [NASA ADS] [CrossRef] [Google Scholar]
 Richardson, J. D., & Kasper, J. C. 2008, J. Atm. SolarTerrestrial Phys., 70, 219 [NASA ADS] [CrossRef] [Google Scholar]
 Schwadron, N. A., & McComas, D. J. 2003, ApJ, 599, 1395 [NASA ADS] [CrossRef] [Google Scholar]
 Smith, E. J. 2011, J. Atm. SolarTerrestrial Phys., 73, 277 [NASA ADS] [CrossRef] [Google Scholar]
 Suzuki, T. K. 2006, ApJ, 640, L75 [NASA ADS] [CrossRef] [Google Scholar]
 Suzuki, T. K., & Inutsuka, S. 2006, Astronomical Herald, 99, 205 [NASA ADS] [Google Scholar]
 Tokumaru, M., Kojima, M., & Fujiki, K. 2010, J. Geophys. Res. (Space Phys.), 115, 04102 [NASA ADS] [CrossRef] [Google Scholar]
 van der Holst, B., Sokolov, I. V., Meng, X., et al. 2014, ApJ, 782, 81 [Google Scholar]
 Wang, Y.M. 1994, ApJ, 437, L67 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, Y.M. 2009, Space Sci. Rev., 144, 383 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, Y.M., & Sheeley, N. R. 1990, ApJ, 355, 726 [Google Scholar]
 Wang, Y.M., & Sheeley, N. R. 1991, ApJ, 372, L45 [Google Scholar]
 Wang, Y.M., & Sheeley, N. R. 2006, ApJ, 653, 708 [Google Scholar]
 Wang, Y.M., Robbrecht, E., Rouillard, A. P., Sheeley, Jr., N. R., & Thernisien, A. F. R. 2010, ApJ, 715, 39 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, Y.M., Grappin, R., Robbrecht, E., & Sheeley, N. R. 2012, ApJ, 749, 182 [NASA ADS] [CrossRef] [Google Scholar]
 Woolsey, L. N., & Cranmer, S. R. 2014, ApJ, 787, 160 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1
Structure of the coronal magnetic field at five representative instants of the modelled solar cycle. The times are, from left to right, t = 0, 2.75 ,3.30, 3.80, 4.40 yr corresponding, respectively, to the activity minimum, the early rising phase, the late rising phase, the maximum, and the declining phase of the cycle. The field lines are coloured according to the terminal speed of the wind (at the outer boundary of the domain, r = 15 R_{⊙}) flowing along each one (dark red and orange represent fast wind flows, and blue slow wind flows; cf. the colourbar in Fig. 2). The grey lines represent closed field lines, for completeness. 

In the text 
Fig. 2
Timelatitude diagram of the terminal wind speed covering 11 yr of the solar cycle (yaxis) and 180 deg of latitude (from north to south pole, xaxis). The colourscheme is the same as in the field lines plotted in Fig. 1. 

In the text 
Fig. 3
Terminal wind speed V_{wind} as a function of the total expansion factor f_{tot} for the whole cycle and for all latitudes. 

In the text 
Fig. 4
Histograms of the terminal wind speed V_{wind} for three different intervals of total expansion factor (1 ≤ f_{tot}< 10 in black, 10 ≤ f_{tot}< 20 in blue, and f_{tot} ≥ 20 in red). The data used covers the whole activity cycle and all latitudes. 

In the text 
Fig. 5
Total expansion factors (f_{tot}; see Eq. (2)) as a function of latitude at five instants of the cycle (from the top to the bottom panel: t = 0, 2.75 ,3.30, 3.85, 4.40 yr). The field lines and points are coloured as a function of the terminal wind speed (with the same colour table as in Fig. 1). The wind speed is anticorrelated with f_{tot} for a large fraction of the simulated activity cycle, but the are exceptions. The blue and orange peaks (corresponding to slow and fast wind) in the fourth panel (t = 3.85 yr) show that fluxtubes with equal values of f_{tot} can bear wind flows with different speeds. 

In the text 
Fig. 6
Profiles of the expansion factors as a function of the distance to the surface corresponding to fluxtubes at different latitudes at five different instants (from the top to the bottom panel: t = 0, 2.75, 3.30, 3.85, 4.40 yr). The curves are coloured according to the asymptotic wind speed, as in the previous figures. The figures show that the terminal wind speed is anticorrelated with the total expansion factor (the value of at r ≫ 1 R_{⊙}), especially when close to the solar minimum, but does not so much on the details of at low altitudes. 

In the text 
Fig. 7
Scatterplots of the total fluxtube expansion rate as a function of the magnetic field strength at the surface. The top panel represent these quantities calculated from our simulations for an 11 yr period. The points are colourcoded according to the speed attained by the wind flow in each fluxtube (orange represents fast wind, blue represents slow wind). The bottom panel represents only the simulation data points relative to the minimum and to the maximum (first and third panels in Fig. 5). 

In the text 
Fig. 8
Terminal wind speed V_{wind} as a function of the ratio B_{0}/f_{tot} for the full data set. The green dashed line corresponds to the curve (B_{0}/f_{tot})^{0.5} proposed by Suzuki (2006) normalised to an arbitrary coefficient for easier representation. The brown dashed and blue continuous lines represent the curves (B_{0}/f_{tot})^{0.1} and (B_{0}/f_{tot})^{2.2}, which are empirical fits to the upper and lower envelopes of the data points (correspondingly, fast and slow wind limits; the curves are also normalised to arbitrary coefficients). 

In the text 
Fig. 9
Profiles of the fieldline inclination angle (angle to the vertical direction, in degrees, from −90 to 90) as a function of the distance to the surface for fluxtubes at different latitudes and at five different instants of the cycles (from the top to the bottom panel: t = 0, 2.75 ,3.30, 3.85, 4.40 yr). The diamonds mark the position of the maximum inclination for each fluxtube. The times and colour scheme are the same as in Figs. 5 and 6. 

In the text 
Fig. 10
Histograms of the terminal wind speed V_{wind} for three different intervals of the maximum fieldline inclination (0° ≤ α< 15° in black, 15° ≤ α< 30° in blue, and α ≥ 30° in red). The data used covers the whole activity cycle and all latitudes. 

In the text 
Fig. 11
Terminal wind speed V_{wind} as a function of , where L is the fieldline length and α is the maximum inclination of each field line with respect to the vertical direction. 

In the text 
Fig. 12
The two terms on the righthand side of Eq. (5) as a function of the distance to the surface for the fluxtubes which carry the fastest and the slowest wind flow at three different instants of the cycles (from the top to the bottom panel: t = 0, 3.30, 3.85 yr), as in Fig. 13. The colour scheme for the wind speed is the same as in Figs. 6 and 9. The continuous lines represent the first term (gravitation and inclination), and the dashed lines represent the second term (expansion). The terms switch sign at the sonic point (cf. the denominator in Eq. (5)) at a height varying between 2 and 6 R_{⊙}. 

In the text 
Fig. 13
Ratio of the absolute value of the two terms on the righthand side of Eq. (5) as a function of the distance to the surface for the fluxtubes which carry the fastest and the slowest wind flow at three different instants of the cycles (the same as in Fig. 12). The black dotted line shows the value of this ratio for a vertically aligned and radially expanding wind flow (a Parker wind) for comparison. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.