Measurements of the large cellular flows on the Sun were made by local correlation tracking of supergranules seen in full-disk Doppler images obtained by the HMI instrument on the NASA SDO satellite. The hourly measurements were averaged over 34 days to produce daily maps of the latitudinal and longitudinal velocities. While flows at all latitudes are largely in the form of vortices with left-handed helicity in the north and right-handed helicity in the south, there are key distinctions between the low latitude and high latitudes cells. The low latitude cells have roughly circular shapes, lifetimes of about one month, rotate nearly rigidly, do not drift in latitude, and do not exhibit any correlation between longitudinal and latitudinal flow. The high latitude cells have long extensions that spiral inward toward the poles. They have lifetimes of several months, rotate differentially with latitude, drift poleward at speeds approaching 2 m s$^{-1}$, and have a strong correlation between prograde and equatorward flows. Spherical harmonic spectral analyses confirm that the flows are dominated by the curl component with RMS velocities of about 12 m s$^{-1}$ at wavenumber $\ell$ = 10. Fourier transforms in time indicate two notable components - an $m = \pm\ell$ feature representing the low latitude component and an $m = \pm1$ feature representing the high latitude component. The dispersion relation for the low latitude component is well represented by that derived for Rossby waves. The high latitude component has a constant temporal frequency for all $\ell$ indicating features advected by differential rotation at rates representative of the base of the convection zone high latitudes. The poleward motions of these features further suggest that the high latitude meridional flow at the base of the convection zone is poleward - not equatorward.

Over the last decade there has been mounting evidence that the strength of the Sun's polar magnetic fields during a solar cycle minimum is the best predictor of the amplitude of the next solar cycle. Surface flux transport models can be used to extend these predictions by evolving the Sun's surface magnetic field to obtain an earlier prediction for the strength of the polar fields, and thus the amplitude of the next cycle. In 2016, our Advective Flux Transport (AFT) model was used to do this, producing an early prediction for Solar Cycle 25. At that time, AFT predicted that Cycle 25 will be similar in strength to the Cycle 24, with an uncertainty of about 15% . AFT also predicted that the polar fields in the southern hemisphere would weaken in late 2016 and into 2017 before recovering. That AFT prediction was based on the magnetic field configuration at the end of January 2016. We now have 2 more years of observations. We examine the accuracy of the 2016 AFT prediction and find that the new observations track well with AFT's predictions for the last two years. We show that the southern relapse did in fact occur, though the timing was off by several months. We propose a possible cause for the southern relapse and discuss the reason for the offset in timing. Finally, we provide an updated AFT prediction for Solar Cycle 25 which includes solar observations through January of 2018.

Evidence strongly indicates that the strength of the Sun's polar fields near the time of a sunspot cycle minimum determines the strength of the following solar activity cycle. We use our Advective Flux Transport (AFT) code, with flows well constrained by observations, to simulate the evolution of the Sun's polar magnetic fields from early 2016 to the end of 2019 --- near the expected time of Cycle 24/25 minimum. We run a series of simulations in which the uncertain conditions (convective motion details, active region tilt, and meridional flow profile) are varied within expected ranges. We find that the average strength of the polar fields near the end of Cycle 24 will be similar to that measured near the end of Cycle 23, indicating that Cycle 25 will be similar in strength to the current cycle. In all cases the polar fields are asymmetric with fields in the south stronger than those in the north. This asymmetry would be more pronounced if not for the predicted weakening of the southern polar fields in late 2016 and through 2017. After just four years of simulation the variability across our ensemble indicates an accumulated uncertainty of about 15\%. This accumulated uncertainty arises from stochastic variations in the convective motion details, the active region tilt, and changes in the meridional flow profile. These variations limit the ultimate predictability of the solar cycle.

