Showing 1–50 of 60 results for author: Miesch, M

Inflows towards Bipolar Magnetic Active Regions and Their Nonlinear Impact on a Three-Dimensional Babcock-Leighton Solar Dynamo Model

Authors: Kinfe Teweldebirhan, Mark Miesch, Sarah Gibson

Abstract: The changing magnetic fields of the Sun are generated and maintained by a solar dynamo, the exact nature of which remains an unsolved fundamental problem in solar physics. Our objective in this paper is to investigate the role and impact of converging flows toward Bipolar Magnetic Regions (BMR inflows) on the Sun's global solar dynamo. These flows are large scale physical phenomena that have been… ▽ More The changing magnetic fields of the Sun are generated and maintained by a solar dynamo, the exact nature of which remains an unsolved fundamental problem in solar physics. Our objective in this paper is to investigate the role and impact of converging flows toward Bipolar Magnetic Regions (BMR inflows) on the Sun's global solar dynamo. These flows are large scale physical phenomena that have been observed and so should be included in any comprehensive solar dynamo model. We have augmented the Surface flux Transport And Babcock LEighton (STABLE) dynamo model to study the nonlinear feedback effect of BMR inflows with magnitudes varying with surface magnetic fields. This fully 3D realistic dynamo model produces the sunspot butterfly diagram and allows a study of the relative roles of dynamo saturation mechanisms such as tilt angle quenching and BMR inflows. The results of our STABLE simulations show that magnetic field dependent BMR inflows significantly affect the evolution of the BMRs themselves and result in a reduced buildup of the global poloidal field due to local flux cancellation within the BMRs, to an extent that is sufficient to saturate the dynamo. As a consequence, for the first time, we have achieved fully 3D solar dynamo solutions in which BMR inflows alone regulate the amplitudes and periods of the magnetic cycles. △ Less

Submitted 1 October, 2023; originally announced October 2023.

Comments: 19 pages 11 figures

The Puzzling Structure of Solar Convection: Window into the Dynamo

Authors: Nicholas A. Featherstone, Kyle C. Augustson, Jonathan M. Aurnou, Catherine Blume, Benjamin P. Brown, Nicholas Brummell, Keaton J. Burns, Michael A. Calkins, Maria Camisassa, Mausumi Dikpati, Yuhong Fan, J. R. Fuentes, Gustavo Guerrero, Bradley W. Hindman, Keith Julien, Irina N. Kitiashvili, Lydia Korre, Daniel Lecoanet, Bhishek Manek, Loren Matilsky, Mark Miesch, Nicholas J. Nelson, Jeffrey S. Oishi, Whitney T. Powers, Matthias Rempel , et al. (3 additional authors not shown)

Abstract: The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective… ▽ More The operation of the solar dynamo, with all of its remarkable spatio-temporal ordering, remains an outstanding problem of modern solar physics. A number of mechanisms that might plausibly contribute to its operation have been proposed, but the relative role played by each remains unclear. This uncertainty stems from continuing questions concerning the speed and structure of deep-seated convective flows. Those flows are in-turn thought to sustain both the Sun's turbulent EMF and the large-scale flows of differential rotation and meridional circulation suspected of influencing the dynamo's organization and timing. Continued progress in this area is complicated by (i) inconsistencies between helioseismic measurements of convective and meridional flow made with different techniques and instruments, and (ii) a lack of high-latitude data for convection, differential rotation, and meridional flow. We suggest that the path forward to resolving these difficulties is twofold. First, the acquisition of long-term helioseismic and emissivity measurements obtained from a polar vantage point is vital to complete our picture of the Sun's outer convection zone. Second, sustained and expanded investment in theory-oriented and combined theory/observational research initiatives will be crucial to fully exploit these new observations and to resolve inconsistencies between existing measurements. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Comments: White paper submitted to the Decadal Survey for Solar and Space Physics 2024-2033. 7 pages; 4 figures

Seismic Monitoring of the Sun's Far Hemisphere: A Crucial Component in Future Space Weather Forecasting (A White Paper Submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) -- SSPH 2024-2033)

Authors: Kiran Jain, C. Lindsey, E. Adamson, C. N. Arge, T. E. Berger, D. C. Braun, R. Chen, Y. M. Collado-Vega, M. Dikpati, T. Felipe, C. J. Henney, J. T. Hoeksema, R. W. Komm, K. D. Leka, A. R. Marble, V. Martinez Pillet, M. Miesch, L. J. Nickisch, A. A. Pevtsov, V. J. Pizzo, W. K. Tobiska, S. C. Tripathy, J. Zhao

