ALBERT

Description: We study the quasi-geostrophic merging dynamics of axisymmetric baroclinic vortices to understand how baroclinicity affects merging rates and the development of the nonlinear cascade of enstrophy. The initial vortices are taken to simulate closely the horizontal' and vertical structure of Gulf Stream rings. A quasigeostrophic model is set with a horizontal resolution of 9 km and 6 vertical levels to resolve the mean stratification of the Gulf Stream region. The results show that the baroclinic merging is slower than the purely barotropic process, The merging is shown to occur in two phases: the tirst, which produces clove-shaped vortices and diffusive mixing of vorticity contours; and the second, which consists of the sliding of the remaining vorticity cores with a second diffusive mixing of the intemal vorticity field. Comparison among Nof, Cushman-Roisin, Polvani et al, and Dewar and Killworth merging events indicates a substantial agreement in the kinematics of the DYOCRSS. Parameter sensitivity experiments show that the decrease of the baroclinicity parameter of the system, Γ^2, [defined as Γ^2 = (D^2 fo^2)/ (No^2 H^2)], increases the speed of merging while its increase slows down the merging. However, the halting elfect of baroclinicity (large Γ^2 or small Rossby radii of deformation) reaches a saturation level where the merging becomes insensitive to larger F2 values. Furthermore, we show that a regime of small Γ^2 exists at which the merged baroclinic vortex is unstable (metastable) and breaks again into two new vortices, Thus, in the baroelinic case the range of Γ^2 detemines the stability of the merged vortex. We analyze these results by local energy and vorticity balances, showing that the horizontal divergence of pressure work term [∇ *(pv)] and the relative-vorticity advection term (v * ∇ (∇ ^2 φ) trigger the merging during the first phase. Due to this horizontal redistribution process, a net kinetic to gravitational energy conversion occurs via buoyancy work in the region external to the cores of the vortices. The second phase of merging is dominated by a direct baroclinic conversion of available gravitational energy into kinetic energy, which in tum triggers a horizontal energy redistribution producing the final fusion of the vortex centers. This energy and vorticity analysis supports the hypothesis that merging is an internal mixing process triggered by a horizontal redistribution of kinetic energy.

Description: The work has been financed by a grant from the Progetto Finalizzato "Calcolo Parallelo"

Description: Published

Description: 1618/1637

Description: 4A. Clima e Oceani

Description: JCR Journal

Description: restricted

Description: 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 2014

Description: This thesis explores the buoyancy-driven circulation in the Red Sea, using a combination of observations, as well as numerical modeling and analytical method. The first part of the thesis investigates the formation mechanism and spreading of Red Sea Overflow Water (RSOW) in the Red Sea. The preconditions required for open-ocean convection, which is suggested to be the formation mechanism of RSOW, are examined. The RSOW is identified and tracked as a layer with minimum potential vorticity and maximum chlorofluorocarbon-12. The pathway of the RSOW is also explored using numerical simulation. If diffusivity is not considered, the production rate of the RSOW is estimated to be 0.63 Sv using Walin’s method. By comparing this 0.63 Sv to the actual RSOW transport at the Strait of Bab el Mandeb, it is implied that the vertical diffusivity is about 3.4 x 10-5m2 s-1 . The second part of the thesis studies buoyancy-forced circulation in an idealized Red Sea. Buoyancy-loss driven circulation in marginal seas is usually dominated by cyclonic boundary currents on f-plane, as suggested by previous observations and numerical modeling. This thesis suggests that by including β-effect and buoyancy loss that increases linearly with latitude, the resultant mean Red Sea circulation consists of an anticyclonic gyre in the south and a cyclonic gyre in the north. In mid-basin, the northward surface flow crosses from the western boundary to the eastern boundary. The observational support is also reviewed. The mechanism that controls the crossover of boundary currents is further explored using an ad hoc analytical model based on PV dynamics. This ad hoc analytical model successfully predicts the crossover latitude of boundary currents. It suggests that the competition between advection of planetary vorticity and buoyancy-loss related term determines the crossover latitude. The third part of the thesis investigates three mechanisms that might account for eddy generation in the Red Sea, by conducting a series of numerical experiments. The three mechanisms are: i) baroclinic instability; ii) meridional structure of surface buoyancy losses; iii) cross-basin wind fields.

Description: This work is supported by Award Nos. USA 00002, KSA 00011 and KSA 00011/02 made by King Abdullah University of Science and Technology (KAUST) , National Science Foundation OCE0927017, and WHOI Academic Program Office.

Description: Author Posting. © American Meteorological Society, 2013. 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 26 (2013): 9839–9859, doi:10.1175/JCLI-D-12-00647.1.

Description: Spatial and temporal covariability between the atmospheric transient eddy heat fluxes (i.e., υ′T′ and υ′q′) in the Northern Hemisphere winter (January–March) and the paths of the Gulf Stream (GS), Kuroshio Extension (KE), and Oyashio Extension (OE) are examined based on an atmospheric reanalyses and ocean observations for 1979–2009. For the climatological winter mean, the northward heat fluxes by the synoptic (2–8 days) transient eddies exhibit canonical storm tracks with their maxima collocated with the GS and KE/OE. The intraseasonal (8 days–3 months) counterpart, while having overall similar amplitude, shows a spatial pattern with more localized maxima near the major orography and blocking regions. Lateral heat flux divergence by transient eddies as the sum of the two frequency bands exhibits very close coupling with the exact locations of the ocean fronts. Linear regression is used to examine the lead–lag relationship between interannual changes in the northward heat fluxes by the transient eddies and the meridional changes in the paths of the GS, KE, and OE, respectively. One to three years prior to the northward shifts of each ocean front, the atmospheric storm tracks shift northward and intensify, which is consistent with wind-driven changes of the ocean. Following the northward shifts of the ocean fronts, the synoptic storm tracks weaken in all three cases. The zonally integrated northward heat transport by the synoptic transient eddies increases by ~5% of its maximum mean value prior to the northward shift of each ocean front and decreases to a similar amplitude afterward.

Description: Support from the National Aeronautics and Space Administration (NASA) Physical Oceanography Program (NNX09AF35G to TJ and Y-OK) and the Department of Energy (DOE) Climate and Environmental Sciences Division (DE-SC0007052 to Y-OK) is gratefully acknowledged.

Description: 2014-06-15

Description: Author Posting. © American Meteorological Society, 2013. 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 43 (2013): 1398–1406, doi:10.1175/JPO-D-13-028.1.

Description: An adiabatic, inertial, and quasigeostrophic model is used to discuss the interaction of surface Ekman transport with an island. The theory extends the recent work of Spall and Pedlosky to include an analytical and nonlinear model for the interaction. The presence of an island that interrupts a uniform Ekman layer transport raises interesting questions about the resulting circulation. The consequential upwelling around the island can lead to a local intake of fluid from the geostrophic region beneath the Ekman layer or to a more complex flow around the island in which the fluid entering the Ekman layer on one portion of the island's perimeter is replaced by a flow along the island's boundary from a downwelling region located elsewhere on the island. This becomes especially pertinent when the flow is quasigeostrophic and adiabatic. The oncoming geostrophic flow that balances the offshore Ekman flux is largely diverted around the island, and the Ekman flux is fed by a transfer of fluid from the western to the eastern side of the island. As opposed to the linear, dissipative model described earlier, this transfer takes place even in the absence of a topographic skirt around the island. The principal effect of topography in the inertial model is to introduce an asymmetry between the circulation on the northern and southern sides of the island. The quasigeostrophic model allows a simple solution to the model problem with topography and yet the resulting three-dimensional circulation is surprisingly complex with streamlines connecting each side of the island.

Description: This research was supported in part by NSF Grant OCE Grant 0925061.

Description: Author Posting. © American Meteorological Society, 2014. 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 44 (2014): 229–245, doi:10.1175/JPO-D-12-0218.1.

Description: Data from a mooring deployed at the edge of the East Greenland shelf south of Denmark Strait from September 2007 to October 2008 are analyzed to investigate the processes by which dense water is transferred off the shelf. It is found that water denser than 27.7 kg m−3—as dense as water previously attributed to the adjacent East Greenland Spill Jet—resides near the bottom of the shelf for most of the year with no discernible seasonality. The mean velocity in the central part of the water column is directed along the isobaths, while the deep flow is bottom intensified and veers offshore. Two mechanisms for driving dense spilling events are investigated, one due to offshore forcing and the other associated with wind forcing. Denmark Strait cyclones propagating southward along the continental slope are shown to drive off-shelf flow at their leading edges and are responsible for much of the triggering of individual spilling events. Northerly barrier winds also force spilling. Local winds generate an Ekman downwelling cell. Nonlocal winds also excite spilling, which is hypothesized to be the result of southward-propagating coastally trapped waves, although definitive confirmation is still required. The combined effect of the eddies and barrier winds results in the strongest spilling events, while in the absence of winds a train of eddies causes enhanced spilling.

Description: The authors wish to thank Paula Fratantoni, Frank Bahr, and Dan Torres for processing the mooring data. The mooring array was capably deployed by the crew of the R/V Arni Fridriksson and recovered by the crew of the R/V Knorr. We thank Hedinn Valdimarsson for his assistance in the field work. Ken Brink provided valuable insights regarding the dynamics of shelf waves. Funding for the study was provided by National Science Foundation Grant OCE-0722694, the Arctic Research Initiative of the Woods Hole Oceanographic Institution. We also wish to thank the Natural Environment Research Council for Ph.D. studentship funding, and the University of East Anglia’s Roberts Fund and Royal Meteorological Society for supporting travel for collaboration.

Description: 2014-07-01

Description: Author Posting. © American Meteorological Society, 2014. 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 44 (2014): 834-849, doi:10.1175/JPO-D-13-0179.1.

Description: A hydrostatic numerical model with alongshore-uniform barotropic M2 tidal boundary forcing and idealized shelfbreak canyon bathymetries is used to study internal-tide generation and onshore propagation. A control simulation with Mid-Atlantic Bight representative bathymetry is supported by other simulations that serve to identify specific processes. The canyons and adjacent slopes are transcritical in steepness with respect to M2 internal wave characteristics. Although the various canyons are symmetrical in structure, barotropic-to-baroclinic energy conversion rates Cυ are typically asymmetrical within them. The resulting onshore-propagating internal waves are the strongest along beams in the horizontal plane, with the stronger beam in the control simulation lying on the side with higher Cυ. Analysis of the simulation results suggests that the cross-canyon asymmetrical Cυ distributions are caused by multiple-scattering effects on one canyon side slope, because the phase variation in the spatially distributed internal-tide sources, governed by variations in the orientation of the bathymetry gradient vector, allows resonant internal-tide generation. A less complex, semianalytical, modal internal wave propagation model with sources placed along the critical-slope locus (where the M2 internal wave characteristic is tangent to the seabed) and variable source phasing is used to diagnose the physics of the horizontal beams of onshore internal wave radiation. Model analysis explains how the cross-canyon phase and amplitude variations in the locally generated internal tides affect parameters of the internal-tide beams. Under the assumption that strong internal tides on continental shelves evolve to include nonlinear wave trains, the asymmetrical internal-tide generation and beam radiation effects may lead to nonlinear internal waves and enhanced mixing occurring preferentially on one side of shelfbreak canyons, in the absence of other influencing factors.

Description: All three authors were supported by Office of Naval Research (ONR) Grant N00014-11-1-0701. WGZ was additionally supported by the National Science Foundation (NSF) Grant OCE-1154575, and TFD was additionally supported by NSF Grant OCE-1060430.

Description: 2014-09-01

Description: Author Posting. © American Meteorological Society, 2014. 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 44 (2014): 2593–2616, doi:10.1175/JPO-D-13-0120.1.

Description: The first direct estimate of the rate at which geostrophic turbulence mixes tracers across the Antarctic Circumpolar Current is presented. The estimate is computed from the spreading of a tracer released upstream of Drake Passage as part of the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES). The meridional eddy diffusivity, a measure of the rate at which the area of the tracer spreads along an isopycnal across the Antarctic Circumpolar Current, is 710 ± 260 m2 s−1 at 1500-m depth. The estimate is based on an extrapolation of the tracer-based diffusivity using output from numerical tracers released in a one-twentieth of a degree model simulation of the circulation and turbulence in the Drake Passage region. The model is shown to reproduce the observed spreading rate of the DIMES tracer and suggests that the meridional eddy diffusivity is weak in the upper kilometer of the water column with values below 500 m2 s−1 and peaks at the steering level, near 2 km, where the eddy phase speed is equal to the mean flow speed. These vertical variations are not captured by ocean models presently used for climate studies, but they significantly affect the ventilation of different water masses.

Description: NSF support through Awards OCE-1233832, OCE-1232962, and OCE-1048926 is gratefully acknowledged.

Description: 2015-04-01

Description: Author Posting. © American Meteorological Society, 2014. 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 27 (2014): 2842–2860, doi:10.1175/JCLI-D-13-00227.1.

Description: Mooring measurements from the Kuroshio Extension System Study (June 2004–June 2006) and from the ongoing Kuroshio Extension Observatory (June 2004–present) are combined with float measurements of the Argo network to study the variability of the North Pacific Subtropical Mode Water (STMW) across the entire gyre, on time scales from days, to seasons, to a decade. The top of the STMW follows a seasonal cycle, although observations reveal that it primarily varies in discrete steps associated with episodic wind events. The variations of the STMW bottom depth are tightly related to the sea surface height (SSH), reflecting mesoscale eddies and large-scale variations of the Kuroshio Extension and recirculation gyre systems. Using the observed relationship between SSH and STMW, gridded SSH products and in situ estimates from floats are used to construct weekly maps of STMW thickness, providing nonbiased estimates of STMW total volume, annual formation and erosion volumes, and seasonal and interannual variability for the past decade. Year-to-year variations are detected, particularly a significant decrease of STMW volume in 2007–10 primarily attributable to a smaller volume formed. Variability of the heat content in the mode water region is dominated by the seasonal cycle and mesoscale eddies; there is only a weak link to STMW on interannual time scales, and no long-term trends in heat content and STMW thickness between 2002 and 2011 are detected. Weak lagged correlations among air–sea fluxes, oceanic heat content, and STMW thickness are found when averaged over the northwestern Pacific recirculation gyre region.

Description: This work was sponsored by the National Science Foundation (Grants OCE-0220161, OCE-0825152, and OCE-0827125).

Description: 2014-10-15

Description: Author Posting. © American Meteorological Society, 2014. 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 27 (2014): 3596–3618, doi:10.1175/JCLI-D-13-00070.1.

Description: Estimates of the recent mean and time varying water mass transformation rates associated with North Atlantic surface-forced overturning are presented. The estimates are derived from heat and freshwater surface fluxes and sea surface temperature fields from six atmospheric reanalyses—the Japanese 25-yr Reanalysis (JRA), the NCEP–NCAR reanalysis (NCEP1), the NCEP–U.S. Department of Energy (DOE) reanalysis (NCEP2), the European Centre for Medium-Range Weather Forecasts (ECMWF) Interim Re-Analysis (ERA-I), the Climate Forecast System Reanalysis (CFSR), and the Modern-Era Reanalysis for Research and Applications (MERRA)—together with sea surface salinity fields from two globally gridded datasets (World Ocean Atlas and Met Office EN3 datasets). The resulting 12 estimates of the 1979–2007 mean surface-forced streamfunction all depict a subpolar cell, with maxima north of 45°N, near σ = 27.5 kg m−3, and a subtropical cell between 20° and 40°N, near σ = 26.1 kg m−3. The mean magnitude of the subpolar cell varies between 12 and 18 Sv (1 Sv ≡ 106 m3 s−1), consistent with estimates of the overturning circulation from subsurface observations. Analysis of the thermal and haline components of the surface density fluxes indicates that large differences in the inferred low-latitude circulation are largely a result of the biases in reanalysis net heat flux fields, which range in the global mean from −13 to 19 W m−2. The different estimates of temporal variability in the subpolar cell are well correlated with each other. This suggests that the uncertainty associated with the choice of reanalysis product does not critically limit the ability of the method to infer the variability in the subpolar overturning. In contrast, the different estimates of subtropical variability are poorly correlated with each other, and only a subset of them captures a significant fraction of the variability in independently estimated North Atlantic Subtropical Mode Water volume.

Description: JPG is funded by UK Natural Environment Research Council New Investigator Grant NE/I001654/1. Y-OK was supported by the U.S. National Science Foundation under Grant OCE-0424492. RJB is supported by a fellowship from the UK National Centre for Earth Observation.

Description: 2014-11-15

Description: Author Posting. © American Meteorological Society, 2014. 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 44 (2014): 2498–2523, doi:10.1175/JPO-D-13-0183.1.

Description: This study examines the observability of a stratified ocean in a square flat basin on a midlatitude beta plane. Here, “observability” means the ability to establish, in a finite interval of time, the time-dependent ocean state given density observations over the same interval and with no regard for errors. The dynamics is linearized and hydrostatic, so that the motion can be decomposed into normal modes and the observability analysis is simplified. An observability Gramian (a symmetric matrix) is determined for the flows in an inviscid interior, in frictional boundary layers, and in a closed basin. Its properties are used to establish the condition for complete observability and to identify optimal data locations for each of these flows. It is found that complete observability of an oceanic interior in time-dependent Sverdrup balance requires that the observations originate from the westernmost location at each considered latitude. The degree of observability increases westward due to westward propagation of long baroclinic Rossby waves: data collected in the west are more informative than data collected in the east. Likewise, the best locations for observing variability in the western (eastern) boundary layer are near (far from) the boundary. The observability of a closed basin is influenced by the westward propagation and the boundaries. Optimal data locations that are identified for different resolutions (0.01 to 1 yr) and lengths of data records (0.2 to 20 yr) show a variable influence of the planetary vorticity gradient. Data collected near the meridional boundaries appear always less informative, from the viewpoint of basin observability, than data collected away from these boundaries.

Description: This work was supported by the U.S. National Science Foundation.

Description: 2015-03-01

Description: 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 1996

Description: This thesis addresses the question of how a highly energetic eddy field could be generated in the interior of the ocean away from the swift boundary currents. The energy radiation due to the temporal growth of non-trapped (radiating) disturbances in such a boundary current is thought to be one of the main sources for the described variability. The problem of stability of an energetic current, such as the Gulf Stream, is formulated. The study then focuses on the ability of the current to support radiating instabilities capable of significant penetration into the far-field and their development with time. The conventional model of the Gulf Stream as a zonal current is extended to allow the jet axis to make an angle to a latitude circle. The linear stability of such a nonzonal flow, uniform in the along-jet direction on a beta-plane, is first studied. The stability computations are performed for piece-wise constant and continuous velocity profiles. New stability properties of nonzonal jets are discussed. In particular, the destabilizing effect of the meridional tilt of the jet axis is demonstrated. The radiating properties of nonzonal currents are found to be very different from those of zonal currents. In particular, purely zonal flows do not support radiating instabilities, whereas flows with a meridional component are capable of radiating long and slowly growing waves. The nonlinear terms are then included in the consideration and the effects of the nonlinear interactions on the radiating properties of the solution are studied in detail. For these purposes, the efficient numerical code for solving equation for the QG potential vorticity with open boundary conditions of Orlanski's type is constructed. The results show that even fast growing linear solutions, which are trapped during the linear stage of developement, can radiate energy in the nonlinear regime if the basic current is nonzonal. The radiation starts as soon as the initial fast exponential growth significantly slows. The initial trapping of those solutions is caused by their fast temporal growth. The new mechanism for radiation is related to the nonzonality of a current.

Description: This work was supported by NSF Grant OCE 9301845.

Description: 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 February 2013

Description: Studying oceanic eddies is important for understanding and predicting ocean circulation and climate variability. The central focus of this dissertation is the energy exchange between eddies and mean flow and banded structures in the low-frequency component of the eddy field. A combination of a realistic eddy-permitting ocean state estimate and simplified theoretical models is used to address the following specific questions. (1) What are the major spatial characteristics of eddy-mean flow interaction from an energy perspective? Is eddy-mean flow interaction a local process in most ocean regions? (2) The banded structures in the low-frequency eddy field are termed striations. How much oceanic variability is associated with striations? How does the time-mean circulation, for example a subtropical gyre or constant mean flow, influence the origin and characteristics of striations? How much do striations contribute to the energy budget and tracer mixing?

Description: This research was supported by the National Aeronautics and Space Administration contracts NNX09AI87G and NNX08AR33G.

Description: Author Posting. © American Meteorological Society, 2013. 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 43 (2013): 283–300, doi:10.1175/JPO-D-11-0240.1.

Description: Motivated by the recent interest in ocean energetics, the widespread use of horizontal eddy viscosity in models, and the promise of high horizontal resolution data from the planned wide-swath satellite altimeter, this paper explores the impacts of horizontal eddy viscosity and horizontal grid resolution on geostrophic turbulence, with a particular focus on spectral kinetic energy fluxes Π(K) computed in the isotropic wavenumber (K) domain. The paper utilizes idealized two-layer quasigeostrophic (QG) models, realistic high-resolution ocean general circulation models, and present-generation gridded satellite altimeter data. Adding horizontal eddy viscosity to the QG model results in a forward cascade at smaller scales, in apparent agreement with results from present-generation altimetry. Eddy viscosity is taken to roughly represent coupling of mesoscale eddies to internal waves or to submesoscale eddies. Filtering the output of either the QG or realistic models before computing Π(K) also greatly increases the forward cascade. Such filtering mimics the smoothing inherent in the construction of present-generation gridded altimeter data. It is therefore difficult to say whether the forward cascades seen in present-generation altimeter data are due to real physics (represented here by eddy viscosity) or to insufficient horizontal resolution. The inverse cascade at larger scales remains in the models even after filtering, suggesting that its existence in the models and in altimeter data is robust. However, the magnitude of the inverse cascade is affected by filtering, suggesting that the wide-swath altimeter will allow a more accurate determination of the inverse cascade at larger scales as well as providing important constraints on smaller-scale dynamics.

Description: BKA received support from Office of Naval Research Grant N00014-11-1-0487, National Science Foundation (NSF) Grants OCE-0924481 and OCE- 09607820, and University of Michigan startup funds. KLP acknowledges support from Woods Hole Oceanographic Institution bridge support funds. RBS acknowledges support from NSF grants OCE-0960834 and OCE-0851457, a contract with the National Oceanography Centre, Southampton, and a NASA subcontract to Boston University. JFS and JGR were supported by the projects ‘‘Global and remote littoral forcing in global ocean models’’ and ‘‘Agesotrophic vorticity dynamics of the ocean,’’ respectively, both sponsored by the Office of Naval Research under program element 601153N.

Description: 2013-08-01

Description: Author Posting. © American Meteorological Society, 2012. 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 42 (2012): 2234–2253, doi:10.1175/JPO-D-12-033.1.

Description: Meridional velocity, mass, and heat transport in the equatorial oceans are difficult to estimate because of the nonapplicability of the geostrophic balance. For this purpose a steady-state model is utilized in the equatorial Indian Ocean using NCEP wind stress and temperature and salinity data from the World Ocean Atlas 2005 (WOA05) and Argo. The results show a Somali Current flowing to the south during the winter monsoon carrying −11.5 ± 1.3 Sv (1 Sv ≡ 106 m3 s−1) and −12.3 ± 0.3 Sv from WOA05 and Argo, respectively. In the summer monsoon the Somali Current reverses to the north transporting 16.8 ± 1.2 Sv and 19.8 ± 0.6 Sv in the WOA05 and Argo results. Transitional periods are considered together and in consequence, there is not a clear Somali Current present in this period. Model results fit with in situ measurements made around the region, although Argo data results are quite more realistic than WOA05 data results.

Description: This study has been partly funded by the MOC Project (CTM 2008- 06438) and the Spanish contribution to the Argo network (AC2009 ACI2009-0998), financed by the Spanish Government and Feder.

Description: 2013-06-01

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 2013

