Bipolar Magnetic Regions in Solar Dynamics

• Bipolar Magnetic Regions are defined by opposite magnetic polarities with systematic tilts and separations, serving as the fundamental units of the Sun’s dynamo-generated magnetic field.
• Advanced detection and tracking methods like BARD and AutoTAB enable precise characterization of BMR emergence, evolution, and their contribution to surface flux transport.
• BMR observations underpin solar dynamo models by revealing nonlinear feedbacks, tilt quenching, and the role of mixed polarity regions in shaping the solar cycle.

Searching arXiv for recent and foundational papers on Bipolar Magnetic Regions to ground the article in the literature. Bipolar Magnetic Regions (BMRs) are pairs of opposite-polarity magnetic concentrations on the solar photosphere, observed in magnetograms as the magnetic counterparts of sunspot pairs and active regions, but extending to weaker spot-free bipoles as well. They are the basic surface manifestations of the Sun’s dynamo-generated toroidal flux, the principal discrete source of poloidal field in Babcock–Leighton theory, and a central organizing unit for surface flux transport, polar-field reversal, and the spatiotemporal structure of the solar cycle (Stenflo et al., 2011, Miesch et al., 2014, Munoz-Jaramillo et al., 2022).

1. Observational definition and physical character

A BMR consists of two main opposite-polarity flux concentrations, usually described as leading and following polarities, with a characteristic east–west ordering, a nonzero polarity separation, and a systematic tilt relative to the equator. In magnetogram-based studies, BMRs span a broad hierarchy of sizes: the SOHO/MDI global survey identified 160,079 BMRs from 1996–2011 spanning nearly four orders of magnitude, with a lower flux limit of about $10^{20}\,\mathrm{Mx}$, just above the size range of ephemeral active regions (Stenflo et al., 2011). In active-region-oriented representations, total unsigned fluxes extend to $10^{23}\,\mathrm{Mx}$, emphasizing that quoted flux ranges depend on the selection function and on whether the sample includes only sunspot-bearing regions or also weaker bipoles (Yeates, 2020).

The observed population is not morphologically uniform. Analyses of MDI and HMI magnetograms found a bimodal distribution of the maximum magnetic field $B_{\rm max}$, with peaks near $600\,\mathrm{G}$ and $2100\,\mathrm{G}$; the low-field peak is dominated by BMRs without white-light sunspot counterparts, whereas the high-field peak is associated with sunspots and pores (Jha et al., 2019). This suggests that “BMR” is a magnetic rather than photometric category: sunspot groups are the strong, high-contrast subset of a broader family of bipolar emergences. In that broader sense, BMRs are the cornerstone of solar variability, tracers of the large-scale magnetic processes that give rise to the solar cycle, shapers of the solar corona, and building blocks of the large-scale solar magnetic field (Munoz-Jaramillo et al., 2022).

2. Detection, tracking, and catalog construction

Systematic BMR studies depend critically on how bipoles are detected and associated across time. In the BARD framework, positive and negative polarities are detected separately from line-of-sight magnetograms, same-polarity patches are merged, and opposite polarities are paired by minimizing a cost function

$$F(d,\Delta A,\Delta \Phi)=d^4\,\Delta A\,\Delta \Phi,$$

which favors close, balanced, similarly sized pairs. Limited human supervision then corrects mispairings, fragmentation, and tracking errors; about $10\%$ of cataloged objects required human adjustment, and the catalog contains more than 10,000 unique BMRs tracked and characterized during every day of their observation across four solar cycles (Munoz-Jaramillo et al., 2022).

AutoTAB shifted the emphasis from snapshot statistics to lifetime-resolved tracking. It starts from pre-identified BMR masks in MDI and HMI line-of-sight magnetograms, refines the masks with a $100\,\mathrm{G}$ threshold and morphological closing, and associates regions across frames by differentially rotating the initial binary mask and testing pixel overlap. In the 2023 description, AutoTAB successfully tracked 9152 BMRs over 1996–2020 and was reported to perform even for small BMRs with flux $\sim 10^{20}\,\mathrm{Mx}$; the 2024 catalog extended this to 9232 tracked BMRs over 1996–2022 (Sreedevi et al., 2023, Sreedevi et al., 2024). This design addresses a major bias of earlier catalogs, namely that long-lived BMRs were counted repeatedly at different evolutionary phases.

Backtracking to true emergence revealed that not every automatically tracked bipole is a fresh emergence. An additional AutoTAB module backtracked about 12000 identified BMRs and successfully recovered 3080 emergence histories. Within that sample, one group showed the expected emergence signature of substantial flux growth, whereas another displayed little or no growth. The latter, classified as non-emerging BMRs, account for about $27\%$ of the backtracked sample and do not show a preferred tilt-angle distribution or systematic latitudinal tilt dependence. Their inclusion can therefore distort statistical inferences about Joy’s law and related emergence physics (Sreedevi et al., 21 Sep 2025).

