Research Article
Towards an algebraic method of solar cycle prediction
I. Calculating the ultimate dipole contributions of individual active regions
^{1}
Department of Astronomy, Eötvös Loránd University, 1053 Budapest, Hungary
^{2}
Department of Mathematical Sciences, Durham University, Science Laboratories, DH1 3LE Durham, UK
^{*} Corresponding author: K.Petrovay@astro.elte.hu
Received:
28
May
2020
Accepted:
31
August
2020
We discuss the potential use of an algebraic method to compute the value of the solar axial dipole moment at solar minimum, widely considered to be the most reliable precursor of the activity level in the next solar cycle. The method consists of summing up the ultimate contributions of individual active regions to the solar axial dipole moment at the end of the cycle. A potential limitation of the approach is its dependence on the underlying surface flux transport (SFT) model details. We demonstrate by both analytical and numerical methods that the factor relating the initial and ultimate dipole moment contributions of an active region displays a Gaussian dependence on latitude with parameters that only depend on details of the SFT model through the parameter η/Δ_{u} where η is supergranular diffusivity and Δ_{u} is the divergence of the meridional flow on the equator. In a comparison with cycles simulated in the 2 × 2D dynamo model we further demonstrate that the inaccuracies associated with the algebraic method are minor and the method may be able to reproduce the dipole moment values in a large majority of cycles.
Key words: solar cycle / rogue sunspots / surface flux transport modeling
© K. Petrovay et al., Published by EDP Sciences 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Predicting the amplitude of an upcoming solar cycle is the central issue of space climate forecasting. It is widely accepted that the best performing, physically well motivated prediction method is based on the good linear correlation between the solar axial dipole moment in solar activity minimum and the amplitude of the next cycle (Schatten et al., 1978; Wang & Sheeley, 2009; MuñozJaramillo et al., 2013; Hathaway & Upton, 2016; Petrovay, 2020). In order to extend the rather short, 3–4 year long time span of these forecasts, in recent years efforts have been made to “predict the predictor”, i.e. to model and forecast the evolution of the magnetic flux distribution over the solar surface (Jiang et al., 2018).
Starting with the pioneering work of DeVore et al. (1985), the standard approach to this problem has been the use of surface flux transport (SFT) simulations. This approach assumes that the lineofsight magnetic field component shown in synoptic maps constructed from solar magnetograms corresponds to the projection of an inherently radial mean photospheric magnetic field with strength B_{r}, the transport of which is governed by advection due to large scale flows and diffusion due to supergranular motions:
(1)where t is time, λ and ϕ are heliographic latitude and longitude, R is the radius of the Sun, Ω is the angular velocity of differential rotation, u(λ) = u_{0} f(λ) is the meridional flow, and η is the supergranular diffusivity. The source term S_{r} represents the emergence of new flux into the atmosphere in active regions, while the term −B_{r}/τ is the simplest (though not the most realistic) form of a sink term representing the effects of radial diffusion that would appear in a consistent derivation of the transport equation from the radial component of the induction equation. For further discussion of the decay term see e.g. Baumann et al. (2006), Whitbread et al. (2017) or Petrovay & Talafha (2019). For general reviews of the SFT modelling approach see Sheeley (2005), Mackay & Yeates (2012), Jiang et al. (2014b), Wang (2017)).
Numerous studies have been performed with the objective to reproduce and predict the evolution of the solar dipole moment (Wang & Sheeley, 1991; Virtanen et al., 2017; Whitbread et al., 2017). Despite some spectacular successes, this approach still suffers from two main limitations:

Ill constrained SFT model ingredients. The model has at least 3 free numerical parameters (u_{0}, η, τ) as well as a free function, the observationally not well mapped meridional flow profile f(λ). Attempts to determine the ingredients from direct observations have dubious relevance and often yield contradictory results. Internal optimizations of the model, in turn, result in parameter combinations that may depend on the choice of merit and still allow a rather wide freedom in the choices (Lemerle et al., 2015; Virtanen et al., 2017; Whitbread et al., 2017; Petrovay & Talafha, 2019). The ill constrained ingredients imply that in order to obtain a realistic estimate of the uncertainities in the predictions an ensemble of models with varying parameters should be studied; with the need to numerically solve the partial differential equation (1) for each of these models the process can get rather lengthy and cumbersome.