Abstract: The purpose of this white paper is to put together a coherent vision for the role of helioseismic monitoring of magnetic activity in the Sun's far hemisphere that will contribute to improving space weather forecasting as well as fundamental research in the coming decade. Our goal fits into the broader context of helioseismology in solar research for any number of endeavors when helioseismic monito… ▽ More The purpose of this white paper is to put together a coherent vision for the role of helioseismic monitoring of magnetic activity in the Sun's far hemisphere that will contribute to improving space weather forecasting as well as fundamental research in the coming decade. Our goal fits into the broader context of helioseismology in solar research for any number of endeavors when helioseismic monitors may be the sole synoptic view of the Sun's far hemisphere. It is intended to foster a growing understanding of solar activity, as realistically monitored in both hemispheres, and its relationship to all known aspects of the near-Earth and terrestrial environment. Some of the questions and goals that can be fruitfully pursued through seismic monitoring of farside solar activity in the coming decade include: What is the relationship between helioseismic signatures and their associated magnetic configurations, and how is this relationship connected to the solar EUV irradiance over the period of a solar rotation?; How can helioseismic monitoring contribute to data-driven global magnetic-field models for precise space weather forecasting?; What can helioseismic monitors tell us about prospects of a flare, CME or high-speed stream that impacts the terrestrial environment over the period of a solar rotation?; How does the inclusion of farside information contribute to forecasts of interplanetary space weather and the environments to be encountered by human crews in interplanetary space? Thus, it is crucial for the development of farside monitoring of the Sun be continued into the next decade either through ground-based or space-borne observations. △ Less

Submitted 3 October, 2022; originally announced October 2022.

Principles Of Heliophysics: a textbook on the universal processes behind planetary habitability

Authors: Karel Schrijver, Fran Bagenal, Tim Bastian, Juerg Beer, Mario Bisi, Tom Bogdan, Steve Bougher, David Boteler, Dave Brain, Guy Brasseur, Don Brownlee, Paul Charbonneau, Ofer Cohen, Uli Christensen, Tom Crowley, Debrah Fischer, Terry Forbes, Tim Fuller-Rowell, Marina Galand, Joe Giacalone, George Gloeckler, Jack Gosling, Janet Green, Nick Gross, Steve Guetersloh , et al. (37 additional authors not shown)

Abstract: Heliophysics is the system science of the physical connections between the Sun and the solar system. As the physics of the local cosmos, it embraces space weather and planetary habitability. The wider view of comparative heliophysics forms a template for conditions in exoplanetary systems and provides a view over time of the aging Sun and its magnetic activity, of the heliosphere in different sett… ▽ More Heliophysics is the system science of the physical connections between the Sun and the solar system. As the physics of the local cosmos, it embraces space weather and planetary habitability. The wider view of comparative heliophysics forms a template for conditions in exoplanetary systems and provides a view over time of the aging Sun and its magnetic activity, of the heliosphere in different settings of the interstellar medium and subject to stellar impacts, of the space physics over evolving planetary dynamos, and of the long-term influence on planetary atmospheres by stellar radiation and wind. Based on a series of NASA-funded summer schools for early-career researchers, this textbook 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 book 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 could thrive. The text includes 200 "Activities" in the form of exercises, explorations, literature readings, "what if" challenges, and group discussion topics; many of the Activities provide additional information complementing the main text. Solutions and discussions are included in an Appendix for a selection of the exercises. △ Less

Submitted 27 March, 2024; v1 submitted 30 October, 2019; originally announced October 2019.

Comments: 429 pages, 121 figures, and 200 "activities" in the form of problems, exercises, literature readings, and "what if" challenges; V2.0: activities modified; added 'solutions' to selected activities; clarified text and corrected typos; added subject index; V2.1 added Activity hyperlinks and corrected typos (see 'Version history' in the manuscript)

Reconstructing Extreme Space Weather from Planet Hosting Stars

