Open AccessArticle

Dynamics of Phase Synchronization between Solar Polar Magnetic Fields Assessed with Van Der Pol and Kuramoto Models

Anton Savostianov

^1,2,*

Alexander Shapoval

and

Mikhail Shnirman

Gran Sasso Science Institute, viale F. Crispi 7, 67100 L’Aquila, Italy

National Research University Higher School of Economics, Myasnitsakaya Ulitsa 20, 101000 Moscow, Russia

Institute of Earthquake Prediction Theory and Mathematical Geophysics, Russian Academy of Sciences, Profsoyuznaya Ulitsa 84/32, 117997 Moscow, Russia

Author to whom correspondence should be addressed.

Entropy 2020, 22(9), 945; https://doi.org/10.3390/e22090945

Submission received: 29 July 2020 / Revised: 23 August 2020 / Accepted: 25 August 2020 / Published: 27 August 2020

(This article belongs to the Collection Feature Papers in Information Theory)

Download

Browse Figures

Review Reports Versions Notes

Abstract

We establish the similarity in two model-based reconstructions of the coupling between the polar magnetic fields of the Sun represented by the solar faculae time series. The reconstructions are inferred from the pair of the coupled oscillators modelled with the Van der Pol and Kuramoto equations. They are associated with the substantial simplification of solar dynamo models and, respectively, a simple ad hoc model reproducing the phenomenon of synchronization. While the polar fields are synchronized, both of the reconstruction procedures restore couplings, which attain moderate values and follow each other rather accurately as the functions of time. We also estimate the evolution of the phase difference between the polar fields and claim that they tend to move apart more quickly than approach each other.

Keywords:

solar hemispheres; solar faculae; synchronisation; coupled oscillators

1. Introduction

Scientists are challenged by existed regular patterns and examples of irregularities in the dynamics of the hemispheric asymmetry in the solar magnetic field [1,2]. The asymmetry in the amplitude is investigated numerically (see, f. e., [3]) and implemented into the flux transport dynamo models [2,4]. The hemispheres contribute non-equally into the solar activity; their dominance is changed inside various solar cycles including 12–15, 17–18, and 21–23 [5,6] The evolution of the phase difference exhibits its own patterns. The hemispheres are typically synchronized [7,8,9], and the phase-leadership of one of them is persistent within several cycles with only two changes in the leadership occurred near 1928 and 1968 during the previous century [10,11].

The northern and southern components of the solar magnetic field can be connected through meridional circulation which transports the field from the equator towards the poles and backward forming the cells [12]. Regular observations of the meridional circulation have been reported only recently. Their empirical investigation and explanation with the flux transport models reveal general mechanisms but still leave open questions regarding the number and geometry of the cells, the influence of the meridional circulation to the solar cycle, and the dynamics of the velocity profile of the meridional circulation [13,14,15,16,17,18].

We are motivated by Blanter et al. [11] who investigated the asymmetry between the hemispheres with the Kuramoto model of the coupling oscillators by reconstructing the dynamics of the phase difference between them. This approach complements other applications of the inverse problem widely exploited to estimate properties of the deep solar dynamo processes that cannot be observed or measured at the solar surface [19,20]. Despite the fact that the Kuramoto model is not derived from the first principles, the range of its applicability includes laser arrays, neural networks, chemical oscillators, and neural activity [21,22,23,24]. Ad hoc modeling, simplifying the underlying physical processes and capturing the effect of synchronization, can give an accurate solution of the inverse problem in just mentioned fields. The choice of an appropriate model is always a challenge. Savostianov et al. [25] linked the efficient reconstruction of the “hidden” features of the solar activity discovered with the Kuramoto model to analogous reconstructions performed with the Van der Pol model of the coupled oscillators, where the latter is considered as a (very) simple descendant of the MHD equations (after [26,27]).

The purpose of this paper is to introduce a measure of similarity in the reconstructions of the coupling performed with the Kuramoto and Van der Pol system of two oscillators. The developed theory is applied to the solar faculae time series, which represent the northern and southern polar fields, to investigate the phase asymmetry in the coupling between the solar hemispheres. Finally, we plan to restore the phase difference between the time series that results in the found level of the similarity in the coupling.

This paper develops our previous study which reveals the similarities in the reconstructions of the coupling between the International Sunspot Numbers ISSN and geomagnetic index aa inferred from the Kuramoto and Van der Pol models of the coupling oscillators [25]. These reconstructions are built on the assumption that the coefficients of the corresponding differential equations vary slowly; in other words, the properties of the solutions are inherited from the equations with the constant coefficients. Now we address the question regarding the dynamics of the coupling and the phase difference when estimating the proximity of the reconstructions.

The applicability of the theory of coupled oscillators to the polar faculae data is worth discussing. A strong correlation between the two time series is observed with a “naked eye”. This correlation is a consequence of the mutual coupling which exists as the chain of the following interactions. A poloidal magnetic field associated with a pole and represented by the solar faculae transits to the near-equator toroidal field. The transition occurs within one of the hemispheres. Further, the poloidal field of this hemisphere is connected to that of the other hemisphere through cross-equatorial links established because of the different polarities of the magnetic fields. Finally, the poloidal field of the other hemisphere is transformed into the toroidal field. The reverse coupling occurs in a symmetric way. Probably, this complex mechanism “trying” to synchronize the polar faculae represents the coupling of the solar hemispheres. As to the transition between the poloidal and toroidal fields, it is performed with the meridional circulation analyzed in empirical way, with ad-hoc models, and in the framework of the solar dynamo modeling [14,28,29,30].

The rest of the paper is organized in the following way. Section 2 introduces the data. Section 3 describes the method of the analysis introducing the two reconstructions. The main findings are placed in Section 4. Section 5 concludes.

