This textbook gives a perspective of heliophysics in a way that emphasizes universal processes from a perspective that draws attention to what provides Earth (and similar (exo-)planets) with a relatively stable setting in which life as we know it can thrive. The book is intended for students in physical sciences in later years of their university training and for beginning graduate students in fields of solar, stellar, (exo-)planetary, and planetary-system sciences.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

During the solar magnetic activity cycle the emergence latitudes of sunspots change, leading to the well-known butterfly diagram. This phenomenon is poorly understood for other stars as starspot latitudes are generally unknown. The related changes in starspot rotation rates caused by latitudinal differential rotation can however be measured. Using the set of 3093 Kepler stars with activity cycles identified by Reinhold et al. (2017), we aim to study the temporal change in starspot rotation rates over magnetic activity cycles, and how this relates to the activity level, mean rotation rate, and effective temperature of the star. We measure the photometric variability as a proxy for the magnetic activity and the spot rotation rate in each quarter over the duration of the Kepler mission. We phase-fold these measurements with the cycle period. We perform averages over stars with comparable mean rotation rates and effective temperature at fixed activity-cycle phases. We detect a clear correlation between the variation of activity level and the variation of the starspot rotation rate. The sign and amplitude of this correlation depends on the mean stellar rotation and, to a lesser extent, on the effective temperature. For slowly rotating stars (with periods between 15-28 days) the starspot rotation rates are clearly anti-correlated with the level of activity during the activity cycles. A transition is seen at periods of 10-15 days, where stars with effective temperature above 4200K instead show positive correlation. Our results can be interpreted in terms of a stellar "butterfly diagram", but these appear different from the Sun's as the starspot rotation rates are either in phase or anti-phase with the activity level. Alternatively, the activity cycles seen by Kepler are short (around 2.5 years) and may therefore be secondary cycles, perhaps analogous to the solar quasi-biennial oscillations.

The observed convective flows on the photosphere (e.g., supergranulation, granulation) play a key role in the Babcock-Leighton (BL) process to generate large-scale polar fields from sunspots fields. In most surface flux transport (SFT) and BL dynamo models, the dispersal and migration of surface fields are modeled as an effective turbulent diffusion. We present the first kinematic 3D FT/BL model to explicitly incorporate realistic convective flows based on solar observations. The results obtained are generally in good agreement with the observed surface flux evolution and with non-convective models that have a turbulent diffusivity on the order of $3 \times 10^{12}$ cm$^2$ s$^{-1}$ (300 km$^2$ s$^{-1}$). However, we find that the use of a turbulent diffusivity underestimates the dynamo efficiency, producing weaker mean fields and shorter cycle.

The Sun occasionally goes through Maunder-like extended grand minima when its magnetic activity drops considerably from the normal activity level for several decades. Many possible theories have been proposed to explain the origin of these minima. However, how the Sun managed to recover from such inactive phases every time is even more enigmatic. The Babcock--Leighton type dynamos, which are successful in explaining many features of the solar cycle remarkably well, are not expected to operate during grand minima due to the lack of a sufficient number of sunspots. In this Letter, we explore the question of how the Sun could recover from grand minima through the Babcock--Leighton dynamo. In our three-dimensional dynamo model, grand minima are produced spontaneously as a result of random variations in the tilt angle of emerging active regions. We find that the Babcock-Leighton process can still operate during grand minima with only a minimal number of sunspots and that the model can emerge from such phases without the need for an additional generation mechanism for the poloidal field. The essential ingredient in our model is a downward magnetic pumping which inhibits the diffusion of the magnetic flux across the solar surface.

We explore the cause of the solar cycle variabilities using a novel 3D Babcock--Leighton dynamo model. In this model, based on the toroidal flux at the base of the convection zone, bipolar magnetic regions (BMRs) are produced with statistical properties obtained from observed distributions. We find that a little quenching in BMR tilt is sufficient to stabilize the dynamo growth. The randomness and nonlinearity in the BMR emergences make the poloidal field unequal and cause some variability in the solar cycle. However, when observed scatter of BMR tilts around Joy's law with a standard deviation of $15^\circ$, is considered, our model produces a variation in the solar cycle, including north-south asymmetry comparable to the observations. The morphology of magnetic fields closely resembles observations, in particular the surface radial field possesses a more mixed polarity field. Observed scatter also produces grand minima. In 11,650 years of simulation, 17 grand minima are detected and $11\%$ of its time the model remained in these grand minima. When we double the tilt scatter, the model produces correct statistics of grand minima. Importantly, the dynamo continues even during grand minima with only a few BMRs, without requiring any additional alpha effect. The reason for this is the downward magnetic pumping which suppresses the diffusion of the magnetic flux across the surface. The magnetic pumping also helps to achieve 11-year magnetic cycle using the observed BMR flux distribution, even at the high diffusivity.

Observations suggest that the large-scale convective velocities obtained by solar convection simulations might be over-estimated (convective conundrum). One plausible solution to this could be the small-scale dynamo which cannot be fully resolved by global simulations. The small-scale Lorentz force suppresses the convective motions and also the turbulent mixing of entropy between upflows and downflows, leading to a large effective Prandtl number ($\Pr$). We explore this idea in three-dimensional global rotating convection simulations at different thermal conductivity ($\kappa$), i.e., at different $\Pr$. In agreement with previous non-rotating simulations, the convective velocity is reduced with the increase of $\Pr$ as long as the thermal conductive flux is negligible. A subadiabatic layer is formed near the base of the convection zone due to continuous deposition of low entropy plumes in low-$\kappa$ simulations. The most interesting result of our low-$\kappa$ simulations is that the convective motions are accompanied by a change in the convection structure that is increasingly influenced by small-scale plumes. These plumes tend to transport angular momentum radially inward and thus establish an anti-solar differential rotation, in striking contrast to the solar rotation profile. If such low diffusive plumes, driven by the radiative-surface cooling, are present in the Sun, then our results cast doubt on the idea that a high effective $\Pr$ may be a viable solution to the solar convective conundrum. Our study also emphasizes that any resolution of the conundrum that relies on the downward plumes must take into account angular momentum transport as well as heat transport.