Authors: V. S. Airapetian, V. Adibekyan, M. Ansdell, D. Alexander, T. Bastian, S. Boro Saikia, A. S. Brun, O. Cohen, M. Cuntz, W. Danchi, J. Davenport, J. DeNolfo, R. DeVore, C. F. Dong, J. J. Drake, K. France, F. Fraschetti, K. Herbst, K. Garcia-Sage, M. Gillon, A. Glocer, J. L. Grenfell, G. Gronoff, N. Gopalswamy, M. Guedel , et al. (58 additional authors not shown)

Abstract: 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 conditi… ▽ More 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. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: White Paper submitted to the Astronomy and Astrophysics Decadal Survey (Astro2020), 8 pages, 1 figure

Starspot rotation rates vs. activity cycle phase: Butterfly diagrams of Kepler stars are unlike the Sun's

Authors: M. B. Nielsen, L. Gizon, R. H. Cameron, M. Miesch

Abstract: 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… ▽ More 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. △ Less

Submitted 16 December, 2018; originally announced December 2018.

Comments: Accepted for publication in Astronomy & Astrophysics. 11 pages, 8 figures

Journal ref: A&A 622, A85 (2019)

A Three- dimensional Babcock-Leighton Solar Dynamo Model with Non-axisymmetric Convective Flows

Authors: Gopal Hazra, Mark Miesch

Abstract: 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 t… ▽ More 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. △ Less

Submitted 18 September, 2018; originally announced September 2018.

Comments: 2 pages,1 figure, to appear in Proceedings of IAU Symposium No. 340, 2018

Long-term variability of the solar cycle in the Babcock-Leighton dynamo framework

Authors: Bidya Binay Karak, Mark Miesch

Abstract: 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… ▽ More 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. △ Less

Submitted 13 May, 2018; originally announced May 2018.

Comments: 4 pages including 4 figures

Incorporating Surface Convection into a 3D Babcock-Leighton Solar Dynamo Model

Authors: Gopal Hazra, Mark Miesch

Abstract: 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, re… ▽ More 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. △ Less

Submitted 9 April, 2018; originally announced April 2018.

Comments: 21 pages, 14 figures, Revised version is submitted to ApJ

Journal ref: The Astrophysical Journal, 864, 110 (2018)

Driving Solar Giant Cells through the Self-Organization of Near-Surface Plumes

Authors: Nicholas J. Nelson, Nicholas A. Featherstone, Mark S. Miesch, Juri Toomre

Abstract: 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 a… ▽ More 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. △ Less

Submitted 3 April, 2018; originally announced April 2018.

Comments: Accepted for publication in ApJ; 21 pages w/ 14 figures

Consequences of high effective Prandtl number on solar differential rotation and convective velocity

Authors: Bidya Binay Karak, Mark Miesch, Yuto Bekki

Abstract: 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 downf… ▽ More 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 ($κ$), 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-$κ$ simulations. The most interesting result of our low-$κ$ 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. △ Less

Submitted 29 April, 2018; v1 submitted 2 January, 2018; originally announced January 2018.

Comments: The published version may be freely available at https://doi.org/10.1063/1.5022034

Journal ref: Physics Of Fluids 30, 046602 (2018)

Recovery from Maunder-like Grand Minima in a Babcock--Leighton Solar Dynamo Model

Authors: Bidya Binay Karak, Mark Miesch

Abstract: 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 succe… ▽ More 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. △ Less

Submitted 20 June, 2018; v1 submitted 29 December, 2017; originally announced December 2017.

Comments: Published in ApJ; 6 pages including 4 figures

Journal ref: The Astrophysical Journal Letters, 860, L26, 2018

Solar Cycle Variability Induced by Tilt Angle Scatter in a Babcock--Leighton Solar Dynamo Model

Authors: Bidya Binay Karak, Mark Miesch

Abstract: 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 t… ▽ More 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 $σ_δ= 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 $σ_δ= 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 $σ_δ= 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. △ Less

Submitted 13 August, 2017; v1 submitted 27 June, 2017; originally announced June 2017.

Comments: Accepted for publication in ApJ

Journal ref: The Astrophysical Journal, 847, 69 (2017)

Meridional Circulation Dynamics in a Cyclic Convective Dynamo

Authors: D. Passos, M. Miesch, G. Guerrero, P. Charbonneau

Abstract: 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… ▽ More 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. △ Less

Submitted 18 July, 2017; v1 submitted 8 February, 2017; originally announced February 2017.

Comments: 20 pages, 13 figures, Accepted in A&A

Journal ref: A&A 607, A120 (2017)

