Breakthroughs in our understanding of physical phenomena have traditionally followed improvements in instrumentation. Studies of the magnetic field of the Sun, and its influence on the solar dynamo and space weather events, have benefited from improvements in resolution and measurement frequency of new instruments. However, in order to fully understand the solar cycle, high-quality data across time-scales longer than the typical lifespan of a solar instrument are required. At the moment, discrepancies between measurement surveys prevent the combined use of all available data. In this work, we show that machine learning can help bridge the gap between measurement surveys by learning to \textbf{super-resolve} low-resolution magnetic field images and \textbf{translate} between characteristics of contemporary instruments in orbit. We also introduce the notion of physics-based metrics and losses for super-resolution to preserve underlying physics and constrain the solution space of possible super-resolution outputs.

Machine learning techniques have been successfully applied to super-resolution tasks on natural images where visually pleasing results are sufficient. However in many scientific domains this is not adequate and estimations of errors and uncertainties are crucial. To address this issue we propose a Bayesian framework that decomposes uncertainties into epistemic and aleatoric uncertainties. We test the validity of our approach by super-resolving images of the Sun's magnetic field and by generating maps measuring the range of possible high resolution explanations compatible with a given low resolution magnetogram.

In this paper we address a discrepancy between the surface flux evolution in a 3D kinematic dynamo model and a 2D surface flux transport model that has been closely calibrated to the real Sun. We demonstrate that the difference is due to the connectivity of active regions to the toroidal field at the base of the convection zone, which is not accounted for in the surface-only model. Initially, we consider the decay of a single active region, firstly in a simplified Cartesian 2D model and subsequently the full 3D model. By varying the turbulent diffusivity profile in the convection zone, we find that increasing the diffusivity - so that active regions are more rapidly disconnected from the base of the convection zone - improves the evolution of the surface field. However, if we simulate a full solar cycle, we find that the dynamo is unable to sustain itself under such an enhanced diffusivity. This suggests that in order to accurately model the solar cycle, we must find an alternative way to disconnect emerging active regions, whilst conserving magnetic flux.

In this paper we present a curated dataset from the NASA Solar Dynamics Observatory (SDO) mission in a format suitable for machine learning research. Beginning from level 1 scientific products we have processed various instrumental corrections, downsampled to manageable spatial and temporal resolutions, and synchronized observations spatially and temporally. We illustrate the use of this dataset with two example applications: forecasting future EVE irradiance from present EVE irradiance and translating HMI observations into AIA observations. For each application we provide metrics and baselines for future model comparison. We anticipate this curated dataset will facilitate machine learning research in heliophysics and the physical sciences generally, increasing the scientific return of the SDO mission. This work is a direct result of the 2018 NASA Frontier Development Laboratory Program. Please see the appendix for access to the dataset.

We test recent claims that the polar field at the end of Cycle 23 was weakened by a small number of large, abnormally oriented regions, and investigate what this means for solar cycle prediction. We isolate the contribution of individual regions from magnetograms for Cycles 21, 22 and 23 using a 2D surface flux transport model, and find that although the top ~10% of contributors tend to define sudden large variations in the axial dipole moment, the cumulative contribution of many weaker regions cannot be ignored. In order to recreate the axial dipole moment to a reasonable degree, many more regions are required in Cycle 23 than in Cycles 21 and 22 when ordered by contribution. We suggest that the negative contribution of the most significant regions of Cycle 23 could indeed be a cause of the weak polar field at the following cycle minimum and the low-amplitude Cycle 24. We also examine the relationship between a region's axial dipole moment contribution and its emergence latitude, flux, and initial axial dipole moment. We find that once the initial dipole moment of a given region has been measured, we can predict the long-term dipole moment contribution using emergence latitude alone.

