Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.izmiran.rssi.ru/~obridko/papers/Obridko_n.pdf
Äàòà èçìåíåíèÿ: Fri Feb 28 17:03:36 2003
Äàòà èíäåêñèðîâàíèÿ: Sun Apr 10 01:05:57 2016
Êîäèðîâêà: IBM-866
Astronomy Reports, Vol. 47, No. 4, 2003, pp. 333í342. Translated from Astronomicheski Zhurnal, Vol. 80, No. 4, 2003, pp. 364í373. i Original Russian Text Copyright c 2003 by Ibragimov, Sule anov, Vikhlinin, Sakhibul lin. im

Meridional Drift of Large-Scale Solar Magnetic Fields
V. N. Obridko and B. D. Shelting
Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation, Troitsk, 142190, Russia
Received August 5, 2002; revised October 10, 2002

Abstract--It is shown that the meridional drift of large-scale fields starts in the equatorial zone and continues over 15í16 yrs (16í17 according to another estimate), i.e., during three fourths of the 22-year cycle. There is an abrupt retardation of the drift at latitudes of 30í50 , and a stagnation region arises, where the drift rate does not exceed several meters per second. The drift becomes rapid again at higher latitudes. The stagnation region coincides with the area in which the radial gradient of the rotational velocity is close to zero in the convective zone. This drift is compared with helio-seismological data on the rotation in the convective zone. A model is proposed taking into account some elements of dynamo theory. c 2003 MAIK "Nauka/Interperiodica".

1. INTRODUCTION The phenomenon of meridional drift, or more precisely, the displacement of zones of solar activity in latitude, was revealed virtually simultaneously with the discovery of solar cyclicity. While we recognize Schwabe's report of 1843 as the first indication of a possible 10-year periodicity, we can find the first statement about the drift of a zone of spots in Carrington's classical treatise of 1863. For a long time, it was considered to be obvious that both phenomena exhibit a single periodicity. We now know that, though the spots indeed appear at a fixed latitude with a period of 11 yrs (or closer to 10 yrs in the twentieth century), the lifetime for the spots of each cycle is significantly longer. The first spots of a new cycle appear at latitudes of 40 and higher several years before the solar minimum, while the last spots are sometimes observed after the solar minimum at latitudes of about 5í10 . Thus, the idea of an extended solar cycle was introduced [1]. The next step, namely the analysis of cyclic variations of the large-scale and background fields, was started only in the 1960s. There is no doubt that the periodicity of the large-scale fields at each latitude is close to the periodicity of local fields and sunspots. However, the latitude dependence is less clear. Some papers [2í14] indicate that the meridional drift is directed from the equator toward the poles, however, other papers [15í17] contest this. Finally, some papers of the former group (see, for example, [7]) admit a possible additional drift from middle latitudes toward the equator, along with the meridional drift from the middle latitudes toward the poles. The inconsistency of these findings is first and foremost connected with the difficulties of eliminating the effects of stronger local fields in the total signal. In addition, the solar

rotation, whose velocity exceeds that of the meridional drift of the large-scale fields by two orders of magnitude, contributes greatly to Doppler measurements. The comparatively short interval for which data on the magnetic fields of the entire solar disk are available represents another difficulty. Therefore, the very thorough works [18í20] analyze data for only 10 yrs, while [9, 21] consider data for 22 yrs. Here, we study the meridional drift of the largescale fields over a long time interval, paying particular attention to the identification of the large-scale fields. For our analysis, we use both magnetic-field data detected using various magnetographs over 40 yrs and unified into a single data set, and magneticfield data recovered from H filament observations obtained over 85 yrs. Since we are analyzing cyclic properties of the large-scale magnetic field, lower limits on the spatial size of the field are imposed by the observational technique used. In our analysis, the magnetograph data limit the spatial size of the field in three ways. (1) The spatial resolution of the Stanford magnetograph is 3 ( 0.2R0 , where R0 is the solar radius). Though other magnetographs provide higher resolution, the data were later averaged. (2) The observations of magnetic fields at various latitudes were averaged over an entire Carrington revolution before the construction of the latitudeítime diagrams. (3) The number of harmonic coefficients used when calculating the magnetic-field components from the longitudinal field using Legendre polynomials was limited. We took into account only the first ten harmonics, which corresponds to sizes no smaller

