Baroclinic Eddies Research Articles

The energy cascade associated with daily variability of the North Atlantic Oscillation

An analysis of anomalies of the energy cascade associated with daily fluctuations of the North Atlantic Oscillation (NAO) is presented here. The analysis consists of a detailed version of the Lorenz energy cycle, which decomposes the energy flows into baroclinic and barotropic terms and into zonal mean and eddy components, and was applied to the 6‐hourly ERA‐Interim reanalysis for the period of 1979–2016. The obtained results show that the positive NAO phase is preceded by a significant increase of synoptic baroclinic eddy activity. The eddy available potential energy is converted into kinetic energy and transferred to barotropic synoptic eddies. Then, the kinetic energy is transferred upscale into the barotropic planetary waves, which produce the NAO pattern. Immediately after the positive NAO peak, there is a significant reduction of baroclinic activity, and no indication of an anomalous increase of the global‐mean baroclinic eddy activity in the 10‐day period after the peak was found. At longer lags, between 10 and 20 days, the kinetic energy of barotropic eddies shows positive anomalies consistent with the eddy–zonal mean flow feedback. We conclude that the synoptic baroclinic eddy activity forces the NAO variability, and no modulating role of the NAO in the global‐mean baroclinic eddy activity (storminess) occurs.

A heuristic dynamical model of the North Atlantic Oscillation with a Lorenz-type chaotic attractor

Despite observational evidence of a distinct regime behaviour in the low-frequency variability of the North Atlantic Oscillation (NAO), the theoretical support for the existence of separate NAO regimes has been limited. Although NAO-like regimes have been detected in numerical models of varying complexity, the role of forcings and feedbacks originated by diabatic processes in the specific geographical setting of the North Atlantic region has not been explicitly advocated as a source of regime behaviour for the NAO. In this study, we develop a minimal theoretical model of the interactions between zonal flow (and associated temperature gradient), planetary waves and surface heat fluxes in the North Atlantic, guided by results from observational diagnostics. Using re-analysis data, we show that interactions between an equivalent barotropic NAO anomaly and climatological stationary waves with a baroclinic structure and a larger meridional scale generate a positive feedback between the NAO-like wave and the strength of the zonal flow. On the other hand, an increased zonal wind generates stronger surface heat fluxes which damp both the zonal and the wave components of the NAO-like anomaly. Using these observational results to guide the choice of empirical parameters, we construct a 3-variable model which is formally equivalent to the Lorenz (J Atmos Sci 20:130–141, 1963) chaotic model for Rayleigh–Benard convection, and possesses two regimes originated by oscillations around two weakly unstable stationary states. From a physical point of view, the two regimes correspond to opposite phases of a near-resonant planetary wave which occur when the zonal flow is either below (negative NAO) or above (positive NAO) the value corresponding to zero phase speed. Finally, we expand the model to a five-variable system by splitting the zonal-mean wind into its barotropic (i.e. height independent) and thermal component, and including a variable representing the area-averaged amplitude of high-frequency baroclinic eddies. Following earlier studies, we assume that high-frequency eddy activity grows by extracting available potential energy from the zonal-mean temperature gradient and decays by feeding kinetic energy into the barotropic component of the zonal wind. The extended model still displays a chaotic, two-regime behaviour, with additional sub-seasonal variability driven by the energy exchanges associated with high-frequency eddy amplification and decay.

Restratification of Abyssal Mixing Layers by Submesoscale Baroclinic Eddies