With multiple vantage points around the Sun, STEREO and SDO imaging observations provide a unique opportunity to view the solar surface continuously. We use He II 304 A data from these observatories to isolate and track ten active regions and study their long-term evolution. We find that active regions typically follow a standard pattern of emergence over several days followed by a slower decay that is proportional in time to the peak intensity in the region. Since STEREO does not make direct observations of the magnetic field, we employ a flux-luminosity relationship to infer the total unsigned magnetic flux evolution. To investigate this magnetic flux decay over several rotations we use a surface flux transport model, the Advective Flux Transport (AFT) model, that simulates convective flows using a time-varying velocity field and find that the model provides realistic predictions when information about the active region's magnetic field strength and distribution at peak flux is available. Finally, we illustrate how 304 \AA\ images can be used as a proxy for magnetic flux measurements when magnetic field data is not accessible.

Spectra of the cellular photospheric flows are determined from full-disk Doppler velocity observations acquired by the Helioseismic and Magnetic Imager (HMI) instrument on the Solar Dynamics Observatory (SDO) spacecraft. Three different analysis methods are used to separately determine spectral coefficients representing the poloidal flows, the toroidal flows, and the radial flows. The amplitudes of these spectral coefficients are constrained by simulated data analyzed with the same procedures as the HMI data. We find that the total velocity spectrum rises smoothly to a peak at a wavenumber of about 120 (wavelength of about 35 Mm), which is typical of supergranules. The spectrum levels off out to wavenumbers of about 400, and then rises again to a peak at a wavenumber of about 3500 (wavelength of about 1200 km), which is typical of granules. The velocity spectrum is dominated by the poloidal flow component (horizontal flows with divergence but no curl) at wavenumbers above 30. The toroidal flow component (horizontal flows with curl but no divergence) dominates at wavenumbers less than 30. The radial flow velocity is only about 3\% of the total flow velocity at the lowest wavenumbers, but increases in strength to become about 50\% at wavenumbers near 4000. The spectrum compares well with the spectrum of giant cell flows at the lowest wavenumbers and with the spectrum of granulation from a 3D radiative-hydrodynamic simulation at the highest wavenumbers.

The Solar Cycle is reviewed. The 11-year cycle of solar activity is characterized by the rise and fall in the numbers and surface area of sunspots. A number of other solar activity indicators also vary in association with the sunspots including; the 10.7cm radio flux, the total solar irradiance, the magnetic field, flares and coronal mass ejections, geomagnetic activity, galactic cosmic ray fluxes, and radioisotopes in tree rings and ice cores. Individual solar cycles are characterized by their maxima and minima, cycle periods and amplitudes, cycle shape, the equatorward drift of the active latitudes, hemispheric asymmetries, and active longitudes. Cycle-to-cycle variability includes the Maunder Minimum, the Gleissberg Cycle, and the Gnevyshev-Ohl (even-odd) Rule. Short-term variability includes the 154-day periodicity, quasi-biennial variations, and double-peaked maxima. We conclude with an examination of prediction techniques for the solar cycle and a closer look at cycles 23 and 24.

After emerging to the solar surface, the Sun's magnetic field displays a complex and intricate evolution. The evolution of the surface field is important for several reasons. One is that the surface field, and its dynamics, sets the boundary condition for the coronal and heliospheric magnetic fields. Another is that the surface evolution gives us insight into the dynamo process. In particular, it plays an essential role in the Babcock-Leighton model of the solar dynamo. Describing this evolution is the aim of the surface flux transport model. The model starts from the emergence of magnetic bipoles. Thereafter, the model is based on the induction equation and the fact that after emergence the magnetic field is observed to evolve as if it were purely radial. The induction equation then describes how the surface flows -- differential rotation, meridional circulation, granular, supergranular flows, and active region inflows -- determine the evolution of the field (now taken to be purely radial). In this paper, we review the modeling of the various processes that determine the evolution of the surface field. We restrict our attention to their role in the surface flux transport model. We also discuss the success of the model and some of the results that have been obtained using this model.