1063-7729/03/4704-0333$24.00 c 2003 MAIK "Nauka/Interperiodica"


334
() 60 40 20 0 í 20 í 40 í 60 (b) 60 40 20 0 í 20 í 40 í 60 1960 1970 1980 Years 1990

OBRIDKO, SHELTING
0í500 í 500í0

2000

Fig. 1. (a) Latitudeítime diagram for the radial magnetic field directly detected by magnetographs. (b) The same diagram combined with the Maunder butterflies.

than 0.2R0 . This provides an additional limitation on the spatial size of the field. For the data, the second and third constraints apply in the same way. In addition, the technique used to recover the field from spectral heliograms reconstructs large-scale fields far from active regions. Since the field strengths are not used in the analysis, the contribution of strong local fields to the reconstructed fields is considerably diminished. Therefore, the reconstructed fields are the most convenient for studies of the meridional drift [22, 23]. In Section 2, we shall compare data on the localization of spots and on the drift of the large-scale magnetic field obtained from direct magnetograph measurements over 40 yrs, as well as with magnetic data reconstructed over 85 yrs. These data are subject to additional filtration in Section 3. A subsidiary analysis to ascertain the phase dependence and drift rate is presented in Section 4. Section 5 presents a discussion of the results. 2. COMPARING THE DRIFTS OF LOCAL AND LARGE-SCALE FIELDS The first set of data obtained by magnetographs covers 40 yrs (1960í2002), and includes the Mount

Wilson magnetic measurements over 1960í1978, Kitt Peak data over 1975í1984, and Stanford measurements from 1976 until the present time. All data were graciously made available by researchers of these observatories. We have unified these data into a single (WSO) system [24]. Using these data, we have constructed latitudeítime diagrams for the longitudinal magnetic fields averaged over a Carrington revolution. Guided by the widespread hypothesis that the large-scale field in the photosphere is radial, we can calculate the radial field and construct the corresponding diagram. The difference between the two diagrams is insignificant, except for a small strengthening of Br in comparison with the longitudinal field in near-polar regions. Therefore, we present only the Br diagram (Fig. 1a). We can immediately see the fine structure of the field: a number of "tongues" and "bays", and, first and foremost, the seeming absence of a clear direction for the meridional drift, especially at low latitudes. The obvious reason is the absence of the third constraint above limiting the size of the fields from below. We observe the total field rather than the sum of its largest-scale components (l 9). The contribution of intense, small-scale, local fields is strong. This is first manifest at low latitudes, where we find a superposition of two meridional drifts: from the equator toward the pole for the large-scale fields and from the middle latitudes toward the equator for the local fields. Figure 1b combines in a single latitudeítime diagram the magnetograph data on the large-scale magnetic fields and the Maunder butterflies for sunspots. We can see that the local fields visibly disrupt the continuous drift of the large-scale fields from the equator toward the pole. The data on the large-scale fields covers 85 yrs (1915í2000). These data are partially taken from the atlas of P. McIntosh (1964í1974), while the remaining are published or presented to us directly by Sivaraman and Makarov (1975í2000). As is discussed above, when working with the data, we used an original technique to derive the magneticfield strength from the synoptic maps, which specify only the sign of the magnetic field. To ensure the values were reasonable, we compared directly measured magnetic-field strengths with those calculated from the observations for overlapping time intervals [24]. A direct study of the polarity of the large-scale fields recovered from data (i.e. without computing the radial field via polynomials) was carried out in [12]. As expected, the meridional drift toward the pole became much clearer in comparison with the magnetograph data. The drift toward the equator became almost invisible, being manifest only in small "tongues"
ASTRONOMY REPORTS Vol. 47 No. 4 2003


