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(6), (2022): 1091–1110, https://doi.org/10.1175/JPO-D-21-0068.1.

Description: Hundreds of full-depth temperature and salinity profiles collected by Deepglider autonomous underwater vehicles (AUVs) in the North Atlantic reveal robust signals in eddy isopycnal vertical displacement and horizontal current throughout the entire water column. In separate glider missions southeast of Bermuda, subsurface-intensified cold, fresh coherent vortices were observed with velocities exceeding 20 cm s−1 at depths greater than 1000 m. With vertical resolution on the order of 20 m or less, these full-depth glider slant profiles newly permit estimation of scaled vertical wavenumber spectra from the barotropic through the 40th baroclinic mode. Geostrophic turbulence theory predictions of spectral slopes associated with the forward enstrophy cascade and proportional to inverse wavenumber cubed generally agree with glider-derived quasi-universal spectra of potential and kinetic energy found at a variety of locations distinguished by a wide range of mean surface eddy kinetic energy. Water-column average spectral estimates merge at high vertical mode number to established descriptions of internal wave spectra. Among glider mission sites, geographic and seasonal variability implicate bottom drag as a mechanism for dissipation, but also the need for more persistent sampling of the deep ocean.

Description: This work was funded by NSF Grant 1736217 and would not have been possible without the help of Kirk O’Donnell, James Bennett, Noel Pelland, and all contributors to Deepglider development. We additionally thank the captain crew of the R/V Atlantic Explorer and the BATS team at the Bermuda Institute of Ocean Sciences, particularly Rod Johnson, as well as Seakeepers International for their professionalism, capability, and generous assistance in deploying and recovering gliders.

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 2023.

Description: The Middle Atlantic Bight (MAB) is a highly productive ecosystem, supporting several economically important commercial fisheries. Chlorophyll enhancement at the MAB shelf-break front has been observed only intermittently, despite numerous studies that suggest persistent upwelling at the front. High resolution cross-frontal transect crossings were collected from three two-week cruises in April 2018, May 2019, and July 2019. Chapter 2 focused on applying a novel method of classifying planktonic images taken by a Video Plankton Recorder to enable processing of the large volumes of data collected with the instrument. Chapter 3 investigated cross-frontal trends by temporally averaging in both Eulerian and frontally-aligned coordinates. For April 2018, transient chlorophyll enhancement was seen at the front in individual transects and within the frontally-aligned mean transect, but not within the Eulerian mean transect. The Eulerian mean for May 2019 showed chlorophyll enhancement as a result of frontal eddies, which were further explored in chapter 4. No frontal enhancement was observed in July 2019. The frontal eddies observed in May 2019 were simulated using an idealized model, which showed that upwelling occurred within both of the frontal eddies, despite having opposite rotational directions. This result was consistent with nutrient enhancement observed within the centers of both eddies. Biological enhancement within each eddy was observed, which may have been a result of advection from source waters and/or a local response to upwelled nutrients. The influence of frontal variability and frontal eddies on nutrients and plankton at the front argues for the necessity for 3-D models to fully explain frontal behavior and its effects on biological responses.

Description: The research in this thesis was supported by the National Science Foundation (OCE-1657803) and the Dalio Foundation. Support from the WHOI Academic Programs Office is gratefully acknowledged.

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(8), (2022): 1593-1611, https://doi.org/10.1175/jpo-d-21-0180.1.

Description: This study presents novel observational estimates of turbulent dissipation and mixing in a standing meander between the Southeast Indian Ridge and the Macquarie Ridge in the Southern Ocean. By applying a finescale parameterization on the temperature, salinity, and velocity profiles collected from Electromagnetic Autonomous Profiling Explorer (EM-APEX) floats in the upper 1600 m, we estimated the intensity and spatial distribution of dissipation rate and diapycnal mixing along the float tracks and investigated the sources. The indirect estimates indicate strong spatial and temporal variability of turbulent mixing varying from O(10−6) to O(10−3) m2 s−1 in the upper 1600 m. Elevated turbulent mixing is mostly associated with the Subantarctic Front (SAF) and mesoscale eddies. In the upper 500 m, enhanced mixing is associated with downward-propagating wind-generated near-inertial waves as well as the interaction between cyclonic eddies and upward-propagating internal waves. In the study region, the local topography does not play a role in turbulent mixing in the upper part of the water column, which has similar values in profiles over rough and smooth topography. However, both remotely generated internal tides and lee waves could contribute to the upward-propagating energy. Our results point strongly to the generation of turbulent mixing through the interaction of internal waves and the intense mesoscale eddy field.

Description: The observations were funded through grants from the Australian Research Council Discovery Project (DP170102162) and Australia’s Marine National Facility. Surface drifters were provided by Dr. Shaun Dolk of the Global Drifter Program. AC was supported by an Australian Research Council Postdoctoral Fellowship. AC, HEP, and NLB acknowledge support from the Australian Government Department of the Environment and Energy National Environmental Science Program and the ARC Centre of Excellence in Climate Extremes. KP acknowledges the support from the National Science Foundation.

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, 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(8), (2022): 1677-1691, https://doi.org/10.1175/jpo-d-21-0269.1.

Description: Oceanic mesoscale motions including eddies, meanders, fronts, and filaments comprise a dominant fraction of oceanic kinetic energy and contribute to the redistribution of tracers in the ocean such as heat, salt, and nutrients. This reservoir of mesoscale energy is regulated by the conversion of potential energy and transfers of kinetic energy across spatial scales. Whether and under what circumstances mesoscale turbulence precipitates forward or inverse cascades, and the rates of these cascades, remain difficult to directly observe and quantify despite their impacts on physical and biological processes. Here we use global observations to investigate the seasonality of surface kinetic energy and upper-ocean potential energy. We apply spatial filters to along-track satellite measurements of sea surface height to diagnose surface eddy kinetic energy across 60–300-km scales. A geographic and scale-dependent seasonal cycle appears throughout much of the midlatitudes, with eddy kinetic energy at scales less than 60 km peaking 1–4 months before that at 60–300-km scales. Spatial patterns in this lag align with geographic regions where an Argo-derived estimate of the conversion of potential to kinetic energy is seasonally varying. In midlatitudes, the conversion rate peaks 0–2 months prior to kinetic energy at scales less than 60 km. The consistent geographic patterns between the seasonality of potential energy conversion and kinetic energy across spatial scale provide observational evidence for the inverse cascade and demonstrate that some component of it is seasonally modulated. Implications for mesoscale parameterizations and numerical modeling are discussed.

Description: This work was generously funded by NSF Grants OCE-1912302, OCE-1912125 (Drushka), and OCE-1912325 (Abernathey) as part of the Ocean Energy and Eddy Transport Climate Process Team.

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, 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: 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 2022.

Description: The ventilation of intermediate waters in the Labrador Sea has important implications for the strength of the Atlantic Meridional Overturning Circulation. Boundary current-interior interactions regulate the exchange of properties between the slope and the basin, which in turn regulates the magnitude of interior convection and the export of ventilated waters from the subpolar gyre. This thesis characterizes theWest Greenland Boundary Current System near Cape Farewell across a range of spatio-temporal scales. The boundary current system is composed of three velocity cores: (1) the West Greenland Coastal Current (WGCC), transporting Greenland and Arctic meltwaters on the shelf; (2) the West Greenland Current (WGC), which advects warm, saline Atlantic-origin water at depth, meltwaters at the surface, and newly-ventilated Labrador Sea Water (LSW); and (3) the Deep Western Boundary Current, which carries dense overflow waters ventilated in the Nordic Seas. The seasonal presence of the LSW and Atlantic-origin water are dictated by air-sea buoyancy forcing, while the seasonality of the WGCC is governed by remote wind forcing and the propagation of coastally trapped waves from East Greenland. Using mooring data and hydrographic surveys, we demonstrate mid-depth intensified cyclones generated at Denmark Strait are found offshore of the WGC and enhance the overflow water transport at synoptic timescales. Using mooring, hydrographic, and satellite data, we demonstrate that the WGC undergoes extensive meandering due to baroclinic instability that is enhanced in winter due to LSW formation adjacent to the current. This leads to the production of small-scale, anticyclonic eddies that can account for the entirety of wintertime heat loss within the Labrador Sea. The meanders are shown to trigger the formation of Irminger Rings downstream. Using mooring, hydrographic, atmospheric, and Lagrangian data, and a mixing model, we find that strong atmospheric storms known as forward tip jets cause upwelling at the shelfbreak that triggers offshore export of freshwater. This freshwater flux can explain the observed lack of ventilation in the eastern Labrador Sea. Together, this thesis documents previously unobserved interannual, seasonal, and synoptic-scale variability and dynamics within the West Greenland boundary current system that must be accounted for in future modeling.

Description: The work in this dissertation was funded by the National Science Foundation grants OCE-1259618 and OCE-1756361.

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): 3267–3294, https://doi.org/10.1175/JPO-D-19-0310.1.

Description: As part of the Flow Encountering Abrupt Topography (FLEAT) program, an array of pressure-sensor equipped inverted echo sounders (PIESs) was deployed north of Palau where the westward-flowing North Equatorial Current encounters the southern end of the Kyushu–Palau Ridge in the tropical North Pacific. Capitalizing on concurrent observations from satellite altimetry, FLEAT Spray gliders, and shipboard hydrography, the PIESs’ 10-month duration hourly bottom pressure p and round-trip acoustic travel time τ records are used to examine the magnitude and predictability of sea level and pycnocline depth changes and to track signal propagations through the array. Sea level and pycnocline depth are found to vary in response to a range of ocean processes, with their magnitude and predictability strongly process dependent. Signals characterized here comprise the barotropic tides, semidiurnal and diurnal internal tides, southeastward-propagating superinertial waves, westward-propagating mesoscale eddies, and a strong signature of sea level increase and pycnocline deepening associated with the region’s relaxation from El Niño to La Niña conditions. The presence of a broad band of superinertial waves just above the inertial frequency was unexpected and the FLEAT observations and output from a numerical model suggest that these waves detected near Palau are forced by remote winds east of the Philippines. The PIES-based estimates of pycnocline displacement are found to have large uncertainties relative to overall variability in pycnocline depth, as localized deep current variations arising from interactions of the large-scale currents with the abrupt topography around Palau have significant travel time variability.

Description: Support for this research was provided by Office of Naval Research Grants N00014-16-1-2668, N00014-18-1-2406, N00014-15-1-2488, and N00014-15-1-2622. R.C.M. was additionally supported by the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Weston Howland Jr. Postdoctoral Scholarship.

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): 3331–3351, https://doi.org/10.1175/JPO-D-20-0035.1.

Description: This study examines the generation of warm spiral structures (referred to as spiral streamers here) over Gulf Stream warm-core rings. Satellite sea surface temperature imagery shows spiral streamers forming after warmer water from the Gulf Stream or newly formed warm-core rings impinges onto old warm-core rings and then intrudes into the old rings. Field measurements in April 2018 capture the vertical structure of a warm spiral streamer as a shallow lens of low-density water winding over an old ring. Observations also show subduction on both sides of the spiral streamer, which carries surface waters downward. Idealized numerical model simulations initialized with observed water-mass densities reproduce spiral streamers over warm-core rings and reveal that their formation is a nonlinear submesoscale process forced by mesoscale dynamics. The negative density anomaly of the intruding water causes a density front at the interface between the intruding water and surface ring water, which, through thermal wind balance, drives a local anticyclonic flow. The pressure gradient and momentum advection of the local interfacial flow push the intruding water toward the ring center. The large-scale anticyclonic flow of the ring and the radial motion of the intruding water together form the spiral streamer. The observed subduction on both sides of the spiral streamer is part of the secondary cross-streamer circulation resulting from frontogenesis on the stretching streamer edges. The surface divergence of the secondary circulation pushes the side edges of the streamer away from each other, widens the warm spiral on the surface, and thus enhances its surface signal.

Description: Authors W. G. Zhang and D. J. McGillicuddy are both supported by the National Science Foundation through Grant OCE 1657803.

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(7), (2021): 2087–2102, https://doi.org/10.1175/JPO-D-20-0255.1.

Description: The boundary current system in the Labrador Sea plays an integral role in modulating convection in the interior basin. Four years of mooring data from the eastern Labrador Sea reveal persistent mesoscale variability in the West Greenland boundary current. Between 2014 and 2018, 197 middepth intensified cyclones were identified that passed the array near the 2000-m isobath. In this study, we quantify these features and show that they are the downstream manifestation of Denmark Strait Overflow Water (DSOW) cyclones. A composite cyclone is constructed revealing an average radius of 9 km, maximum azimuthal speed of 24 cm s−1, and a core propagation velocity of 27 cm s−1. The core propagation velocity is significantly smaller than upstream near Denmark Strait, allowing them to trap more water. The cyclones transport a 200-m-thick lens of dense water at the bottom of the water column and increase the transport of DSOW in the West Greenland boundary current by 17% relative to the background flow. Only a portion of the features generated at Denmark Strait make it to the Labrador Sea, implying that the remainder are shed into the interior Irminger Sea, are retroflected at Cape Farewell, or dissipate. A synoptic shipboard survey east of Cape Farewell, conducted in summer 2020, captured two of these features that shed further light on their structure and timing. This is the first time DSOW cyclones have been observed in the Labrador Sea—a discovery that could have important implications for interior stratification.

Description: A. P. and R. S. P. were funded by National Science Foundation Grants OCE-1259618 and OCE-1756361. I. L. B. and F. S. were funded by National Science Foundation Grants OCE-1258823 and OCE-1756272. N. P. H. was supported by the Natural Environment Research Council U.K. OSNAP program (NE/K010875/1 and NE/K010700/1). M. A. S. was supported by NSF Grants OCE-1558742 and OPP-1822334.

Description: 2021-12-08

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(8),(2020): 2315-2321, doi:10.1175/JPO-D-19-0327.1.

Description: Low-frequency currents and eddies transport sediment, pathogens, larvae, and heat along the coast and between the shoreline and deeper water. Here, low-frequency currents (between 0.1 and 4.0 mHz) observed in shallow surfzone waters for 120 days during a wide range of wave conditions are compared with theories for generation by instabilities of alongshore currents, by ocean-wave-induced sea surface modulations, and by a nonlinear transfer of energy from breaking waves to low-frequency motions via a two-dimensional inverse energy cascade. For these data, the low-frequency currents are not strongly correlated with shear of the alongshore current, with the strength of the alongshore current, or with wave-group statistics. In contrast, on many occasions, the low-frequency currents are consistent with an inverse energy cascade from breaking waves. The energy of the low-frequency surfzone currents increases with the directional spread of the wave field, consistent with vorticity injection by short-crested breaking waves, and structure functions increase with spatial lags, consistent with a cascade of energy from few-meter-scale vortices to larger-scale motions. These results include the first field evidence for the inverse energy cascade in the surfzone and suggest that breaking waves and nonlinear energy transfers should be considered when estimating nearshore transport processes across and along the coast.

Description: Funding was provided by a Vannevar Bush Faculty Fellowship [from OUSD(R&E)] and NSF.

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Cusack, J. M., Brearley, J. A., Garabato, A. C. N., Smeed, D. A., Polzin, K. L., Velzeboer, N., & Shakespeare, C. J. Observed eddy-internal wave interactions in the Southern Ocean. Journal of Physical Oceanography, 50(10), (2020): 3042-3062, doi:10.1175/JPO-D-20-0001.1.

Description: The physical mechanisms that remove energy from the Southern Ocean’s vigorous mesoscale eddy field are not well understood. One proposed mechanism is direct energy transfer to the internal wave field in the ocean interior, via eddy-induced straining and shearing of preexisting internal waves. The magnitude, vertical structure, and temporal variability of the rate of energy transfer between eddies and internal waves is quantified from a 14-month deployment of a mooring cluster in the Scotia Sea. Velocity and buoyancy observations are decomposed into wave and eddy components, and the energy transfer is estimated using the Reynolds-averaged energy equation. We find that eddies gain energy from the internal wave field at a rate of −2.2 ± 0.6 mW m−2, integrated from the bottom to 566 m below the surface. This result can be decomposed into a positive (eddy to wave) component, equal to 0.2 ± 0.1 mW m−2, driven by horizontal straining of internal waves, and a negative (wave to eddy) component, equal to −2.5 ± 0.6 mW m−2, driven by vertical shearing of the wave spectrum. Temporal variability of the transfer rate is much greater than the mean value. Close to topography, large energy transfers are associated with low-frequency buoyancy fluxes, the underpinning physics of which do not conform to linear wave dynamics and are thereby in need of further research. Our work suggests that eddy–internal wave interactions may play a significant role in the energy balance of the Southern Ocean mesoscale eddy and internal wave fields.

Description: Funding for DIMES was provided by U.K. Natural Environment Research Council (NERC) Grants NE/E007058/1 and NE/E005667/1. JMC acknowledges the support of a NERC PhD studentship, and ACNG that of the Royal Society and the Wolfson Foundation. NV acknowledges support from the ARC Centre of Excellence for Climate Extremes (CLEX) Honours Scholarship and the ANU PBSA Partnership - Spotless Scholarship. CJS acknowledges support from an ARC Discovery Early Career Researcher Award DE180100087 and an Australian National University Futures Scheme award. Numerical simulations were conducted on the National Computational Infrastructure (NCI) facility, Canberra, Australia. This study has been conducted using E.U. Copernicus Marine Service Information. We thank two anonymous reviewers for their comments which helped to improve the manuscript significantly. Codes and output files are available online at the project repository (https://github.com/jessecusack/DIMES_eddy_wave_interactions).

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(7), (2020): 1839-1852, https://doi.org/10.1175/JPO-D-19-0209.1.

Description: The Lagrangian characteristics of the surface flow field arising when an idealized, anticyclonic, mesoscale, isolated deep-ocean eddy collides with continental slope and shelf topography are explored. In addition to fluid parcel trajectories, we consider the trajectories of biological organisms that are able to navigate and swim, and for which shallow water is a destination. Of particular interest is the movement of organisms initially located in the offshore eddy, the manner in which the eddy influences the ability of the organisms to reach the shelf break, and the spatial and temporal distributions of organisms that do so. For nonswimmers or very slow swimmers, the organisms arrive at the shelf break in distinct pulses, with different pulses occurring at different locations along the shelf break. This phenomenon is closely related to the episodic formation of trailing vortices that are formed after the eddy collides with the continental slope, turns, and travels parallel to the coast. Analysis based on finite-time Lyapunov exponents reveals initial locations of all successful trajectories reaching the shoreline, and provides maps of the transport pathways showing that much of the cross-shelf-break transport occurs in the lee of the eddy as it moves parallel to the shore. The same analysis shows that the onshore transport is interrupted after a trailing vortex detaches. As the swimming speeds are increased, the organisms are influenced less by the eddy and tend to show up en mass and in a single pulse.

Description: IR and LP were supported by National Science Foundation (NSF) Grant OCE-1558806. DC was supported by NSF U.S. National Science Foundation’s Physical Oceanography program through Grants OCE-1059632 and OCE-1433953 as well as the Academic Programs Office, Woods Hole Oceanographic Institution. We acknowledge high-performance computing support from Yellowstone (http://n2t.net/ark:/85065/d7wd3xhc) provided by NCAR’s Computational and Information Systems Laboratory, sponsored by the National Science Foundation.

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(8),(2021): 2425–2441, https://doi.org/10.1175/JPO-D-20-0317.1.

Description: The frequency and latitudinal dependence of the midlatitude wind-driven meridional overturning circulation (MOC) is studied using theory and linear and nonlinear applications of a quasigeostrophic numerical model. Wind forcing is varied either by changing the strength of the wind or by shifting the meridional location of the wind stress curl pattern. At forcing periods of less than the first-mode baroclinic Rossby wave basin crossing time scale, the linear response in the middepth and deep ocean is in phase and opposite to the Ekman transport. For forcing periods that are close to the Rossby wave basin crossing time scale, the upper and deep MOC are enhanced, and the middepth MOC becomes phase shifted, relative to the Ekman transport. At longer forcing periods the deep MOC weakens and the middepth MOC increases, but eventually for long enough forcing periods (decadal) the entire wind-driven MOC spins down. Nonlinearities and mesoscale eddies are found to be important in two ways. First, baroclinic instability causes the middepth MOC to weaken, lose correlation with the Ekman transport, and lose correlation with the MOC in the opposite gyre. Second, eddy thickness fluxes extend the MOC beyond the latitudes of direct wind forcing. These results are consistent with several recent studies describing the four-dimensional structure of the MOC in the North Atlantic Ocean.

Description: This study was supported by National Science Foundation Grant OCE-1947290.

Description: 2022-01-13

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(6), (2020): 1557-1582, doi:10.1175/JPO-D-19-0239.1.

Description: We examine various contributions to the vertical velocity field within large mesoscale eddies by analyzing multiple solutions to an idealized numerical model of a representative anticyclonic warm core Gulf Stream ring. Initial conditions are constructed to reproduce the observed density and nutrient profiles collected during the Warm Core Rings Program. The contributions to vertical fluxes diagnosed from the numerical simulations are compared against a divergence-based, semidiagnostic equation and a generalized omega equation to better understand the dynamics of the vertical velocity field. Frictional decay alone is found to be ineffective in raising isopycnals and transporting nutrients to the upper ocean. With representative wind forcing, the magnitude of vorticity gradient–induced Ekman pumping is not necessarily larger than the current-induced counterpart on a time scale relevant to ecosystem response. Under realistic forcing conditions, strain deformation can perturb the ring to be noncircular and induce vertical velocities much larger than the Ekman vertical velocities. Nutrient budget diagnosis, together with analysis of the relative magnitudes of the various types of vertical fluxes, allows us to describe the time-scale dependence of nutrient delivery. At time scales that are relevant to individual phytoplankton (from hours to days), the magnitudes of nutrient flux by Ekman velocities and deformation-induced velocities are comparable. Over the life span of a typical warm core ring, which can span multiple seasons, surface current–induced Ekman pumping is the most effective mechanism in upper-ocean nutrient enrichment because of its persistence in the center of anticyclones regardless of the direction of the wind forcing.

Description: This work was supported by the National Science Foundation Ocean Science Division under Grant OCE-1558960. PG also acknowledges support of the NASA Physical Oceanography Program Grant NNX16H59G. KC would like to thank D. McGillicuddy Jr. for inspiring discussions and suggestions during the course of this study. Constructive comments from two anonymous reviewers are appreciated.

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, 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(5),(2020): 1227-1244, doi:10.1175/JPO-D-19-0280.1.

Description: The Nordic seas are commonly described as a single basin to investigate their dynamics and sensitivity to environmental changes when using a theoretical framework. Here, we introduce a conceptual model for a two-basin marginal sea that better represents the Nordic seas geometry. In our conceptual model, the marginal sea is characterized by both a cyclonic boundary current and a front current as a result of different hydrographic properties east and west of the midocean ridge. The theory is compared to idealized model simulations and shows good agreement over a wide range of parameter settings, indicating that the physics in the two-basin marginal sea is well captured by the conceptual model. The balances between the atmospheric buoyancy forcing and the lateral eddy heat fluxes from the boundary current and the front current differ between the Lofoten and the Greenland Basins, since the Lofoten Basin is more strongly eddy dominated. Results show that this asymmetric sensitivity leads to opposing responses depending on the strength of the atmospheric buoyancy forcing. Additionally, the front current plays an essential role for the heat and volume budget of the two basins, by providing an additional pathway for heat toward the interior of both basins via lateral eddy heat fluxes. The variability of the temperature difference between east and west influences the strength of the different flow branches through the marginal sea and provides a dynamical explanation for the observed correlation between the front current and the slope current of the Norwegian Atlantic Current in the Nordic seas.

Description: We thank Ilker Fer and two anonymous reviewers whose comments improved this paper. S. L. Ypma and S. Georgiou were supported by NWO (Netherlands Organisation for Scientific Research) VIDI Grant 864.13.011 awarded to C. A. Katsman. M. A. Spall was supported by National Science Foundation Grants OCE-1558742 and OPP-1822334. E. Lambert is funded by the ERA4CS project INSeaPTION. The model data analyzed in this study are available on request from the corresponding author. This study has been conducted using E.U. Copernicus Marine Service Information. The altimeter products were produced by Ssalto/Duacs and distributed by Aviso+, with support from CNES (https://www.aviso.altimetry.fr).

Description: 2020-10-27

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): 2867-2881, doi: 10.1175/JPO-D-19-0072.1.

Description: The Antarctic Circumpolar Current plays a central role in the ventilation of heat and carbon in the global ocean. In particular, the isopycnal slopes determine where each water mass outcrops and thus how the ocean interacts with the atmosphere. The region-integrated isopycnal slopes have been suggested to be eddy saturated, that is, stay relatively constant as the wind forcing changes, but whether or not the flow is saturated in realistic present day and future parameter regimes is unknown. This study analyzes an idealized two-layer quasigeostrophic channel model forced by a wind stress and a residual overturning generated by a mass flux across the interface between the two layers, with and without a blocking ridge. The sign and strength of the residual overturning set which way the isopycnal slopes change with the wind forcing, leading to an increase in slope with an increase in wind forcing for a positive overturning and a decrease in slope for a negative overturning, following the usual conventions; this behavior is caused by the dominant standing meander weakening as the wind stress weakens causing the isopycnal slopes to become more sensitive to changes in the wind stress and converge with the slopes of a flat-bottomed simulation. Eddy saturation only appears once the wind forcing passes a critical level. These results show that theories for saturation must have both topography and residual overturning in order to be complete and provide a framework for understanding how the isopycnal slopes in the Southern Ocean may change in response to future changes in wind forcing.

Description: MKY and RF acknowledge support through NSF Awards OCE-1536515 and AGS-1835576. MKY acknowledges funding from NDSEG. GRF was supported by NSF OCE-1459702. We are very grateful for conversations with David Marshall, Andrew Stewart, and two anonymous reviewers that greatly improved the manuscript. The code for running the model is found at https://github.com/mkyoungs/JPO-QG-Channel.

Description: 2020-04-30

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(1), (2020): 255-268, doi:10.1175/JPO-D-19-0166.1.

Description: Regional connectivity is important to the global climate salinity response, particularly because salinity anomalies do not have a damping feedback with atmospheric freshwater fluxes and may therefore be advected over long distances by ocean circulation, resulting in nonlocal influences. Climate model intercomparison experiments such as CMIP5 exhibit large uncertainty in some aspects of the salinity response, hypothesized here to be a result of ocean dynamics. We use two types of Lagrangian particle tracking experiments to investigate pathways of exchange for salinity anomalies. The first uses forward trajectories to estimate average transport time scales between water cycle regimes. The second uses reverse trajectories and a freshwater accumulation method to quantitatively identify remote influences in the salinity response. Additionally, we compare velocity fields with both resolved and parameterized eddies to understand the impact of eddy stirring on intergyre exchange. These experiments show that surface anomalies are readily exchanged within the ocean gyres by the mean circulation, but intergyre exchange is slower and largely eddy driven. These dynamics are used to analyze the North Atlantic salinity response to climate warming and water cycle intensification, where the system is broadly forced with fresh surface anomalies in the subpolar gyre and salty surface anomalies in the subtropical gyres. Under these competing forcings, strong intergyre eddy fluxes carry anomalously salty subtropical water into the subpolar gyre which balances out much of the local freshwater input.

Description: We acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups (listed in Table 1 of this paper) for producing and making available their model output. We also thank the creators of the SODA and ECCO reanalysis products. This work was supported by NASA Headquarters under the NASA Earth and Space Science Fellowship Program Award 80NSSC17K0372, and by National Science Foundation Award OCE-1433132. The SODA outputs used here can be accessed at http://www.atmos.umd.edu/~ocean/, and the ECCO outputs at https://ecco.jpl.nasa.gov/. Data from the CMIP5 ensemble is available at https://esgf-node.llnl.gov/projects/esgf-llnl/. The particle tracking code used for these experiments can be found at https://github.com/slevang/particle-tracking.

Description: 2020-07-20

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 (2), (2019): 607-630, doi:10.1175/JPO-D-18-0166.1.

Description: The Lagrangian motion in the eddy field produced from an unstable retrograde jet along the shelf break is studied from idealized numerical experiments with a primitive equation model. The jet is initially in thermal wind balance with a cross-isobath density gradient and is not subjected to any atmospheric forcing. Over the course of the model integration, the jet becomes unstable and produces a quasi-stationary eddy field over a 2-month period. During this period, the cross-slope flow at the shelf break is characterized by along-slope correlation scales of O(10) km and temporal correlation scales of a few days. The relative dispersion of parcels across isobaths is found to increase with time as tb, where 1 〈 b 〈 2. This mixed diffusive–ballistic regime appears to reflect the combined effects of (i) the short length scales of velocity correlation at the shelf break and (ii) the seaward excursion of monopolar and dipolar vortices. Cross-slope dispersion is greater offshore of the front than inshore of the front, as offshore parcels are both subducted onshore below density surfaces and translated offshore with eddies. Nonetheless, the exchange of parcels across the jet remains very limited on the monthly time scale. Particles originating from the bottom experience upward displacements of a few tens of meters and seaward displacements of O(100) km, suggesting that the eddy activity engendered by an unstable along-slope jet provides another mechanism for bottom boundary layer detachment near the shelf edge.

Description: The author expresses his gratitude to the researchers who contributed to the development and public dissemination of POM [for a list of contributors, see Mellor (2002) and comments in the source code]. Discussions with Kenneth Brink, Hyodae Seo, and Weifeng Zhang have been helpful. Comments provided by Kenneth Brink on a draft are gratefully acknowledged. The criticism from two anonymous reviewers allowed us to better focus the manuscript and to significantly improve its clarity. This work has been supported by Grant OCE-1556400 from the U.S. National Science Foundation.

Description: 2020-02-18

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(9), (2019): 2337-2343, doi:10.1175/JPO-D-19-0097.1.

Description: The weakly unstable, two-layer model of baroclinic instability is studied in a configuration in which the flow is perturbed at the inflow section of a channel by a slow and periodic perturbation. In a parameter regime where the governing equation would be the Lorenz equations for chaos if the development occurs only in time, the solution behavior becomes considerably more complex as a function of time and downstream coordinate. In the absence of the beta effect it has earlier been shown that the chaotic behavior along characteristics renders the solution nearly discontinuous in the slow downstream coordinate of the asymptotic model. The additional presence of the beta effect, although expunging the chaos for large enough values of the beta parameter, also provides an additional mechanism for abrupt spatial change.

Description: 2020-02-28

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Oceanography and Applied Ocean Science and Engineering at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution June 2019.

Description: Submesoscale currents, with horizontal length scales of 1-20 km, are an important element of upper ocean dynamics. These currents play a crucial role in the horizontal and vertical redistribution of tracers, the cascade of tracer variance to smaller scales, and in linking the mesoscale circulation with the dissipative scales. This thesis investigates submesoscale flows and their properties using Lagrangian trajectories of observed and modeled drifters. We analyze statistics of observed drifter pairs to characterize turbulent dispersion at submeso-scales. Contrary to theoretical expectations, we find that nonlocal velocity gradients associated with mesoscale eddies dominate the separation of drifters even at the kilometer scale. At submeso-scales, we observe energetic motions, such as near-inertial oscillations, that contribute to the energy spectrum but are inefficient at dispersion. Using trajectories in a model of submesoscale turbulence, we find that, if drifters have a vertical separation, vertical shear dominates the dispersion and conceals horizontal dispersion regimes from drifter observations. Particularly in submesoscale flows, vertical shear is orders of magnitude larger than horizontal gradients in velocity. Since conventional drifters in the ocean are not affected by vertical shear, it is likely that drifter-derived diffusivity underestimates the diffusivity that a tracer would experience. Lastly, we test and apply cluster-based methods, using three or more drifters, to estimate the velocity gradient tensor. Since velocity gradients become large at submesoscales, the divergence, strain, and vorticity control the evolution and deformation of clusters of drifters. Observing the velocity gradients using drifters, enables us to further constrain the governing dynamics and decipher submesoscale motions from inertia-gravity waves. These insights provide a Lagrangian perspective on submesoscale flows that illuminates scales that are challenging to observe from other platforms. We reveal observational and theoretical challenges that need to be overcome in future investigations.

Description: My doctoral studies in the WHOI/MIT Joint Program were funded by the National Science Foundation (OCE-I434788) and the Office of Naval Research (N00014-13-1-0451, Grant N00014-16-1-2470).

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: 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): 773-794, doi:10.1175/JPO-D-17-0205.1.

Description: Fourteen autonomous profiling floats, equipped with CTDs, were deployed in the deep eastern and western basins of the Gulf of Mexico over a four-year interval (July 2011–August 2015), producing a total of 706 casts. This is the first time since the early 1970s that there has been a comprehensive survey of water masses in the deep basins of the Gulf, with better vertical resolution than available from older ship-based surveys. Seven floats had 14-day cycles with parking depths of 1500 m, and the other half from the U.S. Argo program had varying cycle times. Maps of characteristic water masses, including Subtropical Underwater, Antarctic Intermediate Water (AAIW), and North Atlantic Deep Water, showed gradients from east to west, consistent with their sources being within the Loop Current (LC) and the Yucatan Channel waters. Altimeter SSH was used to characterize profiles being in LC or LC eddy water or in cold eddies. The two-layer nature of the deep Gulf shows isotherms being deeper in the warm anticyclonic LC and LC eddies and shallower in the cold cyclones. Mixed layer depths have an average seasonal signal that shows maximum depths (~60 m) in January and a minimum in June–July (~20 m). Basin-mean steric heights from 0–50-m dynamic heights and altimeter SSH show a seasonal range of ~12 cm, with significant interannual variability. The translation of LC eddies across the western basin produces a region of low homogeneous potential vorticity centered over the deepest part of the western basin.

Description: The authors were supported by the Department of the Interior, Bureau of Ocean Energy Management (BOEM), Contract M08PC20043 to Leidos, Inc., Raleigh, North Carolina.

Description: 2018-10-04

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 2018

Description: Highly migratory marine fishes support valuable commercial fisheries worldwide. Yet, many target species have proven difficult to study due to long-distance migrations and regular deep diving. Despite the dominance of oceanographic features, such as fronts and eddies, in the open ocean, the biophysical interactions occurring at the oceanic (sub)mesoscale (〈 100 km) remain poorly understood. This leads to a paucity of knowledge on oceanographic associations of pelagic fishes and hinders management efforts. With ever-improving oceanographic datasets and modeling outputs, we can leverage these tools both to derive better estimates of animal movements and to quantify fish-environment interactions. In this thesis, I developed analytical tools to characterize the biophysical interactions influencing animal behavior and species’ ecology in the open ocean. A novel, observation-based likelihood framework was combined with a Bayesian state-space model to improve geolocation estimates for archival-tagged fishes using oceanographic profile data. Using this approach, I constructed track estimates for a large basking shark tag dataset using a high-resolution oceanographic model and discovered a wide range of movement strategies. I also applied this modeling approach to track archival-tagged swordfish, which revealed affinity for thermal front and eddy habitats throughout the North Atlantic that was further corroborated by synthesizing these results with a fisheries-dependent conventional tag dataset. An additive modeling approach applied to longline catch-per-unit effort data further highlighted the biophysical interactions that characterize variability in swordfish catch. In the final chapter, I designed a synergistic analysis of high-resolution, 3D shark movements and satellite observations to quantify the influence of mesoscale oceanography on blue shark movements and behavior. This work demonstrated the importance of eddies in structuring the pelagic ocean by influencing the movements of an apex predator and governing the connectivity between deep scattering layer communities and deep-diving, epipelagic predators. Together, these studies demonstrate the breadth and depth of information that can be garnered through the integration of traditional animal tagging and oceanographic research with cutting-edge analytical approaches and high-resolution oceanographic model and remote sensing datasets, the product of which provides a transformative view of the biophysical interactions occurring in and governing the structure of the pelagic ocean.

Description: Supported by the NASA Earth and Space Science Fellowship, the MIT John S. Hennessy Fellowship, the MIT Martin Family Society of Fellows for Sustainability Fellowship, the WHOI Ocean Venture, Grassle, and James Stratton Fellowships and the WHOI Academic Programs Office. This research and its dissemination was supported by funds from National Geographic, Amazon Web Services, the Explorers Club, Rolex, Sigma Xi, the MIT Center for International Studies, WHOI Access to the Sea and Coastal Ocean Institute Funds, MIT Graduate Student Council, MIT Student Assistance Fund, WHOI Biology Department, American Fisheries Society, WHOI Academic Programs Office

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): 2703-2719, doi:10.1175/JPO-D-17-0245.1.

Description: A new set of deep float trajectory data collected in the Gulf of Mexico from 2011 to 2015 at 1500- and 2500-m depths is analyzed to describe mesoscale processes, with particular attention paid to the western Gulf. Wavelet analysis is used to identify coherent eddies in the float trajectories, leading to a census of the basinwide coherent eddy population and statistics of the eddies’ kinematic properties. The eddy census reveals a new formation region for anticyclones off the Campeche Escarpment, located northwest of the Yucatan Peninsula. These eddies appear to form locally, with no apparent direct connection to the upper layer. Once formed, the eddies drift westward along the northern edge of the Sigsbee Abyssal Gyre, located in the southwestern Gulf of Mexico over the abyssal plain. The formation mechanism and upstream sources for the Campeche Escarpment eddies are explored: the observational data suggest that eddy formation may be linked to the collision of a Loop Current eddy with the western boundary of the Gulf. Specifically, the disintegration of a deep dipole traveling under the Loop Current eddy Kraken, caused by the interaction with the northwestern continental slope, may lead to the acceleration of the abyssal gyre and the boundary current in the Bay of Campeche region.

Description: The authors were supported by the Department of the Interior, Bureau of Ocean Energy Management (BOEM), Contract M10PC00112 to Leidos, Inc., Raleigh, North Carolina.

Description: 2019-05-07

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): 879-894, doi:10.1175/JPO-D-16-0196.1.

Description: Models show that surface cooling over a sloping continental shelf should give rise to baroclinic instability and thus tend toward gravitationally stable density stratification. Less is known about how alongshore winds affect this process, so the role of surface momentum input is treated here by means of a sequence of idealized, primitive equation numerical model calculations. The effects of cooling rate, wind amplitude and direction, bottom slope, bottom friction, and rotation rate are all considered. All model runs lead to instability and an eddy field. While instability is not strongly affected by upwelling-favorable alongshore winds, wind-driven downwelling substantially reduces eddy kinetic energy, largely because the downwelling circulation plays a similar role to baroclinic instability by flattening isotherms and so reducing available potential energy. Not surprisingly, cross-shelf winds appear to have little effect. Analysis of the model runs leads to quantitative relations for the wind effect on eddy kinetic energy for the equilibrium density stratification (which increases as the cooling rate increases) and for eddy length scale.

Description: This research was supported by the National Science Foundation Physical Oceanography Program through Grant OCE-1433953.

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): 855-866, doi:10.1175/JPO-D-16-0194.1.

Description: Mesoscale eddies shape the Beaufort Gyre response to Ekman pumping, but their transient dynamics are poorly understood. Climate models commonly use the Gent–McWilliams (GM) parameterization, taking the eddy streamfunction to be proportional to an isopycnal slope s and an eddy diffusivity K. This local-in-time parameterization leads to exponential equilibration of currents. Here, an idealized, eddy-resolving Beaufort Gyre model is used to demonstrate that carries a finite memory of past ocean states, violating a key GM assumption. As a consequence, an equilibrating gyre follows a spiral sink trajectory implying the existence of a damped mode of variability—the eddy memory (EM) mode. The EM mode manifests during the spinup as a 15% overshoot in isopycnal slope (2000 km3 freshwater content overshoot) and cannot be explained by the GM parameterization. An improved parameterization is developed, such that is proportional to an effective isopycnal slope , carrying a finite memory γ of past slopes. Introducing eddy memory explains the model results and brings to light an oscillation with a period ≈ 50 yr, where the eddy diffusion time scale TE ~ 10 yr and γ ≈ 6 yr are diagnosed from the eddy-resolving model. The EM mode increases the Ekman-driven gyre variance by γ/TE ≈ 50% ± 15%, a fraction that stays relatively constant despite both time scales decreasing with increased mean forcing. This study suggests that the EM mode is a general property of rotating turbulent flows and highlights the need for better observational constraints on transient eddy field characteristics.

Description: GEM acknowledges the Stanback Postdoctoral Fellowship Fund at Caltech and the Howland Postdoctoral Program Fund at WHOI. MAS was supported by NSF Grants PLR-1415489 and OCE- 1232389. AFT acknowledges support from NSF OCE- 1235488.

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: 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: 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): 2645-2662, doi:10.1175/JPO-D-15-0191.1.

Description: The occurrence, drivers, and implications of small-scale O(2–5) km diameter coherent vortices, referred to as submesoscale eddies, over the inner shelf south of Martha’s Vineyard, Massachusetts, are examined using high-frequency (HF), radar-based, high-resolution (400 m) observations of surface currents. Within the 300 km2 study area, eddies occurred at rates of 1 and 4 day−1 in winter and summer, respectively. Most were less than 5 h in duration, smaller than 4 km in diameter, and rotated less than once over their lifespan; 60% of the eddies formed along the eastern edge of study area, adjacent to Wasque Shoal, and moved westward into the interior, often with relative vorticity greater than f. Eddy generation was linked to vortex stretching on the ebb and flood tide as well as the interaction of the spatially variable tide and the wind-driven currents; however, these features had complex patterns of surface divergence and stretching. Eddies located away from Wasque Shoal were related to the movement of wind-driven surface currents, as wind direction controlled where eddies formed as well as density effects. Using an analysis of particles advected within the radar-based surface currents, the observed eddies were found to be generally leaky, losing 60%–80% of particles over their lifespan, but still more retentive than the background flow. As a result, the combined translation and rotational effects of the observed eddies were an important source of lateral exchange for surface waters over the inner shelf.

Description: The HF radar data utilized here were obtained using internal funding from the Woods Hole Oceanographic Institution. The analysis was supported by NSF OCE Grant 1332646.

Description: 2017-02-19

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: Eddies in the ocean move westwards. Those shed by western boundary currents must then interact with continental shelf-slope topography at the western boundary. The presence of other eddies and mean lows complicates this simple picture, yet satellite images show that mesoscale eddies translating near the shelfbreak routinely affect the continental shelves of the Mid-Atlantic Bight, the Gulf of Mexico etc. The consequent cross-shelfbreak transports are currently of unknown importance to shelf budgets of heat, salt and volume. Thus motivated, this thesis uses idealized continuously stratified numerical experiments to explore eddy-slope interactions under four questions: 1. Can the continental slope prevent an eddy from reaching the shelfbreak? 2. What is the structure of the eddy-driven offshore low? 3. How is the continental shelf affected by an eddy at the shelfbreak? 4. Given surface observations, can one estimate the volume of water transported across the shelfbreak? The experiments show that the efficiency of Rossby wave radiation from the eddy controls whether it can cross isobaths: by radiating energy the eddy becomes shallow enough to move into shallower water. For wide continental slopes, relative to an eddy diameter, a slope can prevent an anticyclone from reaching the shelfbreak by shutting down such radiation. For narrow continental slopes, the interaction repeatedly produces dipoles, whose cyclonic halves contain shelf-slope water stacked over eddy water. The formation of such cyclones is explained. Then, the structure of shelf lows forced by the eddy are studied: their vertical structures are rationalized and scalings derived for their cross-isobath scales; for example, the extent to which the eddy influences the shelf. A recipe for estimating cross-isobath transports based on eddy surface properties is put forward. Finally, the findings are tested against observations in the Middle Atlantic Bight of the northeastern United States.

Description: I acknowledge high-performance computing support from Yellowstone (ark:/85065/d7wd3xhc) provided by NCAR’s Computational and Information Systems Laboratory, sponsored by the National Science Foundation. The research presented here was funded by National Science Foundation grants OCE-1059632 and OCE-1433953. Funding support from the Academic Programs Oice, WHOI is also gratefully acknowledged.

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, 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): 1717-1734, doi:10.1175/JPO-D-15-0124.1.

Description: The contribution of warm-core anticyclones shed by the Irminger Current off West Greenland, known as Irminger rings, to the restratification of the upper layers of the Labrador Sea is investigated in the 1/12° Family of Linked Atlantic Models Experiment (FLAME) model. The model output, covering the 1990–2004 period, shows strong similarities to observations of the Irminger Current as well as ring observations at a mooring located offshore of the eddy formation region in 2007–09. An analysis of fluxes in the model shows that while the majority of heat exchange with the interior indeed occurs at the site of the Irminger Current instability, the contribution of the coherent Irminger rings is modest (18%). Heat is provided to the convective region mainly through noncoherent anomalies and enhanced local mixing by the rings facilitating further exchange between the boundary and interior. The time variability of the eddy kinetic energy and the boundary to interior heat flux in the model are strongly correlated to the density gradient between the dense convective region and the more buoyant boundary current. In FLAME, the density variations of the boundary current are larger than those of the convective region, thereby largely controlling changes in lateral fluxes. Synchronous long-term trends in temperature in the boundary and the interior over the 15-yr simulation suggest that the heat flux relative to the temperature of the interior is largely steady on these time scales.

Description: The authors were supported in this work by the U.S. National Science Foundation.

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): 606–612, doi:10.1175/JPO-D-14-0221.1.

Description: Mesoscale intrathermocline lenses are observed throughout the World Ocean and are commonly attributed to water mass anomalies advected from a distant origin. An alternative mechanism of local generation is offered herein, in which eddy–wind interaction can create lens-shaped disturbances in the thermocline. Numerical simulations illustrate how eddy–wind-driven upwelling in anticyclones can yield a convex lens reminiscent of a mode water eddy, whereas eddy–wind-driven downwelling in cyclones produces a concave lens that thins the mode water layer (a cyclonic “thinny”). Such transformations should be observable with long-term time series in the interiors of mesoscale eddies.

Description: Support of this research by the National Science Foundation and National Aeronautics and Space Administration is gratefully acknowledged.

Description: 2015-08-01

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

Description: The spatial and temporal evolution of stratified shear instabilities is quantified in a highly stratified and energetic estuary. The measurements are made using high-resolution acoustic backscatter from an array composed of six calibrated broadband transducers connected to a six-channel high-frequency (120-600 kHz) broadband acoustic backscatter system. The array was mounted on the bottom of the estuary and looking upward. The spatial and temporal evolution of the waves is described in terms of their wavelength, amplitude and turbulent dissipation as a function of space and time. The observed waves reach an arrested growth stage nearly 10 times faster than laboratory and numerical experiments performed at much lower Reynolds number. High turbulent dissipation rates are observed within the braid regions of the waves, consistent with the rapid transition to arrested growth. Further, it appears that the waves do not undergo periodic doubling and do not collapse once their maximum amplitude is reached. Under some conditions long internal waves may provide the perturbation that decreases the gradient Richardson number so as to initiate shear instability. The initial Richardson number for the observed instabilities is likely between 0.1 and 0.2 based on the slope and growth rate of the shear instabilities.

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): 2820–2835, doi:10.1175/JPO-D-15-0101.1.

Description: The response of a convective ocean basin to variations in atmospheric temperature is explored using numerical models and theory. The results indicate that the general behavior depends strongly on the frequency at which the atmosphere changes relative to the local response time to air–sea heat flux. For high-frequency forcing, the convective region in the basin interior is essentially one-dimensional and responds to the integrated local surface heat flux anomalies. For low-frequency forcing, eddy fluxes from the boundary current into the basin interior become important and act to suppress variability forced by the atmosphere. A theory is developed to quantify this time-dependent response and its influence on various oceanic quantities. The amplitude and phase of the temperature and salinity of the convective water mass, the meridional overturning circulation, the meridional heat flux, and the air–sea heat flux predicted by the theory compare well with that diagnosed from a series of numerical model calculations in both strongly eddying and weakly eddying regimes. Linearized analytic solutions provide direct estimates of each of these quantities and demonstrate their dependence on the nondimensional numbers that characterize the domain and atmospheric forcing. These results highlight the importance of mesoscale eddies in modulating the mean and time-dependent ocean response to atmospheric variability and provide a dynamical framework with which to connect ocean observations with changes in the atmosphere and surface heat flux.

Description: This study was supported by the National Science Foundation under Grant OCE-1232389.

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): 1822–1842, doi:10.1175/JPO-D-14-0147.1.

Description: Influences of time-dependent precipitation on water mass transformation and heat budgets in an idealized marginal sea are examined using theoretical and numerical models. The equations proposed by Spall in 2012 are extended to cases with time-dependent precipitation whose form is either a step function or a sinusoidal function. The theory predicts the differences in temperature and salinity between the convective water and the boundary current as well as the magnitudes of heat fluxes into the marginal sea and across the sea surface. Moreover, the theory reveals that there are three inherent time scales: relaxation time scales for temperature and salinity and a precipitation time scale. The relaxation time scales are determined by a steady solution of the theoretical model with steady precipitation. The relaxation time scale for temperature is always smaller than that for salinity as a result of not only the difference in the form of fluxes at the surface but also the variation in the eddy transport from the boundary current. These three time scales and the precipitation amplitude determine the strength of the ocean response to changes in precipitation and the phase relation between precipitation, changes in salinity and temperature, and changes in heat fluxes. It is demonstrated that the theoretical predictions agree qualitatively well with results from the eddy-resolving numerical model. This demonstrates the fundamental role of mesoscale eddies in the ocean response to time-dependent forcing and provides a framework with which to assess the extent to which observed variability in marginal sea convection and water mass transformation are consistent with an external forcing by variations in precipitation.

Description: This work was initiated at the 2013 WHOI Geophysical Fluid Dynamics Summer Program, which was supported by the National Science Foundation and the Office of Naval Research. This work was also supported by Grant-in-Aid for Research Fellow (25·8466) of the Ministry of Education, Culture, Sports and Technology (MEXT), Japan, the Program for Leading Graduate Schools, MEXT, Japan (YY), and by the National Science Foundation Grant OCE-1232389 (MAS).

Description: 2016-01-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): 1735–1756, doi:10.1175/JPO-D-14-0238.1.

Description: The Lofoten basin of the Nordic Seas is recognized as a crucial component of the meridional overturning circulation in the North Atlantic because of the large horizontal extent of Atlantic Water and winter surface buoyancy loss. In this study, hydrographic and current measurements collected from a mooring deployed in the Lofoten basin from July 2010 to September 2012 are used to describe water mass transformation and the mesoscale eddy field. Winter mixed layer depths (MLDs) are observed to reach approximately 400 m, with larger MLDs and denser properties resulting from the colder 2010 winter. A heat budget of the upper water column requires lateral input, which balances the net annual heat loss of ~80 W m−2. The lateral flux is a result of mesoscale eddies, which dominate the velocity variability. Eddy velocities are enhanced in the upper 1000 m, with a barotropic component that reaches the bottom. Detailed examination of two eddies, from April and August 2012, highlights the variability of the eddy field and eddy properties. Temperature and salinity properties of the April eddy suggest that it originated from the slope current but was ventilated by surface fluxes. The properties within the eddy were similar to those of the mode water, indicating that convection within the eddies may make an important contribution to water mass transformation. A rough estimate of eddy flux per unit boundary current length suggests that fluxes in the Lofoten basin are larger than in the Labrador Sea because of the enhanced boundary current–interior density difference.

Description: The work was supported by NSF OCE 0850416.

Description: 2015-12-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 February 2014

