ALBERT

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 2009

Description: Observations and inverse models suggest that small-scale turbulent mixing is enhanced in the Southern Ocean in regions above rough topography. The enhancement extends 1 km above the topography suggesting that mixing is supported by breaking of gravity waves radiated from the ocean bottom. In other regions, gravity wave radiation by bottom topography has been primarily associated with the barotropic tide. In this study, we explore the alternative hypothesis that the enhanced mixing in the Southern Ocean is sustained by internal waves generated by geostrophic motions flowing over bottom topography. Weakly-nonlinear theory is used to describe the internal wave generation and the feedback of the waves on the zonally averaged flow. A major finding is that the waves generated at the ocean bottom at finite inverse Froude numbers drive vigorous inertial oscillations. The wave radiation and dissipation at equilibrium is therefore the result of both geostrophic flow and inertial oscillations and differs substantially from the classical lee wave problem. The theoretical predictions are tested versus two-dimensional and three-dimensional high resolution numerical simulations with parameters representative of the Drake Passage region. Theory and fully nonlinear numerical simulations are used to estimate internal wave radiation from LADCP, CTD and topography data from two regions in the Southern Ocean: Drake Passage and the Southeast Pacific. The results show that radiation and dissipation of internal waves generated by geostrophic motions reproduce the magnitude and distribution of dissipation measured in the region.

Description: Author Posting. © American Meteorological Society, 2010. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 40 (2010): 789-801, doi:10.1175/2009JPO4039.1.

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

Description: Author Posting. © American Meteorological Society, 2008. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 38 (2008): 2556-2574, doi:10.1175/2008JPO3666.1.

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

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

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

Description: The thesis develops and demonstrates methods of classifying ocean processes using an underwater moving platform such as an Autonomous Underwater Vehicle (AUV). The "mingled spectrum principle" is established which concisely relates observations from a moving platform to the frequency-wavenumber spectrum of the ocean process. It clearly reveals the role of the AUV speed in mingling temporal and spatial information. For classifying different processes, an AUV is not only able to jointly utilize the time-space information, but also at a tunable proportion by adjusting its cruise speed. In this respect, AUVs are advantageous compared with traditional oceanographic platforms. Based on the mingled spectrum principle, a parametric tool for designing an AUVbased spectral classifier is developed. An AUV's controllable speed tunes the separability between the mingled spectra of different processes. This property is the key to optimizing the classifier's performance. As a case study, AUV-based classification is applied to distinguish ocean convection from internal waves. The mingled spectrum templates are derived from the MIT Ocean Convection Model and the Garrett-Munk internal wave spectrum model. To allow for mismatch between modeled templates and real measurements, the AUVbased classifier is designed to be robust to parameter uncertainties. By simulation tests on the classifier, it is demonstrated that at a higher AUV speed, convection's distinct spatial feature is highlighted to the advantage of classification. Experimental data are used to test the AUV-based classifier. An AUV-borne flow measurement system is designed and built, using an Acoustic Doppler Velocimeter (ADV). The system is calibrated in a high-precision tow tank. In February 1998, the AUV acquired field data of flow velocity in the Labrador Sea Convection Experiment. The Earth-referenced vertical flow velocity is extracted from the raw measurements. The classification test result detects convection's occurrence, a finding supported by more traditional oceanographic analyses and observations. The thesis work provides an important foundation for future work in autonomous detection and sampling of oceanographic processes.

Description: This thesis research has been funded by the Office of Naval Research (ONR) under Grants NOOOl4-95-1-1316, NOO0l4-97-1-0470, and by the MIT Sea Grant College Program under Grant NA46RG0434.

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

Description: Steady scale-invariant solutions of a kinetic equation describing the statistics of oceanic internal gravity waves based on wave turbulence theory are investigated. It is shown in the nonrotating scale-invariant limit that the collision integral in the kinetic equation diverges for almost all spectral power-law exponents. These divergences come from resonant interactions with the smallest horizontal wavenumbers and/or the largest horizontal wavenumbers with extreme scale separations. A small domain is identified in which the scale-invariant collision integral converges and numerically find a convergent power-law solution. This numerical solution is close to the Garrett–Munk spectrum. Power-law exponents that potentially permit a balance between the infrared and ultraviolet divergences are investigated. The balanced exponents are generalizations of an exact solution of the scale-invariant kinetic equation, the Pelinovsky–Raevsky spectrum. A small but finite Coriolis parameter representing the effects of rotation is introduced into the kinetic equation to determine solutions over the divergent part of the domain using rigorous asymptotic arguments. This gives rise to the induced diffusion regime. The derivation of the kinetic equation is based on an assumption of weak nonlinearity. Dominance of the nonlocal interactions puts the self-consistency of the kinetic equation at risk. However, these weakly nonlinear stationary states are consistent with much of the observational evidence.

Description: This research is supported by NSF CMG Grants 0417724, 0417732 and 0417466. YL is also supported by NSF DMS Grant 0807871 and ONR Award N00014-09-1-0515.

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: 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): 1524–1547, doi:10.1175/JPO-D-11-0117.1.

Description: Evidence is presented for the transfer of energy from low-frequency inertial–diurnal internal waves to high-frequency waves in the band between 6 cpd and the buoyancy frequency. This transfer links the most energetic waves in the spectrum, those receiving energy directly from the winds, barotropic tides, and parametric subharmonic instability, with those most directly involved in the breaking process. Transfer estimates are based on month-long records of ocean velocity and temperature obtained continuously over 80–800 m from the research platform (R/P) Floating Instrument Platform (FLIP) in the Hawaii Ocean Mixing Experiment (HOME) Nearfield (2002) and Farfield (2001) experiments, in Hawaiian waters. Triple correlations between low-frequency vertical shears and high-frequency Reynolds stresses, uiw∂Ui/∂z, are used to estimate energy transfers. These are supported by bispectral analysis, which show significant energy transfers to pairs of waves with nearly identical frequency. Wavenumber bispectra indicate that the vertical scales of the high-frequency waves are unequal, with one wave of comparable scale to that of the low-frequency parent and the other of much longer scale. The scales of the high-frequency waves contrast with the classical pictures of induced diffusion and elastic scattering interactions and violates the scale-separation assumption of eikonal models of interaction. The possibility that the observed waves are Doppler shifted from intrinsic frequencies near f or N is explored. Peak transfer rates in the Nearfield, an energetic tidal conversion site, are on the order of 2 × 10−7 W kg−1 and are of similar magnitude to estimates of turbulent dissipation that were made near the ridge during HOME. Transfer rates in the Farfield are found to be about half the Nearfield values.

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

Description: 2013-03-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 January 18, 1988

Description: We have contrived a regional model Φ(K, ω, η, φ, λ) for the distribution of low frequency variability energy in horizontal wavenumber, frequency, vertical mode and geography. We assume horizontal isotropy, Φ(K, ω, η, φ, λ) = 2πKψ(k, l, ω, η, φ, λ), with K designating the amplitude of total horizontal wavenumber. The parameters of Φ(K, ω, η, φ, λ) can be derived from observations: (i) satellite altimetry measurements yield the surface eddy kinetic energy wavenumber and frequency spectra and the geographic distribution of surface eddy kinetic energy magnitude, (ii) XBT measurements yield the temperature wavenumber spectra, (iii) current meter and thermistor measurements yield the frequency spectra of kinetic energy and temperature, (iv) tomographic measurements yield the frequency spectra of range— and depth—averaged temperature, and (v) the combination of satellite altimetry and current meter measurements yields the vertical partitioning of kinetic energy among dynamical modes. We assume the form of the geography—independent part of our model Φ(K, ω, η) ∝Kpωq. The observed kinetic energy and temperature wavenumber spectra suggest p = 3/2 at K 〈 K0 and p = —2 at K 〉 K0 for the barotropic mode, and p = —1/2 at K 〈 K0 and p = —3 at K 〉 K0 for the baroclinic mods, where K0 is the transitional wavenumber of the wavenumber spectra. The observed frequency spectra of temperature and kinetic energy suggest that q = —1/2 for ω 〈 ω0 and q = —2 for ω 〉 ω0, where ω0 is the transitional frequency of the frequency spectra. The combination of satellite altimetry and current meter measurements suggests the vertical structure of the low frequency variability is governed by the first few modes. The geography—dependent part of our model is the energy magnitude. Although we have shown analytically that the tomographic measurements behave as a low—pass filter, it is impossible to identify this filtering effect in the real data due to the strong geographic variability of the energy magnitude and the vertical gradient of the mean temperature. The model wavenumber spectrum is appropriate only where the statistical properties are relatively homogeneous in space.

Description: My first year in the Joint Program was supported by the National Science Foundation under grant OCE 92-16628, then were supported by the National Science Foundation under grant OCE 95-29545.

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): 17–28, doi:10.1175/JPO-D-11-0108.1.

Description: Observational evidence is presented for transfer of energy from the internal tide to near-inertial motions near 29°N in the Pacific Ocean. The transfer is accomplished via parametric subharmonic instability (PSI), which involves interaction between a primary wave (the internal tide in this case) and two smaller-scale waves of nearly half the frequency. The internal tide at this location is a complex superposition of a low-mode waves propagating north from Hawaii and higher-mode waves generated at local seamounts, making application of PSI theory challenging. Nevertheless, a statistically significant phase locking is documented between the internal tide and upward- and downward-propagating near-inertial waves. The phase between those three waves is consistent with that expected from PSI theory. Calculated energy transfer rates from the tide to near-inertial motions are modest, consistent with local dissipation rate estimates. The conclusion is that while PSI does befall the tide near a critical latitude of 29°N, it does not do so catastrophically.

Description: This work was sponsored by NSF OCE 04-25283.

Description: 2013-07-01

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

Description: The North Atlantic Oscillation (NAO) is one of the most important modes of variability in the global climate system and is characterized by a meridional dipole in the sea level pressure field, with centers of action near Iceland and the Azores. It has a profound influence on the weather, climate, ecosystems, and economies of Europe, Greenland, eastern North America, and North Africa. It has been proposed that around 1980, there was an eastward secular shift in the NAO’s northern center of action that impacted sea ice export through Fram Strait. Independently, it has also been suggested that the location of its southern center of action is tied to the phase of the NAO. Both of these attributes of the NAO have been linked to anthropogenic climate change. Here the authors use both the one-point correlation map technique as well as empirical orthogonal function (EOF) analysis to show that the meridional dipole that is often seen in the sea level pressure field over the North Atlantic is not purely the result of the NAO (as traditionally defined) but rather arises through an interplay among the NAO and two other leading modes of variability in the North Atlantic region: the East Atlantic (EA) and the Scandinavian (SCA) patterns. This interplay has resulted in multidecadal mobility in the two centers of action of the meridional dipole since the late nineteenth century. In particular, an eastward movement of the dipole has occurred during the 1930s to 1950s as well as more recently. This mobility is not seen in the leading EOF of the sea level pressure field in the region.

Description: GWKM was supported by the Natural Sciences and Engineering Research Council of Canada. IAR was supported in part by NE/C003365/1. RSP was supported by Grant OCE-0959381 from the U.S. National Science Foundation.

Description: 2013-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 April, 1977

Description: A total of four moorings from POLYMODE Array I and II were analyzed in an investigation of internal wavefield-mean flow interactions. In particular, evidence for wave-mean flow interaction was sought by searching for time correlations between the wavefield vertically-acting Reynolds stress (estimated using the temperature and velocity records), and the mean shear. No significant stress-shear correlations were found at the less energetic moorings, indicating that the magnitude of the eddy viscosity was under 200 cm2/sec, with the sign of the energy transfer uncertain. This is considerably below the 0(4500 cm2/sec) predicted by Müller (1976). An extensive error analysis indicates that the large wave stress predicted by the theory should have been clearly observable under the conditions of measurement. Theoretical computations indicate that the wavefield "basic state" may not be independent of the mean flow as assumed by Müller, but can actually be modified by large-scale vertical shear and still remain in equilibrium. In that case, the wavefield does not exchange momentum with a large-scale vertical shear flow, and, excepting critical layer effects, a small vertical eddy viscosity is to be expected. Using the Garrett-Munk (1975) model internal wave spectrum, estimates were made of the maximum momentum flux (stress) expected to be lost to critical layer absorption. Stress was found to increase almost linearly with the velocity difference across the shear zone, corresponding to a vertical eddy viscosity of -100 cm2 s -1. Stresses indicative of this effect were not observed in the data. The only significantly non-zero stress correlations were found at the more energetic moorings. Associated with the 600 m mean velocity and the shear at the thermocline were a positively correlated stress at 600 m, and a negatively correlated stress at 1000 m. These stress correlations were most clearly observable in the frequency range corresponding to 1 to 8 hour wave periods. The internal wavefield kinetic and potential energy were modulated by the mean flow at both levels, increasing by a factor of two with a factor of ten in the mean flow. The observed stress correlations and energy level changes were found to be inconsistent with ideas of a strictly local eddy viscosity, in which the spectrum of waves is only slightly modified by the shear. When Doppler effects in the temperature equation used to estimate vertical velocity were considered, the observations of stress and energy changes were found to be consistent with generation of short (0.4 to 3 km) internal waves at the level of maximum shear, about 800 m. The intensity of the generated waves increases with the shear, resulting in an effective vertical eddy viscosity (based on the main thermocline shear) of about +100 cm2 s-1 The stresses were not observable at the 1500 m level, indicating that the waves were absorbed within 500 m of vertical travel. The tendency for internal wave currents to be horizontally anisotropic in the presence of a mean current was investigated. Using the Garrett- Munk (1975) model internal wave spectrum, it was found that critical layer absorption cannot induce anisotropies as large as observed. A mechanical noise problem was found to be the cause of large anisotropies measured with Geodyne 850 current meters. It could not be decided, however, whether or not the A.M.F. Vector Averaging Current Meter is able to satisfactorily remove the noise with its averaging scheme.