A theoretical study of the build-up of the Sun's polar magnetic field by using a 3D kinematic dynamo model

Authors: Gopal Hazra, Arnab Rai Choudhuri, Mark S. Miesch

Abstract: 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 wit… ▽ More 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. △ Less

Submitted 25 January, 2017; v1 submitted 9 October, 2016; originally announced October 2016.

Comments: 18 pages, 18 figures, published in ApJ

Journal ref: The Astrophysical Journal, Volume 835, Issue 1, article id. 39, 16 pp. (2017)

Magnetic Helicity Reversals in a Cyclic Convective Dynamo

Authors: Mark S. Miesch, Mei Zhang, Kyle C. Augustson

Abstract: 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 ban… ▽ More 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. △ Less

Submitted 31 May, 2016; originally announced May 2016.

Comments: 6 pages, 4 color figures, accepted to Astrophysical Journal Letters

Velocity Amplitudes in Global Convection Simulations: The Role of the Prandtl Number and Near-Surface Driving

Authors: Bridget O'Mara, Mark S. Miesch, Nicholas A. Featherstone, K. C. Augustson

Abstract: 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 $κ$ along two paths in parameter space. If the kinematic viscosi… ▽ More 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 $κ$ along two paths in parameter space. If the kinematic viscosity $ν$ is decreased proportionally with $κ$ (fixing the Prandtl number $P_r = ν/κ$), we find that U increases but asymptotes toward a constant value, as found by Featherstone & Hindman (2016). However, if $ν$ is held fixed while decreasing $κ$ (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 $ν$) 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. △ Less

Submitted 19 March, 2016; originally announced March 2016.

Comments: 34 Pages, 8 Figures, submitted to a special issue of "Advances in Space Research" on "Solar Dynamo Frontiers"

Convection in Oblate Solar-Type Stars

Authors: Junfeng Wang, Mark S. Miesch, Chunlei Liang

Abstract: 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\% larg… ▽ More 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, $ΔΩ$, 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. △ Less

Submitted 6 July, 2016; v1 submitted 16 March, 2016; originally announced March 2016.

A Three-Dimensional Babcock-Leighton Solar Dynamo Model: Initial Results with Axisymmetric Flows

Authors: Mark S. Miesch, Kinfe Teweldebirhan

Abstract: 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… ▽ More 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). △ Less

Submitted 19 March, 2016; v1 submitted 11 November, 2015; originally announced November 2015.

Comments: 41 pages, 13 figures. Accepted for publication in a special issue of Advances in Space Research on "Solar Dynamo Frontiers"

The Wisconsin Plasma Astrophysics Laboratory

Authors: C. B. Forest, K. Flanagan, M. Brookhart, C. M. Cooper, M. Clark, V. Desangles, J. Egedal, D. Endrizzi, M. Miesch, I. V. Khalzov, H. Li, J. Milhone, M. Nornberg, J. Olson, E. Peterson, F. Roesler, A. Schekochihin, O. Schmitz, R. Siller, A. Spitkovsky, A. Stemo, J. Wallace, D. Weisberg, E. Zweibel

Abstract: 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… ▽ More 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. △ Less

Submitted 4 August, 2015; v1 submitted 23 June, 2015; originally announced June 2015.

Comments: 21 pages, 12 figures, 2 tables

Journal ref: J. Plasma Phys. 81 (2015) 345810501

Large-Eddy Simulations of Magnetohydrodynamic Turbulence in Heliophysics and Astrophysics

Authors: Mark S. Miesch, William H. Matthaeus, Axel Brandenburg, Arakel Petrosyan, Annick Pouquet, Claude Cambon, Frank Jenko, Dmitri Uzdensky, James Stone, Steve Tobias, Juri Toomre, Marco Velli

Abstract: 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 i… ▽ More 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. △ Less

Submitted 18 July, 2015; v1 submitted 7 May, 2015; originally announced May 2015.

Comments: 51 pages, 6 figures (Figs 2, 3, and 4 color), accepted to Space Science Reviews (in press). The paper is a product of a workshop on "LES Modeling in MHD Turbulence" held in Boulder, CO in May, 2013, sponsored by the Geophysical Turbulence Program at the National Center for Atmospheric Research

Solar Dynamics, Rotation, Convection and Overshoot

Authors: S. Hanasoge, M. S. Miesch, M. Roth, J. Schou, M. Schuessler, M. J. Thompson