Accurate prediction of solar activity calls for precise calibration of solar cycle models. Consequently we aim to find optimal parameters for models which describe the physical processes on the solar surface, which in turn act as proxies for what occurs in the interior and provide source terms for coronal models. We use a genetic algorithm to optimize surface flux transport models using National Solar Observatory (NSO) magnetogram data for Solar Cycle 23. This is applied to both a 1D model that inserts new magnetic flux in the form of idealized bipolar magnetic regions, and also to a 2D model that assimilates specific shapes of real active regions. The genetic algorithm searches for parameter sets (meridional flow speed and profile, supergranular diffusivity, initial magnetic field, and radial decay time) that produce the best fit between observed and simulated butterfly diagrams, weighted by a latitude-dependent error structure which reflects uncertainty in observations. Due to the easily adaptable nature of the 2D model, the optimization process is repeated for Cycles 21, 22, and 24 in order to analyse cycle-to-cycle variation of the optimal solution. We find that the ranges and optimal solutions for the various regimes are in reasonable agreement with results from the literature, both theoretical and observational. The optimal meridional flow profiles for each regime are almost entirely within observational bounds determined by magnetic feature tracking, with the 2D model being able to accommodate the mean observed profile more successfully. Differences between models appear to be important in deciding values for the diffusive and decay terms. In like fashion, differences in the behaviours of different solar cycles lead to contrasts in parameters defining the meridional flow and initial field strength.

In this work we introduce a new way of binning sunspot group data with the purpose of better understanding the impact of the solar cycle on sunspot properties and how this defined the characteristics of the extended minimum of cycle 23. Our approach assumes that the statistical properties of sunspots are completely determined by the strength of the underlying large-scale field and have no additional time dependencies. We use the amplitude of the cycle at any given moment (something we refer to as activity level) as a proxy for the strength of this deep-seated magnetic field. We find that the sunspot size distribution is composed of two populations: one population of groups and active regions and a second population of pores and ephemeral regions. When fits are performed at periods of different activity level, only the statistical properties of the former population, the active regions, is found to vary. Finally, we study the relative contribution of each component (small-scale versus large-scale) to solar magnetism. We find that when hemispheres are treated separately, almost every one of the past 12 solar minima reaches a point where the main contribution to magnetism comes from the small-scale component. However, due to asymmetries in cycle phase, this state is very rarely reached by both hemispheres at the same time. From this we infer that even though each hemisphere did reach the magnetic baseline, from a heliospheric point of view the minimum of cycle 23 was not as deep as it could have been.

In this work we take advantage of eleven different sunspot group, sunspot, and active region databases to characterize the area and flux distributions of photospheric magnetic structures. We find that, when taken separately, different databases are better fitted by different distributions (as has been reported previously in the literature). However, we find that all our databases can be reconciled by the simple application of a proportionality constant, and that, in reality, different databases are sampling different parts of a composite distribution. This composite distribution is made up by linear combination of Weibull and log-normal distributions -- where a pure Weibull (log-normal) characterizes the distribution of structures with fluxes below (above) $10^{21}$Mx ($10^{22}$Mx). We propose that this is evidence of two separate mechanisms giving rise to visible structures on the photosphere: one directly connected to the global component of the dynamo (and the generation of bipolar active regions), and the other with the small-scale component of the dynamo (and the fragmentation of magnetic structures due to their interaction with turbulent convection). Additionally, we demonstrate that the Weibull distribution shows the expected linear behavior of a power-law distribution (when extended into smaller fluxes), making our results compatible with the results of Parnell et al. (2009).

The solar cycle and its associated magnetic activity are the main drivers behind changes in the interplanetary environment and Earth's upper atmosphere (commonly referred to as space weather and climate). In recent years there has been an effort to develop accurate solar cycle predictions, leading to nearly a hundred widely spread predictions for the amplitude of solar cycle 24. Here we show that cycle predictions can be made more accurate if performed separately for each hemisphere, taking advantage of information about both the dipolar and quadrupolar moments of the solar magnetic field during minimum.

