ALBERT

Description: Author Posting. © American Meteorological Society, 2022. 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 52(12),(2022): 3199-3219, https://doi.org/10.1175/jpo-d-22-0009.1.

Description: The abyssal overturning circulation is thought to be primarily driven by small-scale turbulent mixing. Diagnosed water-mass transformations are dominated by rough topography “hotspots,” where the bottom enhancement of mixing causes the diffusive buoyancy flux to diverge, driving widespread downwelling in the interior—only to be overwhelmed by an even stronger upwelling in a thin bottom boundary layer (BBL). These water-mass transformations are significantly underestimated by one-dimensional (1D) sloping boundary layer solutions, suggesting the importance of three-dimensional physics. Here, we use a hierarchy of models to generalize this 1D boundary layer approach to three-dimensional eddying flows over realistically rough topography. When applied to the Mid-Atlantic Ridge in the Brazil Basin, the idealized simulation results are roughly consistent with available observations. Integral buoyancy budgets isolate the physical processes that contribute to realistically strong BBL upwelling. The downward diffusion of buoyancy is primarily balanced by upwelling along the sloping canyon sidewalls and the surrounding abyssal hills. These flows are strengthened by the restratifying effects of submesoscale baroclinic eddies and by the blocking of along-ridge thermal wind within the canyon. Major topographic sills block along-thalweg flows from restratifying the canyon trough, resulting in the continual erosion of the trough’s stratification. We propose simple modifications to the 1D boundary layer model that approximate each of these three-dimensional effects. These results provide local dynamical insights into mixing-driven abyssal overturning, but a complete theory will also require the nonlocal coupling to the basin-scale circulation.

Description: We acknowledge funding support from National Science Foundation Awards 1536515, 1736109, and 2149080. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship Program under Grant 174530.

Description: 2023-05-18

Description: Author Posting. © American Meteorological Society, 2022. 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 52(6), (2022): 1233-1244, https://doi.org/10.1175/jpo-d-21-0223.1.

Description: The Sverdrup relation is the backbone of wind-driven circulation theory; it is a simple relation between the meridional transport of the wind-driven circulation in the upper ocean and the wind stress curl. However, the relation is valid for steady circulation only. In this study, a time-dependent Sverdrup relation is postulated, in which the meridional transport in a time-dependent circulation is the sum of the local wind stress curl term and a time-delayed term representing the effect of the eastern boundary condition. As an example, this time-dependent Sverdrup relation is evaluated through its application to the equatorial circulation in the Indian Ocean, using reanalysis data and a reduced gravity model. Close examination reveals that the southward Somali Current occurring during boreal winter is due to the combination of the local wind stress curl in the Arabian Sea and delayed signals representing the time change of layer thickness at the eastern boundary.

Description: This work is supported by NSFC (41822602, 41976016, 42005035, 42076021), the Strategic Priority Research Program of Chinese Academy of Sciences (XDB42000000, XDA 20060502), Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) (GML2019ZD0306), Guangdong Basic and Applied Basic Research Foundation (2021A1515011534), Youth Innovation Promotion Association CAS, ISEE2021ZD01, and LTOZZ2002. The numerical simulation is supported by the High-Performance Computing Division in the South China Sea Institute of Oceanology.

Description: 2022-11-27

Description: Author Posting. © American Meteorological Society, 2022. 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 52(7), (2022): 1415–1430. https://doi.org/10.1175/JPO-D-21-0147.1.

Description: Strong subinertial variability near a seamount at the Xisha Islands in the South China Sea was revealed by mooring observations from January 2017 to January 2018. The intraseasonal deep flows presented two significant frequency bands, with periods of 9–20 and 30–120 days, corresponding to topographic Rossby waves (TRWs) and deep eddies, respectively. The TRW and deep eddy signals explained approximately 60% of the kinetic energy of the deep subinertial currents. The TRWs at the Ma, Mb, and Mc moorings had 297, 262, and 274 m vertical trapping lengths, and ∼43, 38, and 55 km wavelengths, respectively. Deep eddies were independent from the upper layer, with the largest temperature anomaly being 〉0.4°C. The generation of the TRWs was induced by mesoscale perturbations in the upper layer. The interaction between the cyclonic–anticyclonic eddy pair and the seamount topography contributed to the generation of deep eddies. Owing to the potential vorticity conservation, the westward-propagating tilted interface across the eddy pair squeezed the deep-water column, thereby giving rise to negative vorticity west of the seamount. The strong front between the eddy pair induced a northward deep flow, thereby generating a strong horizontal velocity shear because of lateral friction and enhanced negative vorticity. Approximately 4 years of observations further confirmed the high occurrence of TRWs and deep eddies. TRWs and deep eddies might be crucial for deep mixing near rough topographies by transferring mesoscale energy to small scales.

Description: This work was supported by the National Natural Science Foundation of China (92158204, 91958202, 42076019, 41776036, 91858203), the Open Project Program of State Key Laboratory of Tropical Oceanography (project LTOZZ2001), and Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) (GML2019ZD0304).

Description: 2022-12-16

Description: Author Posting. © American Meteorological Society, 2022. 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 52(1),(2022): 75–97, https://doi.org/10.1175/JPO-D-21-0099.1.

Description: Mesoscale eddies contain the bulk of the ocean’s kinetic energy (KE), but fundamental questions remain on the cross-scale KE transfers linking eddy generation and dissipation. The role of submesoscale flows represents the key point of discussion, with contrasting views of submesoscales as either a source or a sink of mesoscale KE. Here, the first observational assessment of the annual cycle of the KE transfer between mesoscale and submesoscale motions is performed in the upper layers of a typical open-ocean region. Although these diagnostics have marginal statistical significance and should be regarded cautiously, they are physically plausible and can provide a valuable benchmark for model evaluation. The cross-scale KE transfer exhibits two distinct stages, whereby submesoscales energize mesoscales in winter and drain mesoscales in spring. Despite this seasonal reversal, an inverse KE cascade operates throughout the year across much of the mesoscale range. Our results are not incompatible with recent modeling investigations that place the headwaters of the inverse KE cascade at the submesoscale, and that rationalize the seasonality of mesoscale KE as an inverse cascade-mediated response to the generation of submesoscales in winter. However, our findings may challenge those investigations by suggesting that, in spring, a downscale KE transfer could dampen the inverse KE cascade. An exploratory appraisal of the dynamics governing mesoscale–submesoscale KE exchanges suggests that the upscale KE transfer in winter is underpinned by mixed layer baroclinic instabilities, and that the downscale KE transfer in spring is associated with frontogenesis. Current submesoscale-permitting ocean models may substantially understate this downscale KE transfer, due to the models’ muted representation of frontogenesis.

Description: The OSMOSIS experiment was funded by the U.K. Natural Environment Research Council (NERC) through Grants NE/1019999/1 and NE/101993X/1. ACNG acknowledges the support of the Royal Society and the Wolfson Foundation, and XY that of a China Scholarship Council PhD studentship.

Description: Author Posting. © American Meteorological Society, 2021. 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 51(10), (2021): 3235–3252, https://doi.org/10.1175/JPO-D-20-0288.1.

Description: Recent mooring measurements from the Overturning in the Subpolar North Atlantic Program have revealed abundant cyclonic eddies at both sides of Cape Farewell, the southern tip of Greenland. In this study, we present further observational evidence, from both Eulerian and Lagrangian perspectives, of deep cyclonic eddies with intense rotation (ζ/f 〉 1) around southern Greenland and into the Labrador Sea. Most of the observed cyclones exhibit strongest rotation below the surface at 700–1000 dbar, where maximum azimuthal velocities are ~30 cm s−1 at radii of ~10 km, with rotational periods of 2–3 days. The cyclonic rotation can extend to the deep overflow water layer (below 1800 dbar), albeit with weaker azimuthal velocities (~10 cm s−1) and longer rotational periods of about one week. Within the middepth rotation cores, the cyclones are in near solid-body rotation and have the potential to trap and transport water. The first high-resolution hydrographic transect across such a cyclone indicates that it is characterized by a local (both vertically and horizontally) potential vorticity maximum in its middepth core and cold, fresh anomalies in the deep overflow water layer, suggesting its source as the Denmark Strait outflow. Additionally, the propagation and evolution of the cyclonic eddies are illustrated with deep Lagrangian floats, including their detachments from the boundary currents to the basin interior. Taken together, the combined Eulerian and Lagrangian observations have provided new insights on the boundary current variability and boundary–interior exchange over a geographically large scale near southern Greenland, calling for further investigations on the (sub)mesoscale dynamics in the region.

Description: OOI mooring data are based upon work supported by the National Science Foundation under Cooperative Agreement 1743430. S. Zou, A. Bower, and H. Furey gratefully acknowledge the support from the Physical Oceanography Program of the U.S. National Science Foundation Grant OCE-1756361. R.S. Pickart acknowledges support from National Science Foundation Grants OCE-1259618 and OCE-1756361. N. P. Holliday and L. Houpert were supported by NERC programs U.K. OSNAP (NE/K010875) and U.K. OSNAP-Decade (NE/T00858X/1).

Description: Author Posting. © American Meteorological Society, 2021. 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 34(22), (2021): 8971–8987, https://doi.org/10.1175/JCLI-D-20-0610.1.

Description: The impact of increasing Greenland freshwater discharge on the subpolar North Atlantic (SPNA) remains unknown as there are uncertainties associated with the time scales of the Greenland freshwater anomaly (GFWA) in the SPNA. Results from numerical simulations tracking GFWA and an analytical approach are employed to estimate the response time, suggesting that a decadal time scale (13 years) is required for the SPNA to adjust for increasing GFWA. Analytical solutions obtained for a long-lasting increase of freshwater discharge show a non-steady-state response of the SPNA with increasing content of the GFWA. In contrast, solutions for a short-lived pulse of freshwater demonstrate different responses of the SPNA with a rapid increase of freshwater in the domain followed by an exponential decay after the pulse has passed. The derived theoretical relation between time scales shows that residence time scales are time dependent for a non-steady-state case and asymptote the response time scale with time. The residence time of the GFWA deduced from Lagrangian experiments is close to and smaller than the response time, in agreement with the theory. The Lagrangian analysis shows dependence of the residence time on the entrance route of the GFWA and on the depth. The fraction of the GFWA exported through Davis Strait has limited impact on the interior basins, whereas the fraction entering the SPNA from the southwest Greenland shelf spreads into the interior regions. In both cases, the residence time of the GFWA increases with depth demonstrating long persistence of the freshwater anomaly in the subsurface layers.

Description: D. S. Dukhovskoy and E. P. Chassignet were funded by the DOE (Award DE-SC0014378) and HYCOM NOPP (Award N00014-19-1-2674). The HYCOM-CICE simulations were supported by a grant of computer time from the DoD High-Performance Computing Modernization Program at NRL SSC. G. Platov was funded by the RSF N19-17-00154. P. G. Myers was funded by an NSERC Discovery Grant (Grant RGPIN 04357). A. Proshutinsky was funded by FAMOS project (NSF Grant NSF 14-584).

Description: Author Posting. © American Meteorological Society, 2021. 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 51(12),(2021): 3663–3678, https://doi.org/10.1175/JPO-D-21-0058.1.

Description: The wind-driven exchange through complex ridges and islands between marginal seas and the open ocean is studied using both numerical and analytical models. The models are forced by a steady, spatially uniform northward wind stress intended to represent the large-scale, low-frequency wind patterns typical of the seasonal monsoons in the western Pacific Ocean. There is an eastward surface Ekman transport out of the marginal sea and westward geostrophic inflows into the marginal sea. The interaction between the Ekman transport and an island chain produces strong baroclinic flows along the island boundaries with a vertical depth that scales with the ratio of the inertial boundary layer thickness to the baroclinic deformation radius. The throughflows in the gaps are characterized by maximum transport in the center gap and decreasing transports toward the southern and northern tips of the island chain. An extended island rule theory demonstrates that throughflows are determined by the collective balance between viscosity on the meridional boundaries and the eastern side boundary of the islands. The outflowing transport is balanced primarily by a shallow current that enters the marginal sea along its equatorward boundary. The islands can block some direct exchange and result in a wind-driven overturning cell within the marginal sea, but this is compensated for by eastward zonal jets around the southern and northern tips of the island chain. Topography in the form of a deep slope, a ridge, or shallow shelves around the islands alters the current pathways but ultimately is unable to limit the total wind-driven exchange between the marginal sea and the open ocean.

Description: This research is supported in part by the China Scholarship Council (201906330102). H. G. is financially supported by the China Scholarship Council to study at WHOI for 2 years as a guest student. M. A. S. is supported by the National Science Foundation Grant OCE-1922538.

Description: Author Posting. © American Meteorological Society, 2022. 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 52(11), (2022): 2627-2641, https://doi.org/10.1175/jpo-d-22-0090.1.

Description: Changes in dynamic manometric sea level ζm represent mass-related sea level changes associated with ocean circulation and climate. We use twin model experiments to quantify magnitudes and spatiotemporal scales of ζm variability caused by barometric pressure pa loading at long periods (≳1 month) and large scales (≳300km) relevant to Gravity Recovery and Climate Experiment (GRACE) ocean data. Loading by pa drives basin-scale monthly ζm variability with magnitudes as large as a few centimeters. Largest ζm signals occur over abyssal plains, on the shelf, and in marginal seas. Correlation patterns of modeled ζm are determined by continental coasts and H/f contours (H is ocean depth and f is Coriolis parameter). On average, ζm signals forced by pa represent departures of ≲10% and ≲1% from the inverted-barometer effect ζib on monthly and annual periods, respectively. Basic magnitudes, spatial patterns, and spectral behaviors of ζm from the model are consistent with scaling arguments from barotropic potential vorticity conservation. We also compare ζm from the model driven by pa to ζm from GRACE observations. Modeled and observed ζm are significantly correlated across parts of the tropical and extratropical oceans, on shelf and slope regions, and in marginal seas. Ratios of modeled to observed ζm magnitudes are as large as ∼0.2 (largest in the Arctic Ocean) and qualitatively agree with analytical theory for the gain of the transfer function between ζm forced by pa and wind stress. Results demonstrate that pa loading is a secondary but nevertheless important contributor to monthly mass variability from GRACE over the ocean.

