ALBERT

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 he balance of salinity variance in a partially stratified estuary: Implications for exchange flow, mixing, and stratification. Journal of Physical Oceanography, 48(12), (2018) 2887-2899., doi: 10.1175/JPO-D-18-0032.1.

Description: Salinity variance dissipation is related to exchange flow through the salinity variance balance equation, and meanwhile its magnitude is also proportional to the turbulence production and stratification inside the estuary. As river flow increases, estuarine volume-integrated salinity variance dissipation increases owing to more variance input from the open boundaries driven by exchange flow and river flow. This corresponds to the increased efficient conversion of turbulence production to salinity variance dissipation due to the intensified stratification with higher river flow. Through the spring–neap cycle, the temporal variation of salinity variance dissipation is more dependent on stratification than turbulence production, so it reaches its maximum during the transition from neap to spring tides. During most of the transition time from spring to neap tides, the advective input of salinity variance from the open boundaries is larger than dissipation, resulting in the net increase of variance, which is mainly expressed as vertical variance, that is, stratification. The intensified stratification in turn increases salinity variance dissipation. During neap tides, a large amount of enhanced salinity variance dissipation is induced by the internal shear stress near the halocline. During most of the transition time from neap to spring tides, dissipation becomes larger than the advective input, so salinity variance decreases and the stratification is destroyed.

Description: TW was supported by the National Key R&D Program of China (Grant 2017YFA0604104), National Natural Science Foundation of China (Grant 41706002), Natural Science Foundation of Jiangsu Province (Grant BK20170864), and MEL Visiting Fellowship (MELRS1617). WRG was supported by NSF Grant OCE 1736539. Part of this work is finished during TW’s visit in MEL and WHOI. We would like to acknowledge John Warner for providing the codes of the Hudson estuary model, and Parker MacCready, the editor, and two reviewers for their insightful suggestions on improving the manuscript.

Description: 2019-06-06

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

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

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

Description: 2020-02-18

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

Description: The main source feeding the abyssal circulation of the North Pacific is the deep, northward flow of 5–6 Sverdrups (Sv; 1 Sv ≡ 106 m3 s−1) through the Samoan Passage. A recent field campaign has shown that this flow is hydraulically controlled and that it experiences hydraulic jumps accompanied by strong mixing and dissipation concentrated near several deep sills. By our estimates, the diapycnal density flux associated with this mixing is considerably larger than the diapycnal flux across a typical isopycnal surface extending over the abyssal North Pacific. According to historical hydrographic observations, a second source of abyssal water for the North Pacific is 2.3–2.8 Sv of the dense flow that is diverted around the Manihiki Plateau to the east, bypassing the Samoan Passage. This bypass flow is not confined to a channel and is therefore less likely to experience the strong mixing that is associated with hydraulic transitions. The partitioning of flux between the two branches of the deep flow could therefore be relevant to the distribution of Pacific abyssal mixing. To gain insight into the factors that control the partitioning between these two branches, we develop an abyssal and equator-proximal extension of the “island rule.” Novel features include provisions for the presence of hydraulic jumps as well as identification of an appropriate integration circuit for an abyssal layer to the east of the island. Evaluation of the corresponding circulation integral leads to a prediction of 0.4–2.4 Sv of bypass flow. The circulation integral clearly identifies dissipation and frictional drag effects within the Samoan Passage as crucial elements in partitioning the flow.

Description: This work was supported by the National Science Foundation under Grants OCE-1029268, OCE-1029483, OCE-1657264, OCE-1657870, OCE-1658027, and OCE-1657795. We thank the captain, crew, and engineers at APL/UW for their hard work and skill.

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

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

Description: 2020-02-28

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

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

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

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

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

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

Description: 2020-05-29

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

Description: Many chemical constituents are removed from the ocean by attachment to settling particles, a process referred to as “scavenging.” Radioisotopes of thorium, a highly particle-reactive element, have been used extensively to study scavenging in the ocean. However, this process is complicated by the highly variable chemical composition and concentration of particles in oceanic waters. This thesis focuses on understanding the cycling of thorium as affected by particle concentration and particle composition in the North Atlantic. This objective is addressed using (i) the distributions 228,230,234Th, their radioactive parents, particle composition, and bulk particle concentration, as measured or estimated along the GEOTRACES North Atlantic Transect (GA03) and (ii) a model for the reversible exchange of thorium with particles. Model parameters are either estimated by inversion (chapter 2-4), or prescribed in order to simulate 230Th in a circulation model (chapter 5). The major findings of this thesis follow. In chapters 2 and 3, I find that the rate parameters of the reversible exchange model show systematic variations along GA03. In particular, 𝑘1, the apparent first-order rate "constant" of Th adsorption onto particles, generally presents maxima in the mesopelagic zone and minima below. A positive correlation between 𝑘1 and bulk particle concentration is found, consistent with the notion that the specific rate at which a metal in solution attaches to particles increases with the number of surface sites available for adsorption. In chapter 4, I show that Mn (oxyhydr)oxides and biogenic particles most strongly influence 𝑘1 west of the Mauritanian upwelling, but that biogenic particles dominate 𝑘1 in this region. In chapter 5, I find that dissolved 230Th data are best represented by a model that assumes enhanced values of 𝑘1 near the seafloor. Collectively, my findings suggest that spatial variations in Th radioisotope activities observed in the North Atlantic reflect at least partly variations in the rate at which Th is removed from the water column.

Description: This work was supported by the US National Science Foundation. Two US NSF grants have supported the research in this thesis (OCE-1232578 and OCE-155644).

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): 607-623, doi:10.1175/JPO-D-17-0189.1.

Description: The roles of straining and dissipation in controlling stratification are derived analytically using a vertical salinity variance method. Stratification is produced by converting horizontal variance to vertical variance via straining, that is, differential advection of horizontal salinity gradients, and stratification is destroyed by the dissipation of vertical variance through turbulent mixing. A numerical model is applied to the Changjiang estuary in order to demonstrate the salinity variance balance and how it reveals the factors controlling stratification. The variance analysis reveals that dissipation reaches its maximum during spring tide in the Changjiang estuary, leading to the lowest stratification. Stratification increases from spring tide to neap tide because of the increasing excess of straining over dissipation. Throughout the spring–neap tidal cycle, straining is almost always larger than dissipation, indicating a net excess of production of vertical variance relative to dissipation. This excess is balanced on average by advection, which exports vertical variance out of the estuarine region into the plume. During neap tide, tidal straining shows a general tendency of destratification during the flood tide and restratification during ebb, consistent with the one-dimensional theory of tidal straining. During spring tide, however, positive straining occurs during flood because of the strong baroclinicity induced by the intensified horizontal salinity gradient. These results indicate that the salinity variance method provides a valuable approach for examining the spatial and temporal variability of stratification in estuaries and coastal environments.

Description: X. Li was supported by the China Scholarship Council. W. R. Geyer was supported by NSF Grants OCE 1736539 and OCE 1634480. J. Zhu was supported by the National Natural Science Foundation of China (41476077 and 41676083). H. Wu was supported by the National Natural Science Foundation of China (41576088 and 41776101).

Description: 2018-09-08

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

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

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

Description: 2018-10-04

Description: 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): 905-923, doi:10.1175/JPO-D-17-0133.1.