Description: Between 2002 and 2011 a single mooring was maintained in the core of the Pacific Water boundary current in the Alaskan Beaufort Sea near 152° W. Using velocity and hydrographic data from six year-long deployments during this time period, we examine the interannual variability of the current. It is found that the volume, heat, and freshwater transport have all decreased drastically over the decade, by more than 80%. The most striking changes have occurred during the summer months. Using a combination of weather station data, atmospheric reanalysis fields, and concurrent shipboard and mooring data from the Chukchi Sea, we investigate the physical drivers responsible for these changes. It is demonstrated that an increase in summertime easterly winds along the Beaufort slope is the primary reason for the drop in transport. The intensification of the local winds has in turn been driven by a strengthening of the summer Beaufort High in conjunction with a deepening of the summer Aleutian Low. Since the fluxes of mass, heat, and freshwater through Bering Strait have increased over the same time period, this raises the question as to the fate of the Pacific water during recent years and its impacts. We present evidence that more heat has been fluxed directly into the interior basin from Barrow Canyon rather than entering the Beaufort shelfbreak jet, and this is responsible for a significant portion of the increased ice melt in the Pacific sector of the Arctic Ocean.

Description: The majority of the data for this project was funded by grant # ARC-0856244 from the O ce of Polar Programs of the National Science Foundation. My time at WHOI was funded by the United States Navy, the National Science Foundation Graduate Research Fellowship Program and the WHOI Academic Programs O ffice.

Description: 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 1996

Description: A convection experiment was done with a rotating rectangular tank as a model of oceanic meridional overturning circulation. Heat flux was fixed at one bottom end of the tank using an electrical heater. Temperature was fixed at the other end using a cooling plate. All other boundaries were insulated. The cross sections of temperature field were made at several locations. In equilibrium, the heat input to the fluid H was the same as the meridional heat flux (heat flux from the source to the sink), so it was possible to find a scaling law relating H to the temperature difference across the tank ΔT and rotation rate f. The experimental result suggests that the meridional heat transport in the experiment was mostly due to geostrophic flows with a minor correction caused by the bottom friction. If there was no friction, the scaling law from the experiment resembles the one verified in part in the numerical model by Bryan and Cox (1967). Flow visualization and temperature sections showed that there were meridional geostrophic currents that transported heat. When the typical values of the North Atlantic are introduced, the geostrophic scaling law predicts meridional heat flux comparable to that estimated in the North Atlantic when the vertical eddy diffusivity of heat is about 1cm2s-1. Naturally, this experiment is a only crude model of the oceanic convective circulation. We do not claim that the geostrophic scaling applies in detail to the oceans, however, it may have some important use in climate modeling. For example, almost all existing box models and two-dimensional numerical models of ocean circulation use a frictional scaling law for buoyancy transport. A box model with the geostrophic scaling law is shown to be more robust to a change in the boundary forcing so that it is less likely to have a thermohaline catastrophic transition under the present conditions. It is also shown that a restoring boundary condition for salinity introduces stability to a thermal mode circulation, unless the restoring time for salinity is several orders of magnitude larger than that for temperature.

Description: This study has been funded by NSF grant number OCE92-01464 and Korean Government Overseas Scholarship Grant.

Description: 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 February 2013

Description: Large-scale thermal forcing and freshwater fluxes play an essential role in setting temperature and salinity in the ocean. A number of recent estimates of the global oceanic freshwater balance as well as the global oceanic surface net heat flux are used to investigate the effects of heat- and freshwater forcing at the ocean surface. Such forcing induces changes in both density and density-compensated temperature and salinity changes (’spice’). The ratio of the relative contributions of haline and thermal forcing in the mixed layer is maintained by large-scale surface fluxes, leading to important consequences for mixing in the ocean interior. In a stratified ocean, mixing processes can be either along lines of constant density (isopycnal) or across those lines (diapycnal). The contribution of these processes to the total mixing rate in the ocean can be estimated from the large-scale forcing by evaluating the production of thermal variance, salinity variance and temperature-salinity covariance. Here, I use new estimates of surface fluxes to evaluate these terms and combine them to generate estimates of the production of density and spice variance under the assumption of a linear equation of state. As a consequence, it is possible to estimate the relative importance of isopycnal and diapycnal mixing in the ocean. While isopycnal and diapycnal processes occur on very different length scales, I find that the surface-driven production of density and spice variance requires an approximate equipartition between isopycnal and diapycnal mixing in the ocean interior. In addition, consideration of the full nonlinear equation of state reveals that surface fluxes require an apparent buoyancy gain (expansion) of the ocean, which allows an estimate of the amount of contraction on mixing due to cabbeling in the ocean interior.

Description: The author would like to acknowledge support from the National Aeronautics and Space Administration, grant #NNX12AF59G and the National Science Foundation, grant #OCE-0647949.

Description: Author Posting. © American Meteorological Society, 2012. 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 42 (2012): 2206–2228, doi:10.1175/JPO-D-11-0191.1.

Description: This study investigates the anisotropic properties of the eddy-induced material transport in the near-surface North Atlantic from two independent datasets, one simulated from the sea surface height altimetry and one derived from real-ocean surface drifters, and systematically examines the interactions between the mean- and eddy-induced material transport in the region. The Lagrangian particle dispersion, which is widely used to characterize the eddy-induced tracer fluxes, is quantified by constructing the “spreading ellipses.” The analysis consistently demonstrates that this dispersion is spatially inhomogeneous and strongly anisotropic. The spreading is larger and more anisotropic in the subtropical than in the subpolar gyre, and the largest ellipses occur in the Gulf Stream vicinity. Even at times longer than half a year, the spreading exhibits significant nondiffusive behavior in some parts of the domain. The eddies in this study are defined as deviations from the long-term time-mean. The contributions from the climatological annual cycle, interannual, and subannual (shorter than one year) variability are investigated, and the latter is shown to have the strongest effect on the anisotropy of particle spreading. The influence of the mean advection on the eddy-induced particle spreading is investigated using the “eddy-following-full-trajectories” technique and is found to be significant. The role of the Ekman advection is, however, secondary. The pronounced anisotropy of particle dispersion is expected to have important implications for distributing oceanic tracers, and for parameterizing eddy-induced tracer transfer in non-eddy-resolving models.

Description: IR was supported by Grant NSF-OCE-0725796. IK would like to acknowledge support by the National Science foundation Grant OCE-0842834.

Description: 2013-06-01

Description: Author Posting. © American Meteorological Society, 2013. 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 43 (2013): 905–919, doi:10.1175/JPO-D-12-0150.1.

Description: Interactions between vortices and a shelfbreak current are investigated, with particular attention to the exchange of waters between the continental shelf and slope. The nonlinear, three-dimensional interaction between an anticyclonic vortex and the shelfbreak current is studied in the laboratory while varying the ratio ε of the maximum azimuthal velocity in the vortex to the maximum alongshelf velocity in the shelfbreak current. Strong interactions between the shelfbreak current and the vortex are observed when ε 〉 1; weak interactions are found when ε 〈 1. When the anticyclonic vortex comes in contact with the shelfbreak front during a strong interaction, a streamer of shelf water is drawn offshore and wraps anticyclonically around the vortex. Measurements of the offshore transport and identification of the particle trajectories in the shelfbreak current drawn offshore from the vortex allow quantification of the fraction of the shelfbreak current that is deflected onto the slope; this fraction increases for increasing values of ε. Experimental results in the laboratory are strikingly similar to results obtained from observations in the Middle Atlantic Bight (MAB); after proper scaling, measurements of offshore transport and offshore displacement of shelf water for vortices in the MAB that span a range of values of ε agree well with laboratory predictions.

Description: Laboratory work was supported by the National Science Foundation through Grant OCE- 0081756. Glider observations in March–April 2006 were supported by the National Science Foundation through Grant OCE-0220769. Glider observations in July– October 2007 were supported by a grant from Raytheon. RET was supported by the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Cooperative Institute for the North Atlantic Region. The REMUS observations were funded by the Office of Naval Research. GGG was supported by the National Science Foundation through Grant OCE-1129125 for analysis and writing.

Description: 2013-11-01

Description: 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 December 1996

Description: The transformation of potential vorticity within and stability of nonlinear deep western boundary currents in an idealized tropical ocean are studied using a shallowwater model. Observational evidence indicates that the potential vorticity of fluid parcels in deep western boundary currents must change sign as they cross the equator, but this evidence is otherwise unable to clarify the process. A series of numerical experiments investigate this transformation in a rectangular basin straddling the equator. A mass source located in the northwestern corner feeds fluid into the domain where it is constrained to cross the equator to reach a distributed mass sink. Dissipation is included as momentum diffusion. The Reynolds number, defined as the ratio of the mass source per unit depth to the viscosity, determines the nature of the flow, and a critical value, Rec, divides its possible behavior into two regimes. For Re 〈 Rec, the flow is laminar and well described by linear theory. For Re just above the critical value, the flow is time-dependent, with cyclonic eddies forming in the western boundary current near the equator. For still larger Reynolds number, eddies of both signs emerge and form a complicated, interacting network that extends into the basin several deformation radii from the western boundary, as well as north and south of the equator. The eddy field is established as the mechanism for potential vorticity transformation in nonlinear cross-equatorial flow. The analysis of vorticity fluxes follows from the flux-conservative form of the absolute vorticity equation. It is shown that the zonally integrated meridional flux of vorticity across the equator using no slip boundary conditions is virtually zero even in the strongly nonlinear limit suggesting that the eddies are extremely efficient vorticity transfer agents. A decomposition of the vorticity fluxes into components due to mean advection, eddy transport, and friction, reveals the growth with Reynolds number of a turbulent boundary layer that exchanges vorticity between the inertial portion of the boundary current and a frictional sub-layer where modification is straightforward. A linear stability analysis of the shallow-water system in the tropical ocean examines the initial formation of the eddy field. The formulation assumes that the basic state is purely meridional and on a local f-plane. Realistic western boundary current profiles undergo a horizontal shear instability that is partially stabilized by viscosity. Calculations at several latitudes indicate that the instability is enhanced in the tropics where the internal deformation radius is a maximum. The linear stability analysis predicts a length scale of the disturbance, a location for its origin, and a critical Reynolds number that agree well with numerical results.

Description: Financial support for this research was provided by NSF grant number OCE- 9115915 and ONR ASSERT grant number N00014-94-1-0844.

Description: Author Posting. © American Meteorological Society, 2012. 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 42 (2012): 2283–2296, doi:10.1175/JPO-D-11-0227.1.

Description: The dynamic influence of thermohaline circulation on wind-driven circulation in the South China Sea (SCS) is studied using a simple reduced gravity model, in which the upwelling driven by mixing in the abyssal ocean is treated in terms of an upward pumping distributed at the base of the upper layer. Because of the strong upwelling of deep water, the cyclonic gyre in the northern SCS is weakened, but the anticyclonic gyre in the southern SCS is intensified in summer, while cyclonic gyres in both the southern and northern SCS are weakened in winter. For all seasons, the dynamic influence of thermohaline circulation on wind-driven circulation is larger in the northern SCS than in the southern SCS. Analysis suggests that the upwelling associated with the thermohaline circulation in the deep ocean plays a crucial role in regulating the wind-driven circulation in the upper ocean.

Description: G. Wang is supported by the National Science Foundation of China (NSFC Grants 41125019, 40725017, and 40976017).D.Chen is supported by grants from the Ministry of Science and Technology (2010DFA21012), the State Oceanic Administration (201105018), and the NSFC (91128204).

Description: 2013-06-01

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution May 1996

Description: The water mass distribution in the southwestern Barents Sea, the thermohaline structure of the western Barents Sea Polar Front, and the formation of local water masses are described based on an analysis of historical hydrographic data and a recent process-oriented field experiment. This study concentrated on the frontal region between Bj0rn0ya and Hopen Island where Arctic water is found on the Spitzbergen Bank and Atlantic Water in the Bear Island Trough and Hopen Trench. Distributions of Atlantic, Arctic, and Polar Front waters are consistent with topographic control of Atlantic water circulation. Seasonal buoyancy forcing disrupts the topographic control in the surface layer, altering the frontal structure, and affecting local water mass formation. In the winter, the topographic control is firmly established and both sides of the front are vertically well-mixed. Winter cooling creates sea-ice over Spitzbergen Bank and convectively formed Modified Atlantic Water in the Bear Island Trough and Hopen Trench. In the summer, heating melts the sea-ice, producing a surface meltwater pool that can cross the polar front, disrupting topographic control and substantially increasing the vertical thermohaline gradients in the frontal region. The meltwater pool produces the largest geostrophic shear in the region.

Description: Support for this work was provided by a Department of Defense National Defense Science and Engineering Graduate Fellowship and Office of Naval Research grant N00014- 90-J-1359.

Description: 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 June 2013

Description: This thesis explores the role that the circulation in the Gulf of Maine (GOM) plays in determining the distribution of dense aggregations of copepods. These aggregations are an important part of the marine ecosystem, especially for endangered North Atlantic right whales. Certain ocean processes may generate dense copepod aggregations, while others may destroy them; this thesis looks at how different characteristics of the GOM circulation fit into these two categories. The first part of the thesis investigates a hypothetical aggregation mechanism in which frontal circulation interacts with copepod behavior to generate a dense patch of copepods. The first two chapters of this thesis address this mechanism in the context of coastal river plumes and salinity fronts. One chapter describes the characteristics and variability of coastal freshwater and salinity fronts using a historical dataset and a realistic numerical model. The seasonal variability of freshwater is tied in part to seasonality in river discharge, while variability on shorter time scales in the frontal position is related to wind stress. Another chapter applies the hypothetical mechanism to idealized river plumes using a suite of numerical models. The structure of the plume and plume-relative circulation change the resulting copepod aggregation from what is expected from the hypothetical mechanism. The second part of the thesis discusses the GOM circulation and how it may eliminate copepod patches. The summertime mean surface circulation and eddy kinetic energy are computed from a Lagrangian drifter dataset and an adaptive technique that allows for higher spatial resolution while also keeping uncertainty low. Eddy diffusivity is also computed over different regions of the GOM in an attempt to quantify the spreading of a patch of copepods, and is found to be lower near the coast where right whales are often found feeding on copepod patches. In the next chapter, a numerical drifter dataset is used to understand how the results of the previous chapter depend upon the quantity of observations. It is found that the uncertainty in estimating eddy diffusivity is tightly coupled to the number of drifters in the calculation.

Description: This work was supported by the WHOI Coastal Ocean Institute Graduate Student Fellowship and Student Research Award, the WHOI Academic Programs O ce, the NOAA National Marine Fisheries Service Northeast Fisheries Science Center (NOAA Cooperative Agreement NA09OAR4320129), and the O ce of Naval Research Marine Mammals and Biology Program (Grant N00014-12-1-0208).

Description: Author Posting. © American Meteorological Society, 2013. 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 43 (2013): 1611–1626, doi:10.1175/JPO-D-12-0204.1.

Description: A new method is proposed for extrapolating subsurface velocity and density fields from sea surface density and sea surface height (SSH). In this, the surface density is linked to the subsurface fields via the surface quasigeostrophic (SQG) formalism, as proposed in several recent papers. The subsurface field is augmented by the addition of the barotropic and first baroclinic modes, whose amplitudes are determined by matching to the sea surface height (pressure), after subtracting the SQG contribution. An additional constraint is that the bottom pressure anomaly vanishes. The method is tested for three regions in the North Atlantic using data from a high-resolution numerical simulation. The decomposition yields strikingly realistic subsurface fields. It is particularly successful in energetic regions like the Gulf Stream extension and at high latitudes where the mixed layer is deep, but it also works in less energetic eastern subtropics. The demonstration highlights the possibility of reconstructing three-dimensional oceanic flows using a combination of satellite fields, for example, sea surface temperature (SST) and SSH, and sparse (or climatological) estimates of the regional depth-resolved density. The method could be further elaborated to integrate additional subsurface information, such as mooring measurements.

Description: JW and AM were supported by NASA (NNX12AD47G) and NSF (OCE 0928617). JLM was supported by the Office of Naval Research and the Office of Science (BER), U.S. Department of Energy under DE-GF0205ER64119. GRF is supported by OCE-0752346 and JHL by NORSEE (Nordic Seas Eddy Exchanges) funded by the Norwegian Research Council.

Description: 2014-02-01

Description: Author Posting. © American Meteorological Society, 2013. 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 26 (2013): 1669–1684, doi:10.1175/JCLI-D-12-00246.1.

Description: Climate change west of the Antarctic Peninsula is the most rapid of anywhere in the Southern Hemisphere, with associated changes in the rates and distributions of freshwater inputs to the ocean. Here, results from the first comprehensive survey of oxygen isotopes in seawater in this region are used to quantify spatial patterns of meteoric water (glacial discharge and precipitation) separately from sea ice melt. High levels of meteoric water are found close to the coast, due to orographic effects on precipitation and strong glacial discharge. Concentrations decrease offshore, driving significant southward geostrophic flows (up to ~30 cm s−1). These produce high meteoric water concentrations at the southern end of the sampling grid, where collapse of the Wilkins Ice Shelf may also have contributed. Sea ice melt concentrations are lower than meteoric water and patchier because of the mobile nature of the sea ice itself. Nonetheless, net sea ice production in the northern part of the sampling grid is inferred; combined with net sea ice melt in the south, this indicates an overall southward ice motion. The survey is contextualized temporally using a decade-long series of isotope data from a coastal Antarctic Peninsula site. This shows a temporal decline in meteoric water in the upper ocean, contrary to expectations based on increasing precipitation and accelerating deglaciation. This is driven by the increasing occurrence of deeper winter mixed layers and has potential implications for concentrations of trace metals supplied to the euphotic zone by glacial discharge. As the regional freshwater system evolves, the continuing isotope monitoring described here will elucidate the ongoing impacts on climate and the ecosystem.

Description: The Palmer LTER participants acknowledge Award 0823101 from the Organisms and Ecosystems program in NSF OPP

Description: 2013-09-01

Description: Author Posting. © American Meteorological Society, 2011. 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 Atmospheric and Oceanic Technology 28 (2011): 1539–1553, doi:10.1175/JTECH-D-11-00001.1.

Description: Turbulent Reynolds stresses are now routinely estimated from acoustic Doppler current profiler (ADCP) measurements in estuaries and tidal channels using the variance method, yet biases due to surface gravity waves limit its use in the coastal ocean. Recent modifications to this method, including spatially filtering velocities to isolate the turbulence from wave velocities and fitting a cospectral model to the below-wave band cospectra, have been used to remove this bias. Individually, each modification performed well for the published test datasets, but a comparative analysis over the range of conditions in the coastal ocean has not yet been performed. This work uses ADCP velocity measurements from five previously published coastal ocean and estuarine datasets, which span a range of wave and current conditions as well as instrument configurations, to directly compare methods for estimating stresses in the presence of waves. The computed stresses from each were compared to bottom stress estimates from a quadratic drag law and, where available, estimates of wind stress. These comparisons, along with an analysis of the cospectra, indicated that spectral fitting performs well when the wave climate is wide-banded and/or multidirectional as well as when instrument noise is high. In contrast, spatial filtering performs better when waves are narrow-banded, low frequency, and when wave orbital velocities are strong relative to currents. However, as spatial filtering uses vertically separated velocity bins to remove the wave bias, spectral fitting is able to resolve stresses over a larger fraction of the water column.

Description: J. Rosman acknowledges funding from the National Science Foundation (OCE-1061108).

Description: 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 June 2012

Description: Observations from a three-year field program on the inner shelf south of Martha's Vineyard, MA and a numerical model are used to describe the effect of stratification on inner shelf circulation, transport, and sediment resuspension height. Thermal stratification above the bottom mixed layer is shown to cap the height to which sediment is resuspended. Stratification increases the transport driven by cross-shelf wind stresses, and this effect is larger in the response to offshore winds than onshore winds. However, a one-dimensional view of the dynamics is not sufficient to explain the relationship between circulation and stratification. An idealized, cross-shelf transect in a numerical model (ROMS) is used to isolate the effects of stratification, wind stress magnitude, surface heat flux, cross-shelf density gradient, and wind direction on the inner shelf response to the cross-shelf component of the wind stress. In well mixed and weakly stratified conditions, the cross-shelf density gradient can be used to predict the transport efficiency of the cross-shelf wind stress. In stratified conditions, the presence of an along-shelf wind stress component makes the inner shelf response to cross-shelf wind stress strongly asymmetric.

Description: This work was supported through National Science Foundation grant no. OCE-0548961, the WHOI Academic Programs Office, and the WHOI Coastal Ocean Institute.

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 2012

Description: Interactions between the ocean circulation in sub-ice shelf cavities and the overlying ice shelf have received considerable attention in the context of observed changes in flow speeds of marine ice sheets around Antarctica. Modeling these interactions requires parameterizing the turbulent boundary layer processes to infer melt rates from the oceanic state at the ice-ocean interface. Here we explore two such parameterizations in the context of the MIT ocean general circulation model coupled to the z-coordinates ice shelf cavity model of Losch (2008). We investigate both idealized ice shelf cavity geometries as well as a realistic cavity under Pine Island Ice Shelf (PIIS), West Antarctica. Our starting point is a three-equation melt rate parameterization implemented by Losch (2008), which is based on the work of Hellmer and Olbers (1989). In this form, the transfer coefficients for calculating heat and freshwater fluxes are independent of frictional turbulence induced by the proximity of the moving ocean to the fixed ice interface. More recently, Holland and Jenkins (1999) have proposed a parameterization in which the transfer coefficients do depend on the ocean-induced turbulence and are directly coupled to the speed of currents in the ocean mixed layer underneath the ice shelf through a quadratic drag formulation and a bulk drag coefficient. The melt rate parameterization in the MITgcm is augmented to account for this velocity dependence. First, the effect of the augmented formulation is investigated in terms of its impact on melt rates as well as on its feedback on the wider sub-ice shelf circulation. We find that, over a wide range of drag coefficients, velocity-dependent melt rates are more strongly constrained by the distribution of mixed layer currents than by the temperature gradient between the shelf base and underlying ocean, as opposed to velocity-independent melt rates. This leads to large differences in melt rate patterns under PIIS when including versus not including the velocity dependence. In a second time, the modulating effects of tidal currents on melting at the base of PIIS are examined. We find that the temporal variability of velocity-dependent melt rates under tidal forcing is greater than that of velocity-independent melt rates. Our experiments suggest that because tidal currents under PIIS are weak and buoyancy fluxes are strong, tidal mixing is negligible and tidal rectification is restricted to very steep bathymetric features, such as the ice shelf front. Nonetheless, strong tidally-rectified currents at the ice shelf front significantly increase ablation rates there when the formulation of the transfer coefficients includes the velocity dependence. The enhanced melting then feedbacks positively on the rectified currents, which are susceptible to insulate the cavity interior from changes in open ocean conditions.

Description: National Science and Engineering Research Council of Canada

Description: Author Posting. © American Meteorological Society, 2011. 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 41 (2011): 2168–2186, doi:10.1175/JPO-D-11-08.1.

Description: This paper studies the interaction of an Antarctic Circumpolar Current (ACC)–like wind-driven channel flow with a continental slope and a flat-bottomed bay-shaped shelf near the channel’s southern boundary. Interaction between the model ACC and the topography in the second layer induces local changes of the potential vorticity (PV) flux, which further causes the formation of a first-layer PV front near the base of the topography. Located between the ACC and the first-layer slope, the newly formed PV front is constantly perturbed by the ACC and in turn forces the first-layer slope with its own variability in an intermittent but persistent way. The volume transport of the slope water across the first-layer slope edge is mostly directly driven by eddies and meanders of the new front, and its magnitude is similar to the maximum Ekman transport in the channel. Near the bay’s opening, the effect of the topographic waves, excited by offshore variability, dominates the cross-isobath exchange and induces a mean clockwise shelf circulation. The waves’ propagation is only toward the west and tends to be blocked by the bay’s western boundary in the narrow-shelf region. The ensuing wave–coast interaction amplifies the wave amplitude and the cross-shelf transport. Because the interaction only occurs near the western boundary, the shelf water in the west of the bay is more readily carried offshore than that in the east and the mean shelf circulation is also intensified along the bay’s western boundary.

