Rossby wave instability and apparent phase speeds in large ocean basins

Author Posting. © American Meteorological Society, 2007. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 37 (2007): 1177-1191, doi:10.1175/jpo3054.1.The stability of baroclinic Rossby waves in large ocean basins is examined, and the quasigeostrophic (QG) results of LaCasce and Pedlosky are generalized. First, stability equations are derived for perturbations on large-scale waves, using the two-layer shallow-water system. These equations resemble the QG stability equations, except that they retain the variation of the internal deformation radius with latitude. The equations are solved numerically for different initial conditions through eigenmode calculations and time stepping. The fastest-growing eigenmodes are intensified at high latitudes, and the slower-growing modes are intensified at lower latitudes. All of the modes have meridional scales and growth times that are comparable to the deformation radius in the latitude range where the eigenmode is intensified. This is what one would expect if one had applied QG theory in latitude bands. The evolution of large-scale waves was then simulated using the Regional Ocean Modeling System primitive equation model. The results are consistent with the theoretical predictions, with deformation-scale perturbations growing at rates inversely proportional to the local deformation radius. The waves succumb to the perturbations at the mid- to high latitudes, but are able to cross the basin at low latitudes before doing so. Also, the barotropic waves produced by the instability propagate faster than the baroclinic long-wave speed, which may explain the discrepancy in speeds noted by Chelton and Schlax.PEI was supported by a postdoctoral grant from the Norwegian Research Council, JHL was supported under the Norwegian NOCLIM II program, and JP was partly supported by NSF OCE 0451086

Baroclinic vortices over a sloping bottom

Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 1996Nonlinear quasigeostrophic flows in two layers over a topographic slope are considered. The evolution depends on the size of two parameters which indicate the degree of nonlinearity at depth. The first measures the importance of relative vorticity advection and the second of stretching vorticity. Two types of isolated vortex are used to examine the parameter dependence. An initially barotropic vortex remains barotropic only when the first parameter is large, otherwise topographic waves dominate at depth. An Initially surface-trapped vortex larger than deformation scale is baroclinically unstable when the second is large, but is stabilized by the slope otherwise. Both parameters are also relevant to cascading geostrophic turbulence. If the stretching parameter is large, a "barotropic cascade" occurs at the deformation radius (Rhines, 1977) and the cascade "arrests" when the relative vorticity parameter is order unity. If small, layer coupling is hindered and the cascade is arrested at the deformation scale, with the flow dominated by isotropic surface vortices. In both cases, the distinction between vortices and waves is transparent when viewing potential vorticity. It is more difficult to identify waves and vortices from the streamfunction fields, because the waves are present in both layers.Funding for this research was provided by Office of Naval Research Coastal Science Code, grants N00014-92-J-1643 and N00014-92-J-1528

Baroclinic vortices over a sloping bottom

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences, 1996.Includes bibliographical references (p. 212-220).by Joseph H. LaCasce, Jr.Ph.D

Reconciling float-based and tracer-based estimates of lateral diffusivities

Lateral diffusivities are computed from synthetic particles and tracers advected by a velocity field derived from sea-surface height measurements from the South Pacific, in a region west of Drake Passage. Three different estimates are compared: (1) the tracer-based “effective diffusivity” of Nakamura (1996), (2) the growth of the second moment of a cloud of tracer and (3) the single- and two-particle Lagrangian diffusivities. The effective diffusivity measures the cross-stream component of eddy mixing, so this article focuses on the meridional diffusivities for the others, as the mean flow (the ACC) is zonally oriented in the region. After an initial transient of a few weeks, the effective diffusivity agrees well with the meridional diffusivity estimated both from the tracer cloud and from the particles. This proves that particle- and tracer-based estimates of eddy diffusivities are equivalent, despite recent claims to the contrary. Convergence among the three estimates requires that the Lagrangian diffusivities be estimated using their asymptotic values, not their maximum values. The former are generally much lower than the latter in the presence of a mean flow. Sampling the long-time asymptotic behavior of Lagrangian diffusivities requires very large numbers of floats in field campaigns. For example, it is shown that hundreds of floats would be necessary to estimate the vertical and horizontal variations in eddy diffusivity in a sector of the Pacific Southern Ocean

Southern Ocean biogenic blooms freezing-in Oligocene colder climates

AbstractCrossing a key atmospheric CO2 threshold triggered a fundamental global climate reorganisation ~34 million years ago (Ma) establishing permanent Antarctic ice sheets. Curiously, a more dramatic CO2 decline (~800–400 ppm by the Early Oligocene(~27 Ma)), postdates initial ice sheet expansion but the mechanisms driving this later, rapid drop in atmospheric carbon during the early Oligocene remains elusive and controversial. Here we use marine seismic reflection and borehole data to reveal an unprecedented accumulation of early Oligocene strata (up to 2.2 km thick over 1500 × 500 km) with a major biogenic component in the Australian Southern Ocean. High-resolution ocean simulations demonstrate that a tectonically-driven, one-off reorganisation of ocean currents, caused a unique period where current instability coincided with high nutrient input from the Antarctic continent. This unrepeated and short-lived environment favoured extreme bioproductivity and enhanced sediment burial. The size and rapid accumulation of this sediment package potentially holds ~1.067 × 1015 kg of the ‘missing carbon’ sequestered during the decline from an Eocene high CO2-world to a mid-Oligocene medium CO2-world, highlighting the exceptional role of the Southern Ocean in modulating long-term climate.</jats:p