The observed convective flows on the photosphere (e.g., supergranulation, granulation) play a key role in the Babcock-Leighton (BL) process to generate large-scale polar fields from sunspots fields. In most surface flux transport (SFT) and BL dynamo models, the dispersal and migration of surface fields is modeled as an effective turbulent diffusion. Recent SFT models have incorporated explicit, realistic convective flows in order to improve the fidelity of convective transport but, to our knowledge, this has not yet been implemented in previous BL models. Since most Flux-Transport (FT)/BL models are axisymmetric, they do not have the capacity to include such flows. We present the first kinematic 3D FT/BL model to explicitly incorporate realistic convective flows based on solar observations. Though we describe a means to generalize these flows to 3D, we find that the kinematic small-scale dynamo action they produce disrupts the operation of the cyclic dynamo. Cyclic solution is found by limiting the convective flow to act only on the vertical radial component of the magnetic field. The results obtained are generally in good agreement with the observed surface flux evolution and with non-convective models that have a turbulent diffusivity on the order of $3 \times 10^{12}$ cm$^2$ s$^{-1}$ (300 km$^2$ s$^{-1}$). However, we find that the use of a turbulent diffusivity underestimates the dynamo efficiency, producing weaker mean fields and shorter cycle than in the convective models. Also, the convective models exhibit mixed polarity bands in the polar regions that have no counterpart in solar observations. Also, the explicitly computed turbulent electromotive force (emf) bears little resemblance to a diffusive flux. We also find that the poleward migration speed of poloidal flux is determined mainly by the meridional flow and the vertical diffusion.

Global 3D simulations of solar giant-cell convection have provided significant insight into the processes which yield the Sun's observed differential rotation and cyclic dynamo action. However, as we move to higher resolution simulations a variety of codes have encountered what has been termed the convection conundrum. As these simulations increase in resolution and hence the level of turbulence achieved, they tend to produce weak or even anti-solar differential rotation patterns associated with a weak rotational influence (high Rossby number) due to large convective velocities. One potential culprit for this convection conundrum is the upper boundary condition applied in most simulations which is generally impenetrable. Here we present an alternative stochastic plume boundary condition which imposes small-scale convective plumes designed to mimic near-surface convective downflows, thus allowing convection to carry the majority of the outward solar energy flux up to and through our simulated upper boundary. The use of a plume boundary condition leads to significant changes in the convective driving realized in the simulated domain and thus to the convective energy transport, the dominant scale of the convective enthalpy flux, and the relative strength of the strongest downflows, the downflow network, and the convective upflows. These changes are present even far from the upper boundary layer. Additionally, we demonstrate that in spite of significant changes, giant cell morphology in the convective patterns is still achieved with self-organization of the imposed boundary plumes into downflow lanes, cellular patterns, and even rotationally-aligned banana cells in equatorial regions. This plume boundary presents an alternative pathway for 3D global convection simulations where driving is non-local and may provide a new approach towards addressing the convection conundrum.

We present results from a three-dimensional Babcock--Leighton dynamo model that is sustained by the explicit emergence and dispersal of bipolar magnetic regions (BMRs). On average, each BMR has a systematic tilt given by Joy's law. Randomness and nonlinearity in the BMR emergence of our model produce variable magnetic cycles. However, when we allow for a random scatter in the tilt angle to mimic the observed departures from Joy's law, we find more variability in the magnetic cycles. We find that the observed standard deviation in Joy's law of $\sigma_\delta = 15^\circ$ produces a variability comparable to observed solar cycle variability of $\sim $ 32%, as quantified by the sunspot number maxima between 1755--2008. We also find that tilt angle scatter can promote grand minima and grand maxima. The time spent in grand minima for $\sigma_\delta = 15^\circ$ is somewhat less than that inferred for the Sun from cosmogenic isotopes (about 9% compared to 17%). However, when we double the tilt scatter to $\sigma_\delta = 30^\circ$, the simulation statistics are comparable to the Sun ($\sim $18% of the time in grand minima and $\sim 10$% in grand maxima). Though the Babcock--Leighton mechanism is the only source of poloidal field, we find that our simulations always maintain magnetic cycles even at large fluctuations in the tilt angle. We also demonstrate that tilt quenching is a viable and efficient mechanism for dynamo saturation; a suppression of the tilt by only 1-2$^\circ$ is sufficient to limit the dynamo growth. Thus, any potential observational signatures of tilt quenching in the Sun may be subtle.

Surface observations indicate that the speed of the solar meridional circulation in the photosphere varies in anti-phase with the solar cycle. The current explanation for the source of this variation is that inflows into active regions alter the global surface pattern of the meridional circulation. When these localized inflows are integrated over a full hemisphere, they contribute to the slow down of the axisymmetric poleward horizontal component. The behavior of this large scale flow deep inside the convection zone remains largely unknown. Present helioseismic techniques are not sensitive enough to capture the dynamics of this weak large scale flow. Moreover, the large time of integration needed to map the meridional circulation inside the convection zone, also masks some of the possible dynamics on shorter timescales. In this work we examine the dynamics of the meridional circulation that emerges from a 3D MHD global simulation of the solar convection zone. Our aim is to assess and quantify the behavior of meridional circulation deep inside the convection zone, where the cyclic large-scale magnetic field can reach considerable strength. Our analyses indicate that the meridional circulation morphology and amplitude are both highly influenced by the magnetic field, via the impact of magnetic torques on the global angular momentum distribution. A dynamic feature induced by these magnetic torques is the development of a prominent upward flow at mid latitudes in the lower convection zone that occurs near the equatorward edge of the toroidal bands and that peaks during cycle maximum. Globally, the dynamo-generated large-scale magnetic field drives variations in the meridional flow, in stark contrast to the conventional kinematic flux transport view of the magnetic field being advected passively by the flow.

We develop a three-dimensional kinematic self-sustaining model of the solar dynamo in which the poloidal field generation is from tilted bipolar sunspot pairs placed on the solar surface above regions of strong toroidal field by using the SpotMaker algorithm, and then the transport of this poloidal field to the tachocline is primarily caused by turbulent diffusion. We obtain a dipolar solution within a certain range of parameters. We use this model to study the build-up of the polar magnetic field and show that some insights obtained from surface flux transport (SFT) models have to be revised. We present results obtained by putting a single bipolar sunspot pair in a hemisphere and two symmetrical sunspot pairs in two hemispheres.We find that the polar fields produced by them disappear due to the upward advection of poloidal flux at low latitudes, which emerges as oppositely-signed radial flux and which is then advected poleward by the meridional flow. We also study the effect that a large sunspot pair, violating Hale's polarity law would have on the polar field. We find that there would be some effect---especially if the anti-Hale pair appears at high latitudes in the mid-phase of the cycle---though the effect is not very dramatic.

We present the first global 3D simulations of thermal convection in the oblate envelopes of rapidly-rotating solar-type stars. This has been achieved by exploiting the capabilities of the new Compressible High-ORder Unstructured Spectral difference (CHORUS) code. We consider rotation rates up to 85\% of the critical (breakup) rotation rate, which yields an equatorial radius that is up to 17\% larger than the polar radius. This substantial oblateness enhances the disparity between polar and equatorial modes of convection. We find that the convection redistributes the heat flux emitted from the outer surface, leading to an enhancement of the heat flux in the polar and equatorial regions. This finding implies that lower-mass stars with convective envelopes may not have darker equators as predicted by classical gravity darkening arguments. The vigorous high-latitude convection also establishes elongated axisymmetric circulation cells and zonal jets in the polar regions. Though the overall amplitude of the surface differential rotation, $\Delta \Omega$, is insensitive to the oblateness, the oblateness does limit the fractional kinetic energy contained in the differential rotation to no more than 61\%. Furthermore, we argue that this level of differential rotation is not enough to have a significant impact on the oblateness of the star.