AbstractFor small-scale turbulence to achieve water mass transformation and thus affect the large-scale overturning circulation, it must occur in stratified water. Observations show that abyssal turbulence is strongly enhanced in the bottom few hundred meters in regions with rough topography, and it is thought that these abyssal mixing layers are crucial for closing and shaping the overturning circulation. If it were left unopposed, however, bottom-intensified turbulence would mix away the observed mixing-layer stratification over the course of a few years. It is proposed here that the homogenizing tendency of mixing may be balanced by baroclinic restratification. It is shown that bottom-intensified mixing, if it occurs on a large-scale topographic slope such as a midocean ridge flank, not only erodes stratification but also tilts isopycnals in the bottom few hundred meters. This tilting of isopycnals generates a reservoir of potential energy that can be tapped into by submesoscale baroclinic eddies. The eddies slide dense water under light water and thus restratify the mixing layer, similar to what happens in the surface mixed layer. This restratification is shown to be effective enough to balance the homogenizing tendency of mixing and to maintain the observed mixing-layer stratification. This suggests that submesoscale baroclinic eddies may play a crucial role in providing the stratification mixing can act on, thus allowing sustained water mass transformation. Through their restratification of abyssal mixing layers, submesoscale eddies may therefore directly affect the strength and structure of the abyssal overturning circulation.

Baroclinic Adjustment and Dissipative Control of Storm Tracks

Abstract The steady-state response of a midlatitude storm track to large-scale extratropical thermal forcing and eddy friction is investigated in a dry general circulation model with a zonally symmetric forcing. A two-way equilibration is found between the relative responses of the mean baroclinicity and baroclinic eddy intensity, whereby mean baroclinicity responds more strongly to eddy friction whereas eddy intensity responds more strongly to the thermal forcing of baroclinicity. These seemingly counterintuitive responses are reconciled using the steady state of a predator–prey relationship between baroclinicity and eddy intensity. This relationship provides additional support for the well-studied mechanism of baroclinic adjustment in Earth’s atmosphere, as well as providing a new mechanism whereby eddy dissipation controls the large-scale thermal structure of a baroclinically unstable atmosphere. It is argued that these two mechanisms of baroclinic adjustment and dissipative control should be used in tandem when considering storm-track equilibration.

The Influence of Aerosol Absorption on the Extratropical Circulation

AbstractThis study examines how aerosol absorption affects the extratropical circulation by analyzing the response to a globally uniform increase in black carbon (BC) simulated with an atmospheric general circulation model forced by prescribed sea surface temperatures. The model includes aerosol direct and semidirect effects, but not indirect or cloud-absorption effects. BC-induced heating in the free troposphere stabilizes the midlatitude atmospheric column, which results in less energetic baroclinic eddies and thus reduced meridional energy transport at midlatitudes. Upper-tropospheric BC also decreases the meridional temperature gradient on the equatorward flank of the tropospheric jet and yields a weakening and poleward shift of the jet, while boundary layer BC has no significant influence on the large-scale circulation since most of the heating is diffused by turbulence in the boundary layer. The effectiveness of BC in altering circulation generally increases with height. Dry baroclinic eddy theories can explain most of the extratropical response to free-tropospheric BC. Specifically, the decrease in vertical eddy heat flux related to a more stable atmosphere is the main mechanism for reestablishing atmospheric energy balance in the presence of BC-induced heating. Similar temperature responses are found in a dry idealized model, which further confirms the dominant role of baroclinic eddies in driving the extratropical circulation changes. The strong atmospheric-only response to BC suggests that absorbing aerosols are capable of altering synoptic-scale weather patterns. Its height dependence highlights the importance of better constraining model-simulated aerosol vertical distributions with satellite and field measurements.

Note on the Rate of Restratification in the Baroclinic Spindown of Fronts

AbstractThis paper revisits how the restratifying buoyancy flux generated by baroclinic mixed layer instabilities depends on environmental conditions. The frontal spindown is shown to produce buoyancy fluxes that increase significantly beyond the previously proposed and widely used scaling (f is the Coriolis parameter, Λ is the geostrophic shear, and H is the mixed layer depth), irrespective of whether the initial front is broad or narrow. This increase occurs after the initial phase of the nonlinear evolution, when the baroclinic eddies grow in size and develop velocities significantly in excess of the scaling assumption V ~ ΛH. Implications for parameterizing the restratification caused by baroclinic mixed layer instabilities in coarse-resolution models are discussed.

Computation and Characterization of Local Subfilter-Scale Energy Transfers in Atmospheric Flows