Description: Observations of turbulent kinetic energy, dissipation, and turbulent stress were collected in the middle reaches of Chesapeake Bay and were used to assess second-moment closure predictions of turbulence generated beneath breaking waves. Dissipation scaling indicates that the turbulent flow structure observed during a 10-day wind event was dominated by a three-layer response that consisted of 1) a wave transport layer, 2) a surface log layer, and 3) a tidal, bottom boundary layer limited by stable stratification. Below the wave transport layer, turbulent mixing was limited by stable stratification. Within the wave transport layer, where dissipation was balanced by a divergence in the vertical turbulent kinetic energy flux, the eddy viscosity was significantly underestimated by second-moment turbulence closure models, suggesting that breaking waves homogenized the mixed surface layer to a greater extent than the simple model of TKE diffusing away from a source at the surface. While the turbulent transport of TKE occurred largely downgradient, the intermittent downward sweeps of momentum generated by breaking waves occurred largely independent of the mean shear. The underprediction of stress in the wave transport layer by second-moment closures was likely due to the inability of the eddy viscosity model to capture the nonlocal turbulent transport of the momentum flux beneath breaking waves. Finally, the authors hypothesize that large-scale coherent turbulent eddies played a significant role in transporting momentum generated near the surface to depth.

Description: This work was supported by National Science Foundation Grants OCE-1061609 and OCE-1339032.

Description: 2018-10-19

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): 1815-1830, doi:10.1175/JPO-D-17-0275.1.

Description: Recent progress in direct numerical simulations (DNSs) of stratified turbulent flows has led to increasing attention to the validity of the constancy of the dissipation flux coefficient Γ in the Osborn’s eddy diffusivity model. Motivated by lack of observational estimates of Γ, particularly under weakly stratified deep-ocean conditions, this study estimates Γ using deep microstructure profiles collected in various regions of the North Pacific and Southern Oceans. It is shown that Γ is not constant but varies significantly with the Ozmidov/Thorpe scale ratio ROT in a fashion similar to that obtained by previous DNS studies. Efficient mixing events with Γ ~ O(1) and ROT ~ O(0.1) tend to be frequently observed in the deep ocean (i.e., weak stratification), while moderate mixing events with Γ ~ O(0.1) and ROT ~ O(1) tend to be observed in the upper ocean (i.e., strong stratification). The observed negative relationship between Γ and ROT is consistent with a simple scaling that can be derived from classical turbulence theories. In contrast, the observed results exhibit no definite relationships between Γ and the buoyancy Reynolds number Reb, although Reb has long been thought to be another key parameter that controls Γ.

Description: This study was supported by MEXT KAKENHI Grant JP15H05824 and JSPS KAKENHI Grant JP15H02131.

Description: 2019-02-15

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

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

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

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

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

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

Description: 2019-05-07

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

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

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

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

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

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

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

Description: This thesis documents the origin, distribution, and fate of methane and several of its isotopic forms on Earth. Using observational, experimental, and theoretical approaches, I illustrate how the relative abundances of 12CH4, 13CH4, 12CH3D, and 13CH3D record the formation, transport, and breakdown of methane in selected settings. Chapter 2 reports precise determinations of 13CH3D, a “clumped” isotopologue of methane, in samples collected from various settings representing many of the major sources and reservoirs of methane on Earth. The results show that the information encoded by the abundance of 13CH3D enables differentiation of methane generated by microbial, thermogenic, and abiogenic processes. A strong correlation between clumped- and hydrogen-isotope signatures in microbial methane is identified and quantitatively linked to the availability of H2 and the reversibility of microbially-mediated methanogenesis in the environment. Determination of 13CH3D in combination with hydrogen-isotope ratios of methane and water provides a sensitive indicator of the extent of C–H bond equilibration, enables fingerprinting of methane-generating mechanisms, and in some cases, supplies direct constraints for locating the waters from which migrated gases were sourced. Chapter 3 applies this concept to constrain the origin of methane in hydrothermal fluids from sediment-poor vent fields hosted in mafic and ultramafic rocks on slow- and ultraslow-spreading mid-ocean ridges. The data support a hypogene model whereby methane forms abiotically within plutonic rocks of the oceanic crust at temperatures above ca. 300 C during respeciation of magmatic volatiles, and is subsequently extracted during active, convective hydrothermal circulation. Chapter 4 presents the results of culture experiments in which methane is oxidized in the presence of O2 by the bacterium Methylococcus capsulatus strain Bath. The results show that the clumped isotopologue abundances of partially-oxidized methane can be predicted from knowledge of 13C/12C and D/H isotope fractionation factors alone.

Description: The research activities documented in this thesis were made possible by grants to my advisor from the U.S. National Science Foundation (NSF award EAR-1250394), the National Aeronautics and Space Administration (NASA) Astrobiology Institute (NAI, University of Colorado, Boulder, CAN 7 under Cooperative Agreement NNA15BB02A), the Department of Energy (DOE, Small Business Innovation Research program, contract DE-SC0004575), the Alfred P. Sloan Foundation via the Deep Carbon Observatory, and a Shell Graduate Fellowship through the MIT Energy Initiative. I completed the bulk of the work in this thesis while being supported by a National Defense Science and Engineering Graduate (NDSEG) Fellowship awarded through the Office of Naval Research of the U.S. Department of Defense. The StanleyW.Watson Fellowship Fund provided support during my first summer term at WHOI.The Charles M. Vest Presidential Fellowship at MIT supported me in the first year of my Ph.D. studies. I received additional support that year through NSF award EAR-1159318 (to S. Ono and T. Bosak) and theWalter & Adel Hohenstein Graduate Fellowship of Phi Kappa Phi. The MIT Earth Resources Laboratory and PAOC Houghton Fund funded my attendance at several conferences.

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

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

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

Description: 2017-10-12

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

Description: Dense oceanic overflows descend over the rough topography of the continental slope entraining and mixing with surrounding waters. The associated dilution dictates the fate of these currents and thus is of fundamental importance to the formation of deep water masses. The entrainment in a dense current flowing down a sloping bottom in a rotating homogeneous fluid is investigated using laboratory experiments, focusing on the influence of the bottom roughness on the flow dynamics. The roughness is idealized by an array of vertical rigid cylinders and both their spacing and height are varied as well as the inclination of the sloping bottom. The presence of the roughness is generally observed to decelerate the dense current, with a consequent reduction of the Froude number, when compared to the smooth bottom configuration. However, the dilution of the dense current due to mixing with the ambient fluid is enhanced by the roughness elements, especially for low Froude numbers. When the entrainment due to shear instability at the interface between the dense current and the ambient fluid is low, the additional turbulence and mixing arising at the bottom of the dense current due to the roughness elements strongly affects the dilution of the current. Finally, a strong dependence of the entrainment parameter on the Reynolds number is observed.

Description: Support to C. C. was given by the National Science Foundation Project OCE- 1333174. Support to L. O. during her internship at WHOI was provided by the Lions Club ‘‘Napoli Megaride’’ and the Zoological Station Anton Dohrn through the Paolo Brancaccio fellowship (2012).

Description: 2017-08-20

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): 1205-1220, doi:10.1175/JPO-D-16-0258.1.

Description: The linkage among total exchange flow, entrainment, and diffusive salt flux in estuaries is derived analytically using salinity coordinates, revealing the simple but important relationship between total exchange flow and mixing. Mixing is defined and quantified in this paper as the dissipation of salinity variance. The method uses the conservation of volume and salt to quantify and distinguish the diahaline transport of volume (i.e., entrainment) and diahaline diffusive salt flux. A numerical model of the Hudson estuary is used as an example of the application of the method in a realistic estuary with a persistent but temporally variable exchange flow. A notable finding of this analysis is that the total exchange flow and diahaline salt flux are out of phase with respect to the spring–neap cycle. Total exchange flow reaches its maximum near minimum neap tide, but diahaline salt transport reaches its maximum during the maximum spring tide. This phase shift explains the strong temporal variation of stratification and estuarine salt content through the spring–neap cycle. In addition to quantifying temporal variation, the method reveals the spatial variation of total exchange flow, entrainment, and diffusive salt flux through the estuary. For instance, the analysis of the Hudson estuary indicates that diffusive salt flux is intensified in the wider cross sections. The method also provides a simple means of quantifying numerical mixing in ocean models because it provides an estimate of the total dissipation of salinity variance, which is the sum of mixing due to the turbulence closure and numerical mixing.

Description: T. Wang was supported by the Open Research Fund of State Key Laboratory of Estuarine and Coastal Research (Grant SKLEC-KF201509), the Fundamental Research Funds for the Central Universities (Grant 2017B03514), and the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant XDA11010203). W. R. Geyer was supported by NSF Grant OCE 0926427 and ONR Grant N00014-16-1-2948. P. MacCready was supported by NSF Grant OCE-1634148.

Description: 2017-09-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): 1921-1939, doi:10.1175/JPO-D-16-0146.1.

Description: The role of surface gravity waves in structuring the air–sea momentum flux is examined in the middle reaches of Chesapeake Bay. Observed wave spectra showed that wave direction in Chesapeake Bay is strongly correlated with basin geometry. Waves preferentially developed in the direction of maximum fetch, suggesting that dominant wave frequencies may be commonly and persistently misaligned with local wind forcing. Direct observations from an ultrasonic anemometer and vertical array of ADVs show that the magnitude and direction of stress changed across the air–sea interface, suggesting that a stress divergence occurred at or near the water surface. Using a numerical wave model in combination with direct flux measurements, the air–sea momentum flux was partitioned between the surface wave field and the mean flow. Results indicate that the surface wave field can store or release a significant fraction of the total momentum flux depending on the direction of the wind. When wind blew across dominant fetch axes, the generation of short gravity waves stored as much as 40% of the total wind stress. Accounting for the storage of momentum in the surface wave field closed the air–sea momentum budget. Agreement between the direction of Lagrangian shear and the direction of the stress vector in the mixed surface layer suggests that the observed directional difference was due to the combined effect of breaking waves producing downward sweeps of momentum in the direction of wave propagation and the straining of that vorticity field in a manner similar to Langmuir turbulence.

Description: This work was supported by National Science Foundation Grants OCE-1061609 and OCE-1339032.

Description: 2018-01-13

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

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

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

Description: 2017-06-08

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

Description: Salt marshes are physically, chemically, and biologically dynamic environments found globally at temperate latitudes. Tidal creeks and marshtop ponds may expand at the expense of productive grass-covered marsh platform. It is therefore important to understand the present magnitude and drivers of production and respiration in these submerged environments in order to evaluate the future role of salt marshes as a carbon sink. This thesis describes new methods to apply the triple oxygen isotope tracer of photosynthetic production in a salt marsh. Additionally, noble gases are applied to constrain air-water exchange processes which affect metabolism tracers. These stable, natural abundance tracers complement traditional techniques for measuring metabolism. In particular, they highlight the potential importance of daytime oxygen sinks besides aerobic respiration, such as rising bubbles. In tidal creeks, increasing nutrients may increase both production and respiration, without any apparent change in the net metabolism. In ponds, daytime production and respiration are also tightly coupled, but there is high background respiration regardless of changes in daytime production. Both tidal creeks and ponds have higher respiration rates and lower production rates than the marsh platform, suggesting that expansion of these submerged environments could limit the ability of salt marshes to sequester carbon.

Description: Financial support for my doctoral research was provided by the United States Department of Defense through the National Defense Science and Engineering Graduate Fellowship Program, the National Science Foundation under grant OCE-1233678, and the Woods Hole Oceanographic Institution (WHOI) under grants from the WHOI Coastal Ocean Institute, Ocean and Climate Change Institute, and Ocean Life Institute. WHOI Academic Programs Office also provided funding support for research, through the Ocean Ventures Fund, and for my stipend, as graduate research assistantships including an assistantship from the United States Geological Survey administered by WHOI.

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

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

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

Description: 2017-02-19

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

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

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

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

Description: The behavior of an axisymmetric vertical turbulent jet in an unconfined stratified environment is studied by means of well-resolved, large-eddy simulations. The stratification is two uniform layers separated by a thermocline. This study considers two cases: when the thermocline thickness is small and on the order of the jet diameter at the thermocline entrance. The Froude number of the jet at the thermocline varies from 0.6 to 1.9, corresponding to the class of weak fountains. The mean jet penetration, stratified turbulent entrainment, jet oscillations, and the generation of internal waves are examined. The mean jet penetration is predicted well by a simple model based on the conservation of the source energy in the thermocline. The entrainment coefficient for the thin thermocline is consistent with the theoretical model for a two-layer stratification with a sharp interface, while for the thick thermocline entrainment is larger at low Froude numbers. The data reveal the presence of a secondary horizontal flow in the upper part of the thick thermocline, resulting in the entrainment of fluid from the thermocline rather than from the upper stratification layer. The spectra of the jet oscillations in the thermocline display two peaks, at the same frequencies for both stratifications at fixed Froude number. For the thick thermocline, internal waves are generated only at the lower frequency, since the higher peak exceeds the maximal buoyancy frequency. For the thin thermocline, conversely, the spectra of the internal waves show the two peaks at low Froude numbers, whereas only one peak at the lower frequency is observed at higher Froude numbers.

Description: This work was supported by the Linné FLOW Centre at KTH (E. E.), the European Research Council Grant ERC-2013-CoG-616186, TRITOS (L. B.), and the Swedish Research Council (VR), Outstanding Young Researcher Award (L. B.). Support to C. C. was given by the NSF Project OCE-1434041.

Description: 2017-05-10

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

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

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

Description: 2017-04-20

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

Description: Direct measurements of oceanic turbulent parameters were taken upstream of and across Drake Passage, in the region of the Subantarctic and Polar Fronts. Values of turbulent kinetic energy dissipation rate ε estimated by microstructure are up to two orders of magnitude lower than previously published estimates in the upper 1000 m. Turbulence levels in Drake Passage are systematically higher than values upstream, regardless of season. The dissipation of thermal variance χ is enhanced at middepth throughout the surveys, with the highest values found in northern Drake Passage, where water mass variability is the most pronounced. Using the density ratio, evidence for double-diffusive instability is presented. Subject to double-diffusive physics, the estimates of diffusivity using the Osborn–Cox method are larger than ensemble statistics based on ε and the buoyancy frequency.

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

Description: 2016-10-05

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): 1769-1783, doi:10.1175/JPO-D-15-0193.1.

Description: High-resolution observations of velocity, salinity, and turbulence quantities were collected in a salt wedge estuary to quantify the efficiency of stratified mixing in a high-energy environment. During the ebb tide, a midwater column layer of strong shear and stratification developed, exhibiting near-critical gradient Richardson numbers and turbulent kinetic energy (TKE) dissipation rates greater than 10−4 m2 s−3, based on inertial subrange spectra. Collocated estimates of scalar variance dissipation from microconductivity sensors were used to estimate buoyancy flux and the flux Richardson number Rif. The majority of the samples were outside the boundary layer, based on the ratio of Ozmidov and boundary length scales, and had a mean Rif = 0.23 ± 0.01 (dissipation flux coefficient Γ = 0.30 ± 0.02) and a median gradient Richardson number Rig = 0.25. The boundary-influenced subset of the data had decreased efficiency, with Rif = 0.17 ± 0.02 (Γ = 0.20 ± 0.03) and median Rig = 0.16. The relationship between Rif and Rig was consistent with a turbulent Prandtl number of 1. Acoustic backscatter imagery revealed coherent braids in the mixing layer during the early ebb and a transition to more homogeneous turbulence in the midebb. A temporal trend in efficiency was also visible, with higher efficiency in the early ebb and lower efficiency in the late ebb when the bottom boundary layer had greater influence on the flow. These findings show that mixing efficiency of turbulence in a continuously forced, energetic, free shear layer can be significantly greater than the broadly cited upper bound from Osborn of 0.15–0.17.

Description: Holleman was supported by the Devonshire Scholars program. The field study and the coauthors’ contributions were supported by NSF Grant OCE 0926427.

Description: 2016-11-24

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, 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): 1823-1837, doi:10.1175/JPO-D-15-0165.1.

Description: Measurements just beneath the ocean surface demonstrate that the primary mechanism by which energy from breaking waves is transmitted into the water column is through the work done by the covariance of turbulent pressure and velocity fluctuations. The convergence in the vertical transport of turbulent kinetic energy (TKE) balances the dissipation rate of TKE at first order and is nearly an order of magnitude greater than the sum of the integrated Eulerian and Stokes shear production. The measured TKE transport is consistent with a simple conceptual model that assumes roughly half of the surface flux of TKE by wave breaking is transmitted to depths greater than the significant wave height. During conditions when breaking waves are inferred, the direction of momentum flux is more aligned with the direction of wave propagation than with the wind direction. Both the energy and momentum fluxes occur at frequencies much lower than the wave band, consistent with the time scales associated with wave breaking. The largest instantaneous values of momentum flux are associated with strong downward vertical velocity perturbations, in contrast to the pressure work, which is associated with strong drops in pressure and upward vertical velocity perturbations.

Description: Funding for this research was provided by the National Science Foundation Grants OCE-1339032 and OCE-1338518

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): 2006–2024, doi:10.1175/JPO-D-14-0234.1.

Description: The effects of wind-driven whitecapping on the evolution of the ocean surface boundary layer are examined using an idealized one-dimensional Reynolds-averaged Navier–Stokes numerical model. Whitecapping is parameterized as a flux of turbulent kinetic energy through the sea surface and through an adjustment of the turbulent length scale. Simulations begin with a two-layer configuration and use a wind that ramps to a steady stress. This study finds that the boundary layer begins to thicken sooner in simulations with whitecapping than without because whitecapping introduces energy to the base of the boundary layer sooner than shear production does. Even in the presence of whitecapping, shear production becomes important for several hours, but then inertial oscillations cause shear production and whitecapping to alternate as the dominant energy sources for mixing. Details of these results are sensitive to initial and forcing conditions, particularly to the turbulent length scale imposed by breaking waves and the transfer velocity of energy from waves to turbulence. After 1–2 days of steady wind, the boundary layer in whitecapping simulations has thickened more than the boundary layer in simulations without whitecapping by about 10%–50%, depending on the forcing and initial conditions.

Description: We thank Skidmore College for financial and infrastructure support, and Skidmore and the National Science Foundation for funding travel to meetings where early versions of this work were presented. We also thank the National Science Foundation, Oregon State University, Jonathan Nash, and Joe Jurisa for funding and hosting a workshop on River Plume Mixing in October, 2013, where ideas and context for this paper were developed.

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

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

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

Description: 2015-08-01

Description: Author Posting. © American Meteorological Society, 2015. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 45 (2015): 966–987, doi:10.1175/JPO-D-14-0110.1.

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

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

Description: 2015-10-01

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

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

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

Description: 2016-05-01

Description: Author Posting. © American Meteorological Society, 2015. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 45 (2015): 1822–1842, doi:10.1175/JPO-D-14-0147.1.

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

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

Description: 2016-01-01

Description: Author Posting. © American Meteorological Society, 2015. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 45 (2015): 1735–1756, doi:10.1175/JPO-D-14-0238.1.

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

Description: The work was supported by NSF OCE 0850416.

Description: 2015-12-01

Description: 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): 2497–2521, doi:10.1175/JPO-D-14-0128.1.

Description: Oceanic density overturns are commonly used to parameterize the dissipation rate of turbulent kinetic energy. This method assumes a linear scaling between the Thorpe length scale LT and the Ozmidov length scale LO. Historic evidence supporting LT ~ LO has been shown for relatively weak shear-driven turbulence of the thermocline; however, little support for the method exists in regions of turbulence driven by the convective collapse of topographically influenced overturns that are large by open-ocean standards. This study presents a direct comparison of LT and LO, using vertical profiles of temperature and microstructure shear collected in the Luzon Strait—a site characterized by topographically influenced overturns up to O(100) m in scale. The comparison is also done for open-ocean sites in the Brazil basin and North Atlantic where overturns are generally smaller and due to different processes. A key result is that LT/LO increases with overturn size in a fashion similar to that observed in numerical studies of Kelvin–Helmholtz (K–H) instabilities for all sites but is most clear in data from the Luzon Strait. Resultant bias in parameterized dissipation is mitigated by ensemble averaging; however, a positive bias appears when instantaneous observations are depth and time integrated. For a series of profiles taken during a spring tidal period in the Luzon Strait, the integrated value is nearly an order of magnitude larger than that based on the microstructure observations. Physical arguments supporting LT ~ LO are revisited, and conceptual regimes explaining the relationship between LT/LO and a nondimensional overturn size are proposed. In a companion paper, Scotti obtains similar conclusions from energetics arguments and simulations.

Description: B.D.M. and S.K.V. gratefully acknowledge the support of the Office of Naval Research under Grants N00014-12-1-0279, N00014-12-1-0282, and N00014-12-1-0938 (Program Manager: Dr. Terri Paluszkiewicz). S.K.V. also acknowledges support of the National Science Foundation under Grant OCE-1151838. L.S.L. acknowledges support for BBTRE by the National Science Foundation by Contract OCE94-15589 and NATRE and IWISE by the Office of Naval Research by Contracts N00014-92-1323 and N00014-10-10315. J.N.M. was supported through Grant 1256620 from the National Science Foundation and the Office of Naval Research (IWISE Project).

Description: 2016-04-01

Description: Author Posting. © American Meteorological Society, 2015. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 45 (2015): 2621–2639, doi:10.1175/JPO-D-14-0239.1.