A realistic representation of the active region source. These sources are often represented as simple bipoles instantaneously introduced into the simulation, or, more realistically, by assimilating actual magnetograms taken at the time of their central meridian passage. In reality, however, during the evolution of an active region its magnetic flux distribution keeps changing as a consequence of the emergence of new flux and the proper motions of sunspots and other magnetic flux concentrations driven by subsurface dynamics and/or by random, localized photospheric flows – effects which are not accounted for in the SFT model. The choice of the proper form of the source and the time of its introduction is therefore a highly nontrivial problem. And the very large number (thousands) of active regions arising in a solar cycle makes this task challenging even in the simplest, idealized case of bipole representation. These problems are even further aggravated in the case of historical data, when the objective is to understand and reproduce the course of evolution of solar activity in centuries past.
The objective of this paper is to consider ways to alleviate these difficulties by simplifying the prediction method to its bare essentials. Specifically, in Section 2 we will show that the SFT modelling approach of solving a partial differential equation can be reduced to the calculation of an algebraic sum. Sections 3 and 4 demonstrate that the result of this procedure only depends on two numerical combinations of the model parameters without the need to specify the choice of any unknown function. Then, in Section 5 we validate the approach in a comparison with activity cycles simulated in a dynamo model. These results may pave the way towards a more robust and effective approach to solar cycle prediction.
2 Mathematical formulation of the problem
The axial dipolar moment of the Sun is given by
(2)where is the azimuthally averaged field strength. Throughout this paper, D will denote the global dipole moment of the whole Sun as defined in equation (2), while δD will denote the contribution to D from an individual active region. The evolution of B, and consequently D, is determined by the 1D surface flux transport equation obtained by azimuthally integrating equation (1):
In this azimuthally averaged representation a tilted bipolar region is then a pair of flux rings of opposite polarity, appearing as a bipole source with a finite latitudinal separation in equation (3).
Figure 1 presents solutions of equation (3) for a particular case where the source is a single bipole placed at λ_{0} = +15° latitude at time t = 0. Both polarities were assumed to initially have a Gaussian profile in λ with a halfwidth σ_{0} = 6°. The lefthand panels show the evolution of the axial dipolar moment. The righthand panels display the evolution of B in the timelatitude plane. The cases shown in the top and bottom rows only differ in the value of τ: in the first case τ is effectively infinite (no decay), while in the second case it is shorter than the solar cycle period. As it is easy to understand from the structure of equation (3), the two solutions only differ by the presence of an exponential factor exp(−t/τ) when τ is finite.
Fig. 1 Solutions of the 1D SFT equation with a bipolar source placed at λ_{0} = +15° latitude at time t = 0. Units of the model parameters shown on top of the rows are m/s, km^{2}/s and year, respectively. The meridional flow profile is given by equation (9). Dashed blue lines in the lefthand plots mark the asymptotic solutions D = δD_{∞} (top) and D = δD_{∞}exp(−t/τ) (bottom). 
The evolution is determined by the competition of two processes. The diffusive spreading of the two polarity patches of opposite sign leads to the cancellation of a large part of the flux originally present, yet a small fraction of the flux still manages to reach the Southern hemisphere where it is transported to the South pole by the meridional flow. The “leading” polarity flux patch, situated closer to the equator, gives a larger contribution to the flux ultimately reaching the South pole, so in the final state, a leading polarity patch remains at the South pole, while a corresponding trailing polarity patch remains at the North pole. While the flux in these patches is a fraction of the original flux in the region, their high latitudinal separation gives rise to a dipole moment that, in the case plotted in Figure 1, exceeds the initial value. In the limit τ → ∞ the dipole moment remains very nearly constant at a fixed value δD_{∞}(λ_{0}). For finite τ the dipole moment asymptotically behaves as δD_{∞}exp(−t/τ).
This implies that if, ignoring the complex details of its structure and evolution, the ith active region in cycle n (starting at time t_{n}) is represented by a simple dipole instantaneously introduced into the SFT model at time t_{i} with an initial dipole moment δD_{1,i}, the ultimate contribution of all active regions at the end of the cycle will be given by
(4)where N_{tot} is the total number of ARs in the cycle, δD_{U} is the ultimate contribution of an active region to the global dipole moment at the end of a cycle, and the asymptotic dipole contribution factor is
From equation (2) it is straightforward to show that for a pointlike bipole consisting of a pair of pointlike polarities with a small latitudinal separation d_{λ} the initial dipole moment contribution δD_{1} is given by
(6)where Φ is the magnetic flux in the northern flux patch.
Equations (4)–(6) offer a simple algebraic tool to extend the temporal scope of the polar field precursor method by “predicting the precursor”, i.e. computing the dipole moment built up during a solar activity cycle without the need to solve the partial differential equation (3) or (1). For the case of cycle 23 this approach has been already exploited by Jiang et al. (2019) for one particular SFT model setup.
Generally, however, this approach is subject to a number of limitations. These limitations were already outlined in the Section 1 above. In more detail, they are:

