Yarema Zyelyk

RESEARCH TASKS
At the research of magnetic fields on the Sun on the basis of observations measure a longitudinal magnetic field, i.e. normal component of magnetic field strength. Thus in each point of the Sun surface in some current realization of the field we observed either a field "+" (N) - northern or "-" (S) - southern polarity. The field "+" (N) - northern polarity is directed from the Sun, and a field "-" (S) - southern polarity - is directed to the Sun
The imbalance of fields in each point of the Sun surface can be calculated as the vectorial sum of the N-and S-polarities field strengths. The study of the unbalance of "+" and "-" fields on the Sun is actual as the N-and S-polarities unbalance stream of magnetic fields from all Sun substantially determines a state of an interplanetary magnetic field and influences on the solar-terrestrial relationship
Differential character of the magnetic field rotation is exhibited in the exclusive dependence of the rotation period (velocity) from the latitude, not stacked in the conventional mechanistic representations about the solid-state rotation (on equator rotation is faster, than at high latitudes)
The purpose of the present study is the estimation of the latitudinal drift of the unbalance magnetic field stream as a whole and determination of the differential rotation and the latitudinal drift of the unbalance stream components possessing the certain significant rotation periods, and also comparison of these two drift types

Synoptic maps data of longitudinal vector of the magnetic field strength on the Sun on the 26 years observational interval since January, 1, 1977 till the September, 30, 2003, obtained in the NSO Kitt Peak (USA) and submitted to the author
The strength module intervals of the solar magnetic fields of the N and S polarities, allocated on the basis of the analysis of the histograms of the map pixels amount dependence from the strength
OBSERVATIONAL MATERIAL

Time series constructing
For each of the allocated 7 strength intervals in the range [0; 1500] Gs in everyone of 28 series each term of series represents the magnetic fields unbalance stream in the strength interval Hk (k = 1, 2,…,7) on the central meridian of the Sun from the surface element in a latitude band φn (n = 0, ±5, ±10, ±15, …, ±70) width's 1° longitudes and height's 5° latitudes in an instant time tl (l = 1, 2, …, 128520):
(|N|-|S|)(Hk , φn, tl),
where (|N|-|S|) – signify a difference of the absolute values of the streams of the magnetic fields «+» (N) и «-» (S) polarities.
The sampling time is ∆t=0,0757 day (1° of the solar longitude)

Temporal variations of the magnetic fields unbalance stream in the strength interval [0; 5] Gs for five 5° latitude bands at the latitudes from the interval of [5°; 25°] of the northern Sun hemisphere

LATITUDINAL DRIFT OF THE MAGNETIC FIELDS UNBALANCE STREAM AS A WHOLE
At the estimation of the latitudinal drift of the magnetic fields unbalance stream the method of the time delay definition of the unbalance stream as a whole as a wave train in a latitude band with a latitude |n-1| concerning the unbalance stream in the next band with the greater latitude |n|, based on the definition of a delay of cross correlation function which corresponds to a absolute maximum of this function was used:
Lag = arg{max Rn n-1(τ)}
Latitudes of the bands |n-1| and |n| are specially taken modulo that both in southern, and in northern hemispheres of the Sun the band |n| corresponded to higher next latitude, than a band |n-1|

The latitudinal variations of the cross-correlation function shifts in the next latitude bands for the magnetic fields unbalance stream on the Sun in the strength range of [0; 700] Gs
Lag > 0 – Drift to poles
Lag < 0 – Drift to equator

The latitudinal dependences of the latitudinal drift velocities of the magnetic fields unbalance stream in the strength range of [0; 700] Gs
V > 0 – Drift to poles
V < 0 – Drift to equator

CROSS-CORRELATION DRIFT ESTIMATION RESULTS
The analysis of the obtained latitudinal dependences of the cross-correlation function time shift and velocities of the latitudinal drift of the magnetic fields unbalance stream in the strength range of [0; 700] Gs results in following conclusion
For the fields of strength [0; 200) Gs the pattern of the latitudinal drift velocity variation is qualitatively similar for the all strength ranges which are included in this interval. From the latitudes of ± (15° - 20°) the unbalance stream drift in each hemisphere of the Sun to the areas of high and low latitudes takes place. With move off this latitude interval the drift velocity decreases up to 5-10 m/s at latitudes of ±60° and on equator
For the fields of strength [200-700) Gs the pattern is essentially another. Since at all latitudes Lag < 0 takes place and the latitude drift velocity is negative, the movement from high latitudes to low latitudes is observed only. The greatest velocities are observed at latitudes of ±(15° - 20°) and of ±(40° - 45°)
For the fields of all strength values there is a noticeable difference in velocities of the latitudinal drift between northern and southern hemispheres of the Sun

LATITUDINAL DRIFT OF THE SEPARATE COMPONENTS OF THE UNBALANCE STREAM, POSSESSING CERTAIN SIGNIFICANT ROTATION PERIODS
At the estimating of the latitudinal drift of the unbalance stream separate components, possessing certain significant rotation periods the spectral analysis methods are used.
Using Welch method the power spectral density Gn(P) in each latitude band |n| and the cross power spectral density (CPSD) Gn n-1(P) (P - rotation period) processes in the next latitude bands |n| and |n-1| are estimated
Regarding the significant peaks Gn(P) the significant rotation periods P_peak are determined, and on the bases of Gn n-1(P_peak) the phase angle of CPSD θn n-1(P_peak) is calculated
On the bases of θn n-1(P_peak) the temporal delay τn n-1(P_peak) for harmonic component of the unbalance stream rotating with significant period P_peak is determined
On the known distance along the meridian between medial lines of the next latitude bands and on the estimated temporal delay τn n-1(P_peak) the velocity of the latitudinal drift of these components to equator or to a pole is determined

