ALBERT

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

Description: Inertial terms dominate the single-gyre ocean model and prevent western-intensification when the viscosity is small. This occurs long before the oceanically-appropriate parameter range. It is demonstrated here that the circulation is controlled if a mechanism for ultimate removal of vorticity exists, even if it is active only in a narrow region near the boundary. Vorticity removal is modeled here as a viscosity enhanced very near the solid boundaries to roughly parameterize missing boundary physics like topographic interaction and three dimensional turbulence over the shelf. This boundary-enhanced viscosity allows western-intensified mean flows even when the inertial boundary width, is much wider than the frictional region because eddies flux vorticity from within the interior streamlines to the frictional region for removal. Using boundary-enhanced viscosity, western-intensified calculations are possible with lower interior viscosity than in previous studies. Interesting behaviors result: a boundary-layer balance novel to the model, calculations with promise for eddy parameterization, eddy-driven gyres rotating opposite the wind, and temporal complexity including basin resonances. I also demonstrate that multiple-gyre calculations have weaker mean circulation than single-gyres with the same viscosity and subtropical forcing. Despite traditional understanding, almost no inter-gyre flux occurs if no-slip boundary conditions are used. The inter-gyre eddy flux is in control only with exactly symmetric gyres and free slip boundaries. Even without the inter-gyre flux, the multiple-gyre circulation is weak because of sinuous instabilities on the jet which are not present in the single-gyre model. These modes efficiently flux vorticity to the boundary and reduce the circulation without an inter-gyre flux, postponing inertial domination to much smaller viscosities. Then sinuous modes in combination with boundary-enhanced viscosity can control the circulation.

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

Description: A comparison of monthly biogeochemical measurements made from 1993 to 1995, combined with hydrography and satellite altimetry, was used to observe the impacts of nine eddy events on primary productivity and particle flux in the Sargasso Sea. Measurements of primary production, thorium-234 flux, nitrate+nitrite, and photosynthetic pigments made at the US JGOFS Bermuda Atlantic Time-series Study (BATS) site were used. During the three years of this study, four out of six high thorium- 234 flux events over 1000 dpm/m2/d occurred during the passage of an eddy. Primary production nearly as high as the spring bloom maximum was observed in two modewater eddies (May 1993 and July 1995). The 1994 spring bloom at BATS was suppressed by the passage of an anticyclone. Distinct phytoplankton community shifts were observed in mode-water eddies, which had an increased percentage diatoms and dinoflagelletes, and in cyclones, which had an increased percentage cyanobacteria (excluding Prochlorococcus). The difference in the observations of mode-water eddies and cyclones may result from the age of the eddy, which was very important to the biological response. In general, eddies that were one to two months old elicited a large biological response; eddies that were three months old may show a biological response and were accompanied by high thorium flux measurements; eddies that were four months old or older did not show a biological response or high thorium flux. Our conceptual model depicting the importance of temporal changes during eddy upwelling and decay fit the observations well in all 7 upwelling eddies. Additional information is needed to determine the importance of deeper mixed layers and winter mixing to the magnitude of the eddy impacts. Also, sampling generally captured only the beginning, end, and lor edge of an eddy due to the monthly to semi-monthly frequency of the measurements made at BATS. Lagrangian studies, higher resolution time-series, and/or more spatial coverage is needed to provide additional information for improved C and N budgets in the Sargasso Sea and to complete our understanding of the temporal changes that occur in an eddy.

Description: Funding for this work was provided by NASA and NSF through the JGOFS Synthesis and Modeling Program.

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

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

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

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of 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: 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: Submitted in partial fulfillment of the requirements for the degree of Master of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution June 2006

Description: The continuous supply of heat and fresh water from the boundaries to the interior of the Labrador Sea plays an important role for the dynamics of the region and in particular, for the Labrador Sea Water formation. Thus, it is necessary to understand the factors governing the exchange of properties between the boundary and interior. A significant fraction of heat and fresh water, needed to balance the annual heat loss and to contribute to the seasonal freshening of the Labrador Sea, is thought to be provided by coherent long-lived anticyclonic eddies shed by the Irminger Current. The population, some properties, rates and direction of propagation of these anomalies are known but the evolution and the mechanism of their decay are still far from obvious. In this work I investigated their water mass properties and evolution under the strong wintertime forcing using hydrographic data from 1990-2004 and a 1-dimensional mixed layer model. There were 50 eddies found in the hydrographic data record, 48 of which were identified as anticyclones. Vertical structure of the eddies was investigated, leading to the categorization of all the anticyclones into three classes: 12 - with a fresh surface layer and no mixed layer, 18 - without a fresh layer and at least one mixed layer, and 18 with ambiguous vertical structure. Four eddies of the second group appeared to have cores extending to as deep as 1500 m vertically and an isopycnal displacement of 400-600 m. A number of eddies without a fresh water cap contained Labrador Sea Water from the previous year at mid-depths.

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

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

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

Description: Submitted in partial fulfillment of the requirements for the degree of 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: 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.