Abstract: 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… ▽ More 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. △ Less

Submitted 30 March, 2015; originally announced March 2015.

Comments: 31 pages, 10 figures, Space Science Reviews

Convective Dynamo Simulation with a Grand Minimum

Authors: Kyle Augustson, Sacha Brun, Mark Miesch, Juri Toomre

Abstract: 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… ▽ More 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. △ Less

Submitted 13 March, 2015; originally announced March 2015.

Comments: Cool Stars 18 Conference Proceedings (2014), 15 Pages, 8 Figures, http://www2.lowell.edu/workshops/coolstars18/proceedings.html

A Compressible High-Order Unstructured Spectral Difference Code for Stratified Convection in Rotating Spherical Shells

Authors: Junfeng Wang, Chunlei Liang, Mark S. Miesch

Abstract: 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… ▽ More 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. △ Less

Submitted 6 March, 2015; v1 submitted 3 March, 2015; originally announced March 2015.

Meridional circulation dynamics from 3D MHD global simulations of solar convection

Authors: Dario Passos, Paul Charbonneau, Mark Miesch

Abstract: 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 meridiona… ▽ More 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. △ Less

Submitted 4 February, 2015; originally announced February 2015.

Comments: Accepted in ApJ Letters

Journal ref: ApJ 800 L18, 2015

Meridional Circulation in Solar and Stellar Convection Zones

Authors: Nicholas A. Featherstone, Mark S. Miesch

Abstract: 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… ▽ More 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. △ Less

Submitted 26 January, 2015; originally announced January 2015.

Comments: 23 pages, 14 Figures, Submitted to the Astrophysical Journal

Grand Minima and Equatorward Propagation in a Cycling Stellar Convective Dynamo

Authors: Kyle Augustson, Allan Sacha Brun, Mark Miesch, Juri Toomre

Abstract: 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… ▽ More 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 $α$-$Ω$ type. The time scales that appear to be relevant to the magnetic polarity reversal are also identified. △ Less

Submitted 10 July, 2015; v1 submitted 23 October, 2014; originally announced October 2014.

Comments: 24 pages, 13 figures, ApJ, Accepted

Differential rotation in main-sequence solar-like stars: Qualitative inference from asteroseismic data

Authors: Mikkel N. Lund, Mark S. Miesch, Jørgen Christensen-Dalsgaard

Abstract: 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 mode… ▽ More 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. △ Less

Submitted 28 August, 2014; v1 submitted 30 June, 2014; originally announced June 2014.

Comments: 30 pages, 21 figures, Published in The Astrophysical Journal (Figures 3 and 4 have been changed relative to the published version)

Journal ref: ApJ, 790, 121 (2014)

A 3D Babcock-Leighton Solar Dynamo Model

Authors: Mark S. Miesch, Mausumi Dikpati

Abstract: 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 prescri… ▽ More 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. △ Less

Submitted 25 January, 2014; originally announced January 2014.

Comments: 5 pages, 3 figures (all color). The Figures are low quality to minimize memory; for higher-quality figures contact the authors or obtain the final published version, when available. Submitted to the Astrophysical Journal Letters (Jan. 22, 2014)

Generating buoyant magnetic flux ropes in solar-like convective dynamos

Authors: Nicholas J. Nelson, Mark S. Miesch

Abstract: 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… ▽ More 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 $Ω$-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. △ Less

Submitted 8 January, 2014; originally announced January 2014.

Comments: 16 pages, 6 figures, accepted for publication in Plasma Physics and Controlled Fusion special issue "Self-Organization in Magnetic Flux Ropes"

Report number: LA-UR-13-28610

Cycling Dynamo in a Young Sun: Grand Minima and Equatorward Propagation

Authors: Kyle Augustson, Allan Sacha Brun, Mark Steven Miesch, Juri Toomre

Abstract: 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 zon… ▽ More 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. △ Less

Submitted 31 October, 2013; originally announced October 2013.

Comments: Submitted to ApJL 6 pages, 4 figures

A Theory on the Convective Origins of Active Longitudes on Solar-like Stars

Authors: Maria A. Weber, Yuhong Fan, Mark S. Miesch

Abstract: 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… ▽ More 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. △ Less

Submitted 8 May, 2013; originally announced May 2013.

Comments: This paper was accepted to The Astrophysical Journal on May 6, 2013