MERIDIONAL DRIFT OF LARGE-SCALE SOLAR MAGNETIC FIELDS

335

() 80 60 40

20 0 í 20 í 40 í 60 í80 1960 80 60 40 20 0 í 20 í 40 í 60 í80 1920 1940 1960 Years 1980 2000 1.37¤3 í 2.75¤3 0 í 1.37¤3 í1.37¤3í0 í 2.75.¤ í í1.37¤3 1970 1980 (b) 1990 2000 250 í 500 0 í 250 í250í0 í 500 í í250

Fig. 2. Latitudeítime diagram for the radial magnetic field calculated in the potential approximation from the (a) magnetograph and (b) H data.

of the secondary polarity near the equator. Nevertheless, more careful filtering is obviously necessary. 3. CALCULATIONS OF CYCLIC VARIATIONS OF RADIAL FIELDS IN THE PHOTOSPHERE The calculations of each component of the magnetic field in the photosphere from the detected longiASTRONOMY REPORTS Vol. 47 No. 4 2003

tudinal field are carried out in a potential approximation. The harmonic coefficients are obtained by two different techniques: the method of orthogonal Legendre polynomials [25] and the least squares method [26]. The radius of the source surface is assumed to be 2.5R0 . A "classical" approximation was used; i.e., radial structure for the field was not assumed a priori. Since we were mainly interested in the behavior of


336

OBRIDKO, SHELTING

80 60 40 20 0 í 20 í 40 í 60 í 80 í 100 í 50 0 50 100 í 200 í 100 0 100 200

Correlation shift
Fig. 3. Correlation between the large-scale magnetic fields obtained for various latitudes from the magnetograph (left) and H (right) data.

the radial component Br of the magnetic field, we first calculated the daily radial field values for latitudes from 0 to ‘90 in steps of 1/15 in the sine of the latitude for the whole time interval. The results were then averaged over a Carrington revolution. We subsequently combined all data into a single file and constructed the latitudeítime diagrams. Let us examine these in detail (Figs. 2a, 2b). Figure 2a presents the latitudeítime diagram for the radial magnetic field Br calculated from the direct measurements of the large-scale longitudinal magnetic field. Figure 2b shows the latitudeítime diagram for the Br calculated from the measurements of the large-scale longitudinal magnetic field. Let us analyze these two plots jointly. The dark gray regions that shade into black near the poles correspond to N polarity of the magnetic field (> 0), while the light gray regions that shade into white near the poles correspond to S polarity (< 0). The motion of regions of large-scale magnetic field of each polarity from the equator toward the poles is clearly visible (compare with Fig. 1a). This drift of the large-scale fields differs fundamentally from the drift of the local fields, in particular from the motion of the Maunder butterflies from middle latitudes toward the equator. The drift rate of the large-scale fields depends on the latitude. From the equator to latitudes of 20í 25 , the drift is fairly rapid, and the magnetic field passes through this interval in two to three years. At latitudes of 25í50 , the drift rate sharply decreases (it does not exceed 1 m/s), and the magnetic field passes through this 25 interval in 15 yrs. Further,