2. Data

We use yearly MWO calibrated polar faculae data downloaded from the solar dynamo dataverse (https://dataverse.harvard.edu/dataverse/solardynamo), maintained by A. Munoz-Jaramillo [31]. The faculae data representing the corresponding solar polar magnetic fields share the common solar quasi-periodicity of ∼11 years with the solar cycle.

Despite the distant locations of the data sources, the faculae time series illustrate a strong phase synchronisation (Figure 1). Nevertheless, from time to time the series are de-synchronized. The most evident episode of de-synchronization occurred during cycle 19 in the 1960s. We investigate the synchronization of the series, exploring their dynamics in sliding windows and applying general models which efficiently describe synchronization.

3. Method: Reconstruction of the Coupling with Two Models

Here we describe the Kuramoto and Van der Pol models of two coupled oscillators. The inverse problem intended to reconstruct the coupling between the oscillators is formulated and solved, given the other parameters of the equations and their solution.

3.1. Kuramoto Model

In the case of the Kuramoto model, the oscillators are given by the system of the equations

\begin{matrix} x (t) = sin θ_{x} (t) \end{matrix}

(1)

\begin{matrix} y (t) = sin θ_{y} (t), \end{matrix}

(2)

where the dynamics of the phases is determined by the coupling coefficient

μ

and their natural frequencies

ω_{x}

and

ω_{y}

through the differential equations

\begin{matrix} {\dot{θ}}_{x} = ω_{x} - μ sin (θ_{x} - θ_{y}) \end{matrix}

(3)

\begin{matrix} {\dot{θ}}_{y} = ω_{y} - μ sin (θ_{y} - θ_{x}) . \end{matrix}

(4)

Simplifying comparison between different models, the time normalization is used in such a way that 1 is assigned to the synchronized natural frequency:

(ω_{x} + ω_{y}) / 2 = 1

. We introduce half a difference

Δ ω

between the natural frequencies:

Δ ω = \frac{ω_{x} - ω_{y}}{2} .

Thus, the direct problem is the description of the solutions

θ_{x}

θ_{y}

of Equations (3) and (4), given

Δ ω

and

μ

We also consider the inverse problem for the Kuramoto model formulated in the following way: given series

x (t)

y (t)

, which satisfy Equations (1)–(4), and the model parameter

Δ ω

, reconstruct the coupling coefficient

μ

between the oscillators. Following Blanter et al. [32], the solution of the inverse problem is inferred from the computation of the correlation between the free oscillators. We define the sliding correlation in the window of the length T as:

Corr [X, Y] (t; T) = \frac{E_{t, T} ((X - E_{t, T} X) (Y - E_{t, T} Y))}{\sqrt{E_{t, T} {(X - E_{t, T} X)}^{2} E_{t, T} {(Y - E_{t, T} Y)}^{2}}},

(5)

where

E_{t, T} X = T^{- 1} \int_{t - T / 2}^{t + T / 2} X (τ) d τ

, and put it into the correspondence to the center of this window. One can expect the relationship between the coupling and the correlation. The larger coupling is the stronger it causes the oscillators to follow each other and exhibit bigger values

ρ (t) = Corr [X, Y] (t; T)

of the correlation. When the oscillators are synchronized, so that

\dot{θ_{x}} - \dot{θ_{y}} = 0

, Equations (3) and (4) yield

μ sin (θ_{x} - θ_{y}) = Δ ω .

(6)

The condition

μ > Δ ω

provides the synchronization. Let

T = 2 π

be fixed to the period of the synchronized oscillators. Then the correlation between

x (t)

and

y (t)

computed with (5), where

x (t)

y (t)

are substituted for

X (t)

Y (t)

is given (see [32] for the proof, which consists of the straightforward computation of the integrals) by the equation

ρ = cos (θ_{x} - θ_{y})

(7)

with the autonomous right hand side. Equations (6) and (7) lead to the following reconstruction rule:

μ = \frac{Δ ω}{\sqrt{1 - ρ^{2}}}

(8)

Since the right hand side of (8) equals to

sin arccos (θ_{x} - θ_{y})

up to the sign, this sign describes the leadership among the two oscillators. Assigning only “+” as the sign, one eliminates the information about the leadership. Equation (8) also implies the following simple relation between the correlation and the coupling: the higher the coupling between two oscillators, the more correlated they are (thus, the smaller stationary phase difference they exhibit).

We are going to apply the reconstruction rule (8) when

x (t)

and

y (t)

represent observed oscillators. We stress that solar proxies are not fully synchronized, in contrast to the solutions of the model equations. Nevertheless, Equation (8) is formally applicable, as soon as

Δ ω

is given (the choice of

Δ ω

will be discussed later). Since the correlation between solar proxies varies with time (

ρ = ρ (t)

) it follows that the coupling reconstructed with the rule (8) is also time-dependent (

μ = μ (t)

). The usage of Equation (8) is plausible if the genuine coupling (i.e., that is between solar proxies) varies slowly with time:

\frac{d}{d t} (θ_{x} - θ_{y}) \approx 0 .

(9)

Equation (9) exhibit quasi-synchronization. The assumption of quasi-synchronization underlies our approach. Clearly, if the coupling does vary slowly, then the quasi-synchronization regime is quickly attained. The waiting time depends only on the distance from the initial condition to the attractor.

3.2. Van Der Pol Model

The pair of the Van der Pol (VdP) coupled oscillators is governed by the equations

\begin{matrix} \ddot{x} - (1 - x^{2}) \dot{x} + (1 - Δ ω) x + μ (\dot{x} - \dot{y}) = 0 \end{matrix}

(10)

\begin{matrix} \ddot{y} - (1 - y^{2}) \dot{y} + (1 + Δ ω) y + μ (\dot{y} - \dot{x}) = 0, \end{matrix}

(11)

where

μ

is the symmetrical coupling of the oscillators and the natural frequencies

1 \pm Δ ω

are introduced in such a way that the common synchronized frequency, i. e., their half a sum, is normalized to 1. The normalization is in line with that for the Kuramoto oscillators. Kuznetsov et al. [33] and Astakhov et al. [34] found similar characteristics of the direct problem for systems (3), (4) and (10), (11), whereas Savostianov et al. [25] discussed similarities of the inverse problems posed for these systems.

The formulation of the inverse problem coincides with that for the Kuramoto model. However, the solution is obtained numerically. Namely, Equations (10) and (11) are solved for a broad domain of the pairs

(Δ ω, μ)

and the correlation

ρ = ρ_{VDP} (Δ ω, μ)

of the solutions is computed with the window of the length T, where T is the period of the limit cycle of the synchronized oscillators. In numerical computations, we wait for a sufficiently long time prior to the data processing to avoid the influence of the initial conditions on the correlation. Figure 2 illustrates the relationship

ρ_{VDP}

between

ρ

Δ ω

, and

μ

. The reconstruction of the coupling

μ

is given by the inverse function

ρ_{VDP}^{- 1} (Δ ω, ρ)

(12)

also found numerically. The inverse function is well defined if the correlation

ρ

is positive and

| Δ ω | < μ

. Performing numerical computation, we assign 0 to

ρ^{- 1} (Δ ω, ρ)

when the correlation is negative.

Function

ρ_{VDP}^{- 1} (Δ ω, ρ)

shown in Figure 2 and the relationship described by Equation (8) have much in common. According to them, when the coupling strengthens or the phase difference shrinks, the oscillators evolve more alike each other and their correlation becomes bigger. A fall in the coupling strength causes a drop in the correlation (Figure 2, upward changes along vertical lines with fixed values of

Δ ω

3.3. Reconstruction Scheme

The Kuramoto and VdP models give the two ways of the reconstruction of the coupling exhibited by (8) and (12) respectively. The dynamics of the reconstruction is inferred from the time series

x (t)

and

y (t)

considered as the input when restoring coupling under the assumption of the quasi-stationarity of the solutions. We add into the consideration the dynamics of the differential equations by generating a new series. The full reconstruction procedure is displayed in Figure 3. It consists of the following steps.

(A): Given time series $X (t)$ , $Y (t)$ , and the model parameter $Δ ω$ , we reconstruct the coupling with both models. These series exhibit the solar cycle; $T^{*} = 11$ years is used in the paper to assign a single number to the variable length of the cycle. The time axis is initially stretched by $T^{*} / 2 π$ to transform the estimate of the cycle length into $2 π$ and set the correspondence between the time axis in the models and observations. Clearly, the linear transform does not affect either the correlation between the series $X (t)$ and $Y (t)$ or the variability of the solar cycle. The Kuramoto reconstruction is performed with (8) and denoted $μ_{k} (t)$ . The VdP reconstruction is performed with (12) and denoted $μ_{v} (t)$ . These two procedures are schematically displayed in the left two blocks of Figure 3. The both reconstructions $μ_{k} (t)$ and $μ_{v} (t)$ , in general, depend on time, since the input series represent the observations instead of the solutions of the model equations. The mathematical expectation of the input series is switched into the mean when the correlation is computed.
(B): Following Equation (7), we put

$θ (t) = arcsin \frac{Δ ω}{μ (t)}, x = sin (t + θ (t) / 2), y = sin (t - θ (t) / 2)$

(13)

and define the correlation of the series $ρ (t) = Corr [x, y] (t; T)$ over the period, where $ρ = ρ_{k}$ , $θ = θ_{k}$ , $μ = μ_{k}$ , $x = x_{k}$ , and $y = y_{k}$ when the Kuramoto model is investigated. The subscript v is used for the VdP model. The series $x (t)$ and $y (t)$ introduced by (13) represent free oscillators (the second block from the right in Figure 3).
(C): Finally, we repeat the reconstruction of the coupling from the time series and the phase difference (the first block from the right in Figure 3). Equation (8) is applied for the both types of the input to get ${\hat{μ}}_{k}$ and ${\hat{μ}}_{v}$ from $ρ_{k}$ and $ρ_{v}$ respectively; $Δ ω$ has been fixed during the steps (A)–(C). This part involves the dynamics of the equations into the reconstruction. Namely, the addressed question is how the dynamics of the coupling in the direct problem affects reconstruction. We end up with the reverse transform of the time axis and restore years as the units of the reconstructions found in the paper and displayed on the Figures.

The described procedure greatly benefits from the normalization applied in the modeling. Instead of tracking initial physical variables from the MHD equations (which also implies the existence of uncertainties when choosing the values of the parameters), we turn to the normalized quantities, thus shrinking the parameter space down to the coupling

μ

and the normalized frequency difference

Δ ω

. Baring that in mind, the coupling strength

μ

should be treated as a composite of various physical quantities related to the interactions of the oscillators.

3.4. Comparison of the Reconstructions

In this section, we discuss the proximity between the reconstructions

{\hat{μ}}_{k}

and

{\hat{μ}}_{v}

to each other. By construction, both reconstructions are located above

Δ ω

. Let

r_{Δ ω} [{\hat{μ}}_{k}, {\hat{μ}}_{v}] (t) = \frac{{\hat{μ}}_{k} (t) - Δ ω}{{\hat{μ}}_{v} (t) - Δ ω} .

(14)

This ratio is expected to exhibit larger variations than

{\hat{μ}}_{k} (t) / {\hat{μ}}_{v} (t)

, as both reconstructions can drop to a neighborhood of

Δ ω

. The stability of

r_{Δ ω}

as a function of time and, moreover, the proximity to 1 give evidence that the reconstructions with the Kuramoto and VdP models agree with each other.

The ratio

r_{Δ ω}

is well defined at the points t such that the correlation

ρ (t)

obtained at step (A) is positive and t belongs to the time interval offset by the period T from both ends of the considered time span. The two-time usage of the sliding windows of the length T generates the requirement about offset. The inverse function

ρ_{VDP}^{- 1}

computed through Figure 2 is well posed if the correlation is positive.

We note that there are other quantities besides the pair

({\hat{μ}}_{k} (t), {\hat{μ}}_{v} (t))

that are worth investigating when discussing the reconstruction of the coupling. The proximity of the initial reconstructions

μ_{k} (t), μ_{v} (t)

obtained at the step (A) also evidences in favor of the interchangeability of both models. The transition from

{\hat{μ}}_{*} (t)

μ_{*} (t)

, where the star stands either for k or v, leads to more smooth curves offset to the right (see [35]). Therefore, we’ve decided to restrict themselves here only to the comparison of the quantities

({\hat{μ}}_{k} (t), {\hat{μ}}_{v} (t))

obtained at the end of the reconstruction procedure.

4. Results

4.1. Reconstructed Couplings and Relation between Two Models

Interpreting

Δ ω

, one can say that the “period” of the quasi-periodic series

X (t)

and

Y (t)

belongs to the interval

[T^{*} - Δ ω T^{*}, T^{*} + Δ ω T^{*}]

, where

T^{*}

is a point estimation of the solar cycle period. The setting

Δ ω = 0.2

defines a feasible range of the values attained by the period.

The reconstructions

{\hat{μ}}_{k} (t)

and

{\hat{μ}}_{v} (t)

of the coupling with

Δ ω = 0.2

are displayed on the top panel of Figure 4. The pink vertical stripes indicate the time span where the reconstruction is not well defined. The computation of the correlation is performed with a smaller window when t belongs to the left and right pink time stripes. The middle pink time stripe corresponds to the values which are obtained a negative correlation between the series. In this case,

μ_{v}

is set to

Δ ω

, which is the least feasible value of the coupling. The two curves

{\hat{μ}}_{k} (t)

and

{\hat{μ}}_{v} (t)

follow each other rather accurately (Figure 4, upper panel).

The consistence between these curves conserves when the ratio

r_{Δ ω} (t)

defined by (14) is considered, Figure 4, lower panel. Outside 1958–1975 interval, the ratio

r_{Δ ω} (t)

varies slowly attaining values inside the range

0.83

–

0.91

Δ ω = 0.2

(black curve). Changes of

Δ ω

within the interval

[0.1, 0.3]

slightly affect the average value of

r_{Δ ω} (t)

but not the weak variability of

r_{Δ ω} (t)

, Figure 4 (lower panel).

4.2. Reconstruction of the Frequencies

Now we are going to find a variable frequency

Δ ω (t)

which generates a fixed value of

r_{Δ ω} (t)

which is attained most frequently. To this end, the range of

r_{Δ ω} (t)

, reachable at each time moment t with all possible values of

Δ ω

, is found (numerically). This range is located between the two curves shown in Figure 5. The values within the light red corridor are attained outside the interval 1958–1982. The width of the interval with attainable values of

r_{Δ ω} (t)

varies with time remaining almost unchangeable after 1980.

The above analysis supports the conjecture that the two reconstructions of the coupling are similar in terms of the ratio

r_{Δ ω} (t)

. More precisely, this

r_{Δ ω} (t)

demonstrates a typical range of the values, given fixed

Δ ω

from a broad interval

(0.025, 0.5)

, assuming a uniform distribution of

Δ ω

and focusing on the relevant values.

Then one can assess the phase difference between the oscillators assuming that

r_{Δ ω} (t)

is constant, but

Δ ω (t)

varies. We choose

0.86

as a “typical” value of

r_{Δ ω} (t)

and restore

Δ ω (t)

which results in this

r_{Δ ω} (t)

. The evolution of the restored

1 + Δ ω (t)

interpreted as the normalized phase frequency associated with the northern component of the polar magnetic field is shown in Figure 6. Small values of

Δ ω

in the late 1950s and around 1916 are observed (Figure 6) approximately 7–10 years prior to the changes in the leadership of the oscillators representing the northern and southern solar hemispheres (see Blanter et al. [11] and Deng et al. [36] for the discussion regarding the changes in the leadership and the time of their occurrence). The leadership of the northern hemisphere were ceased to exist in the 1920s and recreated in the 1960s. The second change is known to be much stronger. Since each restored value

Δ ω (t)

is obtained with the 22-years data centered at this t, these small values of

Δ ω

are constructed with the data which represent the changes in leadership. Our rough way to restore the natural frequencies failed at the time of this strong change in the leadership and resumed its functionality in the mid 1980s when the phase difference returned to usual values (Figure 6). A downward trend at the right part of the graph in Figure 6 may also attain values of the 1920s and 1960s, but with the data at hand only the level of the 1940s is reached.

The reconstructed dynamics of the natural frequencies is asymmetrical: the intervals characterized by the outward drift of the two hemispheres (around 1915, 1945, 2000) are noticeably shorter than the intervals with the phases moving to each other (e.g., 1920–1940, 1945–1960, 1990–2000), Figure 6. In general, this asymmetry is not novel and difficult to model. Following Syukuya and Kusano [37], one can attribute the asymmetry to the interactions between dipole- and quadrupole-type solution of the MHD model focused on the properties of the two hemispheres. With this model, Syukuya and Kusano [37] generate the slow phase convergence and quick divergence between the hemispheres (observed in Figure 6) when the corresponding solutions approaches the attractor.

Moreover, the MHD dynamo models admit various patterns in the long-term dynamics of the asymmetry, not limited by phase differences. Passos et al. [38] constructed the range in the feasible dynamics of the N–S asymmetry, including the grand-minima of only one of the hemispheres and almost complete synchronization as possible extreme cases. A few examples of the asymmetry in the phase-difference dynamics derived implicitly in our paper need further empirical justifications.

5. Conclusions

The paper compares the two reconstructions of the coupling between the polar faculae which represent the polar components of the solar magnetic field. The reconstructions are inferred from the Kuramoto and VdP models of the pair of coupled oscillators. Kuznetsov et al. [33] derived the Kuramoto-like equation for phases from the VdP equations Savostianov et al. [25] the established similarity of the inverse problem for these models finding excellent proximity between the two reconstructions of the coupling if the correlation between the oscillators is moderate. Better proximity between the input series makes the reconstructed values more remote. In this paper, we reconstruct the coupling between the solar hemispheres with the two models using the solar faculae as the input series. The reconstructed coupling is weak, in line with results of Norton and Gallagher [39].

The reconstruction with the Kuramoto model is obtained as the rigorous solution of the inverse problem proposed and efficiently applied to solar proxies by Blanter et al. [32]. The inverse problem for the VdP coupled oscillators is solved numerically. As a result, we obtain the coupling between solar hemispheres from the two models. Whilst the overall dynamics of reconstructed couplings qualitatively coincides, their difference exhibits significant variations, Figure 4. The proximity between the two reconstructions is measured with the function

r_{Δ ω} (t)

, defined by (14), which would be equal to 1 if the reconstructions were identical, Figure 4. This function

r_{Δ ω} (t)

varies slowly in time and exhibits a systematic bias towards the values that are a bit lesser than 1. Based on [25], we argue that the level of the correlation between the solar faculae series is larger than that required for the most similar reconstructions.

In this paper, VdP oscillations are considered as a (huge) simplification of the MHD-based models [26,27]. The proximity of both reconstructions links the efficient reconstruction of the coupling and phase differences between solar proxies performed by Blanter et al. [11,32] with the ad hoc Kuramoto models to the solar dynamo models inferred from the MHD-equations [2].

In general, arising long-term asymmetry in the solar hemispheres has been vastly investigated empirically, [5,8,40,41]. Various methods of the time-series analysis, including the Fourier transform, wavelets, singular spectrum analysis, and cross-recurrence plot are used to derive the periodicities and trends of the absolute

(N - S)

or normalized

(N - S) / (N + S)

asymmetries, where N and S represent the magnetic fields associated with the northern and southern hemispheres respectively. Nevertheless, only a few authors focused on the phase difference; the paper by Blanter et al. [11] gives a possible example. Complementing this paper, we end up with the conjecture that the solar hemispheres more rapidly diverge than converge, Figure 6.

We estimate the dynamics of the phase difference between the northern and southern hemispheres of the Sun (Figure 6). It is possible that the oscillators move apart more quickly than tend to each other. A general decrease observed from 1919 to 1958 including ∼2-years abrupt growth in the restored phase difference may be related to the traces of 30-to-40-years periodicities found with harmonic analysis in the amplitude characteristics of solar activity (e.g. [42,43,44]). Blanter et al. [45] found the same 3 cycle quasi-periodicities when modeling coupled oscillations between solar proxies with the Kuramoto equations used also in this paper. Two downward trends of the reconstructed normalized frequencies to moderate minima, which are attained in the 1940s and 1990s, were altered by a subsequent growth (Figure 6). In contrast, a deep minimum assigned to the 1950s occurred prior to the episode of the de-synchronization between the two hemispheres. The last at the right restored values of the normalized frequencies continue a downward trend but still remain larger than that of the 1950s. Only additional values can give evidence of how far the series is from a new episode of the de-synchronization.

We stress that the two principally different models usually agree with each other. This agreement is in favor of the reliability of the models. Roughly speaking, the epoch of the synchronized oscillators can be described not only with the MHD but also with simple ad-hoc models. However, both our models lose their adequacy when tackling cycle 20 characterized by the de-synchronization of the polar faculae. The failure of the models is predictable, as they are built on the assumption of the synchronization between the considered time series. However, the divergence between the models dealing with the de-synchronized data rises important questions regarding the nature of the physical mechanisms which stay beyond the model divergence. As such divergence stems from the different origins of nonlinearity (we work with nonlinearly coupled linear model and linearly coupled nonlinear model), further investigation of such simple models could shed light on the adequate choice of MHD parameters. As a result, a better understanding of the relationship between non-linear characteristics of our models and instabilities in the solar dynamo evolution opens the door to the prediction of these instabilities in advance.

The polar faculae were chosen in this paper as available long daily time series representing the polar magnetic fields. The computational procedure can be repeated with shorter time series exhibiting the components of the solar magnetic fields or sunspot data in order to compare the results.

In the same time, the numerical procedure establishing the relationship between the coupling, the phase difference, and the correlation of the VdP oscillators is useful. It functions efficiently when the variables exhaust a synchronization domain between the boundaries

Δ ω = 0

and

Δ ω = μ

corresponding to the absolute and disappearing synchronization respectively. The numerical analysis is worth being extended to the case of three coupled VdP oscillators and applied to the reconstruction of the velocity profile of the meridional flow. This approach will connect the results by Blanter et al. [28] inferred from the properties of the Kuramoto oscillators to that by Hazra et al. [29], Featherstone and Miesh [46], Choudhuri [47], Cameron et al. [48] obtained with elaborated versions of the solar dynamo model.

Author Contributions

Conceptualization, A.S. (Anton Savostianov), A.S. (Alexander Shapoval) and M.S.; methodology, A.S. (Anton Savostianov), A.S. (Alexander Shapoval) and M.S.; software, A.S. (Anton Savostianov); validation, A.S. (Anton Savostianov), A.S. (Alexander Shapoval) and M.S.; investigation, A.S. (Anton Savostianov), A.S. (Alexander Shapoval) and M.S.; writing—original draft preparation, A.S. (Anton Savostianov); writing—review and editing, A.S. (Alexander Shapoval) and M.S.; visualization, A.S. (Anton Savostianov); supervision, A.S. (Alexander Shapoval) and M.S. All authors have read and agreed to the published version of the manuscript.

Funding

A.Sh and A.S (at the initial stage of this research) were benefited from the support of the Laboratory of of Complex Systems Modeling and Control at the National Research University Higher School of Economics.

Acknowledgments

An earlier version of this paper has been presented during Space Climate 7 meeting organized in Orford, Canada, in July 2019. The authors are thankful to the participants of the meeting who commented the presentation and helped improve the research.

Conflicts of Interest

The authors declare no conflict of interest.

References

Charbonneau, P. Solar Dynamo Theory. Annu. Rev. Astron. Astrophys. 2014, 52, 251–290. [Google Scholar] [CrossRef] [Green Version]
Hathaway, D.H. The Solar Cycle. Living Rev. Sol. Phys. 2015, 12, 4. [Google Scholar] [CrossRef]
Nagovitsyn, Y.A.; Kuleshova, A.I. North–South asymmetry of solar activity on a long timescale. Geomagn. Aeron. 2015, 55, 887–891. [Google Scholar] [CrossRef]
Shetye, J.; Tripathi, D.; Dikpati, M. Observations and modeling of north-south asymmetries using a flux transport dynamo. Astrophys. J. 2015, 799, 220. [Google Scholar] [CrossRef]
Donner, R.; Thiel, M. Scale-resolved phase coherence analysis of hemispheric sunspot activity: A new look at the north-south asymmetry. Astron. Astrophys. 2007, 475, L33–L36. [Google Scholar] [CrossRef]
Temmer, M.; Rybák, J.; Bendík, P.; Veronig, A.; Vogler, F.; Otruba, W.; Pötzi, W.; Hanslmeier, A. Hemispheric sunspot numbers and from 1945–2004: Catalogue and NS asymmetry analysis for solar cycles 18–23. Astron. Astrophys. 2006, 447, 735–743. [Google Scholar] [CrossRef] [Green Version]
Hathaway, D.H.; Wilson, R.M. What the sunspot record tells us about space climate. Sol. Phys. 2004, 224, 5–19. [Google Scholar] [CrossRef] [Green Version]
Zolotova, N.V.; Ponyavin, D.I.; Arlt, R.; Tuominen, I. Secular variation of hemispheric phase differences in the solar cycle. Astron. Nachr. 2010, 331, 765–771. [Google Scholar] [CrossRef] [Green Version]
McIntosh, S.W.; Leamon, R.J.; Gurman, J.B.; Olive, J.P.; Cirtain, J.W.; Hathaway, D.H.; Burkepile, J.; Miesch, M.; Markel, R.S.; Sitongia, L. Hemispheric asymmetries of solar photospheric magnetism: Radiative, particulate, and heliospheric impacts. Astrophys. J. 2013, 765, 146. [Google Scholar] [CrossRef] [Green Version]
Norton, A.A.; Charbonneau, P.; Passos, D. Hemispheric coupling: Comparing dynamo simulations and observations. In The Solar Activity Cycle; Springer: Berlin/Heidelberg, Germany, 2015; pp. 251–283. [Google Scholar] [CrossRef] [Green Version]
Blanter, E.; Le Mouël, J.L.; Shnirman, M.; Courtillot, V. Reconstruction of the North–South Solar Asymmetry with a Kuramoto Model. Sol. Phys. 2017, 292, 54. [Google Scholar] [CrossRef]
Hazra, G.; Choudhuri, A.R. Explaining the variation of the meridional circulation with the solar cycle. Proc. Int. Astron. Union 2018, 13, 313–316. [Google Scholar] [CrossRef] [Green Version]
Zhao, J.; Bogart, R.S.; Kosovichev, A.G.; Duvall, T.L., Jr.; Hartlep, T. Detection of equatorward meridional flow and evidence of double-cell meridional circulation inside the Suna. Astrophys. J. Lett. 2013, 774, L29. [Google Scholar] [CrossRef]
Chen, R.; Zhao, J. A comprehensive method to measure solar meridional circulation and the center-to-limb effect using time–distance helioseismology. Astrophys. J. 2017, 849, 144. [Google Scholar] [CrossRef]
Böning, V.G.; Roth, M.; Jackiewicz, J.; Kholikov, S. Inversions for deep solar meridional flow using spherical born kernels. Astrophys. J. 2017, 845, 2. [Google Scholar] [CrossRef]
Choudhuri, A.R.; Schussler, M.; Dikpati, M. The solar dynamo with meridional circulation. Astron. Astrophys. 1995, 303, L29–L32. [Google Scholar]
Brun, A.S.; Browning, M.K. Magnetism, dynamo action and the solar-stellar connection. Living Rev. Sol. Phys. 2017, 14, 4. [Google Scholar] [CrossRef] [PubMed]
Mandal, K.; Hanasoge, S.M.; Rajaguru, S.P.; Antia, H.M. Helioseismic Inversion to Infer the Depth Profile of Solar Meridional Flow Using Spherical Born Kernels. Astrophys. J. 2018, 863, 39. [Google Scholar] [CrossRef]
Turck-Chieze, S.; Couvidat, S. Solar neutrinos, helioseismology and the solar internal dynamics. Rep. Prog. Phys. 2011, 74, 086901. [Google Scholar] [CrossRef] [Green Version]
Brun, A.S.; Browning, M.K.; Dikpati, M.; Hotta, H.; Strugarek, A. Recent advances on solar global magnetism and variability. Space Sci. Rev. 2015, 196, 101–136. [Google Scholar] [CrossRef]
Acebrón, J.A.; Bonilla, L.L.; Vicente, C.J.P.; Ritort, F.; Spigler, R. The Kuramoto model: A simple paradigm for synchronization phenomena. Rev. Mod. Phys. 2005, 77, 137. [Google Scholar] [CrossRef] [Green Version]
Cumin, D.; Unsworth, C.P. Generalising the Kuramoto model for the study of neuronal synchronisation in the brain. Phys. D Nonlinear Phenom. 2007, 226, 181–196. [Google Scholar] [CrossRef]
Sadilek, M.; Thurner, S. Physiologically motivated multiplex Kuramoto model describes phase diagram of cortical activity. Sci. Rep. 2015, 5, 10015. [Google Scholar] [CrossRef] [Green Version]
Rodrigues, F.A.; Peron, T.K.; Ji, P.; Kurths, J. The Kuramoto model in complex networks. Phys. Rep. 2016, 610, 1–98. [Google Scholar] [CrossRef] [Green Version]
Savostianov, A.; Shapoval, A.; Shnirman, M. Reconstruction of the coupling between solar proxies: When approaches based on Kuramoto and Van der Pol models agree with each other. Commun. Nonlinear Sci. Numer. Simul. 2020, 83, 105149. [Google Scholar] [CrossRef]
Mininni, P.D.; Gomez, D.O.; Mindlin, G.B. Simple model of a stochastically excited solar dynamo. Sol. Phys. 2001, 201, 203–223. [Google Scholar] [CrossRef]
Lopes, I.; Passos, D.; Nagy, M.; Petrovay, K. Oscillator Models of the Solar Cycle. In Space Sciences Series of ISSI; Springer: New York, NY, USA, 2015; pp. 535–559. [Google Scholar] [CrossRef]
Blanter, E.; Le Mouël, J.L.; Shnirman, M.; Courtillot, V. Long Term Evolution of Solar Meridional Circulation and Phase Synchronization Viewed Through a Symmetrical Kuramoto Model. Sol. Phys. 2018, 293, 134. [Google Scholar] [CrossRef]
Hazra, G.; Karak, B.B.; Choudhuri, A.R. Is a deep one-cell meridional circulation essential for the flux transport solar dynamo? Astrophys. J. 2014, 782, 93. [Google Scholar] [CrossRef] [Green Version]
Passos, D.; Charbonneau, P.; Miesch, M. Meridional circulation cynamics from 3D magnetohydrodynamic global simulations of solar convection. Astrophys. J. 2015, 800, L18. [Google Scholar] [CrossRef]
Muñoz Jaramillo, A.; Sheeley, N.R. MWO Polar Faculae Count Calibrated to WSO Polar Fields and SOHO/MDI Polar Flux; Harvard Dataverse: Cambridge, MA, USA, 2016. [Google Scholar]
Blanter, E.M.; Le Mouël, J.L.; Shnirman, M.G.; Courtillot, V. Kuramoto model of nonlinear coupled oscillators as a way for understanding phase synchronization: Application to solar and geomagnetic indices. Sol. Phys. 2014, 289, 4309–4333. [Google Scholar] [CrossRef]
Kuznetsov, A.P.; Stankevich, N.V.; Turukina, L.V. Coupled van der Pol-Duffing oscillators: Phase dynamics and structure of synchronization tongues. Phys. D: Nonlinear Phenom. 2009, 238, 1203–1215. [Google Scholar] [CrossRef]
Astakhov, S.; Gulai, A.; Fujiwara, N.; Kurths, J. The role of asymmetrical and repulsive coupling in the dynamics of two coupled van der Pol oscillators. Chaos Interdiscip. J. Nonlinear Sci. 2016, 26, 023102. [Google Scholar] [CrossRef] [PubMed]
Savostyanov, A.; Shapoval, A.; Shnirman, M. The inverse problem for the Kuramoto model of two nonlinear coupled oscillators driven by applications to solar activity. Phys. D Nonlinear Phenom. 2020, 401, 132160. [Google Scholar] [CrossRef]
Deng, L.H.; Xiang, Y.Y.; Qu, Z.; An, J. Systematic regularity of hemispheric sunspot areas over the past 140 years. Astron. J. 2016, 151, 70. [Google Scholar] [CrossRef]
Syukuya, D.; Kusano, K. Simulation Study of Hemispheric Phase-Asymmetry in the Solar Cycle. arXiv 2016, arXiv:1612.03294. [Google Scholar]
Passos, D.; Nandy, D.; Hazra, S.; Lopes, I. A solar dynamo model driven by mean-field alpha and Babcock-Leighton sources: Fluctuations, grand-minima-maxima, and hemispheric asymmetry in sunspot cycles. Astron. Astrophys. 2014, 563, A18. [Google Scholar] [CrossRef] [Green Version]
Norton, A.A.; Gallagher, J.C. Solar-cycle characteristics examined in separate hemispheres: Phase, gnevyshev gap, and length of minimum. Sol. Phys. 2010, 261, 193. [Google Scholar] [CrossRef] [Green Version]
Gurgenashvili, E.; Zaqarashvili, T.V.; Kukhianidze, V.; Oliver, R.; Ballester, J.L.; Dikpati, M.; McIntosh, S.W. North–South Asymmetry in Rieger-type Periodicity during Solar Cycles 19–23. Astrophys. J. 2017, 845, 137. [Google Scholar] [CrossRef] [Green Version]
Nepomnyashchikh, A.; Mandal, S.; Banerjee, D.; Kitchatinov, L. Can the long-term hemispheric asymmetry of solar activity result from fluctuations in dynamo parameters? Astron. Astrophys. 2019, 625, A37. [Google Scholar] [CrossRef] [Green Version]
Duhau, S.; Chen, C.Y. The sudden increase of solar and geomagnetic activity after 1923 as a manifestation of a non-linear solar dynamo. Geophys. Res. Lett. 2002, 29, 6-1–6-4. [Google Scholar] [CrossRef]
Kane, R.P. Solar cycle predictions based on extrapolation of spectral components: An update. Sol. Phys. 2007, 246, 487–493. [Google Scholar] [CrossRef]
Richards, M.T.; Rogers, M.L.; Richards, D.S.P. Long-term variability in the length of the solar cycle. Publ. Astron. Soc. Pac. 2009, 121, 797. [Google Scholar] [CrossRef] [Green Version]
Blanter, E.; Le Mouël, J.L.; Shnirman, M.; Courtillot, V. Kuramoto model with non-symmetric coupling reconstructs variations of the solar-cycle period. Sol. Phys. 2016, 291, 1003–1023. [Google Scholar] [CrossRef]
Featherstone, N.A.; Miesch, M.S. Meridional circulation in solar and stellar convection zones. Astrophys. J. 2015, 804, 67. [Google Scholar] [CrossRef] [Green Version]
Choudhuri, A.R. A critical assessment of the flux transport dynamo. J. Astrophys. Astron. 2015, 36, 5–14. [Google Scholar] [CrossRef] [Green Version]
Cameron, R.H.; Dikpati, M.; Brandenburg, A. The global solar dynamo. Space Sci. Rev. 2017, 210, 367–395. [Google Scholar] [CrossRef] [Green Version]

Figure 1. Polar faculae data for northern (blue curve) and southern (red curve) hemispheres.

Figure 2. The relationship between coupling

μ

Δ ω

, and the correlation

ρ

between two coupled van der Pol oscillators;

Δ ω < μ

Figure 2. The relationship between coupling

μ

Δ ω

, and the correlation

ρ

between two coupled van der Pol oscillators;

Δ ω < μ

Figure 3. Reconstruction scheme described in details in steps (A)–(C).

Figure 4. Top panel: Reconstructed coupling with the Kuramoto (

{\hat{μ}}_{k} (t)

; green) and van der Pol (

{\hat{μ}}_{v} (t)

; orange) models found with

Δ ω = 0.2

. Bottom panel: ratio

r_{Δ ω} [{\hat{μ}}_{k}, {\hat{μ}}_{v}] (t)

defined by (14) exhibiting the proximity of the reconstructions. Grey figure background means that the reconstruction of the coupling involves either negative correlation between the series (in the middle) or the computation of the correlation on smaller windows (at the left and right);

Δ ω = 0.1, 0.2, 0.3

Figure 4. Top panel: Reconstructed coupling with the Kuramoto (

{\hat{μ}}_{k} (t)

; green) and van der Pol (

{\hat{μ}}_{v} (t)

; orange) models found with

Δ ω = 0.2

. Bottom panel: ratio

r_{Δ ω} [{\hat{μ}}_{k}, {\hat{μ}}_{v}] (t)

Δ ω = 0.1, 0.2, 0.3

Figure 5. The minimal (blue) and maximal (red curve) values of

r_{Δ ω} (t)

obtained with different

Δ ω

vs. time.

Figure 5. The minimal (blue) and maximal (red curve) values of

r_{Δ ω} (t)

obtained with different

Δ ω

vs. time.

Figure 6. Reconstructed normalized natural frequency,

1 + Δ ω

, associated with the signals coming from the northern solar hemisphere. The reconstruction rule is inferred from the assumption that the Kuramoto and VdP reconstruction of the coupling result in the time-independent values of

r_{Δ ω} (t)

fixed to

0.86

Figure 6. Reconstructed normalized natural frequency,

1 + Δ ω

r_{Δ ω} (t)

fixed to

0.86

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Savostianov, A.; Shapoval, A.; Shnirman, M. Dynamics of Phase Synchronization between Solar Polar Magnetic Fields Assessed with Van Der Pol and Kuramoto Models. Entropy 2020, 22, 945. https://doi.org/10.3390/e22090945

AMA Style

Savostianov A, Shapoval A, Shnirman M. Dynamics of Phase Synchronization between Solar Polar Magnetic Fields Assessed with Van Der Pol and Kuramoto Models. Entropy. 2020; 22(9):945. https://doi.org/10.3390/e22090945

Chicago/Turabian Style

Savostianov, Anton, Alexander Shapoval, and Mikhail Shnirman. 2020. "Dynamics of Phase Synchronization between Solar Polar Magnetic Fields Assessed with Van Der Pol and Kuramoto Models" Entropy 22, no. 9: 945. https://doi.org/10.3390/e22090945

APA Style

Savostianov, A., Shapoval, A., & Shnirman, M. (2020). Dynamics of Phase Synchronization between Solar Polar Magnetic Fields Assessed with Van Der Pol and Kuramoto Models. Entropy, 22(9), 945. https://doi.org/10.3390/e22090945

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu