Baroclinic and barotropic instabilities in planetary atmospheres: energetics, equilibration and adjustment

NONLINEAR PROCESSES IN GEOPHYSICS 27:1 (2020) 147-173

Authors:

Peter Read, Daniel Kennedy, Neil Lewis, Helene Scolan, Fachreddin Tabataba-Vakili, Yixiong Wang, Susie Wright, Roland Young

Abstract:

© 2020 BMJ Publishing Group. All rights reserved. <p>Baroclinic and barotropic instabilities are well known as the mechanisms responsible for the production of the dominant energy-containing eddies in the atmospheres of Earth and several other planets, as well as Earth's oceans. Here we consider insights provided by both linear and nonlinear instability theories into the conditions under which such instabilities may occur, with reference to forced and dissipative flows obtainable in the laboratory, in simplified numerical atmospheric circulation models and in the planets of our solar system. The equilibration of such instabilities is also of great importance in understanding the structure and energetics of the observable circulation of atmospheres and oceans. Various ideas have been proposed concerning the ways in which baroclinic and barotropic instabilities grow to a large amplitude and saturate whilst also modifying their background flow and environment. This remains an area that continues to challenge theoreticians and observers, though some progress has been made. The notion that such instabilities may act under some conditions to adjust the background flow towards a critical state is explored here in the context of both laboratory systems and planetary atmospheres. Evidence for such adjustment processes is found relating to baroclinic instabilities under a range of conditions where the efficiency of eddy and zonal-mean heat transport may mutually compensate in maintaining a nearly invariant thermal structure in the zonal mean. In other systems, barotropic instabilities may efficiently mix potential vorticity to result in a flow configuration that is found to approach a marginally unstable state with respect to Arnol'd's second stability theorem. We discuss the implications of these findings and identify some outstanding open questions.</p>.

Thermal versus mechanical topography: an experimental investigation in a rotating baroclinic annulus

Geophysical and Astrophysical Fluid Dynamics Taylor and Francis 114:6 (2020) 763-797

Authors:

SD Marshall, Peter Read

Abstract:

We present a series of experimental investigations in which a differentially-heated annulus was used to investigate the effects of topography on rotating, stratified flows. In particular, we investigate blocking effects via azimuthally varying differential-heating and compare them to previous experiments utilising partial mechanical barriers. The thermal topography used consisted of a flat patch of heating elements covering a small azimuthal extent of the base, forming an equivalent of a partial barrier, to study the difference between blocked and unblocked flow. These azimuthally-varying heating experiments produced results with many similarities to our previous experiments with a mechanical barrier, despite the lack of a physical obstacle or formation of bottom-trapped waves. In particular, a unique flow structure was found when the drifting flow and the topography interacted in the form of an “interference” regime at low Taylor number, but forming an erratic “irregular” regime at higher Taylor number. This suggests that blocking may be induced by either or both of a thermal or mechanical inhomogeneity. Evidence of coherent/persistent resonant wave triads was noted in both kinds of experiment, though the component wavenumbers of the wave-triads and their impact on the flow was found to depend on the topography in question.

Non–adiabatic tidal oscillations induced by a planetary companion

Monthly Notices of the Royal Astronomical Society Oxford University Press (OUP) (2019)

Authors:

Andrew Bunting, John CB Papaloizou, Caroline Terquem

Abstract:

Abstract We calculate the dynamical tides raised by a close planetary companion on non–rotating stars of 1 M⊙ and 1.4 M⊙. Using the Henyey method, we solve the fully non–adiabatic equations throughout the star. The horizontal Lagrangian displacement is found to be 10 to 100 times larger than the equilibrium tide value in a thin region near the surface of the star. This is because non–adiabatic effects dominate in a region that extends from below the outer edge of the convection zone up to the stellar surface, and the equilibrium tide approximation is inconsistent with non–adiabaticity. Although this approximation generally applies in the low frequency limit, it also fails in the parts of the convection zone where the forcing frequency is small but larger than the Brunt-Väisälä frequency. We derive analytical estimates which give a good approximation to the numerical values of the magnitude of the ratio of the horizontal and radial displacements at the surface. The relative surface flux perturbation is also significant, on the order of 0.1% for a system modelled on 51 Pegasi b. Observations affected by the horizontal displacement may therefore be more achievable than previously thought, and brightness perturbations may be the result of flux perturbations rather than due to the radial displacement. We discuss the implication of this on the possibility of detecting such tidally excited oscillations, including the prospect of utilising the large horizontal motion for observations of systems such as 51 Pegasi.

Investigating the semiannual oscillation on Mars using data assimilation

Icarus Elsevier 333 (2019) 404-414 )

Authors:

Tao Ruan, Neil Lewis, S Lewis, Luca Montabone, Peter Read

Abstract:

A Martian semiannual oscillation (SAO), similar to that in the Earths tropical stratosphere, is evident in the Mars Analysis Correction Data Assimilation reanalysis dataset (MACDA) version 1.0, not only in the tropics, but also extending to higher latitudes. Unlike on Earth, the Martian SAO is found not always to reverse its zonal wind direction, but only manifests itself as a deceleration of the dominant wind at certain pressure levels and latitudes. Singular System Analysis (SSA) is further applied on the zonal-mean zonal wind in different latitude bands to reveal the characteristics of SAO phenomena at different latitudes. The second pair of principal components (PCs) is usually dominated by a SAO signal, though the SAO signal can be strong enough to manifest itself also in the first pair of PCs. An analysis of terms in the Transformed Eulerian Mean equation (TEM) is applied in the tropics to further elucidate the forcing processes driving the tendency of the zonal-mean zonal wind. The zonal-mean meridional advection is found to correlate strongly with the observed oscillations of zonal-mean zonal wind, and supplies the majority of the westward (retrograde) forcing in the SAO cycle. The forcing due to various non-zonal waves supplies forcing to the zonal-mean zonal wind that is nearly the opposite of the forcing due to meridional advection above ∼3 Pa altitude, but it also partly supports the SAO between 40 Pa and 3 Pa. Some distinctive features occurring during the period of the Mars year (MY) 25 global-scale dust storm (GDS) are also notable in our diagnostic results with substantially stronger values of eastward and westward momentum in the second half of MY 25 and stronger forcing due to vertical advection, transient waves and thermal tides

Overview of new MAST physics in anticipation of first results from MAST Upgrade

Nuclear Fusion IOP Science 59:11 (2019) 112011

Authors:

Harrison, RJ Akers, SY Allan, JS Allcock, JO Allen, L Appel, M Barnes, N Ben Ben Ayedl, W Boeglin, C Bowman, J Bradley, P Browning, P Bryant, M Carr, M Cecconello, CD Challis, S Chapman, IT Chapman, GJ Colyer, S Conroy, NJ Conway, M Cox, G Cunningham, RO Dendy, W Dorland

Abstract:

The mega amp spherical tokamak (MAST) was a low aspect ratio device (R/a  =  0.85/0.65 ~ 1.3) with similar poloidal cross-section to other medium-size tokamaks. The physics programme concentrates on addressing key physics issues for the operation of ITER, design of DEMO and future spherical tokamaks by utilising high resolution diagnostic measurements closely coupled with theory and modelling to significantly advance our understanding. An empirical scaling of the energy confinement time that favours higher power, lower collisionality devices is consistent with gyrokinetic modelling of electron scale turbulence. Measurements of ion scale turbulence with beam emission spectroscopy and gyrokinetic modelling in up-down symmetric plasmas find that the symmetry of the turbulence is broken by flow shear. Near the non-linear stability threshold, flow shear tilts the density fluctuation correlation function and skews the fluctuation amplitude distribution. Results from fast particle physics studies include the observation that sawteeth are found to redistribute passing and trapped fast particles injected from neutral beam injectors in equal measure, suggesting that resonances between the m  =  1 perturbation and the fast ion orbits may be playing a dominant role in the fast ion transport. Measured D–D fusion products from a neutron camera and a charged fusion product detector are 40% lower than predictions from TRANSP/NUBEAM, highlighting possible deficiencies in the guiding centre approximation. Modelling of fast ion losses in the presence of resonant magnetic perturbations (RMPs) can reproduce trends observed in experiments when the plasma response and charge-exchange losses are accounted for. Measurements with a neutral particle analyser during merging-compression start-up indicate the acceleration of ions and electrons. Transport at the plasma edge has been improved through reciprocating probe measurements that have characterised a geodesic acoustic mode at the edge of an ohmic L-mode plasma and particle-in-cell modelling has improved the interpretation of plasma potential estimates from ball-pen probes. The application of RMPs leads to a reduction in particle confinement in L-mode and H-mode and an increase in the core ionization source. The ejection of secondary filaments following type-I ELMs correlates with interactions with surfaces near the X-point. Simulations of the interaction between pairs of filaments in the scrape-off layer suggest this results in modest changes to their velocity, and in most cases can be treated as moving independently. A stochastic model of scrape-off layer profile formation based on the superposition of non-interacting filaments is in good agreement with measured time-average profiles. Transport in the divertor has been improved through fast camera imaging, indicating the presence of a quiescent region devoid of filament near the X-point, extending from the separatrix to ψ n ~ 1.02. Simulations of turbulent transport in the divertor show that the angle between the divertor leg on the curvature vector strongly influences transport into the private flux region via the interchange mechanism. Coherence imaging measurements show counter-streaming flows of impurities due to gas puffing increasing the pressure on field lines where the gas is ionised. MAST Upgrade is based on the original MAST device, with substantially improved capabilities to operate with a Super-X divertor to test extended divertor leg concepts. SOLPS-ITER modelling predicts the detachment threshold will be reduced by more than a factor of 2, in terms of upstream density, in the Super-X compared with a conventional configuration and that the radiation front movement is passively stabilised before it reaches the X-point. 1D fluid modelling reveals the key role of momentum and power loss mechanisms in governing detachment onset and evolution. Analytic modelling indicates that long legs placed at large major radius, or equivalently low at the target compared with the X-point are more amenable to external control. With MAST Upgrade experiments expected in 2019, a thorough characterisation of the sources of the intrinsic error field has been carried out and a mitigation strategy developed.