The faster meridional flow that preceded the solar cycle 23/24 minimum is thought to have led to weaker polar field strengths, producing the extended solar minimum and the unusually weak cycle 24. To determine the impact of meridional flow variations on the sunspot cycle, we have simulated the Sun's surface magnetic field evolution with our newly developed surface flux transport model. We investigate three different cases: a constant average meridional flow, the observed time-varying meridional flow, and a time-varying meridional flow in which the observed variations from the average have been doubled. Comparison of these simulations shows that the variations in the meridional flow over cycle 23 have a significant impact (~20%) on the polar fields. However, the variations produced polar fields that were stronger than they would have been otherwise. We propose that the primary cause of the extended cycle 23/24 minimum and weak cycle 24 was the weakness of cycle 23 itself - with fewer sunspots, there was insufficient flux to build a big cycle. We also find that any polar counter-cells in the meridional flow (equatorward flow at high latitudes) produce flux concentrations at mid-to-high latitudes that are not consistent with observations.

We have measured the meridional motions of the magnetic elements in the Sun's surface layers since 1996 and find systematic and substantial variations. In general the meridional flow speed is fast at cycle minima and slow at cycle maxima. We find that these systematic variations are characterized by a weakening of the meridional flow on the poleward sides of the active (sunspot) latitudes. This can be interpreted as a inflow toward the sunspot zones superimposed on a more general poleward meridional flow profile. We also find variations in the meridional flow which vary from cycle-to-cycle. The meridional flow was slower at both the minimum and maximum of cycle 23 compared to similar phases of cycles 21, 22, and 24. Models of the magnetic flux transport by a variable meridional flow suggest that it can significantly modulate the size and timing of the following sunspot cycle through its impact on the Sun's polar magnetic fields. We suggest that the meridional flow variations observed in cycle 23 contributed to the weak polar fields at the end of the cycle which then produced a weak cycle 24 and the extraordinary cycle 23/24 minimum.

Heat is transported through the outermost 30% of the Sun's interior by overturning convective motions. These motions are evident at the Sun's surface in the form of two characteristic cellular structures - granules and supergranules (~1000 and ~30,000 km across respectively). The existence of much larger cells has been suggested by both theory and observation for over 45 years. We found evidence for giant cellular flows that persist for months by tracking the motions of supergranules. As expected from the effects of the Sun's rotation, the flows in these cells are clockwise around high pressure in the north, counter-clockwise in the south and transport angular momentum toward the equator, maintaining the Sun's rapid equatorial rotation.

The Sun's polar magnetic fields are directly related to solar cycle variability. The strength of the polar fields at the start (minimum) of a cycle determine the subsequent amplitude of that cycle. In addition, the polar field reversals at cycle maximum alter the propagation of galactic cosmic rays throughout the heliosphere in fundamental ways. We describe a surface magnetic flux transport model that advects the magnetic flux emerging in active regions (sunspots) using detailed observations of the near-surface flows that transport the magnetic elements. These flows include the axisymmetric differential rotation and meridional flow and the non-axisymmetric cellular convective flows (supergranules) all of which vary in time in the model as indicated by direct observations. We use this model with data assimilated from full-disk magnetograms to produce full surface maps of the Sun's magnetic field at 15-minute intervals from 1996 May to 2013 July (all of sunspot cycle 23 and the rise to maximum of cycle 24). We tested the predictability of this model using these maps as initial conditions, but with daily sunspot area data used to give the sources of new magnetic flux. We find that the strength of the polar fields at cycle minimum and the polar field reversals at cycle maximum can be reliably predicted up to three years in advance. We include a prediction for the cycle 24 polar field reversal.

Daily records of sunspot group areas compiled by the Royal Observatory, Greenwich, from May of 1874 through 1976 indicate a curious history for the penumbral areas of the smaller sunspot groups. On average, the ratio of penumbral area to umbral area in a sunspot group increases from 5 to 6 as the total sunspot group area increases from 100 to 2000 microHem (a microHem is a millionth the area of a solar hemisphere). This relationship does not vary substantially with sunspot group latitude or with the phase of the sunspot cycle. However, for the sunspot groups with total areas <100 microHem, this ratio changes dramatically and systematically through this historical record. The ratio for these smallest sunspots is near 5.5 from 1874 to 1900. After a rapid rise to more than 7 in 1905 it drops smoothly to less than 3 by 1930 and then rises smoothly back to more than 7 in 1961. It then returns to near 5.5 from 1965 to 1976. The smooth variation from 1905 to 1961 shows no indication of any step-like changes that might be attributed to changes in equipment or personnel. The overall level of solar activity was increasing monotonically during this time period when the penumbra-to-umbra area ratio dropped to less than half its peak value and then returned. If this history can be confirmed by other observations (e.g., Mt. Wilson or Kodaikanal) it may impact our understanding of penumbra formation, our dynamo models, and our estimates of historical changes in the solar irradiance.