We investigate the role of magnetic helicity in promoting cyclic magnetic activity in a global, 3D, magnetohydrodynamic (MHD) simulation of a convective dynamo. This simulation is characterized by coherent bands of toroidal field that exist within the convection zone, with opposite polarities in the northern and southern hemispheres. Throughout most of the cycle, the magnetic helicity in these bands is negative in the northern hemisphere and positive in the southern hemisphere. However, during the declining phase of each cycle, this hemispheric rule reverses. We attribute this to a global restructuring of the magnetic topology that is induced by the interaction of the bands across the equator. This band interaction appears to be ultimately responsible for, or at least associated with, the decay and subsequent reversal of both the toroidal bands and the polar fields. We briefly discuss the implications of these results within the context of solar observations, which also show some potential evidence for toroidal band interactions and helicity reversals.

The main objective of this paper is to introduce the STABLE (Surface flux Transport And Babcock-LEighton) solar dynamo model. STABLE is a 3D Babcock-Leighton/Flux Transport dynamo model in which the source of poloidal field is the explicit emergence, distortion, and dispersal of bipolar magnetic regions (BMRs). Here we describe the STABLE model in more detail than we have previously and we verify it by reproducing a 2D mean-field benchmark. We also present some representative dynamo simulations, focusing on the special case of kinematic magnetic induction and axisymmetric flow fields. Not all solutions are supercritical; it can be a challenge for the BL mechanism to sustain the dynamo when the turbulent diffusion near the surface is $\geq 10^{12}$ cm$^2$ s$^{-1}$. However, if BMRs are sufficiently large, deep, and numerous, then sustained, cyclic, dynamo solutions can be found that exhibit solar-like features. Furthermore, we find that the shearing of radial magnetic flux by the surface differential rotation can account for most of the net toroidal flux generation in each hemisphere, as has been recently argued for the Sun by Cameron & Schussler (2015).

Several lines of evidence suggest that the velocity amplitude in global simulations of solar convection, U, may be systematically over-estimated. Motivated by these recent results, we explore the factors that determine U and we consider how these might scale to solar parameter regimes. To this end, we decrease the thermal diffusivity $\kappa$ along two paths in parameter space. If the kinematic viscosity $\nu$ is decreased proportionally with $\kappa$ (fixing the Prandtl number $P_r = \nu/\kappa$), we find that U increases but asymptotes toward a constant value, as found by Featherstone & Hindman (2016). However, if $\nu$ is held fixed while decreasing $\kappa$ (increasing $P_r$), we find that U systematically decreases. We attribute this to an enhancement of the thermal content of downflow plumes, which allows them to carry the solar luminosity with slower flow speeds. We contrast this with the case of Rayleigh-Benard convection which is not subject to this luminosity constraint. This dramatic difference in behavior for the two paths in parameter space (fixed $P_r$ or fixed $\nu$) persists whether the heat transport by unresolved, near-surface convection is modeled as a thermal conduction or as a fixed flux. The results suggest that if solar convection can operate in a high-$P_r$ regime, then this might effectively limit the velocity amplitude. Small-scale magnetism is a possible source of enhanced viscosity that may serve to achieve this high-$P_r$ regime.

The Wisconsin Plasma Astrophysics Laboratory (WiPAL) is a flexible user facility designed to study a range of astrophysically relevant plasma processes as well as novel geometries that mimic astrophysical systems. A multi-cusp magnetic bucket constructed from strong samarium cobalt permanent magnets now confines a 10 m$^3$, fully ionized, magnetic-field free plasma in a spherical geometry. Plasma parameters of $ T_{e}\approx5$ to $20$ eV and $n_{e}\approx10^{11}$ to $5\times10^{12}$ cm$^{-3}$ provide an ideal testbed for a range of astrophysical experiments including self-exciting dynamos, collisionless magnetic reconnection, jet stability, stellar winds, and more. This article describes the capabilities of WiPAL along with several experiments, in both operating and planning stages, that illustrate the range of possibilities for future users.

We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions.

The 3-D magnetohydrodynamic (MHD) Anelastic Spherical Harmonic (ASH) code, using slope-limited diffusion, is employed to capture convective and dynamo processes achieved in a global-scale stellar convection simulation for a model solar-mass star rotating at three times the solar rate. The dynamo generated magnetic fields possesses many time scales, with a prominent polarity cycle occurring roughly every 6.2 years. The magnetic field forms large-scale toroidal wreaths, whose formation is tied to the low Rossby number of the convection in this simulation. The polarity reversals are linked to the weakened differential rotation and a resistive collapse of the large-scale magnetic field. An equatorial migration of the magnetic field is seen, which is due to the strong modulation of the differential rotation rather than a dynamo wave. A poleward migration of magnetic flux from the equator eventually leads to the reversal of the polarity of the high-latitude magnetic field. This simulation also enters an interval with reduced magnetic energy at low latitudes lasting roughly 16 years (about 2.5 polarity cycles), during which the polarity cycles are disrupted and after which the dynamo recovers its regular polarity cycles. An analysis of this grand minimum reveals that it likely arises through the interplay of symmetric and antisymmetric dynamo families. This intermittent dynamo state potentially results from the simulations relatively low magnetic Prandtl number. A mean-field-based analysis of this dynamo simulation demonstrates that it is of the $\alpha$-$\Omega$ type. The time scales that appear to be relevant to the magnetic polarity reversal are also identified.

We discuss recent observational, theoretical and modeling progress made in understanding the Sun's internal dynamics, including its rotation, meridional flow, convection and overshoot. Over the past few decades, substantial theoretical and observational effort has gone into appreciating these aspects of solar dynamics. A review of these observations, related helioseismic methodology and inference and computational results in relation to these problems is undertaken here.

The global-scale dynamo action achieved in a simulation of a Sun-like star rotating at thrice the solar rate is assessed. The 3-D MHD Anelastic Spherical Harmonic (ASH) code, augmented with a viscosity minimization scheme, is employed to capture convection and dynamo processes in this G-type star. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the solar convection zone. It is found that dynamo action with a high degree of time variation occurs, with many periodic polarity reversals occurring roughly every 6.2 years. The magnetic energy also rises and falls with a regular period. The magnetic energy cycles arise from a Lorentz-force feedback on the differential rotation, whereas the processes leading to polarity reversals are more complex, appearing to arise from the interaction of convection with the mean toroidal fields. Moreover, an equatorial migration of toroidal field is found, which is linked to the changing differential rotation, and potentially to a nonlinear dynamo wave. This simulation also enters a grand minimum lasting roughly 20~years, after which the dynamo recovers its regular polarity cycles.