Hemispheric Asymmetries of Solar Photospheric Magnetism: Radiative, Particulate, and Heliospheric Impacts

Authors: Scott W. McIntosh, Robert J. Leamon, Joseph B. Gurman, Jean-Philippe Olive, Jonathan W. Cirtain, David H. Hathaway, Joan Burkepile, Mark Miesch, Robert S. Markel, Leonard Sitongia

Abstract: 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… ▽ More 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. △ Less

Submitted 30 January, 2013; originally announced February 2013.

Comments: Accepted to appear in the ApJ - 44 pages in referee format

Buoyant Magnetic Loops Generated by Global Convective Dynamo Action

Authors: Nicholas J. Nelson, Benjamin P. Brown, A. Sacha Brun, Mark S. Miesch, Juri Toomre

Abstract: 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 tim… ▽ More 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. △ Less

Submitted 21 December, 2012; originally announced December 2012.

Comments: Solar Physics format, 22 pages, accepted by Solar Physics for topical issue related NSO Workshop #26: Solar Origins of Space Weather and Space Climate

Magnetic Wreaths and Cycles in Convective Dynamos

Authors: Nicholas J. Nelson, Benjamin P. Brown, A. Sacha Brun, Mark S. Miesch, Juri Toomre

Abstract: 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 ou… ▽ More 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$Ω_\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. △ Less

Submitted 13 November, 2012; originally announced November 2012.

Comments: 21 pages, 16 figures, accepted for publication in ApJ

Comparing Simulations of Rising Flux Tubes Through the Solar Convection Zone with Observations of Solar Active Regions: Constraining the Dynamo Field Strength

Authors: Maria A. Weber, Yuhong Fan, Mark S. Miesch

Abstract: 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 articl… ▽ More 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. △ Less

Submitted 6 August, 2012; originally announced August 2012.

Comments: Accepted to Solar Physics Topical Issue: Solar Dynamics and Magnetism from the Interior to the Atmosphere

On the Amplitude of Convective Velocities in the Deep Solar Interior

Authors: Mark S. Miesch, Nicholas A. Featherstone, Matthias Rempel, Regner Trampedach

Abstract: 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 iso… ▽ More 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. △ Less

Submitted 8 August, 2012; v1 submitted 7 May, 2012; originally announced May 2012.

Comments: 16 pages, 4 figures, accepted to the Astrophysical Journal

Convective Babcock-Leighton Dynamo Models

Authors: Mark S. Miesch, Benjamin P. Brown

Abstract: 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 obser… ▽ More 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. △ Less

Submitted 12 January, 2012; originally announced January 2012.

Comments: 5 pages, 4 figures, submitted to the Astrophysical Journal Letters

Journal ref: The Astrophysical Journal Letters, 746:L26 (5pp) 2012 February 20

The Rise of Active Region Flux Tubes in the Turbulent Solar Convective Envelope

Authors: Maria A. Weber, Yuhong Fan, Mark S. Miesch

Abstract: 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 t… ▽ More 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. △ Less

Submitted 1 September, 2011; originally announced September 2011.

Comments: Accepted for publication in The Astrophysical Journal

Journal ref: 2011, ApJ, 741, 11

Buoyant Magnetic Loops in a Global Dynamo Simulation of a Young Sun

Authors: Nicholas J. Nelson, Benjamin P. Brown, Allan Sacha Brun, Mark S. Miesch, Juri Toomre

Abstract: 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… ▽ More 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. △ Less

Submitted 23 August, 2011; originally announced August 2011.

Comments: 6 pages, 3 figures, accepted in ApJ Letters

Gyroscopic Pumping in the Solar Near-Surface Shear Layer

Authors: Mark S. Miesch, Bradley W. Hindman

Abstract: 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 b… ▽ More 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. △ Less

Submitted 26 August, 2011; v1 submitted 21 June, 2011; originally announced June 2011.

Comments: 26 pages, 11 Figures, Accepted on August 25, 2011 for publication in the Astrophysical Journal

Journal ref: The Astrophysical Journal, 743: 79 (25pp), 2011 Dec 10

Magnetic Cycles in a Convective Dynamo Simulation of a Young Solar-type Star

Authors: Benjamin P. Brown, Mark S. Miesch, Matthew K. Browning, Allen Sacha Brun, Juri Toomre