Description: Y. Zhang acknowledges the support of the MIT-WHOI Joint Program in Physical Oceanography and NSF OCE-9901654 and OCE- 0451086. J. Pedlosky acknowledges the support of NSF OCE-9901654 and OCE-0451086.

Description: Author Posting. © American Meteorological Society, 2012. 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 25 (2012): 343–349, doi:10.1175/JCLI-D-11-00059.1.

Description: The Equatorial Undercurrent (EUC) is a major component of the tropical Pacific Ocean circulation. EUC velocity in most global climate models is sluggish relative to observations. Insufficient ocean resolution slows the EUC in the eastern Pacific where nonlinear terms should dominate the zonal momentum balance. A slow EUC in the east creates a bottleneck for the EUC to the west. However, this bottleneck does not impair other major components of the tropical circulation, including upwelling and poleward transport. In most models, upwelling velocity and poleward transport divergence fall within directly estimated uncertainties. Both of these transports play a critical role in a theory for how the tropical Pacific may change under increased radiative forcing, that is, the ocean dynamical thermostat mechanism. These findings suggest that, in the mean, global climate models may not underrepresent the role of equatorial ocean circulation, nor perhaps bias the balance between competing mechanisms for how the tropical Pacific might change in the future. Implications for model improvement under higher resolution are also discussed.

Description: KBK gratefully acknowledges the J. Lamar Worzel Assistant Scientist Fund. GCJ is supported by NOAA’s Office of Oceanic and Atmospheric Research. RM gratefully acknowledges the generous support and hospitality of the Divecha Centre for Climate Change and CAOS at IISc, Bangalore, and partial support by NASA PO grants.

Description: 2012-07-01

Description: 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 February 1993

Description: Reciprocal acoustic transmissions made in a region just south of the Gulf Stream are analyzed to determine the structure and variability of temperature, current velocity, and vorticity fields at the northern extent of the southern recirculation gyre. For ten months (November, 1988 through August, 1989), a pentagonal array of tomographic transceivers was situated in a region centered at 38°N, 55°W as part of the eastern array of the SYNOP (SYNoptic Ocean Prediction) Experiment. The region of focus is one rich in mesoscale energy, with the influence of local Gulf Stream meandering and cold-core ring activity strikingly evident. Daily-averaged acoustic transmissions yielded travel times which were inverted to obtain estimates of range-averaged temperature and current velocity fields, and area-averaged relative vorticity fields. The acoustically determined estimates are consistent with nearby current meter measurements and satellite infrared imagery. The signature of cold-core rings is clearly evident in the sections. Spectral estimates of the fields are dominated by motions with periodicities ranging from 32-128 days. Second-order statistics, such as eddy kinetic energies, and heat and momentum fluxes, are also estimated. The integrating nature of the tomographic measurement has been exploited to shed some light on the radiation of eddy energy from the Gulf Stream. The Eliassen-Palm flux diagnostic has been applied to an investigation of wave radiation from the Gulf Stream. Results of the diagnosis suggest that the Gulf Stream itself is the source of wave energy radiating into the far field and found in the interior of the North Atlantic subtropical gyre.

Description: This research was carried out under Office of Naval Research (ONR) University Research Initiative contract N00014-86-K-0751 and ONR contract N00014- 90-J-1481. Construction of the tomographic instruments was supported by grants and contracts with MIT: National Science Foundation grant OCE 85-12430 and by ONR. The field work was supported by ONR under contract N00014-85-G-0241 (Secretary of the Navy Professorship (C. Wunsch)).

Description: 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 March 1989

Description: The general theme of this thesis is the study of systematic mathematical techniques for determining the ocean circulation from classical hydrographic data. Two aspects of this theme are analyzed. The first is finding an efficient representation of hydrographic structure so as to make it most useful and informative. The second is application of inverse methods to the data to determine ocean circulation. Both subjects are examined in the North Atlantic Ocean. The efficient representation is examined in terms of empirical orthogonal functions (EOFs) among the variations in vertical hydrographic profiles. The data used are of a new set of high quality hydrography, all obtained in the early 1980s. Common EOFs are examined among temperature, salinity, oxygen, phosphate, silicate, and nitrate. The EOFs identify a fundamental simplicity in the spatial distributions of t hese properties. Although the volume of numbers involved in the raw data is large, the significant degrees of freedom are only six in space and two among the six properties; temperature and salinity are represented by one mode, while the nutrients by another. The modal structure reflects some underlying simplicity in ocean physics. EOFs form a quantitative basis from which models of the ocean's hydrographic structure can be constructed for various degrees of complexities. As for the second aspect, two applications of inverse methods are explored on small regional scales. The first problem addressed concerns the circulation inside a 12° square located in the eastern basin over the axis of the Mediterranean Water tongue. The study is based on an ocean model constructed by mapping the modes identified in the first half of the thesis over the entire North Atlantic Ocean. A combination of box model inverse and β-spiral method is used to determine the geostrophic reference level velocities. The circulation consists of an anticyclonic circulation near the surface, which is part of the eastern half of the wind-driven subtropical gyre. The flow at depth is weak, and is a cyclonic circulation around the core of the Mediterranean Water tongue. In the second inverse problem, we examine a decaying warm-core ring. Observations of a warm-core ring are used to formulate a model for diagnosing the physics of ring change over a two month period. About 30 hydrographic casts and acoustic doppler current measurements are used to generate estimates of an equivalent radially symmetric ring with radial contrasts of stratification, temperature, salinity, azimuthal velocity, angular momentum, and potential vorticity. A series of related models are inverted for the ring circulation and mixing coefficients. The circulation is insensitive to the model details, is well-resolved, and is a radial outflow and upwelling. Eddy coefficients are only partially resolved; determining the mixing with any degree of confidence appears to require a much more elaborate data set than the one available.

Description: This research was funded in part by the Office of Naval Research (Secretary of the Navy Chair) and the National Science Foundation under grant OCE 85-21685.

Description: Author Posting. © American Meteorological Society, 2012. 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 42 (2012): 855–868, doi:10.1175/JPO-D-10-05010.1.

Description: Data from the Hudson River estuary demonstrate that the tidal variations in vertical salinity stratification are not consistent with the patterns associated with along-channel tidal straining. These observations result from three additional processes not accounted for in the traditional tidal straining model: 1) along-channel and 2) lateral advection of horizontal gradients in the vertical salinity gradient and 3) tidal asymmetries in the strength of vertical mixing. As a result, cross-sectionally averaged values of the vertical salinity gradient are shown to increase during the flood tide and decrease during the ebb. Only over a limited portion of the cross section does the observed stratification increase during the ebb and decrease during the flood. These observations highlight the three-dimensional nature of estuarine flows and demonstrate that lateral circulation provides an alternate mechanism that allows for the exchange of materials between surface and bottom waters, even when direct turbulent mixing through the pycnocline is prohibited by strong stratification.

Description: The funding for this research was obtained from NSF Grant OCE-08-25226.

Description: 2012-11-01

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 1989

Description: Eighteen months of sea surface height data from the GEOSAT altimeter along collinear subtracks were analyzed for information on the circulation pattern in the Bering Sea. Seventy subtracks from both ascending and descending orbits, with as many as 35 repeat cycles along each subtrack, were analyzed. Orbit errors were removed from the height data using a least-squares fit to a cubic polynomial, weighted by the inverse of the height variance. Addition of the weights decreased contamination of residual height profiles by the large geoid signal. Composite maps of variability along each track revealed patterns of increased variability in the regions of the documented Bering slope current (BSC) and the proposed western boundary current (WBC); however, no evidence was found of the expected bifurcation of the BSC near the Siberian coast. Past observations of tides in the Bering Sea were reviewed along with a local tide model to detect tidal contributions to the mesoscale sea surface height variability. The tidal analysis suggested that residual tides contributed primarily to the longer wavelengths which were removed in the collinear processing. Examination of the Schwiderski tidal correction proved it to be a sensible correction, reducing the height variance by approximately 60%. Finally, using a Gaussian model for the BSC velocity profile, synthetic residual heights were generated and fit to the actual data to produce estimates of absolute surface geostrophic velocity and transport. Comparisons of mean flow, height fluctuations and seasonal trends across the BSC, the WBC and Bering Strait support the hypothesis that the BSC turns north at Cape Navarin into the WBC which, in turn, is capable of supplying a major part of the transport through the Bering Strait.

Description: Submitted in partial fulfillment of the requirements for the degree of Ocean Engineer at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution August 1993

Description: Travel time perturbations of adiabatic normal modes due to an internal tide and internal mode field in the Barents Sea are examined. A formalism for the travel time perturbation due to a change in sound speed is presented. Internal tide and internal wave amplitude spectra are calculated from Brancker temperature loggers which were deployed on moorings in the Barents Sea during the August 1992 Barents Sea Polar Front Experiment. In particular, the first three internal wave mode amplitudes are estimated from the four Brancker temperature loggers on the southwest mooring of the array. Modal perturbations in acoustic pulse travel time and the travel time covariance are calculated and compared for consistency to a simple ray model. These perturbations are small for the modal arrivals that the vertical acoustic array which was deployed is expected to resolve. The third internal wave mode has the largest impact on the acoustic arrivals, per unit amplitude, but the first internal wave mode dominates the scattering due to having a much larger amplitude overall.

Description: 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 2012

Description: This thesis focuses on ocean circulation and atmospheric forcing in the Atlantic Ocean at the Last Glacial Maximum (LGM, 18-21 thousand years before present). Relative to the pre-industrial climate, LGM atmospheric CO2 concentrations were about 90 ppm lower, ice sheets were much more extensive, and many regions experienced significantly colder temperatures. In this thesis a novel approach to dynamical reconstruction is applied to make estimates of LGM Atlantic Ocean state that are consistent with these proxy records and with known ocean dynamics. Ocean dynamics are described with the MIT General Circulation Model in an Atlantic configuration extending from 35°S to 75°N at 1° resolution. Six LGM proxy types are used to constrain the model: four compilations of near sea surface temperatures from the MARGO project, as well as benthic isotope records of δ18O and δ13C compiled byMarchal and Curry; 629 individual proxy records are used. To improve the fit of the model to the data, a least-squares fit is computed using an algorithm based on the model adjoint (the Lagrange multiplier methodology). The adjoint is used to compute improvements to uncertain initial and boundary conditions (the control variables). As compared to previous model-data syntheses of LGM ocean state, this thesis uses a significantly more realistic model of oceanic physics, and is the first to incorporate such a large number and diversity of proxy records. A major finding is that it is possible to find an ocean state that is consistent with all six LGM proxy compilations and with known ocean dynamics, given reasonable uncertainty estimates. Only relatively modest shifts from modern atmospheric forcing are required to fit the LGM data. The estimates presented herein successfully reproduce regional shifts in conditions at the LGM that have been inferred from proxy records, but which have not been captured in the best available LGM coupled model simulations. In addition, LGM benthic δ18O and δ13C records are shown to be consistent with a shallow but robust Atlantic meridional overturning cell, although other circulations cannot be excluded.

Description: Primary support was provided by a National Defense Science and Engineering Graduate Fellowship and two National Science Foundation awards: Award #OCE-0645936: “Beyond the Instrumental Record: the Case of Circulation at the Last Glacial Maximum” and Award #OCE-1060735: “Collaborative Research: Beyond the Instrumental Record - the Ocean Circulation at the Last Glacial Maximum and the de-Glacial Sequence”. Important secondary support came from the National Ocean Partnership Program and the National Aeronautics and Space Administration via the ECCO effort at MIT.

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

Description: Ocean modellers seek to understand the circulation of the oceans, or portions thereof, by developing models of the ocean they can solve. This tractability constraint forces ocean modellers to make choices. Naturally, they hope to make intelligent choices, but whenever a new model is being developed or an existing one extended, the issue of tractability lurks. The large-scale, basin-wide, circulation of the oceans can be divided into two components, classified by their driving force. The wind-driven circulation, whose flow occurs mainly above the thermocline, was first explained qualitatively by Stommel (1948) with a simple, elegant analytical model. The other component of the oceans' circulation, the density-driven, or thermohaline circulation, flows below the thermocline. Again, the first simple analytical model for the deep thermohaline flow was proposed by Stommel (1958) and developed by Stommel and Aarons (1959) whose basic ideas underlie even the most recent conceptual models of the large-scale circulation. The details of the thermohaline circulation and its interaction with the wind-driven circulation in a realistic ocean basin is a problem which is not tractable analytically. This has driven ocean modellers interested in this aspect of the oceans' circulation to numerical models: ocean circulation models.

Description: 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 1994

Description: This work investigates whether large-scale coherent vortex structures driven by wave-current interaction (Langmuir circulation) are responsible for maintaining the oceanic mixed layer. Langmuir circulations dominate the near-surface vertical transport of momentum and density when the characteristic scale for forcing (defined as the Craik-Leibovich instability parameter γCLS) is stronger than the characteristic scale for diffusive decay γdiff. Since the wave-current forcing is concentrated near the surface both terms depend on the cell geometry. Cells with long wavelengths penetrate more deeply into the water column. These cells grow more slowly than the fastest growing mode for most cases, but always dominate the solution in the absence of Coriolis forces. In the presence of Coriolis forces, the horizontal wavelength and thus the depth of penetration are limited. When a cell geometry is found such that γCLS » γdiff, the current profile produced by small-scale diffusion is unstable to Langmuir cells and the cells replace small-scale diffusion as the dominant vertical transport mechanism for momentum and density. The perturbation crosscell shear is predicted to scale as γCLS. Such a scaling is observed during two field experiments. The observed velocity profile during these experiments is more sheared than predicted by a model which implicitly assumes instantaneous mixing by large eddies, but less sheared than predicted by a model which assumes small-scale mixing by near-isotropic turbulence. The latter profile is unstable to Langmuir cells when waves are present. The inclusion of cells driven by wave-current interaction explains the failure of the mixed layer to restratify on two days with high waves and low wind. Wave-current interaction introduces a small but efficient source of energy for transporting density which goes as the surface stress times the Stokes drift.

Description: The Office of Naval Research supported me throughout graduate school, first as an ONR Graduate Fellow. and later as a research assistant under the Surface Waves Processes Program (ONR Grant N00014-90-J-1495).

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and Woods Hole Oceanographic Institution, February 2012

Description: Trace metal cycling is one of many processes that influence ocean ecosystem dynamics. Cobalt, iron, and manganese are redox active trace metal micronutrients with oceanic distributions that are influenced by both biological and abiotic sources and sinks. Their open ocean concentrations range from picomolar to nanomolar, and their bioavailabilities can impact primary production. Understanding the biogeochemical cycling of these hybrid-type metals with an emphasis on cobalt was the focus of this thesis. This was accomplished by determining the dissolved distributions of these metals in oceanic regions that were characterized by different dominant biogeochemistries. A large subsurface plume of dissolved cobalt, iron, and manganese was found in the Eastern South Atlantic. The cause of this plume is a combination of reductive dissolution in coastal sediments, wind-driven upwelling, advection, biological uptake, and remineralization. Additional processes that are discussed as sources of metals to the regions studied during this thesis include isopycnal uplift within cold-core eddies (Hawaii), ice melt (McMurdo Sound, Antarctica), riverine input (Arctic Ocean), and winter mixing (McMurdo Sound). The biological influence on surface ocean distributions of cobalt was apparent by the observation of linear relationships between cobalt and phosphate in mid to low latitudes. The cobalt:phosphate ratios derived from these correlations changed over orders of magnitude, revealing dynamic variability in the utilization, demand, and sources of this micronutrient. Speciation studies suggest that there may be two classes of cobalt binding ligands, and that organic complexation plays an important role in preventing scavenging of cobalt in the ocean. These datasets provided a basis for comparing the biogeochemical cycles of cobalt, iron, and manganese in three oceanic regimes (Hawaii, South Atlantic, McMurdo Sound). The relative rates of scavenging for these metals show environmental variability: in the South Atlantic, cobalt, iron, and manganese were scavenged at very different rates, but in the Ross Sea, mixing and circulation over the shallow sea was fast, scavenging played a minor role, and the cycles of all three metals were coupled. Studying the distributions of these metals in biogeochemically distinct regions is a step toward a better understanding of their oceanic cycles.

Description: Funding for this research was provided by the the National Science Foundation Chemical Oceanography (Division of Ocean Sciences OCE-0452883, OCE-0752291, OCE-0928414, OCE-0732665, OCE-0440840, OCE-0327225), the Center for Microbial Research and Education, the WHOI Coastal Ocean Institute, and the WHOI Ocean Life Institute, WHOI Academic Programs Office, and a Fye Teaching Fellowship.

Description: 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 May 1989

Description: The relationship between depth-averaged velocity and bottom stress for wind-driven flow in unstratified coastal waters is examined here. The adequacy of traditional linear and quadratic drag laws is addressed by comparison with a 2 1/2-D model. A 2 1/2-D model is one in which a simplified 1-D depth-resolving model (DRM) is used to provide an estimate of the relationship between the flow and bottom stress at each grid point of a depth-averaged model (DAM). Bottom stress information is passed from the DRM to the DAM in the form of drag tensor with two components: one which scales the flow and one which rotates it. This eliminates the problem of traditional drag laws requiring the flow and bottom stress to be collinear. In addition , the drag tensor field is updated periodically so that the relationship between the velocity and bottom stress can be time-dependent. However, simplifications in the 2 1/2-D model that render it computationally efficient also impose restrictions on the time-scale of resolvable processes. Basically, they must be much longer than the vertical diffusion time scale. Four progressively more complicated scenarios are investigated. The important factors governing the importance of bottom friction in each are found to be 1) non-dimensional surface Ekman depth, u.5/fh where u.s is the surface shear velocity, f is the Coriolis parameter and h is the water depth 2) the non-dimensional bottom roughness, zo/h where zo is the roughness length and 3) the angle between the wind stress and the shoreline. Each has significant influence on the drag law. The drag tensor magnitude, r, and the drag sensor angle, θ are functions of all three, while a drag tensor which scales with the square of the depth-averaged velocity has a magnitude, Cd, that only depends on zo/h. The choice of drag Jaw is found to significantly affect the response of a domain. Spin up times and phase relationships vary between models. In general, the 2 1/2-D model responds more quickly than either a constant r or constant Cd model. Steady-state responses are also affected. The two most significant results are that failure to account for θ in the drag law sometimes leads to substantial errors in estimating the sea surface height and to extremely poor resolution of cross-shore bottom stress. The latter implies that cross-shore near-bottom transport is essentially neglected by traditional DAMs.

Description: Financial support during my time in graduate school came from the Woods Hole Oceanographic Institution and grants from the National Science Foundation (OCE84-03249) and the Office of Naval Research (N00014-86-K-0061).

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 1989

Description: This thesis studies mixing and convection in a rectangular basin driven by a specified heat flux at the surface. A numerical model is constructed for this purpose. The main focus of the study is on the density and circulation structure resulting from the thermal forcing. In chapter two, a simple vertical one-dimensional model is developed to examine the mixing processes under a given surface heat flux. In order to simulate strong vertical mixing in the region where stratification is unstable, turbulent processes are modeled by a convective overturning parameterization of eddy viscosity and diffusivity. The results show that the density structure is strongly affected by the convective overturning adjustment as surface cooling prevails, and the resulting density field is nearly depth independent. In chapter three, a more complicated two-dimensional model is constructed to simulate mixing and circulation in a vertical rectangular basin with rigid boundaries. The aspect ratio of the basin ranges from 1 to 0.001 and Rayleigh number from 104 to 2 x 1012. It is found that the circulation pattern is dominated by these two important numbers. The roles of density overturning and density-momentum overturning mixing are further investigated. The results show that the convective overturning not only homogenizes the density field in the unstably stratified region but also contributes to increase the circulation. A crude scale analysis of the system shows that the characteristics of the density and momentum fields from the analysis agree well with the numerical results.

Description: 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 1993

Description: An extensive set of new high-quality hydrographic data is assembled in order to determine the mean circulation in the equatorial Pacific, and thus the pathways for cross-equatorial and cross-gyre exchange. Making up the core of the data set are two onetime transpacific zonal sections nominally at 10°N and 14°S. Supplementing these are repeat surveys of the equatorial currents along the 165°E meridian with direct shear measurements, and repeat surveys of the western boundary current at 8°N including direct velocity measurements. The repeat survey data are crucial for obtaining a good estimate of the mean conditions in the face of strong annual and interannual variability of the near-equatorial flow field. A comparison with historical XBT and hydrographic data shows that the interior thermocline transports in the one-time sections are fortuitously representative of the mean conditions. A detailed study of the water mass distribution along the sections is the basis for choosing reference levels for the thermal wind shear in an initial guess circulation field. Using an inverse model, the initial guess circulation is adjusted such that volume, heat and salt arc conserved in a set of subthermocline layers (δΘ 〉 26.7). Cross-isopycnal diffusion and advection are explicitly accounted for in the inverse model, and the diapycnal diffusivity is constrained to be positive, though its value is allowed to vary with depth and location. Net mass conservation constraints are applied to the enclosed volumes of the North Pacific and eastern Pacific, and essentially require that the Ekman divergence be equal to the geostrophic convergence. The Ekman fluxes as estimated from wind-stress climatologies are an important element of the mass budget, and yet are subject to large uncertainties. The model is therefore given the freedom to determine the Ekman fluxes within the range of error of the wind-stresses. The circulation of the coldest waters (Θ 〈 1.2°C) is dominated by the northward flow of Lower Circumpolar Water (LCPW) in a system of narrow western boundary currents. A net transport of 12.1 Sv of LCPW flows across 14°S, 9.6 Sv of which flows into the North Pacific across 10°N. The bulk of the LCPW flux across the equator appears to occur in the denser part of the western boundary current which follows topography directly across the equator. Dissipation in the boundary layer can thus modify the potential vorticity of the fluid and allow it to cross the equator. The circulation of the upper part of the LCPW is dominated by a strong westward jet at the equator which is supplied both by upwelling from below and the recirculation of modified LCPW from the North Pacific. At mid-depth (4.0 〉 Θ 〉 1.2°C) high silica and low oxygen concentrations mark the North Pacific Deep Water (NPDW) which is present in both the North and South Pacific Oceans. Across both 10°N and 14°S, a net of 11 Sv of NPDW flows southward, returning the northward mass flux associated with the LCPW. In contrast to the LCPW, narrow western boundary currents are not present in this layer, and it is not clear how the deep water flows across the equator. Strong zonal jets on and about the equator may be important in allowing mass to cross the equator by increasing the time available for the cross-equatorial diffusion of potential vorticity to act on a fluid parcel. At intermediate depths equatorward advection is suggested by the presence of intermediate water salinity minima formed in the subpolar latitudes: Antarctic Intermediate Water dominates the 4 to 8°C classes south of the equator, while North Pacific Intermediate Water occupies this range north of the equator. Determination of the mean circulation of the intermediate waters is, however, confounded by the large eddies that dominate the geostrophic transport stream function along the onetime zonal sections. The equatorial thermocline is occupied by waters of subtropical origin: the shallow salinity minimum waters and saline Central Water from both the North and South Pacific Ocean. The equator marks the location of a front between northern and southern subtropical gyre waters, except in the lower thermocline where water from the South Pacific subtropical gyre penetrates to about 4°N to feed the Northern Subsurface Countercurrent at 165°E. All of the equatorward flowing thermocline waters are entrained in the eastward equatorial currents which in turn feed the upwelling system in the eastern Pacific. The upwelled waters largely supply the South Equatorial Current in the eastern Pacific, accounting for its large transport compared to that predicted by Sverdrup dynamics. Northward flow across the equator of the upwelled waters in the thermocline or surface layer in the western Pacific is necessary to supply the Ekman flux into the North Pacific. The analysis indicates that the Pacific Ocean does not convert a large amount of abyssal water to thermocline water, as required by several theories of the global thermohaline circulation. In contrast to the Atlantic Ocean, the thermocline circulation in the Pacific appears decoupled from the abyssal overturning, with little upwelling of abyssal waters occurring in either the North Pacific or the equatorial Pacific. The leakage of Pacific water into the Indian Ocean is deduced to be essentially zero, though an error analysis allows a range of 0-8 x 106m3s-1.