We present a novel and powerful Compressible High-ORder Unstructured Spectral-difference (CHORUS) code for simulating thermal convection and related fluid dynamics in the interiors of stars and planets. The computational geometries are treated as rotating spherical shells filled with stratified gas. The hydrodynamic equations are discretized by a robust and efficient high-order Spectral Difference Method (SDM) on unstructured meshes. The computational stencil of the spectral difference method is compact and advantageous for parallel processing. CHORUS demonstrates excellent parallel performance for all test cases reported in this paper, scaling up to 12,000 cores on the Yellowstone High-Performance Computing cluster at NCAR. The code is verified by defining two benchmark cases for global convection in Jupiter and the Sun. CHORUS results are compared with results from the ASH code and good agreement is found. The CHORUS code creates new opportunities for simulating such varied phenomena as multi-scale solar convection, core convection, and convection in rapidly-rotating, oblate stars.

The form of the solar meridional circulation is a very important ingredient for mean field flux transport dynamo models. Yet a shroud of mystery still surrounds this large-scale flow, given that its measurement using current helioseismic techniques is challenging. In this work we use results from 3D global simulations of solar convection to infer the dynamical behavior of the established meridional circulation. We make a direct comparison between the meridional circulation that arises in these simulations and the latest observations. Based on our results we argue that there should be an equatorward flow at the base of the convection zone at mid latitudes, below the current maximum depth helioseismic measures can probe (0.75 R). We also provide physical arguments to justify this behaviour. The simulations indicate that the meridional circulation undergoes substantial changes in morphology as the magnetic cycle unfolds. We close by discussing the importance of these dynamical changes for current methods of observation that involve long averaging periods of helioseismic data. Also noteworthy is the fact that these topological changes indicate a rich interaction between magnetic fields and plasma flows, which challenges the ubiquitous kinematic approach used in the vast majority of mean field dynamo simulations.

We present a series of 3-D nonlinear simulations of solar-like convection, carried out using the Anelastic Spherical Harmonic (ASH) code, that are designed to isolate those processes that drive and shape meridional circulations within stellar convection zones. These simulations have been constructed so as to span the transition between solar-like differential rotation (fast equator/slow poles) and ``anti-solar' differential rotation (slow equator/fast poles). Solar-like states of differential rotation, arising when convection is rotationally constrained, are characterized by a very different convective Reynolds stress than anti-solar regimes, wherein convection only weakly senses the Coriolis force. We find that the angular momentum transport by convective Reynolds stress plays a central role in establishing the meridional flow profiles in these simulations. We find that the transition from single-celled to multi-celled meridional circulation profiles in strong and weak regimes of rotational constraint is linked to a change in the convective Reynolds stress, a clear demonstration of gyroscopic pumping. Latitudinal thermal variations differ between these different regimes, with those in the solar-like regime conspiring to suppress a single cell of meridional circulation, whereas the cool poles and warm equator established in the anti-solar states tend to promote single-celled circulations. Though the convective angular momentum transport becomes radially inward at mid-latitudes in anti-solar regimes, it is the meridional circulation that is primarily responsible for establishing a rapidly-rotating pole. We conclude with a discussion of how these results relate to the Sun, and suggest that the Sun may lie near the transition between rapidly-rotating and slowly-rotating regimes.

Understanding differential rotation of Sun-like stars is of great importance for insight into the angular momentum transport in these stars. One means of gaining such information is that of asteroseismology. By a forward modeling approach we analyze in a qualitative manner the impact of different differential rotation profiles on the splittings of p-mode oscillation frequencies. The optimum modes for inference on differential rotation are identified along with the best value of the stellar inclination angle. We find that in general it is not likely that asteroseismology can be used to make an unambiguous distinction between a rotation profile such as, e.g., a conical Sun-like profile and a cylindrical profile. In addition, it seems unlikely that asteroseismology of Sun-like stars will result in inferences on the radial profile of the differential rotation, such as can be done for, e.g., red giants. At best one could possibly obtain the sign of the radial differential rotation gradient. Measurements of the extent of the latitudinal differential from frequency splitting are, however, more promising. One very interesting aspect that could likely be tested from frequency splittings is whether the differential rotation is solar-like or anti-solar-like in nature, in the sense that a solar-like profile has an equator rotating faster than the poles.

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.

Our Sun exhibits strong convective dynamo action which results in magnetic flux bundles emerging through the stellar surface as magnetic spots. Global-scale dynamo action is believed to generate large-scale magnetic structures in the deep solar interior through the interplay of convection, rotation, and shear. Portions of these large-scale magnetic structures are then believed to rise through the convective layer, forming magnetic loops which then pierce the photosphere as sunspot pairs. Previous global simulations of 3D MHD convection in rotating spherical shells have demonstrated mechanisms whereby large-scale magnetic wreaths can be generated in the bulk of the convection zone. Our recent simulations have achieved sufficiently high levels of turbulence to permit portions of these wreaths to become magnetically buoyant and rise through the simulated convective layer through a combination of magnetic buoyancy and advection by convective giant cells. These buoyant magnetic loops are created in the bulk of the convective layer as strong Lorentz force feedback in the cores of the magnetic wreaths dampen small-scale convective motions, permitting the amplification of local magnetic energies to over 100 times the local kinetic energy. While the magnetic wreaths are largely generated the shearing of axisymmetric poloidal magnetic fields by axisymmetric rotational shear (the $\Omega$-effect), the loops are amplified to their peak field strengths before beginning to rise by non-axisymmetric processes. This further extends and enhances a new paradigm for the generation of emergent magnetic flux bundles, which we term turbulence-enabled magnetic buoyancy.

We assess the global-scale dynamo action achieved in a simulation of a sun-like star rotating at three times the solar rate. The 3-D MHD Anelastic Spherical Harmonic code, using slope-limited diffusion, is employed to capture convection and dynamo processes in such a young sun. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the solar convection zone. We find that dynamo action with a high degree of time variation occurs, with many periodic polarity reversals every 6.2 years. The magnetic energy also rises and falls with a regular period, with two magnetic energy cycles required to complete a polarity cycle. These magnetic energy cycles arise from a Lorentz-force feedback on the differential rotation, whereas the polarity reversals are present due to the spatial separation of the equatorial and polar dynamos. Moreover, an equatorial migration of toroidal field is found, which is linked to the changing differential rotation and to a near-surface shear layer. This simulation also enters a grand minimum lasting roughly 20 years, after which the dynamo recovers its regular polarity cycles.