Description: The research reported here was provided by Office of Naval Research Contract Numer N00014-76-C-0197 NR 083-400.

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: Intermittent, shoreward propagating packets of high frequency first mode internal waves are common on the continental shelf when the water column is stratified and may induce large fluctuations in near bottom velocity. Simple theoretical considerations here lead to an approximate method for estimating those quantities of most interest for the bottom boundary layer interaction problem. Examination of data from the pilot Coastal Ocean Dynamics Experiment (CODE I) shows that near bottom velocity fluctuations in the high frequency internal wave band were dominated by shoreward propagating, intermittent mode 1 internal events. Predictions of CODE I internal wave characteristics using the above approximate method are shown to be good. A boundary layer model is developed, which allows for the nonlinear interaction of surface waves, internal waves, and a steady current over a rough bottom. Modeling results suggest that internal waves will significantly enhance the stress felt by the steady current, and can increase the variability and decrease the reliability of boundary layer measurements by the "log profile" technique, when the waves are present. Theoretical dissipation of internal wave energy in the bottom boundary layer is found to be significantly enhanced in the presence of surface waves and currents, and may be important to the overall internal wave energy balance on the shelf.

Description: My doctoral work was supported for the first three years by an NSF Graduate Fellowship and has been supported since under NSF grant OCE-8014938.

Description: Author Posting. © American Meteorological Society, 2008. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 38 (2008): 880–895, doi:10.1175/2007JPO3750.1.

Description: The oceanic response to overflows is explored using a two-layer isopycnal model. Overflows enter the open ocean as dense gravity currents that flow along and down the continental slope. While descending the slope, overflows typically double their volume transport by entraining upper oceanic water. The upper oceanic layer must balance this loss of mass, and the resulting convergent flow produces significant vortex stretching. Overflows thus represent an intense and localized mass and vorticity forcing for the upper ocean. In this study, simulations show that the upper ocean responds to the overflow-induced forcing by establishing topographic β plumes that are aligned more or less along isobaths and that have a transport that is typically a few times larger than that of the overflows. For the topographic β plume driven by the Mediterranean overflow, the occurrence of eddies near Cape St. Vincent, Portugal, allows the topographic β plume to flow across isobaths. The modeled topographic β-plume circulation forms two transatlantic zonal jets that are analogous to the Azores Current and the Azores Countercurrent. In other cases (e.g., the Denmark Strait overflow), the same kind of circulation remains trapped along the western boundary and hence would not be readily detected.

Description: SK’s support during the time of his Ph.D. research in the MIT/WHOI Joint Program was provided by the National Science Foundation through Grant OCE04-24741. JP and JY have also received support from the Climate Process Team on Gravity Current Entrainment, NSF Grant OCE-0611530.

Description: Author Posting. © American Meteorological Society, 2008. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 38 (2008): 380-399, doi:10.1175/2007JPO3728.1.

Description: Barotropic to baroclinic conversion and attendant phenomena were recently examined at the Kaena Ridge as an aspect of the Hawaii Ocean Mixing Experiment. Two distinct mixing processes appear to be at work in the waters above the 1100-m-deep ridge crest. At middepths, above 400 m, mixing events resemble their open-ocean counterparts. There is no apparent modulation of mixing rates with the fortnightly cycle, and they are well modeled by standard open-ocean parameterizations. Nearer to the topography, there is quasi-deterministic breaking associated with each baroclinic crest passage. Large-amplitude, small-scale internal waves are triggered by tidal forcing, consistent with lee-wave formation at the ridge break. These waves have vertical wavelengths on the order of 400 m. During spring tides, the waves are nonlinear and exhibit convective instabilities on their leading edge. Dissipation rates exceed those predicted by the open-ocean parameterizations by up to a factor of 100, with the disparity increasing as the seafloor is approached. These observations are based on a set of repeated CTD and microconductivity profiles obtained from the research platform (R/P) Floating Instrument Platform (FLIP), which was trimoored over the southern edge of the ridge crest. Ocean velocity and shear were resolved to a 4-m vertical scale by a suspended Doppler sonar. Dissipation was estimated both by measuring overturn displacements and from microconductivity wavenumber spectra. The methods agreed in water deeper than 200 m, where sensor resolution limitations do not limit the turbulence estimates. At intense mixing sites new phenomena await discovery, and existing parameterizations cannot be expected to apply.

Description: This work was funded by the National Science Foundation as a component of the Hawaii Ocean Mixing Program. Added support for FLIP was provided by the Office of Naval Research.

Description: Author Posting. © American Meteorological Society, 2008. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 38 (2008): 686-701, doi:10.1175/2007JPO3826.1.

Description: The disintegration of a first-mode internal tide into shorter solitary-like waves is considered. Since observations frequently show both tides and waves with amplitudes beyond the restrictions of weakly nonlinear theory, the evolution is studied using a fully nonlinear, weakly nonhydrostatic two-layer theory that includes rotation. In the hydrostatic limit, the governing equations have periodic, nonlinear inertia–gravity solutions that are explored as models of the nonlinear internal tide. These long waves are shown to be robust to weak nonhydrostatic effects. Numerical solutions show that the disintegration of an initial sinusoidal linear internal tide is closely linked to the presence of these nonlinear waves. The initial tide steepens due to nonlinearity and sheds energy into short solitary waves. The disintegration is halted as the longwave part of the solution settles onto a state close to one of the nonlinear hydrostatic solutions, with the short solitary waves superimposed. The degree of disintegration is a function of initial amplitude of the tide and the properties of the underlying nonlinear hydrostatic solutions, which, depending on stratification and tidal frequency, exist only for a finite range of amplitudes (or energies). There is a lower threshold below which no short solitary waves are produced. However, for initial amplitudes above another threshold, given approximately by the energy of the limiting nonlinear hydrostatic inertia–gravity wave, most of the initial tidal energy goes into solitary waves. Recent observations in the South China Sea are briefly discussed.

Description: KRH was supported by a Woods Hole Oceanographic Institution Mellon Independent Study Award and ONR Grant N000140610798.

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

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

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

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

Description: Oceanic profiles of temperature, salinity, horizontal velocity, rate of dissipation of turbulent kinetic energy (ε) and rate of dissipation of thermal variance (χ) are used to examine the parameterization of turbulent mixing in the ocean due to internal waves. Turbulent mixing is quantified through eddy diffusivity parameterizations of the mass (Kρ; Osborn, 1980) and heat fluxes (Kτ; Osborn and Cox, 1972) in turbulent production/dissipation balances. Turbulence in the ocean is generally held to result from the occurrence of shear instability in regions where the Richardson number is locally supercritical (i.e. Ri ≤ 1/4), permitting the growth of small-scale waves which break and result in turbulent mixing. The occurrence of shear instability results from the local intensification of the shear in the internal wave field. The energy dissipated in such events is provided by the energy flux to higher wavenumber due to nonlinear wave/wave interactions on scales of 10's to 100's of meters. In turn, the strength of the wave/wave interactions depends generally on the energy content of the internal wave field, which can vary considerably over even larger scales due to the presence of topography or background flows. The magnitude of turbulent mixing is linked to internal wave dynamics by equating the turbulent dissipation with the energy flux through the vertical wavenumber spectrum under the priviso that the model spectrum which forms the basis for the analysis is statistically stationary with respect to the nonlinear interactions. Dynamical models (McComas and Muller, 1981; Henyey et al., 1986) indicate that the Garrett and Munk (GM; Munk, 1981) spectrum is stationary. Observations from the far field of a seamount in a region of negligible large-scale flow were examined to address the issue of the buoyancy scaling of ε. These data exhibited large variations in background stratification with depth, but the internal wave characteristics were not substantially differentiable from the GM prescription. The magnitude of ε and its functional dependence upon internal wave energy levels (E) and buoyancy frequency (N) was best described by the dynamical model ofHenyey et al. (1986) (ε ~ E2N2). The Richardson number scaling model of Kunze et al. (1990) produced consistent estimates. A second dynamical model, McComas and Muller (1981), predicted an appropriate (E,N) scaling, but overestimated the observed dissipation rates by a factor of five. Two kinematical dissipation parameterizations (Garmett and Holloway (1984) and Munk (1981)) predicted buoyancy scalings of N3/2 which were inconsistent with the observed scaling. Data from an upper-ocean front, a warm core ring and a region of steep topography were analyzed in order to examine the parameter dependence of E in internal wave fields which exhibited potentially nonstationary characteristics. Evidence was provided which implied the internal wave field in an upper ocean front was interacting with and modified by the background flow. Inhomogeneity and anisotropy of the internal wave field were noted in that data set. The model of Gregg (1989), which in turn was based upon the model of Henyey et al., effectively collapsed the observed diffusivity estimates from the front. The warm core ring profiles were noted to be anisotropic, dominated by near-inertial frequencies and to have a peaked vertical wavenumber shear spectrum. The data from a region of steep topography were noted to have a peaked vertical wavenumber spectrum and were characterized by higher than GM frequency motions. For the latter two data sets, application of a frequency based correction to the Henyey et al. model (Henyey, 1991) reduced more than an order of magnitude scatter in the parameterized estimates of E to less than a factor of four. Of the possible non-equilibrium conditions in the internal wave field, the (E,N) scaled dissipation rates were most sensitive to deviations in wave field frequency content. On the basis of a number of theoretical Richardson number probability distributions (Ri = N2/S2, where S2 is the sum of the squared vertical derivatives of horizontal velocity), the nominal dissipation scaling of the Kunze et al. model was determined to be E2N3. This scaling is altered to the observed ε ~ E2N2 scaling by a statistical dependence between N2 and S2 which reduces the occurrence of supercritical Ri values. This statistical dependence is hypothesized to be an effect of the turbulent momentum and buoyancy fluxes on the internal wave shear and strain profiles caused by shear instability. The statistical dependence between N2 and S2 exhibited a buoyancy scaling which was interpreted as resulting from the decreasing ratio between the time scale of the shear instability mechanism [T- 2π/N] and the adiabatic time scale [T - 2π/(Nf)1/2] of the internal wave field (f is the Coriolis parameter). This phenomenology is interpreted in light of saturated spectral theories which suggest that the magnitude and shape of the vertical wavenumber spectrum is controlled by instability mechanisms at large wavenumber ( ≥ .1 cpm). We argue that saturated spectral theories are valid only in the limit where a separation exists between the two time scales, i.e. for large N, low internal wave frequency content, and small f. These results have immediate implications for oceanic mixing driven by internal wave motions. First, background diffusivities are small: at GM energy levels, Kρ - .03x10-4 m2/s (Kρ = .25ε/N2). Secondly, since Kρ is independent of N at constant E, some process or collection of processes must be responsible for heightened E values in the abyss if internal waves cause the 0(1-10x10-4 m2/s) diffusivities generally inferred from deep ocean hydrographic data. We view internal wave reflection and/or internal wave generation associated with topographic features to be likely candidates.

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 1989

Description: Given well known environmental conditions, matched field processing has been shown to be a promising signal processing technique for the localization of acoustic sources. However, when environmental data are incomplete or inaccurate, a 'mismatch' occurs between the measured field and model field which can lead to a severe degradation of the localization estimator. We investigate the possible mismatch effects of surface and internal waves on matched field processing in a shallow water waveguide. We utilize a modified ray theory, based on the work of Tindle, to calculate the acoustic pressure field. This allows us to simply incorporate range dependent environmental conditions as well as to generalize our work to deeper waveguides. In general, the conventional (Bartlett) matched field beamformer does not provide sufficient resolution to unambiguously locate a source, even in a perfectly matched environment. The maximum likelihood method (MLM) matched field beamformer has much better resolution but is extremely susceptible to mismatch. The mismatch due to surface roughness can result in a large reduction of the estimator peak. Part, but not all, of the peak can be regained by 1)using a model which includes incomplete reflection at the surface based on actual sea surface statistics and 2) short time averaging of the measured signal, with times on the order of the period of the surface waves. Mismatch due to internal waves can also result in a large degradation of the estimator. Averaging over the same time period as surface waves provides little improvement and leads one to surmise that internal waves may be a limiting constraint on matched field processing. Finally, we combine the surface and internal wave fields with a slowly moving source. This example highlights the necessity for the development of a beamformer which has a broader mainlobe while maintaining adequate sidelobe suppression, and we address this issue by looking at two such beamformers.

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): 1981–2000, doi:10.1175/JPO-D-12-028.1.

Description: Packets of nonlinear internal waves (NLIWs) in a small area of the Mid-Atlantic Bight were 10 times more energetic during a local neap tide than during the preceding spring tide. This counterintuitive result cannot be explained if the waves are generated near the shelf break by the local barotropic tide since changes in shelfbreak stratification explain only a small fraction of the variability in barotropic to baroclinic conversion. Instead, this study suggests that the occurrence of strong NLIWs was caused by the shoaling of distantly generated internal tides with amplitudes that are uncorrelated with the local spring-neap cycle. An extensive set of moored observations show that NLIWs are correlated with the internal tide but uncorrelated with barotropic tide. Using harmonic analysis of a 40-day record, this study associates steady-phase motions at the shelf break with waves generated by the local barotropic tide and variable-phase motions with the shoaling of distantly generated internal tides. The dual sources of internal tide energy (local or remote) mean that shelf internal tides and NLIWs will be predictable with a local model only if the locally generated internal tides are significantly stronger than shoaling internal tides. Since the depth-integrated internal tide energy in the open ocean can greatly exceed that on the shelf, it is likely that shoaling internal tides control the energetics on shelves that are directly exposed to the open ocean.