(1a) The result depends on details of the SFT model used, mainly through the f_{∞,i}, and for models with a decay term also through the exponential factor in (4).

(1b) As illustrated in Figure 1, the asymptotic solution is approximated after a transitional period of 2–3 years only. The contribution of active regions emerging in the last years of the cycle is therefore not accurately represented by this formula. As, however, activity in this late phase of the cycle is normally rather low, this is not expected to be a major limitation in most cases.

(2a) The assumption that active regions may be represented by the instantaneous introduction of simple bipoles into the simulation is undoubtedly a strong simplification.

(2b) The number N_{tot} of terms in the sum, i.e. the number of bipolar regions contributing to the solar axial dipole moment can be quite high.
In the present work we focus on issues (1a), (1b) and (2a). Issue (2b) will be dealt with in a companion paper.
3 Calculating the ultimate dipole contribution of active regions
The dependence of the asymptotic dipole contribution factor f_{∞} on latitude was first considered by Jiang et al. (2014a). In a series of numerical experiments with one particular SFT model they found a Gaussian dependence on latitude:
In what follows, λ_{R} will be referred to as the dynamo effectivity range of active regions.
Let us consider whether this conclusion holds generally for other SFT model setups. In Figure 2 the results of a similar set of experiments for different model setups are shown: in addition to Jiang et al.’s original setup, the other SFT models used were those of Lemerle et al. (2015), Lemerle & Charbonneau (2017), Whitbread et al. (2017) and Wang (2017). It is apparent that the Gaussian dependence identified by Jiang et al. (2014a) holds generally. Indeed, even in an SFT model where active region shapes are determined by assimilation rather than fitted with bipoles, Whitbread et al. (2018) still found a Gaussian. In all these cases, however, the width and amplitude of the Gaussian are different. It was indeed already found by Nagy et al. (2017) that the dynamo effectivity range in the case of the Lemerle & Charbonneau (2017) setup was significantly wider than expected from the results of Jiang et al. (2014a).
Fig. 2 Dependence of the asymptotic dipole contribution factor f_{∞} of bipolar sources on their latitude in various SFT model setups. Solid lines are Gaussian fits to the numerical results. 
In order to understand the dependence of λ_{R} and A on model parameters we determine these parameters on a large grid of SFT models. Our model grid is essentially identical to the grid considered by Petrovay & Talafha (2019), except that here we limit ourselves to models with effectively no decay term (τ = 1000 years) as the effect of τ has already been separated in the exponential factor in equation (4). Four types of flow profiles are considered (cf. Fig. 3).
Fig. 3 The meridional flow profiles used in the paper. 
Flow 1: a simple sinusoidal profile
Flow 2: a sinusoidal profile with a dead zone around the poles
Flow 3: The Lemerle & Charbonneau (2017) profile peaking at high latitudes,
Flow 4: A profile peaking at very low latitudes, considered by Wang (2017):
For each of the profiles, u_{0} was allowed to vary between 5 and 20 m/s in steps of 2.5, while η varied from 50 to 750 km^{2}/s in steps of 50. From numerical runs like the one plotted in Figure 1 the values of λ_{R} and A were determined for each model.
Experimenting with different simple combinations of the input and output parameters we find the clearest relationship in the case plotted in Figure 4. Here,
(12)is the divergence of the meridional flow at the equator.
Fig. 4 Dependence of the dynamo effectivity range λ_{R} (left) and the amplitude A (right) of the asymptotic dipole contribution factor of bipolar sources on selected parameter combinations in various SFT model setups for σ_{0} = 6°. In the lefthand panel the solid line shows the analytic result (22) in the lowlatitude Cartesian limit, while the dashed line shows a fiducial analytic fit of the type (26). In the righthand panel the shaded region is the range expected for a sin^{ n } λ equilibrium polar field profile with n > 7; the horizontal line marks the value for n = 8. 
The finding that the single parameter combination η/Δ_{u} determines f_{∞} for all but one of our 2parameter model grids, irrespective of the choice of the meridional flow profile is an impressive and somewhat unexpected result which calls for a theoretical interpretation.
4 Analytical derivation of the asymptotic dipole contribution factor
4.1 Lowlatitude Cartesian limit
To make the problem analytically tractable, we limit ourselves to the neighbourhood of the equator (λ ≪ 1 radian) where the spherical coordinate grid may be locally approximated by a Cartesian setup and meridional flow is expressed by the leading term of its Taylor expansion as u = RΔ_{u} λ with Δ_{u} = const. This is formally identical to the cosmological case of a onedimensional Hubble flow in a vacuum dominated universe and the advection of a frozenin magnetic field configuration by the flow is an exponential expansion. This cosmological analogy suggests to consider the evolution in the Lagrangian comoving (coexpanding) frame, where fluid elements are labelled by their latitude λ_{L} at time zero, rather than their current latitude .
We recall the initial condition of the evolution, as illustrated in Figure 1: a pair of opposite polarity flux rings with Gaussian profile of half width σ_{0} ≪ 1 radian, situated at latitude λ_{0}, with a separation d_{λ} ≪ 1 radian between the rings. (In the diffusive case, other assumed initial profiles will also soon approach Gaussian by virtue of the central limit theorem.) What we are looking for is the amount of net transequatorial flux (flux in the other hemisphere) in the limit t → ∞.
4.1.1 Asymptotic magnetic field profile
In the Lagrangian frame flow advection is absent by definition, so the flux transport equation simplifies to a diffusion equation
(13)where, however, the diffusivity η_{L}(t) is now time dependent. Indeed, in the comoving frame the unit of length expands exponentially as exp(Δ_{u} t), so the diffusivity, of dimension length^{2}/time, expressed in these units, will scale as . For the same reason, the Lagrangian flux density B_{L} is related to the Eulerian by .
Consider the evolution of one of the two flux patches comprising the bipole. Our initial condition is
This problem may be solved exactly using Fourier transforms (cf. Mackay et al., 2016). First, we change the time variable from t to to obtain the standard diffusion equation
(15)which may be solved using standard techniques. In particular, if we define the Fourier transform
(16)then equation (15) implies that . The Fourier transform of our initial condition is
Inverting the transform finally gives
Since , we arrive at
Note that this selfsimilar solution might have been anticipated, given that the Gaussian is known to be a selfsimilar solution of the diffusion equation and the meaning of B_{L} is (onedimensional) flux density, so magnetic flux conservation requires the amplitude of the Gaussian to scale with the inverse of σ. Plugging equation (19) back into (13) then returns (20).
4.1.2 Transequatorial flux
In the t → ∞ limit equation (19) gives
Using equation (21), and taking λ_{0} > 0 for concreteness, the fraction of flux of this single polarity in the opposite (Southern) hemisphere is given by
For our pair of flux rings separated by d_{λ}, the net transequatorial flux fraction is then
(24)where the last form is the leading term in a Taylor expansion for small d_{λ}.
4.2 Sphericity effects
To compute the asymptotic dipole contribution δD_{∞} from Φ_{∞} by equation (2) we need to return to spherical geometry. The poleward meridional flow results in a “topknot” equilibrium field distribution strongly peaked near the poles (Sheeley et al., 1989). Indeed, approximating the field profile as B(λ) ∝ sin^{ n } λ, even flows mildly concentrated on the poles will result in n ≃ 7 (cf. Fig. 4 in Wang, 2017). Observational constraints indicate n ≳ 8–9 (Petrie, 2015). With this approximation, using equations (2), (5), (6) and (24), we have the following expression for the asymptotic dipole contribution factor in the lowlatitude limit:
(A factor 1/cosλ_{0} originating from (6) has been omitted as in the lowlatitude limit it becomes unity.) For n = 8 this yields a = 41°.16; values for n > 7 are in the range between 40°.66 and 45°.7. The curves in the righthand panel of Figure 4 are in agreement with this result.
The curves representing flow types 1–3 in the lefthand panel of Figure 4 are well fitted by the solution (22) at low values of λ_{R}. This is to be expected as these flows peak at latitudes above 40° so for low latitudes the profiles are well approximated as linear. For the same reason, while for values λ_{R} ~ 10–20° curvature effects do come into play, the curves still do not diverge as the nonlinearity of the flow profile remains low, hence, with an appropriate planar projection of the spherical surface, the flow may still be transformed out switching to a homologously expanding Lagrangian frame. In this case equation (13) generalizes to the diffusion equation on a spherical surface which has no fluxconserving solution with an exactly Gaussian profile, though Figure 2 indicates that at moderate latitudes the solutions may still be well approximated by a Gaussian. The dynamo effectivity range, however, changes. A good empirical fit to the curves is found to be
The points representing flow type 4 in the lefthand panel of Figure 4 strongly diverge from both equations (22) and (26). The reason is that for this profile peaking at very low latitude, the nonlinearity of the profile becomes important already at λ ~ σ_{0}. In effect, in most of the area covered by the AR flux during its evolution the expansion rate will be far below the nominal equatorial value – at λ > 13° the surface will actually already contract, strengthening the field. Hence, nominal values of the parameter η/Δ_{u} are not really relevant for the determination of λ_{R} in this case.
To close off our analytic discussion we note that the expression (22) for the dynamo effectivity range can also be understood in simple physical terms. The time scale associated with diffusion to the other hemisphere from latitude λ is (Rλ)^{2}/η. The latitude where this equals the advective time scale 1/Δ_{u} is just (η/R ^{2}Δ_{u})^{1/2}: for higher latitudes diffusion cannot compete with the poleward advection and little flux from here can reach the other hemisphere.
5 Comparison to a dynamo model solution
Our suggested algebraic approach to solar cycle prediction consists in using equation (4) to calculate the net dipole moment at the end of a solar cycle, where the dynamo effectivity f_{∞} is given by equation (25), with λ_{R} and a taken from either a direct interpolation of the numerical results plotted in Figure 4, from their analytical approximation (dashed curve) or even from the lowlatitude limit (22).
In order to test the validity and accuracy of the suggested algebraic approach, in Figure 5 we compare the results of a run of the 2 × 2D dynamo model, as described in Lemerle & Charbonneau (2017), with the results of our algebraic approach. For the SFT parameters used in the Lemerle & Charbonneau (2017) model (, u_{0} = 17 m/s, η = 600 km^{2}/s, τ = 10 yr) the numerical results plotted in Figure 4 give λ_{R} = 13°.6 and A = 3.75. The advantage of using the 2 × 2D model is that it explicitly includes a 2D SFT model as one of its components, with the source term represented as idealized bipoles with parameters that can be directly extracted from the models. In this way, factors like an arbitrary choice of model parameters or illspecified sources will not influence the comparison.
Fig. 5 Comparison of total net change in the solar axial dipole moment D_{n+1} − D_{n} in the 2 × 2D dynamo simulation and its approximations with the algebraic method in a segment of 80 cycles. Plotted are absolute values of the quantities shown in the legend. See text for further explanations. 
Plotted in blue is the net change in the axial dipole moment D_{n+1} − D_{n} between subsequent minima of the cycle in the model. The red dashed curve shows the values computed by adding to the actual dipole moment at the start of a cycle the expected total dipole moment contribution calculated by the algebraic method, equations (4) and (25), with the parameter values quoted above; δD_{1,i} values are computed using the BMR properties extracted from the dynamo code. While the overall agreement is quite good, some smaller discrepancies remain. The standard deviation of the relative error is 10.1%. Sources contributing to this unexplained variance include the invalidity of formula (4) for active regions emerging in the last three years of each cycle: as illustrated in Figure 1, in these cases the time elapsed from emergence to solar minimum is too short for the asymptotic solution to set in. Other contributing factors are smaller differences in model details between the 2 × 2D model and the 1D SFT model forming the basis of our algebraic approach.
The green dashed curve, in turn, diplays the result of the algebraic method for the “reduced stochasticity” case. In this case, bipolar magnetic regions (as the active regions are called in the 2 × 2 D model) are substituted with regions of the same flux and latitude, but with tilt and polarity separation values corresponding to the expected values for the given flux and latitude, as given by equations (15) and (16a) in Lemerle et al. (2015). While the agreement is noticeably poorer (standard deviation of the relative error now reaches 21.2%), the net dipole moment change is still reproduced with less than 25% error in about 90% of the cases. This indicates that detailed knowledge of the structure and evolution of each individual active region may not be indispensable for a tolerably good predictive skill in the algebraic method in most (though not all) cycles. This issue will be discussed further in the second paper of this series (Nagy et al., 2020).
6 Conclusions
We have discussed the potential use of an algebraic method to predict the value of the solar axial dipole moment at solar minimum. The method, already applied in the case of one particular SFT model by Jiang et al. (2019), consists in summing up the ultimate contributions of individual active regions the solar axial dipole moment as given by equations (4) and (25).
In Section 2 we listed four potential limitations of the approach. The first of these, (1a) was its dependence on SFT model details. Indeed, the meridional flow profile is still rather uncertain and potentially variable; in addition, systematic inflows towards the activity belts are also expected to impact on flux transport, see Nagy et al. (2020, in this issue). Disregarding time dependent effects here we demonstrated by both analytical and numerical methods that the dynamo effectivity range λ_{R} and the equatorial value of the asymptotic dipole contribution factor A only depend on details of the SFT model through the parameter η/Δ_{u}. This significantly reduces the uncertainty introduced by the choice of model details and makes the algebraic method preferable to the more computationintensive traditional method of numerically solving the SFT equation.
While numerical costs and a more limited freedom in the choice of parameters advocate the use of this algebraic approach, this clearly goes at the cost of accuracy. One source of the inaccuracy, (1b), is inapplicability of the ultimate dipole contribution to active regions appearing in the last years of a solar cycle, as the contributions of such ARs have not yet reached equilibrium at the time of minimum. In a comparison with cycles simulated in the 2 × D dynamo model we demonstrated that the inaccuracy introduced by this effect and by differences in the underlying numerical models is small and does not constitute a great obstacle in the way of correctly reproducing the dipole moment.
Another source of inaccuracy of this approach, (2a), is that the representation of active region sources by idealized bipoles is likely to be far from perfect. Applying our method to the same dynamosimulated cycles but substituting the assumed initial AR dipole moments with their expected values for an AR of the given size and latitude showed that in this case the dipole moment could still be reasonably well reproduced in the large majority of cycles, lending further support to the algebraic approach. Nevertheless, in a small fraction of cycles inaccurate representation of the sources does lead to significant inaccuracies in the resulting dipole moment values. Note that while in the dynamo model used here for comparison ARs were assumed to be be bipolar, in applications to solar data a further source of uncertainty concerns to what extent a bipolar representation reflects the structure of ARs (cf. Iijima et al., 2019; Jiang et al., 2019; Yeates, 2020).
The fourth limitation of the method, (2b) is related to the very high number of terms in the summation that are theoretically needed for the correct representation of the dipole contributions, which exacerbates the issue with the correct representation of the initial AR contributions. On the other hand, the fact that there are many regions might help fluctuations from the Gaussian trend to average out in an overall prediction. This issue will be discussed further in the second part of this series.
Acknowledgments
This research was supported by the Hungarian National Research, Development and Innovation Fund (grant no. NKFI K128384), by the UK STFC (grant no. ST/S000321/1) and by the European Union’s Horizon 2020 research and innovation programme under grant agreement no. 739500. The collaboration of the authors was facilitated by support from the International Space Science Institute in ISSI Team 474. The editor thanks two anonymous reviewers for their assistance in evaluating this paper.
References
 Baumann I, Schmitt D, Schüssler M. 2006. A necessary extension of the surface flux transport model. A&A 446: 307–314. https://doi.org/10.1051/00046361:20053488. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 DeVore CR, Sheeley NR Jr, Boris JP, Young TR Jr, Harvey KL. 1985. Simulations of magneticflux transport in solar active regions. Sol. Phys. 102(1–2): 41–49. https://doi.org/10.1007/BF00154036. [NASA ADS] [CrossRef] [Google Scholar]
 Hathaway DH, Upton LA. 2016. Predicting the amplitude and hemispheric asymmetry of solar cycle 25 with surface flux transport. J Geophys Res (Space Phys) 121(11): 10744–10753. https://doi.org/10.1002/2016JA023190. [CrossRef] [Google Scholar]
 Iijima H, Hotta H, Imada S. 2019. Effect of morphological asymmetry between leading and following sunspots on the prediction of solar cycle activity. Astrophys J 883(1): 24. https://doi.org/10.3847/15384357/ab3b04. [CrossRef] [Google Scholar]
 Jiang J, Cameron RH, Schüssler M. 2014a. Effects of the scatter in sunspot group tilt angles on the largescale magnetic field at the solar surface. Astrophys J 791: 5. https://doi.org/10.1088/0004637X/791/1/5. [NASA ADS] [CrossRef] [Google Scholar]
 Jiang J, Hathaway DH, Cameron RH, Solanki SK, Gizon L, Upton L. 2014b. Magnetic flux transport at the solar surface. Space Sci Rev 186: 491–523. https://doi.org/10.1007/s1121401400831. [NASA ADS] [CrossRef] [Google Scholar]
 Jiang J, Song Q, Wang JX, Baranyi T. 2019. Different contributions to space weather and space climate from different big solar active regions. Astrophys J 871(1): 16. https://doi.org/10.3847/15384357/aaf64a. [CrossRef] [Google Scholar]
 Jiang J, Wang JX, Jiao QR, Cao JB. 2018. Predictability of the solar cycle over one cycle. Astrophys J 863: 159. https://doi.org/10.3847/15384357/aad197. [NASA ADS] [CrossRef] [Google Scholar]
 Lemerle A, Charbonneau P. 2017. A coupled 2 × 2D BabcockLeighton solar dynamo model. II. Reference dynamo solutions. Astrophys J 834: 133. https://doi.org/10.3847/15384357/834/2/133. [NASA ADS] [CrossRef] [Google Scholar]
 Lemerle A, Charbonneau P, CarignanDugas A. 2015. A coupled 2 × 2D BabcockLeighton solar dynamo model. I. Surface magnetic flux evolution. Astrophys J 810(1): 78. https://doi.org/10.1088/0004637X/810/1/78. [NASA ADS] [CrossRef] [Google Scholar]
 Mackay DH, Yeates AR. 2012. The Sun’s global photospheric and coronal magnetic fields: Observations and models. Living Rev Sol Phys 9(1): 6. https://doi.org/10.12942/lrsp20126. [CrossRef] [Google Scholar]
 Mackay DH, Yeates AR, Bocquet FX. 2016. Impact of an L5 magnetograph on nonpotential solar global magnetic field modeling. Astrophys J 825(2): 131. https://doi.org/10.3847/0004637X/825/2/131. [NASA ADS] [CrossRef] [Google Scholar]
 MuñozJaramillo A, DasiEspuig M, Balmaceda LA, DeLuca EE. 2013. Solar cycle propagation, memory, and prediction: Insights from a century of magnetic proxies. Astrophys J Lett 767: L25. https://doi.org/10.1088/20418205/767/2/L25. [NASA ADS] [CrossRef] [Google Scholar]
 Nagy M, Lemerle A, Labonville F, Petrovay K, Charbonneau P. 2017. The effect of “Rogue” active regions on the solar cycle. Sol Phys 292: 167. https://doi.org/10.1007/s1120701711940. [NASA ADS] [CrossRef] [Google Scholar]
 Nagy M, Petrovay K, Lemerle A, Charbonneau P. 2020. Towards an algebraic method of solar cycle prediction II. Reducing the need for detailed input data with ARDOR. J Space Weather Space Clim. This issue. https://doi.org/10.1051/swsc/2020051. [Google Scholar]
 Petrie GJD. 2015. Solar magnetism in the polar regions. Living Rev Sol Phys 12: 5. https://doi.org/10.1007/lrsp20155. [CrossRef] [Google Scholar]
 Petrovay K. 2020. Solar cycle prediction. Living Rev Sol Phys 17(1): 2. https://doi.org/10.1007/s411160200022z. [CrossRef] [Google Scholar]
 Petrovay K, Talafha M. 2019. Optimization of surface flux transport models for the solar polar magnetic field. A&A 632: A87. https://doi.org/10.1051/00046361/201936099. [CrossRef] [EDP Sciences] [Google Scholar]
 Schatten KH, Scherrer PH, Svalgaard L, Wilcox JM. 1978. Using dynamo theory to predict the sunspot number during solar cycle 21. Geophys Res Lett 5: 411–414. https://doi.org/10.1029/GL005i005p00411. [NASA ADS] [CrossRef] [Google Scholar]
 Sheeley NR Jr. 2005. Surface evolution of the sun’s magnetic field: A historical review of the fluxtransport mechanism. Living Rev Sol Phys 2: 5. https://doi.org/10.12942/lrsp20055. [Google Scholar]
 Sheeley NR Jr, Wang YM, DeVore CR. 1989. Implications of a strongly peaked polar magnetic field. Sol Phys 124: 1–13. https://doi.org/10.1007/BF00146515. [NASA ADS] [CrossRef] [Google Scholar]
 Virtanen IOI, Virtanen II, Pevtsov AA, Yeates A, Mursula K. 2017. Reconstructing solar magnetic fields from historical observations. II. Testing the surface flux transport model. A&A 604: A8. https://doi.org/10.1051/00046361/201730415. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wang YM. 2017. Surface flux transport and the evolution of the sun’s polar fields. Space Sci Rev 210: 351–365. https://doi.org/10.1007/s1121401602570. [CrossRef] [Google Scholar]
 Wang YM, Sheeley NR Jr. 1991. Magnetic flux transport and the sun’s dipole moment: New twists to the BabcockLeighton Model. Astrophys J 375: 761. https://doi.org/10.1086/170240. [NASA ADS] [CrossRef] [Google Scholar]
 Wang YM, Sheeley NR. 2009. Understanding the geomagnetic precursor of the solar cycle. Astrophys J Lett 694: L11–L15. https://doi.org/10.1088/0004637X/694/1/L11. [NASA ADS] [CrossRef] [Google Scholar]
 Whitbread T, Yeates AR, MuñozJaramillo A. 2018. How many active regions are necessary to predict the solar dipole moment? Astrophys J 863(2): 116. https://doi.org/10.3847/15384357/aad17e. [NASA ADS] [CrossRef] [Google Scholar]
 Whitbread T, Yeates AR, MuñozJaramillo A, Petrie GJD. 2017. Parameter optimization for surface flux transport models. A&A 607: A76. https://doi.org/10.1051/00046361/201730689. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Yeates AR. 2020. How good is the bipolar approximation of active regions for surface flux transport? Sol Phys 295: 119. https://doi.org/10.1007/s1120702001688y. [CrossRef] [Google Scholar]
