The solar dynamo begins near the surface

The magnetic dynamo cycle of the Sun features a distinct pattern: a propagating region of sunspot emergence appears around 30° latitude and vanishes near the equator every 11 years (ref. 1). Moreover, longitudinal flows called torsional oscillations closely shadow sunspot migration, undoubtedly sharing a common cause2. Contrary to theories suggesting deep origins of these phenomena, helioseismology pinpoints low-latitude torsional oscillations to the outer 5-10% of the Sun, the near-surface shear layer3,4. Within this zone, inwardly increasing differential rotation coupled with a poloidal magnetic field strongly implicates the magneto-rotational instability5,6, prominent in accretion-disk theory and observed in laboratory experiments7. Together, these two facts prompt the general question: whether the solar dynamo is possibly a near-surface instability. Here we report strong affirmative evidence in stark contrast to traditional models8 focusing on the deeper tachocline. Simple analytic estimates show that the near-surface magneto-rotational instability better explains the spatiotemporal scales of the torsional oscillations and inferred subsurface magnetic field amplitudes9. State-of-the-art numerical simulations corroborate these estimates and reproduce hemispherical magnetic current helicity laws10. The dynamo resulting from a well-understood near-surface phenomenon improves prospects for accurate predictions of full magnetic cycles and space weather, affecting the electromagnetic infrastructure of Earth.

Iterative methods for Navier-Stokes inverse problems

Even when the partial differential equation underlying a physical process can be evolved forward in time, the retrospective (backward in time) inverse problem often has its own challenges and applications. Direct adjoint looping (DAL) is the defacto approach for solving retrospective inverse problems, but it has not been applied to deterministic retrospective Navier-Stokes inverse problems in 2D or 3D. In this paper, we demonstrate that DAL is ill-suited for solving retrospective 2D Navier-Stokes inverse problems. Alongside DAL, we study two other iterative methods: simple backward integration (SBI) and the quasireversible method (QRM). As far as we know, our iterative SBI approach is novel, while iterative QRM has previously been used. Using these three iterative methods, we solve two retrospective inverse problems: 1D Korteweg-de Vries-Burgers (decaying nonlinear wave) and 2D Navier-Stokes (unstratified Kelvin-Helmholtz vortex). In both cases, SBI and QRM reproduce the target final states more accurately and in fewer iterations than DAL. We attribute this performance gap to additional terms present in SBI and QRM's respective backward integrations which are absent in DAL.

The photometric variability of massive stars due to gravity waves excited by core convection

Massive stars die in catastrophic explosions that seed the interstellar medium with heavy elements and produce neutron stars and black holes. Predictions of the explosion's character and the remnant mass depend on models of the star's evolutionary history. Models of massive star interiors can be empirically constrained by asteroseismic observations of gravity wave oscillations. Recent photometric observations reveal a ubiquitous red noise signal on massive main sequence stars; a hypothesized source of this noise is gravity waves driven by core convection. We present three-dimensional simulations of massive star convection extending from the star's centre to near its surface, with realistic stellar luminosities. Using these simulations, we predict the photometric variability due to convectively driven gravity waves at the surfaces of massive stars, and find that gravity waves produce photometric variability of a lower amplitude and lower characteristic frequency than the observed red noise. We infer that the photometric signal of gravity waves excited by core convection is below the noise limit of current observations, and thus the red noise must be generated by an alternative process.

Multimodal Excitation to Model the Quasibiennial Oscillation

The quasibiennial oscillation (QBO) of stratospheric winds is the most striking example of mean-flow generation and reversal by the nonlinear interactions of internal waves. Previous studies have used an idealized monochromatic forcing to investigate the QBO. Here we instead force a more realistic continuous wave spectrum. Unexpectedly, spreading the wave energy across a wide frequency range leads to more regular oscillations. We also find that different forcing spectra can yield the same QBO. Multimodal wave forcing is thus essential for understanding wave-mean-flow interactions in nature.

The magnetorotational instability prefers three dimensions

The magnetorotational instability (MRI) occurs when a weak magnetic field destabilizes a rotating, electrically conducting fluid with inwardly increasing angular velocity. The MRI is essential to astrophysical disc theory where the shear is typically Keplerian. Internal shear layers in stars may also be MRI-unstable, and they take a wide range of profiles, including near-critical. We show that the fastest growing modes of an ideal magnetofluid are three-dimensional provided the shear rate, S, is near the two-dimensional onset value, S _c . For a Keplerian shear, three-dimensional modes are unstable above S ≈ 0.10S _c , and dominate the two-dimensional modes until S ≈ 2.05S _c . These three-dimensional modes dominate for shear profiles relevant to stars and at magnetic Prandtl numbers relevant to liquid-metal laboratory experiments. Significant numbers of rapidly growing three-dimensional modes remainy well past 2.05S _c . These finding are significant in three ways. First, weakly nonlinear theory suggests that the MRI saturates by pushing the shear rate to its critical value. This can happen for systems, such as stars and laboratory experiments, that can rearrange their angular velocity profiles. Second, the non-normal character and large transient growth of MRI modes should be important whenever three-dimensionality exists. Finally, three-dimensional growth suggests direct dynamo action driven from the linear instability.