Description: This research was supported by ONR Grants N00014-05-1-0271, N00014-08-1-0991, N00014-04- 1-0146, and N00014-11-1-0194.

Description: 2013-05-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 August 1996

Description: A field experiment was undertaken during July and August of 1995 aimed at understanding the interaction of acoustic signals with the internal wave field off the coast of New Jersey. As part of SWARM (Shallow Water Acoustics in a Random Medium), physical data were collected in 75 m of water near 39°15.34'N, 72°56.59'W with three thermistor strings, a bottom-mounted ADCP, and yo-yo CTDs. These data spanned a two-week period of the month-long study. With the exception of a time following a storm event, during which the generation mechanism near the shelf break was effectively switched off, large-amplitude (up to 20 meters), rank-ordered groups of internal solitons were observed traveling through the region approximately every 12.4 hours. These groups of solitons progressed across the shelf with phase speeds of 61.8 ± 14.9 cm/s with a heading of 280 ± 31° T. Two-layer finite-depth theory was tested on this data and shown to consistently overpredict the phase speed of the internal solitons within each group. Predictions of horizontal scale, particle velocities, and displacements were in qualitative agreement with two-layer finite-depth dynamics.

Description: Support for this work was provided by a National Science Foundation grant OCE-9313670.

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: 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 scattering of low-mode internal tides by ocean-floor topography is extensively studied through analytical models and field observations at the Line Islands Ridge (LIR). An existing Green function method is utilized to examine the generation of internal tides by idealized topographic shapes as well as realistic transects of the LIR. The method is also applied to examine the scattering of a mode-1 internal tide at these topographies to determine the relative high mode energy flux due to generated and scattered internal tides at the realistic transects. A method of determining the modal content of an internal wave field is advanced to account for arbitrary stratification and rotation. It is then adjusted to allow for image loss as is common to oceanographic studies. Its performance is compared to the existing regression method widely used by oceanographers to determine the modal content of internal tides. The results from this comparison are used to inform the analysis of the field observations. This thesis concludes by examining the modal content of the LIR as determined from measurements taken during the 150-day EXperiment on Internal Tide Scattering (EXITS) NSF field study. Motivated by satellite altimetry data and three-dimensional numerical model studies, the EXITS cruise sought to observe the internal tide scattering process in the ocean for the first time. The data from three moorings equipped with moored profilers, spanning total depths of 3000-5000 m is analyzed to determine the modal content of the southward propagating internal tide before and after it encounters the ridge for evidence of topographic scattering.

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

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

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

Description: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 1831-1848, doi:10.1175/JPO-D-18-0068.1.

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

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

Description: 2019-02-15

Description: Part I: It is to be expected that in the future, measurements of pressure fluctuations beneath the ocean surface will provide basic data for solution of many practical problems on the state of the sea. As a result of the wartime impetus, several types of underwater pressure recorders (of similar instrumentation principles) were developed both in this country and abroad. Essentially, the instrument consists of an underwater unit which electrically transmits pressure impulses near the sea bottom to a clockwork recorder installed on the shore. The underwater pressure unit is adjusted for pressure fluctuations resulting from surface waves within the spectrum band of periods set up by winds acting on the sea surface. The resulting records may be scaled for height and period of the pressure fluctuations over known time intervals. An accessory wave analyzer has been constructed for rapid periodogram analysis of the pressure records. This investigation is concerned with a comparative study of observed sea surface waves and recorded sea bottom pressure fluctuations. It was undertaken for the purpose of evaluating sea surface wave heights from sea bottom pressure recordings in the vicnity of Woods Hole and Bermuda. Part II: The results presented here comprise data on the state of the sea surface in the vicinity of Woods Hole (off Cuttyhunk Island, Massachusetts) July 1946 to May 1947, and off Bermuda, British West Indies, February 1947 to May 1947, obtained from analyses of automatic wave recordings. At both locations measuring elements of the wave recording instruments were located on the sea bottom (approximate depth 75 feet at Cuttyhunk and 120 feet at Bermuda) and electrically connected to shore recorders several miles away. The instrument operational schedules were continuous twenty minute recordings every two hours at Bermuda and every six hours at Cuttyhunk.

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

Description: The structure of the inertial peak in deep ocean kinetic energy spectra is studied here. Records were obtained from Polymode arrays deployed in the Western North Atlantic Ocean (40°W to 70°W, 15°N to 42°N). The results are interpreted both in terms of local sources and of turning point effects on internal waves generated at lower latitudes. In most of the data, there is a prominent inertial peak slightly above f; however, the peak height above the background continuum varies with depth and geographical environment. Three classes of environment and their corresponding spectra emerge from peak height variations: class 1 is the 1500 m level near the Mid-Atlantic Ridge, with the greatest peak height of 18 db; class 2 includes (a) the upper ocean (depth less than 2000 m), (b) the deep ocean (depth greater than 2000 m) over rough topography, and (c) the deep ocean underneath the Gulf Stream, with intermediate peak height of 11.5 db; class 3 is the deep ocean over smooth topography, with the lowest peak height of 7.5 db. Near f, the horizontal coherence scale is 0(60 km) at depths from 200 m to 600 m, and the vertical coherence scale is O(200 m) just below the main thermocline. A one turning point model is developed to describe inertial waves at mid-latitudes, based on the assumption that inertial waves are randomly generated at lower latitudes (global generation) where their frequency-wavenumber spectrum is given by the model of Garrett and Munk (1972 a, 1975). Using the globally valid wave functions obtained by Munk and Phillips (1968), various frequency spectra near f are calculated numerically. The model yields a prominent inertial peak of 7 db in the horizontal velocity spectrum but no peaks in the temperature spectrum. The model is latitudinally dependent: the frequency shift and bandwidth of the inertial peak decrease with latitude; energy level near f is minimum at about 30° and higher at low and high latitudes. The observations of class 3 can be well-described by the model; a low zonal wavenumber cutoff is required to produce the observed frequency shift of the inertial peak. The differences between the global generation model and the observations of class 1 and class 2 are interpreted as the effects of local sources. A locally forced model is developed based on the latitudinal modal decomposition of a localized source function. Asymptotic eigensolutions of the Laplace's tidal equation are therefore derived and used as a set of expansion functions. The forcing is through a vertical velocity field specified at the top or bottom boundaries of the ocean. For white noise forcing, the horizontal velocity spectrum of the response has an inertial peak which diminishes in the far-field. With the forcing located at either the surface or the bottom, several properties of the class 2 observations can be described qualitatively by a combination of the global and local models. The reflection of inertial waves from a turbulent benthic boundary layer is studied by a slab model of given depth. Frictional effects are confined to the boundary layer and modelled by a quadratic drag law. For given incident waves, reflection coefficients are found to be greater than 0.9 for the long waves which contain most of the energy. This result suggests that energy-containing inertial waves can propagate over great distance as is required by the validity of the model of global generation.

Description: This work was supported by the National Science Foundation through grant OCE 76-80210 and its continuation OEE 78-19833.

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 1982

Description: The development of nonlinear surface and internal wave groups is investigated. Surface wave evolution was observed in an unusually long wave channel as a function of steepness and group length. Dissipation and frequency downshifting were important characteristics of the long-time evolution. The amplitude and phase modulations were obtained using the Hilbert transform and specified as an initial condition to the cubic nonlinear Schrodinger equation, which was solved numerically. This equation is known to govern the slowly varying complex modulation envelope of gravity waves on deep water. When dissipation was included, the model compared quite well with the observations. Phase modulation was used to interpret the long-time behavior, using the phase evolution of exact asymptotic solutions as a guide. The wave groups exhibited a long-time coherence but not the recurrence predicted by the inviscid theory. An oceanic field study of the generation of groups of large amplitude internal waves by stratified tidal flow over a submarine ridge indicates that the large amplitude and asymmetry of the topography are critical in determining the type of flow response. The calculated Froude numbers response length scale and duration differ markedly between the two phases of the tide due to the asymmetry.

Description: Research assistantship provided by the Office of Naval Research contract no. N00014-80-C-0273

Description: Author Posting. © American Meteorological Society, 2009. 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 39 (2009): 1035-1049, doi:10.1175/2008JPO3920.1.

Description: Seasonal variability of near-inertial horizontal kinetic energy is examined using observations from a series of McLane Moored Profiler moorings located at 39°N, 69°W in the western North Atlantic Ocean in combination with a one-dimensional, depth-integrated kinetic energy model. The time-mean kinetic energy and shear vertical wavenumber spectra of the high-frequency motions at the mooring site are in reasonable agreement with the Garrett–Munk internal wave description. Time series of depth-dependent and depth-integrated near-inertial kinetic energy are calculated from available mooring data after filtering to isolate near-inertial-frequency motions. These data document a pronounced seasonal cycle featuring a wintertime maximum in the depth-integrated near-inertial kinetic energy deriving chiefly from the variability in the upper 500 m of the water column. The seasonal signal in the near-inertial kinetic energy is most prominent for motions with vertical wavelengths greater than 100 m but observable wintertime enhancement is seen down to wavelengths of the order of 10 m. Rotary vertical wavenumber spectra exhibit a dominance of clockwise-with-depth energy, indicative of downward energy propagation and implying a surface energy source. A simple depth-integrated near-inertial kinetic energy model consisting of a wind forcing term and a dissipation term captures the order of magnitude of the observed near-inertial kinetic energy as well as its seasonal cycle.

Description: Funding to initiate the McLane Moored Profiler observations at Line W were provided by grants from the G. Unger Vetlesen Foundation and the Comer Charitable Fund to the Woods Hole Oceanographic Institution’s Ocean and Climate Change Institute. Ongoing moored observations at Line W are supported by the National Science Foundation (NSF Grant OCE-0241354).

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

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

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

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

Description: This thesis describes the physics of fully three-dimensional low frequency acoustic interaction with internal waves, bottom sediment waves and surface swell waves that are often observed in shallow waters and on continental slopes. A simple idealized model of the ocean waveguide is used to analytically study the properties of acoustic normal modes and their perturbations due to waves of each type. The combined approach of a semi-quantitative study based on the geometrical acoustics approximation and on fully three-dimensional coupled mode numerical modeling is used to examine the azimuthal dependence of sound wave horizontal reflection from, transmission through and ducting between straight parallel waves of each type. The impact of the natural crossings of nonlinear internal waves on horizontally ducted sound energy is studied theoretically and modeled numerically using a three-dimensional parabolic equation acoustic propagation code. A realistic sea surface elevation is synthesized from the directional spectrum of long swells and used for three-dimensional numerical modeling of acoustic propagation. As a result, considerable normal mode amplitude scintillations were observed and shown to be strongly dependent on horizontal azimuth, range and mode number. Full field numerical modeling of low frequency sound propagation through large sand waves located on a sloped bottom was performed using the high resolution bathymetry of the mouth of San Francisco Bay. Very strong acoustic ducting is shown to steer acoustic energy beams along the sand wave’s curved crests.

Description: Office of Naval Research for the financial support of this work.

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

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

Description: 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): 259–282, doi:10.1175/JPO-D-11-0194.1.

Description: This study reports on observations of turbulent dissipation and internal wave-scale flow properties in a standing meander of the Antarctic Circumpolar Current (ACC) north of the Kerguelen Plateau. The authors characterize the intensity and spatial distribution of the observed turbulent dissipation and the derived turbulent mixing, and consider underpinning mechanisms in the context of the internal wave field and the processes governing the waves’ generation and evolution. The turbulent dissipation rate and the derived diapycnal diffusivity are highly variable with systematic depth dependence. The dissipation rate is generally enhanced in the upper 1000–1500 m of the water column, and both the dissipation rate and diapycnal diffusivity are enhanced in some places near the seafloor, commonly in regions of rough topography and in the vicinity of strong bottom flows associated with the ACC jets. Turbulent dissipation is high in regions where internal wave energy is high, consistent with the idea that interior dissipation is related to a breaking internal wave field. Elevated turbulence occurs in association with downward-propagating near-inertial waves within 1–2 km of the surface, as well as with upward-propagating, relatively high-frequency waves within 1–2 km of the seafloor. While an interpretation of these near-bottom waves as lee waves generated by ACC jets flowing over small-scale topographic roughness is supported by the qualitative match between the spatial patterns in predicted lee wave radiation and observed near-bottom dissipation, the observed dissipation is found to be only a small percentage of the energy flux predicted by theory. The mismatch suggests an alternative fate to local dissipation for a significant fraction of the radiated energy.

Description: SW acknowledges the support of the Grantham Institute for Climate Change, Imperial College London. ACNG acknowledges the support of a NERC Advanced Research Fellowship (Grant NE/C517633/1). KLP acknowledges support from Woods Hole Oceanographic Institution bridge support funds.

Description: 2013-08-01

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

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

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

Description: 2013-10-01

Description: 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): 222–230, doi:10.1175/JPO-D-12-099.1.

Description: Observations with fine horizontal resolution are used to identify the horizontal scales of variability over the Middle Atlantic Bight (MAB) shelf break and continental rise. Spray gliders collected observations along two alongshelf transects over the continental rise in March–April 2006 and along 16 cross-shelf transects over the shelf break and continental rise during July–October 2007. Horizontal resolution varied from 1 km or finer over the shelf to 6 km in deep water. These observations allow horizontal thermohaline variability offshore of the MAB shelf break to be examined for the first time. Structure functions of temperature and salinity, the mean square difference between observations separated by specified distances, reveal the horizontal spatial scales in the region. Exponential (e-folding) scales of temperature and salinity increase from 8–13 km near the shelf break to about 30 km over the continental rise. Just offshore of the shelf break, alongshelf structure functions exhibit periodicity with a 40–50-km wavelength that matches the wavelength of shelfbreak frontal meanders. Farther offshore, alongshelf structure functions suggest a dominant wavelength of 175–250 km, but these scales are only marginally resolved by the available observations. Examination of structure functions of along-isopycnal salinity (i.e., spice) suggests that interleaving of shelf and slope water masses contributes most of the horizontal variability near the MAB shelf break, but heaving of isopycnals is the primary source of horizontal variability over the continental rise.