Abstract Atmospheric motions are governed by turbulent motions associated to nontrivial energy transfers at small scales (direct cascade) and/or at large scales (inverse cascade). Although it is known that the two cascades coexist, energy fluxes have been previously investigated from the spectral point of view but not on their instantaneous spatial and local structure. Here, we compute local and instantaneous subfilter-scale energy transfers in two sets of reanalyses (NCEP–NCAR and ERA-Interim) in the troposphere and the lower stratosphere for the year 2005. The fluxes are mostly positive (toward subgrid scales) in the troposphere and negative in the stratosphere, reflecting the baroclinic and barotropic nature of the motions, respectively. The most intense positive energy fluxes are found in the troposphere and are associated with baroclinic eddies or tropical cyclones. The computation of such fluxes can be used to characterize the amount of energy lost or missing at the smallest scales in climate and weather models.

The sensitivity of solsticial pauses to atmospheric ice and dust in the MarsWRF General Circulation Model

Mars exhibits less atmospheric variability at the solstices than it does during periods nearer the equinoxes. Much of this variability in air temperature and dust activity is attributable to a significant decrease in eastward traveling transient wave amplitudes in the lower atmosphere near the solstice. Previous versions of the Mars Weather Research and Forecasting (MarsWRF) model using only dust radiative forcing have reproduced the nature but not the magnitude of this ‘solsticial pause’ in atmospheric variability. In this paper, we use a version of MarsWRF that includes a fully-interactive dust and water cycle to simulate winter solsticial pauses under a range of dust and water ice conditions. The upgraded model specifically includes a new hybrid binned/two-moment microphysics model that simulates dust, water ice, and cloud condensation nuclei. The scheme tracks mass and number density for the three particle types throughout the atmosphere and allows advection by resolved winds, mixing by unresolved processes, and sedimentation that depends on particle size and density. Ice and dust particles interact with radiation in the atmosphere using a Mie scattering parameterization that allows for variable particle size and composition. Heterogeneous nucleation and condensation use an adaptive bin size scheme to accurately track the particle size during condensation and sublimation processes. All microphysical processes in the model are calculated within the dynamical timesteps using stability-guaranteed implicit calculations with no sub-timestepping. The impact of the addition of water processes to the model was assessed by comparing simulations with only interactive dust (dry simulations) and ones with a fully-interactive dust and water cycle (wet simulations). In dry simulations with dust storms a solsticial pause occurs in the northern winter with a magnitude (or ‘depth’) that depends on the opacity of the southern summer dust storms. In wet simulations that include water ice and dust particles, deep solsticial pauses are found in both winter hemispheres. In all simulations that reproduce the solsticial pause, energy and instability analysis suggest that a decrease in baroclinic instability and increase in barotropic energy conversion occurs during the solsticial pause. In dry simulations the decrease in baroclinic instability is caused by increased dust opacity leading to increased thermal static stability. In wet simulations, additional opacity from local cap-edge ice clouds reduces the near surface wind shear and further inhibits baroclinic eddy growth. The wet simulations are in better agreement with observations and tend to support results from other models that include ice cloud radiative effects.

Evaluation of the Momentum Closure Schemes in MPAS-Ocean