Description: Measurements made as part of a large-scale experiment to examine wind-driven circulation and mixing in Chesapeake Bay demonstrate that circulations consistent with Langmuir circulation play an important role in surface boundary layer dynamics. Under conditions when the turbulent Langmuir number Lat is low (〈0.5), the surface mixed layer is characterized by 1) elevated vertical turbulent kinetic energy; 2) decreased anisotropy; 3) negative vertical velocity skewness indicative of strong/narrow downwelling and weak/broad upwelling; and 4) strong negative correlations between low-frequency vertical velocity and the velocity in the direction of wave propagation. These characteristics appear to be primarily the result of the vortex force associated with the surface wave field, but convection driven by a destabilizing heat flux is observed and appears to contribute significantly to the observed negative vertical velocity skewness. Conditions that favor convection usually also have strong Langmuir forcing, and these two processes probably both contribute to the surface mixed layer turbulence. Conditions in which traditional stress-driven turbulence is important are limited in this dataset. Unlike other shallow coastal systems where full water column Langmuir circulation has been observed, the salinity stratification in Chesapeake Bay is nearly always strong enough to prevent full-depth circulation from developing.

Description: The funding for this research was provided by the National Science Foundation Grants OCE-1339032 and OCE-1338518.

Description: 2016-04-01

Description: Author Posting. © American Meteorological Society, 2015. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 45 (2015): 778–791, doi:10.1175/JPO-D-14-0164.1.

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

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

Description: 2015-09-01

Description: Author Posting. © American Meteorological Society, 2015. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 45 (2015): 294–312, doi:10.1175/JPO-D-14-0104.1.

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

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

Description: 2015-07-01

Description: Author Posting. © American Meteorological Society, 2015. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 45 (2015): 104–132, doi:10.1175/JPO-D-14-0032.1.

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

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

Description: 2015-07-01

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

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

Description: Published

Description: 1618/1637

Description: 4A. Clima e Oceani

Description: JCR Journal

Description: restricted

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

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

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

Description: 2014-06-15

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

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

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

Description: 2014-07-01

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

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

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

Description: 2015-04-01

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

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

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

Description: 2014-10-15

Description: Author Posting. © American Meteorological Society, 2014. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Atmospheric and Oceanic Technology 31 (2014): 1410–1421, doi:10.1175/JTECH-D-13-00230.1.

Description: Aerial images are used to quantify the concentration of fluorescent Rhodamine water tracing (WT) dye in turbid and optically deep water. Tracer releases near the shoreline of an ocean beach and near a tidal inlet were observed with a two-band multispectral camera and a pushbroom hyperspectral imager, respectively. The aerial observations are compared with near-surface in situ measurements. The ratio of upwelling radiance near the Rhodamine WT excitation and emission peaks varies linearly with the in situ dye concentrations for concentrations 〈20 ppb (r2 = 0.70 and r2 = 0.85–0.88 at the beach and inlet, respectively). The linear relationship allows for relative tracer concentration estimates without in situ calibration. The O(1 m) image pixels resolve complex flow structures on the inner shelf that transport and mix tracer.

Description: We thank ONR and NSF for funding this work.

Description: 2014-12-01

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

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

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

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

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

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

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

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

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

Description: 2013-08-01

Description: Author Posting. © American Meteorological Society, 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 Atmospheric and Oceanic Technology 30 (2013): 1767–1788, doi:10.1175/JTECH-D-12-00140.1.

Description: Seismic images of oceanic thermohaline finestructure record vertical displacements from internal waves and turbulence over large sections at unprecedented horizontal resolution. Where reflections follow isopycnals, their displacements can be used to estimate levels of turbulence dissipation, by applying the Klymak–Moum slope spectrum method. However, many issues must be considered when using seismic images for estimating turbulence dissipation, especially sources of random and harmonic noise. This study examines the utility of seismic images for estimating turbulence dissipation in the ocean, using synthetic modeling and data from two field surveys, from the South China Sea and the eastern Pacific Ocean, including the first comparison of turbulence estimates from seismic images and from vertical shear. Realistic synthetic models that mimic the spectral characteristics of internal waves and turbulence show that reflector slope spectra accurately reproduce isopycnal slope spectra out to horizontal wavenumbers of 0.04 cpm, corresponding to horizontal wavelengths of 25 m. Using seismic reflector slope spectra requires recognition and suppression of shot-generated harmonic noise and restriction of data to frequency bands with signal-to-noise ratios greater than about 4. Calculation of slope spectra directly from Fourier transforms of the seismic data is necessary to determine the suitability of a particular dataset to turbulence estimation from reflector slope spectra. Turbulence dissipation estimated from seismic reflector displacements compares well to those from 10-m shear determined by coincident expendable current profiler (XCP) data, demonstrating that seismic images can produce reliable estimates of turbulence dissipation in the ocean, provided that random noise is minimal and harmonic noise is removed.

Description: This work was funded by NSF Grants 0452744, 0405654, and 0648620, and ONR/DEPSCoR Grant DODONR40027.

Description: 2014-02-01

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

Description: Owing to the larger thermal expansion coefficient at higher temperatures, more buoyancy is put into the ocean by heating than is removed by cooling at low temperatures. The authors show that, even with globally balanced thermal and haline surface forcing at the ocean surface, there is a negative density flux and hence a positive buoyancy flux. As shown by McDougall and Garrett, this must be compensated by interior densification on mixing due to the nonlinearity of the equation of state (cabbeling). Three issues that arise from this are addressed: the estimation of the annual input of density forcing, the effects of the seasonal cycle, and the total cabbeling potential of the ocean upon complete mixing. The annual expansion through surface density forcing in a steady-state ocean driven by balanced evaporation–precipitation–runoff (E–P–R) and net radiative budget at the surface Qnet is estimated as 74 000 m3 s−1 (0.07 Sv; 1 Sv ≡ 106 m3 s−1), which would be equivalent to a sea level rise of 6.3 mm yr−1. This is equivalent to approximately 3 times the estimated rate of sea level rise or 450% of the average Mississippi River discharge. When seasonal variations are included, this density forcing increases by 35% relative to the time-mean case to 101 000 m3 s−1 (0.1 Sv). Likely bounds are established on these numbers by using different Qnet and E–P–R datasets and the estimates are found to be robust to a factor of ~2. These values compare well with the cabbeling-induced contraction inferred from independent thermal dissipation rate estimates. The potential sea level decrease upon complete vertical mixing of the ocean is estimated as 230 mm. When horizontal mixing is included, the sea level drop is estimated as 300 mm.

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

Description: 2013-10-01

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

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

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

Description: 2013-06-01

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

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

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

Description: 2013-11-01

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

Description: The Community Climate System Model, version 4 (CCSM4) is used to assess the climate impact of wind-generated near-inertial waves (NIWs). Even with high-frequency coupling, CCSM4 underestimates the strength of NIWs, so that a parameterization for NIWs is developed and included into CCSM4. Numerous assumptions enter this parameterization, the core of which is that the NIW velocity signal is detected during the model integration, and amplified in the shear computation of the ocean surface boundary layer module. It is found that NIWs deepen the ocean mixed layer by up to 30%, but they contribute little to the ventilation and mixing of the ocean below the thermocline. However, the deepening of the tropical mixed layer by NIWs leads to a change in tropical sea surface temperature and precipitation. Atmospheric teleconnections then change the global sea level pressure fields so that the midlatitude westerlies become weaker. Unfortunately, the magnitude of the real air-sea flux of NIW energy is poorly constrained by observations; this makes the quantitative assessment of their climate impact rather uncertain. Thus, a major result of the present study is that because of its importance for global climate the uncertainty in the observed tropical NIW energy has to be reduced.