3. Tilt laws, polarity rules, and emergence-phase diagnostics

Two empirical regularities organize BMR orientation statistics. Hale’s polarity law specifies the leading-polarity sign in each hemisphere and its reversal from cycle to cycle. Joy’s law specifies a systematic latitudinal dependence of the tilt. In the full SOHO/MDI survey, the mean tilt followed

$10^{23}\,\mathrm{Mx}$0

with no indication of a dependence on region size, flux, or bipolar moment within the accessible range (Stenflo et al., 2011). The same survey also showed that a few percent of all regions violate Hale’s law, and that well-defined medium-size anti-Hale regions can occur side by side with Hale-oriented regions in the same latitude zone, demonstrating that the observed population is not perfectly represented by a single coherent toroidal flux system (Stenflo et al., 2011).

Later MDI/HMI analyses introduced a more field-dependent picture. Using $10^{23}\,\mathrm{Mx}$1 as the organizing variable, the Joy’s-law slope $10^{23}\,\mathrm{Mx}$2 was found to increase slowly with $10^{23}\,\mathrm{Mx}$3 at low field strengths but to decrease once $10^{23}\,\mathrm{Mx}$4; the tilt scatter about Joy’s law decreases systematically with increasing $10^{23}\,\mathrm{Mx}$5 (Jha et al., 2019). This behavior was summarized as a nonlinear tilt-quenching law of the form

$10^{23}\,\mathrm{Mx}$6

with best-fit parameters $10^{23}\,\mathrm{Mx}$7 and $10^{23}\,\mathrm{Mx}$8 (Jha et al., 2019). AutoTAB-based tracking confirmed that the bimodal $10^{23}\,\mathrm{Mx}$9 distribution is not an artifact of multiple counting and found the downturn in $B_{\rm max}$0 for $B_{\rm max}$1, again suggesting tilt quenching in the strong-field regime (Jha et al., 2024).

Emergence-phase analyses sharpened the physical interpretation. In tracked AutoTAB BMRs, polarity separation increases over normalized lifetime, and early-time samples already display Joy’s-law behavior. A more targeted emergence study found that backtracked, genuinely emerging BMRs obey Joy’s law already at appearance, with

$B_{\rm max}$2

while the tilt scatter decreases substantially from emergence to the matured phase (Sreedevi et al., 5 Nov 2025). A related AutoTAB analysis reported that BMR polarity separation increases over normalized lifetime, that tilt increases with flux over much of the observed range, and that lower-flux regions show larger tilt fluctuations than stronger regions (Sreedevi et al., 2024). Taken together, these results support the view that at least part of the systematic tilt is acquired beneath the photosphere, while the large emergence-phase scatter reflects strong interaction with turbulent convection. They also clarify why mixed samples containing non-emerging bipoles can obscure the underlying Joy’s-law signal (Sreedevi et al., 21 Sep 2025).

4. Surface flux transport and the bipolar approximation

In surface flux transport (SFT) theory, BMRs are the source term that determines the buildup of the Sun’s axial dipole and polar field. The key large-scale diagnostic is the axial dipole moment

$B_{\rm max}$3

and for fitted BMRs the initial dipole scales approximately as $B_{\rm max}$4, where $B_{\rm max}$5 is flux, $B_{\rm max}$6 the polarity separation, $B_{\rm max}$7 the tilt, and $B_{\rm max}$8 the emergence latitude (Yeates, 2020). This approximation is foundational but imperfect. Using 1090 BMRs fitted to SHARP active-region patches from 2010–2020, the bipolar approximation was shown to match the initial flux and axial dipole well, yet it overestimates the end-of-cycle net axial dipole moment by $B_{\rm max}$9: the SHARP-driven simulation gave $600\,\mathrm{G}$0, whereas the BMR-driven simulation gave $600\,\mathrm{G}$1 with the same transport parameters (Yeates, 2020). The stated cause is the neglect of multipolar structure and polarity asymmetry in real active regions. A plausible implication is that symmetric single-bipole source terms systematically overstate polar-field production unless complexity-dependent corrections are introduced.

BMR-based SFT models also require a relation between sunspot observables and the total magnetic content of the newly emerged bipole. For sunspot-bearing BMRs, the total magnetic flux inferred from disc-integrated facular and network flux is

$600\,\mathrm{G}$2