In order to compare and evaluate the performances of the Laplacian viscosity closure, the biharmonic viscosity closure, and the Leith closure momentum schemes in the MPAS-Ocean model, a variety of physical quantities, such as the relative reference potential energy (RPE) change, the RPE time change rate (RPETCR), the grid Reynolds number, the root mean square (RMS) of kinetic energy, and the spectra of kinetic energy and enstrophy, are calculated on the basis of results of a 3D baroclinic periodic channel. Results indicate that: 1) The RPETCR demonstrates a saturation phenomenon in baroclinic eddy tests. The critical grid Reynolds number corresponding to RPETCR saturation differs between the three closures: the largest value is in the biharmonic viscosity closure, followed by that in the Laplacian viscosity closure, and that in the Leith closure is the smallest. 2) All three closures can effectively suppress spurious dianeutral mixing by reducing the grid Reynolds number under sub-saturation conditions of the RPETCR, but they can also damage certain physical processes. Generally, the damage to the rotation process is greater than that to the advection process. 3) The dissipation in the biharmonic viscosity closure is strongly dependent on scales. Most dissipation concentrates on small scales, and the energy of small-scale eddies is often transferred to large-scale kinetic energy. The viscous dissipation in the Laplacian viscosity closure is the strongest on various scales, followed by that in the Leith closure. Note that part of the small-scale kinetic energy is also transferred to large-scale kinetic energy in the Leith closure. 4) The characteristic length scale L and the dimensionless parameter Г in the Leith closure are inherently coupled. The RPETCR is inversely proportional to the product of Г and L. When the product of Г and L is constant, both the simulated RPETCR and the inhibition of spurious dianeutral mixing are the same in all tests using the Leith closure. The dissipative scale in the Leith closure depends on the parameter L, and the dissipative intensity depends on the parameter Г. 5) Although optimal results may not be achieved by using the optimal parameters obtained from the 2D barotropic model in the 3D baroclinic simulation, the total energies are dissipative in all three closures. Dissipation is the strongest in the biharmonic viscosity closure, followed by that in the Leith closure, and that in the Laplacian viscosity closure is the weakest. Mesoscale eddies develop the fastest in the biharmonic viscosity closure after the baroclinic adjustment process finishes, and the kinetic energy reaches its maximum, which is attributed to the smallest dissipation of enstrophy in the biharmonic viscosity closure. Mesoscale eddies develop the slowest, and the kinetic energy peak value is the smallest in the Laplacian viscosity closure. Results in the Leith closure are between that in the biharmonic viscosity closure and the Laplacian viscosity closure.

Northern Westerlies during the Last Glacial Maximum: Results from CMIP5 Simulations

Abstract Based upon simulations from phase 5 of the Coupled Model Intercomparison Project (CMIP5), the vertical and regional characteristics of the northern westerlies during the Last Glacial Maximum (LGM) are investigated in this study. At the Northern Hemispheric scale, all nine available models simulate a poleward shift of the 200-hPa jet, and eight models simulate an equatorward shift of the 850-hPa jet compared to the preindustrial period; these shifts are of approximately 2°–3° latitude for the arithmetic multimodel mean. The upper-tropospheric cooling in the tropics, possibly due to reduced latent heat release, is expected to account for the poleward shift of the 200-hPa jet through the thermal wind relationship. Changes in the midlatitude baroclinic instability in response to the amplified polar cooling are associated with the jet stream in the lower troposphere through anomalous eddy activity. In particular, the types of predominant baroclinic eddies are regionally dependent. The behavior of the 850-hPa jet over the North Pacific is steered by transient eddies and characterized by a southward displacement during the LGM. By contrast, the remarkable enhancement of the North Atlantic jet stream throughout the troposphere is associated with the notably increased stationary eddy momentum convergence, presumably due to the presence of the Laurentide ice sheet over North America. In comparison with the proxy records, although there is no observational evidence explicitly indicating changes of the upper-level northern westerlies, the simulated LGM 850-hPa westerly wind field is indirectly concordant with the reconstructed moisture conditions over the Mediterranean region and southwestern North America.

Baroclinic Instability in the Presence of Convection

AbstractBaroclinic mixed-layer instabilities have recently been recognized as an important source of submesoscale energy in deep winter mixed layers. While the focus has so far been on the balanced dynamics of these instabilities, they occur in and depend on an environment shaped by atmospherically forced small-scale turbulence. In this study, idealized numerical simulations are presented that allow the development of both baroclinic instability and convective small-scale turbulence, with simple control over the relative strength. If the convection is only weakly forced, baroclinic instability restratifies the layer and shuts off convection, as expected. With increased forcing, however, it is found that baroclinic instabilities are remarkably resilient to the presence of convection. Even if the instability is too weak to restratify the layer and shut off convection, the instability still grows in the convecting environment and generates baroclinic eddies and fronts. This suggests that despite the vigorous atmospherically forced small-scale turbulence in winter mixed layers, baroclinic instabilities can persistently grow, generate balanced submesoscale turbulence, and modify the bulk properties of the upper ocean.