Description: This research was funded as part of the Climate Process Team on internal wave-driven mixing with NSF Grant Nr E0968771 at NCAR.

Description: 2013-11-01

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

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

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

Description: 2014-02-01

Description: 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): 1841–1861, doi:10.1175/JPO-D-12-0231.1.

Description: In this idealized numerical modeling study, the composition of residual sediment fluxes in energetic (e.g., weakly or periodically stratified) tidal estuaries is investigated by means of one-dimensional water column models, with some focus on the sediment availability. Scaling of the underlying dynamic equations shows dependence of the results on the Simpson number (relative strength of horizontal density gradient) and the Rouse number (relative settling velocity) as well as impacts of the Unsteadiness number (relative tidal frequency). Here, the parameter space given by the Simpson and Rouse numbers is mainly investigated. A simple analytical model based on the assumption of stationarity shows that for small Simpson and Rouse numbers sediment flux is down estuary and vice versa for large Simpson and Rouse numbers. A fully dynamic water column model coupled to a second-moment turbulence closure model allows to decompose the sediment flux profiles into contributions from the transport flux (product of subtidal velocity and sediment concentration profiles) and the fluctuation flux profiles (tidal covariance between current velocity and sediment concentration). Three different types of bottom sediment pools are distinguished to vary the sediment availability, by defining a time scale for complete sediment erosion. For short erosion times scales, the transport sediment flux may dominate, but for larger erosion time scales the fluctuation sediment flux largely dominates the tidal sediment flux. When quarter-diurnal components are added to the tidal forcing, up-estuary sediment fluxes are strongly increased for stronger and shorter flood tides and vice versa. The theoretical results are compared to field observations in a tidally energetic inlet.

Description: Project funding was provided by the German Research Foundation (DFG) in the framework of the Project ECOWS (Role of Estuarine Circulation for Transport of Suspended Particulate Matter in the Wadden Sea, BU 1199/11) and by the German Federal Ministry of Research and Education in the framework of the Project PACE [The future of the Wadden Sea sediment fluxes: still keeping pace with sea level rise? (FKZ 03F0634A)].

Description: 2014-03-01

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

Description: Pacific Water flows across the shallow Chukchi Sea before reaching the Arctic Ocean, where it is a source of heat, freshwater, nutrients, and carbon. A substantial portion of Pacific Water is routed through Barrow Canyon, located in the northeast corner of the Chukchi. Barrow Canyon is a region of complex geometry and forcing where a variety of water masses have been observed to coexist. These factors contribute to a dynamic physical environment, with the potential for significant water mass transformation. The measurements of turbulent kinetic energy dissipation presented here indicate diapycnal mixing is important in the upper canyon. Elevated dissipation rates were observed near the pycnocline, effectively mixing winter and summer water masses, as well as within the bottom boundary layer. The slopes of shear/stratification layers, combined with analysis of rotary spectra, suggest that near-inertial wave activity may be important in modulating dissipation near the bottom. Because the canyon is known to be a hotspot of productivity with an active benthic community, mixing may be an important factor in maintenance of the biological environment.

Description: ELS was supported as a WHOI Postdoctoral Scholar through the WHOI Ocean and Climate Change Institute.

Description: 2012-12-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): 2168–2186, doi:10.1175/JPO-D-11-08.1.

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

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

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

Description: Ventilation, including subduction and obduction, for the global oceans was examined using Simple Ocean Data Assimilation (SODA) outputs. The global subduction rate averaged over the period from 1959 to 2006 is estimated at 505.8 Sv (1 Sv ≡ 106 m3 s−1), while the corresponding global obduction rate is estimated at 482.1 Sv. The annual subduction/obduction rates vary greatly on the interannual and decadal time scales. The global subduction rate is estimated to have increased 7.6% over the past 50 years, while the obduction rate is estimated to have increased 9.8%. Such trends may be insignificant because errors associated with the data generated by ocean data assimilation could be as large as 10%. However, a major physical mechanism that induced these trends is primarily linked to changes in the Southern Ocean. While the Southern Ocean plays a key role in global subduction and obduction rates and their variability, both the Southern Ocean and equatorial regions are critically important sites of water mass formation/erosion.

Description: This work was supported by the Key State Basic Research Program of China under Grant 2012CB417401, the National Natural Science Foundation of China (Grants 40906007, 40890152), and the Open Foundation of Physical Oceanography Laboratory, OUC, under Grant 200902.

Description: 2012-08-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 May 1998

Description: Planktonic protozoan grazers have the potential to significantly affect the chemistry of particle-associated trace metals. This is due both to the importance of protists as consumers of bacterial-sized particles, and to the unique low-pH, enzyme-rich microenvironment of the grazer food vacuole. This thesis examines the role of protozoan grazers in the marine geochemistry of strongly hydrolyzed, particle-reactive trace metals, in particular Th and Fe. A series of tracer experiments was carried out in model systems in order to determine the effect of grazer-mediated transformations on the chemical speciation and partitioning of radioisotopes C9Fe, 234Th, 51Cr) associated with prey cells. Results indicate that protozoan grazers are equally able to mobilize intracellular and extracellular trace metals. In some cases, protozoan regeneration of trace metals appears to lead to the formation of metal-organic complexes. Protozoan grazing may generate colloidal material that can scavenge trace metals and, via aggregation, lead to an increase in the metal/organic carbon ratio of aggregated particles. Model system experiments were also conducted in order to determine the effect of grazers on mineral phases, specifically colloidal iron oxide (ferrihydrite). Several independent techniques were employed, including size fractionation ors9Fe-labeled colloids, competitive ligand exchange, and iron-limited diatoms as "probes" for bioavailable Fe. Experimental evidence strongly suggests that protozoan grazing can affect the surface chemistry and increase the dissolution rate of iron oxide phases through phagotrophic ingestion. In further work on protozoan-mediated dissolution of colloidal Fe oxides, a novel tracer technique was developed based on the synthesis of colloidal ferrihydrite impregnated with 133Ba as an inert tracer. This technique was shown to be a sensitive, quantitative indicator for the extent of ferrihydrite dissolution/alteration by a variety of mechanisms, including photochemical reduction and ligand-mediated dissolution. In field experiments using this technique, grazing by naturally occuring protistan assemblages was shown to significantly enhance the dissolution rate of colloidal ferrihydrite over that in non-grazing controls. Laboratory and field results indicate that, when integrated temporally over the entire euphotic zone, protozoan grazing may equal or exceed photoreduction as a pathway for the dissolution of iron oxides.

Description: This work was financially supported by a Department of Defense ONR-NDSEG Graduate Fellowship, Office ofNaval Research AASERT Award (N00014-94-1-0711), and the National Science Foundation EGB Program (OCE-9523910).

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

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

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

Description: 2012-11-01

Description: 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 Doctor of Philosophy at the Massachusetts Institute of Technology and Woods Hole Oceanographic Institution February 1994

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: 2013-04-01

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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, 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): 2910-2925, doi:10.1175/2009JPO4139.1.