We propose a phenomenological technique for modelling the emergence of active regions within a three-dimensional, kinematic dynamo framework. By imposing localised velocity perturbations, we create emergent flux-tubes out of toroidal magnetic field at the base of the convection zone, leading to the eruption of active regions at the solar surface. The velocity perturbations are calibrated to reproduce observed active region properties (including the size and flux of active regions, and the distribution of tilt angle with latitude), resulting in a more consistent treatment of flux-tube emergence in kinematic dynamo models than artificial flux deposition. We demonstrate how this technique can be used to assimilate observations and drive a kinematic 3D model, and use it to study the characteristics of active region emergence and decay as a source of poloidal field. We find that the poloidal components are strongest not at the solar surface, but in the middle convection zone, in contrast with the common assumption that the poloidal source is located near the solar surface. We also find that, while most of the energy is contained in the lower convection zone, there is a good correlation between the evolution of the surface and interior magnetic fields.

The solar cycle and its associated magnetic activity are the main drivers behind changes in the interplanetary environment and the Earth's upper atmosphere (commonly referred to as space weather). These changes have a direct impact on the lifetime of space-based assets and can create hazards to astronauts in space. In recent years there has been an effort to develop accurate solar cycle predictions (with aims at predicting the long-term evolution of space weather), leading to nearly a hundred widely spread predictions for the amplitude of solar cycle 24. A major contributor to the disagreement is the lack of direct long-term databases covering different components of the solar magnetic field (toroidal vs.\ poloidal). Here we use sunspot area and polar faculae measurements spanning a full century (as our toroidal and poloidal field proxies), to study solar cycle propagation, memory, and prediction. Our results substantiate predictions based on the polar magnetic fields, whereas we find sunspot area to be uncorrelated to cycle amplitude unless multiplied by area-weighted average tilt. This suggests that the joint assimilation of tilt and sunspot area is a better choice (with aims to cycle prediction) than sunspot area alone, and adds to the evidence in favor of active region emergence and decay as the main mechanism of poloidal field generation (i.e. the Babcock-Leighton mechanism). Finally, by looking at the correlation between our poloidal and toroidal proxies across multiple cycles, we find solar cycle memory to be limited to only one cycle.

Direct observations over the past four centuries show that the number of sunspots observed on the Sun's surface vary periodically, going through successive maxima and minima. Following sunspot cycle 23, the Sun went into a prolonged minimum characterized by a very weak polar magnetic field and an unusually large number of days without sunspots. Sunspots are strongly magnetized regions and are generated by a dynamo mechanism which recreates the solar polar field mediated via plasma flows. Here we report results from kinematic dynamo simulations which demonstrate that a fast meridional flow in the early half of a cycle, followed by a slower flow in the latter half, reproduces both the characteristics of the minimum of sunspot cycle 23 - a large number of spotless days and a relatively weak polar field. Our model predicts that, in general, very deep minima are associated with weak polar fields. Sunspots govern the solar radiative energy and radio flux, and in conjunction with the polar field, modulate the solar wind, heliospheric open flux and consequently cosmic ray flux at Earth; our results therefore provide a causal link between solar internal dynamics and the atypical values of these heliospheric forcing parameters during the recently concluded solar minimum.

Although the Sun's polar magnetic fields are thought to provide important clues for understanding the 11-year sunspot cycle, including the observed variations of its amplitude and period, the current database of high-quality polar-field measurements spans relatively few sunspot cycles. In this paper we address this deficiency by consolidating Mount Wilson Observatory polar faculae data from four data reduction campaigns, validating it through a comparison with facular data counted automatically from MDI intensitygrams, and calibrating it against polar field measurements taken by the Wilcox Solar Observatory and average polar field and total polar flux calculated using MDI line-of-sight magnetograms. Our results show that the consolidated polar facular measurements are in excellent agreement with both polar field and polar flux estimates, making them an ideal proxy to study the evolution of the polar magnetic field. Additionally, we combine this database with sunspot area measurements to study the role of the polar magnetic flux in the evolution of the heliospheric magnetic field (HMF). We find that there is a strong correlation between HMF and polar flux at solar minimum and that, taken together, polar flux and sunspot area are better at explaining the evolution of the HMF during the last century than sunspot area alone.