Using a thin flux tube model in a rotating spherical shell of turbulent, solar-like convective flows, we find that the distribution of emerging flux tubes in our simulation is inhomogeneous in longitude, with properties similar to those of active longitudes on the Sun and other solar-like stars. The large-scale pattern of flux emergence our simulations produce exhibits preferred longitudinal modes of low order, drift with respect to a fixed reference system, and alignment across the Equator at low latitudes between 15 degrees. We suggest that these active-longitude-like emergence patterns are the result of columnar, rotationally aligned giant cells present in our convection simulation at low latitudes. If giant convecting cells exist in the bulk of the solar convection zone, this phenomenon, along with differential rotation, could in part provide an explanation for the behavior of active longitudes.

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.

Our global 3D simulations of convection and dynamo action in a Sun-like star reveal that persistent wreaths of strong magnetism can be built within the bulk of the convention zone. Here we examine the characteristics of buoyant magnetic structures that are self-consistently created by dynamo action and turbulent convective motions in a simulation with solar stratification but rotating at three times the current solar rate. These buoyant loops originate within sections of the magnetic wreaths in which turbulent flows amplify the fields to much larger values than is possible through laminar processes. These amplified portions can rise through the convective layer by a combination of magnetic buoyancy and advection by convective giant cells, forming buoyant loops. We measure statistical trends in the polarity, twist, and tilt of these loops. Loops are shown to preferentially arise in longitudinal patches somewhat reminiscent of active longitudes in the Sun, although broader in extent. We show that the strength of the axisymmetric toroidal field is not a good predictor of the production rate for buoyant loops or the amount of magnetic flux in the loops that are produced.

Solar-type stars exhibit a rich variety of magnetic activity. Seeking to explore the convective origins of this activity, we have carried out a series of global 3D magnetohydrodynamic (MHD) simulations with the anelastic spherical harmonic (ASH) code. Here we report on the dynamo mechanisms achieved as the effects of artificial diffusion are systematically decreased. The simulations are carried out at a nominal rotation rate of three times the solar value (3$\Omega_\odot$), but similar dynamics may also apply to the Sun. Our previous simulations demonstrated that convective dynamos can build persistent toroidal flux structures (magnetic wreaths) in the midst of a turbulent convection zone and that high rotation rates promote the cyclic reversal of these wreaths. Here we demonstrate that magnetic cycles can also be achieved by reducing the diffusion, thus increasing the Reynolds and magnetic Reynolds numbers. In these more turbulent models, diffusive processes no longer play a significant role in the key dynamical balances that establish and maintain the differential rotation and magnetic wreaths. Magnetic reversals are attributed to an imbalance in the poloidal magnetic induction by convective motions that is stabilized at higher diffusion levels. Additionally, the enhanced levels of turbulence lead to greater intermittency in the toroidal magnetic wreaths, promoting the generation of buoyant magnetic loops that rise from the deep interior to the upper regions of our simulated domain. The implications of such turbulence-induced magnetic buoyancy for solar and stellar flux emergence are also discussed.

We obtain lower limits on the amplitude of convective velocities in the deep solar convection zone based only on the observed properties of the differential rotation and meridional circulation together with simple and robust dynamical balances obtained from the fundamental MHD equations. The linchpin of the approach is the concept of gyroscopic pumping whereby the meridional circulation across isosurfaces of specific angular momentum is linked to the angular momentum transport by the convective Reynolds stress. We find that the amplitude of the convective velocity must be at least 30 m s$^{-1}$ in the upper CZ ($r \sim 0.95 R$) and at least 8 m s$^{-1}$ in the lower CZ ($r \sim 0.75 R$) in order to be consistent with the observed mean flows. Using the base of the near-surface shear layer as a probe of the rotational influence, we are further able to show that the characteristic length scale of deep convective motions must be no smaller than 5.5--30 Mm. These results are compatible with convection models but suggest that the efficiency of the turbulent transport assumed in advection-dominated flux-transport dynamo models is generally not consistent with the mean flows they employ.

We study how active-region-scale flux tubes rise buoyantly from the base of the convection zone to near the solar surface by embedding a thin flux tube model in a rotating spherical shell of solar-like turbulent convection. These toroidal flux tubes that we simulate range in magnetic field strength from 15 kG to 100 kG at initial latitudes of 1 degree to 40 degrees in both hemispheres. This article expands upon Weber, Fan, and Miesch (Astrophys. J., 741, 11, 2011) (Article 1) with the inclusion of tubes with magnetic flux of 10^20 Mx and 10^21 Mx, and more simulations of the previously investigated case of 10^22 Mx, sampling more convective flows than the previous article, greatly improving statistics. Observed properties of active regions are compared to properties of the simulated emerging flux tubes, including: the tilt of active regions in accordance with Joy's Law as in Article 1, and in addition the scatter of tilt angles about the Joy's Law trend, the most commonly occurring tilt angle, the rotation rate of the emerging loops with respect to the surrounding plasma, and the nature of the magnetic field at the flux tube apex. We discuss how these diagnostic properties constrain the initial field strength of the active region flux tubes at the bottom of the solar convection zone, and suggest that flux tubes of initial magnetic field strengths of \geq 40 kG are good candidates for the progenitors of large (10^21 Mx to 10^22 Mx) solar active regions, which agrees with the results from Article 1 for flux tubes of 10^22 Mx. With the addition of more magnetic flux values and more simulations, we find that for all magnetic field strengths, the emerging tubes show a positive Joy's Law trend, and that this trend does not show a statistically significant dependence on the magnetic flux.

We present the first global, three-dimensional simulations of solar/stellar convection that take into account the influence of magnetic flux emergence by means of the Babcock-Leighton (BL) mechanism. We have shown that the inclusion of a BL poloidal source term in a convection simulation can promote cyclic activity in an otherwise steady dynamo. Some cycle properties are reminiscent of solar observations, such as the equatorward propagation of toroidal flux near the base of the convection zone. However, the cycle period in this young sun (rotating three times faster than the solar rate) is very short ($\sim$ 6 months) and it is unclear whether much longer cycles may be achieved within this modeling framework, given the high efficiency of field generation and transport by the convection. Even so, the incorporation of mean-field parameterizations in 3D convection simulations to account for elusive processes such as flux emergence may well prove useful in the future modeling of solar and stellar activity cycles.

