-
Inertial Waves in a Nonlinear Simulation of the Sun's Convection Zone and Radiative Interior
Authors:
Catherine C. Blume,
Bradley W. Hindman,
Loren I. Matilsky
Abstract:
Recent observations of Rossby waves and other more exotic forms of inertial oscillations in the Sun's convection zone have kindled the hope that such waves might be used as a seismic probe of the Sun's interior. Here we present a 3D numerical simulation in spherical geometry that models the Sun's convection zone and upper radiative interior. This model features a wide variety of inertial oscillati…
▽ More
Recent observations of Rossby waves and other more exotic forms of inertial oscillations in the Sun's convection zone have kindled the hope that such waves might be used as a seismic probe of the Sun's interior. Here we present a 3D numerical simulation in spherical geometry that models the Sun's convection zone and upper radiative interior. This model features a wide variety of inertial oscillations, including both sectoral and tesseral equatorial Rossby waves, retrograde mixed inertial modes, prograde thermal Rossby waves, the recently observed high-frequency retrograde (HFR) vorticity modes, and what may be latitudinal overtones of these HFR modes. With this model, we demonstrate that sectoral and tesseral Rossby waves are ubiquitous within the radiative interior as well as within the convection zone. We suggest that there are two different Rossby-wave families in this simulation that live in different wave cavities: one in the radiative interior and one in the convection zone. Finally, we suggest that many of the retrograde inertial waves that appear in the convection zone, including the HFR modes, are in fact all related, being latitudinal overtones that are mixed modes with the prograde thermal Rossby waves.
△ Less
Submitted 21 December, 2023;
originally announced December 2023.
-
Confinement of the Solar Tachocline by a Non-Axisymmetric Dynamo
Authors:
Loren I. Matilsky,
Nicholas H. Brummell,
Bradley W. Hindman,
Juri Toomre
Abstract:
We recently presented the first 3D numerical simulation of the solar interior for which tachocline confinement was achieved by a dynamo-generated magnetic field. In this followup study, we analyze the degree of confinement as the magnetic field strength changes (controlled by varying the magnetic Prandtl number) in a coupled radiative zone (RZ) and convection zone (CZ) system. We broadly find thre…
▽ More
We recently presented the first 3D numerical simulation of the solar interior for which tachocline confinement was achieved by a dynamo-generated magnetic field. In this followup study, we analyze the degree of confinement as the magnetic field strength changes (controlled by varying the magnetic Prandtl number) in a coupled radiative zone (RZ) and convection zone (CZ) system. We broadly find three solution regimes, corresponding to weak, medium, and strong dynamo magnetic field strengths. In the weak-field regime, the large-scale magnetic field is mostly axisymmetric with regular, periodic polarity reversals (reminiscent of the observed solar cycle), but fails to create a confined tachocline. In the strong-field regime, the large-scale field is mostly non-axisymmetric with irregular, quasi-periodic polarity reversals, and creates a confined tachocline. In the medium-field regime, the large-scale field resembles a strong-field dynamo for extended intervals, but intermittently weakens to allow temporary epochs of strong differential rotation. In all regimes, the amplitude of poloidal field strength in the RZ is very well explained by skin-depth arguments, wherein the oscillating field that gives rise to the skin depth (in the medium- and strong-field cases) is a non-axisymmetric field structure rotating with respect to the RZ. These new simulations reaffirm that tachocline confinement by the solar dynamo (the so-called fast magnetic confinement scenario) is possible, but suggest a new picture in which non-axisymmetric field components rotating with respect to the RZ play the primary role, instead of the regularly reversing axisymmetic field associated with the 22-year cycle.
△ Less
Submitted 8 January, 2024; v1 submitted 16 November, 2023;
originally announced November 2023.
-
Assessing the Observability of Deep Meridional Flow Cells in the Solar Interior
Authors:
J. R. Fuentes,
Bradley W. Hindman,
Junwei Zhao,
Catherine C. Blume,
Maria E. Camisassa,
Nicholas A. Featherstone,
Thomas Hartlep,
Lydia Korre,
Loren I. Matilsky
Abstract:
Meridional circulation regulates the Sun's interior dynamics and magnetism. While it is well accepted that meridional flows are poleward at the Sun's surface, helioseismic observations have yet to provide a definitive answer for the depth at which those flows return to the equator, or the number of circulation cells in depth. Here, we explore the observability of multiple circulation cells stacked…
▽ More
Meridional circulation regulates the Sun's interior dynamics and magnetism. While it is well accepted that meridional flows are poleward at the Sun's surface, helioseismic observations have yet to provide a definitive answer for the depth at which those flows return to the equator, or the number of circulation cells in depth. Here, we explore the observability of multiple circulation cells stacked in radius. Specifically, we examine the seismic signature of several meridional flow profiles by convolving time-distance averaging kernels with mean flows obtained from a suite of 3D hydrodynamic simulations. At mid and high latitudes, we find that weak flow structures in the deep convection zone can be obscured by signals from the much stronger surface flows. This contamination of 1--2 m s$^{-1}$ is caused by extended side lobes in the averaging kernels, which produce a spurious equatorward signal with flow speeds that are one order of magnitude stronger than the original flow speeds in the simulations. At low latitudes, the flows in the deep layers of the simulations are stronger ($> 2$ m s$^{-1}$) and multiple cells across the convection zone can produce a sufficiently strong signal to survive the convolution process. Now that meridional flows can be measured over two decades of data, the uncertainties arising from convective noise have fallen to a level where they are comparable in magnitude to the systematic biases caused by non-local features in the averaging kernels. Hence, these systematic errors are beginning to influence current helioseismic deductions and need broader consideration.
△ Less
Submitted 4 December, 2023; v1 submitted 14 August, 2023;
originally announced August 2023.
-
The stellar thermal wind as a consequence of oblateness
Authors:
Loren I. Matilsky
Abstract:
In many rotating fluids, the lowest-order force balance is between gravity, pressure, and rotational acceleration ('GPR' balance). Terrestrial GPR balance takes the form of geostrophy and hydrostasy, which together yield the terrestrial thermal wind equation. By contrast, stellar GPR balance is an oblateness equation, which determines the departures of the thermal variables from spherical symmetry…
▽ More
In many rotating fluids, the lowest-order force balance is between gravity, pressure, and rotational acceleration ('GPR' balance). Terrestrial GPR balance takes the form of geostrophy and hydrostasy, which together yield the terrestrial thermal wind equation. By contrast, stellar GPR balance is an oblateness equation, which determines the departures of the thermal variables from spherical symmetry; its curl yields the 'stellar thermal wind equation.' In this sense, the stellar thermal wind should be viewed not as a consequence of geostrophy, but of baroclinicity in the oblateness. Here we treat the full stellar oblateness, including the thermal wind, using pressure coordinates. We derive the generalised stellar thermal wind equation and identify the parameter regime for which it holds. In the case of the Sun, not considering the full oblateness has resulted in conflicting calculations of the theoretical aspherical temperature anomaly. We provide new calculation here and find that the baroclinic anomaly is ~3-60 times smaller than the barotropic anomaly. Thus, the anomaly from the thermal wind may not be measurable helioseismically; but if measurement were possible, this would potentially yield a new way to bracket the depth of the solar tachocline.
△ Less
Submitted 25 July, 2023; v1 submitted 18 July, 2023;
originally announced July 2023.
-
Confinement of the Solar Tachocline by Dynamo Action in the Radiative Interior
Authors:
Loren I. Matilsky,
Bradley W. Hindman,
Nicholas A. Featherstone,
Catherine C. Blume,
Juri Toomre
Abstract:
A major outstanding problem in solar physics is the confinement of the solar tachocline, the thin shear layer that separates nearly solid-body rotation in the radiative interior from strong differential rotation in the convection zone. Here, we present the first 3-D, global solar simulation in which a tachocline is confined by a self-excited dynamo. The non-axisymmetric magnetism is initially buil…
▽ More
A major outstanding problem in solar physics is the confinement of the solar tachocline, the thin shear layer that separates nearly solid-body rotation in the radiative interior from strong differential rotation in the convection zone. Here, we present the first 3-D, global solar simulation in which a tachocline is confined by a self-excited dynamo. The non-axisymmetric magnetism is initially built in the convection zone and then diffusively imprints downward. Additionally, the field is locally amplified throughout the radiative interior by vigorous horizontal motions that arise from equatorial Rossby waves and possibly shear instabilities. Our work thus challenges the long-held notion that the Sun's dynamo magnetic field is amplified only as deep as the tachocline and stored in a quiescent radiative interior.
△ Less
Submitted 27 September, 2022; v1 submitted 26 June, 2022;
originally announced June 2022.
-
Constraining Global Solar Models through Helioseismic Analysis
Authors:
Andrey M. Stejko,
Alexander G. Kosovichev,
Nicholas A. Featherstone,
Gustavo Guerrero,
Bradley W. Hindman,
Loren I. Matilsky,
Jörn Warnecke
Abstract:
Global hydrodynamic simulations of internal solar dynamics have focused on replicating the conditions for solar-like differential rotation and meridional circulation using the results of helioseismic inversions as a constraint. Inferences of meridional circulation, however, have provided controversial results showing the possibility of one, two, or multiple cells along the radius. To resolve this…
▽ More
Global hydrodynamic simulations of internal solar dynamics have focused on replicating the conditions for solar-like differential rotation and meridional circulation using the results of helioseismic inversions as a constraint. Inferences of meridional circulation, however, have provided controversial results showing the possibility of one, two, or multiple cells along the radius. To resolve this controversy and develop a more robust understanding of global flow regimes in the solar interior, we apply a "forward-modeling" approach to the analysis of helioseismic signatures of meridional circulation profiles obtained from numerical simulations. We employ the global acoustic modeling code GALE to simulate the propagation of acoustic waves through regimes of mean mass flows generated by global hydrodynamic and magnetohydrodynamic models: EULAG, the Pencil Code, and the Rayleigh code. These models are used to create synthetic dopplergram data products, used as inputs for local time-distance helioseismology techniques. Helioseismic travel-time signals from solutions obtained through global numerical simulations are compared directly with inferences from solar observations, in order to set additional constraints on global model parameters in a direct way. We show that even though these models are able to replicate solar-like differential rotation, the resulting rotationally-constrained convection develops a multi-cell global meridional circulation profile that is measurably inconsistent with local time-distance inferences of solar observations. However, we find that the development of rotationally-unconstrained convection close to the model surface is able to maintain solar-like differential rotation, while having a significant impact on the helioseismic travel-time signal, replicating solar observations within one standard deviation of the error due to noise.
△ Less
Submitted 11 April, 2022;
originally announced April 2022.
-
Building and Maintaining a Solar Tachocline through Convective Dynamo Action
Authors:
Loren I. Matilsky,
Juri Toomre
Abstract:
For more than thirty years, the dynamical maintenance of the thin solar tachocline has remained one of the central outstanding problems of stellar astrophysics. Three main theories have been developed to explain the tachocline's thinness, but so far none of them has been shown to work convincingly in the extreme parameter regime of the solar interior. Here, we present a rotating, 3D, spherical-she…
▽ More
For more than thirty years, the dynamical maintenance of the thin solar tachocline has remained one of the central outstanding problems of stellar astrophysics. Three main theories have been developed to explain the tachocline's thinness, but so far none of them has been shown to work convincingly in the extreme parameter regime of the solar interior. Here, we present a rotating, 3D, spherical-shell simulation of a combined solar-like convection zone and radiative zone that achieves a tachocline built and maintained by convective dynamo action. Because of numerical constraints, the dynamo prevents the viscous spread of the tachocline instead of the Eddington-Sweet-time-scale radiative spread believed to occur in the Sun. Nonetheless, our simulation supports the scenario of tachocline confinement via the cyclic solar dynamo, and is the first time one of the main confinement scenarios has been realized in a global, 3D, spherical-shell geometry including nonlinear fluid motions and a self-consistently generated dynamo.
△ Less
Submitted 11 May, 2021;
originally announced May 2021.
-
Revisiting the Sun's Strong Differential Rotation along Radial Lines
Authors:
Loren I. Matilsky,
Bradley W. Hindman,
Juri Toomre
Abstract:
Current state-of-the-art models of the solar convection zone consist of solutions to the Navier-Stokes equations in rotating, 3D spherical shells. Such models are highly sensitive to the choice of boundary conditions. Here, we present two suites of simulations differing only in their outer thermal boundary condition, which is either one of fixed-entropy or fixed-entropy-gradient. We find that the…
▽ More
Current state-of-the-art models of the solar convection zone consist of solutions to the Navier-Stokes equations in rotating, 3D spherical shells. Such models are highly sensitive to the choice of boundary conditions. Here, we present two suites of simulations differing only in their outer thermal boundary condition, which is either one of fixed-entropy or fixed-entropy-gradient. We find that the resulting differential rotation is markedly different between the two sets. The fixed-entropy-gradient simulations have strong differential rotation contrast and isocontours tilted along radial lines (in good agreement with the Sun's interior rotation revealed by helioseismology), whereas the fixed-entropy simulations have weaker contrast and contours tilted in the opposite sense. We examine in detail the force balances in our models and find that the poleward transport of heat by Busse columns drives a thermal wind responsible for the different rotation profiles. We conclude that the Sun's strong differential rotation along radial lines may result from the solar emissivity being invariant with latitude (which is similar to the fixed-entropy-gradient condition in our models) and the poleward transport of heat by Busse columns. In future work on convection in the solar context, we strongly advise modelers to use a fixed-gradient outer boundary condition.
△ Less
Submitted 16 June, 2020; v1 submitted 31 March, 2020;
originally announced April 2020.
-
Exploring the Influence of Density Contrast on Solar Near-Surface Shear
Authors:
Loren I. Matilsky,
Bradley W. Hindman,
Juri Toomre
Abstract:
The advent of helioseismology has determined in detail the average rotation rate of the Sun as a function of radius and latitude. These data immediately reveal two striking boundary layers of shear in the solar convection zone (CZ): a tachocline at the base, where the differential rotation of the CZ transitions to solid-body rotation in the radiative zone, and a 35-Mm-thick near-surface shear laye…
▽ More
The advent of helioseismology has determined in detail the average rotation rate of the Sun as a function of radius and latitude. These data immediately reveal two striking boundary layers of shear in the solar convection zone (CZ): a tachocline at the base, where the differential rotation of the CZ transitions to solid-body rotation in the radiative zone, and a 35-Mm-thick near-surface shear layer (NSSL) at the top, where the rotation rate slows by about 5% with increasing radius. Though asteroseismology cannot probe the differential rotation of distant stars to the same level of detail that helioseismology can achieve for the Sun, it is possible that many cool stars with outer convective envelopes possess similar differential rotation characteristics, including both a tachocline and a NSSL. Here we present the results of 3D global hydrodynamic simulations of spherical-shell convection for a Sun-like star at different levels of density contrast across the shell. The simulations with high stratification possess characteristics of near-surface shear, especially at low latitudes. We discuss in detail the dynamical balance of torques giving rise to the NSSL in our models and interpret what these balances imply for the real Sun. We further discuss the dynamical causes that may serve to wipe out near-surface shear at high latitudes, and conclude by offering some theories as to how this problem might be tackled in future work.
△ Less
Submitted 1 November, 2018;
originally announced November 2018.
-
The Role of Downflows in Establishing Solar Near-Surface Shear
Authors:
Loren I. Matilsky,
Bradley W. Hindman,
J. Toomre
Abstract:
The dynamical origins of the Sun's tachocline and near-surface shear layer (NSSL) are still not well understood. We have attempted to self-consistently reproduce a NSSL in numerical simulations of a solar-like convection zone by increasing the density contrast across rotating, 3D spherical shells. We explore the hypothesis that high density contrast leads to near-surface shear by creating a rotati…
▽ More
The dynamical origins of the Sun's tachocline and near-surface shear layer (NSSL) are still not well understood. We have attempted to self-consistently reproduce a NSSL in numerical simulations of a solar-like convection zone by increasing the density contrast across rotating, 3D spherical shells. We explore the hypothesis that high density contrast leads to near-surface shear by creating a rotationally unconstrained layer of fast flows near the outer surface. Although our high-contrast models do have near-surface shear, it is confined primarily to low latitudes (between $\pm15^\circ$). Two distinct types of flow structures maintain the shear dynamically: rotationally $\textit{constrained}$ Busse columns aligned with the rotation axis and fast, rotationally $\textit{unconstrained}$ downflow plumes that deplete angular momentum from the outer fluid layers. The plumes form at all latitudes, and in fact are more efficient at transporting angular momentum inward at high latitudes. The presence of Busse columns at low latitudes thus appears essential to creating near-surface shear in our models. We conclude that a solar-like NSSL is unobtainable from a rotationally unconstrained outer fluid layer alone. In numerical models, the shear is eliminated through the advection of angular momentum by the meridional circulation. Therefore, a detailed understanding how the solar meridional circulation is dynamically achieved will be necessary to elucidate the origin of the Sun's NSSL.
△ Less
Submitted 28 September, 2018;
originally announced October 2018.