Description: 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. GGG was supported by the National Science Foundation under Grant OCE-1129125.

Description: 2013-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): 427–444, doi:10.1175/JPO-D-13-070.1.

Description: Between 25 September 2007 and 28 September 2009, a heavily instrumented mooring was deployed in the Labrador Sea, offshore of the location where warm-core, anticyclonic Irminger rings are formed. The 2-year time series offers insight into the vertical and horizontal structure of newly formed Irminger rings and their heat and salt transport into the interior basin. In 2 years, 12 Irminger rings passed by the mooring. Of these, 11 had distinct properties, while 1 anticyclone likely passed the mooring twice. Eddy radii (11–35 km) were estimated using the dynamic height signal of the anticyclones (8–18 cm) together with the observed velocities. The anticyclones show a seasonal cycle in core properties when observed (1.9°C in temperature and 0.07 in salinity at middepth) that has not been described before. The temperature and salinity are highest in fall and lowest in spring. Cold, fresh caps, suggested to be an important source of freshwater, were seen in spring but were almost nonexistent in fall. The heat and freshwater contributions by the Irminger rings show a large spread (from 12 to 108 MJ m−2 and from −0.5 to −4.7 cm, respectively) for two reasons. First, the large range of radii leads to large differences in transported volume. Second, the seasonal cycle leads to changes in heat and salt content per unit volume. This implies that estimates of heat and freshwater transport by eddies should take the distribution of eddy properties into account in order to accurately assess their contribution to the restratification.

Description: This work was supported by the U.S. National Science Foundation and the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Devonshire Foundation.

Description: 2014-08-01

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

Description: Potential vorticity structure in two segments of the North Atlantic’s western boundary current is examined using concurrent, high-resolution measurements of hydrography and velocity from gliders. Spray gliders occupied 40 transects across the Loop Current in the Gulf of Mexico and 11 transects across the Gulf Stream downstream of Cape Hatteras. Cross-stream distributions of the Ertel potential vorticity and its components are calculated for each transect under the assumptions that all flow is in the direction of measured vertically averaged currents and that the flow is geostrophic. Mean cross-stream distributions of hydrographic properties, potential vorticity, and alongstream velocity are calculated for both the Loop Current and the detached Gulf Stream in both depth and density coordinates. Differences between these mean transects highlight the downstream changes in western boundary current structure. As the current increases its transport downstream, upper-layer potential vorticity is generally reduced because of the combined effects of increased anticyclonic relative vorticity, reduced stratification, and increased cross-stream density gradients. The only exception is within the 20-km-wide cyclonic flank of the Gulf Stream, where intense cyclonic relative vorticity results in more positive potential vorticity than in the Loop Current. Cross-stream gradients of mean potential vorticity satisfy necessary conditions for both barotropic and baroclinic instability within the western boundary current. Instances of very low or negative potential vorticity, which predispose the flow to various overturning instabilities, are observed in individual transects across both the Loop Current and the Gulf Stream.

Description: Glider operations in the Gulf Stream were supported by the National Science Foundation under Grant OCE-0220769. Glider operations in the Gulf of Mexico were supported by BP. R.E.T. was supported by the Penzance Endowed Fund in Support of Assistant Scientists and the Independent Research and Development Program at WHOI.

Description: 2016-07-01

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

Description: Measurements of horizontal and vertical current by propeller cluster current meters and temperature by thermistors mounted on a rigid array 8 m high and 20 m long moored in the oceanic main thermocline near Bermuda are interpreted in terms of thermocline-trapped internal wave modes in the presence of temperature and density fine-structure. Two turning-point uniformly valid asymptotic solutions to the internal wave equation are developed to describe the wave functions. Mode decay beyond the turning point in depth or frequency produces a sharp cutoff in vertical current spectra above the local buoyancy frequency N(z). An internal wave wavenumber-frequency spectral model Ε(α,ω) = E(ω/No)-2 (α./α0)-2 describes vertical current spectra and potential energy to horizontal kinetic energy ratios. The red wavenumer shape suppresses peaks in both these quantities at frequencies near N(z). The data are consistent with time-averaged horizontal isotropy of the wave field. A dip in the vertical current spectra at 0.5 cph not predicted by the model appears related to the bottom slope. Temperature fine-structure is modeled as a passive vertical field advected by internal waves. Quasi-permanent fine-scale features of the stratification and vertically small-scale internal waves are indistinguishable in this study. The model of McKean (1974) is generalized to include fine-structure fields specified by their vertical wavenumber spectra as well as different Poisson-distributed layer models. Together with the trapped internal wave model, moored temperature spectra, temperature vertical difference spectra, and coherence over vertical separations are described using a fine-structure vertical wavenumber spectrum PT(k) =ATk-5/2 which agrees with other spectra made using vertical profiling instruments in the range 0.1 to 1.0 cpm. Horizontal current fine-structure is also modeled as a passive field advected vertically by long internal waves. The model describes moored horizontal current spectra (least successfully at frequencies near N(z)) and finite-difference vertical shear spectra. Contours of temperature in depth versus time indicate possible mixing events. These events appear concurrently with high shear and Richardson numbers O. 25≤ R ≤ 1.0. Over 7 m a cutoff in Ri at 0.25 is observed, indicating saturation of the internal wave spectrum. Spectra of finite-difference approximations to shear and buoyancy frequency are dominated by fine-structure contributions over nearly the whole internal wave range, suggesting that breaking is enhanced by fine-structure. Breaking appears equally likely at all frequencies in the internal wave range.

Description: This research was supported by Office of Naval Research contract N00014-67-0204-0047 and continuation contract NOOOl4-75-C-0291.

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 1982

Description: Velocity and temperature time series from Hudson Submarine Canyon and hydrographic surveys of seven canyons of the Middle Atlantic Bight indicate that the effects of storms, tides, and incoming internal waves are intensified in submarine canyons. Storms with strong eastward and westward wind stress were found to cause strong upwelling and downwelling through the upper layers of Hudson Canyon. Storm-forced upwelling also caused strong down-canyon flows at the canyon floor. Internal waves were found to be concentrated in the canyon head and near the floor, in agreement with theoretical predictions. Slope water apparently circulates slowly through the outer part of the canyon and is mixed in near-floor layers which could be caused by breaking internal waves. Internal tides are generated at the floor in the central part of the canyon. Oscillations at tidal frequencies dominate the near-floor velocity field below the thermocline, and are accompanied by high-frequency spikes that may be nonlinear interface waves propagating on the top of the bottom mixed layer. A numerical model was used to calculate mixing in the canyon's bottom boundary layer caused by an unstable density gradient during flood tide. Energetic internal wave activity is apparently responsible for sediment sorting in the canyon head; the internal waves become more energetic as the sediment grain size increases. Below the thermocline, the tidal oscillations vary in amplitude with the phases of the moon; the observed deposition of mud can easily occur during weeks of low velocity.

Description: Foundation graduate fellowship and by the Office of Naval Research under Contracts N00014-75-C-029l and N00014-80-C-0273.

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

Description: Garrett and Munk use linear dynamics to synthesize frequency-wavenumber energy spectra for internal waves (GM72, GM75, GM79). The GM internal wave models are horizontally isotropic, vertically symmetric, purely propagating, and universal in both time and space. This set of properties effectively eliminates all the interesting physics, since such models do not allow localized sources and sinks of energy. Thus an important step in understanding internal wave dynamics is to make measurements of deviations from the simple GM models. This thesis continues the search for deviations from the GM models. It has three advantages over earlier work: extensive data from an equatorial region, long time series (2 years), and relatively sophisticated linear internal wave models. Since the GM models are based on mid-latitude data, having data from an equatorial region which has a strong mean current system offers an opportunity to examine a region with a distinctly different basic state. The longer time series mean there is a larger statistical ensemble of realizations, making it possible to detect smaller internal wave signals. The internal wave models include several important extensions to the GM models: horizontal anisotropy and vertical asymmetry, resolution between standing modes and propagating waves, general vertical structure, and kinematic effects of mean shear flow. Also investigated are the effects of scattering on internal waves, effects that are especially strong on the equator because the buoyancy frequency variability is a factor of ten higher than at mid-latitudes. In the high frequency internal wave field considered (frequencies between .125 cph and .458 cph), several features are found that are not included in the GM models. Both the kinematic effects of a mean shear flow and the phase-locking that distinguishes standing modes from propagating waves are observed. There is a seasonal dependence in energy level of roughly 10% of the mean level. At times the wave field is zonally and vertically asymmetric, with resulting energy fluxes that are a small (4% to 10%) fraction of the maximum energy flux the internal wave field could support. The fluxes are, however, as big as many of the postulated sources of energy for the internal wave field.

Description: This work has been supported under grants from the National Science Foundation and the Office of Naval Research, grants numbered NSF-89076, ONR-88914, NSF-9l002, NSF-94971, and NSF-93661.

Description: Author Posting. © American Meteorological Society, 2009. 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 39 (2009): 3162-3175, doi:10.1175/2009JPO4239.1.

Description: This study analyzes anisotropic properties of the material transport by eddies and eddy-driven zonal jets in a general circulation model of the North Atlantic through the analysis of Lagrangian particle trajectories. Spreading rates—defined here as half the rate of change in the particle dispersion—in the zonal direction systematically exceed the meridional rates by an order of magnitude. Area-averaged values for the upper-ocean zonal and meridional spreading rates are approximately 8100 and 1400 m2 s−1, respectively, and in the deep ocean they are 2400 and 200 m2 s−1. The results demonstrate that this anisotropy is mainly due to the action of the transient eddies and not to the shear dispersion associated with the time-mean jets. This property is consistent with the fact that eddies in this study have zonally elongated shapes. With the exception of the upper-ocean subpolar gyre, eddies also cause the superdiffusive zonal spreading, significant variations in the spreading rate in the vertical and meridional directions, and the difference between the westward and eastward spreading.

Description: Funding for IK was provided by NSF Grants OCE 0346178, 0749722, and 0842834. Funding for PB was provided by NSF Grants OCE 0344094 and OCE 0725796 and by the research grant from the Newton Trust of the University of Cambridge. For JP the acknowledgement is to NSF OCE-0451086.

Description: Author Posting. © American Meteorological Society, 2009. 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 39 (2009): 1361-1379, doi:10.1175/2008JPO4096.1.

Description: Multiple zonal jets are observed in satellite data–based estimates of oceanic velocities, float measurements, and high-resolution numerical simulations of the ocean circulation. This study makes a step toward understanding the dynamics of these jets in the real ocean by analyzing the vertical structure and dynamical balances within multiple zonal jets simulated in an eddy-resolving primitive equation model of the North Atlantic. In particular, the authors focus on the role of eddy flux convergences (“eddy forcing”) in supporting the buoyancy and relative/potential vorticity (PV) anomalies associated with the jets. The results suggest a central role of baroclinic eddies in the barotropic and baroclinic dynamics of the jets, and significant differences in the effects of eddy forcing between the subtropical and subpolar gyres. Additionally, diabatic potential vorticity sources and sinks, associated with vertical diffusion, are shown to play an important role in supporting the potential vorticity anomalies. The resulting potential vorticity profile does not resemble a “PV staircase”—a distinct meridional structure observed in some idealized studies of geostrophic turbulence.

Description: Funding for IK was provided by NSF Grants OCE 0346178 and 0749722. Funding for PB was provided by NSF Grants OCE 0344094 and OCE 0725796 and by the research grant from the Newton Trust of the University of Cambridge. For JP the acknowledgement is to NSF OCE-0451086.

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): 2307–2327, doi:10.1175/JPO-D-10-05004.1.

Description: Results from a high-resolution (~2 km) numerical simulation of the Irminger Basin during summer 2003 are presented. The focus is on the East Greenland Spill Jet, a recently discovered component of the circulation in the basin. The simulation compares well with observations of surface fields, the Denmark Strait overflow (DSO), and the hydrographic structure of typical sections in the basin. The model reveals new aspects of the circulation on scales of O(0.1–10) days and O(1–100) km. The model Spill Jet results from the cascade of dense waters over the East Greenland shelf. Spilling can occur in various locations southwest of the strait, and it is present throughout the simulation but exhibits large variations on periods of O(0.1–10) days. The Spill Jet sometimes cannot be distinguished in the velocity field from surface eddies or from the DSO. The vorticity structure of the jet confirms its unstable nature with peak relative and tilting vorticity terms reaching twice the planetary vorticity term. The average model Spill Jet transport is 4.9 ±1.7 Sv (1 Sv ≡ 106 m3 s−1) equatorward, about 2½ times larger than has been previously reported from a single ship transect in August 2001. Kinematic analysis of the model results suggests two different types of spilling events. In the first case (type I), a local perturbation results in dense waters descending over the shelf break into the Irminger Basin. In the second case (type II), surface cyclones associated with DSO deep domes initiate the spilling process. During summer 2003, more than half of the largest Spill Jet transport values are of type II.

Description: The research is supported by the National Science Foundation Grants OCE-0726393 and OCI-0904640 (MGM and TWNH) and OCE-0726640 (RSP).

Description: 2012-06-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): 2223–2241, doi:10.1175/2011JPO4344.1.

