ALBERT

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution August, 1977

Description: Stimulated by new evidence from both "in situ" oceanic observations and results from numerical modelling, a laboratory study of quasigeostrophic flow and turbulence in a rotating homogeneous fluid has been undertaken. Two dimensional turbulence driven by a uniform distribution of sources and sinks which oscillate in time, can be fairly well reproduced in this context. Inertial time scales are about ten times smaller than Ekman spinup time, and typical Reynolds numbers read 2000. The observations emphasize the spectral tendency of the energy containing eddies. The case of no topography is first discussed. In steadily forced turbulence, it is observed that the energy containing scale is significantly larger than the forcing scale. In the decaying stage the red cascade is observed and rates of interaction are measured. Theoretical arguments for both behaviors are presented; the former concerning the forced turbulence case is believed to be new. The forcing is next applied over various large scale topographies, modelling the geophysical beta effect. The polar beta plane geometry preserves the above spectral characteristics but at the same time introduces anisotropy into the flow pattern. A broad westward mean flow develops in the north and is surrounded by a belt of cyclones lying on its southward side. The calculated second-order Eulerian mean flows induced by steadily and uniformly forced Rossby waves in a long zonal channel, exhibit much of the same momentum distribution in the inertial regime. In contrast, the "sliced cylinder" geometry which possesses no closed geostrophic contours drastically modifies the above picture. Both mean flow production and a large scale tendency for the eddies are inhibited. The geographical distribution of the eddy intensities and scales is now wildly inhomogeneous. The second aspect of this work is a study of the interaction of Rossby waves with mean flows. A zonally traveling, forced wave is generated near the southern boundary of a polar beta plane. Due to energy radiation in the free interior and (or) potential vorticity mixing by the finite amplitude waves, a westward zonal flow develops. The effect of the mean flow upon the forced steady waves is to weaken the anticyclones and intensify the cyclones. Pressure time series reveal a growth of harmonics and general spectral broadening as the waves travel freely inwards, suggesting active nonlinear interactions. An experimental test of Rhines' (1977) potential vorticity mixing theory is also presented at free latitudes. The decay period when the driving is suppressed shows that a net transfer from the waves to the mean flow kinetic energy occurs. Connection with hydrodynamic stability theory is discussed. Interaction of Rossby waves with an externally generated westward mean flow allows one to make a controlled study of the critical layer problem. For small amplitude waves, the mean flow is accelerated in the entire region between the forcing and the critical latitude which acts as a wall for mean wave momentum. In nonlinear runs the steady profile of the westward flow indicates that an accelerating force is acting everywhere, revealing the increasing transmission of wave momentum through the critical layer. At the same time, pressure measurements near the critical point show considerable fine structure developing over a long time scale. The third part deals with steady isolated source-sink flows in the sliced cylinder geometry. The response of the fluid to a meridionally oriented steady dipole extends exclusively westward of the forcing. The viscously balanced solutions are discussed and relevance to oceanic abyssal circulation is emphasized. With strong driving, the combination of a cyclone to the north and an anticyclone to the south is absolutely stable although the reverse configuration is not. A connection with a certain class of free, steady, isolated, inertial solutions developed recently by Stern (1976) is made.

Description: The DGRST . (FRACE) and the Joint Program in Oceanography, Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution offered a fellowship for the first two years. The National Science Foundation under Grant OCE75-2l 674 and the Office of Naval Research under Contract N00014-74-C0262-NR-083-004 supported this study for the final two years.

Description: Author Posting. © American Meteorological Society, 2016. 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 Climate 29 (2016): 1545-1571, doi:10.1175/JCLI-D-15-0509.1.