the drift rate rapidly increases near the pole, and the magnetic field passes through the remaining 40í50 in about one year. Thus, the total time for the transport of the large-scale fields is approximately 17í18 yrs, i.e., there is an extended cycle. It is interesting that the centers of the regions of the large-scale fields of each polarity are always situated near maxima of the local field at the equator and minima of the local field at the poles. Each regular region of largescale field of a certain polarity arises at the equator when there is a reversal of the polarity of the two preceding fields at the poles (the last field, with the same polarity as the current equatorial field, and the next field, with the opposite polarity). The width of the region of each polarity in time is 11 yrs. During this time interval, three regions of large-scale field of alternating polarity corresponding to three different cycles lie one above the other at latitudes from the equator to the poles. Recall that the extended cycle of the local fields lasts 12í13 yrs, and there is a partial overlapping of two neighboring cycles at the latitudes of the wings of the Maunder butterflies. Note also that the maximum strength of the local fields is observed in a narrow latitude interval near 20 , whereas the maximum strength for large-scale fields is observed near the poles. Comparing these large-scale field diagrams with the Maunder butterflies, we arrive at the following conclusions. (1) There is a clear, mutually opposite motion of the large-scale and local fields: the large-scale fields move from the equator toward high latitudes, while
ASTRONOMY REPORTS Vol. 47 No. 4 2003


MERIDIONAL DRIFT OF LARGE-SCALE SOLAR MAGNETIC FIELDS

337

the local fields move from the high latitudes toward the equator. In each activity cycle, the local fields arise before the polarity reversal of the large-scale fields (by two to three years), and disappear near the equator simultaneous with the next polarity reversal of the large-scale fields, covering a time interval of 12í 13 yrs. The narrow intersection of these two regions occurs at latitude 20 . (2) The inclination of the butterflies to the time axis is virtually constant (a linear approximation is good) during their entire motion from the middle latitudes toward the equator. This inclination is about 5 per year. A butterfly moves through its whole path over 10í11 yrs, from the minimum of each cycle to the minimum of the following cycle. The central area of the butterfly (where there is the greatest number of spots) corresponds to the maximum of the local cycle, and coincides in latitude (20 ) with the intersection of the diagrams of the local and large-scale fields. The behavior of the large-scale fields differs from that of the local fields. First, the drift rate is not constant during the whole time interval, as was indicated above. Second, the large-scale fields move through their path from the equator to the middle latitudes, on average, over a longer time of about 17 yrs; i.e., from the maximum of the local cycle nearly to the minimum of the cycle following the next (the time is short by only one to two years). These conclusions agree with the results presented in Fig. 3, which shows the correlations between the large-scale magnetic fields at various latitudes. The cross-correlation of the equatorial radial field with the radial field at other latitudes is calculated. We can see that the correlation region extends from the equator to the poles, and its total length exceeds 200 Carrington revolutions, or 15 yrs. Similar findings are presented in [14]. 4. GENERAL SCHEME FOR THE MERIDIONAL CIRCULATION OF THE LARGE-SCALE MAGNETIC FIELDS To ascertain the drift characteristics obtained from the magnetograph data, we approximate the field distribution in the latitudeítime diagram (see the northern hemisphere in Fig. 2a) using the function Br / () =


Phase, radians 6

4

2

0

20

40

60

80 Latitude

Fig. 4. Phase of the meridional drift model as a function of latitude.

Drift rate, m/s 25 20 5 0 5 0 20 40 60 80 Latitude

Fig. 5. The model rate for the meridional drift of the largescale magnetic fields.

A ()exp(i (t + ())), A2 = 1.


(1)

(2)

We use this form of the Fourier transformation in order to analyze specifically the meridional drift, and to separate it from the latitudinal variations of the strength of large-scale fields. Therefore, we first
ASTRONOMY REPORTS Vol. 47 No. 4 2003

calculate the average Br and variance ( ) over the entire time interval (1960í2002) for the each latitude. The average Br turns out to be close to zero (as was expected), and all the Br (, t) are then normalized to the corresponding ( ). Further, the amplitude A ( ), frequency , and phase ( ) are found for each latitude via a Fourier-series expansion with fractional harmonics. As was mentioned above, the frequency corresponding to the maximum (about 20 yrs) is virtually independent of latitude and equal to 8.706 rad/s (i.e., the period is 20.87 yrs). We shall call this the "22-year cycle" and denote the related parameters by the subscript 22. The phase corresponding to this period is a rather complicated function of , accurate to within an arbitrary initial value. Figure 4 presents the function ( ), which can be