Description: The authors acknowledge support from the National Aeronautics and Space Administration through the GRACE Follow-On Science Team (Grant 80NSSC20K0728) and the Sea Level Change Team (Grant 80NSSC20K1241). The contribution from I. F. and O. W. represents research carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (Grant 80NM0018D0004).

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of [publisher] for personal use, not for redistribution. The definitive version was published in Farrar, J. T., Durland, T., Jayne, S. R., & Price, J. F. Long-distance radiation of Rossby Waves from the equatorial current system. Journal of Physical Oceanography, 51(6), (2021): 1947–1966, https://doi.org/10.1175/JPO-D-20-0048.1.

Description: Measurements from satellite altimetry are used to show that sea surface height (SSH) variability throughout much of the North Pacific Ocean is coherent with the SSH signal of the tropical instability waves (TIWs) that result from instabilities of the equatorial currents. This variability has regular phase patterns consistent with freely propagating barotropic Rossby waves radiating energy away from the unstable equatorial currents, and the waves clearly propagate from the equatorial region to at least 30°N. The pattern of SSH variance at TIW frequencies exhibits remarkable patchiness on scales of hundreds of kilometers, which we interpret as being due to the combined effects of wave reflection, refraction, and interference. North of 40°N, more than 6000 km from the unstable equatorial currents, the SSH field remains coherent with the near-equatorial SSH variability, but it is not as clear whether the variability at the higher latitudes is a simple result of barotropic wave radiation from the tropical instability waves. Even more distant regions, as far north as the Aleutian Islands off of Alaska and the Kamchatka Peninsula of eastern Russia, have SSH variability that is significantly coherent with the near-equatorial instabilities. The variability is not well represented in the widely used gridded SSH data product commonly referred to as the AVISO or DUACS product, and this appears to be a result of spatial variations in the filtering properties of the objective mapping scheme.

Description: This work was supported by NASA Grants NNX13AE46G, NNX14AM71G, and NNX17AH54G.

Description: Author Posting. © American Meteorological Society, 2022. 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 the Atmospheric and Oceanic Technology 39(8), (2022): 1183-1198, https://doi.org/10.1175/jtech-d-21-0068.1.

Description: Horizontal kinematic properties, such as vorticity, divergence, and lateral strain rate, are estimated from drifter clusters using three approaches. At submesoscale horizontal length scales O(1–10)km, kinematic properties become as large as planetary vorticity f, but challenging to observe because they evolve on short time scales O(hourstodays). By simulating surface drifters in a model flow field, we quantify the sources of uncertainty in the kinematic property calculations due to the deformation of cluster shape. Uncertainties arise primarily due to (i) violation of the linear estimation methods and (ii) aliasing of unresolved scales. Systematic uncertainties (iii) due to GPS errors, are secondary but can become as large as (i) and (ii) when aspect ratios are small. Ideal cluster parameters (number of drifters, length scale, and aspect ratio) are determined and error functions estimated empirically and theoretically. The most robust method—a two-dimensional, linear least squares fit—is applied to the first few days of a drifter dataset from the Bay of Bengal. Application of the length scale and aspect-ratio criteria minimizes errors (i) and (ii), and reduces the total number of clusters and so computational cost. The drifter-estimated kinematic properties map out a cyclonic mesoscale eddy with a surface, submesoscale fronts at its perimeter. Our analyses suggest methodological guidance for computing the two-dimensional kinematic properties in submesoscale flows, given the recently increasing quantity and quality of drifter observations, while also highlighting challenges and limitations.

Description: This research was supported by the Office of Naval Research (ONR) Departmental Research Initiative ASIRI under Grant N00014-13-1-0451 (SE and AM) and Grant N00014-13-1-0477 (VH and LC). The authors thank the captain and crew of the R/V Roger Revelle, and Andrew Lucas with the Multiscale Ocean Dynamics group at the Scripps Institution for Oceanography for providing the FastCTD data collected in 2015, which was supported by ONR Grant N00014-13-1-0489, as well as Eric D’Asaro for helpful discussions and Lance Braasch for assistance with the drifter dataset. AM and SE further thank NSF (Grant OCE-I434788) and ONR (Grant N00014-16-1-2470) for support. VH and LC were additionally supported by ONR Grants N00014-15-1-2286, N00014-14-1-0183, N00014-19-1-26-91 and NOAA Global Drifter Program (GDP) Grant NA15OAR4320071.

Description: 2023-02-01

Description: Author Posting. © American Meteorological Society, 2022. 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 35(17), (2022): 5465-5482, https://doi.org/10.1175/jcli-d-21-0671.1.

Description: Understanding the contribution of ocean circulation to glacial–interglacial climate change is a major focus of paleoceanography. Specifically, many have tried to determine whether the volumes and depths of Antarctic- and North Atlantic–sourced waters in the deep ocean changed at the Last Glacial Maximum (LGM; ∼22–18 kyr BP) when atmospheric pCO2 concentrations were 100 ppm lower than the preindustrial. Measurements of sedimentary geochemical proxies are the primary way that these deep ocean structural changes have been reconstructed. However, the main proxies used to reconstruct LGM Atlantic water mass geometry provide conflicting results as to whether North Atlantic–sourced waters shoaled during the LGM. Despite this, a number of idealized modeling studies have been advanced to describe the physical processes resulting in shoaled North Atlantic waters. This paper aims to critically assess the approaches used to determine LGM Atlantic circulation geometry and lay out best practices for future work. We first compile existing proxy data and paleoclimate model output to deduce the processes responsible for setting the ocean distributions of geochemical proxies in the LGM Atlantic Ocean. We highlight how small-scale mixing processes in the ocean interior can decouple tracer distributions from the large-scale circulation, complicating the straightforward interpretation of geochemical tracers as proxies for water mass structure. Finally, we outline promising paths toward ascertaining the LGM circulation structure more clearly and deeply.

Description: S.K.H. was supported by the Investment in Science Fund at WHOI and the John E. and Anne W. Sawyer Endowed Fund in Support of Scientific Staff. F.J.P. was supported by a Stanback Postdoctoral Fellowship at Caltech.

Description: 2023-02-01

Description: Author Posting. © American Meteorological Society, 2020. 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 50(9), (2020): 2491-2506, doi:10.1175/JPO-D-20-0056.1.

Description: An idealized two-layer shallow water model is applied to the study of the dynamics of the Arctic Ocean halocline. The model is forced by a surface stress distribution reflective of the observed wind stress pattern and ice motion and by an inflow representing the flow of Pacific Water through Bering Strait. The model reproduces the main elements of the halocline circulation: an anticyclonic Beaufort Gyre in the western basin (representing the Canada Basin), a cyclonic circulation in the eastern basin (representing the Eurasian Basin), and a Transpolar Drift between the two gyres directed from the upwind side of the basin to the downwind side of the basin. Analysis of the potential vorticity budget shows a basin-averaged balance primarily between potential vorticity input at the surface and dissipation at the lateral boundaries. However, advection is a leading-order term not only within the anticyclonic and cyclonic gyres but also between the gyres. This means that the eastern and western basins are dynamically connected through the advection of potential vorticity. Both eddy and mean fluxes play a role in connecting the regions of potential vorticity input at the surface with the opposite gyre and with the viscous boundary layers. These conclusions are based on a series of model runs in which forcing, topography, straits, and the Coriolis parameter were varied.

Description: This study was supported by National Science Foundation Grant OPP-1822334. Comments and suggestions from two anonymous referees greatly helped to improve the paper.

Description: 2021-02-17

Description: Author Posting. © American Meteorological Society, 2020. 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 50(11), (2020): 3235–3251, https://doi.org/10.1175/JPO-D-20-0095.1.

Description: The dense outflow through Denmark Strait is the largest contributor to the lower limb of the Atlantic meridional overturning circulation, yet a description of the full velocity field across the strait remains incomplete. Here we analyze a set of 22 shipboard hydrographic–velocity sections occupied along the Látrabjarg transect at the Denmark Strait sill, obtained over the time period 1993–2018. The sections provide the first complete view of the kinematic components at the sill: the shelfbreak East Greenland Current (EGC), the combined flow of the separated EGC, and the North Icelandic Jet (NIJ), and the northward-flowing North Icelandic Irminger Current (NIIC). The total mean transport of overflow water is 3.54 ± 0.29 Sv (1 Sv ≡ 106 m3 s−1), comparable to previous estimates. The dense overflow is partitioned in terms of water mass constituents and flow components. The mean transports of the two types of overflow water—Atlantic-origin Overflow Water and Arctic-origin Overflow Water—are comparable in Denmark Strait, while the merged NIJ–separated EGC transports 55% more water than the shelfbreak EGC. A significant degree of water mass exchange takes place between the branches as they converge in Denmark Strait. There are two dominant time-varying configurations of the flow that are characterized as a cyclonic state and a noncyclonic state. These appear to be wind-driven. A potential vorticity analysis indicates that the flow through Denmark Strait is subject to symmetric instability. This occurs at the top of the overflow layer, implying that the mixing/entrainment process that modifies the overflow water begins at the sill.

Description: Funding for the study was provided by National Science Foundation (NSF) Grants OCE-1259618, OCE-1756361, and OCE-1558742. The German research cruises were financially supported through various EU Projects (e.g. THOR, NACLIM) and national projects (most recently TRR 181 “Energy Transfer in Atmosphere and Ocean” funded by the German Research Foundation and RACE II “Regional Atlantic Circulation and Global Change” funded by the German Federal Ministry for Education and Research). GWKM acknowledges the support of the Natural Sciences and Engineering Research Council of Canada. LP was supported by NSF Grant OCE-1657870.

Description: Author Posting. © American Meteorological Society, 2021. 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 51(1), (2021): 247-266, https://doi.org/10.1175/JPO-D-20-0098.1.

Description: This study focuses on mechanisms of shelf valley bathymetry affecting the spread of riverine freshwater in the nearshore region. In the context of Changjiang River, a numerical model is used with different no-tide idealized configurations to simulate development of unforced river plumes over a sloping bottom, with and without a shelf valley off the estuary mouth. All simulated freshwater plumes are surface-trapped with continuously growing bulges near the estuary mouth and narrow coastal currents downstream. The simulations indicate that a shelf valley tends to compress the bulge along the direction of the valley long axis and modify the incident angle of the bulge flow impinging toward the coast, which then affects the strength of the coastal current. The bulge compression results from geostrophic adjustment and isobath-following tendency of the depth-averaged flow in the bulge region. Generally, the resulting change in the direction of the bulge impinging flow enhances down-shelf momentum advection and freshwater delivery into the coastal current. Sensitivity simulations with altered river discharges Q, Coriolis parameter, shelf bottom slope, valley geometry, and ambient stratification show that enhancement of down-shelf freshwater transport in the coastal current, ΔQc, increases with increasing valley depth within the bulge region and decreasing slope Burger number of the ambient shelf. Assuming potential vorticity conservation, a scaling formula of ΔQc/Q is developed, and it agrees well with results of the sensitivity simulations. Mechanisms of valley influences on unforced river plumes revealed here will help future studies of topographic influence on river plumes under more realistic conditions.

Description: This work is conducted by Canbo Xiao and Weifeng (Gordon) Zhang during CX’s one-year visit at Woods Hole Oceanographic Institution (WHOI) in 2018–19. CX was supported by China Scholarship Council.

Description: Author Posting. © American Meteorological Society, 2021. 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 51(3), (2021): 955–973, https://doi.org/10.1175/JPO-D-20-0240.1.

Description: Fresh Arctic waters flowing into the Atlantic are thought to have two primary fates. They may be mixed into the deep ocean as part of the overturning circulation, or flow alongside regions of deep water formation without impacting overturning. Climate models suggest that as increasing amounts of freshwater enter the Atlantic, the overturning circulation will be disrupted, yet we lack an understanding of how much freshwater is mixed into the overturning circulation’s deep limb in the present day. To constrain these freshwater pathways, we build steady-state volume, salt, and heat budgets east of Greenland that are initialized with observations and closed using inverse methods. Freshwater sources are split into oceanic Polar Waters from the Arctic and surface freshwater fluxes, which include net precipitation, runoff, and ice melt, to examine how they imprint the circulation differently. We find that 65 mSv (1 Sv ≡ 106 m3 s−1) of the total 110 mSv of surface freshwater fluxes that enter our domain participate in the overturning circulation, as do 0.6 Sv of the total 1.2 Sv of Polar Waters that flow through Fram Strait. Based on these results, we hypothesize that the overturning circulation is more sensitive to future changes in Arctic freshwater outflow and precipitation, while Greenland runoff and iceberg melt are more likely to stay along the coast of Greenland.

Description: We gratefully acknowledge the U.S. National Science Foundation: this work was supported by Grants OCE-1258823, OCE-1756272, OCE-1948335, and OCE-2038481. L.H.S. thanks the U.S. Norway Fulbright Foundation for the Norwegian Arctic Chair Grant 2019-20 that made the visit to Scripps Institution of Oceanography possible. N.P.H. acknowledges support by the U.K. Natural Environment Research Council (NERC) National Capability program CLASS (NE/R015953/1), and Grants U.K.-OSNAP (NE/K010875/1, NE/K010875/2) and U.K.-OSNAP Decade (NE/T00858X/1). We acknowledge the World Climate Research Programme, which, through its Working Group on Coupled Modelling, coordinated and promoted CMIP6.

