Academic literature on the topic 'Quasigeostrophic formulation'

Abstract This paper examines the sensitivity of the circulation of an idealized primitive equation two-level model on the form and strength of the heating, aiming to understand the qualitatively different sensitivity of the isentropic slope on differential heating reported by previous idealized studies when different model formulations are used. It is argued that this contrasting behavior might arise from differences in the internal determination of the heating. To test this contention, the two-level model is forced using two different heating formulations: a standard Newtonian cooling formulation and a highly simplified formulation in which the net lower-to-upper troposphere heat transport is prescribed by construction. The results are interpreted using quasigeostrophic turbulent closures, which have previously been shown to have predictive power for the model. It is found that the strength of the circulation, as measured by eddy length and velocity scales and by the strength of the energy cycle, scales with the vertical heating (the lower-to-upper troposphere heat transport), with a weak dependence. By contrast, the isentropic slope is only sensitive to the structure of the heating, as measured by the ratio between meridional versus vertical heating, and not to the actual strength of the heating. In general the heating is internally determined, and this ratio may either increase or decrease as the circulation strengthens. It is shown that the sign of the sensitivity depends on the steepness of the relation between vertical heating and stratification for the particular heating formulation used. The quasigeostrophic limit (fixed stratification) and the prescribed heating model constrain the possible range of behaviors and provide bounds of sensitivity for the model. These results may help explain the different sensitivity of the isentropic slope on differential heating for dry and moist models and for quasigeostrophic and primitive equation models.

AbstractA quasigeostrophic model is developed to diagnose the three-dimensional circulation, including the vertical velocity, in the upper ocean from high-resolution observations of sea surface height and buoyancy. The formulation for the adiabatic component departs from the classical surface quasigeostrophic framework considered before since it takes into account the stratification within the surface mixed layer that is usually much weaker than that in the ocean interior. To achieve this, the model approximates the ocean with two constant stratification layers: a finite-thickness surface layer (or the mixed layer) and an infinitely deep interior layer. It is shown that the leading-order adiabatic circulation is entirely determined if both the surface streamfunction and buoyancy anomalies are considered. The surface layer further includes a diabatic dynamical contribution. Parameterization of diabatic vertical velocities is based on their restoring impacts of the thermal wind balance that is perturbed by turbulent vertical mixing of momentum and buoyancy. The model skill in reproducing the three-dimensional circulation in the upper ocean from surface data is checked against the output of a high-resolution primitive equation numerical simulation.

Abstract Motivated by the success of potential vorticity (PV) thinking for Rossby waves and related shear flow phenomena, this work develops a buoyancy–vorticity formulation of gravity waves in stratified shear flow, for which the nonlocality enters in the same way as it does for barotropic/baroclinic shear flows. This formulation provides a time integration scheme that is analogous to the time integration of the quasigeostrophic equations with two, rather than one, prognostic equations, and a diagnostic equation for streamfunction through a vorticity inversion. The invertibility of vorticity allows the development of a gravity wave kernel view, which provides a mechanistic rationalization of many aspects of the linear dynamics of stratified shear flow. The resulting kernel formulation is similar to the Rossby-based one obtained for barotropic and baroclinic instability; however, since there are two independent variables—vorticity and buoyancy—there are also two independent kernels at each level. Though having two kernels complicates the picture, the kernels are constructed so that they do not interact with each other at a given level.

Abstract Some of the advances of the past century in our understanding of the general circulation of the atmosphere are described, starting with a brief summary of some of the key developments from the first half of the twentieth century, but with a primary focus on the period beginning with the midcentury breakthrough in baroclinic instability and quasigeostrophic dynamics. In addition to baroclinic instability, topics touched upon include the following: stationary wave theory, the role played by the two-layer model, scaling arguments for the eddy heat flux, the subtlety of large-scale eddy momentum fluxes, the Eliassen–Palm flux and the transformed Eulerian mean formulation, the structure of storm tracks, and the controls on the Hadley cell.