Description: Motivated by inspection of complex underwater environments, we have developed a system for multi-sensor SLAM utilizing both structured and unstructured environmental features. We present a system for deriving planar constraints from sonar data, and jointly optimizing the vehicle and plane positions as nodes in a factor graph. We also present a system for outlier rejection and smoothing of 3D sonar data, and for generating loop closure constraints based on the alignment of smoothed submaps. Our factor graph SLAM backend combines loop closure constraints from sonar data with detections of visual fiducial markers from camera imagery, and produces an online estimate of the full vehicle trajectory and landmark positions. We evaluate our technique on an inspection of a decomissioned aircraft carrier, as well as synthetic data and controlled indoor experiments, demonstrating improved trajectory estimates and reduced reprojection error in the final 3D map.

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, 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): 104–132, doi:10.1175/JPO-D-14-0032.1.

Description: Three mechanisms for self-induced Ekman pumping in the interiors of mesoscale ocean eddies are investigated. The first arises from the surface stress that occurs because of differences between surface wind and ocean velocities, resulting in Ekman upwelling and downwelling in the cores of anticyclones and cyclones, respectively. The second mechanism arises from the interaction of the surface stress with the surface current vorticity gradient, resulting in dipoles of Ekman upwelling and downwelling. The third mechanism arises from eddy-induced spatial variability of sea surface temperature (SST), which generates a curl of the stress and therefore Ekman pumping in regions of crosswind SST gradients. The spatial structures and relative magnitudes of the three contributions to eddy-induced Ekman pumping are investigated by collocating satellite-based measurements of SST, geostrophic velocity, and surface winds to the interiors of eddies identified from their sea surface height signatures. On average, eddy-induced Ekman pumping velocities approach O(10) cm day−1. SST-induced Ekman pumping is usually secondary to the two current-induced mechanisms for Ekman pumping. Notable exceptions are the midlatitude extensions of western boundary currents and the Antarctic Circumpolar Current, where SST gradients are strong and all three mechanisms for eddy-induced Ekman pumping are comparable in magnitude. Because the polarity of current-induced curl of the surface stress opposes that of the eddy, the associated Ekman pumping attenuates the eddies. The decay time scale of this attenuation is proportional to the vertical scale of the eddy and inversely proportional to the wind speed. For typical values of these parameters, the decay time scale is about 1.3 yr.

Description: This work was funded by NASA Grants NNX08AI80G, NNX08AR37G, NNX13AD78G, NNX10AE91G, NNX13AE47G, and NNX10AO98G.

Description: 2015-07-01

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

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

Description: Published

Description: 1618/1637

Description: 4A. Clima e Oceani

Description: JCR Journal

Description: restricted

Description: Author Posting. © American Meteorological Society, 2013. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 26 (2013): 9839–9859, doi:10.1175/JCLI-D-12-00647.1.

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

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

Description: 2014-06-15

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

Description: Underwater acoustic communication is an extremely complex field that faces many challenges due to the time-varying nature of the ocean environment. Vector sensors are a proven technology that when utilizing their directional sensing capabilities allows us to minimize the effect of interfering noise sources. A traditional pressure sensor array has been the standard for years but suffers at degraded signal to noise ratios (SNR) and requires maneuvers or a lengthly array aperture to direction find. This thesis explores the effect of utilizing a vector sensor array to steer to the direction of signal arrival and the effect it has on equalization of the signal at degraded SNRs. It was demonstrated that utilizing a single vector sensor we were able steer to the direction of arrival and improve the ability of an equalizer to determine the transmitted signal. This improvement was most prominent when the SNR was degraded to levels of 0 and 10 dB where the performance of the vector sensor outperformed that of the pressure sensor in nearly 100% of cases. Finally, this performance improvement occured with a savings in computational expense.

Description: ONR Grants #N00014-11-10426 and #N00014-07-10738

Description: Author Posting. © American Meteorological Society, 2014. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 44 (2014): 229–245, doi:10.1175/JPO-D-12-0218.1.

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

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

Description: 2014-07-01

Description: Author Posting. © American Meteorological Society, 2014. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 44 (2014): 2593–2616, doi:10.1175/JPO-D-13-0120.1.

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

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

Description: 2015-04-01

Description: Author Posting. © American Meteorological Society, 2014. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 27 (2014): 2842–2860, doi:10.1175/JCLI-D-13-00227.1.

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

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

Description: 2014-10-15

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

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

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

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

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

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

Description: Author Posting. © American Meteorological Society, 2013. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 43 (2013): 283–300, doi:10.1175/JPO-D-11-0240.1.

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

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

Description: 2013-08-01

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

Description: There is a growing consensus that the sound generated by breaking waves is responsible for much of the ambient noise level in the ocean. While numerous field measurements have shown a strong correlation between the ambient noise spectrum level (N) in the range 100Hz to 25kHz and wind speed in the ocean, very little has been done to establish a comparable correlation between the ambient noise spectrum level and surface wave field parameters. The difficulty in establishing this relationship is remarkable given that the frequency and intensity of wave breaking are dependent on the characteristics of the wave field. In Fall 1991, an experiment was conducted from the research platform Flip 130 kilometers off the coast of Oregon, where the ambient noise between 2.5 and 25 kHz, the wind speed, and the sea surface elevation using wire wave gauges were measured. The correlation between N and the root mean square wave amplitude a was found to be poor but could be improved if the swell was filtered out from the wave elevation time series. The influence of swell on the value of a was disproportionate to the level of ambient noise since its characteristics were not directly due to the local wind-wave conditions. Observations of the dependence of the high frequency wind waves and the directional wave spectrum under turning winds suggested that the high frequency wave components responded more quickly to changes in the wind speed and wind direction than the energy-containing frequencies. The ambient noise level also correlated well with the root mean square wave slopes. This is consistent with previous laboratory measurements which showed that the steepness of a packet of waves correlates with the strength of wave breaking and with characteristics of breaking waves such as loss of momentum flux, dissipation, initial volume of air entrained, mixing, and sound generation. Comparisons of surface wave dissipation estimates using field measurements and models developed by Phillips (1985) and Hasselmann (1974) show that although the two models have very different forms, they give values that are comparable in magnitude. The relationship between the ambient noise level and log of dissipation give correlation coefficients (0.93-0.95) that are comparable to those between ambient noise and wind speed. The mean square acoustic pressure was shown to vary with the dissipation, with p2 ∝ D0.6-0.8. The results suggest that measurements of ambient sound may prove to be useful in inferring surface wave dissipation.

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 1997

Description: Acoustic propagation in the ocean inevitably encounters inhomogeneities of various types, which give rise to scattering. Acoustic scattering from rough water /bottom interfaces comprised of exposed rocks and sea mountains gives way to volumetric scattering in areas with fiat interfaces and thick sediment cover. The data analysis of the ARSRP backscattering experiment revealed that random inhomogeneities in two irregular layers beneath the seafloor were the primary contributors to oblique backscattering in a sediment pond on the western flank of the Mid-Atlantic Ridge. In this thesis, an attempt has been made to model monostatic backscattering from 3-D volume inhomogeneities in the sediment and to compare the results with the ARSRP backscattering data. A scattering process cannot be modeled correctly without a proper account of the incident field . Several approximate propagation models have been evaluated against the exact solution, while the appropriateness of using the equivalent surface scattering strength in volume scattering characterizations is studied. This study concludes that precautions need to be taken in modeling both the propagation effects and the scattering mechanisms associated with the bottom volume scattering process. A volume scattering model based on perturbation theory and the Born approximation is developed incorporating contributions from both sound speed and density fluctuations. With the propagation part handled accurately by OASES and random fluctuations generated effectively by a new scheme modified from the spectral method, the model is capable of simulating the monostatic backscattered field and time series due to 3-D volumetric sediment inhomogeneities. Both the characteristic length scale and power spectrum descriptions of the random inhomogeneities are shown to have great impact on the backscattered field by parameter studies in a free-space scenario. The important roles played by horizontal anisotropy and the vertical correlation of the random field have been demonstrated. Density fluctuations are further confirmed to be the dominant force in backscattering. The model matches the ARSRP backscattering data very well, with the fluctuations of sound speed and density in the two irregular layers described by a power law type of power spectrum.

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 1997

Description: This work is concerned with coherent communication by means of acoustic signals over underwater communication channels. The estimated scattering functions of real data ranging from the Arctic environment to tropical waters show that underwater communication channels can not be captured by a single, simple channel model. This thesis considers mainly a subset of underwater communication channels where the Doppler spread is more severe than the delay spread. An appropriate representation of the linear time- variant channel is introduced, and the wide sense stationary uncorrelated scattering (WSSUS) channel assumption enables characterization in terms of scattering functions. The concept of Doppler lines, which are frequency domain filters, is used in the derivation of a receiver for Doppler spread channels. The channel is simulated by means of a ray representation for the acoustic field and a time-variant FIR filter. The impact of physical ocean processes on the Doppler spread is demonstrated, and from this modeling explanations for the Doppler spread observed in real data are obtained. A decision feedback equalizer (DFE) adapted with recursive least squares (RLS) is analyzed, and its limit with respect to pure Doppler spread is found. By using the DFE with a phase locked loop (PLL) suboptimal system behavior is found, and this is verified on real data. In the case of a simple Doppler shift the cross-ambiguity function is used to estimate the shift, and the received signal is phase rotated to compensate this before it enters the receiver. A modified RLS called the time updated RLS (TU-RLS) is presented, and it is used in a new receiver. This receiver is initialized by means of the cross-ambiguity function and the performance is characterized by probability of decoding error vs delay spread, Doppler spread and SNR. The receiver uses Doppler lines to compensate both discrete and continuous Doppler spread. The receiver stability depends on the conditioning of the block diagonal correlation matrix propagated by the TU-RLS. The receiver is used to decode both real and simulated data, and some of these data are severely Doppler spread.

Description: The work in this thesis was funded in part by the Norwegian Research Counsil under a general PhD contract, the education office of Woods Hole Oceanographic Institution and the Office of Naval Research under contracts N00014-95-0153 and N00014-95-l-0322.

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 1997

Description: As part of the Shallow Water Acoustics in a Random Medium (SWARM) experiment [1], a sixteen element WHOI vertical line array (WVLA) was moored in 70 meters of water off the New Jersey coast. This array was sampled at 1395 Hz or higher for the seven days it was deployed. Tomography sources with carrier frequencies of 224 and 400 Hz were moored about 32 km shoreward, such that the acoustic path was anti-parallel to the primary propagation direction for shelf generated internal wave solitons. Two models for the propagation of normal modes through a 2-D waveguide with solitary internal wave (soliton) scattering included are developed to help in understanding the very complicated mode arrivals seen at the WVLA. The simplest model uses the Preisig and Duda [2] sharp interface approximation for solitons, allowing for rapid analysis of the effects of various numbers of solitons on mode arrival statistics. The second model, using SWARM thermistor string data to simulate the actual SWARM waveguides, is more realistic, but much slower. The analysis of the actual WVLA data yields spread, bias, wander, and intensity fluctuation signals that are modulated at tidal frequencies. The signals are consistent with predicted relationships to the internal wave distributions in the waveguides.

Description: The funds for my education were provided by the Office of Naval Research through an ONR Fellowship (MIT award 002734-001); the funds for SWARM were also provided by the Office of Naval Research through ONR Grant N00014-95-0051.

Description: Author Posting. © American Meteorological Society, 2012. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 42 (2012): 2206–2228, doi:10.1175/JPO-D-11-0191.1.

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

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

Description: 2013-06-01

Description: Author Posting. © American Meteorological Society, 2013. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 43 (2013): 905–919, doi:10.1175/JPO-D-12-0150.1.

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

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

Description: 2013-11-01

Description: Author Posting. © American Meteorological Society, 2013. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 43 (2013): 1611–1626, doi:10.1175/JPO-D-12-0204.1.

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

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

Description: 2014-02-01

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

Description: Adaptive equalization is an important aspect of communication systems in various environments. It is particularly important in underwater acoustic communication systems, as the channel has a long delay spread and is subject to the effects of time- varying multipath fading and Doppler spreading. The design of the adaptation algorithm has a profound influence on the performance of the system. In this thesis, we explore this aspect of the system. The emphasis of the work presented is on applying concepts from inference and decision theory and information theory to provide an approach to deriving and analyzing adaptation algorithms. Limited work has been done so far on rigorously devising adaptation algorithms to suit a particular situation, and the aim of this thesis is to concretize such efforts and possibly to provide a mathematical basis for expanding it to other applications. We derive an algorithm for the adaptation of the coefficients of an equalizer when the receiver has limited or no information about the transmitted symbols, which we term the Soft-Decision Directed Recursive Least Squares algorithm. We will demonstrate connections between the Expectation-Maximization (EM) algorithm and the Recursive Least Squares algorithm, and show how to derive a computationally efficient, purely recursive algorithm from the optimal EM algorithm. Then, we use our understanding of Markov processes to analyze the performance of the RLS algorithm in hard-decision directed mode, as well as of the Soft-Decision Directed RLS algorithm. We demonstrate scenarios in which the adaptation procedures fail catastrophically, and discuss why this happens. The lessons from the analysis guide us on the choice of models for the adaptation procedure. We then demonstrate how to use the algorithm derived in a practical system for underwater communication using turbo equalization. As the algorithm naturally incorporates soft information into the adaptation process, it becomes easy to fit it into a turbo equalization framework. We thus provide an instance of how to use the information of a turbo equalizer in an adaptation procedure, which has not been very well explored in the past. Experimental data is used to prove the value of the algorithm in a practical context.