Cite this article as: Petrovay K, Nagy M & Yeates A 2020. Towards an algebraic method of solar cycle prediction. J. Space Weather Space Clim. 10, 50. https://doi.org/10.1051/swsc/2020050.
All Figures
Fig. 1 Solutions of the 1D SFT equation with a bipolar source placed at λ_{0} = +15° latitude at time t = 0. Units of the model parameters shown on top of the rows are m/s, km^{2}/s and year, respectively. The meridional flow profile is given by equation (9). Dashed blue lines in the lefthand plots mark the asymptotic solutions D = δD_{∞} (top) and D = δD_{∞}exp(−t/τ) (bottom). 

In the text 
Fig. 2 Dependence of the asymptotic dipole contribution factor f_{∞} of bipolar sources on their latitude in various SFT model setups. Solid lines are Gaussian fits to the numerical results. 

In the text 
Fig. 3 The meridional flow profiles used in the paper. 

In the text 
Fig. 4 Dependence of the dynamo effectivity range λ_{R} (left) and the amplitude A (right) of the asymptotic dipole contribution factor of bipolar sources on selected parameter combinations in various SFT model setups for σ_{0} = 6°. In the lefthand panel the solid line shows the analytic result (22) in the lowlatitude Cartesian limit, while the dashed line shows a fiducial analytic fit of the type (26). In the righthand panel the shaded region is the range expected for a sin^{ n } λ equilibrium polar field profile with n > 7; the horizontal line marks the value for n = 8. 

In the text 
Fig. 5 Comparison of total net change in the solar axial dipole moment D_{n+1} − D_{n} in the 2 × 2D dynamo simulation and its approximations with the algebraic method in a segment of 80 cycles. Plotted are absolute values of the quantities shown in the legend. See text for further explanations. 

In the text 