338

OBRIDKO, SHELTING

(), mu T A22 () 400

rate was obtained in [13, 27, 28]. At higher latitudes from 33 to 53 , the rate is extremely small, 1í3 m/s. Further, the rate rapidly increases to 25 m/s near the pole. The functions ( ) and A22 ( ) (Fig. 6) also indicate these two zones. The function ( ) is the r.m.s magnetic field calculated for a given latitude over the long time interval of two magnetic cycles. We can see that the r.m.s. strength of the large-scale field in the polar zone exceeds that of the equatorial maximum at latitudes 16 í18 by a factor of five or more. The function A22 ( ) describes the contribution of the 22year cycle at each latitude. Here, these two maxima are also seen, but the difference between them is much smaller. The value of A22 ( ) at high latitudes exceeds that at latitude 20 by only 25%. The quantitative results of the foregoing analysis are presented in the Table. Figure 7 shows a scheme for the cyclic variation of the meridional drift obtained by the formula Br = A22 ()exp(i22 (t + 22 ())). This means that we use here two additional types of filtering besides the three previous types presented in the Introduction. First, we add frequency filtering to eliminate every cycle except the 22-year cycle. Second, we remove the strong dependence of the largescale fields on latitude, expressed in Eq. (1) by the denominator ( ). 5. DISCUSSION The drift of the large-scale fields from the equatorial zone toward the pole takes 15í16 yrs (or 16í 17 yrs according to another estimate), i.e., three fourths of the 22-year cycle. The reverse motion, which is not directly observed, takes 5í6 yrs. The drift on the surface during these 5í6 yrs is virtually equal to zero. It then abruptly grows, and its rate exceeds 10 m/s. There is an abrupt retardation of the drift at latitudes 30í50 , and a stagnation region arises, where the rate does not exceed several meters per second. The drift again becomes rapid at higher latitudes. The stagnation region coincides with the region where the radial gradient of the rotational velocity in the convective zone is close to zero (there is a change of sign in / r ). At the same time, the stagnation region is extended in latitude, and its high-latitude and low-latitude boundaries appear as neutral curves in the latitudeítime diagrams. The variation of these boundaries and its secular cycle were studied in [29, 30]. The average rate of the polar drift is 2 m/s.
ASTRONOMY REPORTS Vol. 47 No. 4 2003

300

200

100 A

0

20

40

60

80

Fig. 6. The variance and amplitude of the 22-year cycle as functions of latitude in the model for the meridional drift.

easily recalculated to give the rate of the meridional drift, shown in Fig. 5. Figure 5 indicates that there are two regions of rapid drift. The drift rate grows from almost zero near the equator and reaches its maximum of 13 m/s at a latitude of 20 . Approximately the same
Table 1. Quantitative characteristics of the drift of the large-scale magnetic fields Latitude 90 75 64 56 50 44 39 35 30 26 22 17 13 10 6 2 0


A() 0.6226 0.6227 0.5932 0.4574 0.2950 0.2659 0.2240 0.1384 0.1331 0.3148 0.4302 0.4881 0.4976 0.4865 0.4236 0.2190 0.1108

A2 () 0.3877 0.3878 0.3518 0.2092 0.0871 0.0707 0.0502 0.0191 0.0177 0.0991 0.1851 0.2383 0.2476 0.2366 0.1795 0.0480 0.0123

() 374.8899 251.5153 113.0285 57.6674 47.2941 46.7928 44.6023 41.5526 40.0913 44.7899 53.7202 63.0929 63.6749 58.8488 47.0516 36.1304 34.9017

Phase 6.31828 6.2448 6.06852 5.71597 5.01434 4.62862 4.42442 4.14342 2.11691 1.72671 1.62887 1.56958 1.53136 1.49697 1.41616 0.02248 0.000