Covariate-based stochastic parameterization of baroclinic ocean eddies

AbstractIn this study we investigate a covariate-based stochastic approach to parameterize unresolved turbulent processes within a standard model of the idealised, wind-driven ocean circulation. We focus on vertical instead of horizontal coarse-graining, such that we avoid the subtle difficulties of horizontal coarsegraining. The corresponding eddy forcing is uniquely defined and has a clear physical interpretation related to baroclinic instability.We propose to emulate the baroclinic eddy forcing by sampling from the conditional probability distribution functions of the eddy forcing obtained from the baroclinic reference model data. These conditional probability distribution functions are approximated here by sampling uniformly from discrete reference values. We analyze in detail the different performances of the stochastic parameterization dependent on whether the eddy forcing is conditioned on a suitable flow-dependent covariate or on a timelagged covariate or on both. The results demonstrate that our non-Gaussian, non-linear methodology is able to accurately reproduce the first four statistical moments and spatial/temporal correlations of the stream function, energetics, and enstrophy of the reference baroclinic model.

The Yermak Pass Branch: A Major Pathway for the Atlantic Water North of Svalbard?

AbstractAn upward‐looking Acoustic Doppler Current Profiler deployed from July 2007 to September 2008 in the Yermak Pass, north of Svalbard, gathered velocity data from 570 m up to 90 m at a location covered by sea ice 10 months out of 12. Barotropic diurnal and semidiurnal tides are the dominant signals in the velocity (more than 70% of the velocity variance). In winter, baroclinic eddies at periods between 5 and 15 days and pulses of 1–2 month periodicity are observed in the Atlantic Water layer and are associated with a shoaling of the pycnocline. Mercator‐Ocean global operational model with daily and 1/12° spatial resolution is shown to have skills in representing low‐frequency velocity variations (>1 month) in the West Spitsbergen Current and in the Yermak Pass. Model outputs suggest that the Yermak Pass Branch has had a robust winter pattern over the last 10 years, carrying on average 31% of the Atlantic Water volume transport of the West Spitsbergen Current (36% in autumn/winter). However, those figures have to be considered with caution as the model neither simulates tides nor fully resolves eddies and ignores residual mean currents that could be significant.

Attribution of horizontal and vertical contributions to spurious mixing in an Arbitrary Lagrangian–Eulerian ocean model

We examine the separate contributions to spurious mixing from horizontal and vertical processes in an ALE ocean model, MOM6, using reference potential energy (RPE). The RPE is a global diagnostic which changes only due to mixing between density classes. We extend this diagnostic to a sub-timestep timescale in order to individually separate contributions to spurious mixing through horizontal (tracer advection) and vertical (regridding/remapping) processes within the model. We both evaluate the overall spurious mixing in MOM6 against previously published output from other models (MOM5, MITGCM and MPAS-O), and investigate impacts on the components of spurious mixing in MOM6 across a suite of test cases: a lock exchange, internal wave propagation, and a baroclinically-unstable eddying channel.The split RPE diagnostic demonstrates that the spurious mixing in a lock exchange test case is dominated by horizontal tracer advection, due to the spatial variability in the velocity field. In contrast, the vertical component of spurious mixing dominates in an internal waves test case. MOM6 performs well in this test case owing to its quasi-Lagrangian implementation of ALE. Finally, the effects of model resolution are examined in a baroclinic eddies test case. In particular, the vertical component of spurious mixing dominates as horizontal resolution increases, an important consideration as global models evolve towards higher horizontal resolutions.

The Influence of Atmospheric Cloud Radiative Effects on the Large-Scale Stratospheric Circulation