Description: Author Posting. © American Meteorological Society, 2020. 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 50(9),(2020): 2797-2814, https://doi.org/10.1175/JPO-D-19-0326.1.

Description: Hydrographic measurements recently acquired along the thalweg of the Lifamatola Passage combined with historical moored velocity measurements immediately downstream of the sill are used to study the hydraulics, transport, mixing, and entrainment in the dense overflow. The observations suggest that the mean overflow is nearly critical at the mooring site, suggesting that a weir formula may be appropriate for estimating the overflow transport. Our assessment suggests that the weir formulas corresponding to a rectangular, triangular, or parabolic cross section all result in transports very close to the observation, suggesting their potential usage in long-term monitoring of the overflow transport or parameterizing the transport in numerical models. Analyses also suggest that deep signals within the overflow layer are blocked by the shear flow from propagating upstream, whereas the shallow wave modes of the full-depth continuously stratified flow are able to propagate upstream from the Banda Sea into the Maluku Sea. Strong mixing is found immediately downstream of the sill crest, with Thorpe-scale-based estimates of the mean dissipation rate within the overflow up to 1.1 × 10−7 W kg−1 and the region-averaged diapycnal diffusivity within the downstream overflow in the range of 2.3 × 10−3 to 10.1 × 10−3 m2 s−1. Mixing in the Lifamatola Passage results in 0.6–1.2-Sv (1 Sv ≡ 106 m3 s−1) entrainment transport added to the overflow, enhancing the deep-water renewal in the Banda Sea. A bulk diffusivity coefficient estimated in the deep Banda Sea yields 1.6 × 10−3 ± 5 × 10−4 m2 s−1, with an associated downward turbulent heat flux of 9 W m−2.

Description: This study is supported by NSFC (91858204), the CAS Strategic Priority Research Program (XDB42000000), NSFC(41720104008, 41421005, 41876025), QMSNL (2018SDKJ0104-02), and the Shandong Provincial projects (U1606402). L. Pratt was supported by the U.S. NSF Grant OCE-1657870.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Liang, Y., Kwon, Y., & Frankignoul, C. Autumn Arctic Pacific sea ice dipole as a source of predictability for subsequent spring Barents Sea ice condition. Journal of Climate, 34(2), (2021): 787-804, https://doi.org/10.1175/JCLI-D-20-0172.1.

Description: This study uses observational and reanalysis datasets in 1980–2016 to show a close connection between a boreal autumn sea ice dipole in the Arctic Pacific sector and sea ice anomalies in the Barents Sea (BS) during the following spring. The September–October Arctic Pacific sea ice dipole variations are highly correlated with the subsequent April–May BS sea ice variations (r = 0.71). The strong connection between the regional sea ice variabilities across the Arctic uncovers a new source of predictability for spring BS sea ice prediction at 7-month lead time. A cross-validated linear regression prediction model using the Arctic Pacific sea ice dipole with 7-month lead time is demonstrated to have significant prediction skills with 0.54–0.85 anomaly correlation coefficients. The autumn sea ice dipole, manifested as sea ice retreat in the Beaufort and Chukchi Seas and expansion in the East Siberian and Laptev Seas, is primarily forced by preceding atmospheric shortwave anomalies from late spring to early autumn. The spring BS sea ice increases are mostly driven by an ocean-to-sea ice heat flux reduction in preceding months, associated with reduced horizontal ocean heat transport into the BS. The dynamical linkage between the two regional sea ice anomalies is suggested to involve positive stratospheric polar cap anomalies during autumn and winter, with its center slowly moving toward Greenland. The migration of the stratospheric anomalies is followed in midwinter by a negative North Atlantic Oscillation–like pattern in the troposphere, leading to reduced ocean heat transport into the BS and sea ice extent increase.

Description: This study is supported by NSF’s Office of Polar Programs (Grant 1736738). We also acknowledge support by the Blue-Action project (European Union’s Horizon 2020 research and innovation programme, Grant 727852).

Description: Author Posting. © American Meteorological Society, 2020. 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 33(9), (2020): 3845-3862, doi:10.1175/JCLI-D-19-0215.1.

Description: The latitudinal structure of the Atlantic meridional overturning circulation (AMOC) variability in the North Atlantic is investigated using numerical results from three ocean circulation simulations over the past four to five decades. We show that AMOC variability south of the Labrador Sea (53°N) to 25°N can be decomposed into a latitudinally coherent component and a gyre-opposing component. The latitudinally coherent component contains both decadal and interannual variabilities. The coherent decadal AMOC variability originates in the subpolar region and is reflected by the zonal density gradient in that basin. It is further shown to be linked to persistent North Atlantic Oscillation (NAO) conditions in all three models. The interannual AMOC variability contained in the latitudinally coherent component is shown to be driven by westerlies in the transition region between the subpolar and the subtropical gyre (40°–50°N), through significant responses in Ekman transport. Finally, the gyre-opposing component principally varies on interannual time scales and responds to local wind variability related to the annual NAO. The contribution of these components to the total AMOC variability is latitude-dependent: 1) in the subpolar region, all models show that the latitudinally coherent component dominates AMOC variability on interannual to decadal time scales, with little contribution from the gyre-opposing component, and 2) in the subtropical region, the gyre-opposing component explains a majority of the interannual AMOC variability in two models, while in the other model, the contributions from the coherent and the gyre-opposing components are comparable. These results provide a quantitative decomposition of AMOC variability across latitudes and shed light on the linkage between different AMOC variability components and atmospheric forcing mechanisms.

Description: The authors gratefully acknowledge support from the Physical Oceanography Program of the U.S. National Science Foundation (Awards OCE-1756143 and OCE-1537136) and the Climate Program Office of the National Oceanic and Atmospheric Administration (Award NA15OAR4310088). Gratitude is extended to Claus Böning and Arne Biastoch who shared ORCA025 output. S. Zou thanks F. Li, M. Buckley, and L. Li for helpful discussions. We also thank three anonymous reviewers for helpful suggestions.

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Hahn, L. C., Storelvmo, T., Hofer, S., Parfitt, R., & Ummenhofer, C. C. Importance of Orography for Greenland cloud and melt response to atmospheric blocking. Journal of Climate, 33(10), (2020): 4187-4206, doi:10.1175/JCLI-D-19-0527.1.

Description: More frequent high pressure conditions associated with atmospheric blocking episodes over Greenland in recent decades have been suggested to enhance melt through large-scale subsidence and cloud dissipation, which allows more solar radiation to reach the ice sheet surface. Here we investigate mechanisms linking high pressure circulation anomalies to Greenland cloud changes and resulting cloud radiative effects, with a focus on the previously neglected role of topography. Using reanalysis and satellite data in addition to a regional climate model, we show that anticyclonic circulation anomalies over Greenland during recent extreme blocking summers produce cloud changes dependent on orographic lift and descent. The resulting increased cloud cover over northern Greenland promotes surface longwave warming, while reduced cloud cover in southern and marginal Greenland favors surface shortwave warming. Comparison with an idealized model simulation with flattened topography reveals that orographic effects were necessary to produce area-averaged decreasing cloud cover since the mid-1990s and the extreme melt observed in the summer of 2012. This demonstrates a key role for Greenland topography in mediating the cloud and melt response to large-scale circulation variability. These results suggest that future melt will depend on the pattern of circulation anomalies as well as the shape of the Greenland Ice Sheet.

Description: This research was supported by the Woods Hole Oceanographic Institution Summer Student Fellow program, by the U.S. National Science Foundation under AGS-1355339 to C.C.U., and by the European Research Council through Grant 758005.

Description: Author Posting. © American Meteorological Society, 2020. 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 50(3), (2020): 679-694, doi:10.1175/JPO-D-19-0218.1.

Description: The zonally integrated flow in a basin can be separated into the divergent/nondivergent parts, and a uniquely defined meridional overturning circulation (MOC) can be calculated. For a basin with significant volume exchange at zonal open boundaries, this method is competent in removing the components associated with the nonzero source terms due to zonal transports at open boundaries. This method was applied to the zonally integrated flow in the Indian Ocean basin extended all the way to the Antarctic by virtue of the ECCO dataset. The contributions due to two major zonal flow systems at open boundaries, the Indonesian Throughflow (ITF) and the Antarctic Circumpolar Current (ACC), were well separated from the rotational flow component, and a nondivergent overturning circulation pattern was identified. Comparisons with previous studies on the MOC of the Indian Ocean in different seasons showed overall consistency but with refinements in details to the south of the entry of the ITF, reflecting the influence of ITF on the MOC pattern in the domain. Other options of decomposition are also examined.

Description: LH was supported by the National Basic Research Program of China through Grant 2019YFA0606703 and “The Fundamental Research Funds of Shandong University” (2019GN051). The authors thank the anonymous reviewers and the editor for their constructive comments. Code availability: The Matlab code that performs the decomposition and produces some figures in this paper is available at https://github.com/lei-han-SDU/IMOC/.

Description: 2020-09-02

Description: Author Posting. © American Meteorological Society, 2020. 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 50(4), (2020): 887-905, doi:10.1175/JPO-D-19-0110.1.

Description: The Equatorial Undercurrent (EUC) encounters the Galápagos Archipelago on the equator as it flows eastward across the Pacific. The impact of the Galápagos Archipelago on the EUC in the eastern equatorial Pacific remains largely unknown. In this study, the path of the EUC as it reaches the Galápagos Archipelago is measured directly using high-resolution observations obtained by autonomous underwater gliders. Gliders were deployed along three lines that define a closed region with the Galápagos Archipelago as the eastern boundary and 93°W from 2°S to 2°N as the western boundary. Twelve transects were simultaneously occupied along the three lines during 52 days in April–May 2016. Analysis of individual glider transects and average sections along each line show that the EUC splits around the Galápagos Archipelago. Velocity normal to the transects is used to estimate net horizontal volume transport into the volume. Downward integration of the net horizontal transport profile provides an estimate of the time- and areal-averaged vertical velocity profile over the 52-day time period. Local maxima in vertical velocity occur at depths of 25 and 280 m with magnitudes of (1.7 ± 0.6) × 10−5 m s−1 and (8.0 ± 1.6) × 10−5 m s−1, respectively. Volume transport as a function of salinity indicates that water crossing 93°W south (north) of 0.4°S tends to flow around the south (north) side of the Galápagos Archipelago. Comparisons are made between previous observational and modeling studies with differences attributed to effects of the strong 2015/16 El Niño event, the annual cycle of local winds, and varying longitudes between studies of the equatorial Pacific.

Description: This work was supported by National Science Foundation (Grants OCE-1232971 and OCE-1233282) and the NASA Earth and Space Science Fellowship Program (Grant 80NSSC17K0443).

Description: Author Posting. © American Meteorological Society, 2020. 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 50(4), (2020): 1045-1064, doi:10.1175/JPO-D-19-0137.1.

Description: Three simulations of the circulation in the Gulf of Mexico (the “Gulf”) using different numerical general circulation models are compared with results of recent large-scale observational campaigns conducted throughout the deep (〉1500 m) Gulf. Analyses of these observations have provided new understanding of large-scale mean circulation features and variability throughout the deep Gulf. Important features include cyclonic flow along the continental slope, deep cyclonic circulation in the western Gulf, a counterrotating pair of cells under the Loop Current region, and a cyclonic cell to the south of this pair. These dominant circulation features are represented in each of the ocean model simulations, although with some obvious differences. A striking difference between all the models and the observations is that the simulated deep eddy kinetic energy under the Loop Current region is generally less than one-half of that computed from observations. A multidecadal integration of one of these numerical simulations is used to evaluate the uncertainty of estimates of velocity statistics in the deep Gulf computed from limited-length (4 years) observational or model records. This analysis shows that the main deep circulation features identified from the observational studies appear to be robust and are not substantially impacted by variability on time scales longer than the observational records. Differences in strengths and structures of the circulation features are identified, however, and quantified through standard error analysis of the statistical estimates using the model solutions.