Description: Three sediment records of sea surface temperature (SST) are analyzed that originate from distant locations in the North Atlantic, have centennial-to-multicentennial resolution, are based on the same reconstruction method and chronological assumptions, and span the past 15 000 yr. Using recursive least squares techniques, an estimate of the time-dependent North Atlantic SST field over the last 15 kyr is sought that is consistent with both the SST records and a surface ocean circulation model, given estimates of their respective error (co)variances. Under the authors’ assumptions about data and model errors, it is found that the 10°C mixed layer isotherm, which approximately traces the modern Subpolar Front, would have moved by ~15° of latitude southward (northward) in the eastern North Atlantic at the onset (termination) of the Younger Dryas cold interval (YD), a result significant at the level of two standard deviations in the isotherm position. In contrast, meridional movements of the isotherm in the Newfoundland basin are estimated to be small and not significant. Thus, the isotherm would have pivoted twice around a region southeast of the Grand Banks, with a southwest–northeast orientation during the warm intervals of the Bølling–Allerød and the Holocene and a more zonal orientation and southerly position during the cold interval of the YD. This study provides an assessment of the significance of similar previous inferences and illustrates the potential of recursive least squares in paleoceanography.

Description: OM acknowledges support from the U.S. National Science Foundation. CW acknowledges support from the European Research Council ERC Grant ACCLIMATE 339108.

Description: 2016-08-19

Description: Numerical simulations of a realistic ocean general circulation model forced by prescribed surface fluxes are used to study the origin and structure of intrinsic interdecadal variability of the ocean circulation. When eddy-induced turbulent diffusivities are low enough, spontaneous oscillations of the Atlantic meridional overturning circulation (AMOC) with periods O(20) yr and amplitude O(1) Sv (1 Sv ≡ 106 m3 s−1) emerge. The transition from the steady to the oscillatory regime is shown to be consistent with a supercritical Hopf bifurcation of the horizontal Peclet number. Adding atmospheric thermal damping is shown to have a very limited influence on the domain of existence of intrinsic variability. The spatial structure of the mode consists of a dipole of sea surface temperature (SST)/sea surface height (SSH) anomalies centered at about 50°N with stronger variance in the western part of the subpolar gyre, in agreement with the observed Atlantic multidecadal oscillation (AMO) signature in this region. Specific features include a westward propagation of temperature anomalies from the source region located on the western flank of the North Atlantic Current (NAC) and a one-quarter phase lag between surface and subsurface (800 m) temperature anomalies. Local linear stability calculations including viscous and diffusive effects confirm that the North Atlantic Current is baroclinically unstable on scales of O(1000) km with growth rates of O(1) yr−1. Both the spatial structure of the mode and the period agree in magnitude with in situ measurements in the North Atlantic, suggesting that this intrinsic ocean mode participates in the observed Atlantic bidecadal climate variability.

Description: The turbulent diapycnal mixing in the ocean is currently obtained from microstructure and finestructure measurements, dye experiments, and inverse models. This study presents a new method that infers the diapycnal mixing from low-resolution numerical calculations of the World Ocean whose temperatures and salinities are restored to the climatology. At the difference of robust general circulation ocean models, diapycnal diffusion is not prescribed but inferred. At steady state the buoyancy equation shows an equilibrium between the large-scale diapycnal advection and the restoring terms that take the place of the divergence of eddy buoyancy fluxes. The geography of the diapycnal flow reveals a strong regional variability of water mass transformations. Positive values of the diapycnal flow indicate an erosion of a deep-water mass and negative values indicate a creation. When the diapycnal flow is upward, a diffusion law can be fitted in the vertical and the diapycnal eddy diffusivity is obtained throughout the water column. The basin averages of diapycnal diffusivities are small in the first 1500 m [O(10−5) m2 s−1] and increase downward with bottom values of about 2.5 × 10−4 m2 s−1 in all ocean basins, with the exception of the Southern Ocean (50°–30°S), where they reach 12 × 10−4 m2 s−1. This study confirms the small diffusivity in the thermocline and the robustness of the higher canonical Munk’s value in the abyssal ocean. It indicates that the upward dianeutral transport in the Atlantic mostly takes place in the abyss and the upper ocean, supporting the quasi-adiabatic character of the middepth overturning.