Description: I was supported by the 1986 Caltex Graduate Women Scholarship, and a NASA Scholarship in Global Change Research.

Description: 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 1994

Description: Paleo-tracers such as carbon 13 and cadmium show that the deep Atlantic was enriched in nutrients during the Last Ice Age. The conventionally accepted interpretation of these higher nutrient levels is that a reduction of the rate of formation of nutrient-depleted Lower North Atlantic Deep Water (Lower NADW) allowed nutrient-rich Antarctic Bottom Water (AABW) to push further north during the Last Glacial Maximum (LGM) (Boyle and Keigwin, 1982; 1987; Duplessy et al., 1988). The evidence for this interpretation is re-examined in this work, with an emphasis on the quantitative analysis of the paleo-data. An end-member analysis of the δ13C data indicates a larger volume of AABW and a smaller volume of Lower NADW during the LGM. It is not yet possible, however, to quantify the extent of the volume differences between the modern and the glacial distributions, because the LGM δ13C end-members are poorly known. The second issue examined in this thesis deals with the interpretation of the water mass distribution, inferred from paleo-tracers, in terms of the oceanic circulation. Using a dynamical inverse model of the North Atlantic and a kinematic inverse model of the South Atlantic, it is shown that a tracer distribution corresponding to a significantly reduced volume of Lower NADW does not necessarily correspond to a reduced flux of NADW. Indeed, a circulation almost identical to a modem ocean reference circulation is consistent with the available LGM δ13C and δ18O data A flux of Lower NADW reduced by 50%, though not needed to explain the LGM tracer distribution, is also consistent with the data Thus, the paleo-tracers δ13C and δ18O do not suffice to quantify the flux of NADW in the glacial ocean. The modem ocean circulation is one of many possible circulations consistent with the available δ13C and δ18O data.

Description: This research was funded by the National Science Foundation under grant OCE-9205942.

Description: 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 October 1993

Description: In this thesis, production of dense water that feeds the dense overflows across the Greenland-Scotland Ridge has been considered. A new circulation scheme is developed which is consistent with the water masses, currents and air-sea fluxes in the region, and with the important observation that the dense overflows show little or no seasonal or interannual variability. An inverse box model has been constructed that shows that the new circulation scheme is consistent with conservation statements for mass, heat and salt as well. According to the new circulation scheme the major buoyancy is lost in the North Atlantic Current, which enters the Norwegian Sea between Iceland and Scotland, and flows northward towards the Arctic Ocean and the Barents Sea. The transformation is due to a large net annual heat loss over the North Atlantic Current, combined with a long residence time (2-3 years) and a large surface area. After subduction, one branch of the North Atlantic Current enters the Arctic Ocean, is modified in hydrographic properties into those associated with the Denmark Strait Overflow Waters in the western North Atlantic, exits the Arctic Ocean in the western Fram Strait and flows with the East Greenland Current towards the Denmark Strait Another branch of the North Atlantic Current recirculates directly in the Fram Strait and flows towards the Denmark Strait with the East Greenland Current This branch will not sink to the bottom of the North Atlantic as it is less compressible than the Arctic branch. The third branch of the North Atlantic Current enters the Barents Sea, continues to lose buoyancy, and enters the Arctic Ocean at intermediate depth. This branch exits the Arctic Ocean in the western Fram Strait, circulates around the Greenland Sea, enters the Norwegian Sea, and flows towards the Frer¢-Shetland Channel. The traditional view holds that the major sources of the dense overflows are the Iceland and Greenland gyres, west of the North Atlantic Current. Aside from the finding that the new circulation scheme is more likely in terms of water mass properties, currents etc., one fundamental problem with the old scheme lies with supplying a substantial overflow. There are indications that the production of dense water in the gyres is sensitive to the highly variable surface conditions and that indeed the production tends to shut on and off. The reservoirs in the gyres are so small that they would be drained within a few years if they were to supply the overflows during a shutdown period. Production of dense water within the North Atlantic Current is less sensitive to surface conditions. The density in the gyres is gained at a temperature around freezing, whereas in the North Atlantic Current the density is gained well above freezing. Therefore a freshwater anomaly in the two domains will have different consequences for vertical · overturning: within the North Atlantic Current the freshening can be overcome by further cooling, whereas in the gyres freezing will occur and the vertical overturning will cease. The observed lack of a significant seasonal signal associated with the dense overflows is consistent with the new circulations scheme. The net annual cooling dominates the seasonal oscillation in the atmospheric heat loss for time scales comparable with the residence time of the Atlantic Water within the domain. Thus winter formation of dense water within the North Atlantic Current does not induce a seasonal signal in the transport field of the dense water.

Description: Funding for this work was partly provided by a NASA Global Change Fellowship.

Description: 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 July 1993

Description: In this thesis, the dynamic role of bottom topography in a β-plane channel is systematically studied in both linear homogeneous and stratified layer models in the presence of either wind stress (Chapters 2, 3, 4, and 6) or buoyancy forcing (Chapter 5). In these studies, the structure of the geostrophic contour plays a fundamental role, and the role of bottom topography is looked at from two different angles. It is shown that blocking all the geostrophic contours leads to two different physical processes in which bottom topographic form drag is generated (Chapters 2, 3 and 4) and enables geostrophic flow in a β-plane channel to support a net cross-channel volume transport (Chapters 5 and 6). It is demonstrated that by blocking all the geostrophic contours in the presence of a sufficiently high ridge, the dynamics of both source-sink and wind driven circulations in a β-plane is similar to that in a dosed basin. First, wind-driven circulation in the inviscid limit is discussed in a linear barotropic channel model in the presence of a bottom ridge. There is a critical height of the ridge, above which all geostrophic contours in the channel are blocked. In the subcritical case, the Sverdrupian balance does not apply and there is no solution in the inviscid limit. In the supercritical case, however, the Sverdrupian balance applies. The form drag is generated through two different physical processes: the through-channel recirculating flow and the Sverdrupian gyre flow. These processes are fundamentally different from the nonlinear Rossby wave drag generation. In this linear model, the presence of a supercritical high ridge is essential in the inviscid limit. With this form drag generation determined, an explicit form for the zonal transport in the channel is obtained, which shows what model parameters determine the through-channel transport. In addition, the model demonstrates that most of the potential vorticity dissipation occurs at the northern boundary where the ridge intersects. The result from the homogeneous channel model in Chapter 2 is then extended to a model whose geometry consists of a zonal channel and two partial meridional barriers along each boundary at the same longitude. Both the model transport and especially the model circulation are significantly affected by the presence of the two meridional barriers. The presence of the northern barrier always leads to a decrease in the transport. The presence of the southern barrier, however, increases the transport for a narrow ridge. The northern barrier only has a localized influence on the circulation pattern, while the southern barrier has a global influence in the channel. Then a multi-layer Q-G model is constructed by assuming that potential vorticity in all subsurface layers is homogenized. The circulation is made up of baroclinic and the barotropic part. The barotropic part is same as that in a corresponding barotropic model, and is solely determined by the wind stress, while the baroclinic part is not directly related to the wind stress. It is determined by the potential vorticity homogenization and lateral boundary conditions. The presence of the stratification does not affect the bottom topographic form drag generation. The interfacial form drag is generated by the stationary eddies. Corresponding to the circulation structure, the zonal through-channel transport associated with the barotropic circulation is determined by the wind stress and bottom topography. The other part associated with the baroclinic circulation, however, is not directly related to the wind stress and it is determined by the background stratification. Based upon the discussion on the geostrophic contour, a simple barotropic model of abyssal circulation in a circumpolar ocean basin is constructed. The presence of a supercritically high ridge is both necessary and sufficient for geostrophic flow in a β-plane channel to support a net cross-channel volume flux. In the presence of a sufficiently high ridge, the classical Stommel & Arons theory applies here, but with significant modifications. The major novelty is that a throughchannel recirculation is generated. Both its strength and direction depend critically upon the model parameters. Then, a schematic picture of the abyssal circulation in a rather idealized Southern Ocean is obtained. The most significant feature is the narrow current along the northern boundary of the circumpolar basin, which feeds the deep western boundary currents of the Indian Ocean and Pacific Ocean and connects all the oceanic basins in the Southern Ocean. Finally, the question of how the northward surface Ekman transport out of the circumpolar ocean is returned is discussed in a two-layer model with an infinitesimally thin surface Ekman layer on top of a homogeneous layer of water in a rather idealized Southern Ocean basin. First, the case with a single subtropical ocean basin is discussed. In the case with a sufficiently high ridge connecting the Antarctic and the meridional barrier, an explicit solution is found. The surface Ekman layer sucks water from the lower layer in the circumpolar basin. This same amount of water flows northward as the surface Ekman drift. It downwells in the subtropical gyre, and is carried to the western boundary layer. From there, the same amount of water flows southward as a western boundary current across the inter-gyre boundary between the circumpolar ocean and the subtropical gyre along the west coast to the southern boundary of the meridional barrier. Then, the same amount of water is carried southward and feeds the water loss to the surface Ekman layer due to the Ekman sucking in the interior circumpolar ocean. The case with multiple subtropical ocean basins such as the Southern Ocean is also discussed. It is demonstrated that the surface Ekman drift drives a strong inter-basin water mass exchange.

Description: This thesis was supported by National Science Foundation through grant OCE OCE90-17158. Part of the numerical simulation was performed at NCAR's supercomputer, and was supported by SCD/NCAR.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and Woods Hole Oceanographic Institution February 1994

Description: In this work we study motion of a baroclinic upper-ocean eddy over a large-scale topography which simulates a continental slope. We use a quasigeostrophic f-plane approximation with continuous stratification. To study this problem we develop a new numerical technique which we call "semi-lagrangian contour dynamics". This technique resembles the traditional 2-D contour dynamics method but differs significantly from it in the numerical algorithm. In addition to "Lagrangian" moving contours it includes an underlying "Eulerian" regular grid to which vorticity or density fields are interpolated. To study topographic interactions in a continuously stratified model we use density contours at the bottom in a similar manner as vorticity contours are used in the standard contour dynamics. For the case of a localized upper-ocean vortex moving over a sloping bottom the problem becomes computationally 2-dimensional (we need to follow only bottom density contours and the position of the vortex itself) although the physical domain is still 3-dimensional. Results of the numerical model lndicate importance of baroclinic effects in the vortex-topography interaction. After the initial surge of topographic Rossby waves a vortex moves almost steadily due to the interaction with a bottom density anomaly which is created and supported by a vortex itself. This anomaly is equivalent to a region of opposite-signed vorticity with a total circulation exactly compensating that of a vortex. This results in a vertically aligned dipolar structure with the total barotropic component equal to zero. Analytical considerations explaining this effect are presented and formulated in a more general statement which resembles but does not coincide with the "zero angular momentum theorem" of Flierl, Stern and Whitehead, 1983. In such steady translation the centroid of a bottom density anomaly is displaced horizon tally from the center of an upper-ocean vortex so the whole system moves due to this misalignment, which is known as a "he tonic mechanism". Cyclonic vortices go generally upslope, and anticyclones - in a downslope direction. The along-slope component of their motion depends upon the strength of a vortex, curvature of the bottom slope and background flows. When surrounded by a bowl-shaped topography anticyclonic vortices tend to stay near the deepest center of a basin, even resisting ambient flows which advect them outward. Application of this results to various oceanic examples (particularly to the "Shikmona eddy" in the Eastern Meditenanian) is discussed. Our results show that the behavior of a vortex over a sloping bottom differs significantly from its motion on the planetary beta-plane (but with a flat bottom). To explain this difference we introduce the concept of a "wave-breaking regime" relevant for the case of a planetary beta-effect, and a "wave-gliding regime" which characterizes the interaction of an eddy with a topographic slope.

Description: This work was supported by the NSF grant #OCE 90-12821.

Description: 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 February 2013