AbstractPrevious studies have explored the influence of atmospheric cloud radiative effects (ACRE) on the tropospheric circulation. Here the authors explore the influence of ACRE on the stratospheric circulation. The response of the stratospheric circulation to ACRE is assessed by comparing simulations run with and without ACRE. The stratospheric circulation response to ACRE is reproducible in a range of different GCMs and can be interpreted in the context of both a dynamically driven and a radiatively driven component.The dynamic component is linked to ACRE-induced changes in the vertical and meridional fluxes of wave activity. The ACRE-induced changes in the vertical flux of wave activity into the stratosphere are consistent with the ACRE-induced changes in tropospheric baroclinicity and thus the amplitude of midlatitude baroclinic eddies. They account for a strengthening of the Brewer–Dobson circulation, a cooling of the tropical lower stratosphere, a weakening and warming of the polar vortex, a reduction of static stability near the tropical tropopause transition layer, and a shortening of the time scale of extratropical stratospheric variability. The ACRE-induced changes in the equatorward flux of wave activity in the low-latitude stratosphere account for a strengthening of the zonal wind in the subtropical lower to midstratosphere.The radiative component is linked to ACRE-induced changes in the flux of longwave radiation into the lower stratosphere. The changes in radiative fluxes lead to a cooling of the extratropical lower stratosphere, changes in the static stability and cloud fraction near the extratropical tropopause, and a shortening of the time scales of extratropical stratospheric variability.The results highlight a previously overlooked pathway through which tropospheric climate influences the stratosphere.

Topographic and Mixed Layer Submesoscale Currents in the Near-Surface Southwestern Tropical Pacific

AbstractThe distribution and strength of submesoscale (SM) surface layer fronts and filaments generated through mixed layer baroclinic energy conversion and submesoscale coherent vortices (SCVs) generated by topographic drag are analyzed in numerical simulations of the near-surface southwestern Pacific, north of 16°S. In the Coral Sea a strong seasonal cycle in the surface heat flux leads to a winter SM “soup” consisting of baroclinic mixed layer eddies (MLEs), fronts, and filaments similar to those seen in other regions farther away from the equator. However, a strong wind stress seasonal cycle, largely in sync with the surface heat flux cycle, is also a source of SM processes. SM restratification fluxes show distinctive signatures corresponding to both surface cooling and wind stress. The winter peak in SM activity in the Coral Sea is not in phase with the summer dominance of the mesoscale eddy kinetic energy in the region, implying that local surface layer forcing effects are more important for SM generation than the nonlocal eddy deformation field. In the topographically complex Solomon and Bismarck Seas, a combination of equatorial proximity and boundary drag generates SCVs with large-vorticity Rossby numbers (Ro ~ 10). River outflows in the Bismarck and Solomon Seas make a contribution to SM generation, although they are considerably weaker than the topographic effects. Mean to eddy kinetic energy conversions implicate barotropic instability in SM topographic wakes, with the strongest values seen north of the Vitiaz Strait along the coast of Papua New Guinea.

Vertical resolution of baroclinic modes in global ocean models

Improvements in the horizontal resolution of global ocean models, motivated by the horizontal resolution requirements for specific flow features, has advanced modelling capabilities into the dynamical regime dominated by mesoscale variability. In contrast, the choice of the vertical grid remains a subjective choice, and it is not clear that efforts to improve vertical resolution adequately support their horizontal counterparts. Indeed, considering that the bulk of the vertical ocean dynamics (including convection) are parameterized, it is not immediately obvious what the vertical grid is supposed to resolve. Here, we propose that the primary purpose of the vertical grid in a hydrostatic ocean model is to resolve the vertical structure of horizontal flows, rather than to resolve vertical motion. With this principle we construct vertical grids based on their abilities to represent baroclinic modal structures commensurate with the theoretical capabilities of a given horizontal grid. This approach is designed to ensure that the vertical grids of global ocean models complement (and, importantly, to not undermine) the resolution capabilities of the horizontal grid. We find that for z-coordinate global ocean models, at least 50 well-positioned vertical levels are required to resolve the first baroclinic mode, with an additional 25 levels per subsequent mode. High-resolution ocean-sea ice simulations are used to illustrate some of the dynamical enhancements gained by improving the vertical resolution of a 1/10° global ocean model. These enhancements include substantial increases in the sea surface height variance (∼30% increase south of 40°S), the barotropic and baroclinic eddy kinetic energies (up to 200% increase on and surrounding the Antarctic continental shelf and slopes), and the overturning streamfunction in potential density space (near-tripling of the Antarctic Bottom Water cell at 65°S).