Abstract An eddy-resolving quasigeostrophic model of the Southern Ocean coupled to a dynamic atmospheric mixed layer is used to compare the performance of two different wind stress parameterization schemes. The first is the standard quadratic drag law, based on atmospheric velocity alone, whereas the second (more exact) formulation is based on the difference between ocean and atmosphere velocities. The two different schemes give very similar magnitudes of mean stress; however, the relative velocity scheme has substantially lower power input, resulting in a weaker eddy field, and consequently, greater circumpolar transport. These results are explored in terms of the existing theories of the Antarctic Circumpolar Current (including eddy saturation and eddy damping) and the implications for modeling the Southern Ocean are discussed.

Abstract The double-gyre circulation, formulated in terms of the quasigeostrophic equations, has a symmetry about the basin midlatitude (y → −y, ψ → −ψ), which is absent in a formulation based on the shallow-water equations. As a result, the shallow-water model does not have the pitchfork bifurcation structures that, in the case of the quasigeostrophic model, connect together multiple solution branches with elongated and contracted recirculation gyres. For the shallow-water model, solution branches with elongated recirculation gyres are disconnected, and a one-parameter bifurcation analysis is unable to detect their existence. The deeply penetrating jet solution branches do, however, continue to exist, and can be found using a bifurcation analysis couched in terms of two parameters. An effective pair of parameters is the viscosity and a parameter controlling the symmetry of the wind stress profile. A bifurcation analysis with these parameters reveals the existence of new solution branches that were not found in previous bifurcation analyses of the shallow-water model. The new solutions have a jet extension that penetrates farther eastward and that is more stable than the jet-up and jet-down solutions found in previous studies. Furthermore, the origin of the low-frequency variability at low viscosities is associated with a sequence of bifurcations originating from one of the new steady-state solution branches. In particular, the eigenmode analysis of the new branch reveals that a so-called gyre mode is at the origin of the model’s low-frequency variability at decadal time scales.

AbstractThe quasigeostrophic and the generalized omega equations are the most widely used methods to reconstruct vertical velocity w from in situ data. As observational networks with much higher spatial and temporal resolutions are being designed, the question arises of identifying the approximations and scales at which an accurate estimation of w through the omega equation can be achieved and what critical scales and observables are needed. In this paper we test different adiabatic omega reconstructions of w over several regions representative of main oceanic regimes of the global ocean in a fully eddy-resolving numerical simulation with a 1/60° horizontal resolution. We find that the best reconstructions are observed in conditions characterized by energetic turbulence and/or weak stratification where near-surface frontal processes are felt deep into the ocean interior. The quasigeostrophic omega equation gives satisfactory results for scales larger than ~10 km horizontally while the improvements using a generalized formulation are substantial only in conditions where frontal turbulent processes are important (providing improvements with satisfactory reconstruction skill down to ~5 km in scale). The main sources of uncertainties that could be identified are related to processes responsible for ocean thermal wind imbalance (TWI), which is particularly difficult to account for (especially in observation-based studies) and to the deep flow that is generally improperly accounted for in omega reconstructions through the bottom boundary condition. Nevertheless, the reconstruction of mesoscale vertical velocities may be sufficient to estimate vertical fluxes of oceanic properties in many cases of practical interest.