Description: Efforts to monitor the ocean for signs of climate change are hampered by ever-present noise, in the form of stochastic ocean variability, and detailed knowledge of the character of this noise is necessary for estimating the significance of apparent trends. Typically, uncertainty estimates are made by a variety of ad hoc methods, often based on numerical model results or the variability of the data set being analyzed. We provide a systematic approach based on the four-dimensional frequency-wavenumber spectrum of low-frequency ocean variability. This thesis presents an empirical model of the spectrum of ocean variability for periods between about 20 days and 15 years and wavelengths of about 200{10,000 km, and describes applications to ocean circulation trend detection, observing system design, and satellite data processing. The horizontal wavenumber-frequency part of the model spectrum is based on satellite altimetry, current meter data, moored temperature records, and shipboard ADCP data. The spectrum is dominated by motions along a "nondispersive line". The observations considered are consistent with a universal ω-2 power law at the high end of the frequency range, but inconsistent with a universal wavenumber power law. The model spectrum is globally varying and accounts for changes in dominant phase speed, period, and wavelength with location. The vertical structure of the model spectrum is based on numerical model results, current meter data, and theoretical considerations. We find that the vertical structure of kinetic energy is surface intensified relative to the simplest theoretical predictions. We present a theory for the interaction of linear Rossby waves with rough topography; rough topography can explain both the observed phase speeds and vertical structure of variability. The improved description of low-frequency ocean variability presented here will serve as a useful tool for future oceanographic studies.

Description: This research was supported by NASA under grants NNG06GC28G and NNX08AR33G.

Description: 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 June 1995

Description: Moored time series from the Coastal Ocean Dynamics Experiment (CODE), Shelf Mixed Layer Experiment (SMILE), Sediment Transport Events over the Shelf and Slope (STRESS) study, and Northern California Coastal Circulation Study (NCCCS) are used to study subtidal cross-shelf circulation over the northern California shelf. The northern California shelf, like much of the United States Pacific coast, is subject to strong wind forcing which exhibits characteristic seasonality. In winter and early spring, it is distinguished by poleward and equatorward fluctuations on time scales of days and by weak monthly means. In summer, it is distinguished by periods of equatorward stress lasting several weeks and by relatively strong monthly means. The intensive winter and spring SMILE and STRESS and summer CODE-2 field programs permit the examination of cross-shelf circulation under both types of wind forcing conditions at a mid-shelf site (~90 m) 6 km from the northern California coast. The primary thesis goal is to examine the applicability of a two-dimensional conceptual model of wind-forced cross-shelf circulation. In this conceptual model, surface and bottom cross-shelf flows are forced by along-shelf wind stress and bottom stress, and interior cross-shelf flow compensates such that the depth-averaged flow is zero. A secondary thesis goal is to use the seasonal coverage of available field programs to gain insight into seasonal variability of cross-shelf circulation on the northern California shelf. To accomplish these goals, the observed subtidal cross-shelf circulation is examined in the context of the winter and spring heat and salt balances, an analytic model of wind-forced cross-shelf circulation, and the spatial scales of subtidal velocity. Mean and fluctuating heat and salt balances estimated between December, 1988 and May, 1989 demonstrate the importance of cross-shelf fluxes and their general consistency with the simple conceptual model. Mean fluxes are consistent with the weak mean equatorward wind stress observed during SMILE. The dominant terms in the fluctuating balances are the cross-shelf fluxes and local changes in heat and salt content. These are well correlated with each other and with the local along-shelf wind stress. The along-shelf heat flux divergence is of secondary importance to the fluctuating heat balance. It is uncorrelated with the along-shelf wind stress, and occurrences when it is strong are interpreted as effects of mesoscale features. To examine the applicability of the wind-forced conceptual model in more detail, a simple analytic model incorporating the assumptions of the conceptual model and observed local wind forcing is compared quantitatively to estimates of surface mixed layer, interior, and bottom mixed layer cross-shelf transport for winter SMILE and STRESS and summer CODE-2 observations. This comparison suggests the model is more suited to the transient wind forcing observed during SMILE and STRESS than to the steady wind forcing observed during CODE-2. For 2-3 day wind events between December, 1988 and February, 1989, the model is well correlated with observed depthdependent (total minus depth-averaged) transports throughout the water column and with total surface mixed layer transports. For 2-3 week wind events between April and July, 1982, the model does not work nearly as well below the surface mixed layer. In the absence of other processes, the locally wind-forced model implies that the wind stress sets the horizontal scales of subtidal velocity. Correlation scales estimated for subtidal along-shelf velocity over the northern California shelf are for all field programs longer than the maximum mooring separation (60 km) and are similar to those of the wind stress. However, along-shelf correlation scales of cross-shelf velocity are shorter than minimum mooring separations for CODE. SMILE and NCCCS time series do resolve along-shelf correlation scales for near surface cross-shelf velocity. During this time, along-shelf correlation scales for near surface cross-shelf velocity vary on a monthly time scale. They are generally long (30 km or more) when correlation with wind stress is high and short (15 km or less) when correlation with wind stress is low. On at least one occasion, short along-shelf correlation scales coincide with the intrusion of an offshore mesoscale feature onto the shelf. Results of the three studies show the two-dimensional model offers some insight into the observed subtidal cross-shelf circulation, particularly in winter. During this time, the heat balance, analytical transport model, and correlation scales all provide evidence that the winter wind-forced circulation is quasi-two-dimensional. Threedimensional variability on the shelf, though important on occasion, does not appear to be wind-driven and may result from the influence of offshore mesoscale features. A quite different story emerges for summer when the simple conceptual model of crossshelf circulation fails to describe adequately subsurface cross-shelf flow. Two useful areas of further investigation may be the non-linear response of cross-shelf velocity to wind forcing and its response to other processes such as remotely generated mesoscale features.

Description: In my first year at MIT, Carl Wunsch supported me through NSF grant OCE 88- 23043. At WHOI, I first started work on SMILE with Bob Beardsley under NSF grant OCE 87-16937 and continued working on it with Steve Lentz under NSF grant OCE 91-15713.

Description: 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 August 1995

Description: Data from fifteen globally distributed, modern, high resolution, hydrographic oceanic transects are combined in an inverse calculation using large scale box models. The models provide estimates of the global meridional heat and freshwater budgets and are used to examine the sensitivity of the global circulation, both inter and intra-basin exchange rates, to a variety of external constraints provided by estimates of Ekman, boundary current and throughflow transports. A solution is found which is consistent with both the model physics and the global data set, despite a twenty five year time span and a lack of seasonal consistency among the data. The overall pattern of the global circulation suggested by the models is similar to that proposed in previously published local studies and regional reviews. However, significant qualitative and quantitative differences exist. These differences are due both to the model definition and to the global nature of the data set. The picture of the global circulation which emerges from the models IS a complex, turbulent flow. When integrated across ocean basins not one, but two major cells emerge. The first connects an Atlantic overturning cell (estimated at 18± 4x 109 kg s- 1) to the Southern Ocean where the Antarctic Circumpolar Current carries lower deep waters to the Indian and Pacific basins where they are converted to upper deep and intermediate waters before returning to the Atlantic. The second cell connects the Pacific and Indian Basins to the north and south of Australia. In t his cell deep waters pass into the Pacific and return within the Indian Basin as intermediate waters after passing through the Indonesian Passages. The two cells are found to be independent of one another, i.e. within the models, the Indonesian Passages do not represent a significant element in a net global circulation. While there is ample evidence of westward flow around the southern tip of South Africa which would support a "warm" water path scenario, the variability of flow in this region, rich with eddies makes hydrography a poor estimator of the relative strengths of the controversial "warm" and "cold" water paths. All existing estimates of Indonesian Passage throughflow, including the smallest (O x 106 m3 s-1) and the largest (20 x 106 m3 s-1), are consistent with the model constraints. When the Pacific- Indian throughflow is not constrained, the model produces an estimate of 11 ± 14x 109 kg s-1. The model heat flux estimates are both significantly different from zero and quite robust to changes in initial assumptions, with the exception of the choice of wind field. Although in this work it was not possible to compute freshwater fluxes which were significantly different from zero, future inclusion of salinity anomaly constraints along with terms describing vertical diffusion may yet make it possible to compute significant freshwater :flux estimates from hydrography.

Description: This research was partially funded by a NASA Global Change Fellowship and was also supported by NASA under contract NAGW-1048 and NSF under contract OCE-9205942.

Description: 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 February 1995

Description: Inverse modeling activities in oceanography have recently been intensified, aided by the oncoming observational data stream of WOCE and the advance of computer power. However, interpretations of inverse model results from climatological hydrographic data are far from simple. This thesis examines the behavior of an inverse model in the WOCE CME (Community Modeling Effort) results where the physics and the parameter values are known. The ultimate hypotheses to be tested are whether the inferred circulations from a climatological hydrographic data set (where limited time means and spatial smoothing are usually used) represent the climatological ocean general circulations, and what the inferred "diffusion" coefficients really are. The inverse model is first tested in a non-eddy resolving numerical GCM ocean. Numerical/scale analyses are used to test whether the inverse model properly represents the GCM ocean. Experiments show how biased answers could result from an incorrect model, and how a correct model must produce the right answers. When the inverse model is applied to the time-mean hydrographic data of an eddy-resolving GCM ocean in the fine grid resolution of the GCM, the estimated horizontal circulation is statistically consistent with the EGCM time means in both patterns and values. Although the flow patterns are similar, the uncertainties for the GCM time means and the inverse model estimates are different. The former are very large, such that the GCM time-mean circulation has no significance in the deep ocean. The latter are much smaller, and with them the estimated circulations are well defined. This is consistent with the concept that ocean motions are very energetic, while variations of tracers (temperature, salinity) are low frequency. The inverse model succeeded in extracting the ocean general circulation from the "climatological" hydrographic data. The estimated vertical velocities are also statistically indistinguishable from the GCM time means. However, significant differences between the estimated "diffusion" coefficients and the EGCM eddy diffusion coefficients are found at certain locations. These discrepancies are attributed to the differences in physics of the inverse model and the EGCM ocean. The "diffusion" coefficients from the inversion parameterize not only the eddy fluxes, but also (part of) the temporal variation and biharmonic terms which are not explicitly included in the inverse model. Given the essentially red spectrum of the ocean, it makes sense to look for smooth solutions. Aliasing due to subsampling on a coarse grid and the effects of spatial smoothing are addressed in the last part of this thesis. It is shown that this aliasing could be greatly reduced by spatial smoothing. The estimated horizontal circulation from the spatially smoothed time-mean EGCM hydrographic data with a coarse grid resolution (2.4° longitude by 2.0° latitude) is generally consistent with the spatially smoothed EGCM time means. Significant differences only occur at some grid points at great depths, where the GCM circulations are very weak. The conclusions of this study are different from some previous studies. These discrepancies are explained in the concluding chapter. Finally, it should be pointed out that the issue of properly representing a GCM ocean by an inverse model is not identical to the issue of represent ing the real ocean by the same inverse model, since the GCM ocean is not identical to the real ocean. Numerical calculations show that both the non-eddy resolving and the eddy-resolving GCM oceans used in this work are evolving towards a statistical equilibrium. In the real ocean, the importance of temporal variation terms in the property conservation equations should also be analyzed when a steady mverse model is applied to a limited time-mean (the climatological) data set.

Description: This research was carried out under National Science Foundation grant OCE- 90-04396.

Description: Author Posting. © American Meteorological Society, 2012. 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 42 (2012): 291–305, doi:10.1175/JPO-D-11-043.1.

Description: A number of previous observational studies have found that the waters of the deep Pacific Ocean have an age, or elapsed time since contact with the surface, of 700–1000 yr. Numerical models suggest ages twice as old. Here, the authors present an inverse framework to determine the mean age and its upper and lower bounds given Global Ocean Data Analysis Project (GLODAP) radiocarbon observations, and they show that the potential range of ages increases with the number of constituents or sources that are included in the analysis. The inversion requires decomposing the World Ocean into source waters, which is obtained here using the total matrix intercomparison (TMI) method at up to 2° × 2° horizontal resolution with 11 113 surface sources. The authors find that the North Pacific at 2500-m depth can be no younger than 1100 yr old, which is older than some previous observational estimates. Accounting for the broadness of surface regions where waters originate leads to a reservoir-age correction of almost 100 yr smaller than would be estimated with a two or three water-mass decomposition and explains some of the discrepancy with previous observational studies. A best estimate of mean age is also presented using the mixing history along circulation pathways. Subject to the caveats that inference of the mixing history would benefit from further observations and that radiocarbon cannot rule out the presence of extremely old waters from exotic sources, the deep North Pacific waters are 1200–1500 yr old, which is more in line with existing numerical model results.

Description: GG is supported by the J. Lamar Worzel Assistant Scientist Fund and the Penzance Endowed Fund in Support of Assistant Scientists. PJH is supported by NSF Award 0960787.

Description: 2012-08-01

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 1991

Description: Based on the Levitus atlas, we find that the application of the Montgomery streamfunction to the isopycnal surfaces induces an error which can not be ignored in some regions in the ocean. The error arises from the sloping effect of the specific volume anomaly along isopycnal surfaces. By including the major part of this effect, new streamfunctions, namely the pressure anomaly and main pressure streamfunctions, are suggested for the use in potential density coordinates. By using the newly proposed streamfunction and by including the variations of specific volume anomaly along isopycnal surfaces, the inverse model proposed by Hogg (1987) is modified for increasing accuracy and applied to the Brazil Basin to study the circulation, diffusion and water mass balances. The equations in the model, i.e. the dynamic equation, continuity equation, integrated vorticity equation, and conservation equations for heat, salt and oxygen (in which a consumption sink term is allowed), are written in centered finite difference form with lateral steps of 2 degree latitude and longitude and 8 levels in the vertical. This system of equations with constraints of positive diffusivities and oxygen consumption rates is solved by the inverse method. The results indicate that the circulation in the upper oceans is consistent with previous works, but that in the deep ocean is quite different. In the NADW region, we find a coincidence of the flows with the tongues of water properties. The diffusivities and diapycnal velocities seem stronger in the region near the equator than in the south, with reasonable values. Diffusion plays an important role in the water mass balance. Examples show that similar property fields may results from different processes.

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution July 1991

Description: Six hydrographic basinwide sections, two in each of the three major ocean basins, are employed in a set of inverse calculations to determine the extent of exchange between the Pacific and Indian Oceans through the Indonesian Archipelago and the net global oceanic heat flux at 30°S. Using a model which combines the data for the South Pacific and South Indian Oceans, it is found that even the largest existing estimates of Indonesian Passage through flow (20 Sv) are consistent with the data. However, the available information cannot limit the extent of the exchange, i.e. both smaller and larger through flows produce physically reasonable circulation patterns. The seasonal and interannual variations which have been found by other investigators and which we are incapable of resolving, lead us to conclude that in the long term mean an estimate of ~10 Sv for the through flow is most reasonable. Globally, at 30°S, we find a net oceanic heat flux of -1.1 ± 1.7 PW, which is not significantly different from zero. It is dominated by a large (〉1 PW) southward heat flux in the Indian Ocean. Large equatorward (~0.8 PW) heat flux values in the South Atlantic Basin are not consistent with our data. We therefore conclude that although our data are consistent with some water following the warm water return path for NADW (Gordon 1986), the cold water path must play the dominant role in the maintenance of the global thermohaline cell associated with the formation process of NADW.

Description: 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 March 1992

Description: Rotating baroclinic and barotropic boundary currents flowing around a corner in the laboratory were studied in order to discover the circumstances under which eddies were produced at the corner. Such flows are reminiscent of oceanic coastal flows around capes. When the baroclinic currents, which consisted of surface flows bounded by a density front, encountered a sharp corner, immediately downstream of the corner an anticyclone grew in the surface layer for an angle of greater than 40 degrees. Varying the initial condition of the flow or the depth of the lower layer did not noticeably affect the gyre's properties except for its growth speed, which was greater when the lower layer was shallower. The barotropic currents were pumped along a sloping bottom, and also formed anticyclonic gyres which quickly attained an approximately steady state. For a given topography, the size of the gyre was proportional to the inertial radius u/f. Volume flux calculations based on the surface velocity revealed vertical shear which increased with gyre size. Hydraulic models were also applied to flow around gently curving topography to determine the critical separation curvature as a function of upstream parameters.

Description: This work was supported by the National Science Foundation grant OCE 89-15408.

Description: 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 April 1992

Description: A pair of simple models representing the interaction of a continuously stratified f-plane quasigeostrophic lens with a uniform external shear flow is examined. The study is motivated by the desire to understand the processes that affect Mediterranean Salt Lenses and other mesoscale lenses in the ocean. The first model represents the eddy as a pair of quasigeostrophic 'point potential vortices' in uniform external shear, where the two point vortices are imagined to represent the top and bottom of a baroclinic eddy. While highly idealized, the model succeeds in qualitatively reproducing many aspects of the behavior of more complex models. In the second model the eddy is represented by an isolated three dimensional patch characterized by quasigeostrophic potential vorticity linear in z, in a background flow with constant potential vorticity. The boundary of the lens may be deformed by interactions with a uniform background shear. A family of linearized analytical solutions representing such a vortex is discussed in Chapter 3. These solutions represent lens-like eddies with trapped fluid cores, which may propagate through the surrounding water when there is external vertical shear. The analysis predicts the possible forms of the boundary deformation in a specified external flow, and the precession rate of normal mode boundary perturbations in the absence of external flow. The translation speed of the lens with respect to the surrounding fluid is found to be a simple function of the external vertical shear and the core baroclinicity. A numerical algorithm which is a generalization of the contour dynamics technique to stratified quasigeostrophic flow is used to extend the linear results into the nonlinear regime. This numerical analysis allows a determination of the range of environmental conditions (e.g., the maximum shear and/ or core baroclinicity) in which coherent vortex solutions can be found, and allows the stability of the steadily translating solutions to be examined directly. It is found that the solutions are stable if neither the external shear nor the core baroclinicity is too large, and that the breakdown of the unstable solutions is characterized by the loss of an extrusion of core fluid to the surrounding waters. The translation speeds of the large amplitude numerical solutions are found to have the same functional dependence on the external vertical shear and the core baroclinicity that was found in the linear analysis, and it is demonstrated that the solutions translate at a rate which is equal to the background flow speed at the center of potential vorticity of the lens. As a test of the model results, new data from a recent SOFAR float experiment are presented and compared with the model predictions. The data show that the cores of two different Mediterranean Salt Lenses are tilted, presumably as a result of interactions with external flows. Both the sense of the tilt and its relation to the translation of the lens are in qualitative agreement with the model solutions.

Description: I am happy to acknowledge the support of the Office of Naval Research (grant N00014-89-J-1182) and the National Science Foundation (grants OCE 8916446 and OCE 87-00601).

Description: 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 1992

Description: This study concerns the barotropic interactions between a mesoscale eddy and a straight monotonic bottom topography. Through simple to relatively complicated modeling effort, some of the fundamental properties of the interaction are investigated. In chapter two, the fundamental aspects of the interaction are examined using a simple contour dynamics model. With the simplest model configuration of an ideal vortex and a step topography, the basic dynamical features of the observed oceanic eddy-topography interaction are qualitatively reproduced. The results consist of eddy-induced cross-topography exchange, formation of topographic eddies, eddy propagation and generation of topographic waves. In chapter three, a more complicated primitive equation model is used to investigate a mesoscale eddy interacting with an exponential continental shelf/slope topography on both f and β-planes. The f-plane model recasts the important features of chapter two. The roles of the eddy size and strength and the geometry of topography are studied. It is seen that the multiple anticyclonic eddy-slope interactions strongly affect the total cross-slope volume transport and the evolution of both the original anticyclone and the topographic eddy. Since a cyclone is trapped at the slope and eventually moves on to the slope, it is most effective in causing perturbation on the shelf and slope. The responses on the shelf and slope are mainly wavelike with dispersion relation obeying that of the free shelf-trapped wave modes. On the β-plane, the problem of an eddy colliding onto a continental shelf/slope from a distance with straight or oblique incident angles is investigated. It is found that the straight eddy incident is more effective in achieving large onslope eddy penetration distance than the oblique eddy incident. The formation of a dipole-like eddy pair consisting of the original anticyclone and the topographic cyclone acts to suppress the eddy decay due to long Rossby wave radiation. A weak along-slope current near the edge of the slope is found, which is part of a outer slope circulation cell originated from the Rossby wave wake trailing the propagating eddy. Model-observation comparisons in_chapter four show favorable qualitative agreement of the model results with some of the observed events in the eastern U.S. continental margins and in the Gulf of Mexico. The model results give dynamical interpretations to some observed features of the oceanic eddy-topography interactions and provide enlightening insight into the problem.

Description: Author Posting. © American Meteorological Society, 2012. 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 25 (2012): 1361–1389, doi:10.1175/JCLI-D-11-00091.1.

Description: The ocean component of the Community Climate System Model version 4 (CCSM4) is described, and its solutions from the twentieth-century (20C) simulations are documented in comparison with observations and those of CCSM3. The improvements to the ocean model physical processes include new parameterizations to represent previously missing physics and modifications of existing parameterizations to incorporate recent new developments. In comparison with CCSM3, the new solutions show some significant improvements that can be attributed to these model changes. These include a better equatorial current structure, a sharper thermocline, and elimination of the cold bias of the equatorial cold tongue all in the Pacific Ocean; reduced sea surface temperature (SST) and salinity biases along the North Atlantic Current path; and much smaller potential temperature and salinity biases in the near-surface Pacific Ocean. Other improvements include a global-mean SST that is more consistent with the present-day observations due to a different spinup procedure from that used in CCSM3. Despite these improvements, many of the biases present in CCSM3 still exist in CCSM4. A major concern continues to be the substantial heat content loss in the ocean during the preindustrial control simulation from which the 20C cases start. This heat loss largely reflects the top of the atmospheric model heat loss rate in the coupled system, and it essentially determines the abyssal ocean potential temperature biases in the 20C simulations. There is also a deep salty bias in all basins. As a result of this latter bias in the deep North Atlantic, the parameterized overflow waters cannot penetrate much deeper than in CCSM3.

Description: NCAR is sponsored by the National Science Foundation. The CCSM is also sponsored by the Department of Energy. SGY was supported by the NOAA Climate Program Office under Climate Variability and Predictability Program Grant NA09OAR4310163.

Description: 2012-09-01

Description: 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 November 1990

Description: One and two layer models are used to study flow over axisymmetric isolated topography. Inviscid or nearly inviscid flow in which non-linear effects have order one importance is considered, and both the effects of β and finite topography are included. A one-layer quasi-geostrophic model is used to find the shape of Taylor columns on both the f-plane and the β-plane in the inviscid limit of the frictional problem. In this limit, the boundary of the Taylor column is a streamline, and the velocity in both directions vanishes on the boundary. The fluid within the Taylor column is stagnant, corresponding to the solution that Ingersoll (1969) found for flow over a right circular cylinder on the f-plane. In this case, the Taylor column is circular. An iterative boundary integral technique is used to find the solutions for flow over a cone on the f-plane. In this case the Taylor column has a tear drop shape. Solutions are also found for flow on the β-plane over a cylinder, and the Taylor column is approximately elliptical for westward flow with the major axis in the x direction, while it is slightly elongated in the y direction for eastward flow. The stagnation point of the Taylor column is located on the edge of the topography for all the solutions found. It was not possible to find solutions for smooth topographic shapes. Steady solutions for flow over a right circular cylinder of finite height are studied when the quasi-geostrophic approximation no longer applies. The solution consists of two parts, one which is similar to the quasi-geostrophic solution and is driven by the potential vorticity anomaly over the topography and the other which is similar to the solution of potential flow around an cylinder and is driven by the matching conditions on the edge of the topography. When the effect of β is large, the transport over the topography is enhanced as the streamlines follow lines of constant background potential vorticity. For eastward flow, the Rossby wave drag can be much larger than predicted from quasi-geostrophic theory. A two-layer model over finite topography using the quasi-geostrophic approximation is developed. The topography is a right circular cylinder which goes all of the way through the lower layer and an order Rossby number amount into the upper layer, so that the quasi-geostrophic approximation can be applied consistently. This geometry allows description of flow in which an isopycnal intersects the topography. The model is valid for a different regime than existing models of steady flow over finite topography in a continuously stratified fluid in which the bottom boundary is an isopycnal surface. The solutions contain the two components that are found in the the barotropic model of flow over finite topography. The model breaks down when the interface goes above the topography which occurs more easily when the stratification is weak. Closed streamlines occur more readily over the topography when the stratification is weak, whereas in traditional quasi-geostrophic theory they occur more readily when the stratification is strong. Near the topography, the interface is depressed to the right and raised to the left (looking downstream). A hierarchy of time-dependent models is used to examine the initial value problem of flow initiation over topography on the f-plane. A modified contour dynamics method is developed that extends the range of problems to which contour dynamics can be applied. The method allows boundary and matching conditions to be applied on a circular boundary. A one-layer quasi-geostrophic model is used to show that more fluid that originates over the topography remains there when the flow is turned on slowly than when it is turned on quickly. Flow over finite topography in a one-layer model shows a variety of different behaviors depending on the topographic height. When the topography has moderate height, two cyclonic eddies are created; when the topography fills up most of the water column, the fluid oscillates on and off the topography as it moves around the topography in a clockwise direction, and none of the fluid is shed downstream. Two quasi-geostrophic stratified models are considered, one in which the topography is small, and the other in which it is finite. In the small topography model, an eddy is shed which is cyclonic, warm-core, and bottom-trapped. In contrast, the shed eddy is cyclonic, cold-core, and surfaceintensified in the finite depth model using the geometry described above.

Description: This study was supported by the Office of Naval Research Program N00014- 90-J-1839.

Description: 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 August 1991

Description: The oceanic distributions of tritium 3H), 3He, and the chlorofluorocarbons (CFCs) can be used to constrain the time-scales of the major ventilation pathways for an ocean basin such as the North Atlantic. I present a new global model function, developed from a factor analysis of the WMO/IAEA data set, for predicting the spatial and temporal variability of bomb-tritium in precipitation. Model estimates for the atmospheric 3H delivery to the North Atlantic are recomputed and combined with advective 3H input estimates in a budget for the North Atlantic Basin. Key features of the model budget include refined estimates of the 3H vapor flux and southward advection of 3H in the low salinity, surface flow from the Arctic. Arctic tritium sources contribute about half of the observed increase (40%) in the decay corrected tritium inventory from the 1972 GEOSECS program and the 1981 TTO/NAS program. The 3H concentration in the intermediate and deep waters for the sub-polar North Atlantic increased substantially between 1972 and 1981. A time dependent model for the 3H and 3He inflow to the abyssal Atlantic from the Nordic Seas is developed. The 3H and 3He distributions in the abyssal North Atlantic and Deep Western Boundary Current (DWBC) are also presented. A simple model of abyssal circulation is constructed using the model Nordic Seas overflow curves, the observed tracer gradients in the DWBC, and the GEOSECS and TTO tracer inventories for the deep basins. Although the tracer concentrations in the boundary current are rather insensitive to the velocity of the boundary current, they do place bounds on the magnitude of recirculation between the boundary current and the interior. On average, a volume equal to the boundary current transport is entrained/detrained over a length scale of about 5000 km. About half of the overflow water entering the western basin of North Atlantic since the mid-1960's has been mixed into the deep Labrador Sea and subpolar gyre. The effects of tracer surface boundary conditions on thermocline ventilation and oxygen utilization rate estimates are discussed. Tracers that equilibrate rapidly with the atmosphere, such as 3He and the CFCs, have lower apparent ventilation time scales than t racers, such as tritium and radiocarbon, t hat are reset slowly in the surface layer. The results of a simple box-mixing model are compared with tritium and 3He data from a 1979 survey of the eastern subtropical North Atlantic. On shallow density surfaces, the computed tritium ventilation rates are two to three times slower than those for 3He; deeper in the thermocline, the two tracer ventilation rates converge. This trend may be related to the decreasing effectiveness of 3He gas exchange in equilibrating the deeper winter mixed layers of the more northerly isopycnal outcrops. Box models using limited surface exchange tracers (e.g. tritium and 14C) can under predict oxygen utilization rates (OUR) by up to 3 times due to differences between tracer and oxygen boundary conditions while 3He may overestimate OUR by 10- 20%. I present and discuss the distributions of two chlorofluorocarbons (CFCs) in the eastern North Atlantic measured on a 1988 hydrographic cruise between Iceland and the equator. CFC tagged seawater fills the entire sub-polar water column and subtropical thermocline. Measurable CFC levels are found at the ocean bottom as far south as 35°N; the CFC penetration depth shoals to about 750 meters in the tropics. The CFC data are used to illustrate the ventilation time-scales for the water masses in the eas tern basin and to calculate OUR values in the subtropical thermocline. The CFC data in the tropical oxygen minimum off of Africa are significantly lower than the values on similar density surfaces in the subtropics, providing support for the idea that the tropical oxygen minima are controlled primarily by physical rather than biological mechanisms. The evolution of the tropical and subtropical CFC distributions between the 1972-73 TTO/TAS program and the 1988 cruise are also examined. Other features of the CFC data include a clear signal of Labrador Sea Water mid-depth ventilation, a CFC-enriched overflow water boundary current along the Iceland slope, a northward flowing deep boundary current along the eastern margin of the basin, and a mid-depth equatorial plume of upper North Atlantic Deep Water.

Description: The work carried out in my thesis has been supported in part by a National Science Foundation Graduate Fellowship and by grants OCE-8615289 and OCE-8800957 from the National Science Foundation.

Description: Author Posting. © American Meteorological Society, 2012. 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 42 (2012): 1684–1700, doi:10.1175/JPO-D-11-0230.1.

Description: The influences of precipitation on water mass transformation and the strength of the meridional overturning circulation in marginal seas are studied using theoretical and idealized numerical models. Nondimensional equations are developed for the temperature and salinity anomalies of deep convective water masses, making explicit their dependence on both geometric parameters such as basin area, sill depth, and latitude, as well as on the strength of atmospheric forcing. In addition to the properties of the convective water, the theory also predicts the magnitude of precipitation required to shut down deep convection and switch the circulation into the haline mode. High-resolution numerical model calculations compare well with the theory for the properties of the convective water mass, the strength of the meridional overturning circulation, and also the shutdown of deep convection. However, the numerical model also shows that, for precipitation levels that exceed this critical threshold, the circulation retains downwelling and northward heat transport, even in the absence of deep convection.

Description: This study was supported by the National Science Foundation underGrantsOCE-0850416, OCE-0959381, andOCE-0859381.

Description: 2013-04-01

Description: 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 March 1988

Description: Inverse methods are applied to historical hydrographic data to address two aspects of the general circulation of the Atlantic Ocean. The method allows conservation statements for mass and other properties, along with a variety of other constraints, to be combined in a dynamically consistent way to estimate the absolute velocity field and associated property transports. The method is first used to examine the exchange of mass and heat between the South Atlantic and the neighboring ocean basins. The Antarctic Circumpolar Current (ACC) carries a surplus of intermediate water into the South Atlantic through Drake Passage which is compensated by a surplus of deep and bottom water leaving the basin south of Africa. As a result, the ACC loses .25±.18x1015 W of heat in crossing the Atlantic. At 32°S the meridional flux of heat is .25±.19x1015 W equatorward, consistent in sign but smaller in magnitude than other recent estimates. This heat flux is carried primarily by a meridional overturning cell in which the export of 17 Sv of North Atlantic Deep Water (NADW) is balanced by an equatorward return flow equally split between the surface layers, and the intermediate and bottom water. No "leak" of warm Indian Ocean thermocline water is necessary to account for the equatorward heat flux across 32°S; in fact, a large transfer of warm water from the Indian Ocean to the Atlantic is found to be inconsistent with the present data set. Together these results demonstrate that the Atlantic as a whole acts to convert intermediate water to deep and bottom water, and thus that the global thermohaline cell associated with the formation and export of NADW is closed primarily by a "cold water path," in which deep water leaving the Atlantic ultimately returns as intermediate water entering the basin through Drake Passage. The second problem addressed concerns the circulation and property fluxes across 24°and 36°N in the subtropical North Atlantic. Conservation statements are considered for the nutrients as well as mass, and the nutrients are found to contribute significant information independent of temperature and salinity. Silicate is particularly effective in reducing the indeterminacy of circulation estimates based on mass conservation alone. In turn, the results demonstrate that accurate estimates of the chemical fluxes depend on relatively detailed knowledge of the circulation. The zonal-integral of the circulation consists of an overturning cell at both latitudes, with a net export of 19 Sv of NADW. This cell results in a poleward heat flux of 1.3±.2x1015 Wand an equatorward oxygen flux of 2900±180 kmol S-l across each latitude. The net flux of silicate is also equatorward: 138±38 kmol s-1 and 152±56 kmol s -1 across 36°and 24° N, respectively. However, in contrast to heat and oxygen, the overturning cell is not the only important mechanism responsible for the net silicate transport. A horizontal recirculation consisting of northward flow of silica-rich deep water in the eastern basin balanced by southward flow of low silica water in the western basin results in a significant silicate flux to the north. The net equatorward flux is thus smaller than indicated by the overturning cell alone. The net flux of nitrate across 36°N is n9±35 kmol 8- 1 to the north and is indistinguishable from zero at 24°N (-8±39 kmol 8-1 ), leading to a net divergence of nitrate between these two latitudes. Forcing the system to conserve nitrate leads to an unreasonable circulation. The dominant contribution to the nitrate flux at 36°N results from the correlation of strong northward flow and relatively high nitrate concentrations in the sub-surface waters of the Gulf Stream. The observed nitrate divergence between 24°and 36°N, and convergence north of 36°N, can be accounted for by a shallow cell in which the northward flow of inorganic nitrogen (nitrate) in the Gulf Stream is balanced by a southward flux of dissolved organic nitrogen in the recirculation gyre. Oxidation of the dissolved organic matter during its transit of the subtropical gyre supplies the required source of regenerated nitrate to the Gulf Stream and consumes oxygen, consistent with recent observations of oxygen utilization in the Sargasso Sea.

Description: This research was supported by NASA under contract NAG5-534 and NSF under contract OCE-8521685.

Description: Author Posting. © American Meteorological Society, 2011. 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 41 (2011): 889–910, doi:10.1175/2010JPO4496.1.

Description: This paper examines interaction between a barotropic point vortex and a steplike topography with a bay-shaped shelf. The interaction is governed by two mechanisms: propagation of topographic Rossby waves and advection by the forcing vortex. Topographic waves are supported by the potential vorticity (PV) jump across the topography and propagate along the step only in one direction, having higher PV on the right. Near one side boundary of the bay, which is in the wave propagation direction and has a narrow shelf, waves are blocked by the boundary, inducing strong out-of-bay transport in the form of detached crests. The wave–boundary interaction as well as out-of-bay transport is strengthened as the minimum shelf width is decreased. The two control mechanisms are related differently in anticyclone- and cyclone-induced interactions. In anticyclone-induced interactions, the PV front deformations are moved in opposite directions by the point vortex and topographic waves; a topographic cyclone forms out of the balance between the two opposing mechanisms and is advected by the forcing vortex into the deep ocean. In cyclone-induced interactions, the PV front deformations are moved in the same direction by the two mechanisms; a topographic cyclone forms out of the wave–boundary interaction but is confined to the coast. Therefore, anticyclonic vortices are more capable of driving water off the topography. The anticyclone-induced transport is enhanced for smaller vortex–step distance or smaller topography when the vortex advection is relatively strong compared to the wave propagation mechanism.

Description: Y. Zhang acknowledges the support of theMIT-WHOI Joint Programin Physical Oceanography, NSF OCE-9901654 and OCE-0451086. J. Pedlosky acknowledges the support of NSF OCE- 9901654 and OCE-0451086.