Description: Results are presented from an observational study of stratified, turbulent flow in the bottom boundary layer on the outer southeast Florida shelf. Measurements of momentum and heat fluxes were made using an array of acoustic Doppler velocimeters and fast-response temperature sensors in the bottom 3 m over a rough reef slope. Direct estimates of flux Richardson number Rf confirm previous laboratory, numerical, and observational work, which find mixing efficiency not to be a constant but rather to vary with Frt, Reb, and Rig. These results depart from previous observations in that the highest levels of mixing efficiency occur for Frt 〈 1, suggesting that efficient mixing can also happen in regions of buoyancy-controlled turbulence. Generally, the authors find that turbulence in the reef bottom boundary layer is highly variable in time and modified by near-bed flow, shear, and stratification driven by shoaling internal waves.

Description: Funding was provided by grants from the National Oceanic and Atmospheric Administration’s National Undersea Research Program, National Science Foundation Grants OCE-0622967 and OCE- 0824972 to SGM, and the Singapore Stanford Program. Kristen Davis was supported by a National Defense Science and Engineering Graduate Fellowship and an ARCS Foundation Fellowship.

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 1995

Description: This thesis investigates the amplitude fluctuation effects in acoustic scattering due to shallow water internal waves. Theoretically, it uses the adiabatic approximation and perturbation methods to statistically evaluate acoustic transmission fluctuations caused by internal waves in the ocean; it also investigates acoustic mode coupling effect due to internal waves. Numerically, this thesis simulates the shallow water internal wave (IW) field using the Garrett-Munk internal wave spectrum model and then evaluates acoustic transmission in the simulated internal wave field with the Kraken normal mode program. Theoretical calculations are also performed using the theory developed in this thesis. Comparisons are made between theory and numerical calculations. The results presented and discussed in this thesis are related to the following issues: coherent and incoherent intensity fluctuations for adiabatic approximation, acoustic mode coupling due to IW's, transmission loss difference between adiabatic and coupled mode methods, and their dependence on range, IW amplitude and frequency.

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

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

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

Description: 2015-05-01

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

Description: Salinity and temperature profiles from drifting ice-tethered profilers in the Beaufort gyre region of the Canada Basin are used to characterize and quantify the regional near-inertial internal wave field over one year. Vertical displacements of potential density surfaces from the surface to 750-m depth are tracked from fall 2006 to fall 2007. Because of the time resolution and irregular sampling of the ice-tethered profilers, near-inertial frequency signals are marginally resolved. Complex demodulation is used to determine variations with a time scale of several days in the amplitude and phase of waves at a specified near-inertial frequency. Characteristics and variability of the wave field over the course of the year are investigated quantitatively and related to changes in surface wind forcing and sea ice cover.

Description: The ITP program and J. Toole’s contributions were supported by the National Science Foundation Office of Polar Programs Arctic Observing Network. We acknowledge the support of the Office of Naval Research (Grant N00014-11-1-0454) for this study. Support for H. Dosser was also provided by the Natural Sciences and Engineering Research Council of Canada.

Description: 2014-08-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): 2381–2406, doi:10.1175/JPO-D-14-0086.1.

Description: While near-inertial waves are known to be generated by atmospheric storms, recent observations in the Kuroshio Front find intense near-inertial internal-wave shear along sloping isopycnals, even during calm weather. Recent literature suggests that spontaneous generation of near-inertial waves by frontal instabilities could represent a major sink for the subinertial quasigeostrophic circulation. An unforced three-dimensional 1-km-resolution model, initialized with the observed cross-Kuroshio structure, is used to explore this mechanism. After several weeks, the model exhibits growth of 10–100-km-scale frontal meanders, accompanied by O(10) mW m−2 spontaneous generation of near-inertial waves associated with readjustment of submesoscale fronts forced out of balance by mesoscale confluent flows. These waves have properties resembling those in the observations. However, they are reabsorbed into the model Kuroshio Front with no more than 15% dissipating or radiating away. Thus, spontaneous generation of near-inertial waves represents a redistribution of quasigeostrophic energy rather than a significant sink.

Description: “The Study of Kuroshio Ecosystem Dynamics for Sustainable Fisheries (SKED)” supported by MEXT, MIT-Hayashi Seed Fund, ONR (Awards N000140910196 and N000141210101), NSF (Award OCE 0928617, 0928138) for support.

Description: 2016-03-01

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): 1789-1797, doi:10.1175/JPO-D-16-0240.1.

Description: Internal solitary waves are commonly observed in the coastal ocean where they are known to contribute to mass transport and turbulent mixing. While these waves are often generated by cross-isobath barotropic tidal currents, novel observations are presented suggesting that internal solitary waves result from along-isobath tidal flows over channel-shoal bathymetry. Mooring and ship-based velocity, temperature, and salinity data were collected over a cross-channel section in a stratified estuary. The data show that Ekman forcing on along-channel tidal currents drives lateral circulation, which interacts with the stratified water over the deep channel to generate a supercritical mode-2 internal lee wave. This lee wave propagates onto the shallow shoal and evolves into a group of internal solitary waves of elevation due to nonlinear steepening. These observations highlight the potential importance of three-dimensionality on the conversion of tidal flow to internal waves in the rotating ocean.

Description: National Science Foundation (OCE-1061609)

Description: 2018-01-03

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): 1969-1993, doi:10.1175/JPO-D-18-0031.1.

Description: Upstream mean semidiurnal internal tidal energy flux has been found in the Gulf Stream in hydrodynamical model simulations of the Atlantic Ocean. A major source of the energy in the simulations is the south edge of Georges Bank, where strong and resonant Gulf of Maine tidal currents are found. An explanation of the flux pattern within the Gulf Stream is that internal wave modal rays can be strongly redirected by baroclinic currents and even trapped (ducted) by current jets that feature strong velocities above the thermocline that are directed counter to the modal wavenumber vector (i.e., when the waves travel upstream). This ducting behavior is analyzed and explained here with ray-based wave propagation studies for internal wave modes with anisotropic wavenumbers, as occur in mesoscale background flow fields. Two primary analysis tools are introduced and then used to analyze the strong refraction and ducting: the generalized Jones equation governing modal properties and ray equations that are suitable for studying waves with anisotropic wavenumbers.

Description: The Woods Hole research was supported by National Science Foundation Grant OCE-1060430 and by the Office of Naval Research Grants N00014-11-1-0701 and N00014-17-1-2624. The USM research was supported by ONR Grant N00014-15-1-2288 and National Science Foundation Grant OCE-1537449.

Description: 2019-02-28

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 33(4), (2020): 1535-1545, doi:10.1175/JCLI-D-19-0547.1.

Description: In a transient warming scenario, the North Atlantic is influenced by a complex pattern of surface buoyancy flux changes that ultimately weaken the Atlantic meridional overturning circulation (AMOC). Here we study the AMOC response in the CMIP5 experiment, using the near-geostrophic balance of the AMOC on interannual time scales to identify the role of temperature and salinity changes in altering the circulation. The thermal wind relationship is used to quantify changes in the zonal density gradients that control the strength of the flow. At 40°N, where the overturning cell is at its strongest, weakening of the AMOC is largely driven by warming between 1000- and 2000-m depth along the western margin. Despite significant subpolar surface freshening, salinity changes are small in the deep branch of the circulation. This is likely due to the influence of anomalously salty water in the subpolar intermediate layers, which is carried northward from the subtropics in the upper limb of the AMOC. In the upper 1000 m at 40°N, salty anomalies due to increased evaporation largely cancel the buoyancy increase due to warming. Therefore, in CMIP5, temperature dynamics are responsible for AMOC weakening, while freshwater forcing instead acts to strengthen the circulation in the net. These results indicate that past modeling studies of AMOC weakening, which rely on freshwater hosing in the subpolar gyre, may not be directly applicable to a more complex warming scenario.

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 John Marshall for helpful discussions on the driving mechanisms of the AMOC, and three anonymous reviewers whose comments greatly improved the manuscript. 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.

Description: 2020-07-20

Description: Author Posting. © American Meteorological Society, 2010. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 40 (2010): 2381-2400, doi:10.1175/2010JPO4403.1.

Description: Langmuir circulation (LC) is a turbulent upper-ocean process driven by wind and surface waves that contributes significantly to the transport of momentum, heat, and mass in the oceanic surface layer. The authors have previously performed a direct comparison of large-eddy simulations and observations of the upper-ocean response to a wind event with rapid mixed layer deepening. The evolution of simulated crosswind velocity variance and spatial scales, as well as mixed layer deepening, was only consistent with observations if LC effects are included in the model. Based on an analysis of these validated simulations, in this study the fundamental differences in mixing between purely shear-driven turbulence and turbulence with LC are identified. In the former case, turbulent kinetic energy (TKE) production due to shear instabilities is largest near the surface, gradually decreasing to zero near the base of the mixed layer. This stands in contrast to the LC case in which at middepth range TKE production can be dominated by Stokes drift shear. Furthermore, the Eulerian mean vertical shear peaks near the base of the mixed layer so that TKE production by mean shear flow is elevated there. LC transports horizontal momentum efficiently downward leading to an along-wind velocity jet below LC downwelling regions at the base of the mixed layer. Locally enhanced vertical shear instabilities as a result of this jet efficiently erode the thermocline. In turn, enhanced breaking internal waves inject cold deep water into the mixed layer, where LC currents transport temperature perturbation advectively. Thus, LC and locally generated shear instabilities work intimately together to facilitate strongly the mixed layer deepening process.

Description: This research was supported by the Office of Naval Research through Grants N00014-09-M-0112 (TK) and N00014-06-1-0178 (AP, JT). Author TK also received support from a Woods Hole Oceanographic Institution Cooperative Institute for Climate and Ocean Research Postdoctoral Scholarship.

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

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

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

Description: 2014-09-01

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

Description: Simultaneous full-depth microstructure measurements of turbulence and finestructure measurements of velocity and density are analyzed to investigate the relationship between turbulence and the internal wave field in the Antarctic Circumpolar Current. These data reveal a systematic near-bottom overprediction of the turbulent kinetic energy dissipation rate by finescale parameterization methods in select locations. Sites of near-bottom overprediction are typically characterized by large near-bottom flow speeds and elevated topographic roughness. Further, lower-than-average shear-to-strain ratios indicative of a less near-inertial wave field, rotary spectra suggesting a predominance of upward internal wave energy propagation, and enhanced narrowband variance at vertical wavelengths on the order of 100 m are found at these locations. Finally, finescale overprediction is typically associated with elevated Froude numbers based on the near-bottom shear of the background flow, and a background flow with a systematic backing tendency. Agreement of microstructure- and finestructure-based estimates within the expected uncertainty of the parameterization away from these special sites, the reproducibility of the overprediction signal across various parameterization implementations, and an absence of indications of atypical instrument noise at sites of parameterization overprediction, all suggest that physics not encapsulated by the parameterization play a role in the fate of bottom-generated waves at these locations. Several plausible underpinning mechanisms based on the limited available evidence are discussed that offer guidance for future studies.

Description: The SOFine project is funded by the United Kingdom’s Natural Environmental Research Council (NERC) (Grant NE/G001510/1). SW acknowledges the support of anARCDiscovery Early CareerResearchAward (Grant DE120102927), as well as the Grantham Institute for Climate Change, Imperial College London, and the ARC Centre of Excellence for Climate System Science (Grant CE110001028). ACNG acknowledges the support of a NERC Advanced Research Fellowship (Grant NE/C517633/1).KLP acknowledges support fromWoods Hole Oceanographic Institution bridge support funds.

Description: 2014-11-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): 1854–1872, doi:10.1175/JPO-D-13-0104.1.

Description: The authors present inferences of diapycnal diffusivity from a compilation of over 5200 microstructure profiles. As microstructure observations are sparse, these are supplemented with indirect measurements of mixing obtained from (i) Thorpe-scale overturns from moored profilers, a finescale parameterization applied to (ii) shipboard observations of upper-ocean shear, (iii) strain as measured by profiling floats, and (iv) shear and strain from full-depth lowered acoustic Doppler current profilers (LADCP) and CTD profiles. Vertical profiles of the turbulent dissipation rate are bottom enhanced over rough topography and abrupt, isolated ridges. The geography of depth-integrated dissipation rate shows spatial variability related to internal wave generation, suggesting one direct energy pathway to turbulence. The global-averaged diapycnal diffusivity below 1000-m depth is O(10−4) m2 s−1 and above 1000-m depth is O(10−5) m2 s−1. The compiled microstructure observations sample a wide range of internal wave power inputs and topographic roughness, providing a dataset with which to estimate a representative global-averaged dissipation rate and diffusivity. However, there is strong regional variability in the ratio between local internal wave generation and local dissipation. In some regions, the depth-integrated dissipation rate is comparable to the estimated power input into the local internal wave field. In a few cases, more internal wave power is dissipated than locally generated, suggesting remote internal wave sources. However, at most locations the total power lost through turbulent dissipation is less than the input into the local internal wave field. This suggests dissipation elsewhere, such as continental margins.

Description: This research was funded by the Climate Process Team (CPT) on internal wave–driven mixing throughNSF GrantOCE-0968721. GSC acknowledges support from NSF Grants OCE-0825266 (EXITS), OCE-1029483 (SPAM), and OCE-1029722 (MIXET). LDT and CBW acknowledge support from NSF Grant OCE-0927650. SWand ACNG acknowledge support from NERC Grant NE/G001510/1 (SOFine).

Description: 2015-01-01

Description: © The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Kwon, Y., Seo, H., Ummenhofer, C. C., & Joyce, T. M. Impact of multidecadal variability in Atlantic SST on winter atmospheric blocking. Journal of Climate, 33(3), (2020): 867-892, doi: 10.1175/JCLI-D-19-0324.1.