Description: The time-mean Argo float displacements and the World Ocean Atlas 2009 temperature–salinity climatology are used to obtain the total, top to bottom, mass transports. Outside of an equatorial band, the total transports are the sum of the vertical integrals of geostrophic- and wind-driven Ekman currents. However, these transports are generally divergent, and to obtain a mass conserving circulation, a Poisson equation is solved for the streamfunction with Dirichlet boundary conditions at solid boundaries. The value of the streamfunction on islands is also part of the unknowns. This study presents and discusses an energetic circulation in three basins: the North Atlantic, the North Pacific, and the Southern Ocean. This global method leads to new estimations of the time-mean western Eulerian boundary current transports maxima of 97 Sverdrups (Sv; 1 Sv ≡ 106 m3 s−1) at 60°W for the Gulf Stream, 84 Sv at 157°E for the Kuroshio, 80 Sv for the Agulhas Current between 32° and 36°S, and finally 175 Sv for the Antarctic Circumpolar Current at Drake Passage. Although the large-scale structure and boundary of the interior gyres is well predicted by the Sverdrup relation, the transports derived from the wind stress curl are lower than the observed transports in the interior by roughly a factor of 2, suggesting an important contribution of the bottom torques. With additional Argo displacement data, the errors caused by the presence of remaining transient terms at the 1000-db reference level will continue to decrease, allowing this method to produce increasingly accurate results in the future.

Description: Three sediment records of sea surface temperature (SST) are analyzed that originate from distant locations in the North Atlantic, have centennial-to-multicentennial resolution, are based on the same reconstruction method and chronological assumptions, and span the past 15 000 yr. Using recursive least squares techniques, an estimate of the time-dependent North Atlantic SST field over the last 15 kyr is sought that is consistent with both the SST records and a surface ocean circulation model, given estimates of their respective error (co)variances. Under the authors’ assumptions about data and model errors, it is found that the 10°C mixed layer isotherm, which approximately traces the modern Subpolar Front, would have moved by ~15° of latitude southward (northward) in the eastern North Atlantic at the onset (termination) of the Younger Dryas cold interval (YD), a result significant at the level of two standard deviations in the isotherm position. In contrast, meridional movements of the isotherm in the Newfoundland basin are estimated to be small and not significant. Thus, the isotherm would have pivoted twice around a region southeast of the Grand Banks, with a southwest–northeast orientation during the warm intervals of the Bølling–Allerød and the Holocene and a more zonal orientation and southerly position during the cold interval of the YD. This study provides an assessment of the significance of similar previous inferences and illustrates the potential of recursive least squares in paleoceanography.

Description: The vertically integrated potential energy of an incompressible stratified fluid formulated in density coordinates can be simply written as a weighted vertical sum of the squares of the vertical displacements of density surfaces, a general expression valid for arbitrary displacements. The sum of this form of potential energy and kinetic energy is then a conserved quantity for the multilayer shallow water model. The formulation in density coordinates is a natural one to find the Lorenz reference state of available potential energy (APE). We describe the method to compute the APE of an ocean state and provide two applications. The first is the classical double-gyre, wind-driven circulation simulated by a shallow water model at high resolution. We show that the eddy kinetic and eddy potential energies are localized in regions of large gradients of mean APE. These large gradients surround an APE minimum found between the two gyres. The second is the time-mean World Ocean Circulation reconstructed from hydrography (World Ocean Atlas) and reference velocities at 1000 db from the Argo float program to obtain an absolute circulation. The total available potential energy exceeds the total mean kinetic energy of the World Ocean by three orders of magnitude, pointing out the very small Burger number of the circulation. The Gulf Stream, the Kuroshio, the Agulhas retroflection, and the confluence regions are four examples that confirm the shallow water model results that large gradients of mean available potential energy can be used as predictors for the presence of high eddy kinetic energy (obtained here from satellite altimetry).