Description: This work was supported by the Gulf Research Program of the National Academy of Sciences under Awards 2000006422 and 2000009966. The content is solely the responsibility of the authors and does not necessarily represent the official views of the Gulf Research Program or the National Academy of Sciences. The authors acknowledge the GLORYS project for providing the ocean reanalysis data used in the ROMS simulation. GLORYS is jointly conducted by MERCATOR OCEAN, CORIOLIS, and CNRS/INSU. Installation, recovery, data acquisition, and processing of the CANEK group current-meter moorings were possible because of CICESE-PetróleosMexicanos Grant PEP-CICESE 428229851 and the dedicated work of the crew of the B/O Justo Sierra and scientists of the CANEK group. The authors thank Dr. Aljaz Maslo, CICESE, for assistance with analysis of model data. The Bureau of Ocean Energy Management (BOEM), U.S. Dept. of the Interior, provided funding for the Lagrangian Study of the Deep Circulation in the Gulf of Mexico and the Observations and Dynamics of the Loop Current study. HYCOM simulation data are available from the HYCOM data server (https://www.hycom.org/data/goml0pt04/expt-02pt2), MITgcm data are available from the ECCO data server (http://ecco.ucsd.edu/gom_results2.html), and the ROMS simulation data are available from GRIIDC (NA.x837.000:0001).

Description: Author Posting. © American Meteorological Society, 2019. 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 49(11), (2019): 2781-2797, doi: 10.1175/JPO-D-19-0111.1.

Description: To ground truth the large-scale dynamical balance of the North Atlantic subtropical gyre with observations, a barotropic vorticity budget is constructed in the ECCO state estimate and compared with hydrographic observations and wind stress data products. The hydrographic dataset at the center of this work is the A22 WOCE section, which lies along 66°W and creates a closed volume with the North and South American coasts to its west. The planetary vorticity flux across A22 is quantified, providing a metric for the net meridional flow in the western subtropical gyre. The wind stress forcing over the subtropical gyre to the west and east of the A22 section is calculated from several wind stress data products. These observational budget terms are found to be consistent with an approximate barotropic Sverdrup balance in the eastern subtropical gyre and are on the same order as budget terms in the ECCO state estimate. The ECCO vorticity budget is closed by bottom pressure torques in the western subtropical gyre, which is consistent with previous studies. In sum, the analysis provides observational ground truth for the North Atlantic subtropical vorticity balance and explores the seasonal variability of this balance for the first time using the ECCO state estimate. This balance is found to hold on monthly time scales in ECCO, suggesting that the integrated subtropical gyre responds to forcing through fast barotropic adjustment.

Description: We thank Alonso Hernández-Guerra, M. Dolores Pérez-Hernández, and María Casanova-Masjoan for providing the inverse model results from Casanova-Masjoan et al. (2018). The A22 section is part of the WOCE/CLIVAR observing effort, with all data available at http://cchdo.ucsd.edu/. We thank Carl Wunsch, Patrick Heimbach, Chris Hill, and Diana Lees Spiegel for their assistance with the ECCO fields. The state estimates were provided by the ECCO Consortium for Estimating the Circulation and Climate of the Ocean funded by the National Oceanographic Partnership Program (NOPP) and can be downloaded at http://www.ecco-group.org/products.htm. The citable URL for the ECCO version 4 release 2 product is http://hdl.handle.net/1721.1/102062. We are grateful to Joseph Pedlosky and Glenn Flierl for their comments on an earlier version of this work. IALB and JMT were supported financially by U.S. NSF Grants OCE-0726720, 1332667, and 1332834. MS was supported by the U.S. NASA Sea Level Change Team (Contract NNX14AJ51G) and through the ECCO Consortium funding via the Jet Propulsion Laboratory. We thank two anonymous reviewers, whose thoughtful comments led to improvements.

Description: 2020-04-17

Description: Author Posting. © American Meteorological Society, 2019. 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 49(12), (2019): 3061-3068, doi: 10.1175/JPO-D-18-0172.1.

Description: The calculation of energy flux in coastal trapped wave modes is reviewed in the context of tidal energy pathways near the coast. The significant barotropic pressures and currents associated with coastal trapped wave modes mean that large errors in estimating the wave flux are incurred if only the baroclinic component is considered. A specific example is given showing that baroclinic flux constitutes only 10% of the flux in a mode-1 wave for a reasonable choice of stratification and bathymetry. The interpretation of baroclinic energy flux and barotropic-to-baroclinic conversion at the coast is discussed: in contrast to the open ocean, estimates of baroclinic energy flux do not represent a wave energy flux; neither does conversion represent the scattering of energy from the tidal Kelvin wave to higher modes.

Description: This work was supported by the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Weston Howland Jr. Postdoctoral Scholarship, and by NSF under Grant OCE-1756781. I am grateful to K. Brink for the many useful conversations that contributed to this work and to J. Toole for providing detailed comments on an early version of this paper. The comments of three anonymous reviewers were very helpful in improving this paper.

Description: 2020-06-03

Description: Author Posting. © American Meteorological Society, 2020. 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 50(3), (2020): 715-726, doi:10.1175/JPO-D-19-0021.1.

Description: Closing the overturning circulation of bottom water requires abyssal transformation to lighter densities and upwelling. Where and how buoyancy is gained and water is transported upward remain topics of debate, not least because the available observations generally show downward-increasing turbulence levels in the abyss, apparently implying mean vertical turbulent buoyancy-flux divergence (densification). Here, we synthesize available observations indicating that bottom water is made less dense and upwelled in fracture zone valleys on the flanks of slow-spreading midocean ridges, which cover more than one-half of the seafloor area in some regions. The fracture zones are filled almost completely with water flowing up-valley and gaining buoyancy. Locally, valley water is transformed to lighter densities both in thin boundary layers that are in contact with the seafloor, where the buoyancy flux must vanish to match the no-flux boundary condition, and in thicker layers associated with downward-decreasing turbulence levels below interior maxima associated with hydraulic overflows and critical-layer interactions. Integrated across the valley, the turbulent buoyancy fluxes show maxima near the sidewall crests, consistent with net convergence below, with little sensitivity of this pattern to the vertical structure of the turbulence profiles, which implies that buoyancy flux convergence in the layers with downward-decreasing turbulence levels dominates over the divergence elsewhere, accounting for the net transformation to lighter densities in fracture zone valleys. We conclude that fracture zone topography likely exerts a controlling influence on the transformation and upwelling of bottom water in many areas of the global ocean.

Description: The data used in this study were collected in the context of several projects funded by the U.S. National Science Foundation (NSF), in particular BBTRE (OCE-9415589 and OCE-9415598) and DoMORE (OCE-1235094). Funding for the analysis was provided as part of the NSF DoMORE and DECIMAL (OCE-1735618) projects. Author Ijichi is a Japan Society for the Promotion of Science (JSPS) Overseas Research Fellow. Comments on an early draft of this paper by Jim Ledwell and Bryan Kaiser, as well as topical discussions with Jörn Callies and Trevor McDougall, are gratefully acknowledged. The paper was greatly improved during the review process, in particular because of the critical comments from one of the two anonymous reviewers.

Description: Author Posting. © American Meteorological Society, 2019. 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 32(24), (2019): 8449-8463, doi: 10.1175/JCLI-D-19-0252.1.

Description: A theory for the mean ice thickness and the Transpolar Drift in the Arctic Ocean is developed. Asymptotic expansions of the ice momentum and thickness equations are used to derive analytic expressions for the leading-order ice thickness and velocity fields subject to wind stress forcing and heat loss to the atmosphere. The theory is most appropriate for the eastern and central Arctic, but not for the region of the Beaufort Gyre subject to anticyclonic wind stress curl. The scale analysis reveals two distinct regimes: a thin ice regime in the eastern Arctic and a thick ice regime in the western Arctic. In the eastern Arctic, the ice drift is controlled by a balance between wind and ocean drag, while the ice thickness is controlled by heat loss to the atmosphere. In contrast, in the western Arctic, the ice thickness is determined by a balance between wind and internal ice stress, while the drift is indirectly controlled by heat loss to the atmosphere. The southward flow toward Fram Strait is forced by the across-wind gradient in ice thickness. The basic predictions for ice thickness, heat loss, ice volume, and ice export from the theory compare well with an idealized, coupled ocean–ice numerical model over a wide range of parameter space. The theory indicates that increasing atmospheric temperatures or wind speed result in a decrease in maximum ice thickness and ice volume. Increasing temperatures also result in a decrease in heat loss to the atmosphere and ice export through Fram Strait, while increasing winds drive increased heat loss and ice export.

Description: MAS was supported by the National Science Foundation under Grant OPP-1822334. Comments and suggestions from Michael Steele, Gianluca Meneghello, and an anonymous reviewer helped to clarify the work.

Description: 2020-05-15

Description: Author Posting. © American Meteorological Society, 2020. 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 50(2), (2020): 455-469, doi:10.1175/JPO-D-19-0190.1.

Description: The mechanisms by which time-dependent wind stress anomalies at midlatitudes can force variability in the meridional overturning circulation at low latitudes are explored. It is shown that winds are effective at forcing remote variability in the overturning circulation when forcing periods are near the midlatitude baroclinic Rossby wave basin-crossing time. Remote overturning is required by an imbalance in the midlatitude mass storage and release resulting from the dependence of the Rossby wave phase speed on latitude. A heuristic theory is developed that predicts the strength and frequency dependence of the remote overturning well when compared to a two-layer numerical model. The theory indicates that the variable overturning strength, relative to the anomalous Ekman transport, depends primarily on the ratio of the meridional spatial scale of the anomalous wind stress curl to its latitude. For strongly forced systems, a mean deep western boundary current can also significantly enhance the overturning variability at all latitudes. For sufficiently large thermocline displacements, the deep western boundary current alternates between interior and near-boundary pathways in response to fluctuations in the wind, leading to large anomalies in the volume of North Atlantic Deep Water stored at midlatitudes and in the downstream deep western boundary current transport.

Description: MAS and DN were supported by the National Science Foundation under Grant OCE-1634468.

Description: 2020-11-10

Description: Author Posting. © American Meteorological Society, 2019. 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 49(7), (2019): 1889-1904, doi:10.1175/JPO-D-19-0053.1.

Description: A high-resolution numerical model, together with in situ and satellite observations, is used to explore the nature and dynamics of the dominant high-frequency (from one day to one week) variability in Denmark Strait. Mooring measurements in the center of the strait reveal that warm water “flooding events” occur, whereby the North Icelandic Irminger Current (NIIC) propagates offshore and advects subtropical-origin water northward through the deepest part of the sill. Two other types of mesoscale processes in Denmark Strait have been described previously in the literature, known as “boluses” and “pulses,” associated with a raising and lowering of the overflow water interface. Our measurements reveal that flooding events occur in conjunction with especially pronounced pulses. The model indicates that the NIIC hydrographic front is maintained by a balance between frontogenesis by the large-scale flow and frontolysis by baroclinic instability. Specifically, the temperature and salinity tendency equations demonstrate that the eddies act to relax the front, while the mean flow acts to sharpen it. Furthermore, the model reveals that the two dense water processes—boluses and pulses (and hence flooding events)—are dynamically related to each other and tied to the meandering of the hydrographic front in the strait. Our study thus provides a general framework for interpreting the short-time-scale variability of Denmark Strait Overflow Water entering the Irminger Sea.

Description: MAS was supported by the National Science Foundation (NSF) under Grants OCE-1558742 and OCE-1534618. RSP, PL, and DM were supported by NSF under Grants OCE-1558742 and OCE-1259618. WJvA was supported by the Helmholtz Infrastructure Initiative FRAM. TWNH and MA were supported by NSF under Grants OCE-1633124 and OCE-118123.

Description: 2020-07-01

Description: © The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Semper, S., Vage, K., Pickart, R. S., Valdimarsson, H., Torres, D. J., & Jonsson, S. The emergence of the North Icelandic Jet and its evolution from northeast Iceland to Denmark Strait. Journal of Physical Oceanography, 49(10), (2019): 2499-2521, doi:10.1175/JPO-D-19-0088.1.

Description: The North Icelandic Jet (NIJ) is an important source of dense water to the overflow plume passing through Denmark Strait. The properties, structure, and transport of the NIJ are investigated for the first time along its entire pathway following the continental slope north of Iceland, using 13 hydrographic/velocity surveys of high spatial resolution conducted between 2004 and 2018. The comprehensive dataset reveals that the current originates northeast of Iceland and increases in volume transport by roughly 0.4 Sv (1 Sv ≡ 106 m3 s−1) per 100 km until 300 km upstream of Denmark Strait, at which point the highest transport is reached. The bulk of the NIJ transport is confined to a small area in Θ–S space centered near −0.29° ± 0.16°C in Conservative Temperature and 35.075 ± 0.006 g kg−1 in Absolute Salinity. While the hydrographic properties of this transport mode are not significantly modified along the NIJ’s pathway, the transport estimates vary considerably between and within the surveys. Neither a clear seasonal signal nor a consistent link to atmospheric forcing was found, but barotropic and/or baroclinic instability is likely active in the current. The NIJ displays a double-core structure in roughly 50% of the occupations, with the two cores centered at the 600- and 800-m isobaths, respectively. The transport of overflow water 300 km upstream of Denmark Strait exceeds 1.8 ± 0.3 Sv, which is substantially larger than estimates from a year-long mooring array and hydrographic/velocity surveys closer to the strait, where the NIJ merges with the separated East Greenland Current. This implies a more substantial contribution of the NIJ to the Denmark Strait overflow plume than previously envisaged.

Description: Six different research vessels were involved in the collection of the data used in this study: RRS James Clark Ross, R/V Knorr, R/V Bjarni Sæmundsson, R/V Håkon Mosby, NRV Alliance, and R/V Kristine Bonnevie. We thank the captain and crew of each of these vessels for their hard work as well as the many watch standers who have sailed on the cruises and helped collect the measurements. We also thank Frank Bahr for processing the VMADCP data collected on NRV Alliance and Magnús Danielsen for the processing of the hydrographic data collected on R/V Bjarni Sæmundsson. We acknowledge Leah Trafford McRaven for assistance with Fig. 1 and two anonymous reviewers for their helpful comments, which improved the manuscript. Funding for the project was provided by the Bergen Research Foundation Grant BFS2016REK01 (K. Våge and S. Semper), the Norwegian Research Council under Grant Agreement 231647 (K. Våge), and the U.S. National Science Foundation Grants OCE-1259618 and OCE-1756361 (R. S. Pickart and D. J. Torres), as well as OCE-1558742 (R. S. Pickart). The dataset is available on PANGAEA under https://doi.pangaea.de/10.1594/PANGAEA.903535.

Description: Author Posting. © American Meteorological Society, 2019. 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 49(6), (2019):1463-1483, doi: 10.1175/JPO-D-18-0213.1.

Description: A set of float trajectories, deployed at 1500- and 2500-m depths throughout the deep Gulf of Mexico from 2011 to 2015, are analyzed for mesoscale processes under the Loop Current (LC). In the eastern basin, December 2012–June 2014 had 〉40 floats per month, which was of sufficient density to allow capturing detailed flow patterns of deep eddies and topographic Rossby waves (TRWs), while two LC eddies formed and separated. A northward advance of the LC front compresses the lower water column and generates an anticyclone. For an extended LC, baroclinic instability eddies (of both signs) develop under the southward-propagating large-scale meanders of the upper-layer jet, resulting in a transfer of eddy kinetic energy (EKE) to the lower layer. The increase in lower-layer EKE occurs only over a few months during meander activity and LC eddy detachment events, a relatively short interval compared with the LC intrusion cycle. Deep EKE of these eddies is dispersed to the west and northwest through radiating TRWs, of which examples were found to the west of the LC. Because of this radiation of EKE, the lower layer of the eastern basin becomes relatively quiescent, particularly in the northeastern basin, when the LC is retracted and a LC eddy has departed. A mean west-to-east, anticyclone–cyclone dipole flow under a mean LC was directly comparable to similar results from a previous moored LC array and also showed connections to an anticlockwise boundary current in the southeastern basin.

Description: The authors were supported by the Department of the Interior, Bureau of Ocean Energy Management (BOEM), Contract M08PC20043 to Leidos, Inc., Raleigh, NC. The authors also wish to acknowledge the enthusiastic support of Dr. Alexis Lugo-Fernández, the BOEM Contracting Officer’s Technical Representative, during the study into the Deep Circulation of the Gulf of Mexico, using Lagrangian Methods. Thanks go to the captains and crews of the R/V Pelican and B/O Justo Sierra, J. Malbrough (LUMCON), J. Singer (Leidos), J. Valdes (WHOI), B. Guest (WHOI), and the CANEK group (CICESE).

Description: 2020-05-29

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 June 2019.

Description: Ocean surface transport is at the core of many environmental disasters, including the spread of marine plastic pollution, the Deepwater Horizon oil spill and the Fukushima nuclear contamination. Understanding and predicting flow transport, however, remains a scientific challenge, because it operates on multiple length- and time-scales that are set by the underlying dynamics. Building on the recent emergence of Lagrangian methods, this thesis investigates the present-day abilities to describe and understand the organization of flow transport at the ocean surface, including the abilities to detect the underlying key structures, the regions of stirring and regions of coherence within the flow. Over the past four years, the field of dynamical system theory has adapted several algorithms from unsupervised machine learning for the detection of Lagrangian Coherent Structures (LCS). The robustness and applicability of these tools is yet to be proven, especially for geophysical flows. An updated, parameter-free spectral clustering approach is developed and a noise-based cluster coherence metric is proposed to evaluate the resulting clusters. The method is tested against benchmarks flows of dynamical system theory: the quasi-periodic Bickley jet, the Duffng oscillator and a modified, asymmetric Duffing oscillator. The applicability of this newly developed spectral clustering method, along with several common LCS approaches, such as the Finite-Time Lyapunov Exponent, is tested in several field studies. The focus is on the ability to predict these LCS in submesoscale ocean surface flows, given all the uncertainties of the modeled and observed velocity fields, as well as the sparsity of Lagrangian data. This includes the design and execution of field experiments targeting LCS from predictive models and their subsequent Lagrangian analysis. These experiments took place in Scott Reef, an atoll system in Western Australia, and off the coast of Martha's Vineyard, Massachusetts, two case studies with tidally-driven channel flows. The FTLE and spectral clustering analyses were particularly helpful in describing key transient flow features and how they were impacted by tidal forcing and vertical velocities. This could not have been identified from the Eulerian perspective, showing the utility of the Lagrangian approach in understanding the organization of transport.

Description: Author Posting. © American Meteorological Society, 2018. 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 48 (2018): 2127-2140, doi:10.1175/JPO-D-18-0035.1.

Description: Shipboard hydrographic and velocity measurements collected in summer 2014 are used to study the evolution of the freshwater coastal current in southern Greenland as it encounters Cape Farewell. The velocity structure reveals that the coastal current maintains its identity as it flows around the cape and bifurcates such that most of the flow is diverted to the outer west Greenland shelf, while a small portion remains on the inner shelf. Taking into account this inner branch, the volume transport of the coastal current is conserved, but the freshwater transport decreases on the west side of Cape Farewell. A significant amount of freshwater appears to be transported off the shelf where the outer branch flows adjacent to the shelfbreak circulation. It is argued that the offshore transposition of the coastal current is caused by the flow following the isobaths as they bend offshore because of the widening of the shelf on the west side of Cape Farewell. An analysis of the potential vorticity shows that the subsequent seaward flux of freshwater can be enhanced by instabilities of the current. This set of circumstances provides a pathway for the freshest water originating from the Arctic, as well as runoff from the Greenland ice sheet, to be fluxed into the interior Labrador Sea where it could influence convection in the basin.

Description: Funding for this project was provided by the National Science Foundation under Grant OCE-1259618.

Description: 2019-03-11

Description: Author Posting. © American Meteorological Society, 2018. 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 48 (2018): 573-590, doi:10.1175/JPO-D-17-0206.1.

Description: Motivated by the proximity of the Northern Recirculation Gyre and the deep western boundary current in the North Atlantic, an idealized model is used to investigate how recirculation gyres and a deep flow along a topographic slope interact. In this two-layer quasigeostrophic model, an unstable jet imposed in the upper layer generates barotropic recirculation gyres. These are maintained by an eddy-mean balance of potential vorticity (PV) in steady state. The authors show that the topographic slope can constrain the northern recirculation gyre meridionally and that the gyre’s adjustment to the slope leads to increased eddy PV fluxes at the base of the slope. When a deep current is present along the topographic slope in the lower layer, these eddy PV fluxes stir the deep current and recirculation gyre waters. Increased proximity to the slope dampens the eddy growth rate within the unstable jet, altering the geometry of recirculation gyre forcing and leading to a decrease in overall eddy PV fluxes. These mechanisms may shape the circulation in the western North Atlantic, with potential feedbacks on the climate system.

Description: We gratefully acknowledge an AMS graduate fellowship (IALB) and U.S. National Science Foundation Grants OCE-1332667 and 1332834 (IALB and JMT).

Description: 2018-09-06

Description: Author Posting. © American Meteorological Society, 2018. 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 48 (2018): 739-748, doi:10.1175/JPO-D-17-0089.1.

Description: McDougall and Ferrari have estimated the global deep upward diapycnal flow in the boundary layer overlying continental slopes that must balance both downward diapycnal flow in the deep interior and the formation of bottom water around Antarctica. The decrease of perimeter of isopycnal surfaces with depth and the observed decay with height above bottom of turbulent dissipation in the deep ocean play a key role in their estimate. They argue that because the perimeter of seamounts increases with depth, the net effect of mixing around seamounts is to produce net downward diapycnal flow. While this is true along much of a seamount, it is shown here that diapycnal flow of the densest water around the seamount is upward, with buoyancy being transferred from water just above. The same is true for midocean ridges, whose perimeter is constant with depth. It is argued that mixing around seamounts and especially midocean ridges contributes positively to the global deep overturning circulation, reducing the amount of turbulence demanded over the continental slopes to balance the buoyancy budget for the bottom and deep water.

Description: This work was supported by National Science Foundation Grant OCE- 1232962.

Description: 2018-09-29

Description: Author Posting. © American Meteorological Society, 2018. 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 48 (2018): 643-646, doi:10.1175/JPO-D-17-0240.1.

Description: A simple oceanic model is presented for source–sink flow on the β plane to discuss the pathways from source to sink when transport boundary layers have large enough Reynolds numbers to be inertial in their dynamics. A representation of the flow as a Fofonoff gyre, suggested by prior work on inertial boundary layers and eddy-driven circulations in two-dimensional turbulent flows, indicates that even when the source and sink are aligned along the same western boundary the flow must intrude deep into the interior before exiting at the sink. The existence of interior pathways for the flow is thus an intrinsic property of an inertial circulation and is not dependent on particular geographical basin geometry.

Description: 2018-09-12

Description: Author Posting. © American Meteorological Society, 2018. 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 48 (2018): 1831-1848, doi:10.1175/JPO-D-18-0068.1.

Description: We present a simplified theory using reduced-gravity equations for North Atlantic Deep Water (NADW) and its variation driven by high-latitude deep-water formation. The theory approximates layer thickness on the eastern boundary with domain-averaged layer thickness and, in tandem with a mass conservation argument, retains fundamental physics for cross-equatorial flows on interannual and longer forcing time scales. Layer thickness anomalies are driven by a time-dependent northern boundary condition that imposes a southward volume flux representative of a variable source of NADW and damped by diapycnal mixing throughout the basin. Moreover, an outflowing southern boundary condition imposes a southward volume flux that generally differs from the volume flux at the northern boundary, giving rise to temporal storage of NADW within the Atlantic basin. Closed form analytic solutions for the amplitude and phase are provided when the variable source of NADW is sinusoidal. We provide a nondimensional analysis that demonstrates that solution behavior is primarily controlled by two parameters that characterize the meridional extent of the southern basin and the width of the basin relative to the equatorial deformation radius. Similar scaling applied to the time-lagged equations of Johnson and Marshall provides a clear connection to their results. Numerical simulations of reduced-gravity equations agree with analytic predictions in linear, turbulent, and diabatic regimes. The theory introduces a simple analytic framework for studying idealized buoyancy- and wind-driven cross-equatorial flows on interannual and longer time scales.

Description: This research was supported by the National Science Foundation under Grant OCE- 1634468.

Description: 2019-02-15

Description: Author Posting. © American Meteorological Society, 2018. 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 48 (2018): 883-904, doi:10.1175/JPO-D-17-0084.1.

Description: The dynamics controlling the along-valley (cross shelf) flow in idealized shallow shelf valleys with small to moderate Burger number are investigated, and analytical scales of the along-valley flows are derived. This paper follows Part I, which shows that along-shelf winds in the opposite direction to coastal-trapped wave propagation (upwelling regime) force a strong up-valley flow caused by the formation of a lee wave. In contrast, along-shelf winds in the other direction (downwelling regime) do not generate a lee wave and consequently force a relatively weak net down-valley flow. The valley flows in both regimes are cyclostrophic with 0(1) Rossby number. A major difference between the two regimes is the along-shelf length scales of the along-valley flows L. In the upwelling regime Ls, depends on the valley width W, and the wavelength lambda(1w) of the coastal-trapped lee wave arrested by the along-shelf flow U-s. In the downwelling regime L depends on the inertial length scale U-s|'f and W-c. The along-valley velocity scale in the upwelling regime, given by V-u approximate to root pi H-c/H-s integral W-c lambda(1w)/2 pi L-x (1+L-x(2)/L-c(2))(-1) e(-(pi Wc)/(lambda 1w),) is based on potential vorticity (PV) conservation and lee-wave dynamics (Hs and H, are the shelf and valley depth scales, respectively, and fis the Coriolis parameter). The velocity scale in the downwelling regime, given by |v(d)| approximate to (H-s/H-s)[1 + (L-x(2)/L-x(2))](-1) fL, is based on PV conservation. The velocity scales are validated by the numerical sensitivity simulations and can be useful for observational studies of along -valley transports. The work provides a framework for investigating cross -shelf transport induced by irregular shelf bathymetry and calls for future studies of this type under realistic environmental conditions and over a broader parameter space.

Description: Both WGZ and SJL were supported by the National Science Foundation (NSF) through Grant OCE 1154575.WGZis also supported by the NSF Grant OCE 1634965 and SJL by NSF Grant OCE 1558874.

Description: 2018-10-16

Description: Author Posting. © American Meteorological Society, 2018. 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 31 (2018): 9881-9901, doi:10.1175/JCLI-D-17-0889.1.

Description: The Atlantic meridional overturning circulation and associated poleward heat transport are balanced by northern heat loss to the atmosphere and corresponding water-mass transformation. The circulation of northward-flowing Atlantic Water at the surface and returning overflow water at depth is particularly manifested—and observed—at the Greenland–Scotland Ridge where the water masses are guided through narrow straits. There is, however, a rich variability in the exchange of water masses across the ridge on all time scales. Focusing on seasonal and interannual time scales, and particularly the gateways of the Denmark Strait and between the Faroe Islands and Shetland, we specifically assess to what extent the exchanges of water masses across the Greenland–Scotland Ridge relate to wind forcing. On seasonal time scales, the variance explained of the observed exchanges can largely be related to large-scale wind patterns, and a conceptual model shows how this wind forcing can manifest via a barotropic, cyclonic circulation. On interannual time scales, the wind stress impact is less direct as baroclinic mechanisms gain importance and observations indicate a shift in the overflows from being more barotropically to more baroclinically forced during the observation period. Overall, the observed Greenland–Scotland Ridge exchanges reflect a horizontal (cyclonic) circulation on seasonal time scales, while the interannual variability more represents an overturning circulation.

Description: This research was supported by the Research Council of Norway project NORTH (Grant 229763). Additional support for M. A. Spall was provided by National Science Foundation Grant OCE- 1558742, for T. Eldevik and S. Østerhus by the European Union’s Horizon 2020 research and innovation program project Blue-Action (Grant 727852), and for S. Østerhus by the European Framework Programs under Grant Agreement 308299 (NACLIM).

Description: Author Posting. © American Meteorological Society, 2017. 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 47 (2017): 2927-2947, doi:10.1175/JPO-D-17-0083.1.

Description: Motivated by observations in Hudson shelf valley showing stronger onshore than offshore flows, this study investigates wind-driven flows in idealized shallow shelf valleys. This first part of a two-part sequence focuses on the mechanism of the asymmetrical flow response in a valley to along-shelf winds of opposite directions. Model simulations show that (i) when the wind is in the opposite direction to coastal-trapped wave (CTW) phase propagation, the shelf flow turns onshore in the valley and generates strong up-valley transport and a standing meander on the upstream side (in the sense of CTW phase propagation) of the valley, and (ii) when the wind is in the same direction as CTW phase propagation, the flow forms a symmetric onshore detour pattern over the valley with negligible down-valley transport. Comparison of the modeled upstream meanders in the first scenario with CTW characteristics confirms that the up-valley flow results from CTWs being arrested by the wind-driven shelf flow establishing lee waves. The valley bathymetry generates an initial excessive onshore pressure gradient force that drives the up-valley flow and induces CTW lee waves that sustain the up-valley flow. When the wind-driven shelf flow aligns with CTW phase propagation, the initial disturbance generated in the valley propagates away, allowing the valley flow to adjust to roughly follow isobaths. Because of the similarity in the physical setup, this mechanism of arrested CTWs generating stronger onshore than offshore flow is expected to be applicable to the flow response in slope canyons to along-isobath background flows of opposite directions.

Description: WGZ and SJL were supported by the National Science Foundation through GrantOCE1154575.WGZ is also supported by the NSF Grant OCE 1634965 and SJL by NSF Grant OCE 1558874.

Description: 2018-06-08

Description: Author Posting. © American Meteorological Society, 2017. 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 47 (2017): 339-351, doi:10.1175/JPO-D-16-0165.1.

Description: A novel multi-iteration statistical method for studying tracer spreading using drifter data is introduced. The approach allows for the best use of the available drifter data by making use of a simple iterative procedure, which results in the statistically probable map showing the likelihood that a tracer released at some source location would visit different geographical regions, along with the associated arrival travel times. The technique is tested using real drifter data in the North Atlantic. Two examples are considered corresponding to sources in the western and eastern North Atlantic Ocean, that is, Massachusetts Bay–like and Irish Sea–like sources, respectively. In both examples, the method worked well in estimating the statistics of the tracer transport pathways and travel times throughout the entire North Atlantic. The role of eddies versus mean flow is quantified using the same technique, and eddies are shown to significantly broaden the spread of a tracer. The sensitivity of the results to the size of the source domain is investigated and causes for this sensitivity are discussed.

Description: This work was supported by the Grant OCE-1356630 from the National Science Foundation (NSF). Rypina also acknowledges NSF Grant OCE-1154641 and NASA Grant NNX14AH29G.

Description: 2017-07-31

Description: Author Posting. © American Meteorological Society, 2017. 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 47 (2017): 633-647, doi:10.1175/JPO-D-16-0089.1.

Description: Interannual variability in the volumetric water mass distribution within the North Atlantic Subtropical Gyre is described in relation to variability in the Atlantic meridional overturning circulation. The relative roles of diabatic and adiabatic processes in the volume and heat budgets of the subtropical gyre are investigated by projecting data into temperature coordinates as volumes of water using an Argo-based climatology and an ocean state estimate (ECCO version 4). This highlights that variations in the subtropical gyre volume budget are predominantly set by transport divergence in the gyre. A strong correlation between the volume anomaly due to transport divergence and the variability of both thermocline depth and Ekman pumping over the gyre suggests that wind-driven heave drives transport anomalies at the gyre boundaries. This wind-driven heaving contributes significantly to variations in the heat content of the gyre, as do anomalies in the air–sea fluxes. The analysis presented suggests that wind forcing plays an important role in driving interannual variability in the Atlantic meridional overturning circulation and that this variability can be unraveled from spatially distributed hydrographic observations using the framework presented here.

Description: DGE was supported by a Natural Environment Research Council studentship award at the University of Southampton. JMT’s contribution was supported by the U.S. National Science Foundation (Grant OCE-1332667). GF’s contribution was supported by the U.S. National Science Foundation through Grant OCE-0961713 and by the U.S. National Oceanic and Atmospheric Administration through Grant NA10OAR4310135. The contributions of JDZ and AJGN were supported by the NERC Grant ‘‘Climate scale analysis of air and water masses’’ (NE/ K012932/1). ACNG gratefully acknowledges support from the Leverhulme Trust, the Royal Society, and the Wolfson Foundation. LY was supported by NASA Ocean Vector Wind Science Team (OVWST) activities under Grant NNA10AO86G.

Description: Author Posting. © American Meteorological Society, 2017. 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 47 (2017): 1291-1305, doi:10.1175/JPO-D-16-0160.1.

Description: Along-stream variations in the dynamics of the Antarctic Circumpolar Current (ACC) impact heat and tracer transport, regulate interbasin exchange, and influence closure of the overturning circulation. Topography is primarily responsible for generating deviations from zonal-mean properties, mainly through standing meanders associated with regions of high eddy kinetic energy. Here, an idealized channel model is used to explore the spatial distribution of energy exchange and its relationship to eddy geometry, as characterized by both eddy momentum and eddy buoyancy fluxes. Variations in energy exchange properties occur not only between standing meander and quasi-zonal jet regions, but throughout the meander itself. Both barotropic and baroclinic stability properties, as well as the magnitude of energy exchange terms, undergo abrupt changes along the path of the ACC. These transitions are captured by diagnosing eddy fluxes of energy and by adopting the eddy geometry framework. The latter, typically applied to barotropic stability properties, is applied here in the depth–along-stream plane to include information about both barotropic and baroclinic stability properties of the flow. These simulations reveal that eddy momentum fluxes, and thus barotropic instability, play a leading role in the energy budget within a standing meander. This result suggests that baroclinic instability alone cannot capture the dynamics of ACC standing meanders, a challenge for models where eddy fluxes are parameterized.

Description: The authors all acknowledge support from NSF OCE-1235488. MKY also acknowledges support from the AMS Graduate Student Fellowship.

Description: 2017-10-12

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, 1982

Description: Oceanic fluctuations are dependent on geographical location. Near intense currents, the eddy field is highly energetic and has broad meridional extent. It is likely that the energy arises from instabilities of the intense current. However, the meridional extent of the linearly most unstable modes of such intense jets is much narrower than the observed region of energetic fluctuations. It is proposed here that weaker instabilities, in the linear sense, which are very weakly trapped to the current, may be the dominant waves in the far field. As a preliminary problem, the (barotropic) instability of parallel shear flow on the beta plane is discussed. An infinite zonal flow with a continuous cross-stream velocity gradient is approximated with segments of uniform flow, joined together by segments of uniform potential vorticity. This simplification allows an exact dispersion relation to be found. There are two classes of linearly unstable solutions. One type is trapped to the source of energy and has large growth rates. The second type are weaker instabilities of the shear flow which excite Rossby waves in the far field: the influence of these weaker instabilities extends far beyond that of the most unstable waves. The central focus of the thesis i: the linear stability of thin, twolayer, zonal jets on the beta plane, with both horizontal and vertical shear. The method used for the parallel shear flow is extended to the two-layer flow. Each layer of the jet has uniform velocity in the center, bordered by shear zones with zero potential vorticity gradient. The velocity in each layer outside the jet is constant in latitude. Separate linearly unstable modes arise from horizontal and vertical shear. The energy source for the vertical shear modes is nearly all potential while the source for the horizontal shear modes is both kinetic and potential. The most unstable waves are tightly trapped to the jet, within two or three deformation radii for small but nonzero beta. Rossby waves and baroclinically unstable waves (in the presence of vertical shear) exist outside the jet because of a nonzero potential vorticity gradient there. Weakly growing jet instabilities can force these waves when their phase speeds and wavelengths match. In particular, westward jets and any jets with vertical shear exterior to the jet can radiate in this sense. The radiating modes influence a large region, their decay scales inversely proportional to the growth rate. Two types of radiating instability are found: (1) a subset of the main unstable modes near marginal stability and (2) modes which appear to be destabilized neutral modes. Westward jets have more vigorously unstable radiating modes. Applications of the model are made to the eddy field south of the Gulf Stream, using data from the POLYMODE settings along 55°W and farther into the gyre at MODE. The energy decay scale and the variation of vertical structure with latitude in different frequency bands can be roughly explained by the model. The lower frequency disturbances decay more slowly and become more surface intensified in the far field. These disturbances are identified with the weak, radiating instabilities of the model. The higher frequency disturbances are more trapped and retain their vertical structure as they decay, and are identified with the trapped, strongly unstable modes of the jet.

Description: This work was supported by a grant from the National Science Foundation, Office of Atmospheric Science.

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 2017

Description: This thesis investigates the evolution of the oceanic lithosphere in a broad sense from formation to subduction, in a focused case at the ridge, and in a focused case proximal to subduction. In general, alteration of the oceanic lithosphere begins at the ridge through focused and diffuse hydrothermal flow, continues off axis through low temperature circulation, and may occur approaching subduction zones as bending related faulting provides fluid pathways. In Chapter 2 I use a dataset of thousands of microearthquakes recorded at the Rainbow massif on the Mid-Atlantic Ridge to characterize the processes which are responsible for the long-term, high-temperature, hydrothermal discharge found hosted in this oceanic core complex. I find that the detachment fault responsible for the uplift of the massif is inactive and that the axial valleys show no evidence for faulting or active magma intrusion. I conclude that the continuous, low-magnitude seismicity located in diffuse pattern in a region with seismic velocities indicating ultramafic host rock suggests that serpentinization may play a role in microearthquake generation but the seismic network was not capable of providing robust focal mechanism solutions to constrain the source characteristics. In Chapter 3 I find that the Juan de Fuca plate, which represents the young/hot end-member of oceanic plates, is lightly hydrated at upper crustal levels except in regions affected by propagator wakes where hydration of lower crust and upper mantle is evident. I conclude that at the subduction zone the plate is nearly dry at upper mantle levels with the majority of water contained in the crust. Finally, in Chapter 4 I examine samples of cretaceous age serpentinite sampled just before subduction at the Puerto Rico Trench. I show that these upper mantle rocks were completely serpentinized under static conditions at the Mid-Atlantic Ridge. Further, they subsequently underwent 100 Ma of seafloor weathering wherein the alteration products of serpentinization themselves continue to be altered. I conclude that complete hydration of the upper mantle is not the end point in the evolution of oceanic lithosphere as it spreads from the axis to subduction.

Description: Funding was provided by the National Science Foundation through grants OCE-1029305 and OCE-0961680, the Deep Ocean Exploration Institute - Ocean Ridge Initiative, and by the WHOI Academic Programs office

Description: Author Posting. © American Meteorological Society, 2017. 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 34 (2017): 1679-1691, doi:10.1175/JTECH-D-16-0162.1.

Description: For direction-finding high-frequency (HF) radar systems, the correct separation of backscattered spectral energy due to Bragg resonant waves from that due to more complex double-scattering represents a critical first step toward attaining accurate estimates of surface currents from the range-dependent radar backscatter. Existing methods to identify this “first order” region of the spectra, generally sufficient for lower-frequency radars and low-velocity or low-surface gravity wave conditions, are more likely to fail in higher-frequency systems or locations with more variable current, wave, or noise regimes, leading to elevated velocity errors. An alternative methodology is presented that uses a single and globally relevant smoothing length scale, careful pretreatment of the spectra, and marker-controlled watershed segmentation, an image processing technique, to separate areas of spectral energy due to surface currents from areas of spectral energy due to more complex scattering by the wave field or background noise present. Applied to a number of HF radar datasets with a range of operating frequencies and characteristic issues, the new methodology attains a higher percentage of successful first-order identification, particularly during complex current and wave conditions. As operational radar systems continue to expand to more systematically cover areas of high marine traffic, close approaches to ports and harbors, or offshore energy installations, use of this type of updated methodology will become increasingly important to attain accurate current estimates that serve both research and operational interests.

Description: This analysis was supported by internal funds from the Woods Hole Oceanographic Institution.

Description: 2018-02-11

Description: Author Posting. © American Meteorological Society, 2016. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 46 (2016): 3599-3621, doi:10.1175/JPO-D-16-0085.1.

Description: At continental margins, energetic deep-ocean eddies can transport shelf water offshore in filaments that wrap around the eddy. One example is that of Gulf Stream warm-core rings interacting with the Mid-Atlantic Bight shelf. The rate at which shelf water is exported in these filaments is a major unknown in regional budgets of volume, heat, and salt. This unknown transport is constrained using a series of idealized primitive equation numerical experiments wherein a surface-intensified anticyclonic eddy interacts with idealized shelf–slope topography. There is no shelfbreak front in these experiments, and shelf water is tracked using a passive tracer. When anticyclones interact with shelf–slope topography, they suffer apparent intrusions of shelf–slope water, resulting in a subsurface maximum in offshore transport. The simulations help construct an approximate model for the filament of exported water that originates inshore of any given isobath. This model is then used to derive an expression for the total volume of shelf–slope water transported by the eddy across that isobath. The transport scales with water depth, radius, and azimuthal velocity scale of the eddy. The resulting expression can be used with satellite-derived eddy properties to estimate approximate real-world transports ignoring the presence of a shelfbreak front. The expression assumes that the eddy’s edge is at the shelf break, a condition not always satisfied by real eddies.

Description: The research presented here was funded by NSF Grants OCE-1059632 and OCE-1433953. Funding support from the Academic Programs Office, and WHOI is also gratefully acknowledged.

Description: 2017-06-08

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 2017

Description: The daily heating of the ocean by the sun can create a stably stratified near-surface layer when the winds are slight and solar insolation is strong. This type of shallow stable layer is called a Diurnal Warm Layer (DWL). This thesis examines the physics and dynamics of DWLs from observations of the subtropical North Atlantic Ocean associated with the Salinity Processes in the Upper ocean Regional Study (SPURS-I). Momentum transferred from the atmosphere to the ocean through wind stress becomes trapped within the DWL, generating shear across the layer. During SPURS-I, strong diurnal shear across the DWL was coincident with enhanced turbulent kinetic energy (TKE) dissipation (𝜖, 𝜖 〉 10−5 W/kg) observed from glider microstructure profiles of the near-surface. However, a scale analysis demonstrated that surface forcing, including diurnal shear, could not be the sole mechanism for the enhanced TKE dissipation. High-frequency internal waves (𝜔 ≫ 𝑓) were observed in the upper ocean during the daytime within the DWL. Internal waves are able to transfer energy from the deep ocean into the DWL through the unstratified remnant mixed layer, which is the intervening layer between the DWL and seasonal thermocline. As the strength of the stratification of the DWL increases, so does the shear caused by the tunneling internal waves. The analysis demonstrates that internal waves can generate strong enough shear to cause a shear-induced instability, and are a plausible source of the observed enhanced TKE dissipation. Vertically-varying horizontal transport across the upper ocean occurs because a diurnal current exists within the DWL, but not in the unstratified remnant mixed layer below. Therefore, when a DWL is present, the water within DWL is horizontally transported a different distance than the water below. Coupled with nocturnal convection that mixes the DWL with the unstratified layer at night, this cycle is a mechanism for submesoscale (1-10 km) lateral diffusion across the upper ocean. Estimates of a horizontal diffusion coefficient are similar in magnitude to current estimates of submesoscale diffusion based on observations, and are likely an important source of horizontal diffusion in the upper ocean.

Description: Supported by the Department of Defense (DoD) through the National Defense Science & Engineering Graduate Fellowship (NDSEG) Program and the National Science Foundation under Grant No. OCE-1129646. The collection and analysis of data from the SPURS-I central mooring were supported under National Aeronautics and Space Administration (NASA) Grant No. NNX11AE84G and NNX14AH38G.

Description: Author Posting. © American Meteorological Society, 2015. 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 45 (2015): 2913–2932, doi:10.1175/JPO-D-14-0179.1.

Description: The oceanic deep circulation is shared between concentrated deep western boundary currents (DWBCs) and broader interior pathways, a process that is sensitive to seafloor topography. This study investigates the spreading and deepening of Denmark Strait overflow water (DSOW) in the western subpolar North Atlantic using two ° eddy-resolving Atlantic simulations, including a passive tracer injected into the DSOW. The deepest layers of DSOW transit from a narrow DWBC in the southern Irminger Sea into widespread westward flow across the central Labrador Sea, which remerges along the Labrador coast. This abyssal circulation, in contrast to the upper levels of overflow water that remain as a boundary current, blankets the deep Labrador Sea with DSOW. Farther downstream after being steered around the abrupt topography of Orphan Knoll, DSOW again leaves the boundary, forming cyclonic recirculation cells in the deep Newfoundland basin. The deep recirculation, mostly driven by the meandering pathway of the upper North Atlantic Current, leads to accumulation of tracer offshore of Orphan Knoll, precisely where a local maximum of chlorofluorocarbon (CFC) inventory is observed. At Flemish Cap, eddy fluxes carry ~20% of the tracer transport from the boundary current into the interior. Potential vorticity is conserved as the flow of DSOW broadens at the transition from steep to less steep continental rise into the Labrador Sea, while around the abrupt topography of Orphan Knoll, potential vorticity is not conserved and the DSOW deepens significantly.

Description: This work is supported by ONR Award N00014-09-1-0587, the NSF Physical Oceanography Program, and NASA Ocean Surface Topography Science Team Program.

Description: 2016-06-01

Description: Author Posting. © American Meteorological Society, 2016. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 46 (2016): 3011-3029, doi:10.1175/JPO-D-15-0248.1.

Description: Seasonal variability of the tropical Atlantic circulation is dominated by the annual cycle, but semiannual variability is also pronounced, despite weak forcing at that period. This study uses multiyear, full-depth velocity measurements from the central equatorial Atlantic to analyze the vertical structure of annual and semiannual variations of zonal velocity. A baroclinic modal decomposition finds that the annual cycle is dominated by the fourth mode and the semiannual cycle is dominated by the second mode. Similar local behavior is found in a high-resolution general circulation model. This simulation reveals that the annual and semiannual cycles of the respective dominant baroclinic modes are associated with characteristic basinwide structures. Using an idealized, linear, reduced-gravity model to simulate the dynamics of individual baroclinic modes, it is shown that the observed circulation variability can be explained by resonant equatorial basin modes. Corollary simulations of the reduced-gravity model with varying basin geometry (i.e., square basin vs realistic coastlines) or forcing (i.e., spatially uniform vs spatially variable wind) show a structural robustness of the simulated basin modes. A main focus of this study is the seasonal variability of the Equatorial Undercurrent (EUC) as identified in recent observational studies. Main characteristics of the observed EUC including seasonal variability of transport, core depth, and maximum core velocity can be explained by the linear superposition of the dominant equatorial basin modes as obtained from the reduced-gravity model.

Description: This study was supported by the Deutsche Forschungsgemeinschaft as part of the Sonderforschungsbereich 754 (SFB754) ‘‘Climate–Biogeochemistry Interactions in the Tropical Ocean’’ and through several research cruises with R/V Meteor, R/V Maria S. Merian, andR/VL’Atalante by the German Federal Ministry of Education and Research as part of the cooperative projects RACE (03F0605B) and SACUS (03G0837A) and by European Union 7th Framework Programme (FP7 2007–13) under Grant Agreement 603521 PREFACE project.

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 2016

Description: Since the Last Glacial Maximum (LGM, ~ 20,000 years ago) air temperatures warmed, sea level rose roughly 130 meters, and atmospheric concentrations of carbon dioxide increased. This thesis combines global models and paleoceanographic observations to constrain the ocean’s role in storing and transporting heat, salt, and other tracers during this time, with implications for understanding how the modern ocean works and how it might change in the future. • By combining a kinematic ocean model with “upstream” and “downstream” deglacial oxygen isotope time series from benthic and planktonic foraminifera, I show that the data are in agreement with the modern circulation, quantify their power to infer circulation changes, and propose new data locations. • An ocean general circulation model (the MITgcm) constrained to fit LGM sea surface temperature proxy observations reveals colder ocean temperatures, greater sea ice extent, and changes in ocean mixed layer depth, and suggests that some features in the data are not robust. • A sensitivity analysis in the MITgcm demonstrates that changes in winds or in ocean turbulent transport can explain the hypothesis that the boundary between deep Atlantic waters originating from Northern and Southern Hemispheres was shallower at the LGM than it is today.

Description: Support for this work came from an MIT Presidential Fellowship, an NSF Graduate Research Fellowship, and grants NASA NNX12AJ93G – Gravity data for ocean circulation and climate studies, NSF OCE-0961713 – Collaborative Research: The Physics and Statistics of Global Sea Level Change, NSF OCE-1060735 – Collaborative Research: Beyond the Instrumental Record - the Ocean Circulation at the last Glacial maximum and the deglacial sequence, and NASA NNX08AR33G – Application of Satellite Altimetry Gravity Winds and in Situ Data to Problems of the Ocean Circulation.

Description: Author Posting. © American Meteorological Society, 2016. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 46 (2016): 3263-3278, doi:10.1175/JPO-D-16-0091.1.

Description: The halocline of the Beaufort Gyre varies significantly on interannual to decadal time scales, affecting the freshwater content (FWC) of the Arctic Ocean. This study explores the role of eddies in the Ekman-driven gyre variability. Following the transformed Eulerian-mean paradigm, the authors develop a theory that links the FWC variability to the stability of the large-scale gyre, defined as the inverse of its equilibration time. The theory, verified with eddy-resolving numerical simulations, demonstrates that the gyre stability is explicitly controlled by the mesoscale eddy diffusivity. An accurate representation of the halocline dynamics requires the eddy diffusivity of 300 ± 200 m2 s−1, which is lower than what is used in most low-resolution climate models. In particular, on interannual and longer time scales the eddy fluxes and the Ekman pumping provide equally important contributions to the FWC variability. However, only large-scale Ekman pumping patterns can significantly alter the FWC, with spatially localized perturbations being an order of magnitude less efficient. Lastly, the authors introduce a novel FWC tendency diagnostic—the Gyre Index—that can be conveniently calculated using observations located only along the gyre boundaries. Its strong predictive capabilities, assessed in the eddy-resolving model forced by stochastic winds, suggest that the Gyre Index would be of use in interpreting FWC evolution in observations as well as in numerical models.

Description: GEMacknowledges the support from theHowland Postdoctoral Program Fund at WHOI and the Stanback Fellowship Fund at Caltech.MAS was supported by NSF Grants PLR-1415489 and OCE-1232389. AFT acknowledges support from NASA Award NNN12AA01C.

Description: 2017-04-20

Description: Author Posting. © American Meteorological Society, 2015. 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 28 (2015): 8574–8584, doi:10.1175/JCLI-D-14-00809.1.

Description: The subsurface ocean response to anthropogenic climate forcing remains poorly characterized. From the Coupled Model Intercomparison Project (CMIP), a robust response of the lower thermocline is identified, where the warming is considerably weaker in the subtropics than in the tropics and high latitudes. The lower thermocline change is inversely proportional to the thermocline depth in the present climatology. Ocean general circulation model (OGCM) experiments show that sea surface warming is the dominant forcing for the subtropical gyre change in contrast to natural variability for which wind dominates, and the ocean response is insensitive to the spatial pattern of surface warming. An analysis based on a ventilated thermocline model shows that the pattern of the lower thermocline change can be interpreted in terms of the dynamic response to the strengthened stratification and downward heat mixing. Consequently, the subtropical gyres become intensified at the surface but weakened in the lower thermcline, consistent with results from CMIP experiments.

Description: The work was supported by the National Basic Research Program of China (2012CB955600), the National Natural Science Foundation of China (41125019, 41206021), and the U.S. National Science Foundation (AGS 1249145, 1305719).

Description: 2016-05-01

Description: Author Posting. © American Meteorological Society, 2015. 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 45 (2015): 966–987, doi:10.1175/JPO-D-14-0110.1.

Description: A key remaining challenge in oceanography is the understanding and parameterization of small-scale mixing. Evidence suggests that topographic features play a significant role in enhancing mixing in the Southern Ocean. This study uses 914 high-resolution hydrographic profiles from novel EM-APEX profiling floats to investigate turbulent mixing north of the Kerguelen Plateau, a major topographic feature in the Southern Ocean. A shear–strain finescale parameterization is applied to estimate diapycnal diffusivity in the upper 1600 m of the ocean. The indirect estimates of mixing match direct microstructure profiler observations made simultaneously. It is found that mixing intensities have strong spatial and temporal variability, ranging from O(10−6) to O(10−3) m2 s−1. This study identifies topographic roughness, current speed, and wind speed as the main factors controlling mixing intensity. Additionally, the authors find strong regional variability in mixing dynamics and enhanced mixing in the Antarctic Circumpolar Current frontal region. This enhanced mixing is attributed to dissipating internal waves generated by the interaction of the Antarctic Circumpolar Current and the topography of the Kerguelen Plateau. Extending the mixing observations from the Kerguelen region to the entire Southern Ocean, this study infers a large water mass transformation rate of 17 Sverdrups (Sv; 1 Sv ≡ 106 m3 s−1) across the boundary of Antarctic Intermediate Water and Upper Circumpolar Deep Water in the Antarctic Circumpolar Current. This work suggests that the contribution of mixing to the Southern Ocean overturning circulation budget is particularly significant in fronts.

Description: AM was supported by the joint CSIRO–University of Tasmania Quantitative Marine Science (QMS) program and the 2009 CSIRO Wealth from Ocean Flagship Collaborative Fund. BMS was supported by the Australian Climate Change Science Program, jointly funded by the Department of the Environment and CSIRO. KLPs salary support was provided by Woods Hole Oceanographic Institution bridge support funds.

Description: 2015-10-01

Description: Author Posting. © American Meteorological Society, 2015. 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 45 (2015): 2598–2620, doi:10.1175/JPO-D-14-0249.1.

Description: Through combining analytical arguments and numerical models, this study investigates the finite-amplitude meanders of shelfbreak fronts characterized by sloping isopycnals outcropping at both the surface and the shelfbreak bottom. The objective is to provide a formula for the meander length scale that can explain observed frontal length scale variability and also be verified with observations. Considering the frontal instability to be a mixture of barotropic and baroclinic instability, the derived along-shelf meander length scale formula is [b1/(1 + a1S1/2)]NH/f, where N is the buoyancy frequency; H is the depth of the front; f is the Coriolis parameter; S is the Burger number measuring the ratio of energy conversion associated with barotropic and baroclinic instability; and a1 and b1 are empirical constants. Initial growth rate of the frontal instability is formulated as [b2(1 + a1S1/2)/(1 + a2αS1/2)]NH/L, where α is the bottom slope at the foot of the front, and a2 and b2 are empirical constants. The formulas are verified using numerical sensitivity simulations, and fitting of the simulated and formulated results gives a1 = 2.69, b1 = 14.65, a2 = 5.1 × 103, and b2 = 6.2 × 10−2. The numerical simulations also show development of fast-growing frontal symmetric instability when the minimum initial potential vorticity is negative. Although frontal symmetric instability leads to faster development of barotropic and baroclinic instability at later times, it does not significantly influence the meander length scale. The derived meander length scale provides a framework for future studies of the influences of external forces on shelfbreak frontal circulation and cross-frontal exchange.

Description: WGZ and GGG were supported by the National Science Foundation through Grant OCE-1129125.

Description: 2016-04-01

Description: Author Posting. © American Meteorological Society, 2015. 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 45 (2015):1189–1204, doi:10.1175/JPO-D-14-0122.1.

Description: Winter outcropping of the Eighteen Degree Water (EDW) and its subsequent dispersion are studied using a ° eddy-resolving simulation of the Family of Linked Atlantic Modeling Experiments (FLAME). Outcropped EDW columns in the model simulations are detected in each winter from 1990 to 1999, and particles are deployed in the center of each outcropped EDW column. Subsequently, the trajectories of these particles are calculated for the following 5 yr. The particles slowly spread away from the outcropping region into the nonoutcropping/subducted EDW region south of ~30°N and eventually to the non-EDW region in the greater subtropical gyre. Approximately 30% of the particles are found in non-EDW waters 1 yr after deployment; after 5 yr, only 25% of the particles are found within EDW. The reoutcropping time is defined as the number of years between when a particle is originally deployed in an outcropping EDW column and when that particle is next found in an outcropping EDW column. Of the particles, 66% are found to reoutcrop as EDW in 1 yr, and less than 5% of the particles outcrop in each of the subsequent 4 yr. While the individual trajectories exhibit significant eddy-like motions, the time scale of reoutcropping is primarily set by the mean circulation. The dominance of reoutcropping in 1 yr suggests that EDW outcropping contributes considerably to the persistence of surface temperature anomalies from one winter to the next, that is, the reemergence of winter sea surface temperature anomalies.

Description: We gratefully acknowledge the support from the NSF OCE Physical Oceanography program (NSF OCE-0961090 to Y-OK and J-JP; NSF OCE-0960776 to MSL and SFG; and NSF OCE-1242989 to Y-OK).

Description: 2015-10-01

Description: Author Posting. © American Meteorological Society, 2015. 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 45 (2015): 1410–1425, doi:10.1175/JPO-D-14-0192.1.

Description: The west-to-east crossover of boundary currents has been seen in mean circulation schemes from several past models of the Red Sea. This study investigates the mechanisms that produce and control the crossover in an idealized, eddy-resolving numerical model of the Red Sea. The authors also review the observational evidence and derive an analytical estimate for the crossover latitude. The surface buoyancy loss increases northward in the idealized model, and the resultant mean circulation consists of an anticyclonic gyre in the south and a cyclonic gyre in the north. In the midbasin, the northward surface flow crosses from the western boundary to the eastern boundary. Numerical experiments with different parameters indicate that the crossover latitude of the boundary currents changes with f0, β, and the meridional gradient of surface buoyancy forcing. In the analytical estimate, which is based on quasigeostrophic, β-plane dynamics, the crossover is predicted to lie at the latitude where the net potential vorticity advection (including an eddy component) is zero. Various terms in the potential vorticity budget can be estimated using a buoyancy budget, a thermal wind balance, and a parameterization of baroclinic instability.

Description: This work is supported by Award USA 00002, KSA 00011, and KSA 00011/02 made by King Abdullah University of Science and Technology (KAUST), by National Science Foundation Grants OCE0927017, OCE1154641, and OCE85464100, and by the Woods Hole Oceanographic Institution Academic Program Office.

Description: 2015-11-01

Description: Author Posting. © American Meteorological Society, 2015. 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 45 (2015): 2806–2819, doi:10.1175/JPO-D-15-0061.1.

Description: An eastward-flowing current of a homogeneous fluid with velocity U, contained in a channel of width L, impinges on an island of width of O(L), and the resulting interaction and dynamics are studied for values of the supercriticality parameter, b = βL2/U, both larger and smaller than π2. The former case is subcritical with respect to Rossby waves, and the latter is supercritical. The nature of the flow field depends strongly on b, and in particular, the nature of the flow around the island and the proportion of the flow passing to the north or south of the island are sensitive to b and to the position of the island in the channel. The problem is studied analytically in a relatively simple, nonlinear quasigeostrophic and adiabatic framework and numerically with a shallow-water model that allows a qualitative extension of the results to the equator. Although the issues involved are motivated by the interaction of the Equatorial Undercurrent and the Galapagos Islands, the analysis presented here focuses on the fundamental issue of the distinctive nature of the flow as a function of Rossby wave criticality.

Description: Supported by the National Science Foundation Grant OCE-0959381.

Description: 2016-05-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): 3033–3053, doi:10.1175/JPO-D-13-0227.1.