Description: 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 June 1998

Description: A freshwater plume often forms when a river or an estuary discharges water onto the continental shelf. Freshwater plumes are ubiquitous features of the coastal ocean and usually leave a striking signature in the coastal hydrography. The present study combines both hydrographic data and idealized numerical simulations to examine how ambient currents and winds influence the transport and mixing of plume waters. The first portion of the thesis considers the alongshore transport of freshwater using idealized numerical simulations. In the absence of any ambient current, the downstream coastal current only carries a fraction of the discharged fresh water; the remaining fraction recirculates in a continually growing "bulge" of fresh water in the vicinity of the river mouth. The fraction of fresh water transported in the coastal current is dependent on the source conditions at the river mouth. The presence of an ambient current augments the transport in the plume so that its freshwater transport matches the freshwater source. For any ambient current in the same direction as the geostrophic coastal current, the plume will evolve to a steady-state width. A key result is that an external forcing agent is required in order for the entire freshwater volume discharged by a river to be transported as a coastal current. The next section of the thesis addresses the wind-induced advection of a river plume, using hydrographic data collected in the western Gulf of Maine. The observations suggest that the plume's cross-shore structure varies markedly as a function of fluctuations in alongshore wind forcing. Consistent with Ekman dynamics, upwelling favorable winds spread the plume offshore, at times widening it to over 50 km in offshore extent, while downwelling favorable winds narrow the plume width to a few Rossby radii. Near-surface current meters show significant correlations between cross-shore currents and alongshore wind stress, consistent with Ekman theory. Estimates of the terms in the alongshore momentum equation calculated from moored current meter arrays also indicate an approximate Ekman balance within the plume. A significant correlation between alongshore currents and alongshore wind stress suggests that interfacial drag may be important. The final section of the thesis is an investigation of the advection and mixing of a surface-trapped river plume in the presence of an upwelling favorable wind stress, using a three-dimensional model in a simple, rectangular domain. Model simulations demonstrate that the plume thins and is advected offshore by the crossshore Ekman transport. The thinned plume is susceptible to significant mixing due to the vertically sheared horizontal currents. The first order plume response is explained by Ekman dynamics and a Richardson number mixing criterion. Under a sustained wind event, the plume evolves to a quasi-steady, uniform thickness. The rate of mixing slowly decreases for longer times as the stratification in the plume weakens, but mixing persists under a sustained upwelling wind until the plume is destroyed. Mixing is most intense at the seaward plume front due to an Ekman straining mechanism in which the advection of cross-shore salinity gradients balances vertical mixing. The mean mixing rate observed in the plume is consistent with the mixing power law suggested by previous studies of I-D mixing, in spite of the two-dimensional dynamics driving the mixing in the plume.

Description: This research was funded by a National Science Foundation graduate fellowship, and Gulf of Maine Regional Marine Research Program grants UM-S227 and UM-S276.

Description: 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 June 2011

Description: Eastern oceanic boundary currents are subject to hydrodynamic instability, generate small scale features that are visible in satellite images and may radiate westward into the interior, where they can be modified by the large-scale circulations. This thesis studies the stability of an eastern boundary current with and without the large-scale flow influence in an idealized framework represented by barotropic quasi-geostrophic dynamics. The linear stability analysis of a meridional current with a continuous velocity profile shows that meridional eastern and western boundary currents support a limited number of radiating modes with long meridional and zonal wavelengths and small growth rates. However, the linearly stable, long radiating modes of an eastern boundary current can become nonlinearly unstable by resonating with short trapped unstable modes. This phenomenon is clearly demonstrated in the weakly nonlinear simulations. Results suggest that linearly stable longwave modes deserve more attention when the radiating instability of a meridional boundary current is considered. A large-scale flow affects the short trapped unstable mode and long radiating mode through different mechanisms. The large-scale flow modifies the structure of the boundary current to stabilize or destabilize the unstable modes, leading to a meridionally localized maximum in the perturbation kinetic energy field. The shortwave mode is accelerated or decelerated by the meridional velocity adjustment of the large-scale flow to have an elongated or a squeezed meridional structure, which is confirmed both in a linear WKB analysis and in nonlinear simulations. The squeezed or elongated unstable mode detunes the nonlinear resonance with the longwave modes, which then become less energetic. These two modes show different meridional structures in kinetic energy field because of the different mechanisms. In spite of the model simplicity, these results can potentially explain the formation of the zonal jets observed in altimeter data, and indicate the influence of the large-scale wind-driven circulation on eastern boundary upwelling systems in the real ocean. Studies with more realistic configurations remain future challenges.

Description: 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 February 1999

Description: This thesis studies the problems of generation and maintenance of recirculations by Gulf Stream instabilities. Observations show that the horizontal structure of the jet and its recirculations suffer significant changes in time. Here, the role of internal dynamics of the jet is isolated as one of the possible sources of such variability, and the differences between barotropic and baroclinic instabilities are investigated. The problem of recirculation development is considered in a framework of a free spin down of the 2-layer and the 1-layer, zonally symmetric, quasi-geostrophic jets. Linear stability analysis shows that in strongly baroclinic basic flows, eddies are capable of driving recirculations in the lower layer through the residual meridional circulation. In strongly barotropic jets, the linearly most unstable wave simply diffuses the jet. Nonlinear stability analysis indicates that recirculations are robust features of the 2-layer model. The strength of recirculations is a function of the model’s parameters. It increases with a decrease in the value of the nondimensional /3 due to potential vorticity homogenization constrained by enstrophy conservation. The recirculation strength is a non-monotonic function of the baroclinic velocity parameter; it is the strongest for strongly baroclinic basic flows, weakest for flows with intermediate baroclinic structure and of medium strength for strongly barotropic basic flows. Such non-monotonic behavior is the result of two different processes responsible for the recirculation development: linear eddy-mean flow interactions for strongly baroclinic basic flows and strongly nonlinear eddy-eddy and eddy-mean flow interaction for strongly barotropic flows. In the case of the reduced-gravity model, recirculations develop only for infinite deformation raduis. Basic flows with finite deformation radius are only weakly supercritical and therefore produced negligible recirculations after equilibration. The problem of maintenance of the recirculations is coupled to the questions of existence of low frequency variability and of multiple dynamical regimes of a system consisting of a quasi-geostrophic jet and its recirculations. The problem is studied in a framework of a 2-layer or a reduced-gravity colliding jets model which has no windforcing. Instead, it is forced by inflows and outflows through the open boundaries. Oniy the western boundary of the domain is closed, and a free slip boundary condition is used there. The results of the numerical experiments show that when oniy the mechanism of barotropic instability is present, the model has two energy states for a wide range of interfacial friction coefficients. The high energy state is characterized by well-developed recirculations and displays strong variability associated with either large recirculating gyres and a weak eddy field or small recirculations and a strong eddy field. The iow energy state is characterized by large meridional excursions in the separation point and large amplitude, westward propagating meanders that produce strong rings after interacting with the western wall. For physically relevant bottom friction values, the presence of baroclinic in stability in the recirculation regions of the 2-layer model allows for a unique dynamical regime characterized by well-developed recirculations in both layers. The low-frequency variability associated with the regime is weak and is related to meridional shifts in the position of the jet, to wrapping of the recirculations around each other, and to pulsations in their zonal extent. For strong bottom friction, the 2-layer model has only the mechanism of barotropic instability which reduces it to a 1 1/2-layer configuration; the model displays two dynamical regimes and strong low frequency variability in the upper layer, while the lower layer is strongly frictional.

Description: Financial support for this research was provided by NSF grant number OCE 9617848.

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution January 1999

Description: Today, deep waters produced in the North Atlantic are exported through the western South Atlantic. Antarctic intermediate water AAJW also enters the Atlantic in this region. Circumpolar deep water (CDW) fills the depths below AAIW and above and below northern source waters. A depth transect of cores from 1567-3909 m water depth in the western South Atlantic are ideally located to monitor inter-ocean exchange of deep water, and variations in the relative strength of northern versus southern source water production. Last glacial maximum (LGM) Cd/Ca and δ13C data indicate a nutrient-depleted intermediate-depth water mass. In the mid-depth western South Atlantic, a simple conversion of LGM δ13C data suggests significantly less nutrient enrichment than LGM Cd/Ca ratios, but Cd/Ca and δ13C data can be reconciled when plotted in CdW/δ13C space. Paired LGM Cd/Ca and δ13C data from mid-depth cores suggest increasingly nutrient rich waters below 2000 m, but do not require an increase in Southern Ocean water contribution relative to today. Cd/Ca data suggest no glacial-interglacial change in the hydrography of the deepest waters ofthe region. To maintain relatively low Cd/Ca ratios low nutrients in the deepest western South Atlantic waters, and in CDW in general, during the LGM requires an increased supply ofnutrient-depleted glacial North Atlantic intermediate water (GNA1W) and/or nutrient-depleted glacial Subantarctic surface waters to CDW to balance reduced NADW contribution to CDW. LGM Cd/Ca and δ13C data suggest strong GNA1W influence in the western South Atlantic which in turn implies export of GNAIW from the Atlantic, and entrainment of GNA1W into the Antarctic Circumpolar current.

Description: Author Posting. © American Meteorological Society, 2011. 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 24 (2011): 4844–4858, doi:10.1175/2011JCLI4130.1.

Description: The factors that determine the heat transport and overturning circulation in marginal seas subject to wind forcing and heat loss to the atmosphere are explored using a combination of a high-resolution ocean circulation model and a simple conceptual model. The study is motivated by the exchange between the subpolar North Atlantic Ocean and the Nordic Seas, a region that is of central importance to the oceanic thermohaline circulation. It is shown that mesoscale eddies formed in the marginal sea play a major role in determining the mean meridional heat transport and meridional overturning circulation across the sill. The balance between the oceanic eddy heat flux and atmospheric cooling, as characterized by a nondimensional number, is shown to be the primary factor in determining the properties of the exchange. Results from a series of eddy-resolving primitive equation model calculations for the meridional heat transport, overturning circulation, density of convective waters, and density of exported waters compare well with predictions from the conceptual model over a wide range of parameter space. Scaling and model results indicate that wind effects are small and the mean exchange is primarily buoyancy forced. These results imply that one must accurately resolve or parameterize eddy fluxes in order to properly represent the mean exchange between the North Atlantic and the Nordic Seas, and thus between the Nordic Seas and the atmosphere, in climate models.

Description: This study was supported by the National Science Foundation under Grants OCE-0726339 and OCE-0850416.

Description: 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 February 1999

Description: This thesis addresses the parameterization of the heat and momentum transporting properties of eddy motions for use in three-dimensional, primitive equation, z-coordinate atmosphere and ocean models. Determining the transport characteristics of these eddies is fundamental to understanding their effect on the large-scale ocean circulation and global climate. The approach is to transform the primitive equations to yield the altered 'transformed Eulerian mean' (TEM) equations. The assumption is made that the eddy motions obey quasigeostrophic dynamics while the mean flow obeys the primitive equations. With this assumption, the TEM framework leads to the eddies appearing as one term, which acts as a body force in the momentum equations. This force manifests itself as a flux of potential vorticity (PV) - a quantity that incorporates both eddy momentum and heat transporting properties. Moreover, the dynamic velocities are those of the residual mean circulation, a much more relevant velocity for understanding heat and tracer transport. Closure for the eddy PV flux is achieved through a flux-gradient relationship, which directs the flux down the large scale PV gradient. For zonal flows, care is taken to ensure that the resulting force does not generate any net momentum, acting only to redistribute it. Neglect of relative vorticity fluxes in the PV flux yields the parameterization scheme of Gent and McWilliams. The approach is investigated by comparing a zonally-averaged parameterized model with a three dimensional eddy-resolving calculation of flow in a stress-driven channel. The stress at the upper surface is communicated down the water column to the bottom by eddy form drag. Moreover, lateral eddy momentum fluxes act to strengthen and sharpen the mean flow, transporting eastward momentum up its large scale gradient. Both the vertical momentum transfer and lateral, upgradient momentum transfer by eddies, are captured in the parameterized model. The advantages of this approach are demonstrated in two further zonal cases: 1) the spin-down of a baroclinic zone, and 2) the atmospheric jet stream. The time mean TEM approach and the eddy PV flux closure are explored in the context of an eddy-resolving closed basin flow which breaks the zonal symmetry. Decomposition of eddy PV fluxes into components associated with advective and dissipative effects suggest that the component associated with eddy flux divergence, and therefore forcing of the mean flow, is mainly directed down the large scale gradient and can be parameterized as before. Thus, the approach can be used to capture eddy transport properties for both zonal mean and time mean flows. The PV flux embodies both the eddy heat and momentum fluxes and so presents a more unified picture of their transferring properties. It therefore provides a powerful conceptual and practical framework for representing eddies in numerical models of the atmsophere and ocean.

Description: The work in this thesis was funded by grants from NSF, (OCE-9634331, OCE- 9503895), ONR (NOOOI4-95-1-0967), and by a fellowship from the Joint Program on the Science and Policy of Global Change at MIT.

Description: Author Posting. © American Meteorological Society, 2011. 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 41 (2011): 1874–1893, doi:10.1175/2011JPO4604.1.

Description: A two-dimensional cross-shelf model of the New England continental shelf and slope is used to investigate the mean cross-shelf and vertical circulation at the shelf break and their seasonal variation. The model temperature and salinity fields are nudged toward climatology. Annual and seasonal mean wind stresses are applied on the surface in separate equilibrium simulations. The along-shelf pressure gradient force associated with the along-shelf sea level tilt is tuned to match the modeled and observed depth-averaged along-shelf velocity. Steady-state model solutions show strong seasonal variation in along-shelf and cross-shelf velocity, with the strongest along-shelf jet and interior onshore flow in winter, consistent with observations. Along-shelf sea level tilt associated with the tuned along-shelf pressure gradient increases shoreward because of decreasing water depth. The along-shelf sea level tilt varies seasonally with the wind and is the strongest in winter and weakest in summer. A persistent upwelling is generated at the shelf break with a maximum strength of 2 m day−1 at 50-m depth in winter. The modeled shelfbreak upwelling differs from the traditional view in that most of the upwelled water is from the upper continental slope instead of from the shelf in the form of a detached bottom boundary layer.

Description: WGZ was supported by the Woods Hole Oceanographic Institution postdoctoral scholarship program. GGGandDJMwere supported byONRGrant N-00014- 06-1-0739.

Description: 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 June 1998

Description: The evolution of a coastal ocean undergoing uniform surface heat loss is examined. The dynamics of this ocean are initially modulated by the intense vertical mixing driven by surface cooling. The strong vertical mixing prevents the formation of geostrophic flows and inhibits the cross-shelf flux of heat. The vertical mixing is eventually suppressed by the advective transport of cold, dense water offshore. Once this happens, alongshore geostrophic flows form, and become baroclinically unstable. The surface heat flux is then balanced by a cross-shelf eddy heat flux. Scales are found for the cross-shelf density gradient which results from this balance. Solutions for linear internal waves are found for a wedge-shaped bathymetry with bottom friction. Bottom friction is capable of entirely dissipating the waves before they reach the coast, and waves traveling obliquely offshore are reflected back to the coast from a caustic. The internal wave climate near two moorings of the Coastal Ocean Dynamics Experiment observation program is analyzed. The high frequency internal wave energy levels were elevated above the Garrett and Munk spectrum, and the spectrum becomes less red as one moves to the shore. The wave field is dominated by vertical-mode one waves, and internal wave energy propagates shoreward.

Description: This work was funded by an Office of Naval Research fellowship and and Office of Naval Research AASERT fellowship, N00014-95:-1-0746.

Description: 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 January 1988

Description: This study focuses on the zonal weakening, eastern termination and seasonal variations of the Atlantic equatorial undercurrent (EUC). The main and most original contribution of the dissertation is a detailed analysis of the Atlantic EUC simulated by Philander and Pacanowski's (1986) general circulation model (GCM), which provides a novel description of the dynamical regimes governing various regions of a nonlinear stratified undercurrent. Only in a narrow deep western region of the simulation does one find an approximately inertial regime corresponding to zonal acceleration. Elsewhere frictional processes cannot be ignored. The bulk of the mid-basin model EUC terminates in the overlying westward surface flow while only a small fraction (the deeper more inertial layers) terminates at the eastern coast. In agreement with observations, a robust feature of the GCM not present in simpler models is the apparent migration of the EUC core from above the thermocline in the west to below it in the east. In the GCM, this happens because the eastward flow is eroded more efficiently by vertical friction above the base of the thermocline than by lateral friction at greater depths. This mechanism is a plausible one for the observed EUC. A scale analysis using a depth scale which decreases with distance eastwards predicts the model zonal transition between western inertial and eastern inertio-frictional regimes. Historical and recent observations and simple models of the equatorial and coastal eastern undercurrents are reviewed, and a new analysis of current measurements in the eastern equatorial Atlantic is presented. Although the measurements are inadequate for definitive conclusions, they suggest that Lukas' (1981) claim of a spring surge of the Pacific EUC to the eastern coast and a seasonal branching of the EUC into a coastal southeastward undercurrent may also hold for the Atlantic Ocean. To improve the agreement between observed and modelled strength of the eastern undercurrent, it is suggested that the eddy coefficient of horizontal mixing should be reduced in future GCM simulations.

Description: This work was supported by NSF grants OCE82-14771, OCE82-08744 and OCE85-14885.

Description: 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 June 1998

Description: The steady states of two models of the double-gyre wind-driven ocean circulation are studied. The link between the steady state solutions of the models and their time-mean and low-frequency variability is explored to test the hypothesis that both stable and unstable fixed points influence shape the model's attractor in phase space. The steady state solutions of a barotropic double-gyre ocean model in which the wind-stress curl input of vorticity is balanced primarily by bottom friction are studied. The bifurcations away from a unique and stable steady state are mapped as a function of two nondimensional parameters, (δI,δS), which can be thought of as measuring respectively the relative importance of the nonlinear advection and bottom damping of relative vorticity to the advection of planetary vorticity. A highly inertial branch characterized by a circulation with transports far in excess of those predicted by Sverdrup balance is present over a wide range of parameters including regions of parameter space where other solutions give more realistic flows. For the range of parameters investigated, in the limit of large Reynolds number, δI,δS → ∞, the inertial branch is stable and appears to be unique. This branch is anti-symmetric with respect to the mid-basin latitude like the prescribed wind-stress curl. For intermediate values of δI,δS, additional pairs of mirror image non-symmetric equilibria come into existence. These additional equilibria have currents which redistribute relative vorticity across the line of zero wind-stress curl. This internal redist~ibution of vorticity prevents the solution from developing the large transports that are necessary for the anti-symmetric solution to achieve a global vorticity balance. Beyond some critical Reynolds number, the nonsymmetric solutions are unstable to time-dependent perturbations. Time-averaged solutions in' this parameter regime have transports comparable in magnitude to those of the non-symmetric steady state branch. Beyond a turning point, where the non-symmetric steady state solutions cease to exist, all the computed time-dependent model trajectories converge to the anti-symmetric inertial runaway solution. The internal compensation mechanism which acts through explicitly simulated eddies is itself dependent explicit dissipation parameter. Using the reduced-gravity quasigeostrophic model an investigation of the link between the steady state solutions and the model's low-frequency variability is conducted. If the wind-stress curl is kept anti-symmetric, successive pairs of non-symmetric equilibria come into existence via symmetry-breaking pitchfork bifurcations as the model's biharmonic viscosity is reduced. Succesive pairs of mirror image equilibria have an additional half meander in the jet. The distinct energy levels of the steady state solutiOris can be understood in part by there different inter-gyre fluxes of vorticity. Those solutions with weak inter-gyre fluxes of vorticity have large and energetic recirculation cells which remove excess vorticity through bottom friction. Those solutions with strong inter-gyre fluxes of vorticity have much smaller and ·less energetic recirculation cells. A significant fraction of the variance (30%) of the interface height anomaly can be accounted by four coherent structures which point away from the time-mean state and towards four steady state solutions in phase space. After removing the variance which projects onto the four modes, the remaining variance is reduced predominantly at low-frequencies, showing that these modes are linked to the low-frequency variability of the model. Furthermore, the time-averaged flow fields within distinct energy ranges show distinct patterns which are in turn similar to the distinct steady state solutions.

Description: 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 August 1998

Description: Double-diffusive processes are studied and parameterized, and their impacts on the oceanic thermohaline circulation are investigated by single-hemisphere numerical models and scaling analysis. Scaling analysis on the thermohaline circulation has been done under three types of surface boundary conditions. (a) Under "relaxation" conditions, there is a two-thirds power law dependence of the meridional overturning rate (and the poleward heat transport) on the diapycnal diffusivity. For any given external forcing, there is only one equilibrium state for the thermohaline circulation. (b) Under "flux" boundary conditions, there is a half power law dependence of the meridional overturning rate on the diapycnal diffusivity. Only one mode is possible for given external forcing. (c) Under "mixed" boundary conditions, multiple equilibria become possible. For given thermal forcing, the existence of multiple equilibria depends on the relative contributions of diapycnal diffusivity and the hydrologic forcing. Numerical experiments are implemented to test the above scaling arguments. Consistent results have been obtained under the above three types of boundary conditions. These provide a basis for understanding how the thermohaline circulation depends on the diapycnal diffusivity, which we know is influenced by the double-diffusive processes of "salt fingering" and "diffusive layering" in some parts of the ocean. In order to examine this issue, the double-diffusive processes are parameterized by diapycnal eddy diffusivities for heat and salt that are different and depend on the local density ratio, Rp= αTz/βSz. A background diffusivity is applied to represent turbulent mixing in the stratified environment. The implementation of this double-diffusive - parameterization in numerical models has significant impacts on the thermohaline circulation. (a) Under "relaxation" boundary conditions, the meridional overturning rate and the poleward heat transport are reduced, and water mass properties are also changed. Similar results are obtained under "flux" boundary conditions. (b) Under "mixed" boundary conditions, the critical freshwater flux for the existence of the thermal mode becomes smaller with the double-diffusive parameterization. The extent to which the thermohaline circulation is affected by double-diffusive processes depends on the magnitude of the freshwater forcing.

Description: This thesis is supported by a grant from the Ocean Sciences Division of the National Science Foundation, OCE94-155S9.

Description: 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 June 1998

Description: Efforts to understand the Arctic system have recently focused on the role in local and global circulation of waters from the Arctic shelf seas. In this study, steady-state exchanges between the Arctic shelves and the central basins are estimated using an inverse box model. The model accounts for data uncertainty in the estimates, and quantifies the solution uncertainty. Other features include resolution of the two-basin Arctic hydrographic structure two-way shelf-basin exchange in the surface mixed layer, the capacity for shelfbreak upwelling, and recognition that most inflows enter the Arctic via the shelves. Aggregate estimates of all fluxes across the Arctic boundary, with their uncertainties, are generated from flux estimates published between 1975 and 1997. From the aggregate estimates, mass-, heat-, and salt-conserving boundary flux estimates are derived, which imply a net flux of water from the shelves to the basins of 1.2±0.4 Sv. Due primarily to boundary flux data uncertainty, constraints of mass, heat, and salt conservation alone cannot determine how much shelf-basin exchange occurs via dense overflows, and how much via the surface mixed layer. Adding δ180 constraints, however, greatly reduces the uncertainty. Dense water flux from the shelves to the basins is necessary for maintaining steady state, but shelfbreak upwelling is not required. Proper representation of external sources feeding the shelves, rather than the basins, is important to obtain the full range of plausible steady solutions. Implications of the results for the study of Arctic change are discussed.

Description: This work was supported by National Science Foundation grant OPP-9422292 as part of the Arctic System Science ARCSS program, administered by the Office of Polar Programs.

Description: Author Posting. © American Meteorological Society, 2011. 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 41 (2011): 1182–1208, doi:10.1175/2010JPO4564.1.