where $600\,\mathrm{G}$3 is the total sunspot area in $600\,\mathrm{G}$4Hem (Yeo et al., 2021). Implementing this relation in an established SFTM preserved its ability to reproduce open flux, polar-field reversals, and the correlation between the polar field at cycle end and the amplitude of the following cycle (Yeo et al., 2021). This places BMRs at the interface between magnetogram-based source assimilation and long-baseline sunspot reconstructions.

5. BMRs in dynamo models and nonlinear feedbacks

In three-dimensional Babcock–Leighton dynamo models, BMRs are not merely empirical input but explicit dynamical agents. The 3D kinematic model of Miesch and Dikpati inserts tilted sunspot pairs, more generally BMRs, on the surface in response to the dynamo-generated toroidal field. The surface radial field of a BMR is represented as

$600\,\mathrm{G}$5

with the tilt imposed by a Joy’s-law relation

$600\,\mathrm{G}$6

thereby unifying the axisymmetric Babcock–Leighton source with the non-axisymmetric surface-flux-transport evolution of individual bipoles (Miesch et al., 2014).

Subsequent models turned BMR statistics into explicit saturation physics. In a 3D Babcock–Leighton dynamo, the eruption threshold

$600\,\mathrm{G}$7

forces stronger cycles to produce BMRs at higher mean latitudes; because high-latitude BMRs are less efficient at generating axial dipole moment, the latitudinal variation of BMR emergence acts as a nonlinear saturation mechanism (Karak, 2020). In a different 3D setting, converging inflows toward BMRs with observed amplitudes of $600\,\mathrm{G}$8–$600\,\mathrm{G}$9 were modeled as magnetic-field-dependent flows proportional to the gradients of the smoothed unsigned radial field. These inflows increase local flux cancellation within BMRs, reduce the buildup of the global poloidal field, and by themselves can saturate a fully 3D solar dynamo (Teweldebirhan et al., 2023).

Stochastic BMR properties are equally important. Using the STABLE framework in both SFT and dynamo modes, one study found that randomness in emergence delay or latitude produced negligible or modest cycle variability under fixed distributions, whereas randomness in BMR flux produced substantial effects and scatter in the tilt around Joy’s law produced the largest variation in the polar field and the solar cycle (Kumar et al., 2024). Mean-field models that include BMR-related buoyancy and tilt terms similarly conclude that BMRs substantially affect cycle strength, but that the observed weak Cycle 24 still required a decrease of turbulent helicity in the bulk of the convection zone, rather than BMR effects alone (Pipin et al., 2022). More recent non-linear mean-field work connects BMR tilt and twist to magnetic helicity transport, finding that BMR-related tilt/twist dominates the helicity flux at the beginning of a BMR’s evolution, whereas differential rotation dominates at later stages (Pipin et al., 1 Aug 2025).

6. Open problems, competing interpretations, and current directions

A central debate concerns the physical origin of BMR tilt. One line of work, based on the full MDI survey, emphasized the absence of size dependence in the mean Joy’s-law amplitude and the coexistence of Hale and anti-Hale medium-size regions in the same latitude zone; it argued that such observations are incompatible with a paradigm based on coherent subsurface toroidal flux ropes as the sole source of sunspots and that there is no observational support for a separation between a global and a local dynamo (Stenflo et al., 2011). Another line, based on tracked and backtracked AutoTAB emergence samples, finds Joy’s law already at appearance, increasing polarity separation over normalized lifetime, and larger tilt fluctuations in weaker BMRs, all of which support the thin-flux-tube picture in which the Coriolis force and helical convection act on rising flux tubes beneath the surface (Sreedevi et al., 2024, Sreedevi et al., 5 Nov 2025). This suggests that the observational tension may depend strongly on which population is measured, which evolutionary phase is isolated, and whether non-emerging remnants are removed from the sample.

A second unresolved issue is representational fidelity. Symmetric single-bipole models are indispensable for large-scale transport and dynamo calculations, yet they miss complexity, asymmetry, and mixed populations. The SHARP comparison showed a $2100\,\mathrm{G}$0 overestimate of the end-of-cycle dipole when real active regions are reduced to idealized BMRs (Yeates, 2020). Emergence backtracking showed that non-emerging bipoles lack Joy’s law and can mask the physical properties of genuine emergers if they remain in the sample (Sreedevi et al., 21 Sep 2025). In parallel, long-baseline catalog work argues for homogeneous, multi-scale magnetic databases and explicitly envisions combining BARD with SWAMIS to build an extended multi-scale magnetic catalog (Munoz-Jaramillo et al., 2022). A plausible implication is that future progress will depend less on ever larger raw BMR counts than on cleaner emergence classification, cross-instrument homogeneity, and systematic treatment of multipolar active-region structure.