MERIDIONAL DRIFT OF LARGE-SCALE SOLAR MAGNETIC FIELDS

339

90 0.600 í 0.800 0.400 í 0.600 0.200 í 0.400 0 í 0.200 í 0.200íí0 í 0.400 í í 0.600 í 0.600 í í 0.400 í 0.800 í í 0.600

60

30

0 1960

1970

1980 Year

1990

2000

Fig. 7. Scheme for the cyclic variation of the meridional drift.

Let us consider the agreement between these findings and classical dynamo models. In most models (see, for example, [31] using the model of the convective zone [32]), the meridional drift naturally arises in the generation process due to the negative dynamo number at the base of the convective zone. This drift is directed toward the equator, and explains the Maunder butterflies for the local fields. These models neglect drift motions in higher layers of the convective zone [33]. The meridional drift of the largescale fields toward the pole is frequently explained as a consequence of the diffusion of tail polarities of the local fields due to differential rotation (see [10, 11] and references therein). However, in this case, it is difficult to understand why the equatorial fields follow the evolution of the polar fields with a time lag of 5í6 yrs (see [29] and references therein). There are also other direct objections to this scheme with the diffusion transport of some portion of the local fields toward the pole. It is unclear how the largescale, deep, poloidal magnetic field can originate from surface elements in the form of sunspots [34, 35]. Finally, it remains unclear why the large-scale fields in the polar zone significantly (by a factor of five or more) exceed the equatorial large-scale fields in the case of simple transport of the tails of active regions toward the pole. A flux of material directed toward the pole is postulated in [18, 19], where it is shown that
ASTRONOMY REPORTS Vol. 47 No. 4 2003

the calculated drift can agree with observations for appropriately chosen parameters. However, the time interval used for the comparison with observations is short, and the polar drift covers only 10 yrs. One qualitative scheme for the generation of large-scale fields in the solar cycle is the following. We will describe it in terms of the concentrated fields (CF) and diffusive fields (DF) introduced in [18], which we consider to be more convenient from the physical point of view than the terms toroidal and poloidal fields that are usually used in dynamo theory. We shall use the structure of motions in the convective zone described in [32]. Figures 8a and 8b show the important d/dr and d/d diagrams. The strongest radial gradient takes place within a narrow layer of thickness 0.05R0 [36] situated directly beneath the base of the convective zone at latitudes lower than 36 (sin < 0.6). This is the tachoclinic region where the - dynamo operates [37]. Here, the dynamo creates strong concentrated fields. Since < 0 in this region [31], this region of CF drifts from the middle latitudes toward the equator, while the fields rapidly (over several days) rise to the surface radially due to the buoyancy force. These fields are tubes of CF, with their meridional drift during their ascent being insignificant. For a meridional drift rate of 10 m/s, the tube drift toward the pole over


340

OBRIDKO, SHELTING

, radians 1.2 1.0 0.8 0.6 0.4

41í60 22í 41 25í 22 í 17í 25 í36íí17 í55íí36 í 74í í55 í93íí74 í 1.1E2í í93 í 1.3E2í í1.1E2 í 1.5E2í í 1.3E2 í 1.7E2í í1.5E2

0.2 0 0.6 1.2 1.0 0.8 0.6 0.4

0.7

0.8

0.9

1.0 r/r

700 í 800 600 í 700 500 í 600 400 í 500 300 í 400 200 í 300 100 í 200 0 í 100 í100í0 í 200í í 100 í 300í í 200 í 400í í 300 í 500í í 400 í 600í í 500 í 700í í 600 í 800í í 700

0.2 0

0.7

0.8

0.9

r/r

Fig. 8. d/dr (top) and d/d (bottom) diagrams.

15 days does not exceed 104 km, remaining less than the size of the tube itself. As a result, we observe on the surface Maunder butterflies directed toward the equator, reflecting the drift of the "spot creating" zone in the tachoclinic region. However, the spot tubes are not directly connected to the tachoclinic region, and are immersed below the solar surface no deeper than 6í12 000 km [34, 35]. The sunspots thus remain surface phenomena. The (or 2 ) dynamo operates simultaneously