Connection between nonlinear energy optimization and instantons

How systems transit between different stable states under external perturbation is an important practical issue. We discuss here how a recently developed energy optimization method for identifying the minimal disturbance necessary to reach the basin boundary of a stable state is connected to the instanton trajectory from large deviation theory of noisy systems. In the context of the one-dimensional Swift-Hohenberg equation, which has multiple stable equilibria, we first show how the energy optimization method can be straightforwardly used to identify minimal disturbances-minimal seeds-for transition to specific attractors from the ground state. Then, after generalizing the technique to consider multiple, equally spaced-in-time perturbations, it is shown that the instanton trajectory is indeed the solution of the energy optimization method in the limit of infinitely many perturbations provided a specific norm is used to measure the set of discrete perturbations. Importantly, we find that the key features of the instanton can be captured by a low number of discrete perturbations (typically one perturbation per basin of attraction crossed). This suggests a promising new diagnostic for systems for which it may be impractical to calculate the instanton.

Order Out of Chaos: Slowly Reversing Mean Flows Emerge from Turbulently Generated Internal Waves

We demonstrate via direct numerical simulations that a periodic, oscillating mean flow spontaneously develops from turbulently generated internal waves. We consider a minimal physical model where the fluid self-organizes in a convective layer adjacent to a stably stratified one. Internal waves are excited by turbulent convective motions, then nonlinearly interact to produce a mean flow reversing on timescales much longer than the waves' period. Our results demonstrate for the first time that the three-scale dynamics due to convection, waves, and mean flow is generic and hence can occur in many astrophysical and geophysical fluids. We discuss efforts to reproduce the mean flow in reduced models, where the turbulence is bypassed. We demonstrate that wave intermittency, resulting from the chaotic nature of convection, plays a key role in the mean-flow dynamics, which thus cannot be captured using only second-order statistics of the turbulent motions.

Anomalous Chained Turbulence in Actively Driven Flows on Spheres

Recent experiments demonstrate the importance of substrate curvature for actively forced fluid dynamics. Yet, the covariant formulation and analysis of continuum models for nonequilibrium flows on curved surfaces still poses theoretical challenges. Here, we introduce and study a generalized covariant Navier-Stokes model for fluid flows driven by active stresses in nonplanar geometries. The analytical tractability of the theory is demonstrated through exact stationary solutions for the case of a spherical bubble geometry. Direct numerical simulations reveal a curvature-induced transition from a burst phase to an anomalous turbulent phase that differs distinctly from externally forced classical 2D Kolmogorov turbulence. This new type of active turbulence is characterized by the self-assembly of finite-size vortices into linked chains of antiferromagnetic order, which percolate through the entire fluid domain, forming an active dynamic network. The coherent motion of the vortex chain network provides an efficient mechanism for upward energy transfer from smaller to larger scales, presenting an alternative to the conventional energy cascade in classical 2D turbulence.

Numerical simulations of internal wave generation by convection in water

Water's density maximum at 4°C makes it well suited to study internal gravity wave excitation by convection: an increasing temperature profile is unstable to convection below 4°C, but stably stratified above 4°C. We present numerical simulations of a waterlike fluid near its density maximum in a two-dimensional domain. We successfully model the damping of waves in the simulations using linear theory, provided we do not take the weak damping limit typically used in the literature. To isolate the physical mechanism exciting internal waves, we use the spectral code dedalus to run several simplified model simulations of our more detailed simulation. We use data from the full simulation as source terms in two simplified models of internal-wave excitation by convection: bulk excitation by convective Reynolds stresses, and interface forcing via the mechanical oscillator effect. We find excellent agreement between the waves generated in the full simulation and the simplified simulation implementing the bulk excitation mechanism. The interface forcing simulations overexcite high-frequency waves because they assume the excitation is by the "impulsive" penetration of plumes, which spreads energy to high frequencies. However, we find that the real excitation is instead by the "sweeping" motion of plumes parallel to the interface. Our results imply that the bulk excitation mechanism is a very accurate heuristic for internal-wave generation by convection.