We use a thin flux tube model in a rotating spherical shell of turbulent convective flows to study how active region scale flux tubes rise buoyantly from the bottom of the convection zone to near the solar surface. We investigate toroidal flux tubes at the base of the convection zone with field strengths ranging from 15 kG to 100 kG at initial latitudes ranging from 1 degree to 40 degrees with a total flux of 10^22 Mx. We find that the dynamic evolution of the flux tube changes from convection dominated to magnetic buoyancy dominated as the initial field strength increases from 15 kG to 100 kG. At 100 kG, the development of Omega-shaped rising loops is mainly controlled by the growth of the magnetic buoyancy instability. However, at low field strengths of 15 kG, the development of rising Omega-shaped loops is largely controlled by convective flows, and properties of the emerging loops are significantly changed compared to previous results in the absence of convection. With convection, rise times are drastically reduced (from years to a few months), loops are able to emerge at low latitudes, and tilt angles of emerging loops are consistent with Joy's Law for initial field strengths of greater than or equal to 40 kG. We also examine other asymmetries that develop between the leading and following legs of the emerging loops. Taking all the results together, we find that mid-range field strengths of approximately 40 - 50 kG produce emerging loops that best match the observed properties of solar active regions.

We use global and local helioseismic inversions to explore the prevailing dynamical balances in the solar Near-Surface Shear Layer (NSSL). The differential rotation and meridional circulation are intimately linked, with a common origin in the turbulent stresses of the upper solar convection zone. The existence and structure of the NSSL cannot be attributed to the conservation of angular momentum by solar surface convection, as is often supposed. Rather, the turbulent angular momentum transport accounts for the poleward meridional flow while the often overlooked meridional force balance maintains the mid-latitude rotational shear. We suggest that the base of the NSSL is marked by a transition from baroclinic to turbulent stresses in the meridional plane which suppress Coriolis-induced circulations that would otherwise establish a cylindrical rotation profile. The turbulent angular momentum transport must be non-diffusive and directed radially inward. Inferred mean flows are consistent with the idea that turbulent convection tends to mix angular momentum but only if the mixing efficiency is inhomogeneous and/or anisotropic. The latitudinal and longitudinal components of the estimated turbulent transport are comparable in amplitude and about an order of magnitude larger than the vertical component. We estimate that it requires 2--4% of the solar luminosity to maintain the solar NSSL against the inertia of the mean flow. Most of this energy is associated with the turbulent transport of angular momentum out of the layer, with a spin-down time scale of $\sim$ 600 days. We also address implications of these results for numerical modeling of the NSSL.

The current dynamo paradigm for the Sun and sun-like stars places the generation site for strong toroidal magnetic structures deep in the solar interior. Sunspots and star-spots on sun-like stars are believed to arise when sections of these magnetic structures become buoyantly unstable and rise from the deep interior to the photosphere. Here we present the first 3-D global magnetohydrodynamic (MHD) simulation in which turbulent convection, stratification, and rotation combine to yield a dynamo that self-consistently generates buoyant magnetic loops. We simulate stellar convection and dynamo action in a spherical shell with solar stratification, but rotating three times faster than the current solar rate. Strong wreaths of toroidal magnetic field are realized by dynamo action in the convection zone. By turning to a dynamic Smagorinsky model for subgrid-scale turbulence, we here attain considerably reduced diffusion in our simulation. This permits the regions of strongest magnetic field in these wreaths to rise toward the top of the convection zone via a combination of magnetic buoyancy instabilities and advection by convective giant cells. Such a global simulation yielding buoyant loops represents a significant step forward in combining numerical models of dynamo action and flux emergence.

Young solar-type stars rotate rapidly and many are magnetically active; some undergo magnetic cycles similar to the 22-year solar activity cycle. We conduct simulations of dynamo action in rapidly rotating suns with the 3D MHD anelastic spherical harmonic (ASH) code to explore dynamo action achieved in the convective envelope of a solar-type star rotating at 5 times the current solar rotation rate. Striking global-scale magnetic wreaths appear in the midst of the turbulent convection zone and show rich time-dependence. The dynamo exhibits cyclic activity and undergoes quasi-periodic polarity reversals where both the global-scale poloidal and toroidal fields change in sense on a roughly 1500 day time scale. These magnetic activity patterns emerge spontaneously from the turbulent flow and are more organized temporally and spatially than those realized in our previous simulations of the solar dynamo. We assess in detail the competing processes of magnetic field creation and destruction within our simulations that contribute to the global-scale reversals. We find that the mean toroidal fields are built primarily through an $\Omega$-effect, while the mean poloidal fields are built by turbulent correlations which are not necessarily well represented by a simple $\alpha$-effect. During a reversal the magnetic wreaths propagate towards the polar regions, and this appears to arise from a poleward propagating dynamo wave. The primary response in the convective flows involves the axisymmetric differential rotation which shows variations associated with the poleward propagating magnetic wreaths. In the Sun, similar patterns are observed in the poleward branch of the torsional oscillations, and these may represent poleward propagating magnetic fields deep below the solar surface. [abridged]

Young solar-type stars rotate rapidly and are very magnetically active. The magnetic fields at their surfaces likely originate in their convective envelopes where convection and rotation can drive strong dynamo action. Here we explore simulations of global-scale stellar convection in rapidly rotating suns using the 3-D MHD anelastic spherical harmonic (ASH) code. The magnetic fields built in these dynamos are organized on global-scales into wreath-like structures that span the convection zone. We explore one case rotates five times faster than the Sun in detail. This dynamo simulation, called case D5, has repeated quasi-cyclic reversals of global-scale polarity. We compare this case D5 to the broader family of simulations we have been able to explore and discuss how future simulations and observations can advance our understanding of stellar dynamos and magnetism.

When our Sun was young it rotated much more rapidly than now. Observations of young, rapidly rotating stars indicate that many possess substantial magnetic activity and strong axisymmetric magnetic fields. We conduct simulations of dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code to explore the complex coupling between rotation, convection and magnetism. Here we study dynamo action realized in the bulk of the convection zone for a system rotating at three times the current solar rotation rate. We find that substantial organized global-scale magnetic fields are achieved by dynamo action in this system. Striking wreaths of magnetism are built in the midst of the convection zone, coexisting with the turbulent convection. This is a surprise, for it has been widely believed that such magnetic structures should be disrupted by magnetic buoyancy or turbulent pumping. Thus, many solar dynamo theories have suggested that a tachocline of penetration and shear at the base of the convection zone is a crucial ingredient for organized dynamo action, whereas these simulations do not include such tachoclines. We examine how these persistent magnetic wreaths are maintained by dynamo processes and explore whether a classical mean-field $\alpha$-effect explains the regeneration of poloidal field.

When stars like our Sun are young they rotate rapidly and are very magnetically active. We explore dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code. The magnetic fields built in these dynamos are organized on global-scales into wreath-like structures that span the convection zone. Wreath-building dynamos can undergo quasi-cyclic reversals of polarity and such behavior is common in the parameter space we have been able to explore. These dynamos do not appear to require tachoclines to achieve their spatial or temporal organization. Wreath-building dynamos are present to some degree at all rotation rates, but are most evident in the more rapidly rotating simulations.