The emergence of tilted bipolar active regions and the dispersal of their flux, mediated via processes such as diffusion, differential rotation and meridional circulation is believed to be responsible for the reversal of the Sun's polar field. This process (commonly known as the Babcock-Leighton mechanism) is usually modeled as a near-surface, spatially distributed $\alpha$-effect in kinematic mean-field dynamo models. However, this formulation leads to a relationship between polar field strength and meridional flow speed which is opposite to that suggested by physical insight and predicted by surface flux-transport simulations. With this in mind, we present an improved double-ring algorithm for modeling the Babcock-Leighton mechanism based on active region eruption, within the framework of an axisymmetric dynamo model. Using surface flux-transport simulations we first show that an axisymmetric formulation -- which is usually invoked in kinematic dynamo models -- can reasonably approximate the surface flux dynamics. Finally, we demonstrate that our treatment of the Babcock-Leighton mechanism through double-ring eruption leads to an inverse relationship between polar field strength and meridional flow speed as expected, reconciling the discrepancy between surface flux-transport simulations and kinematic dynamo models.

The turbulent magnetic diffusivity in the solar convection zone is one of the most poorly constrained ingredients of mean-field dynamo models. This lack of constraint has previously led to controversy regarding the most appropriate set of parameters, as different assumptions on the value of turbulent diffusivity lead to radically different solar cycle predictions. Typically, the dynamo community uses double step diffusivity profiles characterized by low values of diffusivity in the bulk of the convection zone. However, these low diffusivity values are not consistent with theoretical estimates based on mixing-length theory -- which suggest much higher values for turbulent diffusivity. To make matters worse, kinematic dynamo simulations cannot yield sustainable magnetic cycles using these theoretical estimates. In this work we show that magnetic cycles become viable if we combine the theoretically estimated diffusivity profile with magnetic quenching of the diffusivity. Furthermore, we find that the main features of this solution can be reproduced by a dynamo simulation using a prescribed (kinematic) diffusivity profile that is based on the spatiotemporal geometric-average of the dynamically quenched diffusivity. Here, we provide an analytic fit to the dynamically quenched diffusivity profile, which can be used in kinematic dynamo simulations. Having successfully reconciled the mixing-length theory estimated diffusivity profile with kinematic dynamo models, we argue that they remain a viable tool for understanding the solar magnetic cycle.

An essential ingredient in kinematic dynamo models is the velocity field within the solar convection zone. In particular, the differential rotation is now well constrained by helioseismic observations. Helioseismology also gives us information about the depth-dependence of the meridional circulation in the near-surface layers. The typical velocity inputs used in solar dynamo models, however, continue to be an analytic fit to the observed differential rotation and a theoretically constructed meridional flow profile that matches only the peak flow speed at the surface. Here we take the first steps towards realistic helioseismic data assimilation, by presenting methodologies for constructing differential rotation and meridional circulation profiles that more closely conform to the observational constraints currently available. We also present simulations driven by the assimilated rotation and four plausible profiles for the internal meridional circulation -- all of which match the helioseismically inferred near-surface depth-dependence, but whose magnitudes are made to vary. We discuss how the results compare with those that are driven by purely analytic fits. Our results indicate that the latitudinal shear of the rotation in the bulk of the solar convection zone plays a more important role, than either the tachocline or surface radial shear, in the induction of toroidal field. We also find that it is the speed of the equatorward counter-flow in the meridional flow at the base of the convection zone, and not how far into the radiative interior it penetrates, that primarily determines the dynamo cycle period. Given that improved helioseismic constraints are expected to be available in the future, our analysis lays the basis for assimilating these data within dynamo models.