Description: Recent studies have suggested that coherent multidecadal variability exists between North Atlantic atmospheric blocking frequency and the Atlantic multidecadal variability (AMV). However, the role of AMV in modulating blocking variability on multidecadal times scales is not fully understood. This study examines this issue primarily using the NOAA Twentieth Century Reanalysis for 1901–2010. The second mode of the empirical orthogonal function for winter (December–March) atmospheric blocking variability in the North Atlantic exhibits oppositely signed anomalies of blocking frequency over Greenland and the Azores. Furthermore, its principal component time series shows a dominant multidecadal variability lagging AMV by several years. Composite analyses show that this lag is due to the slow evolution of the AMV sea surface temperature (SST) anomalies, which is likely driven by the ocean circulation. Following the warm phase of AMV, the warm SST anomalies emerge in the western subpolar gyre over 3–7 years. The ocean–atmosphere interaction over these 3–7-yr periods is characterized by the damping of the warm SST anomalies by the surface heat flux anomalies, which in turn reduce the overall meridional gradient of the air temperature and thus weaken the meridional transient eddy heat flux in the lower troposphere. The anomalous transient eddy forcing then shifts the eddy-driven jet equatorward, resulting in enhanced Rossby wave breaking and blocking on the northern flank of the jet over Greenland. The opposite is true with the AMV cold phases but with much shorter lags, as the evolution of SST anomalies differs in the warm and cold phases.

Description: We gratefully acknowledge support from the NSF Climate and Large-scale Dynamics Program (AGS-1355339) to Y-OK, HS, CCU, and TMJ, the NASA Physical Oceanography Program (NNX13AM59G) to Y-OK, HS, and TMJ, NOAA CPO Climate Variability and Predictability Program (NA13OAR4310139) and DOE CESD Regional and Global Model Analysis Program (DE-SC0019492) to Y-OK, and NSF Physical Oceanography Program (OCE-1419235) to HS. We are very grateful to the three anonymous reviewers and editor Dr. Mingfang Ting, for their thorough and insightful suggestions. The NOAA 20CR dataset was downloaded from the NOAA Earth System Research Laboratory Physical Science Division webpage (https://www.esrl.noaa.gov/psd/data/20thC_Rean/). Support for the 20CR Project version 2c dataset is provided by the U.S. Department of Energy, Office of Science Biological and Environmental Research (BER), and by the National Oceanic and Atmospheric Administration Climate Program Office. The HadISST dataset was downloaded from the U.K. Met Office Hadley Centre webpage (https://www.metoffice.gov.uk/hadobs/hadisst/). The ERA-20C dataset was downloaded from the ECMWF webpage (https://apps.ecmwf.int/datasets/data/era20c-daily/). The ERSST5 dataset was provided by the NOAA Earth System Research Laboratory Physical Science Division (https://www.esrl.noaa.gov/psd/data/gridded/data.noaa.ersst.v5.html).

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 51(1), (2021): 19-35, https://doi.org/10.1175/JPO-D-19-0233.1.

Description: In the Beaufort Sea in September of 2015, concurrent mooring and microstructure observations were used to assess dissipation rates in the vicinity of 72°35′N, 145°1′W. Microstructure measurements from a free-falling profiler survey showed very low [O(10−10) W kg−1] turbulent kinetic energy dissipation rates ε. A finescale parameterization based on both shear and strain measurements was applied to estimate the ratio of shear to strain Rω and ε at the mooring location, and a strain-based parameterization was applied to the microstructure survey (which occurred approximately 100 km away from the mooring site) for direct comparison with microstructure results. The finescale parameterization worked well, with discrepancies ranging from a factor of 1–2.5 depending on depth. The largest discrepancies occurred at depths with high shear. Mean Rω was 17, and Rω showed high variability with values ranging from 3 to 50 over 8 days. Observed ε was slightly elevated (factor of 2–3 compared with a later survey of 11 profiles taken over 3 h) from 25 to 125 m following a wind event which occurred at the beginning of the mooring deployment, reaching a maximum of ε= 6 × 10−10 W kg−1 at 30-m depth. Velocity signals associated with near-inertial waves (NIWs) were observed at depths greater than 200 m, where the Atlantic Water mass represents a reservoir of oceanic heat. However, no evidence of elevated ε or heat fluxes was observed in association with NIWs at these depths in either the microstructure survey or the finescale parameterization estimates.

Description: This work was supported by NSF Grants PLR 14-56705 and PLR-1303791 and by NSF Graduate Research Fellowship Grant DGE-1650112.

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, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 51(3), (2021): 955–973, https://doi.org/10.1175/JPO-D-20-0240.1.

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

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

Description: This paper contains twelve charts of the North Atlantic Ocean west of the 60th meridian which show the average temperature of the surface water for each month of the year. Two additional charts show respectively the maximum and minimum values reached by the average temperature for any month. The next chart presents the yearly range, that is, the difference between the average temperature of the warmest and coldest months. The last chart of the series shows the depth of the virtually isothermal water during the summer and winter months.

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

Description: This thesis reports on an investigation into the structure, energetics and propagation of tidal frequency internal waves. Data from Site D, near the New England continental slope, Muir Seamount northeast of Bermuda, and the Mid-Ocean Dynamics Experiment in the deep Sargasso Sea were used. Site D, in the near-field of a near-critical semidiurnal generation region, shows variable tidal currents and a marked surface intensification of M2 energy at the southern Site, related to the beam-like nature of the internal tide. The M2 tide dominates the semidiurnal band, with about 3 times more energy than at adjacent frequencies at 1/15 cpd separation. There is a significant phase locking between the M2 baroclinic currents and the equilibrium tide, and evidence for southward propagation of internal wave energy, suggesting generation at the slope to the north. The M2 baroclinic energy density is about 40% as great as the total barotropic energy density, but the internal tides have more horizontal kinetic energy. A seaward energy flux of .6 x 106 erg/s cm in the first three baroclinic M2 modes is much less than the .2 x 1010 erg/s cm shoreward energy flux in the surface tide. Difficul ties in interpreting the measurements are ascribed to the near-singular generation case. The MODE-l semidiurnal internal tides are also dominated by the M2 frequency, with a 3-fold energy increase over adjacent frequencies at 1/15 cpd separation. MODE-l is far from any major source of internal tides, but the measurements are much less variable than those from Site D. The extensive temperature measurements defining the MODE-l M2 internal tide are significantly coherent (phase locked) with the equilibrium tide, with about 80% of the coherent energy deriving from the first baroclinic mode, typical thermocline displacements being 3 m. A horizontal wavenumber spectrum estimate for the first mode M2 displacement fluctuations gives a peak at 160 km wavelength, in excellent agreement with the theoretical dispersion relation. The coherent first mode propagates on a bearing of 125°T, with a horizontal energy flux of .3 x 108 erg/s cm. Use of the weaker S2 internal tide and the dispersive nature of oceanic internal waves yields an estimate of 700 km to a common semidiurnal source region. The inferred range and bearing are consistent with generation at the Blake Escarpment and the continental slope to the northwest of the experiment. In one special case current and temperature measurements are combined in a local demonstration of the first mode M2 propagation, and the less extensive current data gives estimates of the barotropic tidal currents. Mooring motion, measured by pressure recorders on the mooring lines, accounts for about 15% of the semidiurnal temperature variance, but it is incoherent with the equilibrium tide. Diurnal tides were examined at all three locations. At the MODE-1 site - near the critical latitude for diurnal period internal waves - the current and temperture fields are dominated by high mode, incoherent, inertial-character morions which mask the tidal currents. About 25% of the diurnal band temperature variance is related to mooring motion. Muir Seamount provides a clear example of diurnal period internal tides trapped to their source region north of the critical latitude. A simple analytical model is developed for the diurnal period flow adjustment in a seamount geometry. Site D shows some evidence for diurnal period internal tides, but most of the energy in the diurnal currents is not simply related to the tidal forcing. Diurnal barotropic currents measured at Site D are combined with currents on the New England shelf, showing that the diurnal tidal wave behaves as a Kelvin-Stokes mode trapped to the slope, propagating along the depth contours to the west. Some aspects of simple generation models are considered. The slope north of Site D is not at all well described as an abrupt step for the M2 generations problem, but a more realistic model of Baines (1974) predicts the coherent fields observed. But the relatively small energy conversion from the surface tide to internal modes suggests that the globally near-critical slope north of Site D is a poor generator of internal tides in the deep sea, although the local energy density is high. The step shelf generation model is well suited for the steep Blake Escarpment, and predicts a seaward energy flux of .4 x 108 erg/s cm in the first mode, comparable to the measurements at MODE-1. This confirms theoretical expectations that the first baroclinic mode is not significantly damped by turbulent diffusion after propagating through the 700 km of ocean between the generation region and the MODE-1 deep ocean site.

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

Description: A set of vertical profiles of horizontal ocean currents, obtained by electro-magnetic profilers in the Atlantic Ocean southwest of Bermuda in the spring of 1973, has been analyzed in order to study the vertical structure and temporal behavior of internal waves, particularly those with periods near the local inertial period. An important feature of the observed structure is the polarization of horizontal velocity components in the vertical. This polarization, along with temporal changes of the vertical wave structure seen in a time series of profiles made at one location, has been related to the direction of vertical energy flux due to the observed waves. Whereas the observed vertical phase propagation can be affected by horizontal advection of waves past the point of observation, the use of wave polarization to infer the direction of vertical energy propagation has the advantage that it is not influenced by horizontal advection. The result shows that at a location where profiles were obtained over smooth topography, the net energy flux was downward, indicating that the energy sources for these waves were located at or near the sea surface. An estimate of the net, downward energy flux (~ .2 - .3 erg/cm2/sec) has been obtained. Calculations have been made which show that a frictional bottom boundary layer can be an important energy sink for near-inertial waves. A rough estimate suggests that the observed, net, downward energy flux coul d be accounted for by energy losses in this frictional boundary layer. A reflection coefficient for the observed waves as they reflect off the bottom has been estimated. In contrast, some profiles made over a region of rough topography indicate that the rough bottom may also be acting to generate near-inertial waves which propagate energy upward. Ca1culations of vertical flux of horizontal kinetic energy, using an empirical form for the energy spectrum of internal waves, show that this vertical flux reaches a maximum for frequencies 10% - 20% greater than the local inertial frequency. Comparison with profiler velocity data and frequency spectra supports the conclusion that the dominant waves had frequencies 10% - 20% greater than the inertial frequency. The fact that the waves were propagating energy in the vertical is proposed as the reason for the observed frequency shift. Finally, energy spectra in vertical wave number have been calculated from the profiles in order to compare the data with an empirical model of the energy density spectrum for internal waves proposed by C. Garrett and W. Munk (1975). The result shows that although the general shape and magnitude of the observed spectrum compares well with the empirica1 model, the two-sided spectrum is not symmetric in vertical wave number. This asymmetry has been used to infer that more energy was propagating downward than upward. These calculations have also been used to obtain the coherence between profiles made at the same location, but separated in time (the so-called dropped, lagged, rotary coherence). This coherence is compared with the aforementioned empirical model. The coherence results show that the contribution of the semidiurnal tide to the energy of the profiles is restricted to long vertical wave lengths.

Description: Support for the experiment which is described in this report was provided by the Office of Naval Research under contracts N00014-66-C-0241, NR 083-004 and N00014-74-C-0262, NR 083-004.

Description: Author Posting. © American Meteorological Society, 2009. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Atmospheric and Oceanic Technology 26 (2009): 2228-2242, doi:10.1175/2009JTECHO652.1.

Description: The performance of pressure sensor–equipped inverted echo sounders for monitoring nonlinear internal waves is examined. The inverted echo sounder measures the round-trip acoustic travel time from the sea floor to the sea surface and thus acquires vertically integrated information on the thermal structure, from which the first baroclinic mode of thermocline motion may be inferred. This application of the technology differs from previous uses in that the wave period (30 min) is short, requiring a more rapid transmission rate and a different approach to the analysis. Sources of error affecting instrument performance include tidal effects, barotropic adjustment to internal waves, ambient acoustic noise, and sea surface roughness. The latter two effects are explored with a simulation that includes surface wave reconstruction, acoustic scattering based on the Kirchhoff approximation, wind-generated noise, sound propagation, and the instrument’s signal processing circuitry. Bias is introduced as a function of wind speed, but the simulation provides a basis for bias correction. The assumption that the waves do not significantly affect the mean stratification allows for a focus on the dynamic response. Model calculations are compared with observations in the South China Sea by using nearby temperature measurements to provide a test of instrument performance. After applying corrections for ambient noise and surface roughness effects, the inverted echo sounder exhibits an RMS variability of approximately 4 m in the estimated depth of the eigenfunction maximum in the wind speed range 0 ≤ U10 ≤ 10 m s−1. This uncertainty may be compared with isopycnal excursions for nonlinear internal waves of 100 m, showing that the observational approach is effective for measurements of nonlinear internal waves in this environment.

Description: This project was supported by the ONR Nonlinear Wave Program under Contract N0014-05-1-0286.

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): 602–615, doi:10.1175/JPO-D-12-055.1.