Observations of sun-like stars rotating faster than our current sun tend to exhibit increased magnetic activity as well as magnetic cycles spanning multiple years. Using global simulations in spherical shells to study the coupling of large-scale convection, rotation, and magnetism in a younger sun, we have probed effects of rotation on stellar dynamos and the nature of magnetic cycles. Major 3-D MHD simulations carried out at three times the current solar rotation rate reveal hydromagnetic dynamo action that yields wreaths of strong toroidal magnetic field at low latitudes, often with opposite polarity in the two hemispheres. Our recent simulations have explored behavior in systems with considerably lower diffusivities, achieved with sub-grid scale models including a dynamic Smagorinsky treatment of unresolved turbulence. The lower diffusion promotes the generation of magnetic wreaths that undergo prominent temporal variations in field strength, exhibiting global magnetic cycles that involve polarity reversals. In our least diffusive simulation, we find that magnetic buoyancy coupled with advection by convective giant cells can lead to the rise of coherent loops of magnetic field toward the top of the simulated domain.

We review recent insights into the dynamics of the solar convection zone obtained from global numerical simulations, focusing on two recent developments in particular. The first is quasi-cyclic magnetic activity in a long-duration dynamo simulation. Although mean fields comprise only a few percent of the total magnetic energy they exhibit remarkable order, with multiple polarity reversals and systematic variability on time scales of 6-15 years. The second development concerns the maintenance of the meridional circulation. Recent high-resolution simulations have captured the subtle nonlinear dynamical balances with more fidelity than previous, more laminar models, yielding more coherent circulation patterns. These patterns are dominated by a single cell in each hemisphere, with poleward and equatorward flow in the upper and lower convection zone respectively. We briefly address the implications of and future of these modeling efforts.

When our Sun was young it rotated much more rapidly than now. Observations of young, rapidly rotating stars indicate that many possess substantial magnetic activity and strong axisymmetric magnetic fields. We conduct simulations of dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code to explore the complex coupling between rotation, convection and magnetism. Here we study dynamo action realized in the bulk of the convection zone for two systems, rotating at three and five times the current solar rate. We find that substantial organized global-scale magnetic fields are achieved by dynamo action in these systems. Striking wreathes of magnetism are built in the midst of the convection zone, coexisting with the turbulent convection. This is a great surprise, for many solar dynamo theories have suggested that a tachocline of penetration and shear at the base of the convection zone is a crucial ingredient for organized dynamo action, whereas these simulations do not include such tachoclines. Some dynamos achieved in these rapidly rotating states build persistent global-scale fields which maintain amplitude and polarity for thousands of days. In the case at five times the solar rate, the dynamo can undergo cycles of activity, with fields varying in strength and even changing polarity. As the magnetic fields wax and wane in strength, the primary response in the convective flows involves the axisymmetric differential rotation, which begins to vary on similar time scales. Bands of relatively fast and slow fluid propagate toward the poles on time scales of roughly 500 days. In the Sun, similar patterns are observed in the poleward branch of the torsional oscillations, and these may represent a response to poleward propagating magnetic field deep below the solar surface.

In the outer envelope of the Sun and in other stars, differential rotation and meridional circulation are maintained via the redistribution of momentum and energy by convective motions. In order to properly capture such processes in a numerical model, the correct spherical geometry is essential. In this paper I review recent insights into the maintenance of mean flows in the solar interior obtained from high-resolution simulations of solar convection in rotating spherical shells. The Coriolis force induces a Reynolds stress which transports angular momentum equatorward and also yields latitudinal variations in the convective heat flux. Meridional circulations induced by baroclinicity and rotational shear further redistribute angular momentum and alter the mean stratification. This gives rise to a complex nonlinear interplay between turbulent convection, differential rotation, meridional circulation, and the mean specific entropy profile. I will describe how this drama plays out in our simulations as well as in solar and stellar convection zones.

We review recent advances in modeling global-scale convection and dynamo processes with the Anelastic Spherical Harmonic (ASH) code. In particular, we have recently achieved the first global-scale solar convection simulations that exhibit turbulent pumping of magnetic flux into a simulated tachocline and the subsequent organization and amplification of toroidal field structures by rotational shear. The presence of a tachocline not only promotes the generation of mean toroidal flux, but it also enhances and stabilizes the mean poloidal field throughout the convection zone, promoting dipolar structure with less frequent polarity reversals. The magnetic field generated by a convective dynamo with a tachocline and overshoot region is also more helical overall, with a sign reversal in the northern and southern hemispheres. Toroidal tachocline fields exhibit little indication of magnetic buoyancy instabilities but may be undergoing magneto-shear instabilities.

In the solar convection zone, rotation couples with intensely turbulent convection to drive a strong differential rotation and achieve complex magnetic dynamo action. Our sun must have rotated more rapidly in its past, as is suggested by observations of many rapidly rotating young solar-type stars. Here we explore the effects of more rapid rotation on the global-scale patterns of convection in such stars and the flows of differential rotation and meridional circulation which are self-consistently established. The convection in these systems is richly time dependent and in our most rapidly rotating suns a striking pattern of localized convection emerges. Convection near the equator in these systems is dominated by one or two nests in longitude of locally enhanced convection, with quiescent streaming flow in between at the highest rotation rates. These active nests of convection maintain a strong differential rotation despite their small size. The structure of differential rotation is similar in all of our more rapidly rotating suns, with fast equators and slower poles. We find that the total shear in differential rotation Delta Omega grows with more rapid rotation while the relative shear Delta Omega/Omega_0 decreases. In contrast, at more rapid rotation the meridional circulations decrease in energy and peak velocities and break into multiple cells of circulation in both radius and latitude.

This paper reports on efforts to simulate wave propagation in the solar interior. Presented is work on extending a numerical code for constant entropy acoustic waves in the absence of magnetic fields to the case where magnetic fields are present. A set of linearized magnetohydrodynamic (MHD) perturbation equations has been derived and implemented.

Stellar dynamos are driven by complex couplings between rotation and turbulent convection, which drive global-scale flows and build and rebuild stellar magnetic fields. When stars like our sun are young, they rotate much more rapidly than the current solar rate. Observations generally indicate that more rapid rotation is correlated with stronger magnetic activity and perhaps more effective dynamo action. Here we examine the effects of more rapid rotation on dynamo action in a star like our sun. We find that vigorous dynamo action is realized, with magnetic field generated throughout the bulk of the convection zone. These simulations do not possess a penetrative tachocline of shear where global-scale fields are thought to be organized in our sun, but despite this we find strikingly ordered fields, much like sea-snakes of toroidal field, which are organized on global scales. We believe this to be a novel finding.