Description: The authors use data collected by a line of tall current meter moorings deployed across the axis of the Kuroshio Extension (KE) jet at the location of maximum time-mean eddy kinetic energy to characterize the mean jet structure, the eddy variability, and the nature of eddy–mean flow interactions observed during the Kuroshio Extension System Study (KESS). A picture of the 2-yr record mean jet structure is presented in both geographical and stream coordinates, revealing important contrasts in jet strength, width, vertical structure, and flanking recirculation structure. Eddy variability observed is discussed in the context of some of its various sources: jet meandering, rings, waves, and jet instability. Finally, various scenarios for eddy–mean flow interaction consistent with the observations are explored. It is shown that the observed cross-jet distributions of Reynolds stresses at the KESS location are consistent with wave radiation away from the jet, with the sense of the eddy feedback effect on the mean consistent with eddy driving of the observed recirculations. The authors consider these results in the context of a broader description of eddy–mean flow interactions in the larger KE region using KESS data in combination with in situ measurements from past programs in the region and satellite altimetry. This demonstrates important consistencies in the along-stream development of time-mean and eddy properties in the KE with features of an idealized model of a western boundary current (WBC) jet used to understand the nature and importance of eddy–mean flow interactions in WBC jet systems.

Description: This work was supported by National Science Foundation funding for the KESS program under Grants OCE-0220161 (SW, NGH, and SRJ), OCE- 0825550 (SW), OCE-0850744 (NGH), and OCE-0849808 (SRJ). SW was also supported by the MIT Presidential Fellowship. The financial assistance of the Houghton Fund, the MIT Student Assistance Fund, and WHOI Academic Programs is also gratefully acknowledged.

Description: 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 August 1999

Description: The thesis investigates the circulation at a 76-m deep study site on the southern flank of Georges Bank, a shallow submarine bank located between the deeper Gulf of Maine and the continental slope. Emphasis is placed on the vertical structure of the bottom boundary layer driven by the semi diurnal tides and the flow field's response to wind forcing. The observational analysis presented here is based on a combination of moored array and bottom tripod-mounted current, temperature, conductivity, and meteorological data taken between February and August 1995. Results from the bottom boundary layer analysis are compared to numerical model predictions for tidal flow over rough bottom topography. The flow response to wind forcing is examined and brought into context with the existing understanding of the wind-induced circulation in the Georges Bank region. Particular attention is given to the vertical distribution of the wind-driven currents, whose variation with background stratification is discussed and compared to observations from open ocean studies.

Description: The research presented in this thesis was generously supported by the National Science Foundation under grants OCE 93-13671 and OCE 96-32357 as part of the U.S. GLOBEC/Georges Bank Program.

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 February 1999

Description: A state-of-the-art, high-resolution ocean general circulation model is used to estimate the time-dependent global ocean heat transport and investigate its dynamics. The north-south heat transport is the prime manifestation of the ocean’s role in global climate, but understanding of its variability has been fragmentary owing to uncertainties in observational analyses, limitations in models, and the lack of a convincing mechanism. These issues are addressed in this thesis. Technical problems associated with the forcing and sampling of the model, and the impact of high-frequency motions are discussed. Numerical schemes are suggested to remove the inertial energy to prevent aliasing when the model fields are stored for later analysis. Globally, the cross-equatorial, seasonal heat transport fluctuations are close to +4.5 x 1015 watts, the same amplitude as the seasonal, cross-equatorial atmospheric energy transport. The variability is concentrated within 200 of the equator and dominated by the annual cycle. The majority of it is due to wind-induced current fluctuations in which the time-varying wind drives Ekman layer mass transports that are compensated by depth-independent return flows. The temperature difference between the mass transports gives rise to the time-dependent heat transport. The rectified eddy heat transport is calculated from the model. It is weak in the central gyres, and strong in the western boundary currents, the Antarctic Circumpolar Current, and the equatorial region. It is largely confined to the upper 1000 meters of the ocean. The rotational component of the eddy heat transport is strong in the oceanic jets, while the divergent component is strongest in the equatorial region and Antarctic Circumpolar Current. The method of estimating the eddy heat transport from an eddy diffusivity derived from mixing length arguments and altimetry data, and the climatological temperature field, is tested and shown not to reproduce the model’s directly evaluated eddy heat transport. Possible reasons for the discrepancy are explored.

Description: Funding for this research came from the Department of Defense under a National Defense Science and Engineering Graduate Fellowship. Financial support was also contributed by the National Science Foundation through grants #OCE-9617570 and #OCE-9730071, and the Tokyo Electric Power Company through the TEPCO/MIT Environmental Research Program. The author received partial support from an MIT Climate Modeling Fellowship, made possible by a gift from the American Automobile Manufacturers Association.

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

Description: A tomographic array is placed in a 2-layer, flat bottom, steady-wind driven quasi-geostrophic circulation model to investigate whether the analysis of acoustic travel time changes can detect large-scale barotropic oscillations. Time series of sea surface elevation and upper and lower layer meridional currents are generated for comparison against a series of acoustic travel times. The spectra of these time series exhibit a broad mesoscale peak near a period of 40 days. The spectrum of the acoustic travel time contains a significant peak due to a resonant barotropic oscillation with a period of 28.6 days which is not present in the spectra of the point measurements. In this numerical model, basin-scale tomographic measurements are a better method of sensing the large-scale resonant barotropic oscillations than are conventional point measurements because the tomographic system attenuates the "noise" from the mesoscale.

Description: Author Posting. © American Meteorological Society, 2011. 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 41 (2011): 241-246, doi:10.1175/2010JPO4557.1.

Description: The vertical dispersion of a tracer released on a density surface near 1500-m depth in the Antarctic Circumpolar Current west of Drake Passage indicates that the diapycnal diffusivity, averaged over 1 yr and over tens of thousands of square kilometers, is (1.3 ± 0.2) × 10−5 m2 s−1. Diapycnal diffusivity estimated from turbulent kinetic energy dissipation measurements about the area occupied by the tracer in austral summer 2010 was somewhat less, but still within a factor of 2, at (0.75 ± 0.07) × 10−5 m2 s−1. Turbulent diapycnal mixing of this intensity is characteristic of the midlatitude ocean interior, where the energy for mixing is believed to derive from internal wave breaking. Indeed, despite the frequent and intense atmospheric forcing experienced by the Southern Ocean, the amplitude of finescale velocity shear sampled about the tracer was similar to background amplitudes in the midlatitude ocean, with levels elevated to only 20%–50% above the Garrett–Munk reference spectrum. These results add to a long line of evidence that diapycnal mixing in the interior middepth ocean is weak and is likely too small to dictate the middepth meridional overturning circulation of the ocean.

Description: This material is based upon work supported by the National Science Foundation Grants OCE-0622825,OCE-0622670, OCE-0622630, and OCE-0623177.

Description: 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 February 2000

Description: Benthic foraminiferal δ13C, Cd/Ca, and Ba/Ca are important tools for reconstructing nutrient distributions, and thus ocean circulation, on glacial-interglacial timescales. However, each tracer has its own "artifacts" that can complicate paleoceanographic interpretations. It is therefore advantageous to measure multiple nutrient proxies with the aim of separating the various complicating effects. Zn/Ca is introduced as an important aid toward this goal. Benthic (Hoeglundina elegans) Cd/Ca ratios from the Bahama Banks indicate that the North Atlantic subtropical gyre was greatly depleted in nutrients during the last glacial maximum (LGM). A high-resolution Cd/Ca record from 965 m water depth suggests that Glacial North Atlantic Intermediate Water formation was strong during the LGM, weakened during the deglaciation, and strengthened again during the Younger Dryas cold period. Comparison of Cd/Ca and δ13C data reveals apparent short-term changes in carbon isotopic air-sea signatures. Benthic foraminiferal Zn/Ca could be a sensitive paleoceanographic tracer because deep water masses have characteristic Zn concentrations that increase about ten-fold from the deep North Atlantic to the deep North Pacific. A "core top calibration" shows that Zn/Ca is controlled by bottom water dissolved Zn concentration and, like Cd/Ca and BalCa, by bottom water saturation state with respect to calcite Since Zn/Ca responds to a different range of saturation states than Cd/Ca, the two may be used together to evaluate changes in deep water carbonate ion (CO32-) concentration. Zn/Ca and Cd/Ca ratios in the benthic foraminifer Cibicidoides wuellerstorfi exhibit large fluctuations over the past 100,000 years in a deep (3851 m) eastern equatorial Pacific sediment core. The data imply that bottom water CO32- concentrations were lowest during glacial Marine Isotope Stage 4 and highest during the last deglaciation. LGM CO32- concentrations appear to have been within a few μmol kg-1 of modern values. Deep North Atlantic Cd/Ca ratios imply much higher nutrient concentrations during the LGM. Although such data have usually been explained by a northward penetration of Southern Ocean Water (SOW), it has been suggested that they could result from increased preformed nutrient levels in the high-latitude North Atlantic or by increased aging of lower North Atlantic Deep Water (NADW). Glacial Zn/Ca data, however, require a substantially increased mixing with SOW and thus a reduction in NADW formation. Large changes in carbon isotopic air-sea exchange are invoked to reconcile benthic δ13C and trace metal data.

Description: This work was supported by a JOIlUSSAC Ocean Drilling Fellowship (subgrant JSG-CY 12-4), the R. H. Cole Ocean Ventures Fund, the Joint Program Education Office, and the National Science Foundation (grants OCE-9402804 and OCE-9503135 to W. Curry, and grant OCE-9633499 to D. Oppo).

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution April 2000

Description: A simple quasi-two dimensional dynamical model of Thermohaline circulation (THC) is developed, assuming that the mixing only occurs near western and eastern boundary layers. When the surface density is prescribed, the climatically important quantities, such as the strength of overturning and meridional heat transport, are related to the zonal integral over the vigorously mixing regions and scaled as (KvΔx )2/3. The numerical results suggest that the density difference between eastern and western boundaries play an important role in the meridional overturning. The eastern boundary is characterized by the upwelling on top of downwelling. The western boundary layer is featured by the universal upwelling. The inefficiency of diffusion heat transport accounts for the narrowness of sinking region and shallowness of overturning cell in one-hemisphere. The experiments with other surface boundary conditions are also explored. The circulation patterns obtained are similar under various surface temperature distributions, suggesting these are very robust features of THC. The role of boundary mixing is further explored in global ocean. The 2 1/2 dynamical model is extended to two-hemisphere ocean. Additional dynamics such as Rayleigh friction and abyssal water properties are taken into account. A set of complicated governing equations are derived and numerically solved to obtain steady state solution. The basic circulation features are revealed in our dynamical model. An equtorially asymmetric meridional circulation is observed due to small perturbation at the surface temperature in the high latitude. The density differences between eastern and western boundaries are distinct in both hemispheres. This is achieved during the spin-up process. Although the dynamical model results agree well with OGCM results in one-hemisphere, several important dynamics are lacking and exposed in two-hemisphere experiments. We need to consider horizontal advection terms which will effectively advect positive density anomalies across the equator and form the deep water for the entire system.

Description: 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 1999

Description: Climate simulation with numerical oceanic models requires a proper parameterization scheme in order to represent the effects of unresolved mesoscale eddies. Even though a munber of schemes have been proposed and some have led to improvements in the simulation of the bulk climatological properties, the success of the parameterizations in representing the mesoscale eddies has not been investigated in detail. This thesis examines the role of eddies in a 105-years long basin scale eddy resolving simulation with the MIT General Circulation Model (GCM) forced by idealized wind stress and relaxation to prescribed meridional temperature; this thesis also evaluates the Fickian diffusive, the diabatic Green-Stone (GS) and the quasi-adiabatic Gent-McWilliams (GM) parameterizations in a diagnostic study and a series of coarse resolution experiments with the same model in the same configuration. The mesoscale eddies in the reference experiment provide a significant contribution to the thermal balance in limited areas of the domain associated with the upper 1000M of the boundary regions. Specifically designed diagnostic tests of the schemes show that the horizontal and vertical components of the parameterized flux are not simultaneously downgradient to the eddy heat flux. The transfer vectors are more closely aligned with the isopycnal surfaces for deeper layers, thus demonstrating the adiabatic nature of the eddy heat flux for deeper layers. The magnitude of the coefficients is estimated to be consistent with traditionally used values. However, the transfer of heat associated with timedependent motions is identified as a complicated process that cannot be fully explained with any of the local parameterization schemes considered. The eddy parameterization schemes are implemented in the coarse resolution configuration with the same model. A series of experiments exploring the schemes' parameter space demonstrate that Fickian diffusion has the least skill in the climatological simulations because it overestimates the temperature of the deep ocean and underestimates the total heat transport. The GS and GM schemes perform better in the simulation of the bulk climatological properties of the reference solution, although the GM scheme in particular produces an ocean that is consistently colder than the reference state. Comparison of the eddy heat flux divergence with the parameterized divergences for typical parameter values demonstrates that the success of the schemes in the climatological simulation is not related to the representation of the eddy heat flux but to the representation of the overall internal mixing processes.

Description: The financial support for this research was provided by ONR grant number NOOOl4- 98-1-0881, Alliance for Global Sustainability and American Automobile Manufactures Association.

Description: 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 June 1999

Description: Data assimilation methods, such as the Kalman filter, are routinely used in oceanography. The statistics of the model and measurement errors need to be specified a priori. In this study we address the problem of estimating model and measurement error statistics from observations. We start by testing the Myers and Tapley (1976, MT) method of adaptive error estimation with low-dimensional models. We then apply the MT method in the North Pacific (5°-60° N, 132°-252° E) to TOPEX/POSEIDON sea level anomaly data, acoustic tomography data from the ATOC project, and the MIT General Circulation Model (GCM). A reduced state linear model that describes large scale internal (baroclinic) error dynamics is used. The MT method, closely related to the maximum likelihood methods of Belanger (1974) and Dee (1995), is shown to be sensitive to the initial guess for the error statistics and the type of observations. It does not provide information about the uncertainty of the estimates nor does it provide information about which structures of the error statistics can be estimated and which cannot. A new off-line approach is developed, the covariance matching approach (CMA), where covariance matrices of model-data residuals are "matched" to their theoretical expectations using familiar least squares methods. This method uses observations directly instead of the innovations sequence and is shown to be related to the MT method and the method of Fu et al. (1993). The CMA is both a powerful diagnostic tool for addressing theoretical questions and an efficient estimator for real data assimilation studies. It can be extended to estimate other statistics of the errors, trends, annual cycles, etc. Twin experiments using the same linearized MIT GCM suggest that altimetric data are ill-suited to the estimation of internal GCM errors, but that such estimates can in theory be obtained using acoustic data. After removal of trends and annual cycles, the low frequency /wavenumber (periods〉 2 months, wavelengths〉 16°) TOPEX/POSEIDON sea level anomaly is of the order 6 cm2. The GCM explains about 40% of that variance. By covariance matching, it is estimated that 60% of the GCM-TOPEX/POSEIDON residual variance is consistent with the reduced state linear model. The CMA is then applied to TOPEX/POSEIDON sea level anomaly data and a linearization of a global GFDL GCM. The linearization, done in Fukumori et al.(1999), uses two vertical mode, the barotropic and the first baroclinic modes. We show that the CMA method can be used with a global model and a global data set, and that the estimates of the error statistics are robust. We show that the fraction of the GCMTOPEX/ POSEIDON residual variance explained by the model error is larger than that derived in Fukumori et al.(1999) with the method of Fu et al.(1993). Most of the model error is explained by the barotropic mode. However, we find that impact of the change in the error statistics on the data assimilation estimates is very small. This is explained by the large representation error, i.e. the dominance of the mesoscale eddies in the TIP signal, which are not part of the 20 by 10 GCM. Therefore, the impact of the observations on the assimilation is very small even after the adjustment of the error statistics. This work demonstrates that simultaneous estimation of the model and measurement error statistics for data assimilation with global ocean data sets and linearized GCMs is possible. However, the error covariance estimation problem is in general highly underdetermined, much more so than the state estimation problem. In other words there exist a very large number of statistical models that can be made consistent with the available data. Therefore, methods for obtaining quantitative error estimates, powerful though they may be, cannot replace physical insight. Used in the right context, as a tool for guiding the choice of a small number of model error parameters, covariance matching can be a useful addition to the repertory of tools available to oceanographers.

Description: This research was supported by NASA through Earth Science Fellowship under contract number NGTS-30084, Global Change Science Fellowship under contract number NGT-30309, contract NAG 5-3274, and NASA-JPL under contract number 958125.

Description: 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 2009

Description: Interaction between the Antarctic Circumpolar Current and the continental slope/shelf in the Marguerite Bay and west Antarctic Peninsula is examined as interaction between a wind-driven channel flow and a zonally uniform slope with a bay-shaped shelf to the south. Two control mechanisms, eddy advection and propagation of topographic waves, are identified in barotropic vortex-escarpment interactions. The two mechanisms advect the potential vorticity (PV) perturbations in opposite directions in anticyclone-induced interactions but in the same direction in cyclone-induced interactions, resulting in dramatic differences in the two kinds of interactions. The topographic waves become more nonlinear near the western(eastern if in the Northern Hemisphere) boundary of the bay, where strong cross-escarpment motion occurs. In the interaction between a surface anticyclone and a slope penetrating into the upper layer in a two-layer isopycnal model, the eddy advection decays on length scales on the order of the internal deformation radius, so shoreward over a slope that is wider than the deformation radius, the wave mechanism becomes noticeably significant. It acts to spread the cross-isobath transport in a much wider range while the transport directly driven by the anticyclone is concentrated in space. A two-layer wind-driven channel flow is constructed to the north of the slope in the Southern Hemisphere, spontaneously generating eddies through baroclinic instability. A PV front forms in the first layer shoreward of the base of the topography due to the lower-layer eddy-slope interactions. Perturbed by the jet in the center of the channel, the front interacts with the slope/shelf persistently yet episodically, driving a clockwise mean circulation within the bay as well as crossisobath transport. Both the transports across the slope edge and out of the bay are comparable with the maximum Ekman transport in the channel, indicative of the significance of the examined mechanism. The wave-boundary interaction identified in the barotropic model is found essential for the out-of-bay transport and responsible for the heterogeneity of the transport within the bay. Much more water is transported out of the bay from the west than from the east, and the southeastern area is the most isolated region. These results suggest that strong out-of-bay transport may be found near the western boundary of the Marguerite Bay while the southeastern region is a retention area where high population of Antarctic krill may be found.

Description: 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 August 1984

Description: This thesis addresses several aspects of the problem of determining the effect of the low-frequency eddy variability on the mean circulation of the Western North Atlantic. A framework for this study is first established by scale analysis of the eddy and mean terms in the mean momentum, vorticity, and heat balances in three regions of the Western North Atlantic -- the northern recirculation, the southern recirculation, and the mid-ocean. The data from the last decade of field experiments suggest somewhat different conclusions from the earlier analysis of Harrison (1980). In the momentum balance we confirm that the eddy terms are negligible compared to the lowest order mean geostrophic balance. The eddy term may be an 0(1) term in the vorticity balance only in the northern recirculation region where the mean flow is anisotropic. In the mean heat balance, if the mean temperature advection is scaled using the thermal wind relation, then the eddy heat flux is negligible in the mid-ocean, but it may be important in the recirculation areas. For all the balances the eddy terms are comparable to or an order of magnitude larger than the mean advective terms. We conclude from the scale analysis that the eddy field is most likely to be important in the Gulf Stream recirculation region. These balances are subsequently examined in more detail using data from the Local Dynamics Experiment (LDE). Several inconsistencies are first shown in McWilliams' (1983) model for the mean dynamical balances in the LDE. The sampling uncertainties do not allow us to draw conclusions about the long-term dynamical balances. However, it is shown that if we assume that the linear vorticity balance holds between the surface and the thermocline for a finite record, then the vertical velocity induced by the eddy heat flux divergence is non-zero. The local effect of the mesoscale eddy field on the mean potential vorticity distribution of the Gulf Stream recirculation region is determined from the quasigeostrophic eddy potential vorticity flux. This flux is calculated by finite difference of current and temperature time series from the Local Dynamics Experiment. This long-term array of moorings is the only experimental data from which the complete eddy flux can be calculated. The total eddy flux is dominated by the term due to the time variation in the thickness of isopycnal layers. This thickness flux is an order of magnitude larger than the relative vorticity flux. The total flux is statistically significant and directed 217° T to the southwest with a magnitude of 1.57 x 10 -5 cm/2s. The direction of the eddy flux with respect to the mean large scale potential vorticity gradient from hydrographic data indicates that eddies in this region tend to reduce the mean potential vorticity gradient. The results are qualitatively consistent with numerical model results and with other data from the Gulf Stream recirculation region. We find that the strength of the eddy transfer in the enstrophy cascade is comparable to the source terms in the mean enstrophy balance. The Austauch coefficient for potential vorticity mixing is estimated to be 0(107cm2/sec). An order of magnitude estimate of the enstrophy dissipation due only to the internal wave field shows that other processes must be important in enstrophy dissipation. The measured eddy potential vorticity fluxes are compared to the linear stability model of Gill, Green, and Simmons (1974). An earlier study (Hogg, 1984) has shown agreement between the empirical orthogonal modes of the data and the predicted wavenumbers, growth rates, and phase speeds of the most unstable waves. However, we show substantial disagreement in a comparison of the higher moments the eddy heat and potential vorticity fluxes. Because the critical layer of the model is located near the surface, the model predicts that most of the eddy potential vorticity and eddy heat flux should occur within about 300 meters of the surface. The data show much greater deep eddy heat flux than predicted by the model. It is suggested that the unstable modes in the ocean have a longer vertical scale because of the reduction in the buoyancy frequency near the surface. The evidence for in situ instability is also examined in the decay region of the Gulf Stream from an array of current and temperature recorders. Although there is vertical phase propagation in the empirical orthogonal modes for some of the variables at some of the moorings, there is not much evidence for a strong ongoing process of wave generation.

Description: This research has been conducted under NSF contract numbers OCE 77-19403, ATM 79-21431, and OCE 82-00154.

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 1987

Description: Several problems connected by the theme of thermal forcing are addressed herein. The main topic is the stratification and flow field resulting from imposing a specified heat flux on a fluid that is otherwise confined to a rigid insulating basin. In addition to the traditional eddy viscosity and diffusivity, turbulent processes are also included by a convective overturning adjustment at locations where the local density field is unstable. Two classes of problems are treated. The first is the large scale meridional pattern of a fluid in an annulus. The detailed treatment is carried out in two steps. In the beginning (chapter 2) it is assumed that the fluid is very diffusive, hence, to first approximation no flow field is present. It is found that the convective overturning adjustment changes the character of the stratification in all the regions that are cooled from the top, resulting in a temperature field that is nearly depth independent in the northernmost latitudes. The response to a seasonal cycle in the forcing, and the differences between averaging the results from the end of each season compared to driving the fluid by a mean forcing are analyzed. In particular, the resulting sea surface temperature is warmer in the former procedure. This observation is important in models where the heat flux is sensitive to the gradient of air to sea surface temperatures. The analysis of the problem continues in chapter 5 where the contribution of the flow field is included in the same configuration. The dimensionless parameter controlling the circulation is now the Rayleigh number, which is a measure of the relative importance of gravitational and viscous forces. The effects of the convective overturning adjustment is investigated at different Rayleigh numbers. It is shown that not only is the stratification now always stable, but also that the vigorous vertical mixing reduces the effective Rayleigh number; thereby the flow field is more moderate, the thermocline deepens, and the horizontal surface temperature gradients are weaker. The interior of the fluid is colder compared to cases without convective overturning, and, because the amount of heat in the system is assumed to be fixed, the surface temperature is warmer. The fluid is not only forced by a mean heat flux, or a seasonally varying one, but its behavior under permanent winter and summer conditions is also investigated. A steady state for the experiments where the net heat flux does not vanish is defined as that state where the flow field and temperature structure are not changing with time except for an almost uniform temperature decrease or increase everywhere. It is found that when winter conditions prevail the circulation is very strong, while it is rather weak for continuous summer forcing. In contrast to those results, if a yearly cycle is imposed, the circulation tends to reach a minimum in the winter time and a maximum in the summer. This suggests that, depending on the Rayleigh number, there is a phase leg of several months between the response of the ocean and the imposed forcing. Differences between the two averaging procedures mentioned before are also observed when the flow field is present, especially for large Rayleigh numbers. The circulation is found to be weaker and the sea surface temperature colder in the mean of the seasonal realizations compared to the steady state derived by the mean forcing. As an extension to the numerical results, an analytic model is presented in chapter 4 for a similar annular configuration. The assumed dynamics is a bit different, with a mixed layer on top of a potential vorticity conserving interior. It is demonstrated that the addition of the thermal wind balance to the conservation of potential vorticity in the axially symmetric problem leads to the result that typical fluid trajectories in the interior are straight lines pointing downward going north to south. The passage of information in the system is surprisingly in the opposite sense to the clockwise direction of the flow. A model for water mass formation by buoyancy loss in the absence of a flow field is introduced in chapter 3. The idea behind it is to use the turbulent mixing parameterization to generate chimney-like structures in open water, followed by along-isopycnal advection and diffusion. This model can be applied to many observations of mode water. In particular, in this work it is related to the chimneys observed by the MEDOC Group (1970), and the Levantine Intermediate Water in the Eastern Mediterranean Basin. An analytic prediction of the depth of the water mass is derived and depends on the forcing and initial stratification. It suggests that the depth of shallow mode water like the 18°C water or the Levantine Intermediate Water would not be very sensitive to reasonable changes in atmospheric forcing. Similar conclusions were also reached by Warren (1972) by assuming that the temperature in the thermocline decreases linearly with depth, and by approximating the energy balance in a water column by a Newtonian cooling law.