Description: The ocean interior stratification and meridional overturning circulation are largely sustained by diapycnal mixing. The breaking of internal tides is a major source of diapycnal mixing. Many recent climate models parameterize internal-tide breaking using the scheme of St. Laurent et al. While this parameterization dynamically accounts for internal-tide generation, the vertical distribution of the resultant mixing is ad hoc, prescribing energy dissipation to decay exponentially above the ocean bottom with a fixed-length scale. Recently, Polzin formulated a dynamically based parameterization, in which the vertical profile of dissipation decays algebraically with a varying decay scale, accounting for variable stratification using Wentzel–Kramers–Brillouin (WKB) stretching. This study compares two simulations using the St. Laurent and Polzin formulations in the Climate Model, version 2G (CM2G), ocean–ice–atmosphere coupled model, with the same formulation for internal-tide energy input. Focusing mainly on the Pacific Ocean, where the deep low-frequency variability is relatively small, the authors show that the ocean state shows modest but robust and significant sensitivity to the vertical profile of internal-tide-driven mixing. Therefore, not only the energy input to the internal tides matters, but also where in the vertical it is dissipated.

Description: This work is a component of the Internal- Wave Driven Mixing Climate Process Team funded by the National Science Foundation Grant OCE-0968721 and the National Oceanic and Atmospheric Administration, U.S. Department of Commerce, Award NA08OAR4320752.

Description: 2013-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 46 (2016): 417-437, doi:10.1175/JPO-D-15-0055.1.

Description: In the stratified ocean, turbulent mixing is primarily attributed to the breaking of internal waves. As such, internal waves provide a link between large-scale forcing and small-scale mixing. The internal wave field north of the Kerguelen Plateau is characterized using 914 high-resolution hydrographic profiles from novel Electromagnetic Autonomous Profiling Explorer (EM-APEX) floats. Altogether, 46 coherent features are identified in the EM-APEX velocity profiles and interpreted in terms of internal wave kinematics. The large number of internal waves analyzed provides a quantitative framework for characterizing spatial variations in the internal wave field and for resolving generation versus propagation dynamics. Internal waves observed near the Kerguelen Plateau have a mean vertical wavelength of 200 m, a mean horizontal wavelength of 15 km, a mean period of 16 h, and a mean horizontal group velocity of 3 cm s−1. The internal wave characteristics are dependent on regional dynamics, suggesting that different generation mechanisms of internal waves dominate in different dynamical zones. The wave fields in the Subantarctic/Subtropical Front and the Polar Front Zone are influenced by the local small-scale topography and flow strength. The eddy-wave field is influenced by the large-scale flow structure, while the internal wave field in the Subantarctic Zone is controlled by atmospheric forcing. More importantly, the local generation of internal waves not only drives large-scale dissipation in the frontal region but also downstream from the plateau. Some internal waves in the frontal region are advected away from the plateau, contributing to mixing and stratification budgets elsewhere.

Description: A.M. was supported by the joint CSIRO-University of Tasmania Quantitative Marine Science (QMS) program and the 2009 CSIRO Wealth from Ocean Flagship Collaborative Fund. K.L.P.’s salary support was provided by Woods Hole Oceanographic Institution bridge support funds. B.M.S. was supported by the Australian Climate Change Science Program.

Description: 2016-06-07

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, 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): 1306–1328, doi:10.1175/JPO-D-12-0191.1.

Description: The ice–ocean system is investigated on inertial to monthly time scales using winter 2009–10 observations from the first ice-tethered profiler (ITP) equipped with a velocity sensor (ITP-V). Fluctuations in surface winds, ice velocity, and ocean velocity at 7-m depth were correlated. Observed ocean velocity was primarily directed to the right of the ice velocity and spiraled clockwise while decaying with depth through the mixed layer. Inertial and tidal motions of the ice and in the underlying ocean were observed throughout the record. Just below the ice–ocean interface, direct estimates of the turbulent vertical heat, salt, and momentum fluxes and the turbulent dissipation rate were obtained. Periods of elevated internal wave activity were associated with changes to the turbulent heat and salt fluxes as well as stratification primarily within the mixed layer. Turbulent heat and salt fluxes were correlated particularly when the mixed layer was closest to the freezing temperature. Momentum flux is adequately related to velocity shear using a constant ice–ocean drag coefficient, mixing length based on the planetary and geometric scales, or Rossby similarity theory. Ekman viscosity described velocity shear over the mixed layer. The ice–ocean drag coefficient was elevated for certain directions of the ice–ocean shear, implying an ice topography that was characterized by linear ridges. Mixing length was best estimated using the wavenumber of the beginning of the inertial subrange or a variable drag coefficient. Analyses of this and future ITP-V datasets will advance understanding of ice–ocean interactions and their parameterizations in numerical models.

Description: Support for this study and the overall ITP program was provided by the National Science Foundation and Woods Hole Oceanographic Institution. Support for S. Cole was partially though the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Devonshire Foundation.

Description: 2014-11-01

Description: Author Posting. © American Meteorological Society, 2016. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 29 (2016): 1545-1571, doi:10.1175/JCLI-D-15-0509.1.

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

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

Description: 2016-08-19

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 51(10), (2021): 3235–3252, https://doi.org/10.1175/JPO-D-20-0288.1.

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

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

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 51(9), (2021): 2721–2733, https://doi.org/10.1175/JPO-D-20-0298.1.

Description: A linear numerical model of an island or a tall seamount is used to explore superinertial leaky resonances forced by ambient vertically and horizontally uniform current fluctuations. The model assumes a circularly symmetric topography (including a shallow reef) and allows realistic stratification and bottom friction. As long as there is substantial stratification, a number of leaky resonances are found, and when the island’s flanks are narrow relative to the internal Rossby radius, some of the near-resonant modes resemble leaky internal Kelvin waves. Other “resonances” resemble higher radial mode long gravity waves as explored by Chambers. The near-resonances amplify the cross-reef velocities that help fuel biological activity. Results for cases with the central island replaced by a lagoon do not differ greatly from the island case which has land at the center. As an aside, insight is provided on the question of offshore boundary conditions for superinertial nearly trapped waves along a straight coast.

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 34(22), (2021): 8971–8987, https://doi.org/10.1175/JCLI-D-20-0610.1.

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

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

Description: Author Posting. © American Meteorological Society, 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): 2251-2270, doi:10.1175/JPO-D-19-0303.1.

Description: The Gulf Stream affects global climate by transporting water and heat poleward. The current’s volume transport increases markedly along the U.S. East Coast. An extensive observing program using autonomous underwater gliders provides finescale, subsurface observations of hydrography and velocity spanning more than 15° of latitude along the path of the Gulf Stream, thereby filling a 1500-km-long gap between long-term transport measurements in the Florida Strait and downstream of Cape Hatteras. Here, the glider-based observations are combined with shipboard measurements along Line W near 68°W to provide a detailed picture of the along-stream transport increase. To account for the influences of Gulf Stream curvature and adjacent circulation (e.g., corotating eddies) on transport estimates, upper- and lower-bound transports are constructed for each cross–Gulf Stream transect. The upper-bound estimate for time-averaged volume transport above 1000 m is 32.9 ± 1.2 Sv (1 Sv ≡ 106 m3 s−1) in the Florida Strait, 57.3 ± 1.9 Sv at Cape Hatteras, and 75.6 ± 4.7 Sv at Line W. Corresponding lower-bound estimates are 32.3 ± 1.1 Sv in the Florida Strait, 54.5 ± 1.7 Sv at Cape Hatteras, and 69.9 ± 4.2 Sv at Line W. Using the temperature and salinity observations from gliders and Line W, waters are divided into seven classes to investigate the properties of waters that are transported by and entrained into the Gulf Stream. Most of the increase in overall Gulf Stream volume transport above 1000 m stems from the entrainment of subthermocline waters, including upper Labrador Sea Water and Eighteen Degree Water.

Description: We gratefully acknowledge funding from the Office of Naval Research (N000141713040), the National Science Foundation (OCE-0220769, OCE-1633911, OCE-1923362), NOAA’s Global Ocean Monitoring and Observing Program (NA14OAR4320158, NA19OAR4320074), WHOI’s Oceans and Climate Change Institute, Eastman Chemical Company, and the W. Van Alan Clark, Jr. Chair for Excellence in Oceanography at WHOI (awarded to Breck Owens).

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(9),(2020): 2797-2814, https://doi.org/10.1175/JPO-D-19-0326.1.

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

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

Description: Author Posting. © American Meteorological Society, 2019. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 32(13), (2019): 3883-3898, doi:10.1175/JCLI-D-18-0735.1.

Description: While it has generally been understood that the production of Labrador Sea Water (LSW) impacts the Atlantic meridional overturning circulation (MOC), this relationship has not been explored extensively or validated against observations. To explore this relationship, a suite of global ocean–sea ice models forced by the same interannually varying atmospheric dataset, varying in resolution from non-eddy-permitting to eddy-permitting (1°–1/4°), is analyzed to investigate the local and downstream relationships between LSW formation and the MOC on interannual to decadal time scales. While all models display a strong relationship between changes in the LSW volume and the MOC in the Labrador Sea, this relationship degrades considerably downstream of the Labrador Sea. In particular, there is no consistent pattern among the models in the North Atlantic subtropical basin over interannual to decadal time scales. Furthermore, the strong response of the MOC in the Labrador Sea to LSW volume changes in that basin may be biased by the overproduction of LSW in many models compared to observations. This analysis shows that changes in LSW volume in the Labrador Sea cannot be clearly and consistently linked to a coherent MOC response across latitudes over interannual to decadal time scales in ocean hindcast simulations of the last half century. Similarly, no coherent relationships are identified between the MOC and the Labrador Sea mixed layer depth or the density of newly formed LSW across latitudes or across models over interannual to decadal time scales.

Description: FL and MSL are thankful for the financial support from the National Science Foundation (NSF) Physical Oceanography Program (NSF-OCE-12-59102, NSF-OCE-12-59103). The NCAR contribution was supported by the National Oceanic and Atmospheric Administration (NOAA) Climate Program Office (CPO) under Climate Variability and Predictability Program (CVP) Grant NA13OAR4310138 and by the NSF Collaborative Research EaSM2 Grant OCE-1243015. NCAR is sponsored by the NSF. NPH is supported by NERC programs U.K. OSNAP (NE/K010875) and ACSIS (National Capability, NE/N018044/1). Y-OK is supported by NOAA CPO CVP (NA17OAR4310111) and NSF EaSM2 grant (OCE-1242989). AR is supported by NASA-ROSES Modeling, Analysis and Prediction 2016 NNX16AC93G-MAP. RZ is supported by NOAA/OAR. Argo data were collected and made freely available by the International Argo Program and the national programs that contribute to it (http://www.argo.ucsd.edu, http://argo.jcommops.org). The Argo Program is part of the Global Ocean Observing System (http://doi.org/10.17882/42182). Data from the RAPID-MOCHA-WBTS array funded by NERC, NSF and NOAA are freely available from www.rapid.ac.uk/rapidmoc. We thank Stephen Griffies for providing access to the GFDL-MOM025 COREII simulation output and Matthew Harrison and Xiaoqin Yan for their comments on the manuscript. We also thank the anonymous reviewers for their valuable comments.

Description: 2020-06-11

Description: Author Posting. © American Meteorological Society, 2019. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 49(11), (2019): 2781-2797, doi: 10.1175/JPO-D-19-0111.1.

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

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

Description: 2020-04-17

Description: Author Posting. © American Meteorological Society, 2019. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 49(12), (2019): 3061-3068, doi: 10.1175/JPO-D-18-0172.1.

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

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

Description: 2020-06-03

Description: Author Posting. © American Meteorological Society, 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): 2611-2630, doi:10.1175/JPO-D-16-0259.1.

Description: This study reports the results of large-eddy simulations of an axisymmetric turbulent buoyant plume in a stratified fluid. The configuration used is an idealized model of the plume generated by a subglacial discharge at the base of a tidewater glacier with an ambient stratification typical of Greenland fjords. The plume is discharged from a round source of various diameters and characteristic stratifications for summer and winter are considered. The classical theory for the integral parameters of a turbulent plume in a homogeneous fluid gives accurate predictions in the weakly stratified lower layer up to the pycnocline, and the plume dynamics are not sensitive to changes in the source diameter. In winter, when the stratification is similar to an idealized two-layer case, turbulent entrainment and generation of internal waves by the plume top are in agreement with the theoretical and numerical results obtained for turbulent jets in a two-layer stratification. In summer, instead, the stratification is more complex and turbulent entrainment by the plume top is significantly reduced. The subsurface layer in summer is characterized by a strong density gradient and the oscillating plume generates internal waves that might serve as an indicator of submerged plumes not penetrating to the surface.

Description: This work was supported by Linné FLOW Centre at KTH and the Academy of Finland Centre of Excellence program (Grant 307331) (E. E.) and VR Swedish Research Council, Outstanding Young Researcher Award, Grant VR 2014-5001 (L. B.). Support to C. C. was given by the NSF Project OCE-1434041.

Description: 2018-04-26

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): 2479-2498, doi:10.1175/JPO-D-16-0167.1.

Description: The generation of trapped and radiating internal tides around Izu‐Oshima Island located off Sagami Bay, Japan, is investigated using the three-dimensional Stanford Unstructured Nonhydrostatic Terrain-following Adaptive Navier–Stokes Simulator (SUNTANS) that is validated with observations of isotherm displacements in shallow water. The model is forced by barotropic tides, which generate strong baroclinic internal tides in the study region. Model results showed that when diurnal K1 barotropic tides dominate, resonance of a trapped internal Kelvin wave leads to large-amplitude internal tides in shallow waters on the coast. This resonance produces diurnal motions that are much stronger than the semidiurnal motions. The weaker, freely propagating, semidiurnal internal tides are generated on the western side of the island, where the M2 internal tide beam angle matches the topographic slope. The internal wave energy flux due to the diurnal internal tides is much higher than that of the semidiurnal tides in the study region. Although the diurnal internal tide energy is trapped, this study shows that steepening of the Kelvin waves produces high-frequency internal tides that radiate from the island, thus acting as a mechanism to extract energy from the diurnal motions.

