next up previous
Next: North-South Asymmetries Up: Electron Density and Core Previous: Electron Density and Core

Radial Gradients at High Latitudes

The two sets of solar wind data obtained in the steady high-latitude regions determined above, (> tex2html_wrap_inline1026 south and north, respectively), a priori depend on the heliocentric distance and latitude of Ulysses and possibly on temporal variations in solar activity. Nevertheless, if we expect that the expansion of the solar wind in these regions has a spherical symmetry and has no large-scale temporal or latitudinal variations in the observing time range (of about 4 months each), we should find that the density decreases as a power law: tex2html_wrap_inline1166 since the radial variation in bulk velocity is very small in high-speed wind (see, e.g., [Arya & Freeman, 1991]). To analyze such a basic trend of the density variation with heliocentric distance, we used a linear regression method using the two full density data sets in logarithmic scales as shown in the lower part of Figures 5a and 5b. As compared with a linear regression on binned data, this method has the advantage of avoiding an arbitrary choice of the bin size and locations, while averaging as well short term variations of the data.

   figure116
Figure 5: Radial variations of the electron density and core temperature at high latitudes in both solar hemispheres. The solid lines are the best fit power laws to the data whose parameters are shown in Table 1. The vertical bars are the standard deviations tex2html_wrap_inline1028 and tex2html_wrap_inline1030 , determined by the specific fitting used in QTN spectroscopy, on the electron density and core temperature respectively (see section 3.1). Note that southward of tex2html_wrap_inline1026 there are about 53,000 data points, whereas there are about 46,000 data points northward of tex2html_wrap_inline1026 .

Let us note that such a method can be performed because our data sets have no significant measurement gaps and have not been selected on any arbitrary criterion among the raw data.

We processed the linear regression by minimizing the corresponding tex2html_wrap_inline1176 merit function defined by

eqnarray122

with tex2html_wrap_inline1178 , where tex2html_wrap_inline1028 is the uncertainty on each electron density logarithm provided by the specific measurement method of QTN spectroscopy reviewed in section 2. Indeed, we do not compute the tex2html_wrap_inline1176 with tex2html_wrap_inline1184 because we need to take into account the actual short scale variations of the density (with respect to the observing interval), associated for instance to the solar rotation; these variations are represented by the standard deviation of the density (shown in Figure 4) which is about three times larger than the measurement uncertainties tex2html_wrap_inline1028 .

This method thus provides the slope - tex2html_wrap_inline1188 , which is the power law index of the variation of tex2html_wrap_inline1190 , and tex2html_wrap_inline1192 , where tex2html_wrap_inline1194 is the density extrapolated to 1 AU; it also provides uncertainties tex2html_wrap_inline1196 and tex2html_wrap_inline1198 depending on the tex2html_wrap_inline1176 variation ( tex2html_wrap_inline1202 ) around its minimum. Indeed, the tex2html_wrap_inline1202 analysis shows a strongly anticorrelated error-ellipse (not shown) between the intercept and the slope of the log-log fit; in order to give uncertainties with a maximum confidence level, we have derived tex2html_wrap_inline1196 and tex2html_wrap_inline1198 by setting the tex2html_wrap_inline1202 to 18.4, which corresponds to 99.99% of the probability distribution for two fitted parameters varying jointly [Press et al., 1992]

The result for the southern hemisphere data is very close to that we expected: tex2html_wrap_inline1212 ; it is thus compatible with our initial assumption of independence on latitude or time. In contrast, the power law index obtained in the northern hemisphere, tex2html_wrap_inline1214 , suggests that this northern density is dependent on latitude and/or time.

We also processed the ditto linear regression method to our temperature measurements as shown in the top part of Figures 5a and 5b, in order to estimate the core electron temperature variation with heliocentric distance in both hemispheres as a power law: tex2html_wrap_inline1216 . In this case, the standard deviation of the core temperature is of the same order as the measurements uncertainties tex2html_wrap_inline1030 , which are thus used to evaluate the uncertainty on the slope tex2html_wrap_inline1220 in the same manner as above; i.e., tex2html_wrap_inline1222 . Southward of tex2html_wrap_inline1026 we obtain a power index tex2html_wrap_inline1226 , while tex2html_wrap_inline1228 0.83 northward of tex2html_wrap_inline1026 , showing an asymmetry already noted for the density index values (see Table 1 for a summary of the results). Since we can consider from our linear regression on the density that the southern data set is a good sample of stationary high-speed solar wind in spherical expansion, we conclude that the core electron temperature variation with heliocentric distance in such a steady wind can be described by the power law tex2html_wrap_inline1232 between 1.52 and 2.31 AU. This result is obtained with tex2html_wrap_inline1234 and N=52,644 data points, which gives a goodness-of-fit probability very close to 1; (such a good value suggests that the measurement errors might have been slightly overestimated). Hereafter, we will use this law to scale our measured core electron temperatures, and the tex2html_wrap_inline1048 law to scale the densities.


next up previous
Next: North-South Asymmetries Up: Electron Density and Core Previous: Electron Density and Core

Karine Issautier
Fri Nov 27 18:47:01 MET 1998