Among many other measurable quantities the summer of 2009 saw a considerable low in the radiative output of the Sun that was temporally coincident with the largest cosmic ray flux ever measured at 1AU. A hemispheric asymmetry in magnetic activity is clearly observed and its evolution monitored and the resulting (prolonged) magnetic imbalance must have had a considerable impact on the structure and energetics of the heliosphere. While we cannot uniquely tie the variance and scale of the surface magnetism to the dwindling radiative and particulate output of the star, or the increased cosmic ray flux through the 2009 minimum, the timing of the decline and rapid recovery in early 2010 would appear to inextricably link them. These observations support a picture where the Sun's hemispheres are significantly out of phase with each other. Studying historical sunspot records with this picture in mind shows that the northern hemisphere has been leading since the middle of the last century and that the hemispheric "dominance" has changed twice in the past 130 years. The observations presented give clear cause for concern, especially with respect to our present understanding of the processes that produce the surface magnetism in the (hidden) solar interior - hemispheric asymmetry is the normal state - the strong symmetry shown in 1996 was abnormal. Further, these observations show that the mechanism(s) which create and transport the magnetic flux magnetic flux are slowly changing with time and, it appears, with only loose coupling across the equator such that those asymmetries can persist for a considerable time. As the current asymmetry persists and the basal energetics of the system continue to dwindle we anticipate new radiative and particulate lows coupled with increased cosmic ray fluxes heading into the next solar minimum.

The meridional circulation at high latitudes is crucial to the build-up and reversal of the Sun's polar magnetic fields. Here we characterize the axisymmetric flows by applying a magnetic feature cross-correlation procedure to high resolution magnetograms obtained by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We focus on Carrington Rotations 2096-2107 (April 2010 to March 2011) - the overlap interval between HMI and the Michelson Doppler Investigation (MDI). HMI magnetograms averaged over 720 seconds are first mapped into heliographic coordinates. Strips from these maps are then cross-correlated to determine the distances in latitude and longitude that the magnetic element pattern has moved, thus providing meridional flow and differential rotation velocities for each rotation of the Sun. Flow velocities were averaged for the overlap interval and compared to results obtained from MDI data. This comparison indicates that these HMI images are rotated counter-clockwise by 0.075 degrees with respect to the Sun's rotation axis. The profiles indicate that HMI data can be used to reliably measure these axisymmetric flow velocities to at least within 5 degrees of the poles. Unlike the noisier MDI measurements, no evidence of a meridional flow counter-cell is seen in either hemisphere with the HMI measurements: poleward flow continues all the way to the poles. Slight North-South asymmetries are observed in the meridional flow. These asymmetries should contribute to the observed asymmetries in the polar fields and the timing of their reversals.