Description: The East Greenland Current (EGC) had long been considered the main pathway for the Denmark Strait overflow (DSO). Recent observations, however, indicate that the north Icelandic jet (NIJ), which flows westward along the north coast of Iceland, is a major separate pathway for the DSO. In this study a two-layer numerical model and complementary integral constraints are used to examine various pathways that lead to the DSO and to explore plausible mechanisms for the NIJ’s existence. In these simulations, a westward and NIJ-like current emerges as a robust feature and a main pathway for the Denmark Strait overflow. Its existence can be explained through circulation integrals around advantageous contours. One such constraint spells out the consequences of overflow water as a source of low potential vorticity. A stronger constraint can be added when the outflow occurs through two outlets: it takes the form of a circulation integral around the Iceland–Faroe Ridge. In either case, the direction of overall circulation about the contour can be deduced from the required frictional torques. Some effects of wind stress forcing are also examined. The overall positive curl of the wind forces cyclonic gyres in both layers, enhancing the East Greenland Current. The wind stress forcing weakens but does not eliminate the NIJ. It also modifies the sign of the deep circulation in various subbasins and alters the path by which overflow water is brought to the Faroe Bank Channel, all in ways that bring the idealized model more in line with observations. The sequence of numerical experiments separates the effects of wind and buoyancy forcing and shows how each is important.