In the solar convection zone, rotation couples with intensely turbulent convection to build global-scale flows of differential rotation and meridional circulation. Our sun must have rotated more rapidly in its past, as is suggested by observations of many rapidly rotating young solar-type stars. Here we explore the effects of more rapid rotation on the patterns of convection in such stars and the global-scale flows which are self-consistently established. The convection in these systems is richly time dependent and in our most rapidly rotating suns a striking pattern of spatially localized convection emerges. Convection near the equator in these systems is dominated by one or two patches of locally enhanced convection, with nearly quiescent streaming flow in between at the highest rotation rates. These active nests of convection maintain a strong differential rotation despite their small size. The structure of differential rotation is similar in all of our more rapidly rotating suns, with fast equators and slower poles. We find that the total shear in differential rotation, as measured by latitudinal angular velocity contrast, Delta_Omega, increases with more rapid rotation while the relative shear, Delta_Omega/Omega, decreases. In contrast, at more rapid rotation the meridional circulations decrease in both energy and peak velocities and break into multiple cells of circulation in both radius and latitude.

The global scales of solar convection are studied through three-dimensional simulations of compressible convection carried out in spherical shells of rotating fluid which extend from the base of the convection zone to within 15 Mm of the photosphere. Such modelling at the highest spatial resolution to date allows study of distinctly turbulent convection, revealing that coherent downflow structures associated with giant cells continue to play a significant role in maintaining the strong differential rotation that is achieved. These giant cells at lower latitudes exhibit prograde propagation relative to the mean zonal flow, or differential rotation, that they establish, and retrograde propagation of more isotropic structures with vortical character at mid and high latitudes. The interstices of the downflow networks often possess strong and compact cyclonic flows. The evolving giant-cell downflow systems can be partly masked by the intense smaller scales of convection driven closer to the surface, yet they are likely to be detectable with the helioseismic probing that is now becoming available. Indeed, the meandering streams and varying cellular subsurface flows revealed by helioseismology must be sampling contributions from the giant cells, yet it is difficult to separate out these signals from those attributed to the faster horizontal flows of supergranulation. To aid in such detection, we use our simulations to describe how the properties of giant cells may be expected to vary with depth, how their patterns evolve in time, and analyze the statistical features of correlations within these complex flow fields.

We present nonlinear three-dimensional simulations of the stably-stratified portion of the solar tachocline in which the rotational shear is maintained by mechanical forcing. When a broad toroidal field profile is specified as an initial condition, a clam-shell instability ensues which is similar to the freely-evolving cases studied previously. After the initial nonlinear saturation, the residual mean fields are apparently too weak to sustain the instability indefinitely. However, when a mean poloidal field is imposed in addition to the rotational shear, a statistically-steady state is achieved in which the clam-shell instability is operating continually. This state is characterized by a quasi-periodic exchange of energy between the mean toroidal field and the instability mode with a longitudinal wavenumber m=1. This quasi-periodic behavior has a timescale of several years and may have implications for tachocline dynamics and field emergence patterns throughout the solar activity cycle.

The operation of the solar global dynamo appears to involve many dynamical elements. Self-consistent MHD simulations which realistically incorporate all of these processes are not yet computationally feasible, though some elements can now be studied with reasonable fidelity. Here we consider the manner in which turbulent compressible convection within the bulk of the solar convection zone can generate large-scale magnetic fields through dynamo action. We accomplish this through a series of three-dimensional numerical simulations of MHD convection within rotating spherical shells using our ASH code on massively parallel supercomputers. Since differential rotation is a key ingredient in all dynamo models, we also examine here the nature of the rotation profiles that can be sustained within the deep convection zone as strong magnetic fields are built and maintained. We find that the convection is able to maintain a solar-like angular velocity profile despite the influence of Maxwell stresses which tend to oppose Reynolds stresses and thus reduce the latitudinal angular velocity contrast throughout the convection zone. The dynamo-generated magnetic fields exhibit a complex structure and evolution, with radial fields concentrated in downflow lanes and toroidal fields organized into twisted ribbons which are extended in longitude and which achieve field strengths of up to 5000 G. The flows and fields exhibit substantial kinetic and magnetic helicity although systematic hemispherical patterns are only apparent in the former. Fluctuating fields dominate the magnetic energy and account for most of the back-reaction on the flow via Lorentz forces. Mean fields are relatively weak and do not exhibit systematic latitudinal propagation or periodic polarity reversals as in the sun. This may be attributed to the absence of a tachocline.

We present the first results from three-dimensional spherical shell simulations of magnetic dynamo action realized by turbulent convection penetrating downward into a tachocline of rotational shear. This permits us to assess several dynamical elements believed to be crucial to the operation of the solar global dynamo, variously involving differential rotation resulting from convection, magnetic pumping, and amplification of fields by stretching within the tachocline. The simulations reveal that strong axisymmetric toroidal magnetic fields (about 3000 G in strength) are realized within the lower stable layer, unlike in the convection zone where fluctuating fields are predominant. The toroidal fields in the stable layer possess a striking persistent antisymmetric parity, with fields in the northern hemisphere largely of opposite polarity to those in the southern hemisphere. The associated mean poloidal magnetic fields there have a clear dipolar geometry, but we have not yet observed any distinctive reversals or latitudinal propagation. The presence of these deep magnetic fields appears to stabilize the sense of mean fields produced by vigorous dynamo action in the bulk of the convection zone.

The probability density functions (pdfs) of molecular line centroid velocity fluctuations and fluctuation differences at different spatial lags are estimated for several nearby molecular clouds with active internal star formation. The data consist of over 75,000 $^{13}$CO line profiles divided among twelve spatially and/or kinematically distinct regions. Although three regions (all in Mon R2) appear nearly Gaussian, the others show strong evidence for non-Gaussian, often nearly exponential, centroid velocity pdfs, possibly with power law contributions in the far tails. Evidence for nearly exponential centroid pdfs in the neutral HI component of the ISM is also presented, based on older optical and radio observations. These results are in contrast to pdfs found in isotropic incompressible turbulence experiments and simulations. Furthermore, no evidence is found for the scaling of difference pdf kurtosis with Reynolds number which is seen in incompressible turbulence, and the spatial distribution of high-amplitude velocity differences shows little indication of the filamentary appearance predicted by decay simulations dominated by vortical interactions. The variation with lag of the difference pdf moments is presented as a constraint on future simulations.

Probability density functions (pdfs) of $^{13}CO$ emission line centroid (line-of-sight, intensity-weighted average) velocities are presented for several densely sampled molecular clouds as quantitative descriptors of their underlying dynamics. Although some are approximately Gaussian in form, most of the pdfs exhibit relatively broader, often nearly exponential, tails, similar to the pdfs of velocity {\em differences} and {\em derivatives} (but not the velocity field itself) found in experiments and numerical simulations of incompressible turbulence. The broad pdf tails found in the present work are also similar to those found in decades-old measurements of interstellar velocity pdfs using atomic line centroids, and to the excess wing emission recently found in individual molecular line profiles. Some possible interpretations of the observed deviations are briefly discussed, although none of these account for the nearly exponential tails.