Description: The propagation of Rossby waves on a midlatitude β plane is investigated in the presence of density diffusion with the aid of linear hydrostatic theory. The search for wave solutions in a vertically bounded medium subject to horizontal (vertical) diffusion leads to an eigenvalue problem of second (fourth) order. Exact solutions of the problem are obtained for uniform background stratification (N), and approximate solutions are constructed for variable N using the Wentzel–Kramers–Brillouin method. Roots of the eigenvalue relations for free waves are found and discussed. The barotropic wave of adiabatic theory is also a solution of the eigenvalue problem as this is augmented with density diffusion in the horizontal or vertical direction. The barotropic wave is undamped as fluid parcels in the wave move only horizontally and are therefore insensitive to the vortex stretching induced by mixing. On the other hand, density diffusion modifies the properties of baroclinic waves of adiabatic theory. In the presence of horizontal diffusion the baroclinic modes are damped but their vertical structure remains unaltered. The ability of horizontal diffusion to damp baroclinic waves stems from its tendency to counteract the deformation of isopycnal surfaces caused by the passage of these waves. The damping rate increases (i) linearly with horizontal diffusivity and (ii) nonlinearly with horizontal wavenumber and mode number. In the presence of vertical diffusion the baroclinic waves suffer both damping and a change in vertical structure. In the long-wave limit the damping is critical (wave decay rate numerically equal to wave frequency) and increases as the square roots of vertical diffusivity and zonal wavenumber. Density diffusion in the horizontal or vertical direction reduces the amplitude of the phase speed of westward-propagating waves. Observational estimates of eddy diffusivities suggest that horizontal and vertical mixing strongly attenuates baroclinic waves in the ocean but that vertical mixing is too weak to notably modify the vertical structure of the gravest modes.

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

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

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

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

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

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

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

Description: 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): 909–917, doi:10.1175/2007JPO3535.1.

Description: The classical two-box model of Stommel is extended in two directions: replacing the buoyancy constraint with an energy constraint and including the wind-driven gyre. Stommel postulated a buoyancy constraint for the thermohaline circulation, and his basic idea has evolved into the dominating theory of thermohaline circulation; however, recently, it is argued that the thermohaline circulation is maintained by mechanical energy from wind stress and tides. The major difference between these two types of models is the bifurcation structure: the Stommel-like model has two thermal modes (one stable and another one unstable) and one stable haline mode, whereas the energy-constraint model has one stable thermal mode and two saline modes (one stable and another one unstable). Adding the wind-driven gyre changes the threshold value of thermohaline bifurcation greatly; thus, the inclusion of the wind-driven gyre is a vital step in completely modeling the physical processes related to thermohaline circulation.

Description: YPG was supported by the National Science Foundation of China (NSFC, 40676022), the National Basic Research Program of China (2006CB403605), and the Guangdong Natural Science Foundation (5003672). RXH was supported by the National Oceanic and Atmospheric Administration through CICOR Cooperative Agreement NA17RJ1223 to the Woods Hole Oceanographic Institution.

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: A series of dye releases in the Hudson River estuary elucidated diapycnal mixing rates and temporal variability over tidal and fortnightly time scales. Dye was injected in the bottom boundary layer for each of four releases during different phases of the tide and of the spring–neap cycle. Diapycnal mixing occurs primarily through entrainment that is driven by shear production in the bottom boundary layer. On flood the dye extended vertically through the bottom mixed layer, and its concentration decreased abruptly near the base of the pycnocline, usually at a height corresponding to a velocity maximum. Boundary layer growth is consistent with a one-dimensional, stress-driven entrainment model. A model was developed for the vertical structure of the vertical eddy viscosity in the flood tide boundary layer that is proportional to u2*/N∞, where u* and N∞ are the bottom friction velocity and buoyancy frequency above the boundary layer. The model also predicts that the buoyancy flux averaged over the bottom boundary layer is equal to 0.06N∞u2* or, based on the structure of the boundary layer equal to 0.1NBLu2*, where NBL is the buoyancy frequency across the flood-tide boundary layer. Estimates of shear production and buoyancy flux indicate that the flux Richardson number in the flood-tide boundary layer is 0.1–0.18, consistent with the model indicating that the flux Richardson number is between 0.1 and 0.14. During ebb, the boundary layer was more stratified, and its vertical extent was not as sharply delineated as in the flood. During neap tide the rate of mixing during ebb was significantly weaker than on flood, owing to reduced bottom stress and stabilization by stratification. As tidal amplitude increased ebb mixing increased and more closely resembled the boundary layer entrainment process observed during the flood. Tidal straining modestly increased the entrainment rate during the flood, and it restratified the boundary layer and inhibited mixing during the ebb.

Description: The work was supported by the National Science Foundation Grant OCE00-95972 (W. Geyer, J. Lerczak), OCE00-99310 (R. Houghton), and OCE00-95913 (R. Chant, E. Hunter).

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): 1551-1573, doi:10.1175/2008JPO4152.1.

Description: A conceptually simple model is presented for predicting the amplitude and periodicity of eddies generated by a steady poleward outflow in a 1½-layer β-plane formulation. The prediction model is rooted in linear quasigeostrophic dynamics but is capable of predicting the amplitude of the β plume generated by outflows in the nonlinear range. Oscillations in the plume amplitude are seen to represent a near-zero group velocity response to an adjustment process that can be traced back to linear dynamics. When the plume-amplitude oscillations become large enough so that the coherent β plume is replaced by a robust eddy field, the eddy amplitude is still constrained by the plume-amplitude prediction model. The eddy periodicity remains close to that of the predictable, near-zero group-velocity linear oscillations. Striking similarities between the patterns of variability in the model and observations south of Indonesia’s Lombok Strait suggest that the processes investigated in this study may play an important role in the generation of the observed eddy field of the Indo-Australian Basin.

Description: This work was completed at the Woods Hole Oceanographic Institution while TS Durland was supported by the Ocean and Climate Change Institute. MA Spall was supported by NSF Grant OCE-0423975 and J Pedlosky by NSF Grant OCE-0451086. TS Durland acknowledges additional report preparation support from NASA Grant NNG05GN98G.

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): 418-434, doi:10.1175/2007JPO3372.1.

Description: Stratification and turbulent mixing exhibit a flood–ebb tidal asymmetry in estuaries and continental shelf regions affected by horizontal density gradients. The authors use a large-eddy simulation (LES) model to investigate the penetration of a tidally driven bottom boundary layer into stratified water in the presence of a horizontal density gradient. Turbulence in the bottom boundary layer is driven by bottom stress during flood tides, with low-gradient (Ri) and flux (Rf) Richardson numbers, but by localized shear during ebb tides, with Ri = ¼ and Rf = 0.2 in the upper half of the boundary layer. If the water column is unstratified initially, the LES model reproduces periodic stratification associated with tidal straining. The model results show that the energetics criterion based on the competition between tidal straining and tidal stirring provides a good prediction for the onset of periodic stratification, but the tidally averaged horizontal Richardson number Rix has a threshold value of about 0.2, which is lower than the 3 suggested in a recent study. Although the tidal straining leads to negative buoyancy flux on flood tides, the authors find that for typical values of the horizontal density gradient and tidal currents in estuaries and shelf regions, buoyancy production is much smaller than shear production in generating turbulent kinetic energy.

Description: This work is supported by Grants OCE-0451699 and OCE-0451740 from the National Science Foundation.

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 2000