Description: The last glacial period was punctuated by rapid climate shifts, known as Dansgaard–Oeschger events, with strong imprint in the North Atlantic sector suggesting that they were linked with the Atlantic meridional overturning circulation. Here an idealized single-hemisphere three-dimensional ocean–atmosphere–sea ice coupled model is used to explore the possible origin of the instability driving these abrupt events and to provide a plausible explanation for the relative stability of the Holocene. Focusing on the physics of noise-free millennial oscillations under steady external (solar) forcing, it was shown that cold climates become unstable, that is, exhibit abrupt millennial-scale transitions, for significantly lower freshwater fluxes than warm climates, in agreement with previous studies making use of zonally averaged coupled models. This fundamental difference is a direct consequence of the weaker stratification of the glacial ocean, mainly caused by upper-ocean cooling. Using a two-hemisphere configuration of a coupled climate model of intermediate complexity, it is shown that this result is robust to the added presence of a bottom water mass of southern origin. The analysis reveals that under particular conditions, a pronounced interdecadal variability develops during warm interstadials. While the nature of the instability driving the millennial oscillations is identical to that found in ocean models under mixed boundary conditions, these interstadial–interdecadal oscillations share the same characteristics as those previously found in ocean models forced by fixed surface fluxes. The wind stress forcing is shown to profoundly affect both the properties and bifurcation structure of thermohaline millennial oscillations across a wide range of variation of freshwater forcing. In particular, it is shown that the wind stress forcing favors the maintenance of thermally direct meridional overturning circulations during the cold stadial phases of Dansgaard–Oeschger cycles.

Description: The stationary atmospheric response to an idealized sea surface temperature anomaly (SSTA) is studied with a quasigeostrophic atmospheric model of the Southern Hemisphere. Sensitivity of the stationary response to the midlatitude SSTA location is determined and responses are decomposed on vertical modes. The SSTA almost directly forces baroclinic responses, inducing warm-air anomalies 40°–50° downstream, eastward, to the SSTA. These baroclinic responses arise from an equilibrium between the SSTA-induced anomalous vortex stretching and (i) advection by the quasi-stationary flow and (ii) dissipation by high-frequency eddies. The barotropic response consists of a midlatitude ridge (trough) and a South Pole trough (ridge) for SSTAs localized from the Drake Passage to the western Indian Ocean (from south of Australia to the center of the Pacific Ocean). This response can be further decomposed into (i) a zonally asymmetric component, a quasi-stationary wave train forced by a barotropic ridge downstream of the SSTA; and (ii) a zonal-mean component similar to a meridional shift of westerlies and hence a southern annular mode (SAM)-like pattern. The former component is phase locked with the SSTA position, while the latter has a phase that depends on the relative SSTA position with regard to the background quasi-stationary wave pattern. The study shows that the barotropic downstream ridge response is responsible for modifying the low-frequency eddy–mean flow interactions through relative vorticity fluxes and inducing the bipolar projection of the zonal-mean response.

Description: The interaction of internal waves with geostrophic flows is found to be strongly dependent upon the background stratification. Under the traditional approximation of neglecting the horizontal component of the earth’s rotation vector, the well-known inertial and symmetric instabilities highlight the asymmetry between positive and negative vertical components of relative vorticity (horizontal shear) of the mean flow, the former being stable. This is a strong stratification limit but, if it becomes too low, the traditional approximation cannot be made and the Coriolis terms caused by the earth’s rotation vector must be kept in full. A new asymmetry then appears between positive and negative horizontal components of relative vorticity (vertical shear) of the mean flow, the latter becoming more unstable. Particularly conspicuous at low latitudes, this new asymmetry does not require vanishing stratification to occur as it operates readily for rotation/stratification ratios 2Ω/N as small as 0.25 (the stratification still dominates over rotation) for realistic vertical shears. Given that such ratios are easily found in ocean–atmosphere boundary layers or in the deep ocean, such ageostrophic instabilities may be important for the routes to dissipation of the energy of the large-scale motions. The energetics show that, depending on the orientation of the internal wave crests with respect to the mean isopycnal surfaces, the unstable motions can draw their energy either from the kinetic energy or from the available potential energy of the mean flow. The kinetic energy source is usually the leading contribution when the growth rates reach their maxima.