State Dependence of Atmospheric Response to Extratropical North Pacific SST Anomalies

AbstractBy performing two sets of high-resolution atmospheric general circulation model (AGCM) experiments, the authors find that the atmospheric response to a sea surface temperature (SST) anomaly in the extratropical North Pacific is sensitive to decadal variations of the background SST on which the SST anomaly is superimposed. The response in the first set of experiments, in which the SST anomaly is superimposed on the observed daily SST of 1981–90, strongly differs from the response in the second experiment, in which the same SST anomaly is superimposed on the observed daily SST of 1991–2000. The atmospheric response over the North Pacific during 1981–90 is eddy mediated, equivalent barotropic, and concentrated in the east. In contrast, the atmospheric response during 1991–2000 is weaker and strongest in the west. The results are discussed in terms of Rossby wave dynamics, with the proposed primary wave source switching from baroclinic eddy vorticity forcing over the eastern North Pacific in 1981–90 to mean-flow divergence over the western North Pacific in 1991–2000. The wave source changes are linked to the decadal reduction of daily SST variability over the eastern North Pacific and strengthening of the Oyashio Extension front over the western North Pacific. Thus, both daily and frontal aspects of the background SST variability in determining the atmospheric response to extratropical North Pacific SST anomalies are emphasized by these AGCM experiments.

Barotropic and Baroclinic Eddy Feedbacks in the Midlatitude Jet Variability and Responses to Climate Change–Like Thermal Forcings

Abstract Studies have suggested that the persistence in the meridional vacillation of the midlatitude jet (i.e., annular mode time scale) in comprehensive climate models is related to the model biases in climatological jet latitude, with important implications for projections of future climates and midlatitude weather events. Through the use of the recently developed finite-amplitude wave activity formalism and feedback quantifying techniques, this paper has quantified the role of barotropic and baroclinic eddy feedbacks in annular mode time scales using an idealized dry atmospheric model. The eddy–mean flow interaction that characterizes the persistent anomalous state of the midlatitude jet depends on processes associated with the lower-tropospheric source of vertically propagating Rossby waves, baroclinic mechanisms, and processes associated with upper-tropospheric wave propagation and breaking, barotropic mechanisms. A variety of climate change–like thermal forcings are used to generate a range of meridional shifts in the midlatitude eddy-driven jet. The idealized model shows a reduction in annular mode time scale associated with an increase in jet latitude, similar to comprehensive climate models. This decrease in time scale can be attributed to a similar decrease in the strength of the barotropic eddy feedback, which, in the positive phase of the annular mode, is characterized by anomalous potential vorticity (PV) mixing on the equatorward flank of the climatological jet. The decrease in subtropical PV mixing is, in turn, associated with a stronger subtropical jet as the eddy-driven jet is more distant from the subtropics. These results highlight the importance of subtropical eddy–mean flow interactions for the persistence of an eddy-driven jet.

Quantifying spatial distribution of spurious mixing in ocean models

Numerical mixing is inevitable for ocean models due to tracer advection schemes. Until now, there is no robust way to identify the regions of spurious mixing in ocean models. We propose a new method to compute the spatial distribution of the spurious diapycnic mixing in an ocean model. This new method is an extension of available potential energy density method proposed by Winters and Barkan (2013). We test the new method in lock-exchange and baroclinic eddies test cases. We can quantify the amount and the location of numerical mixing. We find high-shear areas are the main regions which are susceptible to numerical truncation errors. We also test the new method to quantify the numerical mixing in different horizontal momentum closures. We conclude that Smagorinsky viscosity has less numerical mixing than the Leith viscosity using the same non-dimensional constant.