Description: Author Posting. © American Meteorological Society, 2010. 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 40 (2010): 789-801, doi:10.1175/2009JPO4039.1.

Description: The issue of internal wave–mesoscale eddy interactions is revisited. Previous observational work identified the mesoscale eddy field as a possible source of internal wave energy. Characterization of the coupling as a viscous process provides a smaller horizontal transfer coefficient than previously obtained, with vh 50 m2 s−1 in contrast to νh 200–400 m2 s−1, and a vertical transfer coefficient bounded away from zero, with νυ + (f2/N2)Kh 2.5 ± 0.3 × 10−3 m2 s−1 in contrast to νυ + (f2/N2)Kh = 0 ± 2 × 10−2 m2 s−1. Current meter data from the Local Dynamics Experiment of the PolyMode field program indicate mesoscale eddy–internal wave coupling through horizontal interactions (i) is a significant sink of eddy energy and (ii) plays an O(1) role in the energy budget of the internal wave field.

Description: Author Posting. © American Meteorological Society, 2008. 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 38 (2008): 1486–1500, doi:10.1175/2007JPO3767.1.

Description: Fits of an annual harmonic to depth-average along-shelf current time series longer than 200 days from 27 sites over the Middle Atlantic Bight (MAB) continental shelf have amplitudes of a few centimeters per second. These seasonal variations are forced by seasonal variations in the wind stress and the cross-shelf density gradient. The component of wind stress that drives the along-shelf flow over most of the MAB mid- and outer shelf is oriented northeast–southwest, perpendicular to the major axis of the seasonal variation in the wind stress. Consequently, there is not a significant seasonal variation in the wind-driven along-shelf flow, except over the southern MAB shelf and the inner shelf of New England where the wind stress components forcing the along-shelf flow are north–south and east–west, respectively. The seasonal variation in the residual along-shelf flow, after removing the wind-driven component, has an amplitude of a few centimeters per second with maximum southwestward flow in spring onshore of the 60-m isobath and autumn offshore of the 60-m isobath. The spring maximum onshore of the 60-m isobath is consistent with the maximum river discharges in spring enhancing cross-shelf salinity gradients. The autumn maximum offshore of the 60-m isobath and a steady phase increase with water depth offshore of Cape Cod are both consistent with the seasonal variation in the cross-shelf temperature gradient associated with the development and destruction of a near-bottom pool of cold water over the mid and outer shelf (“cold pool”) due to seasonal variations in surface heat flux and wind stress.

Description: This research was funded by the Ocean Sciences Division of the National Science Foundation under Grants OCE-820773, OCE-841292, and OCE- 848961.

Description: Author Posting. © American Meteorological Society, 2008. 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 38 (2008): 1091-1106, doi:10.1175/2007JPO3805.1.

Description: A model of deep ocean circulation driven by turbulent mixing is produced in a long, rectangular laboratory tank. The salinity difference is substituted for the thermal difference between tropical and polar regions. Freshwater gently flows in at the top of one end, dense water enters at the same rate at the top of the other end, and an overflow in the middle removes the same amount of surface water as is pumped in. Mixing is provided by a rod extending from top to bottom of the tank and traveling back and forth at constant speed with Reynolds numbers 〉500. A stratified upper layer (“thermocline”) deepens from the mixing and spreads across the entire tank. Simultaneously, a turbulent plume (“deep ocean overflow”) from a dense-water source descends through the layer and supplies bottom water, which spreads over the entire tank floor and rises into the upper layer to arrest the upper-layer deepening. Data are taken over a wide range of parameters and compared to scaling theory, energetic considerations, and simple models of turbulently mixed fluid. There is approximate agreement with a simple theory for Reynolds number 〉1000 in experiments with a tank depth less than the thermocline depth. A simple argument shows that mixing and plume potential energy flux rates are equal in magnitude, and it is suggested that the same is approximately true for the ocean.

Description: The research was supported by the Ocean Climate Change Institute of Woods Hole Oceanographic Institution.

Description: Author Posting. © American Meteorological Society, 2008. 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 38 (2008): 2556-2574, doi:10.1175/2008JPO3666.1.

Description: Vertical profiles of horizontal velocity obtained during the Mid-Ocean Dynamics Experiment (MODE) provided the first published estimates of the high vertical wavenumber structure of horizontal velocity. The data were interpreted as being representative of the background internal wave field, and thus, despite some evidence of excess downward energy propagation associated with coherent near-inertial features that was interpreted in terms of atmospheric generation, these data provided the basis for a revision to the Garrett and Munk spectral model. These data are reinterpreted through the lens of 30 years of research. Rather than representing the background wave field, atmospheric generation, or even near-inertial wave trapping, the coherent high wavenumber features are characteristic of internal wave capture in a mesoscale strain field. Wave capture represents a generalization of critical layer events for flows lacking the spatial symmetry inherent in a parallel shear flow or isolated vortex.

Description: Salary support for this analysis was provided by Woods Hole Oceanographic Institution bridge support funds.

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 December 1999

Description: A new, global inversion is used to estimate the large scale oceanic circulation based on the World Ocean Circulation Experiment and Java Australia Dynamic Experiment hydrographic data. A linear inverse "box" model is used to combine consistently the transoceanic sections. The circulation is geostrophic with an Ekman layer at the surface and oceanic layers defined by neutral surfaces. Near-conservation of mass, salt and top-to-bottom silica is required and, in addition, heat and the phosphate-oxygen combination (170[P04]+[02]) are conserved in layers that are not in contact with the surface. A globally-consistent solution is obtained for a depth-independent adjustment to the thermal wind field, freshwater flux divergenees, the Ekman transport, and the advective and diffusive dianeutral fluxes between layers. A detailed error budget permits calculation of statistical uncertainties, taking into account both the non-resolved part of the solution and the systematic errors due to the temporal oceanic variability. The estimated water mass transports during the WOCE period (1985-1996) are generally similar to previous published estimates. However, important differences are found. In particular, the inflow of bottom waters into the Pacific Ocean is smaller than in most previous estimates. Utilization of property anomaly conservation constraints allows the estimation of significant dianeutral diffusivities in deep layers, with a global average of 3 ± lcm2s- 1 north of 30°S. Dianeutral transfers indicate that about 20 Sv of bottom water is formed in the Southern Ocean. Significant ocean-atmosphere heat fluxes are found, with a global heating of 2.3 ± 0.4PW in the tropical band and a corresponding cooling at high latitudes. The signature of a large-scale average export production is found for nutrients in several temperate regions. Despite the large uncertainties, the production magnitudes are consistent with independent measurements from sediment traps and isotopic data. Net nutrient sources or sinks are found in several regions, suggesting either transport of dissolved organic matter or a seasonal alias. Oxygen indicates large exchanges with the atmosphere, with intake at high latitudes and outgassing/remineralization at low latitudes.

Description: This work was supported in part by the Jet Propulsion Laboratory/CALTECH (contract #958125), and by gifts from Ford, General Motors, and Daimler-Chrysler to MIT's Climate Modelling Initiative.

Description: 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 February 1981

Description: The equatorial Pacific heat flow low, a major oceanic geothermal anomaly centered on the equatorial sediment bulge, was investigated using deeply penetrating heat flow probes (6-11 meters penetration) within three detailed surveys (400 km2) and along over 10,000 km of continuous seismic profiles (CSP). Previous heat flow measurements in this region defined a broad region characterized by a heat flux well below 1 HFU. We report 98 new measurements collected during cruises PLEIADES 3 and KNORR 73-4 that verify the anomalous nature of the heat flux and also define non-linear temperature gradients (concave down). Temperature field disturbances due to perturbations of a purely conductive heat transport regime are incapable of suitably explaining either of these observations . A simple model incorporating heat transport by both conduction and fluid convection through the sediments fits the observations. A volume flux of (hydrothermal) fluid in the range of 10-6 to 10-5 cm3/sec/cm2 (0.1 liter/yr/cm2) is required. The sense of the flow for all measurements exhibiting non-linear gradients is upward out of the sediment column; no evidence for the recharging of the system was observed. Investigation of a well-defined boundary of the low zone at 4°N and 114°W showed a transition from low and variable heat flow to values compatible with thermal models that correlated with a change in the nature of the basement from rough to smooth. A few outcrops occur in the area of rough basement, but otherwise the region is well-sedimented (greater than 200 meters). Measurements within a detailed survey centered at this transition showed a dramatic increase in heat flow from 1.21 HFU to values greater than 3 HFU over a horizontal distance of 10km. A similar transition from non-linear to linear temperature gradients was not observed as nearly every measurement was non-linear. Heat flow measurements located in well-sedimented, outcrop-free areas (A environments) were associated with linear gradients and a heat flux greater than 1 HFU, however, several of these values were well below the theoretical heat flow for the appropriate age crust. Values measured in environments other than A exhibited variable heat flow and non-linear gradients. The average value of measurements located in A environments within the equatorial Pacific heat flow low was 1.37±0.27 HFU. The previously reported average was 0.92±0.48 HFU based on several measurements from L-DGO cruise VEMA 24-3. The average heat flow measured at a survey located outside the low heat flow zone on crust of 55 ±5 m.a. was 1.76 ±0.30 HFU which is in good agreement with the theoretical value of 1.60. The measurements in this survey were not located in A environments suggesting that crustal convection has ceased or is greatly attenuated within crust of this age. Error analysis of the geothermal data reduction using the convective/conductive heat transport model suggests that the volume flux parameter is sensitive to temperature measurement errors greater than a few millidegrees. Volume fluxes less than 10-7 cm/sec are difficult to distinguish from the purely conductive case assuming instrumental accuracies of 0.001°C. Resolution of the volume flux deteriorates as heat flow decreases and is poor for values less than 0.5 HFU. A detailed survey located within the low zone confirmed previous measurements of low heat flow, however, due to the low value of heat flow (about 0.5 HFU) the small-scale variability could not be clearly defined.

Description: Author Posting. © American Meteorological Society, 2008. 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 38 (2008): 1203–1221, doi:10.1175/2007JPO3768.1.

Description: Analyses of current time series longer than 200 days from 33 sites over the Middle Atlantic Bight continental shelf reveal a consistent mean circulation pattern. The mean depth-averaged flow is equatorward, alongshelf, and increases with increasing water depth from 3 cm s−1 at the 15-m isobath to 10 cm s−1 at the 100-m isobath. The mean cross-shelf circulation exhibits a consistent cross-shelf and vertical structure. The near-surface flow is typically offshore (positive, range −3 to 6 cm s−1). The interior flow is onshore and remarkably constant (−0.2 to −1.4 cm s−1). The near-bottom flow increases linearly with increasing water depth from −1 cm s−1 (onshore) in shallow water to 4 cm s−1 (offshore) at the 250-m isobath over the slope, with the direction reversal near the 50-m isobath. A steady, two-dimensional model (no along-isobath variations in the flow) reproduces the main features of the observed circulation pattern. The depth-averaged alongshelf flow is primarily driven by an alongshelf pressure gradient (sea surface slope of 3.7 × 10−8 increasing to the north) and an opposing mean wind stress that also drives the near-surface offshore flow. The alongshelf pressure gradient accounts for both the increase in the alongshelf flow with water depth and the geostrophic balance onshore flow in the interior. The increase in the near-bottom offshore flow with water depth is due to the change in the relative magnitude of the contributions from the geostrophic onshore flow that dominates in shallow water and the offshore flow driven by the bottom stress that dominates in deeper water.

Description: This research was funded by Ocean Sciences Division of the National Science Foundation under Grants OCE-820773, OCE-841292, and OCE-848961.

Description: Author Posting. © American Meteorological Society, 2008. 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 38 (2008): 133–145, doi:10.1175/2007JPO3782.1.

Description: Five ice-tethered profilers (ITPs), deployed between 2004 and 2006, have provided detailed potential temperature θ and salinity S profiles from 21 anticyclonic eddy encounters in the central Canada Basin of the Arctic Ocean. The 12–35-m-thick eddies have center depths between 42 and 69 m in the Arctic halocline, and are shallower and less dense than the majority of eddies observed previously in the central Canada Basin. They are characterized by anomalously cold θ and low stratification, and have horizontal scales on the order of, or less than, the Rossby radius of deformation (about 10 km). Maximum azimuthal speeds estimated from dynamic heights (assuming cyclogeostrophic balance) are between 9 and 26 cm s−1, an order of magnitude larger than typical ambient flow speeds in the central basin. Eddy θ–S and potential vorticity properties, as well as horizontal and vertical scales, are consistent with their formation by instability of a surface front at about 80°N that appears in historical CTD and expendable CTD (XCTD) measurements. This would suggest eddy lifetimes longer than 6 months. While the baroclinic instability of boundary currents cannot be ruled out as a generation mechanism, it is less likely since deeper eddies that would originate from the deeper-reaching boundary flows are not observed in the survey region.

Description: The engineering design work for the ITP was initiated by the Cecil H. and Ida M. Green Technology Innovation Program (an internal program at the Woods Hole Oceanographic Institution). Prototype development and construction were funded jointly by the U.S. National Science Foundation (NSF) Oceanographic Technology and Interdisciplinary Coordination Program and Office of Polar Programs (OPP) under Award OCE-0324233. Continued support has been provided by the OPP Arctic Sciences Section under Award ARC-0519899 and internal WHOI funding.

Description: Author Posting. © American Meteorological Society, 2008. 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 38 (2008): 1644-1668, doi:10.1175/2007JPO3829.1.

Description: The mean structure and time-dependent behavior of the shelfbreak jet along the southern Beaufort Sea, and its ability to transport properties into the basin interior via eddies are explored using high-resolution mooring data and an idealized numerical model. The analysis focuses on springtime, when weakly stratified winter-transformed Pacific water is being advected out of the Chukchi Sea. When winds are weak, the observed jet is bottom trapped with a low potential vorticity core and has maximum mean velocities of O(25 cm s−1) and an eastward transport of 0.42 Sv (1 Sv ≡ 106 m3 s−1). Despite the absence of winds, the current is highly time dependent, with relative vorticity and twisting vorticity often important components of the Ertel potential vorticity. An idealized primitive equation model forced by dense, weakly stratified waters flowing off a shelf produces a mean middepth boundary current similar in structure to that observed at the mooring site. The model boundary current is also highly variable, and produces numerous strong, small anticyclonic eddies that transport the shelf water into the basin interior. Analysis of the energy conversion terms in both the mooring data and the numerical model indicates that the eddies are formed via baroclinic instability of the boundary current. The structure of the eddies in the basin interior compares well with observations from drifting ice platforms. The results suggest that eddies shed from the shelfbreak jet contribute significantly to the offshore flux of heat, salt, and other properties, and are likely important for the ventilation of the halocline in the western Arctic Ocean. Interaction with an anticyclonic basin-scale circulation, meant to represent the Beaufort gyre, enhances the offshore transport of shelf water and results in a loss of mass transport from the shelfbreak jet.

Description: This study was supported by the National Science Foundation Office of Polar Programs under Grants 0421904 and 035268 (MS), and by the Office of Naval Research Grant N00014-02-1-0317 (RP and PF). Analysis by AJP was supported by the Office of Naval Research under Grant N00014-97-1-0135 and by the National Science Foundation under Grant OPP-9815303.

Description: Author Posting. © American Meteorological Society, 2008. 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 38 (2008): 1992-2002, doi:10.1175/2008JPO3669.1.

Description: This paper extends A. Bracco and J. Pedlosky’s investigation of the eddy-formation mechanism in the eastern Labrador Sea by including a more realistic depiction of the boundary current. The quasigeostrophic model consists of a meridional, coastally trapped current with three vertical layers. The current configuration and topographic domain are chosen to match, as closely as possible, the observations of the boundary current and the varying topographic slope along the West Greenland coast. The role played by the bottom-intensified component of the boundary current on the formation of the Labrador Sea Irminger Rings is explored. Consistent with the earlier study, a short, localized bottom-trapped wave is responsible for most of the perturbation energy growth. However, for the instability to occur in the three-layer model, the deepest component of the boundary current must be sufficiently strong, highlighting the importance of the near-bottom flow. The model is able to reproduce important features of the observed vortices in the eastern Labrador Sea, including the polarity, radius, rate of formation, and vertical structure. At the time of formation, the eddies have a surface signature as well as a strong circulation at depth, possibly allowing for the transport of both surface and near-bottom water from the boundary current into the interior basin. This work also supports the idea that changes in the current structure could be responsible for the observed interannual variability in the number of Irminger Rings formed.

Description: AB is supported by WHOI unrestricted funds, JP by the National Science Foundation OCE 85108600, and RP by 0450658.

Description: 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): 1103-1121, doi:10.1175/jpo3041.1.

Description: The role of mesoscale oceanic eddies is analyzed in a quasigeostrophic coupled ocean–atmosphere model operating at a large Reynolds number. The model dynamics are characterized by decadal variability that involves nonlinear adjustment of the ocean to coherent north–south shifts of the atmosphere. The oceanic eddy effects are diagnosed by the dynamical decomposition method adapted for nonstationary external forcing. The main effects of the eddies are an enhancement of the oceanic eastward jet separating the subpolar and subtropical gyres and a weakening of the gyres. The flow-enhancing effect is due to nonlinear rectification driven by fluctuations of the eddy forcing. This is a nonlocal process involving generation of the eddies by the flow instabilities in the western boundary current and the upstream part of the eastward jet. The eddies are advected by the mean current to the east, where they backscatter into the rectified enhancement of the eastward jet. The gyre-weakening effect, which is due to the time-mean buoyancy component of the eddy forcing, is a result of the baroclinic instability of the westward return currents. The diagnosed eddy forcing is parameterized in a non-eddy-resolving ocean model, as a nonstationary random process, in which the corresponding parameters are derived from the control coupled simulation. The key parameter of the random process—its variance—is related to the large-scale flow baroclinicity index. It is shown that the coupled model with the non-eddy-resolving ocean component and the parameterized eddies correctly simulates climatology and low-frequency variability of the control eddy-resolving coupled solution.

Description: Funding for this work came from NSF Grants OCE 02-221066 and OCE 03-44094. Additional funding for PB was provided by the U.K. Royal Society Fellowship and by WHOI Grants 27100056 and 52990035.

Description: 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 2009

Description: Internal tides are internal waves of tidal period generated by tidal currents flowing over submarine topography. Tall ridges that are nominally two-dimensional (2-D) are sites of particularly strong generation. The subsequent dissipation of internal tides contributes to ocean mixing, thereby playing an important role in the circulation of the ocean. Strong internal tides can also evolve into internal wave solitons, which affect acoustic communication, offshore structures and submarine navigation. This thesis addresses the generation of internal tides by tall submarine ridges using a combined analytical and experimental approach. The first part of the thesis is an experimental investigation of a pre-existing Green function formulation for internal tide generation by a tall symmetric ridge in a uniform density stratification. A modal decomposition technique was developed to characterize the structure of the experimental wave fields generated by 2D model topographies in a specially configured wave tank. The theory accurately predicts the low mode structure of internal tides, and reasonably predicts the conversion rate of internal tides infinite tidal excursion regimes, for which the emergence of non-linearities was notable in the laboratory. In the second part of the thesis, the Green function method is advanced for asymmetric and multiple ridges in weakly non-uniform stratifications akin to realistic ocean situations. A preliminary investigation in uniform stratification with canonical asymmetric and double ridges reveals asymmetry in the internal tide that can be very sensitive to the geometric configuration. This approach is then used with realistic topography and stratification data to predict the internal tide generated by the ridges at Hawaii and at the Luzon Strait. Despite the assumption of two-dimensionality, there is remarkably good agreement between field data, simulations and the new theory for the magnitude, asymmetry and modal content of the internal tide at these sites. The final part of the thesis investigates the possibility of internal wave attractors in the valley of double-ridge configurations. A one-dimensional map is developed to identify the existence and stability of attractors as a function of the ridge geometry. The Green function method is further advanced to include a viscous correction to balance energy focusing and dissipation along an attracting orbit of internal wave rays, and very good agreement is obtained between experiment and theory, even in the presence of an attractor.

Description: My Ph. D. and the work in this thesis have been generously funded by the National Science Foundation under grants OCE 0645529 and OCE 04-25283 and the Office of Naval Research under grants N00014-08-0390, N00014-05-1-0573 and N00014-09-0282.

Description: 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 August 1987

Description: This thesis consists of two loosely related problems, both of which analyze some consequences of the failure of Sverdrup relation. In the first part, Chapters 2 and 3, the Sverdrup relation is invalidated because substantial flow is obtained at the bottom where topography exists. The eddies play an essential role in transfering momentum vertically from the surface, where the forcing is applied, to the bottom, which is otherwise unforced. If the topography has a structure in the longitudinal direction, then the inviscid theory predicts the occurence of strong jets in the interior of the model ocean. According to the structure of the topography these internal jets can occur in both vertically homeogenous and baroclinic oceans. If the topographic slope changes sign, then one kind of jets is observed both in stratified and in homogeneous oceans. This phenomenon is robust to moderate amounts of dissipation and is not disturbed by the occurrence of recirculating gyres within the basin. If the topographic slope is constant, then another kind of internal jets is observed, and it occurs in stratified models only. I was unable to observe this kind of jets in the presence of weak dissipation. The reason for this failure is twofold: on one hand friction, especially interfacial friction, tends to make the flow more barotropic (and we believe that indeed this is one of the processes that the eddies accomplish in a stratified fluid) and therefore the phenomena that rely strongly on baroclinicity are discouraged. On the other hand, reduction of the dissipation leads to the onset of a strong recirculating, inertial gyre which, although confined in space, affects the global properties of the flow. In the second part of the thesis (Chapters 4 and 5) I developed a simple model of the recirculating, inertial gyre. Again the dynamics of this feature are far from being in Sverdrup balance. In this case inertia is responsible for the failure of Sverdrup relation, together with the eddy field which provides a mean for transfering momentum vertically and laterally into regions away from where the forcing is applied. In this model there is no direct forcing in the recirculation region, and the input of momentum is confined to the boundary currents surrounding the gyre, for example the separated Gulf Stream. One of the results of the recirculation model is the prediction of its transport. It is shown that most of the transport is depth independent, i.e. it can be calculated without detailed knowledge of the density structure of the ocean. It is also shown that the "barotropic" part of the transport increases as the cube of the meridional extent of the gyre.

Description: The thesis work has been supported by a National Foundation grant from the Office of Atmospheric Sciences.