Spectral methods of the rotation estimation implementation: Shuster Periodogram, Welch, Multitaper

Number of samples in each time series – N = 357x360 = 128520
Sampling time – dt = 27,2753/ 360 = 0,07576572 (day)
Sampling frequency – fs = 1/ dt = 13,19875491 (1/day)
Nyquist frequency – nyquist = 1/2 fs = 6,599377 (1/day)
Number of time samples subjected to FFT – nfft = 262144
Number of frequency samples
(Fourier coefficients) – M= nfft/2 + 1 = 131073
Sampling and spectral analysis parameters

In the same latitude band can exist from one up to four different rotation periods corresponding to significant PSD peaks.
These different rotation periods presented in the some current band, are spread to different on the width intervals of the next latitude bands everyone.
An example of the power spectral density dependences from the period in the period interval of [20; 35] day for the time series of the magnetic fields unbalance stream in the strength interval of [10; 20] Gs in the latitude band interval of [5°; 35°] in the northern Sun hemisphere is represent in the figure

Latitudinal intervals (in degree) of the determined periods of significant PSD peaks in the 5 ranges of magnetic field strength

Phase angle, time shift and velocity of the latitudinal drift determination
The phase angle θn n-1(P) of cross spectral density Gn n-1(P) at each value of period P is determined as follows:
θn n-1(P)=arctg {-Im[Gn n-1(P)]/ Re[Gn n-1(P)]} * 180/ π (degree)
The positive phase angle θn n-1(P) > 0 in a latitude band n show, on how many degrees on a longitude the harmonic component with period P of the rotation process of the magnetic fields unbalance stream in a latitude band (n-1) lags behind from harmonic component with period P in a band n. At θn n-1(P) < 0 the corresponding harmonic component with period P advancing takes place
With the phase angle θn n-1(P) is unequivocally related the time shift τn n-1(P) at the corresponding rotation period P which is determined per phase angle as follows:
τn n-1 (P) = θn n-1 (P) * P / (2 π) (day)
After determination time shift τn n-1 (P) for any value of significant rotation period P it is possible to estimate the velocity Vn n-1 (P) of the latitudinal drift of harmonic component with period P of rotation process of the magnetic fields unbalance stream:
Vn n-1 (P) = dn n-1/ [τn n-1 (P)* 86400] (m/s),
where dn n-1= 2 π * RSun* 5/ 360 (m) - distance on the meridian between medial lines of the next 5° latitude bands (meter), and RSun= 696000000 m is radius of the Sun
The sign of the latitudinal drift velocity Vn n-1 (P) coincides with the sign of the time shift τn n-1(P) and, accordingly, phase angle θn n-1(P)

The significant peaks PSD, phase angle, time shift and latitudinal drift velocity estimation procedure
Magnetic fields unbalance stream PSD. Significant period P = 28,1 day in the latitudinal interval [-40 °; -20 °]. The magnetic field strength range is [20; 50] Gs

The phase angle dependence on period and the monotonic character of its variation in neighborhood of P=28,1 day for the magnetic field strength range of [20; 50] Gs

The time shift on period dependence and the monotonic character of its variation in neighborhood of P=28,1 day for the magnetic field strength range of [20; 50] Gs

The latitudinal drift velocity on period dependence in the latitudinal interval of [-40°;-20°] in which there is a significant rotation period of P=28,1 day
Magnetic field strength range is [20; 50] Gs
The phase angle on period dependence in the latitudinal interval of [-40°;-20°] in which there is a significant rotation period of P=28,1 day
Magnetic field strength range is [20; 50] Gs

Latitudinal dependences of the phase angle of the magnetic fields unbalance stream components rotating with periods in the interval of [27; 28) day for five magnetic field strength ranges

Latitudinal dependences of the latitudinal drift velocity of the magnetic fields unbalance stream components rotating with periods in the interval of [27; 28) day for five magnetic field strength ranges

The dependences of significant rotation periods of the unbalance stream structures from the emersion latitude for different magnetic field strength ranges
DIFFERENTIAL ROTATION AT THE LEVEL OF THE BACKGROUND FIELD SOURCE

The correlation and spectral models of the stationary random processes describing the fields rotation in the assigned latitude bands of the Sun have been effectively applied for the estimation of stationary differential rotation and the latitudinal drift of magnetic fields on the Sun
The comparison of the estimates of the cross correlation function time shifts with phase angles of the cross spectral density in the next latitude bands and the latitudinal drift velocities of the unbalance stream as a whole with the latitudinal drift velocities of the unbalance stream structures possessing the significant rotation periods has shown qualitative coincidence of the latitudinal drift of these two types
Thus, at the magnetic field strength of [0; 200) Gs the latitudinal drift of the unbalance stream as a whole is determined generally a latitude drift of its structures rotated with significant periods in some latitude intervals, enveloping from three and more five-degree latitude bands
CONCLUSION