Description: This study has been supported by National Science Foundation (OCE0927017 and ARC1107412).

Description: 2015-06-01

Description: Author Posting. © American Meteorological Society, 2015. 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 45 (2015): 778–791, doi:10.1175/JPO-D-14-0164.1.

Description: This study examines anisotropic transport properties of the eddying North Atlantic flow, using an idealized model of the double-gyre oceanic circulation and altimetry-derived velocities. The material transport by the time-dependent flow (quantified by the eddy diffusivity tensor) varies geographically and is anisotropic, that is, it has a well-defined direction of the maximum transport. One component of the time-dependent flow, zonally elongated large-scale transients, is particularly important for the anisotropy, as it corresponds to primarily zonal material transport and long correlation time scales. The importance of these large-scale zonal transients in the material distribution is further confirmed with simulations of idealized color dye tracers, which has implications for parameterizations of the eddy transport in non-eddy-resolving models.

Description: IK would like to acknowledge support through the NSF Grant OCE-1154923. IR was supported by the NSF OCE-1154641 and NASA Grant NNX14AH29G.

Description: 2015-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 Climate 27 (2014): 9359–9376, doi:10.1175/JCLI-D-14-00228.1.

Description: Multidecadal variability of the Atlantic meridional overturning circulation (AMOC) is examined based on a comparison of the AMOC streamfunctions in depth and in density space, in a 700-yr present-day control integration of the fully coupled Community Climate System Model, version 3. The commonly used depth-coordinate AMOC primarily exhibits the variability associated with the deep equatorward transport that follows the changes in the Labrador Sea deep water formation. On the other hand, the density-based AMOC emphasizes the variability associated with the subpolar gyre circulation in the upper ocean leading to the changes in the Labrador Sea convection. Combining the two representations indicates that the ~20-yr periodicity of the AMOC variability in the first half of the simulation is primarily due to an ocean-only mode resulting from the coupling of the deep equatorward flow and the upper ocean gyre circulation near the Gulf Stream and North Atlantic Current. In addition, the density-based AMOC reveals a gradual change in the deep ocean associated with cooling and increased density, which is likely responsible for the transition of AMOC variability from strong ~20-yr oscillations to a weaker red noise–like multidecadal variability.

