Thin-Shell Magnetohydrodynamic Equations for the Solar Tachocline

We present a new system of equations designed to study global-scale dynamics in the stably-stratified portion of the solar tachocline. This system is derived from the 3D equations of magnetohydrodynamics in a rotating spherical shell under the assumption that the shell is thin and stably-stratified (subadiabatic). The resulting thin-shell model can be regarded as a magnetic generalization of the hydrostatic primitive equations often used in meteorology. It is simpler in form than the more general anelastic or Boussinesq equations, making it more amenable to analysis and interpretation and more computationally efficient. However, the thin-shell system is still three-dimensional and as such represents an important extension to previous 2D and shallow-water approaches. In this paper we derive the governing equations for our thin-shell model and discuss its underlying assumptions, its context relative to other models, and its application to the solar tachocline. We also demonstrate that the dissipationless thin-shell system conserves energy, angular momentum and magnetic helicity.     

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

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^° to 40^° 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²⁰ Mx and 10²¹ Mx, and more simulations of the previously investigated case of 10²² 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 ≥?40 kG are good candidates for the progenitors of large (10²¹ Mx to 10²² Mx) solar active regions, which agrees with the results from Article 1 for flux tubes of 10²² 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.     

Rapid rotation,active nests of convection and global-scale flows in solar-likestars

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

Dynamo action in simulations of penetrative solar convection with an imposed tachocline

We summarize new and continuing three-dimensional spherical shell simulations of dynamo action by convection allowed to penetrate downward into a tachocline of rotational shear. The inclusion of an imposed tachocline allows us to examine several processes believed to be essential in the operation of the global solar dynamo, including differential rotation, magnetic pumping, and the stretching and organization of fields within the tachocline. In the stably stratified core, our simulations reveal that strong axisymmetric magnetic fields (of ∼ 3000 G strength) can be built, and that those fields generally exhibit a striking antisymmetric parity, with fields in the northern hemisphere largely of opposite polarity to those in the southern hemisphere. In the convection zone above, fluctuating fields dominate over weaker mean fields. New calculations indicate that the tendency toward toroidal fields of antisymmetric parity is relatively insensitive to initial magnetic field configurations; they also reveal that on decade-long timescales, the magnetic fields can briefly enter (and subsequently emerge from) states of symmetric parity.We have not yet observed any overall reversals of the field polarity, nor systematic latitudinal propagation. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Anelastic convection-driven dynamo benchmarks

Benchmark solutions for fully nonlinear anelastic compressible convection and dynamo action in a rotating spherical shell are proposed. Three benchmarks are specified. The first is a purely hydrodynamic case, which is steady in a uniformly drifting frame. The second is a self-excited saturated dynamo solution, also steady in a drifting frame. The third is again a self-excited dynamo but is unsteady in time, and it has a higher Rayleigh number than the steady dynamo benchmark. Four independent codes have been tested against these benchmarks, and very satisfactory agreement has been found. This provides an accurate reference standard against which new anelastic codes can be tested.     

Buoyant Magnetic Loops Generated by Global Convective Dynamo Action

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 higher 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.     

Scaling variables and interpretation of eigenvalues in principal component analysis of geologic data

The dominant feature distinguishing one method of principal components analysis from another is the manner in which the original data are transformed prior to the other computations. The only other distinguishing feature of any importance is whether the eigenvectors of the inner product-moment of the transformed data matrix are taken directly as the Q-mode scores or scaled by the square roots of their associated eigenvalues and called the R-mode loadings. If the eigenvectors are extracted from the product-moment correlation matrix, the variables, in effect, were transformed by column standardization (zero means and unit variances), and the sum of the p-largest eigenvalues divided by the sum of all the eigenvalues indicates the degree to which a model containing pcomponents will account for the total variance in the original data. However, if the data were transformed in any manner other than column standardization, the eigenvalues cannot be used in this manner, but can only be used to determine the degree to which the model will account for the transformed data. Regardless of the type of principal components analysis that is performed—even whether it is Ror Q-mode—the goodness-of-fit of the model to the original data is given better by the eigenvalues of the correlation matrix than by those of the matrix that was actually factored.