[38, 39]. It is more efficient in higher layers, where the parameter changes sign for R 0.72R0 [31]. The evolution of diffuse fields is independent of the CFs [18], since the equation for weak, large-scale DFs branches off.The dynamo wave becomes the lower boundary condition for the dynamo. Thus, the DFs arise in a wider region of the convective zone, rather than in the narrow tachoclinic region. There are no sharp boundaries for these fields, and the magneticpressure gradients are weak. Therefore, these fields
ASTRONOMY REPORTS Vol. 47 No. 4 2003


MERIDIONAL DRIFT OF LARGE-SCALE SOLAR MAGNETIC FIELDS

341

move along with the medium. They rise slowly (they have low buoyancy) and are strongly dragged by the meridional drift in the convective zone. These fields can be additionally strengthened in the polar region due to the strong negative gradient d/d (Fig. 8b). The rates of the meridional drift obtained in [33] from helio-seismological data agree well with our results. Though the DFs arise above the CFs, they must be considered as being deeper when interpreting the observations. Therefore, we usually treat the largescale and global magnetic fields as deep phenomena (see the discussion in [14]). Stix [40] has calculated a similar pattern for the drift of the large-scale fields toward the pole using an dynamo model. It is extremely difficult to reveal the reverse drift in both our data and the helio-seismological data. The analysis of the cross-correlation of the largescale field strengths in [14] does not indicate any direct drift, but shows that the equatorial structure is similar to the polar structure with a time lag of 5í 6 yrs. This implies a rate for the reverse drift of 5í 6 m/s. Giles [33] likewise was not able to demonstrate the presence of a reverse drift directly using helioseismological data, but estimated from the continuity equation the reverse rate to be 3 m/s and the level in the convective zone where the meridional flow toward the pole is changed to the reverse flow to be 0.80R0 . In our discussion, we have used a model for the convective zone constructed especially for the rising branch of the 23rd cycle. If this model is changed, some of our results will require refinement. All our calculations have used a "classical" model in which the source surface is 2.5R0 from the Sun's center, and the field at the surface of the photosphere is potential. The calculations in an alternative model in which the photospheric field is taken to be radial and the source surface is at 3.25R0 give no significantly new results. The diagram obtained is virtually the same as that in Fig. 2a.