This thesis presents a numerical exploration of the dynamics governing rotating flow driven by a surface stress in the " sliced cylinder " model of Pedlosky & Greenspan (1967) and Beardsley (1969), and its close relative, the " sliced cone " model introduced by Griffiths & Veronis (1997). The sliced cylinder model simulates the barotropic wind-driven circulation in a circular basin with vertical sidewalls, using a depth gradient to mimic the effects of a gradient in Coriolis parameter. In the sliced cone the vertical sidewalls are replaced by an azimuthally uniform slope around the perimeter of the basin to simulate a continental slope. Since these models can be implemented in the laboratory, their dynamics can be explored by a complementary interplay of analysis and numerical and laboratory experiments. ¶ In this thesis a derivation is presented of a generalised quasigeostrophic formulation which is valid for linear and moderately nonlinear barotropic flows over large-amplitude topography on an f-plane, yet retains the simplicity and conservation properties of the standard quasigeostrophic vorticity equation (which is valid only for small depth variations). This formulation is implemented in a numerical model based on a code developed by Page (1982) and Becker & Page (1990). ¶ The accuracy of the formulation and its implementation are confirmed by detailed comparisons with the laboratory sliced cylinder and sliced cone results of Griffiths (Griffiths & Kiss, 1999) and Griffiths & Veronis (1997), respectively. The numerical model is then used to provide insight into the dynamics responsible for the observed laboratory flows. In the linear limit the numerical model reveals shortcomings in the sliced cone analysis by Griffiths & Veronis (1998) in the region where the slope and interior join, and shows that the potential vorticity is dissipated in an extended region at the bottom of the slope rather than a localised region at the east as suggested by Griffiths & Veronis (1997, 1998). Welander's thermal analogy (Welander, 1968) is used to explain the linear circulation pattern, and demonstrates that the broadly distributed potential vorticity dissipation is due to the closure of geostrophic contours in this geometry. ¶ The numerical results also provide insight into features of the flow at finite Rossby number. It is demonstrated that separation of the western boundary current in the sliced cylinder is closely associated with a " crisis " due to excessive potential vorticity dissipation in the viscous sublayer, rather than insufficient dissipation in the outer western boundary current as suggested by Holland & Lin (1975) and Pedlosky (1987). The stability boundaries in both models are refined using the numerical results, clarifying in particular the way in which the western boundary current instability in the sliced cone disappears at large Rossby and/or Ekman number. A flow regime is also revealed in the sliced cylinder in which the boundary current separates without reversed flow, consistent with the potential vorticity " crisis " mechanism. In addition the location of the stability boundary is determined as a function of the aspect ratio of the sliced cylinder, which demonstrates that the flow is stabilised in narrow basins such as those used by Beardsley (1969, 1972, 1973) and Becker & Page (1990) relative to the much wider basin used by Griffiths & Kiss (1999). ¶ Laboratory studies of the sliced cone by Griffiths & Veronis (1997) showed that the flow became unstable only under anticyclonic forcing. It is shown in this thesis that the contrast between flow under cyclonic and anticyclonic forcing is due to the combined effects of the relative vorticity and topography in determining the shape of the potential vorticity contours. The vorticity at the bottom of the sidewall smooths out the potential vorticity contours under cyclonic forcing, but distorts them into highly contorted shapes under anticyclonic forcing. In addition, the flow is dominated by inertial boundary layers under cyclonic forcing and by standing Rossby waves under anticyclonic forcing due to the differing flow direction relative to the direction of Rossby wave phase propagation. The changes to the potential vorticity structure under strong cyclonic forcing reduce the potential vorticity changes experienced by fluid columns, and the flow approaches a steady free inertial circulation. In contrast, the complexity of the flow structure under anticyclonic forcing results in strong potential vorticity changes and also leads to barotropic instability under strong forcing. ¶ The numerical results indicate that the instabilities in both models arise through supercritical Hopf bifurcations. The two types of instability observed by Griffiths & Veronis (1997) in the sliced cone are shown to be related to the western boundary current instability and " interior instability " identified by Meacham & Berloff (1997). The western boundary current instability is trapped at the western side of the interior because its northward phase speed exceeds that of the fastest interior Rossby wave with the same meridional wavenumber, as discussed by Ierley & Young (1991). ¶ Numerical experiments with different lateral boundary conditions are also undertaken. These show that the flow in the sliced cylinder is dramatically altered when the free-slip boundary condition is used instead of the no-slip condition, as expected from the work of Blandford (1971). There is no separated jet, because the flow cannot experience a potential vorticity " crisis " with this boundary condition, so the western boundary current overshoots and enters the interior from the east. In contrast, the flow in the sliced cone is identical whether no-slip, free-slip or super-slip boundary conditions are applied to the horizontal flow at the top of the sloping sidewall, except in the immediate vicinity of this region. This insensitivity results from the extremely strong topographic steering near the edge of the basin due to the vanishing depth, which demands a balance between wind forcing and Ekman pumping on the upper slope, regardless of the lateral boundary condition. The sensitivity to the lateral boundary condition is related to the importance of lateral friction in the global vorticity balance. The integrated vorticity must vanish under the no-slip condition, so in the sliced cylinder the overall vorticity budget is dominated by lateral viscosity and Ekman friction is negligible. Under the free-slip condition the Ekman friction assumes a dominant role in the dissipation, leading to a dramatic change in the flow structure. In contrast, the much larger depth variation in the sliced cone leads to a global vorticity balance in which Ekman friction is always dominant, regardless of the boundary condition.