We observe the abrupt end of solar activity cycles at the Sun's equator by combining almost 140 years of observations from ground and space. These "terminator" events appear to be very closely related to the onset of magnetic activity belonging to the next sunspot cycle at mid-latitudes and the polar-reversal process at high-latitudes. Using multi-scale tracers of solar activity we examine the timing of these events in relation to the excitation of new activity and find that the time taken for the solar plasma to communicate this transition is of the order of one solar rotation, but could be shorter. Utilizing uniquely comprehensive solar observations from the Solar Terrestrial Relations Observatory (STEREO), and Solar Dynamics Observatory (SDO) we see that this transitional event is strongly longitudinal in nature. Combined, these characteristics imply that magnetic information is communicated through the solar interior rapidly. A range of possibilities exist to explain such behavior: the presence of magnetic reconnection in the deep interior, internal gravity waves on the solar tachocline, or that the magnetic fields present in the Sun's convection zone could be very large, with a poloidal field strengths reaching 50k - considerably larger than conventional explorations of solar and stellar dynamos estimate. Regardless of mechanism responsible, the rapid timescales demonstrated by the Sun's global magnetic field reconfiguration present strong constraints on first-principles numerical simulations of the solar interior and, by extension, other stars.

In the recent years there has been a resurgence of the study of Extended Solar Cycles (ESCs) through observational proxies mainly in Extreme Ultraviolet. But most of them are limited only to space-based era covering only about two solar cycles. Long-term historical data-sets are worth in examining the consistency of ESCs. Kodaikanal Solar Observatory (KSO) and Mount Wilson Observatory (MWO) are the two major sources of long-term Ca II K digitised spectroheliograms covering the temporal spans 1907-2007 and 1915-1985 respectively. In this study, we detected supergranule boundaries, commonly known as networks, using the Carrington maps from both KSO and MWO datasets. Subsequently we excluded the plage areas to consider only quiet sun (QS) and detected small scale bright features through intensity thresholding over the QS network. Latitudinal density of those features, which we named as `Network Bright Elements' (NBEs), could clearly depict the existence of overlapping cycles with equator-ward branches starting at latitude $\approx 55^{\circ}$ and taking about $15\pm1$ years to reach the equator. We performed superposed epoch analysis to depict the similarity of those extended cycles. Knowledge of such equator-ward band interaction, for several cycles, may provide critical constraints on solar dynamo models.

In 1844 Schwabe discovered that the number of sunspots increased and decreased over a period of about 11 years, that variation became known as the sunspot cycle. Almost eighty years later, Hale described the nature of the Sun's magnetic field, identifying that it takes about 22 years for the Sun's magnetic polarity to cycle. It was also identified that the latitudinal distribution of sunspots resembles the wings of a butterfly showing migration of sunspots in each hemisphere that abruptly start at mid-latitudes towards the Sun's equator over the next 11 years. These sunspot patterns were shown to be asymmetric across the equator. In intervening years, it was deduced that the Sun (and sun-like stars) possess magnetic activity cycles that are assumed to be the physical manifestation of a dynamo process that results from complex circulatory transport processes in the star's interior. Understanding the Sun's magnetism, its origin and its variation, has become a fundamental scientific objective \-- the distribution of magnetism, and its interaction with convective processes, drives various plasma processes in the outer atmosphere. In the past few decades, a range of diagnostic techniques have been employed to systematically study finer scale magnetized objects, and associated phenomena. The patterns discerned became known as the ``Extended Solar Cycle'' (ESC). The patterns of the ESC appeared to extend the wings of the activity butterfly back in time, nearly a decade before the formation of the sunspot pattern, and to much higher solar latitudes. In this short review, we describe their observational patterns of the ESC and discuss possible connections to the solar dynamo as we depart on a multi-national collaboration to investigate the origins of solar magnetism through a blend of archived and contemporary data analysis with the goal of improving solar dynamo understanding and modeling.

Rieger-type periodicity has been detected in different activity indices over many solar cycles. It was recently shown that the periodicity correlates with solar activity having a shorter period during stronger cycles. Solar activity level is generally asymmetric between northern and southern hemispheres, which could suggest the presence of a similar behavior in the Rieger-type periodicity. We analyse the sunspot area/number and the total magnetic flux data for northern and southern hemispheres during solar cycles 19-23 which had remarkable north-south asymmetry. Using wavelet analysis of sunspot area and number during the north-dominated cycles (19-20) we obtained the periodicity of 160-165 days in the stronger northern hemisphere and 180-190 days in the weaker southern hemisphere. On the other hand, south-dominated cycles (21-23) display the periodicity of 155-160 days in the stronger southern hemisphere and 175-188 days in the weaker northern hemisphere. Therefore, the Rieger-type periodicity has the north-south asymmetry in sunspot area/number data during solar cycles with strong hemispheric asymmetry. We suggest that the periodicity is caused by magnetic Rossby waves in the internal dynamo layer. Using the dispersion relation of magnetic Rossby waves and observed Rieger periodicity we estimated the magnetic field strength in the layer as 45-49 kG in more active hemispheres (north during the cycles 19-20 and south during the cycles 21-23) and 33-40 kG in weaker hemispheres. The estimated difference in the hemispheric field strength is around 10 kG, which provides a challenge for dynamo models. Total magnetic flux data during the cycle 20-23 reveals no clear north-south asymmetry which needs to be explained in the future.

A brief summary of the various observations and constraints that underlie solar dynamo research are presented. The arguments that indicate that the solar dynamo is an alpha-omega dynamo of the Babcock-Leighton type are then shortly reviewed. The main open questions that remain are concerned with the subsurface dynamics, including why sunspots emerge at preferred latitudes as seen in the familiar butterfly wings, why the cycle is about 11 years long, and why the sunspot groups emerge tilted with respect to the equator (Joy's law). Next, we turn to magnetic helicity, whose conservation property has been identified with the decline of large-scale magnetic fields found in direct numerical simulations at large magnetic Reynolds numbers. However, magnetic helicity fluxes through the solar surface can alleviate this problem and connect theory with observations, as will be discussed.

Babcock-Leighton type solar dynamo models with single-celled meridional circulation are successful in reproducing many solar cycle features. Recent observations and theoretical models of meridional circulation do not indicate a single-celled flow pattern. We examine the role of complex multi-cellular circulation patterns in a Babcock-Leighton solar dynamo in advection- and diffusion-dominated regimes. We show from simulations that presence of a weak, second, high-latitude reverse cell speeds up the cycle and slightly enhances the poleward branch in butterfly diagram, whereas the presence of a second cell in depth reverses the tilt of butterfly wing to an anti-solar type. A butterfly diagram constructed from middle of convection zone yields a solar-like pattern, but this may be difficult to realize in the Sun because of magnetic buoyancy effects. Each of the above cases behaves similarly in higher and lower magnetic diffusivity regimes. However, our dynamo with a meridional circulation containing four cells in latitude behaves distinctly differently in the two regimes, producing solar-like butterfly diagrams with fast cycles in the higher diffusivity regime, and complex branches in butterfly diagrams in the lower diffusivity regime. We also find that dynamo solutions for a four-celled pattern, two in radius and two in latitude, prefer to quickly relax to quadrupolar parity if the bottom flow-speed is strong enough, of similar order of magnitude as the surface flow-speed.

The peculiar behaviour of the solar cycle 23 and its prolonged minima has been one of the most studied problems over the last few years. In the present paper, we study the asymmetries in active region magnetic flux in the northern and southern hemispheres during complete solar cycle 23 and rising phase of solar cycle 24. During the declining phase of solar cycle 23, we find that the magnetic flux in the southern hemisphere is about 10 times stronger than that in the northern hemisphere during the declining phase of the solar cycle 23 and during the rising phase of cycle 24, however, this trend reversed. The magnetic flux becomes about a factor of 4 stronger in the northern hemisphere to that of southern hemisphere. Additionally, we find that there was significant delay (about 5 months) in change of the polarity in the southern hemisphere in comparison with the northern hemisphere. These results provide us with hints of how the toroidal fluxes have contributed to the solar dynamo during the prolonged minima in the solar cycle 23 and in the rising phase of the solar cycle 24. Using a solar flux-transport dynamo model, we demonstrate that persistently stronger sunspot cycles in one hemisphere could be caused by the effect of greater inflows into active region belts in that hemisphere. Observations indicate that greater inflows are associated with stronger activity. Some other change or difference in meridional circulation between hemispheres could cause the weaker hemisphere to become the stronger one.

Accurate knowledge of time-variation in meridional flow-speed and profile is crucial for estimating a solar cycle's features, which are ultimately responsible for causing space climate variations. However, no consensus has been reached yet about the Sun's meridional circulation pattern observations and theories. By implementing an Ensemble Kalman Filter (EnKF) data assimilation in a Babcock-Leighton solar dynamo model using Data Assimilation Research Testbed (DART) framework, we find that the best reconstruction of time-variation in meridional flow-speed can be obtained when ten or more observations are used with an updating time of 15 days and a $\le 10\%$ observational error. Increasing ensemble-size from 16 to 160 improves reconstruction. Comparison of reconstructed flow-speed with "true-state" reveals that EnKF data assimilation is very powerful for reconstructing meridional flow-speeds and suggests that it can be implemented for reconstructing spatio-temporal patterns of meridional circulation.

We present a 3D kinematic solar dynamo model in which poloidal field is generated by the emergence and dispersal of tilted sunspot pairs (more generally Bipolar Magnetic Regions, or BMRs). The axisymmetric component of this model functions similarly to previous 2D Babcock-Leighton (BL) dynamo models that employ a double-ring prescription for poloidal field generation but we generalize this prescription into a 3D flux emergence algorithm that places BMRs on the surface in response to the dynamo-generated toroidal field. In this way, the model can be regarded as a unification of BL dynamo models (2D in radius/latitude) and surface flux transport models (2D in latitude/longitude) into a more self-consistent framework that captures the full 3D structure of the evolving magnetic field. The model reproduces some basic features of the solar cycle including an 11-yr periodicity, equatorward migration of toroidal flux in the deep convection zone, and poleward propagation of poloidal flux at the surface. The poleward-propagating surface flux originates as trailing flux in BMRs, migrates poleward in multiple non-axisymmetric streams (made axisymmetric by differential rotation and turbulent diffusion), and eventually reverses the polar field, thus sustaining the dynamo. In this letter we briefly describe the model, initial results, and future plans.

Meridional circulation in stellar convection zones is not generally well observed, but may be critical for MHD dynamos. Coriolis forces from differential rotation (DR) play a large role in determining what the meridional circulation is. Here we consider whether a stellar DR that is constant on cylinders concentric with the rotation axis can drive a meridional circulation.Conventional wisdom says that it can not. Using two related forms of governing equations that respectively estimate the longitudinal components of the curl of meridional mass flux and the vorticity, we show that such DR will drive a meridional flow. This is because to satisfy anelastic mass conservation, non-spherically symmetric pressure contours must be present for all DRs, not just ones that depart from constancy on cylinders concentric with the rotation axis. Therefore the fluid is always baroclinic if DR is present, because, in anelastic systems, the perturbation pressure must satisfy a Poisson type equation, as well as an equation of state and a thermodynamic equation. We support our qualitative reasoning with numerical examples, and show that meridional circulation is sensitive to the magnitude and form of departures from rotation constant on cylinders. The effect should be present in 3D global anelastic convection simulations, particularly those for which the DR driven by global convection is nearly cylindrical in profile. For solar-like DR, Coriolis forces generally drive a two-celled circulation in each hemisphere, with a second, reversed flow at high latitudes. For solar like turbulent viscosities, the meridional circulation produced by Coriolis forces is much larger than observed on the Sun. Therefore there must be at least one additional force, probably a buoyancy force, which opposes the meridional flow to bring its amplitude down to observed values.

We estimate here a flux-transport dynamo model's response time to changes in meridional flow speed. Time-variation in meridional flow primarily determines the shape of a cycle in this class of dynamo models. In order to simultaneously predict the shape, amplitude and timing of a solar cycle by implementing an Ensemble Kalman Filter in the framework of Data Assimilation Research Testbed (DART), it is important to know the model's sensitivity to flow variation. Guided by observations we consider a smooth increase or decrease in meridional flow speed for a specified time (a few months to a few years), after which the flow speed comes back to the steady speed, and implement that time-varying meridional flow at different phases of solar cycle. We find that the model's response time to change in flow speed peaks at four to six months if the flow change lasts for one year. The longer the changed flow lasts, the longer the model takes to respond. Magnetic diffusivity has no influence in model's response to flow variation as long as the dynamo operates in the advection-dominated regime. Experiments with more complex flow variations indicate that the shape and amplitude of flow-perturbation have no influence in the estimate of model's response time.

We build a hydrodynamical model for computing and understanding the Sun's large-scale high latitude flows, including Coriolis forces, turbulent diffusion of momentum and gyroscopic pumping. Side boundaries of the spherical 'polar cap', our computational domain, are located at latitudes $\geq 60^{\circ}$. Implementing observed low latitude flows as side boundary conditions, we solve the flow equations for a cartesian analog of the polar cap. The key parameter that determines whether there are nodes in the high latitude meridional flow is $\epsilon=2 \Omega n \pi H^2/\nu$, in which $\Omega$ is the interior rotation rate, n the radial wavenumber of the meridional flow, $H$ the depth of the convection zone and $\nu$ the turbulent viscosity. The smaller the $\epsilon$ (larger turbulent viscosity), the fewer the number of nodes in high latitudes. For all latitudes within the polar cap, we find three nodes for $\nu=10^{12}{\rm cm}^2{\rm s}^{-1}$, two for $10^{13}$, and one or none for $10^{15}$ or higher. For $\nu$ near $10^{14}$ our model exhibits 'node merging': as the meridional flow speed is increased, two nodes cancel each other, leaving no nodes. On the other hand, for fixed flow speed at the boundary, as $\nu$ is increased the poleward most node migrates to the pole and disappears, ultimately for high enough $\nu$ leaving no nodes. These results suggest that primary poleward surface meridional flow can extend from $60^{\circ}$ to the pole either by node-merging or by node migration and disappearance.

We present a fully nonlinear hydrodynamic 'shallow water' model of the solar tachocline. The model consists of a global spherical shell of differentially rotating fluid, which has a deformable top, thus allowing motions in radial directions along with latitude and longitude directions. When the system is perturbed, in the course of its nonlinear evolution it can generate unstable low-frequency shallow-water shear modes from the differential rotation, high-frequency gravity waves, and their interactions. Radiative and overshoot tachoclines are characterized in this model by high and low effective gravity values respectively. Building a semi-implicit spectral scheme containing very low numerical diffusion, we perform nonlinear evolution of shallow-water modes. Our first results show that, (i) high latitude jets or polar spin-up occurs due to nonlinear evolution of unstable hydrodynamic shallow-water disturbances and differential rotation, (ii) Reynolds stresses in the disturbances together with changing shell-thickness and meridional flow are responsible for the evolution of differential rotation, (iii) disturbance energy primarily remains concentrated in the lowest longitudinal wavenumbers, (iv) an oscillation in energy between perturbed and unperturbed states occurs due to evolution of these modes in a nearly dissipation-free system, and (v) disturbances are geostrophic, but occasional nonadjustment in geostrophic balance can occur, particularly in the case of high effective gravity, leading to generation of gravity waves. We also find that a linearly stable differential rotation profile remains nonlinearly stable.

Recent observations revealed that small magnetic elements abundant at the solar surface move poleward with a velocity which seems to be lower than the plasma velocity $U_{\theta}$. Guerrero et al. (2011) explained this discrepancy as a consequence of diffusive spreading of the magnetic elements due to a positive radial gradient of $|U_{\theta}|$. As the gradient's sign (inferred by local helioseismology) is still unclear, cases with a negative gradient are studied in this paper. Under this condition, the velocity of the magnetic tracers turns out to be larger than the plasma velocity, in disagreement with the observations. Alternative mechanisms for explaining them independently are proposed. For the turbulent magnetic pumping it is shown that it has to be unrealistically strong to reconcile the model with the observations.

We simulate the magnetic feature tracking (MFT) speed using advective-diffusive transport models in both one and two dimensions. By depositing magnetic bipolar regions at different latitudes at the Sun's surface and following their evolution for a prescribed meridional circulation and magnetic diffusivity profiles, we derive the MFT speed as a function of latitude. We find that in a one dimensional surface-transport model the simulated MFT speed at the surface is always the same as the meridional flow-speed used as input to the model, but is different in a two-dimensional transport model in the meridional (r,theta) plane. The difference depends on the value of the magnetic diffusivity and on the radial gradient of the latitudinal velocity. We have confirmed our results with two different codes in spherical and Cartesian coordinates.

Motivated by observations that only a very few stars have been found to have antisolar differential rotation, much weaker in amplitude than that of the Sun, we analyze the stability of antisolar and solar type latitudinal differential rotations in the tachoclines of typical F,G and K stars. We employ two three-dimensional thin shell models, one for a Boussinesq but non-hydrostatic system, the other for a hydrostatic but non-Boussinesq system. We find that, in general, the combination of toroidal field band and differential rotation is more unstable, and unstable for lower toroidal fields, for antisolar than for solar-type differential rotation. In the antisolar case, the instability is always found to weaken the differential rotation, even if the primary energy source for the instability is the magnetic field. This favors surface antisolar differential rotations in stars being weaker than solar types, if the instability in the tachocline is felt at the surface of the star. This is most likely to happen in F stars, whose convection zones are much thinner than they are in G and K stars. This effect could help explain why the antisolar differential rotations that have been found are very weak compared to the Sun.

We show that simple 2 and 3-layer flux-transport dynamos, when forced at the top by a poloidal source term, can produce a widely varying amplitude of toroidal field at the bottom, depending on how close the meridional flow speed of the bottom layer is to the propagation speed of the forcing applied above the top layer, and how close the amplitude of the $\alpha$-effect is to two values that give rise to a resonant response. This effect should be present in this class of dynamo model no matter how many layers are included. This result could have implications for the prediction of future solar cycles from the surface magnetic fields of prior cycles. It could be looked for in flux-transport dynamos that are more realistic for the Sun, done in spherical geometry with differential rotation, meridional flow and $\alpha$-effect that vary with latitude and time as well as radius. Because of these variations, if resonance occurs, it should be more localized in time, latitude and radius.

Polar fields in solar cycle 23 were about 50% weaker than those in cycle 22. The only theoretical models which have addressed this puzzle are surface transport models and flux-transport dynamo models. Comparing polar fields obtained from numerical simulations using surface flux transport models and flux-transport dynamo models, we show that both classes of models can explain the polar field features within the scope of the physics included in the respective models. In both models, how polar fields change as a result of changes in meridional circulation depends on the details of meridional circulation profile used. Using physical reasoning and schematics as well as numerical solutions from a flux-transport dynamo model, we demonstrate that polar fields are determined mostly by the strength of surface poloidal source provided by the decay of tilted, bipolar active regions. Profile of meridional flow with latitude and its changes with time have much less effect in flux-transport dynamo models than in surface transport models.

We show that systematic differences between surface Doppler and magnetic element tracking measures of solar meridional flow can be explained by the effects of surface turbulent magnetic diffusion. Feature-tracking speeds are lower than plasma speeds in low and mid-latitudes, because magnetic diffusion opposes poleward plasma flow in low-latitudes whereas it adds to plasma flow at high latitudes. Flux transport dynamo models must input plasma flow; the model-outputs yield estimates of the surface magnetic feature tracking speed. We demonstrate that the differences between plasma speed and magnetic pattern speed in a flux-transport dynamo are consistent with the observed difference between these speeds.

In the non-linear phase of a dynamo process, the back-reaction of the magnetic field upon the turbulent motion results in a decrease of the turbulence level and therefore in a suppression of both the magnetic field amplification (the alpha-quenching effect) and the turbulent magnetic diffusivity (the eta-quenching effect). While the former has been widely explored, the effects of eta-quenching in the magnetic field evolution have rarely been considered. In this work we investigate the role of the suppression of diffusivity in a flux-transport solar dynamo model that also includes a non-linear alpha-quenching term. Our results indicate that, although for alpha-quenching the dependence of the magnetic field amplification with the quenching factor is nearly linear, the magnetic field response to eta-quenching is non-linear and spatially non-uniform. We have found that the magnetic field can be locally amplified in this case, forming long-lived structures whose maximum amplitude can be up to ~2.5 times larger at the tachocline and up to ~2 times larger at the center of the convection zone than in models without quenching. However, this amplification leads to unobservable effects and to a worse distribution of the magnetic field in the butterfly diagram. Since the dynamo cycle period increases when the efficiency of the quenching increases, we have also explored whether the eta-quenching can cause a diffusion-dominated model to drift into an advection-dominated regime. We have found that models undergoing a large suppression of the diffusivity produce a strong segregation of magnetic fields that may lead to unsteady dynamo-oscillations. On the other hand, an initially diffusion-dominated model undergoing a small suppression in eta remains in the diffusion-dominated regime.

The solar tachocline, located at the interface between the latitude-dependent rotation of the convection zone and the rigid radiative interior, presents high gradients of angular velocity which are of particular interest for the models of the solar dynamo and angular momentum transport. Furthermore, latitudinal and temporal variations of the tachocline parameters, if any, are also of particular interest in order to constrain models. We present a review of some of the theories of the tachocline and their predictions that may be tested by helioseismology. We describethe methods for inferring the tachocline parameters from observations and the associated difficulties. A review of results previously obtained is given and an analysis of the new 6 years database of LOWL observations is presented which yields no compelling evidence of variations or general trend of the tachocline parameters during the ascending phase of the current solar cycle (1994-2000).