Circulation and Stirring in the Southeast Pacific Ocean and the Scotia Sea Sectors of the Antarctic Circumpolar Current

The large-scale middepth circulation and eddy diffusivities in the southeast Pacific Ocean and Scotia Sea sectors between 110° and 45°W of the Antarctic Circumpolar Current (ACC) are described based on a subsurface quasi-isobaric RAFOS-float-based Lagrangian dataset. These RAFOS float data were collected during the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES). The mean flow, adjusted to a common 1400-m depth, shows the presence of jets in the time-averaged sense with speeds of 6 cm s⁻¹ in the southeast Pacific Ocean and upward of 13 cm s−¹ in the Scotia Sea. These jets appear to be locked to topography in the Scotia Sea but, aside from negotiating a seamount chain, are mostly free of local topographic constraints in the southeast Pacific Ocean. The eddy kinetic energy (EKE) is higher than the mean kinetic energy everywhere in the sampled domain by about 50%. The magnitude of the EKE increases drastically (by a factor of 2 or more) as the current crosses over the Hero and Shackleton fracture zones into the Scotia Sea. The meridional isopycnal stirring shows lateral and vertical variations with local eddy diffusivities as high as 2800 ± 600 m2 s⁻¹ at 700 m decreasing to 990 ± 200 m² s⁻¹ at 1800 m in the southeast Pacific Ocean. However, the cross-ACC diffusivity in the southeast Pacific Ocean is significantly lower, with values of 690 ± 150 and 1000 ± 200 m² s⁻¹ at shallow and deep levels, respectively, due to the action of jets. The cross-ACC diffusivity in the Scotia Sea is about 1200 ± 500 m² s⁻¹

Deep-Reaching Global Ocean Overturning Circulation Generated by Surface Buoyancy Forcing

In contrast with the atmosphere, which is heated from below by solar radiation, the ocean is both heated and cooled from above. To drive a deep-reaching overturning circulation in this context, it is generally assumed that either intense interior mixing by winds and internal tides, or wind-driven upwelling is required; in their absence, the circulation is thought to collapse to a shallow surface cell. We demonstrate, using a primitive equation model with an idealized domain and no wind forcing, that the surface temperature forcing can in fact drive an interhemispheric overturning provided that there is an open channel unblocked in the zonal direction, such as in the Southern Ocean. With this geometry, rotating horizontal convection, in combination with asymmetric surface cooling between the north and south, drives a deep-reaching two-cell overturning circulation. The resulting vertical mid-depth stratification closely resembles that of the real ocean, suggesting that wind-driven pumping is not necessary to produce a deep-reaching overturning circulation, and that buoyancy forcing plays a more important role than is usually assumed

Deep-Reaching Global Ocean Overturning Circulation Generated by Surface Buoyancy Forcing

In contrast with the atmosphere, which is heated from below by solar radiation, the ocean is both heated and cooled from above. To drive a deep-reaching overturning circulation in this context, it is generally assumed that either intense interior mixing by winds and internal tides, or wind-driven upwelling is required; in their absence, the circulation is thought to collapse to a shallow surface cell. We demonstrate, using a primitive equation model with an idealized domain and no wind forcing, that the surface temperature forcing can in fact drive an interhemispheric overturning provided that there is an open channel unblocked in the zonal direction, such as in the Southern Ocean. With this geometry, rotating horizontal convection, in combination with asymmetric surface cooling between the north and south, drives a deep-reaching two-cell overturning circulation. The resulting vertical mid-depth stratification closely resembles that of the real ocean, suggesting that wind-driven pumping is not necessary to produce a deep-reaching overturning circulation, and that buoyancy forcing plays a more important role than is usually assumed

The Richardson's Law in Large-Eddy Simulations of Boundary Layer flows

Relative dispersion in a neutrally stratified planetary boundary layer (PBL) is investigated by means of Large-Eddy Simulations (LES). Despite the small extension of the inertial range of scales in the simulated PBL, our Lagrangian statistics turns out to be compatible with the Richardson $t^3$ law for the average of square particle separation. This emerges from the application of nonstandard methods of analysis through which a precise measure of the Richardson constant was also possible. Its values is estimated as $C_2\sim 0.5$ in close agreement with recent experiments and three-dimensional direct numerical simulations.Comment: 15 LaTex pages, 4 PS figure

Oceanic three-dimensional Lagrangian Coherent Structures: A study of a mesoscale eddy in the Benguela ocean region

We study three dimensional oceanic Lagrangian Coherent Structures (LCSs) in the Benguela region, as obtained from an output of the ROMS model. To do that we first compute Finite-Size Lyapunov exponent (FSLE) fields in the region volume, characterizing mesoscale stirring and mixing. Average FSLE values show a general decreasing trend with depth, but there is a local maximum at about 100 m depth. LCSs are extracted as ridges of the calculated FSLE fields. They present a "curtain-like" geometry in which the strongest attracting and repelling structures appear as quasivertical surfaces. LCSs around a particular cyclonic eddy, pinched off from the upwelling front are also calculated. The LCSs are confirmed to provide pathways and barriers to transport in and out of the eddy