REFERENCES
1. K. Harvey, in The Solar Cycle , Ed. by K. L. Harvey, Astron. Soc. Pac. Conf. Ser. 27, 335 (1992). 2. V. Bumba and R. Howard, Astrophys. J. 141, 1502 (1965). 3. T. L. Duvall, Jr., Sol. Phys. 63, 3 (1979). 4. R. Howard and B. J. LaBonte, Sol. Phys. 74, 131 (1981). 5. B. J. LaBonte and R. F. Howard, Sol. Phys. 80, 361 (1982). 6. V. I. Makarov and K. R. Sivaraman, Sol. Phys. 119, 35, (1989). 7. V. I. Makarov, M. P. Fatianov and K. P. Sivaraman, Sol. Phys. 85, 215 (1983). 8. R. K. Ulrich, J. E. Boyden, L. Webster, et al., Sol. Phys. 117, 291 (1988). 9. Y.-M. Wang, J. Lean, and N. R. Sheeley, Geophys. Res. Lett. 27, 505 (2000). 10. Y.-M. Wang, A. G. Nash, and N. R. Sheeley, Astrophys. J. 347, 529 (1989). 11. Y.-M. Wang, A. G. Nash, and N. R. Sheeley, Science 245, 712 (1989). 12. V. N. Obridko and G. Gaziev, in The Solar Cycle ,Ed. by K. L. Harvey, Astron. Soc. Pac. Conf. Ser. 27, 410 (1992). 13. R. W. Komm, R. F. Howard, and J. W. Harvey, Sol. Phys. 147, 207(1993). 14. E. V. Ivanov and V. N. Obridko, Sol. Phys. 206, 1 (2002). 15. M. Perez Garde, M. Vazquez, H. Schwan, and H. Wohl, Astron. Astrophys. 93, 67 (1981). 16. B. N. Anderson, Sol. Phys. 94, 49. (1984). 17. G. Lustig and H. Wohl, Astron. Astrophys. 229, 224. (1990). 18. M. Dikpati and A. R. Choudhuri, Astron. Astrophys. 291, 975. (1994). 19. M. Dikpati and A. R. Choudhuri, Sol. Phys. 161, 9 (1995). 20. A. R. Choudhuri and M. Dikpati, Sol. Phys. 184, 61 (1999). 21. Y.-M. Wang, N. R. Sheeley, Jr., and J. Lean, Geophys. Res. Lett., 27, 621 (2000). 22. V. N. Obridko and B. D. Shelting , Sol. Phys. 201, 1 (2001). 23. B. D. Shelting and V.N.Obridko, Astron. Astrophys. Trans. 20(3), 491 (2001). 24. V. N. Obridko and B. D. Shelting , Sol. Phys. 184,187 (1999). 25. J. T. Hoeksema and P. H. Sherrer, Solar Magnetic Field "-- 1976 through 1985 (WCDA, Boulder, 1996). 26. K. G. Ivanov and A. P. Kharshiladze, Geomagn. Aeron. 34, 22 (1994). 27. R. K. Ulrich, in IAU Colloq. 137: Inside the Stars (A93-53126 23-90, 1993), p. 25. 28. H. B. Snodgrass, Sol. Phys. 94, 13 (1984). 29. V. I. Makarov, A. G. Tlatov, D. K. Callebaut, et al., Sol. Phys. ,198, 409 (2001). 30. V. I. Makarov, A. G. Tlatov and K. R. Sivaraman, Sol. Phys. 202, 11 (2001).

6. ACKNOWLEDGEMENTS The authors are grateful to E. V. Ivanov, K. M. Kuzanyan, and D. D. Sokolov for helpful discussions, to V. I. Makarov and A. G. Tlatov for providing us with unpublished data, and to the Wilcox Solar Observatory staff for the data received via INTERNET. This work was supported by the Russian Foundation for Basic Research (project codes 00-15-96661 and 0202-16199) and by INTAS (project 2000-840).
ASTRONOMY REPORTS Vol. 47 No. 4 2003


342

OBRIDKO, SHELTING 37. E. A. Spiegel and J. P. Zahn, Astron. Astrophys. 265, 106 (1992). 38. H. K. Moffat, in Magnetic Field Generation in Electrically Conducting Fluids (Cambridge Univ., Cambridge, 1978). 39. F. Krause and K.-H. Radler, in Mean Field Magnetogydrodynamics and Dynamo Theory (Pergamon, Oxford, 1980). 40. M. Stix, Astron. Astrophys. Trans. 20(3), 417 (2001).

31. G. Belvedere, K. M. Kuzanyan, and D. D. Sokoloff, Mon. Not. R. Astron. Soc. 315, 778 (2000). 32. J. Shou et al., Astrophys. J. 505, 390 (1998) . 33. P. M. Giles, Time-Distance Measurements of Large-scale Flows in the Solar Convection Zone, Dissertation, Stanford University (1999). 34. V. N. Obridko, Sun Spots and Activity Complexes [in Russian] (Nauka, Moscow, 1985). 35. J. Zhao, A. G. Kosovichev, and T. L. Duvall, Jr., Astrophys. J. 557, 384 (2001). 36. A. G. Kosovichev, Astrophys. J. 469, L61 (1996).

Translated by V. Badin

ASTRONOMY REPORTS Vol. 47

No. 4

2003