Description: Support from the agencies that funded this research- the Academic Programs Office at WHOI and the Office of Naval Research (through ONR Grant #N00014-07-10738 and #N00014-10-10259).

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 1998

Description: Work on the forward problem in zooplankton bioacoustics has resulted in the identification of three categories of acoustic scatterers: elastic-shelled (e.g. pteropods), fluid-like (e.g. euphausiids), and gas-bearing (e.g. siphonophores). The relationship between backscattered energy and animal biomass has been shown to vary by a factor of —19,000 across these categories, so that to make accurate estimates of zooplankton biomass from acoustic backscatter measurements of the ocean, the acoustic characteristics of the species of interest must be well-understood. This thesis describes the development of both feature based and model based classification techniques to invert broadband acoustic echoes from individual zooplankton for scatterer type, as well as for particular parameters such as animal orientation. The feature based Empirical Orthogonal Function Classifier (EOFC) discriminates scatterer types by identifying characteristic modes of variability in the echo spectra, exploiting only the inherent characteristic structure of the acoustic signatures. The model based Model Parameterisation Classifier (MPC) classifies based on correlation of observed echo spectra with simplified parameterisations of theoretical scattering models for the three classes. The Covariance Mean Variance Classifiers (CMVC) are a set of advanced model based techniques which exploit the full complexity of the theoretical models by searching the entire physical model parameter space without employing simplifying parameterisations. Three different CMVC algorithms were developed: the Integrated Score Classifier (ISC), the Pairwise Score Classifier (PSC) and the Bayesian Probability Classifier (BPC); these classifiers assign observations to a class based on similarities in covariance, mean, and variance, while accounting for model ambiguity and validity. These feature based and model based inversion techniques were successfully applied to several thousand echoes acquired from broadband (-350 kHz - 750 kHz) insonifications of live zooplankton collected on Georges Bank and the Gulf of Maine to determine scatterer class. CMVC techniques were also applied to echoes from fluid-like zooplankton (Antarctic krill) to invert for angle of orientation using generic and animal-specific theoretical and empirical models. Application of these inversion techniques in situ will allow correct apportionment of backscattered energy to animal biomass, significantly improving estimates of zooplankton biomass based on acoustic surveys.

Description: Thanks to the WHOI/MIT Joint Program Education Office for partial funding. Other sources of funding for my thesis work include the Ocean Acoustics, Oceanic Biology and URIP programs of the Office of Naval Research grant numbers N00014-89- J-1729, N00014-95-1-0287 and N00014-924-1527, and the Biological Oceanography program of the National Science Foundation grant number OCE-9201264.

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 1990

Description: Travel-times of acoustic signals were measured between a bottom-mounted source near Oahu and four bottom-mounted receivers located near Washington, Oregon, and California in 1988 and 1989. This paper discusses the observed tidal signals. At three out of four receivers, observed travel times at M2 and S2 periods agree with predictions from barotropic tide models to within ±30° in phase and a factor of 1.6 in amplitude. The discrepancy at the fourth receiver can be removed by including predicted effects of phase-locked baroclinic tides generated by seamounts. Our estimates of barotropic M2 tidal dissipation by seamounts vary between 2 x 1016 and 1 X 1018 erg·s-1. The variation by two orders of magnitude is due to uncertainties in the numbers and sizes of seamounts. The larger dissipation (1 x 1018 erg·s-1) is the same order as previous estimates and amounts to 4% of the total dissipation at M2.

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

Description: In 1991 the Heard Island Feasibility Test demonstrated that it is possible to transmit coded acoustic signals nearly half way around the world. One of the key issues in the feasibility test was to determine the spatial structure of the received transmissions. In this thesis, data from the Canadian Defense Research Establishment Pacific horizontal line array is used to form an estimate of the directional power spectrum. This spectrum determines if any horizontal multipath is detectable. The preliminary signal conditioning, including frequency spectrum estimation and demodulation required before beamforming is described. Conventional and adaptive beamforming methods are examined with synthetic data to demonstrate the limitations on the directional spectrum results. The principle result of this work is that no stable horizontal multipath is evident. The mean arrival angle for the five hours of data analyzed is 212° ± 1.5°.

Description: The Office of Naval Research provided funding for the author under the ONR fellowship program.

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: 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 problem of a low-frequency acoustic plane wave incident upon a free surface coupled to a semi-infinite elastic plate surface, is solved using an analytic approach based on the Wiener-Hopf method. By low-frequency it is meant that the elastic properties of the plate are adequately described by the thin plate equation (kH ≲ 1). The diffraction problem relates to issues in long range sound propagation through partially ice-covered Arctic waters, where open leads or polynya on the surface represent features from which acoustic energy can be diffracted or scattered. This work focusses on ice as the material for the elastic plate surface, and, though the solution methods presented here have applicability to general edge diffraction problems, the results and conclusions are directed toward the ice lead diffraction process. The work begins with the derivation of an exact solution to a canonical problem: a plane wave incident upon a free surface (Dirichlet boundary condition) coupled to a perfectly rigid surface (Neumann boundary condition). Important features of the general edge diffraction problem are included here, with the solution serving as a guideline to the more complicated solutions presented later involving material properties of the boundary. The ice material properties are first addressed using the locally reacting approximation for the input impedance of an ice plate, wherein the effects of elasticity are ignored. This is followed by use of the thin plate equation to describe the input impedance, which incorporates elements of elastic wave propagation. An important issue in working with the thin plate equation is the fluid loading pertaining to sea ice and low-frequency acoustics, which cannot be characterized by simplifying heavy or light fluid loading limits. An approximation to the exact kernel of the Wiener-Hopf functional equation is used here, which is valid in this mid-range fluid loading regime. Use of this approximate kernel allows one to proceed to a complete and readily interpretable solution for the far field diffracted pressure, which includes a subsonic flexural wave in the ice plate. By using Green's theorem, in conjunction with the behavior of the diffracted field along the two-part planar boundary, the functional dependence of ∏D (total diffracted power) in terms of k (wavenumber), H (ice thickness), α (grazing angle) and the combined elastic properties of the ice sheet and ambient medium, is determined. A means to convert ∏D into an estimate of dB loss per bounce is developed using ray theoretical methods, in order to demonstrate a mechanism for acoustic propagation loss attributed directly to ice lead diffraction effects. Data from the 1984 MIZEX (Marginal Ice Zone Experiments) narrow-band acoustic transmission experiments are presented and discussed in this context.

Description: I also gratefully acknowledge financial support provided by the WHOI Education Office and the Office of Naval Research.

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: Variability of the Florida Current has been monitored via acoustic tomography. A reciprocal tomography experiment was conducted in the eastern half of the Florida Straits during mid October and November, 1983. A triangular array of transceivers, with leg separations of approximately 45 kilometers, was deployed at 27°N. The presence of a surface mixed layer in the region allowed for the ducted propagation of acoustic energy in the surface layer. A deeper layer was sampled by an unresolved group of refracted, bottom reflected ray arrivals. Incorporating the complete set of arrivals, we are able to obtain depth dependent estimates of the temperature field, current velocity, and relative vorticity. The oceanography of the region has been shown to be dominated by the lateral shifting of the surface and subsurface core of the Florida Current. The influx of westward flowing water through the Northwest Providence Channel at 26°N also appears as a large scale signal in the eastern Florida Straits. Low frequency fluctuations of temperature, current velocity, and vorticity occur at periods ranging from several days to nearly two weeks, and are intimately related to meandering of the Florida Current system.

Description: This research was carried out under ONR contract N00014-86-K-0751.

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 September 1990

Description: The propagation of low frequency seismo-acoustic waves in the Arctic Ocean ice canopy is examined through the analysis of hydrophone and geophone data sets collected in 1987 at an ice camp designated PRUDEX in the Beaufort Sea. Study of the geophone time series generated by under-ice explosive detonations reveals not only the expected longitudinal and flexural waves in the ice plate, but also an unexpected horizontally-polarized transverse (SH) wave arriving at a higher amplitude than the other wave types. The travel paths of all three observed wave types are found to be refracted in the horizontal plane along a line coincident with a known ridge separating the ice canopy locally into two distinct half-plates, the first of thin first year ice and the second of thicker multi-year ice. The origin of the SH wave appears to be near the detonation and not associated with the interaction of longitudinal, flexural or waterborne waves with the ridge line. The need to determine the exact location of each detonation from the received time series highlights the dramatic superiority of geophones over hydrophones in this application, as does the ability to detect the anomalous SH waves and the refracted ray paths, neither of which are visible in the hydrophone data. Inversion of the geophone data sets for the low frequency elastic parameters of the ice is conducted initially by treating the ice as a single homogeneous isotropic plate to demonstrate the power of SAFARI numerical modeling in this application. A modified stationary phase approach is then used to extend SAFARI modeling to invert the data sets for the elastic parameters of the two ice half-plates simultaneously. The compressional/shear bulk wave speeds estimated in the half-plates, 3500/1750 m/s in the multi-year ice and 3000/1590 m/ s in the new ice, are comparable to previously obtained values; however, the compressional/shear attenuation values in the two half-plates, 1.0/2.99 d/Bλ. and 1.0/2.67 dB/λ, respectively, are somewhat greater than previously measured values and four times greater than estimates extrapolated from high frequency data.

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 September 1990

Description: Pulse-like acoustic signals are transmitted from an acoustic source near Oahu to seven receivers off the west coast of the United States for a 124-day period in 1988. Acoustic travel-time oscillations are observed in the received signal at periods between 15 and 23 hours, which are caused by barotropic (or first or second mode baroclinic) flu ctuations in the ocean. It is shown that these fluctuations cannot be local processes isolated to either the source or to the receivers. It is further shown that resonant barotropic gravity wave modes (Platzman et al., 1981) are not consistent with the data. The cause of these flu ctuations remains unresolved, but the data and other oceanographic measurements put many constraints on the process causing these fluctuations.

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 1996

Description: The shallow water acoustic channel supports far-field propagation in a discrete set of modes. Ocean experiments have confirmed the modal nature of acoustic propagation, but no experiment has successfully excited only one of the suite of mid-frequency propagating modes propagating in a coastal environment. The ability to excite a single mode would be a powerful tool for investigating shallow water ocean processes. A feedback control algorithm incorporating elements of adaptive estimation, underwater acoustics, array processing and control theory to generate a high-fidelity single mode is presented. This approach also yields a cohesive framework for evaluating the feasibility of generating a single mode with given array geometries, noise characteristics and source power limitations. Simulations and laboratory waveguide experiments indicate the proposed algorithm holds promise for ocean experiments.

Description: Josko Catipovic funded my research for summer of 1992 on the Office of Naval Research Grant Number N00014-92-J-1661 and from June 1993 through August 1995 on Defense Advanced Research Projects Agency Grant Number MDA972-92-J- 1041. The Office of Naval Research Grant N00014-95-1-0362 to MIT supported the computer facilities used to do much of this work.

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

Description: The goal of this thesis is to develop a methodology to interpret sound scattered from the seafloor in terms of seafloor structure and subseafloor geological properties. Specifically, this work has been directed towards the interpretation of matched-filtered, beamformed monostatic acoustic reverberation data acquired on the west flank of the Mid-Atlantic Ridge when the seafloor is insonified by a band-limited, low-grazing-angle acoustic pulse. This research is based on the hypothesis that observed backscatter signals are produced by a combination of seafloor (interface) scattering and subseafloor (volume) scattering from structure having variations at scale lengths similar to the wavelength of the insonifying acoustic field. Analysis of monostatic reverberation data acquired during the Site A experiment (Run 1) of the Acoustic Reverberation Special Research Program 1993 Acoustics Cruise suggests that the scattered signals cannot be accounted for quantitatively in terms of large-scale slope, even though a strong correspondence between high intensity backscatter and seafloor ridges is observed. In order to investigate and quantify the actual sources of seafloor scattering, a numerical modeling study of seafloor models is undertaken using a finite-difference solution to the elastic wave equation. Geological data available at Site A and published reports describing geological properties of similar deep ocean crustal regions are used to develop a realistic seafloor model for the study area with realistic constraints on elastic parameters. Wavelength-scale heterogeneity in each model, in the form of seafloor roughness and subseafloor volume heterogeneity is defined using stochastic distributions with Gaussian autocorrelations. These distributions are quantified by their correlation lengths and standard deviation in amplitude. In order to incorporate all seafloor structure in a single parameterization of seafloor scattering, large-scale slope and wavelength-scale seafloor spatial parameters (rms height and correlation length), are included, along with the acoustic beam grazing-angle relative to a horizontal seafloor, in the definition of an 'effective grazing angle'. The Rayleigh roughness parameter, which depends on grazing angle of the insonification, is then redefined using the effective grazing angle and calculated for a variety of seafloor models. Scattering strengths are shown to vary systematically but nonlinearly with the 'effective Rayleigh roughness parameters' of horizontal rough seafloor models. This leads to an approximate interpretation scheme for backscatter intensity. In general, variation in backscattering is found to be dominated by the scattering from rough seafloor. If the seafloor is smooth or very low velocity (e.g., sediment), then scattering from volume heterogeneity becomes an important factor in the backscattered field. Both wavelength-scale seafloor roughness and volume heterogeneity are shown to be capable of producing the levels of variation in intensity observed in monostatic reverberation experiments. Variations in large-scale seafloor slope and subseafloor average velocity are shown to influence the backscatter response of seafloor models.

Description: My first year of study at MIT was supported by a Chevron Fellowship. The work in this thesis was funded by the Office of Naval Research under grant numbers N00014-93-1-1352, N00014-90-J-1493, and N00014-95-1-0506.

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 May 1995

Description: The data analysis of a deep-sea bottom backscattering experiment, carried out over a sediment pond on the western flank of the Mid-Atlantic Ridge in July 1993 with a 250- 650 Hz chirp source and a vertical receiving array suspended near the fiat seafloor, is presented in this thesis. Reflected signals in the normal incidence direction as the output of endfire beamforming are used to determine the sediment structure. The sediment is found to be horizontally stratified, except for two irregular regions, each about 20 m t hick, located around 18 m and 60 m beneath the water-sediment interface. Multiple constraints beamforming is shown to be effective in removing coherent reflections from internal stratified layers, which is critical to the analysis of bottom backscattering. With backscattered signals obtained by beamforming, the above-mentioned two inhomogeneous regions are found to be the dominant factors on the bottom backscattered field, both in the normal incidence and oblique directions. The backscattering strength as a function of grazing angle is estimated for each of the two regions.

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 1990

Description: This thesis introduces a new procedure for the enhancement of acoustic images of the bottom of the sea produced by side-scan sonars. Specifically, it addresses the problem of estimating and correcting geometric distortions frequently observed in such images as a consequence of motion instabilities of the sonar array. This procedure estimates the geometric distortions from the image itself, without requiring any navigational or attitude measurements. A mathematical model for the distortions is derived from the geometry of the problem, and is applied to estimates of the local degree of geometric distortion obtained by cross-correlating segments of adjacent lines of the image. The model parameters are then recursively estimated through deterministic least-squares estimation. An alternative approach based on adaptive Kalman filtering is also proposed, providing a natural framework in which a priori information about the array dynamics may be easily incorporated. The estimates of the parameters of the distortion model are used to rectify the image, and may also be used for estimating the attitude parameters of the array. A simulation is employed to evaluate the effectiveness of this technique and examples of its application to high-resolution side-scan sonar images are provided.

Description: This work was produced under sponsorship of the Conselho Nacional de Desenvolvimento Cient{jico e Tecnol6gico (CNPq), an agency of the Government of the Federative Republic of Brazil, and was supported in part by the Defense Advanced Research Projects Agency monitored by the Office of Naval Research under Grant No. N00014-89-J-1489, in part by the National Science Foundation under Grant No. MIP 87-14969, and in part by Sanders Associates, Incorporated.

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 September 1991

Description: A model is developed for the prediction of the seismo-acoustic noise spectrum in the microseism peak region (0.1 to 0.7 Hz). The model uses a theory developed by Cato [J. Acoust. Soc. Am., 89 , 1096-1112 (1991)] for an infinite depth ocean in which the surface orbital motion caused by gravity waves may produce acoustic waves at twice the gravity wave frequency. Using directional wave spectra as inputs, acoustic source levels are computed and incorporated into a more realistic environment consisting of a horizontally stratified ocean with an elastic bottom. Noise predictions are made using directional wave spectra obtained from the SWADE surface buoys moored off the coast of Virginia and the SAFARI sound propagation code, with a bottom model derived using wave speeds measured in the EDGE deep seismic reflection survey. The predictions are analyzed for noise level variations with frequency, wave height, wind direction, and receiver depth. These predictions are compared to noise measurements made in ECONOMEX using near-bottom receivers located close to the surface buoys. Good agreement is found between the predictions and observations under a variety of environmental conditions.

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 1992

Description: An ultra-short baseline acoustic navigation system has been developed which is capable of determining bearing to a sound source with an error of less than 1° in typical operational conditions. The system has a demonstrated ability to operate in an environment in which multipath interference is significant. A DSP microprocessor is used to process signals received by two hydrophones from a 26 kHz toneburst sound source. This processing power is used to implement features not commonly available with commercial systems. The system has the capability to make an on-line measurement of the signal-tonoise ratio, which can be used to estimate the confidence which should be placed in the data. Estimates of phase difference and signal power are generated many times within each received pulse, so the effects of multi path interference throughout the pulse can be observed. Results of tests at several ranges are presented, and compared to performance models developed in the thesis. System performance is quantified, and an effort is made to understand the effects of multipath arrivals.

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 1993

Description: This thesis examines the use of a single, omnidirectional hydrophone as a receiving sensor to passively localize an acoustic beacon. The localization problem is presented as a constrained, nonlinear parameter estimation problem, and Lagrange multipliers are introduced to solve for the maximum likelihood estimate of the acoustic beacon's position. An iterative algorithm is developed using range difference measurements to solve for the maximum likelihood estimate of a stationary acoustic beacon's position. This algorithm is then _extended to include linear, constant velocity motion of the acoustic beacon. Finally, design specifications for a receiver to implement the maximum likelihood estimation algorithms are developed. To test the maximum likelihood estimate algorithms, Monte Carlo simulations are conducted. Results from six representative scenarios are presented. Test results show that as the number of range differences used increases, or the distance that the observer travels between received beacon signals increases, the accuracy of the estimated position improves. Also, tests show that accuracy of the estimated beacon position is directly related to the accuracy in which the observer's position is measured. To test the receiver's design specifications, a prototype receiver is built using commonly available components. It is then shown that the prototype receiver meets or exceeds the design specifications.

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

Description: This thesis presents the design and simulation of an acoustic listen-only navigation system for use in a time-varying high-multipath environment. The system uses a ray tracing model of sound propagation in a horizontally stratified ocean to associate each multipath arrival with a particular ray path through the ocean. This ray path identification allows an inversion for both vehicle position and sound speed profile parameters. Tracking changes in the sound speed profile allows the vehicle to maintain navigation accuracy in a time-varying acoustic environment. Simulation results for typical Arctic conditions indicate the potential for accurate navigation from acoustic beacons at ranges up to tens of kilometers.

Description: I also wish to acknowledge the support received from a National Science Foundation Graduate Fellowship under Grant #RCD-9154652 and from the Office of Naval Research Arctic Programs Office, which made possible this work.

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

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

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

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

Description: Ambient noise in the sea has been observed for over 100 years. Previous studies conclude that the primary source of microseisms is nonlinear interaction of surface gravity waves at the sea surface. Though this source relationship is generally accepted, the actual processes by which the wave generated acoustic noise in the water column couples and propagates to and along the sea floor are not well understood. In this thesis, the sources and propagation of sea floor and sub-sea floor microseismic noise between 0.2 and 10 Hz are investigated. This thesis involves a combination of theoretical, observational and numerical analysis to probe the nature of the microseismic field in the Blake Bahama Basin. Surface waves are the primary mechanism for noise propagation in the crust and fall into two separate groups depending on the relative wavelength/water depth ratio. Asymptotic analysis of the Sommerfeld integral in the complex ray parameter plane shows results that agree with previous findings by Strick (1959) and reveal two fundamental interface wave modes for short wavelength noise propagation in the crust: the Stoneley and pseudo-Rayleigh wave. For ocean sediments, where the shear wave velocity is less than the acoustic wave velocity of water, only the Stoneley interface wave can exist. For well consolidated sediments and basalt, the shear velocity exceeds the acoustic wave velocity of water and the pseudo-Rayleigh wave can also exist. Both interface waves propagate with retrograde elliptic motion at the sea floor and attenuate with depth into the crust, however the pseudo-Rayleigh wave travels along the interface with dispersion and attenuation and "leaks" energy into the water column for a half-space ocean over elastic crust model. For finite depth ocean models, the pseudo-Rayleigh wave is no longer leaky and approaches the Rayleigh wave velocity of the crust. The analysis shows that longer wavelength noise propagates as Rayleigh and Stoneley modes of the ocean+crust waveguide. These long wavelength modes are the fundamental mechanism for long range noise propagation. During the Low Frequency Acoustic Seismic Experiment (LFASE) a four-node, 12- channel borehole array (SEABASS) was deployed in the Blake Bahama Basin off the coast of eastern Florida (DSDP Hole 534B). This experiment is unique and is the first use of a borehole array to measure microseismic noise below the sea floor. Ambient background noise from a one week period is compared between an Ocean Bottom Seismometer (OBS) and SEABASS at sub-bottom depths of 10, 40, 70 and 100 meters below the sea floor. The 0.3 H z microseism peak is found to be nearly invariant with depth and has a power level of 65 and 75 dB rel 1 (nm/ s2)2)/ H z for the vertical and horizontal components respectively. At 100 m depth, the mean microseismic noise levels above 0.7 Hz are 10 dB and 15-20 dB quieter for the vertical and horizontal components respectively. Most of this attenuation occurs in the upper 10 m above 1.0 Hz, however higher modes in the spectra show narrow bandwidth variability in the noise field that is not monotonic with depth. Dispersion calculations show normal mode Stoneley waves below 0.7 Hz and evidence of higher modes above 0.8 Hz. A strong correlation between noise levels in the borehole and local sea state conditions is observed along with clear observation of the nonlinear frequency doubling effect between ocean surface waves and microseisms. Particle motion analysis further verifies that noise propagates through the array as Rayleigh/Stoneley waves. Polarization direction indicates at least two sources; distant westerly swell during quiescent times and local surface waves due to a passing storm. Above 1.0 Hz the LFASE data shows little coherence and displays random polarization. Because of this, we believe scattered energy is a significant component of the noise field in the Blake Bahama Basin. A fully 3-D finite difference algorithm is used to model both surface and volume heterogeneities in the ocean crust. Numerical modeling of wave propagation for hard and soft bottom environments shows that heterogeneities on the order of a seismic wavelength radiate energy into the water column and convert acoustic waves in the water into small wavelength Stoneley waves observed at the borehole. Sea floor roughness is the most important elastic scattering feature of the ocean crust. Comparisons of 2D and 3D rough sea floor models show that out-of-plane effects necessitate the use of 3D methods. The out-of-plane energy that is present in the LFASE data comes from either heterogeneities in the source field (i.e. mixed gravity wave directions) or, equally likely, scattering of the source field from surface or volume heterogeneities in the sea floor.

Description: This research was supported by Office of Naval Research grants N00014-89-C-0018, N00014-89-J-1012, N00014-90-C-0098, N00014-90-J-1493 and N00014-93-1-1352.

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

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

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

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

Description: This paper develops a beam pattern design procedure for general multidimensional irregular sonar arrays that incorporates the not well understood effects of array geometry into the design process. The procedure is implemented by generating a "penalty function" in a spectral covariance function form. Processing the penalty function causes beam pattern high sidelobes to be penalized and the main lobe to be emphasized. This is accomplished by forming the penalty function in terms of an isotropic noise field of specified strength modified with a finite sector of low coherent energy and stabilized with incoherent sensor noise. By inputting the penalty function into a minimum variance beamformer, the beam pattern and aperture weights are calculated based on the given array geometry. The beamformer used is Capon's Maximum Likelihood Method. The array used to test the procedure is located on a sixty degree sector of a cylindrical surface. The procedure is implemented by two different methods, each with some desirable characteristics. One method suppresses sidelobes directly by the placement of nulls. The other method suppresses sidelobes indirectly by the enhancement of the main lobe with anti-nulls. Both methods are evaluated in terms of a sensitivity factor which constrains the maximum white noise array gain. Results show that both methods result in sidelobe levels that range from 20 to 35 dB lower compared to a conventional beam pattern with uniform aperture weighting and that the design procedure is applicable to beam patterns steered to both true broadside and to off-broadside directions.

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

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

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

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

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

Description: 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: Signal detection and synchronization in the time varying ocean environment is a difficult endeavor. The current common methods include using a linear frequency modulated chirped pulse or maximal length sequence as a detection pulse, then match filtering to that signal. In higher signal to noise ratio (SNR) environments (~0 dB and higher) this has been a suitable solution. As the SNR drops lower however, this solution no longer provides an acceptable probability of detection for a given tolerable probability of false alarm. The issue derives from the inherent coherence issues in the ocean environment which limit the useful matched filter length. This thesis proposes an alternative method of detection based on a recursive least squares linearly adaptive equalizer which we term the Adaptive Linear Equalizer Detector (ALED). This detectors performance has demonstrated reliable probability of detection with minimal interfering false alarms with SNR as low as -20 dB. Additionally this thesis puts forth a computationally feasible method for implementing the detector.

Description: Support from the Office of Naval Research (through ONR grant #N00014-07-10738 and #N00014-11-10426).

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

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

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

Description: Author Posting. © American Meteorological Society, 2012. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 42 (2012): 1684–1700, doi:10.1175/JPO-D-11-0230.1.

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

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

Description: 2013-04-01

Description: Author Posting. © American Meteorological Society, 2011. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 41 (2011): 889–910, doi:10.1175/2010JPO4496.1.

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

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

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

Description: During July and August of 1996, a large acoustics/physical oceanography experiment was fielded in the Mid-Atlantic Bight, south of Nantucket Island, MA. Known as the Shelfbreak Front PRIMER Experiment, the study combined acoustic data from a moored array of sources and receivers with very high resolution physical oceanographic measurements. This thesis addresses two of the primary goals of the experiment, explaining the properties of acoustic propagation in the region, and tomographic inversion of the acoustic data. In addition, this thesis develops a new method for predicting acoustic coherence in such regions. Receptions from two 400 Hz tomography sources, transmitting from the continental slope onto the shelf, are analyzed. This data, along with forward propagation modeling utilizing SeaSoar thermohaline measurements, reveal that both the shelfbreak front and tidally-generated soliton packets produce stronger coupling between the acoustic waveguide modes than expected. Arrival time wander and signal spread show variability attributable to the presence of a shelf water meander, changes in frontal configuration, and variability in the soliton field. The highly-coupled nature of the acoustic mode propagation prevents detailed tomographic inversion. Instead, methods based on only the wander of the mode arrivals are used to estimate path-averaged temperatures and internal tide "strength". The modal phase structure function is introduced as a useful proxy for acoustic coherence, and is related via an integral transform to the environmental sound speed correlation function. Advantages of the method are its flexibility and division of the problem into independent contributions, such as from the water column and seabed.

Description: Office of Naval Research for providing funding for this thesis through AASERT Grant N00014-96-1-0918, and through ONR Grant N00014-98-1-0059.

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

Description: In this thesis the analysis of natural ice events is carried out based on direct measurements of ice-borne seismo-acoustic waves generated by ice fracturing processes. A major reason for studying this phenomenon is that this acoustic emission is a significant contributor to Arctic ocean ambient noise. Also the Arctic contains rich mineral and oil resources and in order to design mining facilities able to withstand the harsh environmental conditions, we need to have a better understanding of the processes of sea ice mechanics. The data analyzed in this thesis were collected during the Sea Ice Mechanics Initiative SIMI’94 experiment which was carried out in the spring of 1994 in the Central Arctic. One of the contributions of this thesis was the determination of the polarization characteristics of elastic waves using multicomponent geophone data. Polarization methods are well known in seismology, but they have never been used for ice event data processing. In this work one of the polarization methods so called Motion Product Detector method has been successfully applied for localization of ice events and determination of polarization characteristics of elastic waves generated by fracturing events. This application demonstrates the feasibility of the polarization method for ice event data processing because it allows one to identify areas of high stress concentration and "hot spots" in ridge building process. The identification of source mechanisms is based on the radiation patterns of the events. This identification was carried out through the analysis of the seismo-acoustic emission of natural ice events in the ice sheet. Previous work on natural ice event identification was done indirectly by analyzing the acoustic energy radiated into the water through coupling from elastic energy in the ice sheet. After identification of the events, the estimation of the parameters of fault processes in Arctic ice is carried out. Stress drop, seismic moment and the type of ice fracture are determined using direct near-field measurements of seismo-acoustic signals generated by ice events. Estimated values of fracture parameters were in good agreement with previous work for marginal ice zone. During data processing the new phenomenon was discovered: "edge waves", which are waves propagating back and forth along a newly opened ice lead. These waves exhibit a quasi-periodic behavior suggesting some kind of stick-slip generation mechanism somewhere along the length of the lead. The propagation characteristics of these waves were determined using seismic wavenumber estimation techniques. In the low frequency limit the dispersion can be modeled approximately by an interaction at the lead edges of the lowest order, antisymmetric modes of the infinite plate.

Description: Support for this thesis was provided by Office of Naval Research.

Description: Author Posting. © American Meteorological Society, 2011. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 24 (2011): 4844–4858, doi:10.1175/2011JCLI4130.1.

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

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

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

Description: An algorithm is developed for underwater acoustic signal processing with an array of hydrophones. With various acoustic signals coming from different directions, the maximum likelihood approach is used to estimate the source bearings and time series. Simulated annealing is used to implement the resulting time-domain beamformer. Broadband signals in spatially correlated noise are treated. Previous time-domain beamformers did not consider the correlation between random noise, and they did not use the concept of maximum likelihood, which is asymptotically optimal. We show that improved resolution can be achieved using this new method.

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 June 1988

Description: Two approaches to determining the ocean sound speed profile using measured acoustic modal eigenvalues are examined. Both methods use measured eigenvalues and mode dependent assumed values of the WKB phase integral as input data and use the WKB phase integral as a starting point for relating the index of refraction to depth. Inversion method one is restricted to monotonic or symmetric sound speed profiles and requires a measurement of the sound speed at one depth to convert the index of refraction profile to a sound speed profile. Inversion method two assumes that the sound speed at the surface and the minimum sound speed in the profile are known and is applicable to monotonic profiles and to general single duct sound speed profiles. For asymmetric profiles, inversion method two gives the depth difference between two points of equal sound speed in the portion of the profile having two turning points, and in the remainder of the profile it gives sound speed versus depth directly. A numerical implementation of the methods is demonstrated using idealized ocean sound speed profiles numerical experiments used to test the performance of the inversions using noisy data. The two methods are used to determine the sediment sound speed profiles in two shallow water waveguide models, and inversion method one is used to find the sediment sound speed profile using data from an experiment performed in the Gulf of Mexico.

Description: Funding through the Office of Naval Research Fellowship Program.

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

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

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

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

Description: Current underwater acoustic channel estimation techniques generally apply linear MMSE estimation. This approach is optimal in a mean square error sense under the assumption that the impulse response fluctuations are well characterized by Gaussian statistics, leading to a Rayleigh distributed envelope. However, the envelope statistics of the underwater acoustic communication channel are often better modeled by the K-distribution. In this thesis, by presenting and analyzing field data to support this claim, I demonstrate the need to investigate channel estimation algorithms that exploit K-distributed fading statistics. The impact that environmental conditions and system parameters have on the resulting distribution are analyzed. In doing so, the shape parameter of the K-distribution is found to be correlated with the source-to-receiver distance, bandwidth, and wave height. Next, simulations of the scattering behavior are carried out in order to gain insight into the physical mechanism that cause these statistics to arise. Finally, MAP and MMSE based algorithms are derived assuming K-distributed fading models. The implementation of these estimation algorithms on simulated data demonstrates an improvement in performance over linear MMSE estimation.

Description: This work was supported by the Office of Naval Research grant #N00014-05-10085 and the National Science Foundation grant #OCE-0519903.

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

Description: Several experiments are presented in this thesis which examine methods to measure and monitor fluid flow from hydrothermal vent fields. Simultaneous velocity, temperature, and conductivity data were collected in the convective flow emanating from a hydrothermal vent field located at 10°S6'N, 103° 41'W on the East Pacific rise. The horizontal profiles obtained indicate that the flow field approaches an ideal plume in the temperature and velocity distribution. Such parameters as total heat flow and maximum plume height can be estimated using either the velocity or the temperature information. The results of these independent calculations are in close agreement, yielding a total heat flow from this vent site of 3.7 ± 0.8 MW and a maximum height of 150±10 m. The nonlinear effects of large temperature variations on heat capacity and volume changes slightly alter the calculations applied to obtain these values. In Guaymas Basin, a twelve day time series of temperature data was collected from a point three centimeters above a diffuse hydrothermal flow area. Using concurrent tidal gauge data from the town of Guaymas it is shown that the effects of tidar currents can be strong enough to dominate the time variability of a temperature signal at a fixed point in hydrothermal flow and are a plausible explanation for the variations seen in the Guaymas Basin temperature data. Theoretical examination of hot, turbulent, buoyant jets exiting from hydrothermal chimneys revealed acoustic source mechanisms capable of producing sound at levels higher than ambient ocean noise. Pressure levels and frequency generated by hydrothermal jets are dependent on chimney dimensions, fluid velocity and temperature and therefore can be used to monitor changes in these parameters over time. A laboratory study of low Mach number jet noise and amplification by flow inhomogeneities confirmed theoretical predictions for homogeneous jet noise power and frequency. The increase in power due to convected flow inhomogeneities, however, was lower in the near field than expected. Indirect evidence of hydrothermal sound fields (Reidesel et al., 1982; Bibee and Jacobson, 1986) showing anomalous high power and low frequency noise associated with vents is due to processes other than jet noise. On Axial Seamount, Juan de Fuca Ridge, high quality acoustic noise measurements were obtained by two hydrophones located 3 m and 40 m from an active hydrothermal vent, in an effort to determine the feasibility of monitoring hydrothermal vent activity through floW noise generation. Most of the noise field could be attributed to ambient ocean noise sources of microseisms, distant shipping and weather, punctuated by local ships and biological sources. Water/rock interface waves of local origin, were detected which showed high pressure amplitudes near the seafloor and, decaying with vertical distance, produced low pressures at 40 m above the bottom. Detection of vent signals was hampered by unexpected spatial non- stationarity due to shadowing effects of the caldera wall. No continuous vent signals were deemed significant based on a criterion of 90% probability of detection and 5% probability of false alarm. However, a small signal near 40 Hz, with a power level of 1x10-4 Pa2/Hz was noticed on two records taken near the Inferno black smoker. The frequency of this signal is consistent with predictions and the power level suggests the occurrence of jet noise amplification due to convected density inhomogeneities. Ambient noise from the TAG (Trans-Atlantic Geotraverse) hydrothermal area on the Mid-Atlantic Ridge near 26°N, in the frequency band 1-30 Hz at a range of 0.75-14 km from the site of an extremely active high temperature hydrothermal vent field (Rona, 1986) was examined. The ambient noise field exhibits great temporal and spatial variations attributed in part to typical ocean noise sources such as distant shipping and microseisms. Power spectral levels as measured at each of six ocean bottom hydrophones (OBH) were used to estimate the location of point sources of sound in the area, if any. The hydrothermal vent did not produce enough sound to be located as a point source using data from the OBH array. The only consistently identifiable point source found with the data set was generating sound in a 0.8-3.5 Hz bandwidth and located outside the median valley. It appears to be harmonic tremor associated with the tip of a ridge on the western side of the spreading axis and may be volcanic in origin.

Description: This work was supported by the WHOI/MIT Education Office, the Center for Analysis of Marine Systems, the National Science Foundation (grant OCE83-l0l75), NOAA National Sea Grant College Program Office, Dept. of Commerce under grant #NA86-AA-D-SG090, WHOI Sea Grant (R/6-l4), the Office of Naval Research grant #N0014-87-K-0007, and the NOAA Vents Program.

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

Description: In practical applications with bathymetric sidescan sonars, the multipath reflections and other directional interferences are the key limiting factors for a better performance. This thesis proposes a new scheme to deal with the interferences using a multiple-row bathymetric sidescan sonar. Instead of smoothing the measurements over some time or angle intervals, which was previously widely investigated, we resolve the multipath interferences from the direct signal. Two approaches on signal direction-of-arrival DOA and amplitude estimation are developed, the correlated signal direction estimate CSDE for three-row systems and the ESPRIT-based method. These approaches are compared using different sonar data models, including a stochastic model from the statistical analysis on bottom scattering and a coherent model from the analysis on interference field; the simulations show the ESPRIT-based approach is quite robust at the angular separation of 100 between two sources and at the signal-to-noise ratio above 10dB except for highly coherent or temporally correlated signals, for which CSDE works very well. The computer simulation results and the discussions on practical algorithm implementation indicate the proposed scheme can be applied to a real multiple-row bathymetric sidescan sonar. With the capability to simultaneously resolve two or more directional signals, the new sonar model should work better for a wider variety of practical situations in shallow water with out significant increase of the system cost.

Description: Funding supporting my thesis research project was provided by the Office of Naval Research ONR.

Description: Author Posting. © American Meteorological Society, 2011. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 41 (2011): 1182–1208, doi:10.1175/2010JPO4564.1.

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

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

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

Description: In order to determine the efficacy of tomographic reconstructions of the ocean sound speed structure in improving acoustic field predictions for source localization, a 150 km by 350 km volume of ocean 3000 meters deep was synthetically modeled to be similar to the Gulf Stream system, including an eddy and a front. The features were Gaussian, with the eddy's maximum sound speed perturbation being 10ms-1 and the front's maximum perturbation 15ms-1. Two vertical slices through this system were inverted in a synthetic tomography experiment using linear optimal estimation theory. Inversions were also performed using XSV and satellite sea surface temperature data. Gaussian fits to the reconstructed features were constructed for use with a three dimensional raytrace program (HARPO). Three dimensional rays were propagated both through the reconstructions and the original model. Travel time versus intensity (transmission loss) for the eigenrays was used as a basis for intercomparison. Tomographic results showed good reconstruction for a first iteration of the inversion, but inadequate vertical resolution. Iterations and the use of more refractive eigenrays are needed for improvement of the reconstruction, especially for the front. Reconstructed results for the acoustic field should improve conventional beamforming, but are probably inadequate for matched field processing.

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

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

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

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

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

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