Description: This study was supported by JST CREST Grant Number JPRMJCR12A6.

Description: 2018-04-12

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): 2531-2543, doi:10.1175/JPO-D-17-0051.1.

Description: Argo floats are used to investigate Labrador Sea overturning and its variability on seasonal time scales. This is the first application of Argo floats to estimate overturning in a deep-water formation region in the North Atlantic. Unlike hydrographic measurements, which are typically confined to the summer season, floats offer the advantage of collecting data in all seasons. Seasonal composite potential density and absolute geostrophic velocity sections across the mouth of the Labrador Sea assembled from float profiles and trajectories at 1000 m are used to calculate the horizontal and overturning circulations. The overturning exhibits a pronounced seasonal cycle; in depth space the overturning doubles throughout the course of the year, and in density space it triples. The largest overturning [1.2 Sv (1 Sv ≡ 106 m3 s−1) in depth space and 3.9 Sv in density space] occurs in spring and corresponds to the outflow of recently formed Labrador Sea Water. The overturning decreases through summer and reaches a minimum in winter (0.6 Sv in depth space and 1.2 Sv in density space). The robustness of the Argo seasonal overturning is supported by a comparison to an overturning estimate based on hydrographic data from the AR7W line.

Description: NSF OCE-1459474 supported this work.

Description: 2018-04-17

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(5), (2022): 595–617, https://doi.org/10.1175/jtech-d-21-0039.1.

Description: The future Surface Water and Ocean Topography (SWOT) mission aims to map sea surface height (SSH) in wide swaths with an unprecedented spatial resolution and subcentimeter accuracy. The instrument performance needs to be verified using independent measurements in a process known as calibration and validation (Cal/Val). The SWOT Cal/Val needs in situ measurements that can make synoptic observations of SSH field over an O(100) km distance with an accuracy matching the SWOT requirements specified in terms of the along-track wavenumber spectrum of SSH error. No existing in situ observing system has been demonstrated to meet this challenge. A field campaign was conducted during September 2019–January 2020 to assess the potential of various instruments and platforms to meet the SWOT Cal/Val requirement. These instruments include two GPS buoys, two bottom pressure recorders (BPR), three moorings with fixed conductivity–temperature–depth (CTD) and CTD profilers, and a glider. The observations demonstrated that 1) the SSH (hydrostatic) equation can be closed with 1–3 cm RMS residual using BPR, CTD mooring and GPS SSH, and 2) using the upper-ocean steric height derived from CTD moorings enable subcentimeter accuracy in the California Current region during the 2019/20 winter. Given that the three moorings are separated at 10–20–30 km distance, the observations provide valuable information about the small-scale SSH variability associated with the ocean circulation at frequencies ranging from hourly to monthly in the region. The combined analysis sheds light on the design of the SWOT mission postlaunch Cal/Val field campaign.

Description: The research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (80NM0018D0004). All authors are supported by the SWOT project. J. T. Farrar was partially supported by NASA NNX16AH76G.

Description: 2022-11-01

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

Description: A hydrostatic, coupled-mode, shallow-water model (CSW) is described and used to diagnose and simulate tidal dynamics in the greater Mid-Atlantic Bight region. The reduced-physics model incorporates realistic stratification and topography, internal tide forcing from a priori estimates of the surface tide, and advection terms that describe first-order interactions of internal tides with slowly varying mean flow and mean buoyancy fields and their respective shear. The model is validated via comparisons with semianalytic models and nonlinear primitive equation models in several idealized and realistic simulations that include internal tide interactions with topography and mean flows. Then, 24 simulations of internal tide generation and propagation in the greater Mid-Atlantic Bight region are used to diagnose significant internal tide interactions with the Gulf Stream. The simulations indicate that locally generated mode-one internal tides refract and/or reflect at the Gulf Stream. The redirected internal tides often reappear at the shelf break, where their onshore energy fluxes are intermittent (i.e., noncoherent with surface tide) because meanders in the Gulf Stream alter their precise location, phase, and amplitude. These results provide an explanation for anomalous onshore energy fluxes that were previously observed at the New Jersey shelf break and linked to the irregular generation of nonlinear internal waves.

Description: We thank the National Science Foundation for support under Grant OCE-1061160 (ShelfIT) to the Massachusetts Institute of Technology (MIT) and under Grant OCE-1060430 to the Woods Hole Oceanographic Institution. PFJL and PJH also thank the Office of Naval Research for research support under Grants N00014-11-1-0701 (MURI-IODA), N00014-12-1-0944 (ONR6.2), and N00014-13-1-0518 (Multi-DA) to MIT.

Description: 2017-06-14

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): 85-100, doi:10.1175/JPO-D-15-0234.1.

Description: Observations and analyses of two tidally recurring, oblique, internal hydraulic jumps at a stratified estuary mouth (Columbia River, Oregon/Washington) are presented. These hydraulic features have not previously been studied due to the challenges of both horizontally resolving the sharp gradients and temporally resolving their evolution in numerical models and traditional observation platforms. The jumps, both of which recurred during ebb, formed adjacent to two engineered lateral channel constrictions and were identified in marine radar image time series. Jump occurrence was corroborated by (i) a collocated sharp gradient in the surface currents measured via airborne along-track interferometric synthetic aperture radar and (ii) the transition from supercritical to subcritical flow in the cross-jump direction via shipborne velocity and density measurements. Using a two-layer approximation, observed jump angles at both lateral constrictions are shown to lie within the theoretical bounds given by the critical internal long-wave (Froude) angle and the arrested maximum-amplitude internal bore angle, respectively. Also, intratidal and intertidal variability of the jump angles are shown to be consistent with that expected from the two-layer model, applied to varying stratification and current speed over a range of tidal and river discharge conditions. Intratidal variability of the upchannel jump angle is similar under all observed conditions, whereas the downchannel jump angle shows an additional association with stratification and ebb velocity during the low discharge periods. The observations additionally indicate that the upchannel jump achieves a stable position that is collocated with a similarly oblique bathymetric slope.

Description: We acknowledge the financial support of the Office of Naval Research under Awards N00014-10-1-0932 and N00014-13-1-0364.

Description: 2017-07-04

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): 2999-3013, doi:10.1175/JPO-D-17-0129.1.

Description: Initial results are presented from a yearlong, high-resolution (~2 km) numerical simulation covering the east Greenland shelf and the Iceland and Irminger Seas. The model hydrography and circulation in the vicinity of Denmark Strait show good agreement with available observational datasets. This study focuses on the variability of the Denmark Strait overflow (DSO) by detecting and characterizing boluses and pulses, which are the two dominant mesoscale features in the strait. The authors estimate that the yearly mean southward volume flux of the DSO is about 30% greater in the presence of boluses and pulses. On average, boluses (pulses) are 57.1 (27.5) h long, occur every 3.2 (5.5) days, and are more frequent during the summer (winter). Boluses (pulses) increase (decrease) the overflow cross-sectional area, and temperatures around the overflow interface are colder (warmer) by about 2.6°C (1.8°C). The lateral extent of the boluses is much greater than that of the pulses. In both cases the along-strait equatorward flow of dense water is enhanced but more so for pulses. The sea surface height (SSH) rises by 4–10 cm during boluses and by up to 5 cm during pulses. The SSH anomaly contours form a bowl (dome) during boluses (pulses), and the two features cross the strait with a slightly different orientation. The cross streamflow changes direction; boluses (pulses) are associated with veering (backing) of the horizontal current. The model indicates that boluses and pulses play a major role in controlling the variability of the DSO transport into the Irminger Sea.

Description: This work was supported by the NSF Grants OCE-1433448, OCE-1633124, and OCE- 1259618 and the Institute for Data Intensive Engineering and Science (IDIES) seed grant funding.

Description: 2018-06-13

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

Description: Global patterns of ocean salinity arise from the exchange of freshwater between the sea surface and the atmosphere. For a quasi-steady state system, these surface fluxes are balanced by compensating transports of salt in the ocean interior. In a warming climate, the atmosphere holds additional water vapor which acts to intensify the global water cycle. Amplified freshwater fluxes are then absorbed at the surface and propagate along ocean circulation pathways. Here, we use coupled model results from the CMIP5 experiment to identify coherent responses in the atmospheric water cycle and in ocean salinity patterns. Some aspects of the response are consistent across models, while other regions show large inter-model spread. In particular, the salinity response in the North Atlantic subpolar gyre, where the mean salinity plays a role in maintaining high surface density for deep-water formation, has low confidence in CMIP5 models. To understand how differences in ocean circulation may affect this response, we use two techniques to diagnose the role of salt transports in the present-day climate. The first is a salt budget within the surface mixed layer, which identifies major transport processes. The second is a Lagrangian particle tracking tool, used to understand the regional connectivity of water masses. From this analysis, we find that anomalous freshwater signals become well mixed within the ocean gyres, but can be isolated on larger scales. The subpolar Atlantic salinity response generally shows freshening at the surface, but is sensitive to the transport of anomalously salty water from the subtropics, a largely eddy-driven process. As CMIP5 models use a range of eddy parameterizations, this is likely a source of uncertainty in the salinity response. Finally, we investigate the effect of salinity changes on the deep overturning cells and other circulations, and find a complex influence that also depends on the details of advective pathways. In a warming scenario, water cycle amplification actually works to strengthen the Atlantic meridional overturning circulation due to the influence of enhanced subtropical evaporation.

Description: Funding for this thesis was provided by NASA grant NNX12AF59GS03, a NASA Earth and Space Science Fellowship award 80NSSC17K0372, and the WHOI Academic Programs Office.

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): 1639-1649, doi: 10.1175/JPO-D-18-0154.1.

Description: Using a recently developed asymptotic theory of internal solitary wave propagation over a sloping bottom in a rotating ocean, some new qualitative and quantitative features of this process are analyzed for internal waves in a two-layer ocean. The interplay between different singularities—terminal damping due to radiation and disappearing quadratic nonlinearity, and reaching an “internal beach” (e.g., zero lower-layer depth)—is discussed. Examples of the adiabatic evolution of a single solitary wave over a uniformly sloping bottom under realistic conditions are considered in more detail and compared with numerical solutions of the variable-coefficient, rotation-modified Korteweg–de Vries (rKdV) equation.

Description: LAO is thankful to Yu. Stepanyants for broad discussions of mutual benefit. KRH was supported by Grant N00014-18-1-2542 from the Office of Naval Research.

Description: 2020-06-13

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 Climate 30 (2017): 9679-9702, doi:10.1175/JCLI-D-16-0707.1.

Description: The North Atlantic atmospheric circulation response to the meridional shifts of the Gulf Stream (GS) path is examined using a large ensemble of high-resolution hemispheric-scale Weather Research and Forecasting Model simulations. The model is forced with a broad range of wintertime sea surface temperature (SST) anomalies derived from a lag regression on a GS index. The primary result of the model experiments, supported in part by an independent analysis of a reanalysis dataset, is that the large-scale quasi-steady North Atlantic circulation response is remarkably nonlinear about the sign and amplitude of the SST anomaly chosen over a wide range of GS shift scenarios. The nonlinear response prevails over the weak linear response and resembles the negative North Atlantic Oscillation (NAO), the leading intrinsic mode of variability in the model and the observations. Further analysis of the associated dynamics reveals that the nonlinear responses are accompanied by the shift of the North Atlantic eddy-driven jet, which is reinforced, with nearly equal importance, by the high-frequency transient eddy feedback and the low-frequency wave-breaking events. Additional sensitivity simulations confirm that the nonlinearity of the circulation response is a robust feature found over the broad parameter space encompassing not only the varied SST but also the absence/presence of tropical influence, the varying lateral boundary conditions, and the initialization scheme. The result highlights the fundamental importance of the intrinsically nonlinear transient eddy dynamics and the eddy–mean flow interactions in generating the nonlinear downstream response to the meridional shifts in the Gulf Stream.

Description: The authors are grateful for the support from NASA (NNX13AM59G) and the NSF (AGS-1355339, OCE-1419235).

Description: 2018-05-07

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

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

Description: 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): 1116–1132, doi:10.1175/JPO-D-13-0194.1.

Description: Internal solitary waves commonly observed in the coastal ocean are often modeled by a nonlinear evolution equation of the Korteweg–de Vries type. Because these waves often propagate for long distances over several inertial periods, the effect of Earth’s background rotation is potentially significant. The relevant extension of the Kortweg–de Vries is then the Ostrovsky equation, which for internal waves does not support a steady solitary wave solution. Recent studies using a combination of asymptotic theory, numerical simulations, and laboratory experiments have shown that the long time effect of rotation is the destruction of the initial internal solitary wave by the radiation of small-amplitude inertia–gravity waves, and the eventual emergence of a coherent, steadily propagating, nonlinear wave packet. However, in the ocean, internal solitary waves are often propagating over variable topography, and this alone can cause quite dramatic deformation and transformation of an internal solitary wave. Hence, the combined effects of background rotation and variable topography are examined. Then the Ostrovsky equation is replaced by a variable coefficient Ostrovsky equation whose coefficients depend explicitly on the spatial coordinate. Some numerical simulations of this equation, together with analogous simulations using the Massachusetts Institute of Technology General Circulation Model (MITgcm), for a certain cross section of the South China Sea are presented. These demonstrate that the combined effect of shoaling and rotation is to induce a secondary trailing wave packet, induced by enhanced radiation from the leading wave.

Description: KH was supported by Grants N00014-09-10227 and N00014-11-0701 from the Office of Naval Research.

Description: 2014-10-01

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

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

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

Description: 2015-06-15

Description: Author Posting. © American Meteorological Society, 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, 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.