Description: The plausibility of local baroclinic instability as a generation mechanism for midocean mesoscale eddies is examined with a two-layer, quasi-geostrophic (QG) model forced by an imposed, horizontally homogeneous, vertically sheared mean flow and dissipated through bottom Ekman friction, Explanations are sought for two observed features of mid-ocean eddies: 1) substantial energy is retained in the baroclinic mode and in the associated deformation radius (Rd) scale, and 2) the ratio of eddy to mean kinetic energy is much larger than one, The tendency of QG to cascade energy into the barotropic mode and into scales larger than Rd can be counteracted when stratification is surface-trapped, for then the baroclinic mode is weakly damped, and hence enhanced, Numerical experiments are performed with both surface-trapped and uniform stratification to quantify this, Experiments with equal Ekman frictions in the two layers are also performed for purposes of contrast, Interpretation is aided with an inequality derived from the energy and enstrophy equations, The inequality forbids the simultaneous retention of substantial energy in the baroclinic mode and in scales near Rd when Ekman friction is symmetric, but points towards surface-trapped stratification and bottomtrapped friction as an environment in which both of these can be achieved, The dissertation also contains a systematic study of geostrophic turbulence forced by nonzonal flows, Narrow zonal jets emerge when shear-induced mean potential vorticity (PV) gradients are small compared to the planetary gradient (β), and energy is a strong function of the angle shear presents to the east-west direction, When shear-induced PV gradients are comparable to β, and the mean shear has a westward component, fields of monopolar vortices form and persist, Energy is asymmetric between fields of cyclones and anticyclones, Such asymmetry was commonly thought not to occur in QG, but is shown here to be introduced by the nonzonal basic state, In both jet and vortex regimes, eddy energy can be much larger than mean kinetic energy, contrary to the expectation that β stabilizes weak shear flows,

Description: My first three years here were funded by an Office of Naval Research/National Defense Science and Engineering Graduate Fellowship administered by Jeff Jarocz at the American Society for Engineering Education. During the last three years, most of my support has come the National Science Foundation via grant OCE-9617848, with some additional support coming from the Office of Naval Research via grant N00014-95-1-0824.

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

Description: The circulation in the equatorial Pacific Ocean is studied in a series of numerical experiments based on an isopycnal coordinate model. The model is subject to monthly mean climatology of wind stress and surface thermohaline forcing. In response to decadal variability in the diapycnal mixing coefficient, sea surface temperature and other properties of the circulation system oscillate periodically. The strongest sea surface temperature anomaly appears in the geographic location of Niño-3 region with the amplitude on the order of 0.5°C, if the model is subject to a 30-yr sinusoidal oscillation in diapycnal mixing coefficient that varies between 0.03 × 10−4 and 0.27 × 10−4 m2 s−1. Changes in diapycnal mixing coefficient of this amplitude are within the bulk range consistent with the external mechanical energy input in the global ocean, especially when considering the great changes of tropical cyclones during the past decades. Thus, time-varying diapycnal mixing associated with changes in wind energy input into the ocean may play a nonnegligible role in decadal climate variability in the equatorial circulation and climate.

Description: CJH and WW were supported by The National Natural Science Foundation of China through Grant 40476010 and National Basic Research Priorities Programmer of China through Grant 2005CB422302. RXH was supported by the National Oceanic and Atmospheric Administration through CICOR Cooperative Agreement NA17RJ1223 to the Woods Hole Oceanographic Institution. This study is also supported through the Chinese 111 Project under Contract B07036.

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): 1253-1266, doi:10.1175/2007JPO3786.1.

Description: Wind stress and tidal dissipation are the most important sources of mechanical energy for maintaining the oceanic general circulation. The contribution of mechanical energy due to tropical cyclones can be a vitally important factor in regulating the oceanic general circulation and its variability. However, previous estimates of wind stress energy input were based on low-resolution wind stress data in which strong nonlinear events, such as tropical cyclones, were smoothed out. Using a hurricane–ocean coupled model constructed from an axisymmetric hurricane model and a three-layer ocean model, the rate of energy input to the world’s oceans induced by tropical cyclones over the period from 1984 to 2003 was estimated. The energy input is estimated as follows: 1.62 TW to the surface waves and 0.10 TW to the surface currents (including 0.03 TW to the near-inertial motions). The rate of gravitational potential energy increase due to tropical cyclones is 0.05 TW. Both the energy input from tropical cyclones and the increase of gravitational potential energy of the ocean show strong interannual and decadal variability with an increasing rate of 16% over the past 20 years. The annual mean diapycnal upwelling induced by tropical cyclones over the past 20 years is estimated as 39 Sv (Sv ≡ 106 m3 s−1). Owing to tropical cyclones, diapycnal mixing in the upper ocean (below the mixed layer) is greatly enhanced. Within the regimes of strong activity of tropical cyclones, the increase of diapycnal diffusivity is on the order of (1 − 6) × 10−4 m2 s−1. The tropical cyclone–related energy input and diapycnal mixing may play an important role in climate variability, ecology, fishery, and environments.

Description: LLL and WW were supported by the National Basic Research Priorities Programmer of China through Grant 2007CB816004 and National Outstanding Youth Natural Science Foundation of China FIG. 15. Annual-mean vertical diffusivity induced by tropical cyclones from 1984 to 2003 (units: 10 4 m2 s 1): (right) the horizontal distribution and (left) the zonally averaged vertical diffusivity. JUNE 2008 L IU ET AL . 1265 under Grant 40725017. RXH was supported by the W. Alan Clark Chair from Woods Hole Oceanographic Institution.

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): 2341–2347, doi:10.1175/2010JPO4465.1.

Description: The mean downwelling in an eddy-resolving model of a convective basin is concentrated near the boundary where eddies are shed from the cyclonic boundary current into the interior. It is suggested that the buoyancy-forced downwelling in the Labrador Sea and the Lofoten Basin is similarly concentrated in analogous eddy formation regions along their eastern boundaries. Use of a transformed Eulerian mean depiction of the density transport reveals the central role eddy fluxes play in maintaining the adiabatic nature of the flow in a nonperiodic region where heat is lost from the boundary current. The vorticity balance in the downwelling region is primarily between stretching of planetary vorticity and eddy flux divergence of relative vorticity, although a narrow viscous boundary layer is ultimately important in closing the regional vorticity budget. This overall balance is similar in some ways to the diffusive–viscous balance represented in previous boundary layer theories, and suggests that the downwelling in convective basins may be properly represented in low-resolution climate models if eddy flux parameterizations are adiabatic, identify localized regions of eddy formations, and allow density to be transported far from the region of eddy formations.

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

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

Description: With Lagrangian and hydrographic data taken in the deep Brazil Basin, we identify several submesoscale coherent vortices (SCVs). These features contrast with SCV paradigms in that float data indicate approximately equal populations of cyclonic and anticyclonic vortices, and hydrographic data suggest that roughly half exhibit the convex lens shape generally associated with SCVs, while half are instead shaped like a concave lens, with isopycnal surfaces pinched together. There is some evidence that the vortex cores may be enriched in warm, salty North Atlantic Deep Water, suggesting formation in the north or northwest regions of the basin. Data is available from 153 floats which were ballasted for 2500 and 4000 db pressures. They tracked 34 eddies, which are believed to be roughly 30 km in diameter, and rotate with apparent periods of about 30 days. Many floats experienced formation or entrainment events, and destruction or detrainment events, near seamounts.

Description: This work was supported by the National Science Foundation, grants OCE 99-11148 and OCE 94-15509.