Recent analysis revealed that supergranules (convection cells seen at the Sun's surface) are advected by the zonal flows at depths equal to the widths of the cells themselves. Here we probe the structure of the meridional circulation by cross-correlating maps of the Doppler velocity signal using a series of successively longer time lags between maps. We find that the poleward meridional flow decreases in amplitude with time lag and reverses direction to become an equatorward return flow at time lags > 24 hours. These cross-correlation results are dominated by larger and deeper cells at longer time lags. (The smaller cells have shorter lifetimes and do not contribute to the correlated signal at longer time lags.) We determine the characteristic cell size associated with each time lag by comparing the equatorial zonal flows measured at different time lags with the zonal flows associated with different cell sizes from a Fourier analysis. This association gives a characteristic cell size of ~50 Mm at a 24 hour time lag. This indicates that the poleward meridional flow returns equatorward at depths > 50 Mm -- just below the base of the surface shear layer. A substantial and highly significant equatorward flow (4.6 +/- 0.4 m/s) is found at a time lag of 28 hours corresponding to a depth of ~70 Mm. This represents one of the first positive detections of the Sun's meridional return flow and illustrates the power of using supergranules to probe the Sun's internal dynamics.

Using magnetic and microwave butterfly diagrams, we compare the behavior of solar polar regions to show that (i) the polar magnetic field and the microwave brightness temperature during the solar minimum substantially diminished during the cycle 23/24 minimum compared to the 22/23 minimum. (ii) The polar microwave brightness temperature (b) seems to be a good proxy for the underlying magnetic field strength (B). The analysis indicates a relationship, B = 0.0067Tb - 70, where B is in G and Tb in K. (iii) Both the brightness temperature and the magnetic field strength show north-south asymmetry most of the time except for a short period during the maximum phase. (iv) The rush-to-the-pole phenomenon observed in the prominence eruption activity seems to be complete in the northern hemisphere as of March 2012. (v) The decline of the microwave brightness temperature in the north polar region to the quiet-Sun levels and the sustained prominence eruption activity poleward of 60oN suggest that solar maximum conditions have arrived at the northern hemisphere. The southern hemisphere continues to exhibit conditions corresponding to the rise phase of solar cycle 24.

Supergranules are convection cells seen at the Sun's surface as a space filling pattern of horizontal flows. While typical supergranules have diameters of about 35 Mm, they exhibit a broad spectrum of sizes from ~10 Mm to ~100 Mm. Here we show that supergranules of different sizes can be used to probe the rotation rate in the Sun's outer convection zone. We find that the equatorial rotation rate as a function of depth as measured by global helioseismology matches the equatorial rotation as a function of wavelength for the supergranules. This suggests that supergranules are advected by flows at depths equal to their wavelengths and thus can be used to probe flows at those depths. The supergranule rotation profiles show that the surface shear layer, through which the rotation rate increases inward, extends to depths of ~50 Mm and to latitudes of at least 70 degrees. Typical supergranules are well observed at high latitudes and have a range of sizes that extend to greater depths than those typically available for measuring subsurface flows with local helioseismology. These characteristics indicate that probing the solar convection zone dynamics with supergranules can complement the results of helioseismology.

The latitudinal location of the sunspot zones in each hemisphere is determined by calculating the centroid position of sunspot areas for each solar rotation from May 1874 to June 2011. When these centroid positions are plotted and analyzed as functions of time from each sunspot cycle maximum there appears to be systematic differences in the positions and equatorward drift rates as a function of sunspot cycle amplitude. If, instead, these centroid positions are plotted and analyzed as functions of time from each sunspot cycle minimum then most of the differences in the positions and equatorward drift rates disappear. The differences that remain disappear entirely if curve fitting is used to determine the starting times (which vary by as much as 8 months from the times of minima). The sunspot zone latitudes and equatorward drift measured relative to this starting time follow a standard path for all cycles with no dependence upon cycle strength or hemispheric dominance. Although Cycle 23 was peculiar in its length and the strength of the polar fields it produced, it too shows no significant variation from this standard. This standard law, and the lack of variation with sunspot cycle characteristics, is consistent with Dynamo Wave mechanisms but not consistent with current Flux Transport Dynamo models for the equatorward drift of the sunspot zones.

The Sun's global meridional circulation is evident as a slow poleward flow at its surface. This flow is observed to carry magnetic elements poleward - producing the Sun's polar magnetic fields as a key part of the 11-year sunspot cycle. Current theories for the sunspot cycle assume that this surface flow is part of a circulation which sinks inward at the poles and turns equatorward at depths below 100 Mm. Here we use the advection of the Sun's convection cells by the meridional flow to map the flow velocity in latitude and depth. Our measurements show the largest cells clearly moving equatorward at depths below 35 Mm - the base of the Sun's surface shear layer. This surprisingly shallow return flow indicates the need for substantial revisions to solar/stellar dynamo theory.

We measure the axisymmetric transport of magnetic flux on the Sun by cross-correlating narrow strips of data from line-of-sight magnetograms obtained at a 96-minute cadence by the MDI instrument on the ESA/NASA SOHO spacecraft and then averaging the flow measurements over each synodic rotation of the Sun. Our measurements indicate that the axisymmetric flows vary systematically over the solar cycle. The differential rotation is weaker at maximum than at minimum. The meridional flow is faster at minimum and slower at maximum. The meridional flow speed on the approach to the Cycle 23/24 minimum was substantially faster than it was at the Cycle 22/23 minimum. The average latitudinal profile is largely a simple sinusoid that extends to the poles and peaks at about $35\degr$ latitude. As the cycle progresses a pattern of in-flows toward the sunspot zones develops and moves equatorward in step with the sunspot zones. These in-flows are accompanied by the torsional oscillations. This association is consistent with the effects of the Coriolis force acting on the in-flows. The equatorward motions associated with these in-flows are identified as the source of the decrease in net poleward flow at cycle maxima. We also find polar counter-cells (equatorward flow at high latitudes) in the south from 1996 to 2000 and in the north from 2002 to 2010. We show that these measurements of the flows are not affected by the non-axisymmetric diffusive motions produced by supergranulation.

We show that the motions of supergranules are consistent with a model in which they are simply advected by the axisymmetric flows in the Sun's surface shear layer. We produce a 10-day series of simulated Doppler images at a 15-minute cadence that reproduces most spatial and temporal characteristics seen in the SOHO/MDI Doppler data. Our simulated data have a spectrum of cellular flows with just two components -- a granule component that peaks at spherical wavenumbers of about 4000 and a supergranule component that peaks at wavenumbers of about 110. We include the advection of these cellular components by the axisymmetric flows -- differential rotation and meridional flow -- whose variations with latitude and depth (wavenumber) are consistent with observations. We mimic the evolution of the cellular pattern by introducing random variations to the phases of the spectral components at rates that reproduce the levels of cross-correlation as functions of time and latitude. Our simulated data do not include any wave-like characteristics for the supergranules yet can reproduce the rotation characteristics previously attributed to wave-like behavior. We find rotation rates which appear faster than the actual rotation rates and attribute this to the projection effects. We find that the measured meridional flow does accurately represent the actual flow and that the observations indicate poleward flow to $65\degr-70\degr$ latitude with equatorward counter cells in the polar regions.

We produce a 10-day series of simulated Doppler images at a 15-minute cadence that reproduces the spatial and temporal characteristics seen in the SOHO/MDI Doppler data. Our simulated data contains a spectrum of cellular flows with but two necessary components --- a granule component that peaks at wavenumbers of about 4000 and a supergranule component that peaks at wavenumbers of about 110. We include the advection of these cellular components by a differential rotation profile that depends on latitude and wavenumber (depth). We further mimic the evolution of the cellular pattern by introducing random variations to the amplitudes and phases of the spectral components at rates that reproduce the level of cross-correlation as a function of time and latitude. Our simulated data do not include any wave-like characteristics for the supergranules yet can accurately reproduce the rotation characteristics previously attributed to wave-like characteristics.

This deep, extended solar minimum and the slow start to Cycle 24 strongly suggest that Cycle 24 will be a small cycle. A wide array of solar cycle prediction techniques have been applied to predicting the amplitude of Cycle 24 with widely different results. Current conditions and new observations indicate that some highly regarded techniques now appear to have doubtful utility. Geomagnetic precursors have been reliable in the past and can be tested with 12 cycles of data. Of the three primary geomagnetic precursors only one (the minimum level of geomagnetic activity) suggests a small cycle. The Sun's polar field strength has also been used to successfully predict the last three cycles. The current weak polar fields are indicative of a small cycle. For the first time, dynamo models have been used to predict the size of a solar cycle but with opposite predictions depending on the model and the data assimilation. However, new measurements of the surface meridional flow indicate that the flow was substantially faster on the approach to Cycle 24 minimum than at Cycle 23 minimum. In both dynamo predictions a faster meridional flow should have given a shorter cycle 23 with stronger polar fields. This suggests that these dynamo models are not yet ready for solar cycle prediction.