Description: Support from the NOAA Climate Program Office (Grant NA10OAR4310202 and NA13OAR4310139) and NSF EaSM2 (OCE1242989) is gratefully acknowledged.

Description: 2015-06-15

Description: Author Posting. © American Meteorological Society, 2015. 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 45 (2015): 294–312, doi:10.1175/JPO-D-14-0104.1.

Description: Model analyses of an alongshelf flow over a continental shelf and slope reveal upwelling near the shelf break. A stratified, initially uniform, alongshelf flow undergoes a rapid adjustment with notable differences onshore and offshore of the shelf break. Over the shelf, a bottom boundary layer and an offshore bottom Ekman transport develop within an inertial period. Over the slope, the bottom offshore transport is reduced from the shelf’s bottom transport by two processes. First, advection of buoyancy downslope induces vertical mixing, destratifying, and thickening the bottom boundary layer. The downward-tilting isopycnals reduce the geostrophic speed near the bottom. The reduced bottom stress weakens the offshore Ekman transport, a process known as buoyancy shutdown of the Ekman transport. Second, the thickening bottom boundary layer and weakening near-bottom speeds are balanced by an upslope ageostrophic transport. The convergence in the bottom transport induces adiabatic upwelling offshore of the shelf break. For a time period after the initial adjustment, scalings are identified for the upwelling speed and the length scale over which it occurs. Numerical experiments are used to test the scalings for a range of initial speeds and stratifications. Upwelling occurs within an inertial period, reaching values of up to 10 m day−1 within 2 to 7 km offshore of the shelf break. Upwelling drives an interior secondary circulation that accelerates the alongshelf flow over the slope, forming a shelfbreak jet. The model results are compared with upwelling estimates from other models and observations near the Middle Atlantic Bight shelf break.