Abstract: 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… ▽ More 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 $Ω$-effect, while the mean poloidal fields are built by turbulent correlations which are not necessarily well represented by a simple $α$-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] △ Less

Submitted 9 February, 2011; originally announced February 2011.

Comments: 20 pages, 14 figures, emulateapj format; accepted for publication in ApJ. Expanded and published version of sections 5-6 from http://arxiv.org/abs/0906.2407

Journal ref: ApJ, 731 (2011), 69 (19pp)

Global-scale Magnetism (and Cycles) in Dynamo Simulations of Stellar Convection Zones

Authors: Benjamin P Brown, Matthew K Browning, Allan Sacha Brun, Mark S Miesch, Juri Toomre

Abstract: 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… ▽ More 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. △ Less

Submitted 30 December, 2010; originally announced January 2011.

Comments: 8 pages, 3 figures; to appear in proceedings of Cool Stars 16, Seattle Aug 29 - Sep 2, 2010

Persistent Magnetic Wreaths in a Rapidly Rotating Sun

Authors: Benjamin P. Brown, Matthew K. Browning, Allen Sacha Brun, Mark S. Miesch, Juri Toomre

Abstract: 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 magneti… ▽ More 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 $α$-effect explains the regeneration of poloidal field. △ Less

Submitted 12 November, 2010; originally announced November 2010.

Comments: 17 pages, 9 figures, 1 appendix, emulateapj format; published version of sections 3-4, 7 and appendix from arXiv:0906.2407

Journal ref: 2010, ApJ, 711, 424-438

Global-scale wreath-building dynamos in stellar convection zones

Authors: Benjamin P Brown, Matthew K Browning, Allan S Brun, Mark S Miesch, Juri Toomre

Abstract: 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… ▽ More 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. △ Less

Submitted 1 November, 2010; originally announced November 2010.

Comments: 8 pages, 4 figures. To appear in IAU 271: "Astrophysical Dynamics: from Stars to Galaxies"

Global magnetic cycles in rapidly rotating younger suns

Authors: Nicholas J. Nelson, Benjamin P. Brown, Matthew K. Browning, Allan Sacha Brun, Mark S. Miesch, Juri Toomre

Abstract: 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 M… ▽ More 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. △ Less

Submitted 28 October, 2010; originally announced October 2010.

Comments: 4 pages, 3 figures, from IAU 273: The Physics of Sun and Star Spots

Magnetic Cycles and Meridional Circulation in Global Models of Solar Convection

Authors: Mark S. Miesch, Benjamin P. Brown, Matthew K. Browning, Allan Sacha Brun, Juri Toomre

Abstract: 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 syst… ▽ More 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. △ Less

Submitted 2 October, 2010; v1 submitted 30 September, 2010; originally announced September 2010.

Comments: 8 pages, 2 color figures. To appear in Proc. IAU Symposium 271, "Astrophysical Dynamics, from Stars to Galaxies", ed. A.S. Brun, N.H. Brummell & M.S. Miesch, in press Second version: corrected an error in the caption of Figure 2

Wreathes of Magnetism in Rapidly Rotating Suns

Authors: Benjamin P. Brown, Matthew K. Browning, Mark S. Miesch, Allan Sacha Brun, Juri Toomre

Abstract: 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 magne… ▽ More 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. △ Less

Submitted 12 June, 2009; originally announced June 2009.

Comments: 22 pages, 16 color figures in emulateapj format. Images heavily degraded, but legible. Submitted to ApJ

Differential rotation and meridional circulation in global models of solar convection

Authors: M. S. Miesch

Abstract: 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 obtai… ▽ More 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. △ Less

Submitted 19 November, 2008; originally announced November 2008.

Comments: 5 pages, 2 figures (one color)

Journal ref: Astron.Nachr.328:998-1001,2007

Three-Dimensional Simulations of Solar and Stellar Dynamos: The Influence of a Tachocline

Authors: M. S. Miesch, M. K. Browning, A. S. Brun, J. Toomre, B. P. Brown

Abstract: 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 she… ▽ More 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. △ Less

Submitted 18 November, 2008; originally announced November 2008.

Comments: 14 pages, 5 color figures, to appear in Proc. GONG 2008/SOHO XXI Meeting on Solar-Stellar Dynamos as Revealed by Helio and Asteroseismology, held August 15-18, 2008, Boulder, CO, Astronomical Soc. Pac. Conf. Series, volume TBD