Description: J. Benthuysen acknowledges support from the ARC Centre of Excellence for Climate System Science (CE110001028) and the MIT/WHOI Joint Program, where this work was initiated.

Description: 2015-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 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 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): 2938–2950, doi:10.1175/JPO-D-13-0201.1.

Description: Direct observations in the Southern Ocean report enhanced internal wave activity and turbulence in a kilometer-thick layer above rough bottom topography collocated with the deep-reaching fronts of the Antarctic Circumpolar Current. Linear theory, corrected for finite-amplitude topography based on idealized, two-dimensional numerical simulations, has been recently used to estimate the global distribution of internal wave generation by oceanic currents and eddies. The global estimate shows that the topographic wave generation is a significant sink of energy for geostrophic flows and a source of energy for turbulent mixing in the deep ocean. However, comparison with recent observations from the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean shows that the linear theory predictions and idealized two-dimensional simulations grossly overestimate the observed levels of turbulent energy dissipation. This study presents two- and three-dimensional, realistic topography simulations of internal lee-wave generation from a steady flow interacting with topography with parameters typical of Drake Passage. The results demonstrate that internal wave generation at three-dimensional, finite bottom topography is reduced compared to the two-dimensional case. The reduction is primarily associated with finite-amplitude bottom topography effects that suppress vertical motions and thus reduce the amplitude of the internal waves radiated from topography. The implication of these results for the global lee-wave generation is discussed.

Description: This research was supported by the National Science Foundation under Award CMG-1024198.

Description: 2015-05-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): 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): 418–431, doi:10.1175/JPO-D-12-087.1.

Description: The overflow of the dense water mass across the Greenland–Scotland Ridge (GSR) from the Nordic Seas drives the Atlantic meridional overturning circulation (AMOC). The Nordic Seas is a large basin with an enormous reservoir capacity. The volume of the dense water above the GSR sill depth in the Nordic Seas, according to previous estimates, is sufficient to supply decades of overflow transport. This large capacity buffers overflow’s responses to atmospheric variations and prevents an abrupt shutdown of the AMOC. In this study, the authors use a numerical and an analytical model to show that the effective reservoir capacity of the Nordic Seas is actually much smaller than what was estimated previously. Basin-scale oceanic circulation is nearly geostrophic and its streamlines are basically the same as the isobaths. The vast majority of the dense water is stored inside closed geostrophic contours in the deep basin and thus is not freely available to the overflow. The positive wind stress curl in the Nordic Seas forces a convergence of the dense water toward the deep basin and makes the interior water even more removed from the overflow-feeding boundary current. Eddies generated by the baroclinic instability help transport the interior water mass to the boundary current. But in absence of a robust renewal of deep water, the boundary current weakens rapidly and the eddy-generating mechanism becomes less effective. This study indicates that the Nordic Seas has a relatively small capacity as a dense water reservoir and thus the overflow transport is sensitive to climate changes.

Description: This study has been supported by National Science Foundation (OCE0927017,ARC1107412).

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: 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: 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): 766–789, doi:10.1175/JPO-D-12-0141.1.

Description: Nonlinear energy transfers from the semidiurnal internal tide to high-mode, near-diurnal motions are documented near Kaena Ridge, Hawaii, an energetic generation site for the baroclinic tide. Data were collected aboard the Research Floating Instrument Platform (FLIP) over a 35-day period during the fall of 2002, as part of the Hawaii Ocean Mixing Experiment (HOME) Nearfield program. Energy transfer terms for a PSI resonant interaction at midlatitude are identified and compared to those for near-inertial PSI close to the M2 critical latitude. Bispectral techniques are used to demonstrate significant energy transfers in the Nearfield, between the low-mode M2 internal tide and subharmonic waves with frequencies near M2/2 and vertical wavelengths of O(120 m). A novel prefilter is used to test the PSI wavenumber resonance condition, which requires the subharmonic waves to propagate in opposite vertical directions. Depth–time maps of the interactions, formed by directly estimating the energy transfer terms, show that energy is transferred predominantly from the tide to subharmonic waves, but numerous reverse energy transfers are also found. A net forward energy transfer rate of 2 × 10−9 W kg−1 is found below 400 m. The suggestion is that the HOME observations of energy transfer from the tide to subharmonic waves represent a first step in the open-ocean energy cascade. Observed PSI transfer rates could account for a small but significant fraction of the turbulent dissipation of the tide within 60 km of Kaena Ridge. Further extrapolation suggests that integrated PSI energy transfers equatorward of the M2 critical latitude may be comparable to PSI energy transfers previously observed near 28.8°N.

Description: This work was supported by the National Science Foundation and the Office of Naval Research.

Description: 2013-10-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 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.