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

Description: This thesis explores the buoyancy-driven circulation in the Red Sea, using a combination of observations, as well as numerical modeling and analytical method. The first part of the thesis investigates the formation mechanism and spreading of Red Sea Overflow Water (RSOW) in the Red Sea. The preconditions required for open-ocean convection, which is suggested to be the formation mechanism of RSOW, are examined. The RSOW is identified and tracked as a layer with minimum potential vorticity and maximum chlorofluorocarbon-12. The pathway of the RSOW is also explored using numerical simulation. If diffusivity is not considered, the production rate of the RSOW is estimated to be 0.63 Sv using Walin’s method. By comparing this 0.63 Sv to the actual RSOW transport at the Strait of Bab el Mandeb, it is implied that the vertical diffusivity is about 3.4 x 10-5m2 s-1 . The second part of the thesis studies buoyancy-forced circulation in an idealized Red Sea. Buoyancy-loss driven circulation in marginal seas is usually dominated by cyclonic boundary currents on f-plane, as suggested by previous observations and numerical modeling. This thesis suggests that by including β-effect and buoyancy loss that increases linearly with latitude, the resultant mean Red Sea circulation consists of an anticyclonic gyre in the south and a cyclonic gyre in the north. In mid-basin, the northward surface flow crosses from the western boundary to the eastern boundary. The observational support is also reviewed. The mechanism that controls the crossover of boundary currents is further explored using an ad hoc analytical model based on PV dynamics. This ad hoc analytical model successfully predicts the crossover latitude of boundary currents. It suggests that the competition between advection of planetary vorticity and buoyancy-loss related term determines the crossover latitude. The third part of the thesis investigates three mechanisms that might account for eddy generation in the Red Sea, by conducting a series of numerical experiments. The three mechanisms are: i) baroclinic instability; ii) meridional structure of surface buoyancy losses; iii) cross-basin wind fields.

Description: This work is supported by Award Nos. USA 00002, KSA 00011 and KSA 00011/02 made by King Abdullah University of Science and Technology (KAUST) , National Science Foundation OCE0927017, and WHOI Academic Program Office.

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

Description: The nitrogen fixation and abundance of Trichodesmium colonies and their connections with physical processes were investigated through Video Plankton Recorder (VPR) and other observations collected in fall 2010 and spring 2011 in the western subtropical–tropical North Atlantic. A data processing procedure for estimating rare taxon abundance was devised to leverage the accuracy of manual classification and the effort savings of automatic classification. In fall 2010, local maxima in colony abundance were observed in a series of cyclones. We hypothesized Ekman transport convergence/divergence in cyclones/anticyclones as a driving mechanism and investigated the process using idealized three-dimensional models. Elevated abundances in anticyclones in spring 2011 were correlated with anomalously fresh water connected to river outflow. A bio-optical model based on carbon-normalized nitrogen fixation rates measured in fall 2010 and spring 2011 was used to estimate nitrogen fixation over the VPR transects. Mean VPR-based estimates of abundance and volume-specific nitrogen fixation rates at depth in the tropical North Atlantic were not inconsistent with estimates derived from conventional sampling methods compiled in a database by Luo et al. (2012). These findings did not reveal the systematic underestimation of deep colony populations and nitrogen fixation hypothesized by Davis and McGillicuddy (2006).

Description: This work was supported through a NASA Earth and Space Science Fellowship (NASA NNX11AL59H Understanding the role of the nitrogen- xing cyanobacterium Trichodesmium in the oceanic nitrogen and carbon cycles: in situ measurement, satellite observation, and biogeochemical modeling) as well as fellowship support from the Ocean Life Institute and Academic Programs O ce at WHOI. Additional grant support was provided by NSF OCE-0925284 Quanti cation of Trichodesmium spp. Vertical and Horizontal Abundance Patterns and Nitrogen Fixation in theWestern North Atlantic, NSF OCE-1048897 MOBY: Modeling Ocean Variability and Biogeochemical Cycles, NASA NNX13AE47G Physical and Biological Dynamics of Nonlinear Mesoscale Eddies: Satellite Observations, in situ Measurements, and Numerical Simulations on a Global Scale, and NASA NNX08AL71G Carbon cycling in the North Atlantic from regional to basin scales: satellite data, in situ observations, and numerical models.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 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 December 1996

Description: Two-thirds of the surface of the Earth is created at mid-ocean ridges where magmas rise from the mantle and cool to form the oceanic crust. The objective of this Thesis is to examine the influence of magma supply and eruptive processes on axial morphology, crustal construction, and the properties of crustal magma chambers at intermediate and fast spreading ridges. Variations in magma supply on time scales of ~100 Kyr generate along-axis changes in crustal thickness and temperature. Magma sill properties and hydrothermal activity are closely linked to spreading events which occur on much shorter time scales (ca. 10-100 yr) than the longer-term variations in magma supply reflected in along-axis changes in ridge morphology. The seismically constrained depths of ridge crest magma sills (〉1-2 km) are considerably deeper than the level of neutral buoyancy (100-400 m). The apparent inverse relationship between magma sill depth and spreading rate suggests that a thermally controlled permeability boundary, such as the solidus horizon, controls the depth at which magma ponds beneath mid-ocean ridges. Recent thermo-mechanical models predict that, at intermediate spreading rates, rift valley and magma sill formation are sensitive to small changes in crustal thickness and mantle temperature. Analysis of gravity at an intermediate spreading ridge shows that small differences in crustal thickness (300-700 m) and mantle temperature (10-15°C) are indeed sufficient to produce major changes in lithospheric strength and axial morphology. A stochastic model for the emplacement of dikes and lava flows with a bimodal distribution of lava flows is required to satisfy geological and geophysical constraints on the construction of the extrusive section. Most dikes are intruded within a narrow zone at the ridge axis. Short flows build up approximately half the extrusive volume. Occasional flows that pond at a considerable distance off-axis build up the remainder of the extrusive section. This Thesis underlines the importance of eruption dynamics in the emplacement of the uppermost volcanic layer of the crust and of the crustal thermal structure in controlling local variations in magma sill depth and ridge morphology.

Description: Author Posting. © American Geophysical Union, 2012. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Global Biogeochemical Cycles 26 (2012): GB4018, doi:10.1029/2011GB004192.

Description: A series of seasonally distributed measurements from the six largest Arctic rivers (the Ob', Yenisey, Lena, Kolyma, Yukon and Mackenzie) was used to examine the magnitude and significance of Arctic riverine DIC flux to larger scale C dynamics within the Arctic system. DIC concentration showed considerable, and synchronous, seasonal variation across these six large Arctic rivers, which have an estimated combined annual DIC flux of 30 Tg C yr−1. By examining the relationship between DIC flux and landscape variables known to regulate riverine DIC, we extrapolate to a DIC flux of 57 ± 9.9 Tg C yr−1for the full pan-arctic basin, and show that DIC export increases with runoff, the extent of carbonate rocks and glacial coverage, but decreases with permafrost extent. This pan-arctic riverine DIC estimate represents 13–15% of the total global DIC flux. The annual flux of selected ions (HCO3−, Na+, Ca2+, Mg2+, Sr2+, and Cl−) from the six largest Arctic rivers confirms that chemical weathering is dominated by inputs from carbonate rocks in the North American watersheds, but points to a more important role for silicate rocks in Siberian watersheds. In the coastal ocean, river water-induced decreases in aragonite saturation (i.e., an ocean acidification effect) appears to be much more pronounced in Siberia than in the North American Arctic, and stronger in the winter and spring than in the late summer. Accounting for seasonal variation in the flux of DIC and other major ions gives a much clearer understanding of the importance of riverine DIC within the broader pan-arctic C cycle.

Description: Funding for this work was provided through NSF-OPP-0229302 and NSF-OPP-0732985. Additional support to SET was provided by an NSERC Postdoctoral Fellowship.

Description: 2013-06-14

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: Two high-resolution mooring arrays extending from the outer shelf to the mid continental slope are used to elucidate shelf-basin exchange at the inflow to and the outflow from the Arctic Ocean. Pacific Water entering the Arctic Ocean forms the Western Arctic shelfbreak current along the Beaufort Sea slope. Data from the mooring array at 152°W—approximately 150 km east of Pt. Barrow, AK—reveals that this current has two distinct states in summer depending on the water mass it transports. When advecting Alaskan Coastal Water it is surface-intensified and both baroclinically and barotropically unstable. This configuration lasts about a month with an average transport of 0.7 Sv. When advecting Chukchi Summer Water the current is bottom-intensified and is only baroclinically unstable. This state also exists for approximately a month with an average transport of 0.6 Sv. The strong mean-to-eddy energy conversion causes both configurations of the current to spin down over a distance of a few hundred kilometers, suggesting that warm Pacific Water does not enter the Canadian Arctic Archipelago via this route. Dense water formed in the Nordic Seas overflows the Denmark Strait and undergoes vortex stretching, forming intense cyclones that propagate along the East Greenland slope. Data from the mooring array at 65°N—roughly 300 km downstream of Denmark Strait—was used to determine the full water column structure of the cyclones. On average a cyclone passes the array every other day in the vicinity of the 900 m isobath, although the depth range of individual cyclones ranges between the 500 m and 1600 m isobaths. The cyclones self-propagate at 0.45 m/s and are also advected by the mean flow of 0.27 m/s, resulting in a total propagation speed of 0.72 m/s. They have a peak azimuthal speed of 0.22 m/s at a radius of 7.8 km and contain overflow water in their core. In the absence of the cyclones, the background flow is dominated by the East Greenland Spill Jet. This is shown to be a year-round feature transporting 2–4 Sv of dense water equatorward along the upper continental slope.

Description: Financial support for this work was provided by National Science Foundation grants OCE-0726640 and OCE-0612143, by the Arctic Research Initiative at WHOI, by the Y-S Anonymous Fellowship from the Office of the Dean of Graduate Education at MIT, and by WHOI Academic Programs Office funds.

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

Description: The transformation of potential vorticity within and stability of nonlinear deep western boundary currents in an idealized tropical ocean are studied using a shallowwater model. Observational evidence indicates that the potential vorticity of fluid parcels in deep western boundary currents must change sign as they cross the equator, but this evidence is otherwise unable to clarify the process. A series of numerical experiments investigate this transformation in a rectangular basin straddling the equator. A mass source located in the northwestern corner feeds fluid into the domain where it is constrained to cross the equator to reach a distributed mass sink. Dissipation is included as momentum diffusion. The Reynolds number, defined as the ratio of the mass source per unit depth to the viscosity, determines the nature of the flow, and a critical value, Rec, divides its possible behavior into two regimes. For Re 〈 Rec, the flow is laminar and well described by linear theory. For Re just above the critical value, the flow is time-dependent, with cyclonic eddies forming in the western boundary current near the equator. For still larger Reynolds number, eddies of both signs emerge and form a complicated, interacting network that extends into the basin several deformation radii from the western boundary, as well as north and south of the equator. The eddy field is established as the mechanism for potential vorticity transformation in nonlinear cross-equatorial flow. The analysis of vorticity fluxes follows from the flux-conservative form of the absolute vorticity equation. It is shown that the zonally integrated meridional flux of vorticity across the equator using no slip boundary conditions is virtually zero even in the strongly nonlinear limit suggesting that the eddies are extremely efficient vorticity transfer agents. A decomposition of the vorticity fluxes into components due to mean advection, eddy transport, and friction, reveals the growth with Reynolds number of a turbulent boundary layer that exchanges vorticity between the inertial portion of the boundary current and a frictional sub-layer where modification is straightforward. A linear stability analysis of the shallow-water system in the tropical ocean examines the initial formation of the eddy field. The formulation assumes that the basic state is purely meridional and on a local f-plane. Realistic western boundary current profiles undergo a horizontal shear instability that is partially stabilized by viscosity. Calculations at several latitudes indicate that the instability is enhanced in the tropics where the internal deformation radius is a maximum. The linear stability analysis predicts a length scale of the disturbance, a location for its origin, and a critical Reynolds number that agree well with numerical results.

Description: Financial support for this research was provided by NSF grant number OCE- 9115915 and ONR ASSERT grant number N00014-94-1-0844.

Description: 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: We analyze bathymetric and gravity anomalies at five plume-ridge systems to constrain crustal and mantle density structure at these prominent oceanic features. Numerical models are then used to explore the physical mechanisms controlling plume-ridge interaction and to place theoretical constraints on the temperature anomalies, dimensions, and fluxes of the Icelandic and Galapagos plumes. In Chapter 1 we analyze bathymetric and gravity anomalies along the hotspot-influenced Galapagos Spreading Center. We find that the Galapagos plume generates along-axis bathymetric and mantle-Bouguer gravity anomalies (MBA) that extend 〉500 km east and west of the Galapagos Islands. The along-axis MBA becomes increasingly negative towards the plume center, reaching a minimum of ~-90 mGal near 91°W, and axial topography shallows by ~1.1 km toward the plume. These variations in MBA and bathymetry are attributed to the combined effects of crustal thickening and anomalously low mantle densities, both of which are due to a mantle temperature anomaly imposed beneath the ridge by the Galapagos plume. Passive mantle flow models predict a temperature anomaly of 50±25°C is sufficient to produce the 2-4 km excess crust required to explain the along-axis anomalies. 70-75% of the along-axis bathymetric and MBA variations are estimated to arise from the crust with the remaining 25-30% generated by the anomalously hot, thus low-density mantle. Along Cocos-plate isochrons, bathymetric and MBA variations increase with increasing isochron age, suggesting the subaxial mantle temperature anomaly was greater in the past when the plume was closer, to the ridge axis. In addition to the Galapagos plume-ridge system, in Chapter 2 we examine alongisochron bathymetric and MBA variations at four other plume-ridge systems associated with the Iceland, Azores, Easter and Tristan hotspots. We show that residual bathymetry (up to 4.7 km) and mantle-Bouguer gravity anomalies (up to -340 mGal) are greatest at on-axis plumes and decreases with increasing ridge-hotspot separation distance, until becoming insignificant at a plume-ridge separation of ~500 km. Along-isochron widths of bathymetric anomalies (up to 2700 km) decrease with increasing paleo-spreading rate, reflecting the extent to which plume material flows along-axis before being swept away by the spreading lithosphere. Scaling arguments suggest an average ridgeward plume flux of -2.2x106 km/my. Assuming that the amplitudes of the MBA and bathymetric anomalies reflect crustal thickness and mantle density variations, passive mantle flow models predict maximum subaxial mantle temperature anomalies to be 150-225°C for ridge-center plumes, which decrease as the ridges migrate away from the plumes. The dynamics of mantle flow and melting at ridge-centered plumes are investigated in Chapters 3 using three-dimensional, variable-viscosity, numerical models. Three buoyancy sources are examined: temperature, melt depletion, and melt retention. The width W to which a plume spreads along a ridge axis depends on plume volume flux Q, full spreading rate U, buoyancy number B = (QΔρg)/(48η0U2), and ambient/plume viscosity contrast ϒ according to W=2.37(Q/U)l/2(Bϒ)0.04. Thermal buoyancy is first order in controlling along-axis plume spreading while latent heat loss due to melting, and depletion and retention buoyancy forces contribute second order effects. Two end-member models of the Iceland-Mid-Atlantic Ridge (MAR) system are examined. The first endmember model has a broad plume source of radius 300 km, temperature anomaly of 75°C, and volume flux of 1.2xl07 km3/my. The second model has a narrower plume source of radius 60 km, temperature anomaly of l70°C, and flux of 2.1 x106 km3/my. The first model predicts successfully the observed crustal thickness, topographic, and MBA variations along the MAR, but the second model requires substantial along-axis melt transport in order to explain the observed along-axis variations in crustal thickness, bathymetry, and gravity. We favor this second model because it predicts a mantle P-wave velocity reduction in the plume of ~2% as consistent with recent seismic observations beneath Iceland. Finally in Chapter 4 we use three-dimensional numerical models to investigate the interaction of plumes and migrating midocean ridges. Scaling laws of axial plume spreading width Ware derived first for stationary ridges and off-axis plumes, which yield results consistent with those obtained from independent studies of Ribe [1996]. Wand the maximum plume-ridge interaction distance Xmax again scale with (Q/U)l/2 as in the case of ridge-centered plumes and increase with ϒ and buoyancy number. In the case of a migrating ridge, Xmax is reduced when a ridge migrates toward the plume due to excess drag of the faster-moving leading plate, and enhanced when a ridge migrates away from the plume due to reduced drag of the slower-moving trailing plate. Thermal erosion of the lithospheric boundary layer by the previously ridge-centered plume further enhances Wand Xmax but to a degree that is secondary to the differential migration rates of the two plates. Model predictions are compared with observed along-isochron bathymetric and MBA variations at the Galapagos plume-ridge system. The anomaly amplitudes and widths, as well as the increase in anomaly amplitude with age are predicted with a plume source temperature anomaly of 80-120°C, radius of 80-100 km, and volume flux of 4.5x106 km3/m.y. Our numerical models also predict crustal production rates of the Galapagos Islands consistent with those estimated independently using the observed island topography. Predictions of the geochemical signature of the plume along the present-day ridge suggest that mixing between the plume and ambient mantle sources is unlikely to occur in the asthenosphere or shallow crust, but most likely deeper in the mantle possibly by entrainment of ambient mantle as the plume ascends through the depleted portion of the mantle from its deep source reservoir.

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 1997

Description: The formation of new oceanic crust is the result of a complex geodynamic system in which mantle rises beneath spreading centers and undergoes decompression melting. The melt segregates from the matrix and is focused to the rise axis, where it is eventually intruded and/or erupted to form the oceanic crust. This thesis combines surface observations with laboratory studies and geodynamic modeling to study this crustal-production system. Quantitative modeling of the crustal and mantle contributions to the axial gravity and topography observed at the East Pacific Rise shows that the retained melt fraction in the mantle is small (〈3%) and is focused into a narrow column extending up to 70 km beneath the ridge axis. Consistent with geochemical constraints, the extraction of melt from the mantle therefore appears to be efficiently focus melt toward the ridge axis. A combination of laboratory and numerical studies are used to constrain the pattern of mantle flow beneath highly-segmented ridges. Even when the buoyant component of mantle flow is constrained to be two-dimensional, laboratory studies show that a segmented ridge will drive three-dimensional mantle upwelling. However, using reasonable mantle parameters in numerical models, it is difficult to induce large-amplitude three-dimensional mantle upwelling at the relatively short wavelengths of individual segments (~50 km). Instead, a simple model of three-dimensional melt migration shows that the observed segment-scale variations in crustal thickness can be explained by focusing of melt as it upwells through a more two-dimensional mantle flow field. At the Reykjanes Ridge, the melt appears to accumulate in small crustal magma chambers, before erupting in small batches to form numerous overlapping hummocky lava flows and small volcanoes. This suggests that crustal accretion, particularly at slow-spreading centers, may be a highly discontinuous process. Long-wavelength variations in crustal accretion may be dominated by variations in mantle upwelling while short-wavelength, segment-scale variations are more likely controlled by a complex three-dimensional processes of melt extraction and magma eruption.

Description: During my first three years in the Joint Program, I was supported by an National Science Foundation Graduate Student Fellowship. Other support has been derived from National Science Foundation grants OCE-9296017, OCE-9224738, OCE-9215544, and EAR grant 93-07400.

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

Description: This thesis addresses the issue, "Which approach to instabilities-temporal, spatial or pulse theory- is the most appropriate model for the Gulf Stream?" I also address the question of how the observations might be compared to theory. This thesis consists of two closely related parts: analytical studies that compare the three types of instability using the same realistic velocity and topography profiles; and numerical modeling that uses a continuous forcing function to examine the three types of theory in the direct context of the Gulf Stream.

Description: My first three years in the Joint Program were supported by the National Science Foundation under grant OCE-9011066 and last two and half years under NSF grant OCE-9314140.

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

Description: When ocean waves in deep water interact with a current, the direction of propagation and characteristics of the waves such as height and length are affected. Swell in the open ocean can undergo significant refraction as it passes through major current systems like the Gulf Stream or Antarctic Circumpolar Current. Remote sensing techniques such as synthetic aperture radars (SAR) have the potential to detect wave systems over a wide geographical area. Combining a model for wave refraction in the presence of currents with SAR measurements, the inverse problem of using the measured wave data can be solved to determine the direction and magnitude of the intervening currents. In this study the behavior of swell measured by SAR on a satellite pass over the Gulf Stream is examined. The refraction predicted by a numerical model under conditions of varying current profiles and velocities is compared to SAR generated wave spectra. By matching the current profile which results in the best correlation of wave refraction to the SAR data, the tomographic problem of measuring the Gulf Stream current is solved. The best correlation between the model and SAR data is obtained when a current is modeled by a top hat velocity profile with a direction of 75° and a current speed of 2 m/s. The direction agrees with that visually observed from the SAR images, and the direction and speeds are close to the Coast Guard estimates for the Gulf Stream at the time of the SEASAT,pass. The current profiles used did not take into account a possible widening of the Gulf Stream at the position of the satellite overpass. There is a great deal of scatter in the SAR data, both before and in the Gulf Stream, so it is difficult to correlate every point with specific current behavior, but the increase in wave length and change in wave angle in the center of the Gulf Stream seem to indicate that there may be a non-uniform feature such as the formation of an eddy or other lateral variability near the current's edge.

Description: I was supported by the U. S. Navy.

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

Description: Interactions between the ocean circulation in sub-ice shelf cavities and the overlying ice shelf have received considerable attention in the context of observed changes in flow speeds of marine ice sheets around Antarctica. Modeling these interactions requires parameterizing the turbulent boundary layer processes to infer melt rates from the oceanic state at the ice-ocean interface. Here we explore two such parameterizations in the context of the MIT ocean general circulation model coupled to the z-coordinates ice shelf cavity model of Losch (2008). We investigate both idealized ice shelf cavity geometries as well as a realistic cavity under Pine Island Ice Shelf (PIIS), West Antarctica. Our starting point is a three-equation melt rate parameterization implemented by Losch (2008), which is based on the work of Hellmer and Olbers (1989). In this form, the transfer coefficients for calculating heat and freshwater fluxes are independent of frictional turbulence induced by the proximity of the moving ocean to the fixed ice interface. More recently, Holland and Jenkins (1999) have proposed a parameterization in which the transfer coefficients do depend on the ocean-induced turbulence and are directly coupled to the speed of currents in the ocean mixed layer underneath the ice shelf through a quadratic drag formulation and a bulk drag coefficient. The melt rate parameterization in the MITgcm is augmented to account for this velocity dependence. First, the effect of the augmented formulation is investigated in terms of its impact on melt rates as well as on its feedback on the wider sub-ice shelf circulation. We find that, over a wide range of drag coefficients, velocity-dependent melt rates are more strongly constrained by the distribution of mixed layer currents than by the temperature gradient between the shelf base and underlying ocean, as opposed to velocity-independent melt rates. This leads to large differences in melt rate patterns under PIIS when including versus not including the velocity dependence. In a second time, the modulating effects of tidal currents on melting at the base of PIIS are examined. We find that the temporal variability of velocity-dependent melt rates under tidal forcing is greater than that of velocity-independent melt rates. Our experiments suggest that because tidal currents under PIIS are weak and buoyancy fluxes are strong, tidal mixing is negligible and tidal rectification is restricted to very steep bathymetric features, such as the ice shelf front. Nonetheless, strong tidally-rectified currents at the ice shelf front significantly increase ablation rates there when the formulation of the transfer coefficients includes the velocity dependence. The enhanced melting then feedbacks positively on the rectified currents, which are susceptible to insulate the cavity interior from changes in open ocean conditions.

Description: National Science and Engineering Research Council of Canada

Description: Author Posting. © American Geophysical Union, 2012. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 13 (2012): Q08014, doi:10.1029/2012GC004163.

Description: Mapping and sampling of 18 eruptive units in two study areas along the Galápagos Spreading Center (GSC) provide insight into how magma supply affects mid-ocean ridge (MOR) volcanic eruptions. The two study areas have similar spreading rates (53 versus 55 mm/yr), but differ by 30% in the time-averaged rate of magma supply (0.3 × 106 versus 0.4 × 106 m3/yr/km). Detailed geologic maps of each study area incorporate observations of flow contacts and sediment thickness, in addition to sample petrology, geomagnetic paleointensity, and inferences from high-resolution bathymetry data. At the lower-magma-supply study area, eruptions typically produce irregularly shaped clusters of pillow mounds with total eruptive volumes ranging from 0.09 to 1.3 km3. At the higher-magma-supply study area, lava morphologies characteristic of higher effusion rates are more common, eruptions typically occur along elongated fissures, and eruptive volumes are an order of magnitude smaller (0.002–0.13 km3). At this site, glass MgO contents (2.7–8.4 wt. %) and corresponding liquidus temperatures are lower on average, and more variable, than those at the lower-magma-supply study area (6.2–9.1 wt. % MgO). The differences in eruptive volume, lava temperature, morphology, and inferred eruption rates observed between the two areas along the GSC are similar to those that have previously been related to variable spreading rates on the global MOR system. Importantly, the documentation of multiple sequences of eruptions at each study area, representing hundreds to thousands of years, provides constraints on the variability in eruptive style at a given magma supply and spreading rate.

Description: This work was supported by the National Science Foundation grants OCE08–49813, OCE08–50052, and OCE08– 49711.

Description: 2013-02-25

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 2012

Description: The structure of the oceanic lithosphere results from magmatic and extensional processes taking place at mid-ocean ridges (MORs). The temporal and spatial scales of the variability of these two processes control the degree of heterogeneity of the oceanic lithosphere, represented by two end-member models: the classical Penrose Model exemplified by layered magmatic crust formed along fast-spreading MORs, e.g., East Pacific Rise (EPR); and the recently defined Chapman Model describing heterogeneous mafic and ultramafic lithosphere formed in settings of oceanic detachment faulting common along slow-spreading MORs, e.g., Mid-Atlantic Ridge (MAR). This thesis is using advanced marine geophysical methods (including finite-difference wave propagation modeling, 3D multi-channel seismic reflection imaging, waveform inversion, streamer tomography, and near-bottom magnetics) to study lithospheric accretion processes in MORs characterized by contrasting tectono-magmatic settings: the magmatically dominated EPR axis between 9°30'-10°00’N, and the Kane Oceanic Core Complex (KOCC), a section of MAR lithosphere (23°20’-23°38’N) formed by detachment faulting. At the EPR study area, I found that the axial magma chamber (AMC) melt sill is segmented into four prominent 2-4-km-long sections spaced every ~5- 10 km along the ridge axis characterized by high melt content (〉95%). In contrast, within the intervening sections, the AMC sill has a lower melt content (41-46%). The total magma volume extracted from the AMC sill was estimated of ~46 × 106 m3, with ~24 × 106 m3 left unerupted in the upper crust as dikes after 2005-06 eruption. At the KOCC, I used streamer tomography to constrain the shallow seismic velocity structure. Lithological interpretation of the seismic tomographic models provides insights into the temporal and spatial evolution of the melt supply at the spreading axis as the KOCC formed and evolved. Investigation of a magnetic polarity reversal boundary in crosssection at the northern boundary of KOCC suggests that the boundary (representing both a frozen isotherm and an isochron) dips away from the ridge axis along the Kane transform fault scarp, with a west-dipping angle of ~45° in the shallow (〈1 km) crust and 〈20° in the deeper crust.

Description: This thesis was funded by National Science Foundation grants OCE-9987004, OCE- 0621660 and OCE-0327885, WHOI Academic Program Office and WHOI Deep Ocean Exploration Institute.

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

Description: A triangular CTD/ADCP survey was made across the Kuroshio west of Kyushu aboard the R/V Thompson during January, 1986 in order to investigate the water properties and flow field in the Kuroshio. A similar CTD survey was made in July, 1986 aboard the R/V Washington to study the seasonal variability in the Kuroshio. The Kuroshio in this region exhibited a marked seasonal change in its near-surface stratification and water properties. In January, the Kuroshio water was separated from the vertically well-mixed coastal water over the shelf by a strong front located near the shelf break. Horizontal mixing between the Kuroshio and coastal water was observed but was limited near the shelf break. In July, surface coastal water extended far past the shelf break over the Kuroshio region near the surface, and in turn, Kuroshio water intruded onto the shelf near the bottom. Mixing between the Kuroshio and coastal water was found over much of the mid and outer shelf and upper slope, spanning a cross-stream distance of 75 km. In addition, evidence of deep vertical mixing within the Kuroshio itself was found near 32.0°N and 128.2°E, most likely due to internal tidal mixing over the slope. Since Loran C navigation coverage in the study region was poor during the R/V Thompson cruise, a simple averaging technique has been used to convert the ADCP data into an absolute velocity. An error analysis shows that the total error in the absolute ADCP velocity was about ±5 cm/s. The absolute geostrophic velocity using the absolute Doppler velocity at 60 m as the reference velocity was then calculated for the sides of the triangle. The results show that the ADCP velocity shear was in good agreement with the geostrophic shear in the Kuroshio. The Kuroshio flowed through the western section as a coherent current, then split into two streams around a tall seamount as it left through the eastern section. Some recirculation also occurred between the core of the Kuroshio and the slope as well as near the seamount. The geostrophic velocity field calculated relative to the bottom missed some of the important features of the true flow field such as splitting of the Kuroshio and the recirculation in the slope region. The volume, salt and heat transports of the Kuroshio during the January 1986 survey have been cakulated using the absolute geostrophic velocity and CTD data. The volume transport of the Kuroshio west of Kyushu in January 1986 was 31.7± 2.0 Sv, which is comparable to that of the Gulf Stream in the Florida Strait. The volume transport through the triangle was conserved within measurement uncertainty, so that a streamfuction field can be defined by the transport. The resulting streamlines clearly show the structure of the flow field in the Kuroshio and its adjacent currents during the survey. The advective heat transport of the Kuroshio west of Kyushu in January 1986 was 28.2 ± 1.8 x 1014 W. The salt transport in January 1986 was about 108.0 ± 7.3 x 1010 kg/s, and the net salt flux was zero within measurement error. Analysis of the potential vorticity based on the January 1986 absolute geostrophic velocity field shows that the total potential vorticity in the Kuroshio may be approximately given by the product of the vertical gradient of the potential density and the sum of the planetary and relative vorticities. The distribution of relative vorticity plays a significant role in determining the structure of the potential vorticity in the Kuroshio. The path of the Kuroshio can be traced in the field of potential vorticity. Facing in the direction of the current, the axis of the maximum velocity is located to the right of the core of maximum potential vorticity. Finally, the Kuroshio was potentially unstable since the gradient of potential vorticity changed its sign on potential density surfaces across the Kuroshio.

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 1992

Description: This thesis consists of two parts: (I) variability of currents and water properties in late spring in the northern Great South Channel and (II) numerical study of stratified tidal rectification over Georges Bank. In part I, the combined analysis of CTD, ADCP, satellite-track drifters, and zooplankton distributions in the northern Great South Channel clearly shows (1) the seasonal evolution of the surface low salinity plume, (2) the threedimensional structure of residual flow, and (3) a coherent relationship between the surface low salinity plume and high concentration of zooplankton. In part II, the numerical model of Georges Bank shows that as the fluid becomes stratified , tidal mixing and rectification intensify both along- and cross-bank residual currents and modify the vertical structure of the flow. Along- and cross-bank residual currents increase as either stratification increases or the depth of the bank decreases. Model results over Georges Bank are in good agreement with observation, particularly in the position of the tidal mixing front and residual currents on the northern flank of the Bank.

Description: This research was supported by the National Science Foundation under grants OCE 87-13988 and OCE 91-01034 and by the National Center for Atmospheric Research (NCAR) under computer time grants Nos. 35781029 and 35781035.

Description: Author Posting. © American Geophysical Union, 2012. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 117 (2012): C06010, doi:10.1029/2011JC007652.

Description: We propose a conceptual model for an Arctic sea that is driven by river runoff, atmospheric fluxes, sea ice melt/growth, and winds. The model domain is divided into two areas, the interior and boundary regions, that are coupled through Ekman and eddy fluxes of buoyancy. The model is applied to Hudson and James Bays (HJB, a large inland basin in northeastern Canada) for the period 1979–2007. Several yearlong records from instruments moored within HJB show that the model results are consistent with the real system. The model notably reproduces the seasonal migration of the halocline, the baroclinic boundary current, spatial variability of freshwater content, and the fall maximum in freshwater export. The simulations clarify the important differences in the freshwater balance of the western and eastern sides of HJB. The significant role played by the boundary current in the freshwater budget of the system, and its sensitivity to the wind-forcing, are also highlighted by the simulations and new data analyses. We conclude that the model proposed is useful for the interpretation of observed data from Arctic seas and model outputs from more complex coupled/climate models.

Description: We thank NSERC and the Canada Research Chairs program for funding. FS acknowledges support from NSF OCE–0927797 and ONR N00014-08-10490.

Description: 2012-12-20

Description: Author Posting. © American Geophysical Union, 2012. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Global Biogeochemical Cycles 26 (2012): GB0E02, doi:10.1029/2012GB004299.

Description: While much of the dissolved organic carbon (DOC) within rivers is destined for mineralization to CO2, a substantial fraction of riverine bicarbonate (HCO3−) flux represents a CO2 sink, as a result of weathering processes that sequester CO2 as HCO3−. We explored landscape-level controls on DOC and HCO3− flux in subcatchments of the boreal, with a specific focus on the effect of permafrost on riverine dissolved C flux. To do this, we undertook a multivariate analysis that partitioned the variance attributable to known, key regulators of dissolved C flux (runoff, lithology, and vegetation) prior to examining the effect of permafrost, using riverine biogeochemistry data from a suite of subcatchments drawn from the Mackenzie, Yukon, East, and West Siberian regions of the circumboreal. Across the diverse catchments that we study, controls on HCO3− flux were near-universal: runoff and an increased carbonate rock contribution to weathering (assessed as riverwater Ca:Na) increased HCO3− yields, while increasing permafrost extent was associated with decreases in HCO3−. In contrast, permafrost had contrasting and region-specific effects on DOC yield, even after the variation caused by other key drivers of its flux had been accounted for. We used ionic ratios and SO4 yields to calculate the potential range of CO2 sequestered via weathering across these boreal subcatchments, and show that decreasing permafrost extent is associated with increases in weathering-mediated CO2 fixation across broad spatial scales, an effect that could counterbalance some of the organic C mineralization that is predicted with declining permafrost.

Description: Funding for this work was provided through NSF-OPP-0229302 and NSF-OPP-0732985. Additional support to S.E.T. was provided by an NSERC Postdoctoral Fellowship.

Description: 2013-02-21

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

Description: The dynamical aspects involved in the assimilation of altimeter data in a numerical ocean model have been investigated. The model used for this study is a quasi-geostrophic model of the Gulf Stream region. The data that have been assimilated are maps of sea surface height which have been obtained as the superposition of sea surface height variability deduced from the Geosat altimeter measurements and a mean field constructed from historical hydrographic data. The method used for assimilating the data is the nudging technique. Nudging has been implemented in such a way as to achieve a high degree of convergence of the surface model fields toward the observations. We have analyzed the mechanisms of the model adjustment, and the final statistical equilibrium characteristics of the model simulation when the surface data are assimilated. Since the surface data are the superposition of a mean component and an eddy component, in order to understand the relative role of these two components in determining the characteristics of the final st atistical steady state, we have considered two different experiments: in the first experiment only the climatological mean field is assimilated, while in the second experiment the total surface streamfunction field (mean + eddies) has been used. We have found that the mean component of the surface data determines, to a large extent, the structure of the flow field in the subsurface layers, while the eddy field, as well as the inflow/outflow conditions at the open boundaries, affect its intensity. In particular, if surface eddies are not assimilated only a weak flow develops in the two deeper model layers where no inflow/ outflow is prescribed at the boundaries. Comparisons of the assimilation results with available in situ observations show a considerable improvement in the degree of realism of the climatological model behavior, with respect to the model in which no data are assimilated. In particular, the possibility of building into the model more realistic eddy characteristics, through the assimilation of the surface eddy field, proves very successful in driving components of the mean model circulation that are in good agreement with the available observations.

Description: This research was carried out with the support of the National Aeronaut ics Space Administration, through a contract to MIT from the Jet Propulsion Laboratory, # 958208, as a part of the TOPEX-Poseidon investigation.

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

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

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

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

Description: Hydrographic and expendable current profiler (XCP) data taken during the Gulf of Cadiz Expedition in September 1988 are analyzed to diagnose the mixing and dynamics of the Mediterranean outflow. The overall structure of the outflow is consistent with that described in the historical literature (Heezen and Johnson, 1969). This data shows that the overflow transport doubles from .85 Sv to 1.9 Sv, and that the velocity weighted salinity decreases from 37.8 pss to 36.7 pss in the first 60 km of the path. The core salinity of the neutrally buoyant outflow near Cape St. Vincent is 36.6 pss, which indicates that most of the mixing has taken place close to the Strait in the initial descent of the outflow. Cross stream variations in the overflow T/S properties increase as the flow spreads from 10 km to 90 km wide. The outflow begins with less than a 0.5°C across-stream variation in temperature in the Strait with the saltiest, coldest water to the south and slightly fresher and warmer outflow to the north. As the outflow spreads, the northern near-shelf flow follows a path higher in the water column and mixes with warmer North Atlantic water than does the deeper offshore flow. Within the first 100 km, the cross stream variation in temperature on an isopycnal becomes more than a 2°C. The flow eventually settles along two preferred isopycnals: 27.5 and 27.8 (Zenk 1975b). The spreading of the flow contains both a barotropic and baroclinic character. The average change in angle above and below the maximum velocity of the outflow is 8°while at the edges of the flow the average direction of the outflow diverges by as much as 50°. Gradient Richardson numbers less than 1/4 are found in the interface (up to 50 m thick) between westward flowing Mediterranean water and eastward flowing North Atlantic water, even though there is a strong stabilizing stratification present. Bulk Froude numbers greater than 1 are found near the Strait coincident with the vigorous mixing noted above. Lower bulk Froude numbers were observed in regions where less entrainment was taking place. The momentum balances are diagnosed using hydrographic and XCP data. Evaluation of the cross stream momentum balance shows the importance of advection as the flow makes a 90 degree inertial turn upon entering the Gulf of Cadiz. A form of the Bernoulli function can be evaluated to infer the total stress (entrainment and bottom drag) acting on the outflow. This stress is as large as 5 Pa within 20 km of the Strait, while further downstream the stress decreases to about 1/2 Pa. The entrainment stress estimated from the property fluxes reaches a maximum of about 0.8 Pa near section C, indicating that bottom stress is dominant. Near the Strait, advection, bottom drag and the Coriolis force are all critical to the dynamics of the outflow. Further downstream, the outflow becomes a damped geostrophic current. A simple geostrophic adjustment model is used to show that in the absence of frictional stresses, the outflow would very quickly become geostrophically balanced and descend only about 10 m down the continental slope. Thus, friction is critical for the outflow to cross isobaths. A simple numerical model that uses a Froude number dependent entrainment and a quadratic bottom friction law is used to simulate the outflow (Price and Baringer, 1993). Some of the properties of the outflow including localized entrainment, large stresses and high Rossby number of the flow (initially as high as 0.6), are simulated rather well, though the model overestimates the magnitude of the outflow current. We suspect that this is a consequence of assuming a passive ocean. Two different methods for specifying the broadening of the flow are compared: one using the highly parameterized concept of Ekman spreading, the other using the conservation of potential vorticity. The potential vorticity broadening more accurately reproduces the observed width of the flow near Cape St. Vincent where the width varies inversely with the bottom slope. However, both methods produce essentially the same equilibrium temperature, salinity and transport of the outflow which is a testament to the robustness of the model solution. with the formation process of NADW.

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

Description: During January-March, Scotian Shelf water has been observed to flow episodically from the southwestern Scotian Shelf directly across the Northeast Channel to Georges Bank. The possible factors that allow Scotian Shelf water to break the topographic constraint presented by the Northeast Channel and flow directly to Georges Bank are considered. As a simple analog to the flow over the southwestern Scotian Shelf near the Northeast Channel, the adjustment of a barotropic current near a shelf-break to a sharp bend in the shelf topography is studied numerically. For parameters within the oceanographic range, the adjustment to the bend is smooth and steady with no eddies shed at the corner. The vorticity dynamics allow a balance between the vortex stretching in the flow and the curvature in the flow. This is possible since the bend is a right-hand one facing downstream, a similar balance not being possible for a left-hand bend, in which case eddy formation is likely. A simple model of this balance clarifies the vorticity dynamics and provides the scaling rc = √eL/0.765 for any streamline in the flow, where rc is the radius of curvature at the corner, E = u0/fL and L = h0/b, where uo is the initial speed, f the coriolis parameter, h0 the initial depth and b the bottom slope. These results show that other factors such as stratification, wind stress, and time-dependent inflow must play a role in any flow across the Northeast Channel.

Description: I am very grateful to the US-GLOBEC program for providing the funding for this study (N.S.F. grant OCE-9313671).

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 1995

Description: Geosat altimeter data and numerical model output are used to examine the circulation and dynamics of the Antarctic Circumpolar Current (ACC). The mean sea surface height across the ACC has been reconstructed from height variability measured by the Geosat altimeter, without assuming prior knowledge of the geoid. For this study, an automated technique has been developed to estimate mean sea surface height for each satellite ground track using a meandering Gaussian jet model, and errors have been estimated using Monte Carlo simulation. The results are objectively mapped to produce a picture of the mean Subantarctic and Polar Fronts, which together comprise the major components of the ACC. The locations of the fronts are consistent with in situ observations and indicate that the fronts are substantially steered by bathymetry. The jets have an average Gaussian width of about 44 km in the meridional direction and meander about 75 km to either side of their mean locations. The width of the fronts is proportional to 1/f, indicating that with constant stratification, the width is proportional to the baroclinic. Rossby radius. The average height difference across the Subantarctic Front (SAF) is 0.7 m and across the Polar Front (PF) 0.6 m. The mean widths of the fronts are correlated with the size of the baroclinic Rossby radius. The meandering jet model explains between 40% and 70% of the height variance along the jet axes. Bathymetric constrictions are associated with increased eddy variability, a smaller percentage of which may be explained by the meandering of the ACC fronts, indicating that propagating eddies and rings may be spawned at topographic features. Detailed examination of spatial and temporal variability in the altimeter data indicates a spatial decorrelation scale of 85 km and a temporal e-folding scale of 34 days. The sea surface height variability is objectively mapped using these scales to define autocovariance functions. The resulting maps indicate substantial evidence of mesoscale eddy activity. Over 17-day time intervals, meanders of the PF and SAF appear to elongate, break off as rings, and propagate. Statistical analysis of ACC variability from altimeter data is conducted using empirical orthogonal functions (EOFs ). The first mode EOF describes 16% of the variance in total sea surface height across the ACC; reducing the domain into basin scales does not significantly increase the variance represented by the first EOF, suggesting that the scales of motion are relatively short, and may be determined by local instability mechanisms rather than larger basin scale processes. Likewise, frequency domain EOFs indicate no statistically significant traveling wave modes. The momentum balance of the ACC has been investigated using both output from a high resolution primitive equation model and sea surface height measurements from the Geosat altimeter. In the Semtner-Chervin general circulation model, run with approximately quarter-degree resolution and time varying ECMWF winds, topographic form stress is the dominant process balancing the surface wind forcing. Detailed examination of form stress in the model indicates that it is due to three large topographic obstructions located at Kerguelen Island, Campbell Plateau, and Drake Passage. In order to reduce the effects of standing eddies, the model momentum balance is considered in stream coordinates; vertically integrated through the entire water column, topographic form drag is the dominant balance for wind stress. However, at mid-depth the cross-stream momentum transfer is dominated by horizontal biharmonic friction. In the upper ocean, horizontal friction, mean momentum flux divergence, transient momentum flux divergence, and mean vertical flux divergence all contribute significantly to the momentum balance. Although the relative importance of individual terms in the momentum balance does not vary substantially along streamlines, elevated levels of eddy kinetic energy are associated with the three major topographic features. In contrast, altimeter data show elevated energy levels at many more topographic features of intermediate scales, suggesting that smaller topographic effects are better able to communicate with the surface in the real ocean than in the model. Transient Reynolds stress terms play a small role in the the overall momentum balance; nonetheless, altimeter and model measurements closely agree, and suggest that transient eddies tend to accelerate the mean flow, except in the region between the major fronts which comprise the ACC. Potential vorticity is considered in the model output along Montgomery streamfunction. Even at about 1000 m depth, it varies in response to wind forcing, largely as a result of changes in vertical stratification, indicating that forcing and dissipation do not locally balance in the Southern Ocean. In order to compare model and altimeter potential vorticity estimates, two different proxies for potential vorticity on surface streamlines are considered. Both proxies show very similar results for model and altimeter, suggesting that differences in surface streamlines estimated by the altimeter and the model are not significant in explaining the Southern Ocean flow. The proxies are both roughly conserved along surface height contours but undergo substantial jumps near topographic features. However, they cannot capture stratification changes which may be critically important to the overall potential vorticity balance.

Description: Funding for this research was provided by an Office of Naval Research graduate student fellowship and National Aeronautics and Space Administration contract NAGW-1666.

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 1991

Description: A set of four hydrographic sections through the Brazil Current are analyzed to identify downstream changes in the Brazil Current. The data, from the Thomas Washington Marathon Cruise, Leg 9, are at 27, 31, 34 and 36°S. The region they span details the change of the current from a relatively small near surface feature to a large, deep current. While the Brazil Current does not appear to develop transports as large as those found in the Gulf Stream, the calculated transports greatly exceed previous estimates. At 27°S the current extends down to approximately 700 m and transports 12 Sv southward; this value is consistent with previous estimates farther north. Downstream, surface layer transport increases, the current deepens, and the transport reaches a maximum of 80 Sv at 36°S. Part of the growth comes from the tight recirculation found just offshore of the Brazil Current. The recirculation strengthens and deepens to the south, with a minimum transport of 4 Sv north at 27°S and a maximum of 33 Sv at 36°S. The change in the current is also reflected in its shear profiles. At 27°S Brazil Current shear is found only in the upper portion of the water column, over the continental slope. Downstream, the current moves off the slope into deeper water and develops top-to-bottom, monotonic shear. To obtain velocity from the shear profiles, zero velocity surfaces are chosen based on conservative use of tracer information. A simple basin-wide model is used at 31°S to tie limits on the size of the Brazil Current and recirculation to various limits on layer-to-layer exchanges south of the section. The multi-layer model - based on changes in depth of several isotherms is used to extend the interpretation of the current beyond that of an isolated ocean feature. The model is required to conserve mass in each layer, either by applying barotropic transports or by allowing layer-to-layer exchanges south of the section. Solutions are deemed acceptable if the sense, or direction, of the various layer-to-layer conversions are consistent with accepted ideas of water mass formation. Initially, a two layer model is employed. Governed by the conservation of mass in each layer, the two layer model has only one constraint on the resulting solutions: a conversion of cold-to-warm water in the south (or the surface layer flowing north and the deep layer flowing south). Such a meridional flow pattern is consistent with the equatorward heat flux in the South Atlantic. The single constraint, however, is not strong enough to limit the solution region in any significant way. The final model presented has four layers, and acceptable solutions have the net transports of the surface layer and the bottom water northward and form intermediate water from North Atlantic Deep Water in the south. The resulting solution set has a fairly small range of transports for the Brazil Current, with surface layer transports between 20 and 35 Sv; this range is consistent with the value calculated from hydrographic data at 31°S. Given the complex interleavings of the South Atlantic water masses, the four layer model performs remarkably well. Finally, total potential vorticity is calculated from the hydrographic sections. Contrary to what one might expect, the reference level choice is not a significant problem: where currents are large, most of the signal in relative potential vorticity comes from the measured shear, and where currents are small, the relative potential vorticity is not significant compared to the planetary vorticity. Unfortunately, the process of taking two horizontal derivatives of the density field results in a jittery relative potential vorticity signal. As a result, a clear potential vorticity profile could not be constructed for the current. This variablitiy may be real -the ocean is frequently much noisier than one imagines. It may also be possible, though, to smooth the data sufficiently so that a cleaner picture emerges. Despite the problems involved in obtaining a quantitative profile of the potential vorticity, qualitative changes are useful in detecting different flow regimes. By comparing the downstream changes in total and planetary potential vorticity, one can deduce frictional and inertial regimes in the different layers. The presence of a frictional regime at the inshore edge suggests that care should be taken in assuming that potential vorticity is conserved in western boundary currents. In addition the potential vorticity sections trace a pattern of the recirculation feeding into the Brazil Current in the upper layers; other tracers did not provide a clear image. The final picture which emerges is not of a small, surface-trapped Brazil Current; rather, it is that of a classic western boundary current, increasing in strength and depth before turning east into the interior ocean.

Description: Financial support for the data collection and initial analysis was provided through the Office of Naval Research South Atlantic Accelerated Research proposal under contract N00014-82-C-0019. Continued analysis was supported by the National Science Foundation under grant OCE86-14486.

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 1991

Description: Using a semi-geostrophic, reduced gravity thin jet model, we analytically study the evolution of initial meanders into pinched-off rings. The model used is similar to the path equation developed by Flierl and Robinson {1984) for vertically coherent meanders. However, in the present model, the meanders are baroclinic, and a stretching term arises due to the motion of the interface. It can be shown that the equation governing the time-dependent meander of this jet (Pratt, 1988) can be transformed into the Modified Korteweg- deVries (MKdV) equation in intrinsic coordinates. The MKdV equation admits two types of solitary wave solutions, loop solitons and breathers. The breathers are permanent meanders which propagate on the path , and some are able to form rings. Using the inverse scattering transform , we can predict breather and ring formation for simple initial meanders. The inverse scattering t ransform is applied to S and Ω shaped meanders with piecewise constant and continuous curvature. S shaped meanders, or steps, must be multi-valued to form breathers, and must have very steep angles in order to form rings. Due to integral constraints, Ω shaped meanders, or lobes, are unable to pinch together to form rings unless they are wide enough so that the two side flanks of the lobe act as two independent steps. The numerical solutions indicate that the breathers predicted by the inverse scattering is a very good approximation to the full solution.

Description: This work has been supported by NSF contract number OCE 89-16446 and ONR contract number N00014-S9-J-1182.

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

Description: This thesis examines the effect of mean large-scale currents on the vertical structure of the upper ocean during two recent observational programs: the Long Term Upper Ocean Study (LOTUS) and the TROPIC HEAT experiments. The LOTUS experiment took place in the northwest Atlantic Ocean, a mid-latitude region away from strong mean currents, and extended over one entire seasonal cycle. The TROPIC HEAT experiments took place in the central equatorial Pacific Ocean during two 12-day periods in 1984 and 1987, at opposite extremes of the seasonal cycle. We use observations from these field experiments as well as one-dimensional numerical models of the upper ocean to analyze the dynamics of the vertical structure of the upper ocean at the equator and in mid-latitudes. Due to the different nature of the observations, we focus on the long term mean structure of the upper ocean in the LOTUS observations (Chapters 2 and 3), and on the diurnal cycle in the equatorial upper ocean in our analysis of the TROPIC HEAT observations (Chapters 4 and 5). In the LOTUS observations, we find that the observed current is coherent with the wind over low frequencies (greater than an inertial period). Using a wind-relative averaging method we find good agreement with Ekman transport throughout the first summer and winter of the LOTUS experiment, with the exception of a downwind component in the wintertime. The mean current spiral is flat compared to the classic Ekman spiral, in that it rotates less with depth than does the Ekman spiral. The mean current has an e-folding depth scale of 12m in the summer and 25 min the winter. Diurnal cycling is the dominant variability in the summer and determines the vertical structure of the spiral. In the winter, diurnal cycling is almost non-existent due to greatly reduced solar insolation. There is a persistent downwind shear in the upper 15 m during the winter which may be partially due to a bias induced by surface wave motion but which is also consistent with a logarithmic boundary layer. The Price et al. (1986) model is reasonably successful in simulating the current structure during the summer, capturing both the mean and the diurnal variation. The model is less successful in the winter, though it does capture the overall depth scale of the current spiral. In our analysis of the TROPIC HEAT observations, we extend the Price et al. (1986) model to the equatorial upper ocean. The model is initialized with the stratification and shear of the Equatorial Undercurrent (EUC), and is driven with heating and wind stress. A surface mixed layer is determined by bulk stability requirements, and a transition layer below the mixed layer is simulated by requiring that the gradient Richardson number be no less than 1/4. A principal result is that the nighttime phase of the diurnal cycle is strongly affected by the EUC, resulting in deep mixing and large dissipation at night consistent with observations of the equatorial upper ocean during TROPIC HEAT. Other features of the equatorial circulation (upwelling and the zonal pressure gradient) are of little direct importance to the diurnal cycle. The daytime (heating) phase of the simulated diurnal cycle is unaffected by equatorial circulation and is very similar to its mid-latitude counterpart. Solar heating produces a stably stratified surface layer roughly 10 m thick within which there is little, 0(3 x 10-8 W kg-1), turbulent dissipation. The diurnal stratification, though small compared to the EUC, is sufficient to insulate the EUC from wind stress during the day. For the typical range of conditions at the equator, diurnal warming of the sea surface is 0.2-0.5°C, and the diurnal variation of surface current (diurnal jet) is 0.1-0.2 m s-1, consistent with observations. The nighttime (cooling) phase of the simulated diurnal cycle is quite different from that seen at mid-latitudes. As cooling removes the warm, stable surface layer, the wind stress can work directly against the shear of the EUC. This produces a transition layer that can reach to 80 m depth, or nearly to the core of the EUC. Within this layer the turbulent dissipation is quite large, 0(2 x 10-7 W kg-1). Thus, the simulated dissipation has a diurnal range of more than a factor of five, as observed in the 1984 TROPIC HEAT experiment, though the diurnal cycle of stratification and current are fairly modest. Dissipation estimated from the model is due to wind working directly against EUC, and is similar to observed values of dissipation in both magnitude and depth range. Overall dissipation values in the model are set by the strength of the wind stress rather than the structure of the EUC, and rise approximately like u*3 for a given Undercurrent. This suggests that the lower values of dissipation observed in the 1987 TROPIC HEAT experiment were due to the lower wind stress values rather than the relatively weak Undercurrent. The main findings of this thesis are: 1) When the diurnal cycle in solar heating is strong, it determines the local vertical structure of the upper ocean (in both the LOTUS and TROPIC HEAT observations). The Price et al. (1986) model and its extension to the equator simulate the upper ocean fairly well when the diurnal cycle is strong. Under these conditions it is necessary to make measurements very near the surface ( 〈 10 m depth) to fully resolve the wind-driven flow. 2) When surface waves are strong, surface-moored measurements of current may have a significant wave bias. To accurately estimate this bias, simultaneous measurements of current, current meter motion, and surface waves are needed. 3) Mean currents strongly amplify the nighttime phase of the diurnal cycle in the equatorial upper ocean, and therefore alter the mean structure of the equatorial upper ocean.

Description: The Office of Naval Research supported this work under contract N00014-89-J-1053.

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: Author Posting. © American Geophysical Union, 2012. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 13 (2012): Q04007, doi:10.1029/2011GC003990.

Description: Here we demonstrate with a study of the Lucky Strike hydrothermal field that image mosaicing over large seafloor areas is feasible with new image processing techniques, and that repeated surveys allow temporal studies of active processes. Lucky Strike mosaics, generated from 〉56,000 images acquired in 1996, 2006, 2008 and 2009, reveal the distribution and types of diffuse outflow throughout the field, and their association with high-temperature vents. In detail, the zones of outflow are largely controlled by faults, and we suggest that the spatial clustering of active zones likely reflects the geometry of the underlying plumbing system. Imagery also provides constraints on temporal variability at two time-scales. First, based upon changes in individual outflow features identified in mosaics acquired in different years, we document a general decline of diffuse outflow throughout the vent field over time-scales up to 13 years. Second, the image mosaics reveal broad patches of seafloor that we interpret as fossil outflow zones, owing to their association with extinct chimneys and hydrothermal deposits. These areas encompass the entire region of present-day hydrothermal activity, suggesting that the plumbing system has persisted over long periods of time, loosely constrained to hundreds to thousands of years. The coupling of mosaic interpretation and available field measurements allow us to independently estimate the heat flux of the Lucky Strike system at ~200 to 1000 MW, with 75% to 〉90% of this flux taken up by diffuse hydrothermal outflow. Based on these heat flux estimates, we propose that the temporal decline of the system at short and long time scales may be explained by the progressive cooling of the AMC, without replenishment. The results at Lucky Strike demonstrate that repeated image surveys can be routinely performed to characterize and study the temporal variability of a broad range of vent sites hosting active processes (e.g., cold seeps, hydrothermal fields, gas outflows, etc.), allowing a better understanding of fluid flow dynamics from the sub-seafloor, and a quantification of fluxes.

Description: This project was funded by CNRS/IFREMER through the 2006, 2008, 2009 and 2010 cruises within the MoMAR program (France), by ANR (France) Mothseim Project NT05-3 42213 to J. Escartín, and by grant CTM2010-15216/MAR from the Spanish Ministry of Science to R. Garcia and J. Escartín. T. Barreyre was supported by University Paris Diderot (Paris 7– France) and Institut de Physique du Globe de Paris (IPGP, France). E. Mittelstaedt was supported by the International Research Fellowship Program of the U.S. National Science Foundation (OISE-0757920).

Description: 2012-10-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 1990

Description: The contribution of tropical instability waves to the momentum and energy balances of the Pacific Equatorial Undercurrent is investigated using velocity and temperature time series from the three-dimensional Equatorial Pacific Ocean Climate Study mooring array at 110°W. Tropical instability waves are an energetic band of variability typically with periods between 14 and 36 days which are thought to be generated by instability of the equatorial currents. They are frequently observed as meanders of the equatorial front in satellite sea surface temperature maps. Here, they are observed as large oscillations in the meridional velocity records at l10°W with an energy peak at 21 days. Westward phase propagation is observed in this band with a phase speed of -0.9 (±0.3) m s-1 and a wavelength of 1660 km. Upward phase propagation is observed which is consistent with downward energy propagation. The observed propagation characteristics are compared with those of the mixed Rossby-gravity wave. The variability in this band produces large northward fluxes of eastward momentum and southward fluxes of temperature which affect the dynamics of the mean Undercurrent through the Reynolds stress divergence, and the Eliassen-Palm flux divergence. The waves produce a northward flux of eastward momentum, uv, which is largest at the northern mooring in the upper part of the array. The meridional divergence of eastward momentum, -δ(uv)/δy, decelerates the Undercurrent core down to 150 m. This implies a coupling between the Undercurrent and the South Equatorial Current with the eastward momentum of the Undercurrent transferred to the westward flowing South Equatorial Current. To estimate the vertical momentum flux divergence, the vertical eddy flux of eastward momentum, uw, is inferred using the eddy temperature equation. The vertical eddy momentum flux is positive and largest at the core of the Undercurrent, implying an acceleration of the eastward flow above the core and a deceleration below. The Eliassen-Palm flux divergence is small above the core of the Undercurrent at 75 m, but below the core, is sufficient to balance the deeply penetrating eastward pressure gradient force. The instability waves are important to the energetics of the mean Undercurrent. An exchange of kinetic energy from the mean Undercurrent to the waves through shear production is estimated. A local exchange is suggested since the rate at which the mean Undercurrent loses kinetic energy through instability is comparable to the rate at which the waves gain energy through shear production. The conversion from mean to eddy potential energy is an order of magnitude smaller with the waves gaining potential energy through conversion of mean available potential energy. The observations of upward phase propagation and downward Eliassen-Palm flux suggest that the waves propagate energy downward into the deep ocean. The energetics and momentum balance of the mean Undercurrent is investigated further by analyzing the downstream change in the Bernoulli function on the equator along isentropes or potential density surfaces using mean hydrographic sections at 150°W and 110°W. A downstream decrease in the Bernoulli function is observed which is due to a decrease in the Acceleration Potential since the mean kinetic energy of the Undercurrent changes little from 150°W to 110°W. The lateral divergence of eddy momentum fluxes calculated on isotherms is sufficient to balance the observed decrease in the Acceleration Potential. The downstream decrease in the Acceleration potential has further implications for the mean energetics since this "downhill" flow releases mean available potential energy stored in the east-west sloping thermocline. The rate at which the Undercurrent releases available potential energy, is shown to be comparable to the rate at which the mean flow loses kinetic energy by interaction with the waves, with the waves gaining kinetic energy in the process. Thus, it is hypothesized that in the eastern Pacific this downstream release of available potential energy is ultimately converted into a downstream increase in the kinetic energy of the waves rather than the kinetic energy of the mean flow as occurs in the western Pacific. To maintain an equilibrium, the waves radiate energy into the deep ocean as is suggested by the upward phase propagation and the downward Eliassen-Palm flux.

Description: Financial support of the National Science Foundation under contracts OCE 82-14955 and OCE 85-19551, for participation in Tropic Heat, and OCE 85-04125.

Description: Author Posting. © American Geophysical Union, 2011. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 12 (2011): Q05008, doi:10.1029/2010GC003439.

Description: Variations in topography and seismic structure are observed along the Juan de Fuca (JdF) Ridge axis and in the vicinity of pseudofaults on the ridge flanks left by former episodes of ridge propagation. Here we analyze gravity data coregistered with multichannel seismic data from the JdF Ridge and flanks in order to better understand the origin of crustal structure variations in this area. The data were collected along the ridge axis and along three ridge-perpendicular transects at the Endeavor, Northern Symmetric, and Cleft segments. Negative Mantle Bouguer anomalies of −21 to −28 mGal are observed at the axis of the three segments. Thicker crust at the Endeavor and Cleft segments is inferred from seismic data and can account for the small differences in axial gravity anomalies (3–7 mGal). Additional low densities/elevated temperatures within and/or below the axial crust are required to explain the remaining axial MBA low at all segments. Gravity models indicate that the region of low densities is wider beneath the Cleft segment. Gravity models for pseudofaults crossed along the three transects support the presence of thinner and denser crust within the pseudofault zones that we attribute to iron-enriched crust. On the young crust side of the pseudofaults, a 10–20 km wide zone of thicker crust is found. Reflection events interpreted as subcrustal sills underlie the zones of thicker crust and are the presumed source for the iron enrichment.

Description: This work was supported by the National Science Foundation grants OCE‐0648303 to Lamont‐Doherty Earth Observatory, OCE‐0648923 to Woods Hole Oceanographic Institution.

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 studies the problems of generation and maintenance of recirculations by Gulf Stream instabilities. Observations show that the horizontal structure of the jet and its recirculations suffer significant changes in time. Here, the role of internal dynamics of the jet is isolated as one of the possible sources of such variability, and the differences between barotropic and baroclinic instabilities are investigated. The problem of recirculation development is considered in a framework of a free spin down of the 2-layer and the 1-layer, zonally symmetric, quasi-geostrophic jets. Linear stability analysis shows that in strongly baroclinic basic flows, eddies are capable of driving recirculations in the lower layer through the residual meridional circulation. In strongly barotropic jets, the linearly most unstable wave simply diffuses the jet. Nonlinear stability analysis indicates that recirculations are robust features of the 2-layer model. The strength of recirculations is a function of the model’s parameters. It increases with a decrease in the value of the nondimensional /3 due to potential vorticity homogenization constrained by enstrophy conservation. The recirculation strength is a non-monotonic function of the baroclinic velocity parameter; it is the strongest for strongly baroclinic basic flows, weakest for flows with intermediate baroclinic structure and of medium strength for strongly barotropic basic flows. Such non-monotonic behavior is the result of two different processes responsible for the recirculation development: linear eddy-mean flow interactions for strongly baroclinic basic flows and strongly nonlinear eddy-eddy and eddy-mean flow interaction for strongly barotropic flows. In the case of the reduced-gravity model, recirculations develop only for infinite deformation raduis. Basic flows with finite deformation radius are only weakly supercritical and therefore produced negligible recirculations after equilibration. The problem of maintenance of the recirculations is coupled to the questions of existence of low frequency variability and of multiple dynamical regimes of a system consisting of a quasi-geostrophic jet and its recirculations. The problem is studied in a framework of a 2-layer or a reduced-gravity colliding jets model which has no windforcing. Instead, it is forced by inflows and outflows through the open boundaries. Oniy the western boundary of the domain is closed, and a free slip boundary condition is used there. The results of the numerical experiments show that when oniy the mechanism of barotropic instability is present, the model has two energy states for a wide range of interfacial friction coefficients. The high energy state is characterized by well-developed recirculations and displays strong variability associated with either large recirculating gyres and a weak eddy field or small recirculations and a strong eddy field. The iow energy state is characterized by large meridional excursions in the separation point and large amplitude, westward propagating meanders that produce strong rings after interacting with the western wall. For physically relevant bottom friction values, the presence of baroclinic in stability in the recirculation regions of the 2-layer model allows for a unique dynamical regime characterized by well-developed recirculations in both layers. The low-frequency variability associated with the regime is weak and is related to meridional shifts in the position of the jet, to wrapping of the recirculations around each other, and to pulsations in their zonal extent. For strong bottom friction, the 2-layer model has only the mechanism of barotropic instability which reduces it to a 1 1/2-layer configuration; the model displays two dynamical regimes and strong low frequency variability in the upper layer, while the lower layer is strongly frictional.

Description: Financial support for this research was provided by NSF grant number OCE 9617848.

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

Description: Remote sensing and in situ observations are used to investigate the ocean response to the Tokar Wind Jet in the Red Sea. The wind jet blows down the pressure gradient through the Tokar Gap on the Sudanese coast, at about 18°N, during the summer monsoon season. It disturbs the prevailing along-sea (southeastward) winds with strong cross-sea (northeastward) winds that can last from days to weeks and reach amplitudes of 20-25 m/s. By comparing scatterometer winds with along-track and gridded sea level anomaly observations, it is shown that an intense dipolar eddy spins up in less than seven days in response to the wind jet. The eddy pair has a horizontal scale of 140 km. Maximum ocean surface velocities can reach 1 m/s and eddy currents extend at least 200 m into the water column. The eddy currents appear to cover the width of the sea, providing a pathway for rapid transport of marine organisms and other drifting material from one coast to the other. Interannual variability in the strength of the dipole is closely matched with variability in the strength of the wind jet. The dipole is observed to be quasi-stationary, although there is some evidence for slow eastward propagation—simulation of the dipole in an idealized high-resolution numerical model suggests that this is the result of self-advection. These and other recent in situ observations in the Red Sea show that the upper ocean currents are dominated by mesoscale eddies rather than by a slow overturning circulation.

Description: This work is supported by Award Nos. USA 00002, KSA 00011 and KSA 00011/02 made by King Abdullah University of Science and Technology (KAUST).

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: Observations of diapycnal mixing rates are examined and related to diapycnal advection for both double-diffusive and turbulent regimes. The role of double-diffusive mixing at the site of the North Atlantic Tracer Release Experiment is considered. The strength of salt-finger mixing is analyzed in terms of the stability parameters for shear and double-diffusive convection, and a nondimensional ratio of the thermal and energy dissipation rates. While the model for turbulence describes most dissipation occurring in high shear, dissipation in low shear is better described by the salt-finger model, and a method for estimating the salt-finger enhancement of the diapycnal haline diffusivity over the thermal diffusivity is proposed. Best agreement between tracer-inferred mixing rates and microstructure based estimates is achieved when the salt-finger enhancement of haline flux is taken into account. The role of turbulence occurring above rough bathymetry in the abyssal Brazil Basin is also considered. The mixing levels along sloping bathymetry exceed the levels observed on ridge crests and canyon floors. Additionally, mixing levels modulate in phase with the spring-neap tidal cycle. A model of the dissipation rate is derived and used to specify the turbulent mixing rate and constrain the diapycnal advection in an inverse model for the steady circulation. The inverse model solution reveals the presence of a secondary circulation with zonal character. These results suggest that mixing in abyssal canyons plays an important role in the mass budget of Antarctic Bottom Water.

Description: This work was supported by contracts N00014-92-1323 and N00014-97-10087 of the Office of Naval Research and grant OCE94-15589 of 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 January 1988

Description: Analysis of vertical profiles of absolute horizontal velocity collected in January 1981, February 1982 and April 1982 in the central equatorial Pacific as part of the Pacific Equatorial Ocean Dynamics (PEQUOD) program, revealed two significant narrow band spectral peaks in the zonal velocity records, centered at vertical wavelengths of 560 and 350 stretched meters (sm). Both signals were present in all three cruises, but the 350 sm peak showed a more steady character in amplitude and a higher signal-to-noise ratio. In addition, its vertical scales corresponded to the scales of the conspicuous alternating flows generically called the equatorial deep jets in the past (the same terminology will be used here). Meridional velocity and vertical displacement spectra did not show any such energetic features. Energy in the 560 sm band roughly doubled between January 1981 and April 1982. Time lagged coherence results suggested upward phase propagation at time scales of about 4 years. East-west phase lines computed from zonally lagged coherences, tilted downward towards the west, implying westward phase propagation. Estimates of zonal wavelength (on the order of 10000 km) and period based on these coherence calculations, and the observed energy meridional structure at this vertical wavenumber band, seem consistent, within experimental errors, with the presence of a first meridional mode long Rossby wave packet, weakly modulated in the zonal direction. The equatorial deep jets, identified with the peak centered at 350 sm, are best defined as a finite narrow band process in vertical wavenumber (311-400 sm), accounting for only 20% of the total variance present in the broad band energetic background. At the jets wavenumber band, latitudinal energy scaling compared well with Kelvin wave theoretical values and a general tilt of phase lines downward towards the east yielded estimates of 10000-16000 km for the zonal wavelengths. Time-lagged coherence calculations revealed evidence for vertical shifting of the jets on interannual time scales. Interpretation of results in terms of single frequency linear wave processes led to inconsistencies, but finite bandwidth (in frequency and wavenumber) Kelvin wave processes of periods on the order of three to five years could account for the observations. Thus, the records do not preclude equatorial waves as a reasonable kinematic description of the jets.

Description: This research was supported by grant OCE-8600052 from the National Science Foundation, through the Woods Hole Oceanographic Institution.

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

Description: Hurricane activity in the Northeastern Gulf of Mexico and its relationship to regional and large-scale climate variability during the Late Holocene is explored. A 4500-year record of hurricane-induced storm surges is developed from sediment cores collected from a coastal sinkhole near Apalachee Bay, Florida. Reconstructed hurricane frequency is shown to exhibit statistically significant variability with the greatest activity occurring between 2700 and 2400 years ago and the least activity between 1900 to 1600 years ago and after 600 years ago. Proxy records of stormrelevant climate variables contain similar timescales of variability and suggest both regional and large-scale mechanisms have influenced hurricane activity on centennial to millennial timescales. In particular, low-frequency migrations of the Loop Current may exercise control over regional hurricane activity by changing the thermal structure of the upper ocean and influencing the role of storm-induced upwelling on hurricane intensification. A new method for estimating the frequency of hurricanegenerated storm surges is presented and applied to Apalachee Bay, Florida. Multisite paleohurricane reconstructions from this region are developed, and the effects of geographic boundary conditions and temporal resolution on estimates of paleohurricane frequency are explored.

Description: Financial support provided by the American Meteorological Society, the National Science Foundation, the Bermuda Risk Prediction Initiative, the National Center for Airborne Laser Mapping, and the Coastal Ocean Institute.

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 1998

Description: The upper ocean crust contains a comprehensive record of the shallow geological processes active along the world's mid-ocean ridge system. This thesis examines the magnetic and seismic structure of the upper crust at two contrasting ridges-the fast spreading East Pacific Rise (EPR) and the slow spreading Mid-Atlantic Ridge (MAR)-to build a more complete understanding about the roles of volcanic emplacement, tectonic disruption and hydrothermal alteration in the near-ridge environment. A technique that inverts potential field measurements. directly from an uneven observation track is developed and applied to near-bottom magnetic data from the spreading segments north of the Kane transform on the MAR. It is concluded that the central anomaly magnetization high marks the locus of focused volcanic emplacement. A cyclic faulting model is proposed to explain the oscillatory magnetization pattern associated with discrete blocks of crust being transported out of the rift valley between intensely altered fault zones. Seismic waveform and amplitude analyses of the magma sill along the EPR reveal it to be a thin (〈100 m) body of partial melt. These characteristics have important implications for melt availability and transport within the cycle of eruption and replenishment. A genetic algorithm-based seismic waveform inversion technique is developed and applied to on- and near-axis multichannel data from 17°20'S on the EPR and the spreading segment south of the Oceanographer transform (MAR) to map and compare for the first time the detailed velocity structure of the upper crust at two different spreading rates. Combined with conventionally processed seismic profiles, our results show that, while final extrusive thickness is comparable at all spreading ridges (300-500 m), the style of thickening may vary. While a thin (≤100 m) extrusive carapace quadruples in thickness within 1-4 km of the EPR crest, the extrusive section at the MAR achieves its final thickness within the inner valley. Both show evidence for a narrow zone of volcanic emplacement. Vigorous hydrothermalism at the EPR may produce a more rapid increase in basement velocities relative to the MAR. Rapid modification of the extrusive/dike transition at both ridges indicates that hydrothermalism is enhanced in this interval. Along-axis transport of lavas may thicken the extrusive pile at slow spreading segment ends, strengthening the magnetic highs generated by lava chemistry.

Description: ONR graduate fellow. The magnetics portion of this thesis was also supported by NSF grants OCE-9204141, OCE-9200905 (M. A. Tivey & H. Schouten) and NERC GR3/7702 (R. C. Searle). Research for Chapter 4 was partially funded by NSF grant OCE-9402933 (R. S. Detrick). The final two science chapters were supported by NSF grants OCE-9012707, OCE-9300450 (R. S. Detrick), OCE-9401717 (G. M. Kent & R. S. Detrick), OCE-9400623 (M. A. Tivey) and the Education office.

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

Description: This study focuses on the zonal weakening, eastern termination and seasonal variations of the Atlantic equatorial undercurrent (EUC). The main and most original contribution of the dissertation is a detailed analysis of the Atlantic EUC simulated by Philander and Pacanowski's (1986) general circulation model (GCM), which provides a novel description of the dynamical regimes governing various regions of a nonlinear stratified undercurrent. Only in a narrow deep western region of the simulation does one find an approximately inertial regime corresponding to zonal acceleration. Elsewhere frictional processes cannot be ignored. The bulk of the mid-basin model EUC terminates in the overlying westward surface flow while only a small fraction (the deeper more inertial layers) terminates at the eastern coast. In agreement with observations, a robust feature of the GCM not present in simpler models is the apparent migration of the EUC core from above the thermocline in the west to below it in the east. In the GCM, this happens because the eastward flow is eroded more efficiently by vertical friction above the base of the thermocline than by lateral friction at greater depths. This mechanism is a plausible one for the observed EUC. A scale analysis using a depth scale which decreases with distance eastwards predicts the model zonal transition between western inertial and eastern inertio-frictional regimes. Historical and recent observations and simple models of the equatorial and coastal eastern undercurrents are reviewed, and a new analysis of current measurements in the eastern equatorial Atlantic is presented. Although the measurements are inadequate for definitive conclusions, they suggest that Lukas' (1981) claim of a spring surge of the Pacific EUC to the eastern coast and a seasonal branching of the EUC into a coastal southeastward undercurrent may also hold for the Atlantic Ocean. To improve the agreement between observed and modelled strength of the eastern undercurrent, it is suggested that the eddy coefficient of horizontal mixing should be reduced in future GCM simulations.

Description: This work was supported by NSF grants OCE82-14771, OCE82-08744 and OCE85-14885.

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

Description: In previous Gulf Stream work (Hall and Bryden, 1985, Hall, 1985, 198GA, 198GB), a decomposition of multiple depth current records was developed which produced along- and cross-stream components. The cross-stream component was found to occasionally match lateral displacements of the Stream, as determined by temperature changes measured at the current meters. This study determined where within the meander pattern of the Gulf Stream the cross-stream velocity calculated from current meters at depth correctly predicted translations of the Gulf Stream as measured by satellite data. Additionally, the effects of recently quantified cross-stream velocities associated with the curvature of Gulf Stream meanders were analyzed.

Description: Funds for this work were provided by ONR contracts N00014-86-K-0751 and N00014-87-K-0007.

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 1988

Description: The highly energetic Agulhas Retroflection region south of the African continent lies at the junction of the South Indian, South Atlantic, and Circumpolar Oceans. A new survey of the Agulhas Retroflection taken in March 1985, plus historical hydrographic data, allow its dynamical and water-mass characteristics, and its role in exchanging mass, tracers, and vorticity between the three oceans, to be extensively characterized. The 1985 survey is composed of three independent, synoptic elements: a grid of closely-spaced, full-water-depth hydrographic stations (the first entirely full-water-column survey in this area), including several transects of the Agulhas and Agulhas Return Currents; a continuous survey of the path of the currents (the first such survey in the Agulhas); and a contemporaneous and relatively cloud-free sea surface temperature image derived from satellite infrared measurements. Mass transport balances within the closed grid boxes of the 1985 hydrographic survey provide information about current transport, recirculation (transport in excess of estimated returning interior ocean transport), and the overall Retroflection transport pattern. The current transport values exceed by as much as a factor of 1.5 the maximum interior transport computed from observed wind-stress curl and linear theory. Agulhas Current transports ranged from 56 to 95 x 106 m s-l at four 1985 transects crossing the current. Agulhas Return Current transports at the two 1985 transects were 54 and 65 x 106 m s-l. These transports are computed relative to 2400 dbar, which lies below the deep oxygen minimum emanating from the South Indian Ocean, and above the North Atlantic Deep Water salinity maximum. The current retroflected in two distinct branches in 1985, with a cold ring and a partially isolated warm recirculation cell found between the two branches. The satellite-derived sea surface temperature (SST) image, in agreement with the in situ measurements, showed that the cold ring lacked a cold SST anomaly; that the subsurface current path, as represented by a survey of the 15 C isotherm and 200 dbar surface intersection, was closely followed by a sharp front in sea surface temperature; and that most of the Agulhas's surface warm core retroflected upstream of the second retroflection branch. Anticyclonic curvature vorticity at sharp turns in the subsurface current path was found to exceed the maximum allowed by gradient wind balance, indicating that at these locations time-dependence and cross-frontal flow are important. The current's density field is found to meet necessary conditions for baroclinic and barotropic instability. These instability mechanisms may play a role in ring formation and current meandering. Top-to-bottom cross-stream spatial and isopycnal water-mass layering in the Agulhas Current, Agulhas Return Current, and associated rings are presented in two sets of sections, one contoured with pressure and the other with potential density as vertical coordinate. Temperature, salinity, oxygen, potential density and velocity sections are shown contoured versus pressure; and pressure, salinity, oxygen, and planetary potential vorticity are shown contoured versus potential density. These sections clearly illustrate water-mass structure both in space and relative to isopycnal surfaces. Strong salt, oxygen, and potential vorticity fronts on isopycnals in the upper -300 m across the Agulhas and Agulhas Return Current are observed, as are deep western boundary filaments of (i) salty, low oxygen water at intermediate depths traceable to Red Sea Water influences, and (ii) salty North Atlantic Deep Water close round the tip of Africa. The 1985 cold-core ring is the first cold-cored isolated feature to be observed within the Retroflection itself. Its transport was 64 x 106 m s-1, its integrated kinetic and available potential energy anomalies were 8.3 and 61 x 1015 J respectively, and its integrated planetary potential vorticity anomaly was 2.8 x 10-12 m-1 s-1. The potential vorticity flux associated with the exchange of 25 warm ring/cold ring pairs per year between the South Indian and Southern Oceans would balance the potential vorticity input by the wind to the entire South Indian Ocean. Interbasin flow of warm thermocline water (warmer than 8 C) from the South Indian to the South Atlantic Ocean is reconsidered in light of the 1985 hydrographic data. Thermocline water flow from the South Indian Ocean into the South Atlantic in the 1985 and historical observations is found to range from 2.8 to 〈9.6 x 106 m s-I. These values are less than the S;10 x 106 m s·1 needed to balance the Atlantic Ocean export of deep water, and implies that the deep water export is balanced in part by water colder than 8 C.

Description: Funding was provided by the Office of Naval Research under contract numbers NOOOI4-84-C-OI34 (NR083-400), NOOOI4-85-C-OOOl (NR083-004), and NOOOI4-87-K-0007 (NR083-004).

Description: Author Posting. © American Geophysical Union, 2011. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 38 (2011): L14307, doi:10.1029/2011GL047798.

Description: The heat output and thermal regime of fast and intermediate spreading centers are strongly controlled by boundary layer processes between the hydrothermal system and the underlying crustal magma chamber (AMC), which remain to be fully understood. Here, we model the interactions between a shallow two-dimensional cellular hydrothermal system at temperatures 〈700°C, and a deeper AMC at temperatures up to 1200°C. We show that hydrothermal cooling can freeze the AMC in years to decades, unless melt injections occur on commensurate timescales. Moreover, the differential cooling between upflow and downflow zones can segment the AMC into mush and melt regions that alternate on sub-kilometric length scales. These predictions are consistent with along-axis variations in AMC roof depth observed in ophiolites and oceanic settings. In this respect, fine-scale geophysical investigations of the structure of AMCs may help constrain hydrothermal recharge locations associated with active hydrothermal sites.

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 2011

Description: Throughout this thesis we will discuss the theoretical background and empirical observation of a swell band shore normal flux divergence reversal. Specifically, we will demonstrate the existence and persistence of the energy flux divergence reversal in the nearshore region of Atchafalaya Bay, Gulf of Mexico, across storms during the March through April 2010 deployment. We will show that the swell band offshore component of energy flux is rather insignificant during the periods of interest, and as such we will neglect it during the ensuing analysis. The data presented will verify that the greatest flux divergence reversal is seen with winds from the East to Southeast, which is consistent with theories which suggest shoreward energy flux as well as estuarine sediment transport and resuspension prior to passage of a cold front. Employing the results of theoretical calculations and numerical modeling we will confirm that a plausible explanation for this phenomena can be found in situations where temporally varying wind input may locally balance or overpower bottom induced dissipation, which may also contravene the hypothesis that dissipation need increase shoreward due to nonlinear wave-wave interactions and maturation of the spectrum. Lastly, we will verify that the data presented is consistent with other measures collected during the same deployment in the Atchafalaya Bay during March - April 2010.

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 1999

Description: The thesis investigates the circulation at a 76-m deep study site on the southern flank of Georges Bank, a shallow submarine bank located between the deeper Gulf of Maine and the continental slope. Emphasis is placed on the vertical structure of the bottom boundary layer driven by the semi diurnal tides and the flow field's response to wind forcing. The observational analysis presented here is based on a combination of moored array and bottom tripod-mounted current, temperature, conductivity, and meteorological data taken between February and August 1995. Results from the bottom boundary layer analysis are compared to numerical model predictions for tidal flow over rough bottom topography. The flow response to wind forcing is examined and brought into context with the existing understanding of the wind-induced circulation in the Georges Bank region. Particular attention is given to the vertical distribution of the wind-driven currents, whose variation with background stratification is discussed and compared to observations from open ocean studies.

Description: The research presented in this thesis was generously supported by the National Science Foundation under grants OCE 93-13671 and OCE 96-32357 as part of the U.S. GLOBEC/Georges Bank Program.

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

Description: A state-of-the-art, high-resolution ocean general circulation model is used to estimate the time-dependent global ocean heat transport and investigate its dynamics. The north-south heat transport is the prime manifestation of the ocean’s role in global climate, but understanding of its variability has been fragmentary owing to uncertainties in observational analyses, limitations in models, and the lack of a convincing mechanism. These issues are addressed in this thesis. Technical problems associated with the forcing and sampling of the model, and the impact of high-frequency motions are discussed. Numerical schemes are suggested to remove the inertial energy to prevent aliasing when the model fields are stored for later analysis. Globally, the cross-equatorial, seasonal heat transport fluctuations are close to +4.5 x 1015 watts, the same amplitude as the seasonal, cross-equatorial atmospheric energy transport. The variability is concentrated within 200 of the equator and dominated by the annual cycle. The majority of it is due to wind-induced current fluctuations in which the time-varying wind drives Ekman layer mass transports that are compensated by depth-independent return flows. The temperature difference between the mass transports gives rise to the time-dependent heat transport. The rectified eddy heat transport is calculated from the model. It is weak in the central gyres, and strong in the western boundary currents, the Antarctic Circumpolar Current, and the equatorial region. It is largely confined to the upper 1000 meters of the ocean. The rotational component of the eddy heat transport is strong in the oceanic jets, while the divergent component is strongest in the equatorial region and Antarctic Circumpolar Current. The method of estimating the eddy heat transport from an eddy diffusivity derived from mixing length arguments and altimetry data, and the climatological temperature field, is tested and shown not to reproduce the model’s directly evaluated eddy heat transport. Possible reasons for the discrepancy are explored.

Description: Funding for this research came from the Department of Defense under a National Defense Science and Engineering Graduate Fellowship. Financial support was also contributed by the National Science Foundation through grants #OCE-9617570 and #OCE-9730071, and the Tokyo Electric Power Company through the TEPCO/MIT Environmental Research Program. The author received partial support from an MIT Climate Modeling Fellowship, made possible by a gift from the American Automobile Manufacturers Association.

Description: Author Posting. © American Geophysical Union, 2011. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 116 (2011): B10102, doi:10.1029/2011JB008259.

Description: Crustal thickness variations at the ultraslow spreading 10–16°E region of the Southwest Indian Ridge are used to constrain melt migration processes. In the study area, ridge morphology correlates with the obliquity of the ridge axis with respect to the spreading direction. A long oblique “supersegment”, nearly devoid of magmatism, is flanked at either end by robust magmatic centers (Joseph Mayes Seamount and Narrowgate segment) of much lesser obliquity. Plate-driven mantle flow and temperature structure are calculated in 3-D based on the observed ridge segmentation. Melt extraction is assumed to occur in three steps: (1) vertical migration out of the melting region, (2) focusing along an inclined permeability barrier, and (3) extraction when the melt enters a region shallower than ∼35 km within 5 km of the ridge axis. No crust is predicted in our model along the oblique supersegment. The formation of Joseph Mayes Seamount is consistent with an on-axis melt anomaly induced by the local orthogonal spreading. The crustal thickness anomaly at Narrowgate results from melt extracted at a tectonic damage zone as it travels along the axis toward regions of lesser obliquity. Orthogonal spreading enhances the Narrowgate crustal thickness anomaly but is not necessary for it. The lack of a residual mantle Bouguer gravity high along the oblique supersegment can be explained by deep serpentization of the upper mantle permissible by the thermal structure of this ridge segment. Buoyancy-driven upwelling and/or mantle heterogeneities are not required to explain the extreme focusing of melt in the study area.

Description: This work was supported by grants OCE‐ 0623188 and OCE‐0937277 from the National Science Foundation.

Description: Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 35 (2008): L03402, doi:10.1029/2007GL032837.

Description: Arctic rivers transport huge quantities of dissolved organic carbon (DOC) to the Arctic Ocean. The prevailing paradigm is that DOC in arctic rivers is refractory and therefore of little significance for the biogeochemistry of the Arctic Ocean. We show that there is substantial seasonal variability in the lability of DOC transported by Alaskan rivers to the Arctic Ocean: little DOC is lost during incubations of samples collected during summer, but substantial losses (20–40%) occur during incubations of samples collected during the spring freshet when the majority of the annual DOC flux occurs. We speculate that restricting sampling to summer may have biased past studies. If so, then fluvial inputs of DOC to the Arctic Ocean may have a much larger influence on coastal ocean biogeochemistry than previously realized, and reconsideration of the role of terrigenous DOC on carbon, microbial, and food-web dynamics on the arctic shelf will be warranted.

Description: This material is based on work supported by the National Science Foundation under grant numbers OPP-0436106, OPP- 0519840, and EAR-0403962, and is a contribution to the Study of Environmental Arctic Change (SEARCH).

Description: Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 35 (2008): L08606, doi:10.1029/2008GL033532.

Description: Turbulent-scale temperature and conductivity were measured during the pan-arctic Beringia 2005 Expedition. The rates of dissipation of thermal variance and diapycnal diffusivities are calculated along a section from Alaska to the North Pole, across deep flat basins (Canada and Makarov Basins) and steep ridges (Alpha-Mendeleev and Lomonosov Ridges). The mixing rates are observed to be small relative to lower latitudes but also remarkably non-uniform. Relatively elevated turbulence is found over deep topography, confirming the dominant role of bottom-generated internal waves. Measured patterns of mixing in the Arctic are also associated with other mechanisms, such as double-diffusive structures and deep overflows. A better knowledge of the distribution of mixing is essential to understand the dynamics of the changing Arctic environment.

Description: This work was funded by the National Science Foundation through a Small Grant for Exploratory Research (ARC-0527874) and grant ARC-0612342 with additional support from the Doherty Foundation and internal WHOI Funds.

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

Description: Interaction between the Antarctic Circumpolar Current and the continental slope/shelf in the Marguerite Bay and west Antarctic Peninsula is examined as interaction between a wind-driven channel flow and a zonally uniform slope with a bay-shaped shelf to the south. Two control mechanisms, eddy advection and propagation of topographic waves, are identified in barotropic vortex-escarpment interactions. The two mechanisms advect the potential vorticity (PV) perturbations in opposite directions in anticyclone-induced interactions but in the same direction in cyclone-induced interactions, resulting in dramatic differences in the two kinds of interactions. The topographic waves become more nonlinear near the western(eastern if in the Northern Hemisphere) boundary of the bay, where strong cross-escarpment motion occurs. In the interaction between a surface anticyclone and a slope penetrating into the upper layer in a two-layer isopycnal model, the eddy advection decays on length scales on the order of the internal deformation radius, so shoreward over a slope that is wider than the deformation radius, the wave mechanism becomes noticeably significant. It acts to spread the cross-isobath transport in a much wider range while the transport directly driven by the anticyclone is concentrated in space. A two-layer wind-driven channel flow is constructed to the north of the slope in the Southern Hemisphere, spontaneously generating eddies through baroclinic instability. A PV front forms in the first layer shoreward of the base of the topography due to the lower-layer eddy-slope interactions. Perturbed by the jet in the center of the channel, the front interacts with the slope/shelf persistently yet episodically, driving a clockwise mean circulation within the bay as well as crossisobath transport. Both the transports across the slope edge and out of the bay are comparable with the maximum Ekman transport in the channel, indicative of the significance of the examined mechanism. The wave-boundary interaction identified in the barotropic model is found essential for the out-of-bay transport and responsible for the heterogeneity of the transport within the bay. Much more water is transported out of the bay from the west than from the east, and the southeastern area is the most isolated region. These results suggest that strong out-of-bay transport may be found near the western boundary of the Marguerite Bay while the southeastern region is a retention area where high population of Antarctic krill may be found.

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 Geophysical Union, 2007. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 112 (2007): C05011, doi:10.1029/2006JC003899.

Description: In September 2004 a detailed physical and chemical survey was conducted on an anticyclonic, cold-core eddy located seaward of the Chukchi Shelf in the western Arctic Ocean. The eddy had a diameter of ∼16 km and was centered at a depth of ∼160 m between the 1000 and 1500 m isobaths over the continental slope. The water in the core of the eddy (total volume of 25 km3) was of Pacific origin, and contained elevated concentrations of nutrients, organic carbon, and suspended particles. The feature, which likely formed from the boundary current along the edge of the Chukchi Shelf, provides a mechanism for transport of carbon, oxygen, and nutrients directly into the upper halocline of the Canada Basin. Nutrient concentrations in the eddy core were elevated compared to waters of similar density in the deep Canada Basin: silicate (+20 μmol L−1), nitrate (+5 μmol L−1), and phosphate (+0.4 μmol L−1). Organic carbon in the eddy core was also elevated: POC (+3.8 μmol L−1) and DOC (+11 μmol L−1). From these observations, the eddy contained 1.25 × 109 moles Si, 4.5 × 108 moles NO3 −, 5.5 × 107 moles PO3 −, 1.2 × 108 moles POC, and 1.9 × 109 moles DOC, all available for transport to the interior of the Canada Basin. This suggests that such eddies likely play a significant role in maintaining the nutrient maxima observed in the upper halocline. Assuming that shelf-to-basin eddy transport is the dominant renewal mechanism for waters of the upper halocline, remineralization of the excess organic carbon transported into the interior would consume 6.70 × 1010 moles of O2, or one half the total oxygen consumption anticipated arising from all export processes impacting the upper halocline.

Description: This work was supported by the National Science Foundation, and office of Naval Research; DH OPP-0124900, NB OPP-0124868, DK OPP 0124872, RP N00014-02-1-0317.

Description: Author Posting. © American Geophysical Union, 2007. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 112 (2007): G04S60, doi:10.1029/2006JG000371.

Description: Export of nitrate and dissolved organic carbon (DOC) from the upper Kuparuk River between the late 1970s and early 2000s was evaluated using long-term ecological research (LTER) data in combination with solute flux and catchment hydrology models. The USGS Load Estimator (LOADEST) was used to calculate June–August export from 1978 forward. LOADEST was then coupled with a catchment-based land surface model (CLSM) to estimate total annual export from 1991 to 2001. Simulations using the LOADEST/CLSM combination indicate that annual nitrate export from the upper Kuparuk River increased by ~5 fold and annual DOC export decreased by about one half from 1991 to 2001. The decrease in DOC export was focused in May and was primarily attributed to a decrease in river discharge. In contrast, increased nitrate export was evident from May to September and was primarily attributed to increased nitrate concentrations. Increased nitrate concentrations are evident across a wide range of discharge conditions, indicating that higher values do not simply reflect lower discharge in recent years but a significant shift to higher concentration per unit discharge. Nitrate concentrations remained elevated after 2001. However, extraordinarily low discharge during June 2004 and June–August 2005 outweighed the influence of higher concentrations in determining export during these years. The mechanism responsible for the recent increase in nitrate concentrations is uncertain but may relate to changes in soils and vegetation associated with regional warming. While changes in nitrate and DOC export from arctic rivers reflect changes in terrestrial ecosystems, they also have significant implications for Arctic Ocean ecosystems.

Description: This work was supported by the Arctic System Science Program of the National Science Foundation (OPP- 0436118) and by NSF funding for the Arctic LTER through a series of grants from 1987 to present.

Description: Author Posting. © American Geophysical Union, 2007. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 112 (2007): C04S06, doi:10.1029/2006JC003643.

Description: A three-dimensional coupled ocean/ice model, intended for long-term Arctic climate studies, is extended to include tidal effects. From saved output of an Arctic tides model, we introduce parameterizations for (1) enhanced ocean mixing associated with tides and (2) the role of tides fracturing and mobilizing sea ice. Results show tides enhancing loss of heat from Atlantic waters. The impact of tides on sea ice is more subtle as thinning due to enhanced ocean heat flux competes with net ice growth during rapid openings and closings of tidal leads. Present model results are compared with an ensemble of nine models under the Arctic Ocean Model Intercomparison Project (AOMIP). Among results from AOMIP is a tendency for models to accumulate excessive Arctic Ocean heat throughout the intercomparison period 1950 to 2000 which is contrary to observations. Tidally induced ventilation of ocean heat reduces this discrepancy.

Description: This research is supported by the National Science Foundation Office of Polar Programs under cooperative agreements OPP-0002239 and OPP-0327664 with the International Arctic Research Center, University of Alaska Fairbanks.

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 2001

Description: This thesis studies interactions between mid-ocean ridges and mantle plumes using geophysics, geochemistry, and geodynamical modeling. Chapter 1 investigates the effects of the Marion and Bouvet hotspots on the ultra-slow spreading, highly-segmented Southwest Indian Ridge (SWIR). Gravity data indicate that both Marion and Bouvet impart high-amplitude mantle Bouguer anomaly lows to the ridge axis, and suggest that long-offset transforms may diminish along-axis plume flow. Building upon this observation, Chapter 2 presents a series of 3D numerical models designed to quantify the sensitivity of along-axis plume-driven mantle flow to transform offset length, spreading rate, and mantle viscosity structure. The calculations illustrate that long-offset transforms in ultra-slow spreading environments may significantly curtail plume dispersion. Chapter 3 investigates helium isotope systematics along the western SWIR as well as near a global array of hotspots. The first part of this study reports uniformly low 3HetHe ratios of 6.3-7.3 RlRa along the SWIR from 9°-24°E, compared to values of 8±1 Ra for normal mid-ocean ridge basalt. The favored explanation for these low values is addition of (U+ Th) into the mantle source by crustal and/or lithospheric recycling. Although high HetHe values have been observed along the SWIR near Bouvet Island to the west, there is no evidence for elevated 3HetHe ratios along this section of the SWIR. The second part of Chapter 3 investigates the relationship between 3HetHe ratios and geophysical indicators of plume robustness for nine hotspots. A close correlation between a plume's flux and maximum 3HetHe ratio suggests a link between plume upwelling strength and origination in the deep, relatively undegassed mantle. Chapter 4 studies 3D mantle flow and temperature patterns beneath oceanic ridge-ridge-ridge triple junctions (TJs). In non-hotspot- affected TJs with geometry similar to the Rodrigues TJ, temperature and upwelling velocity along the slowest-spreading of the three ridges are predicted to increase within a few hundred kilometers of the TJ, to approach those of the fastest-spreading ridge. Along the slowest-spreading branch in hotspot-affected TJs such as the Azores, a strong component of along-axis flow directed away from the TJ is predicted to advect a hotspot thermal anomaly away from its deep-seated source.

Description: Funding was provided by the National Science Foundation through grants OCE- 9811924 and OCE-9907630, and by a National Defense Science and Engineering Graduate Fellowship.

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

Description: A new, global inversion is used to estimate the large scale oceanic circulation based on the World Ocean Circulation Experiment and Java Australia Dynamic Experiment hydrographic data. A linear inverse "box" model is used to combine consistently the transoceanic sections. The circulation is geostrophic with an Ekman layer at the surface and oceanic layers defined by neutral surfaces. Near-conservation of mass, salt and top-to-bottom silica is required and, in addition, heat and the phosphate-oxygen combination (170[P04]+[02]) are conserved in layers that are not in contact with the surface. A globally-consistent solution is obtained for a depth-independent adjustment to the thermal wind field, freshwater flux divergenees, the Ekman transport, and the advective and diffusive dianeutral fluxes between layers. A detailed error budget permits calculation of statistical uncertainties, taking into account both the non-resolved part of the solution and the systematic errors due to the temporal oceanic variability. The estimated water mass transports during the WOCE period (1985-1996) are generally similar to previous published estimates. However, important differences are found. In particular, the inflow of bottom waters into the Pacific Ocean is smaller than in most previous estimates. Utilization of property anomaly conservation constraints allows the estimation of significant dianeutral diffusivities in deep layers, with a global average of 3 ± lcm2s- 1 north of 30°S. Dianeutral transfers indicate that about 20 Sv of bottom water is formed in the Southern Ocean. Significant ocean-atmosphere heat fluxes are found, with a global heating of 2.3 ± 0.4PW in the tropical band and a corresponding cooling at high latitudes. The signature of a large-scale average export production is found for nutrients in several temperate regions. Despite the large uncertainties, the production magnitudes are consistent with independent measurements from sediment traps and isotopic data. Net nutrient sources or sinks are found in several regions, suggesting either transport of dissolved organic matter or a seasonal alias. Oxygen indicates large exchanges with the atmosphere, with intake at high latitudes and outgassing/remineralization at low latitudes.

Description: This work was supported in part by the Jet Propulsion Laboratory/CALTECH (contract #958125), and by gifts from Ford, General Motors, and Daimler-Chrysler to MIT's Climate Modelling Initiative.

Description: Author Posting. © American Geophysical Union 2003. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 4 (2003): 8515, doi:10.1029/2003GC000609.

Description: Complete multibeam bathymetric coverage of the western Galápagos Spreading Center (GSC) between 90.5°W and 98°W reveals the fine-scale morphology, segmentation and influence of the Galápagos hot spot on this intermediate spreading ridge. The western GSC comprises three morphologically defined provinces: A Western Province, located farthest from the Galápagos hot spot west of 95°30′W, is characterized by an axial deep, rift valley morphology with individual, overlapping, E-W striking segments separated by non-transform offsets; A Middle Province, between the propagating rift tips at 93°15′W and 95°30′W, with transitional axial morphology strikes ∼276°; An Eastern Province, closest to the Galápagos hot spot between the ∼90°50′W Galápagos Transform and 93°15′W, with an axial high morphology generally less than 1800 m deep, strikes ∼280°. At a finer scale, the axial region consists of 32 individual segments defined on the basis of smaller, mainly 〈2 km, offsets. These offsets mainly step left in the Western and Middle Provinces, and right in the Eastern Province. Glass compositions indicate that the GSC is segmented magmatically into 8 broad regions, with Mg # generally decreasing to the west within each region. Striking differences in bathymetric and lava fractionation patterns between the propagating rifts with tips at 93°15′W and 95°30′W reflect lower overall magma supply and larger offset distance at the latter. The structure of the Eastern Province is complicated by the intersection of a series of volcanic lineaments that appear to radiate away from a point located on the northern edge of the Galápagos platform, close to the southern limit of the Galápagos Fracture Zone. Where these lineaments intersect the GSC, the ridge axis is displaced to the south through a series of overlapping spreading centers (OSCs); abandoned OSC limbs lie even farther south. We propose that southward displacement of the axis is promoted during intermittent times of increased plume activity, when lithospheric zones of weakness become volcanically active. Following cessation of the increased plume activity, the axis straightens by decapitating southernmost OSC limbs during short-lived propagation events. This process contributes to the number of right stepping offsets in the Eastern Province.

Description: This work was supported by NSF grants OCE98- 18632 to the University of Hawai’i and OCE98-19117 to the Woods Hole Oceanographic Institution; support was provided to M. B. by a CIW/DTM Postdoctoral Fellowship

Description: Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 9 (2008): Q04015, doi:10.1029/2007GC001611.

Description: Near-bottom magnetic data collected along the crest of the East Pacific Rise between 9°55′ and 9°25′N identify the Central Anomaly Magnetization High (CAMH), a geomagnetic anomaly modulated by crustal accretionary processes over timescales of ∼104 years. A significant decrease in CAMH amplitude is observed along-axis from north to south, with the steepest gradient between 9°42′ and 9°36′N. The source of this variation is neither a systematic change in geochemistry nor varying paleointensity at the time of lava eruption. Instead, magnetic moment models show that it can be accounted for by an observed ∼50% decrease in seismic Layer 2A thickness along-axis. Layer 2A is assumed to be the extrusive volcanic layer, and we propose that this composes most of the magnetic source layer along the ridge axis. The 9°37′N overlapping spreading center (OSC) is located at the southern end of the steep CAMH gradient, and the 9°42′–9°36′N ridge segment is interpreted to be a transition zone in crustal accretion processes, with robust magmatism north of 9°42′N and relatively low magmatism at present south of 9°36′N. The 9°37′N OSC is also the only bathymetric discontinuity associated with a shift in the CAMH peak, which deviates ∼0.7 km to the west of the axial summit trough, indicating southward migration of the OSC. CAMH boundaries (defined from the maximum gradients) lie within or overlie the neovolcanic zone (NVZ) boundaries throughout our survey area, implying a systematic relationship between recent volcanic activity and CAMH source. Maximum flow distances and minimum lava dip angles are inferred on the basis of the lateral distance between the NVZ and CAMH boundaries. Lava dip angles average ∼14° toward the ridge axis, which agrees well with previous observations and offers a new method for estimating lava dip angles along fast spreading ridges where volcanic sequences are not exposed.

Description: The research project was funded by National Science Foundation under grants OCE-9819261 and OCE- 0096468.

Description: Author Posting. © American Geophysical Union, 2007. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 8 (2007): Q08013, doi:10.1029/2007GC001652.

Description: We report first evidence for hydrothermal activity from the southern Knipovich Ridge, an ultra-slow spreading ridge segment in the Norwegian-Greenland Sea. Evidence comes from optical backscatter anomalies collected during a systematic side-scan sonar survey of the ridge axis, augmented by the identification of biogeochemical tracers in the overlying water column that are diagnostic of hydrothermal plume discharge (Mn, CH4, ATP). Analysis of coregistered geologic and oceanographic data reveals that the signals we have identified are consistent with a single high-temperature hydrothermal source, located distant from any of the axial volcanic centers that define second-order segmentation along this oblique ridge system. Rather, our data indicate a hydrothermal source associated with highly tectonized seafloor that may be indicative of serpentinizing ultramafic outcrops. Consistent with this hypothesis, the hydrothermal plume signals we have detected exhibit a high methane to manganese ratio of 2–3:1. This is higher than that typical of volcanically hosted vent sites and provides further evidence that the source of the plume signals reported here is most probably a high-temperature hydrothermal field that experiences some ultramafic influence (compare to Rainbow and Logachev sites, Mid-Atlantic Ridge). While such sites have previously been invoked to be common on the SW Indian Ridge, this may be the first such site to be located along the Arctic ultra-slow spreading ridge system.

Description: Connelly and German were funded by NERC grant NER/B/S/ 2000/00755, NERC Core Strategic Funding at NOC, and the ChEss project of the Census of Marine Life.

Description: Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 9 (2008): Q08O10, doi:10.1029/2008GC001965.

Description: We use 2-D numerical models to explore the thermal and mechanical effects of magma intrusion on fault initiation and growth at slow and intermediate spreading ridges. Magma intrusion is simulated by widening a vertical column of model elements located within the lithosphere at a rate equal to a fraction, M, of the total spreading rate (i.e., M = 1 for fully magmatic spreading). Heat is added in proportion to the rate of intrusion to simulate the thermal effects of magma crystallization and the injection of hot magma into the crust. We examine a range of intrusion rates and axial thermal structures by varying M, spreading rate, and the efficiency of crustal cooling by conduction and hydrothermal circulation. Fault development proceeds in a sequential manner, with deformation focused on a single active normal fault whose location alternates between the two sides of the ridge axis. Fault spacing and heave are primarily sensitive to M and secondarily sensitive to axial lithosphere thickness and the rate that the lithosphere thickens with distance from the axis. Contrary to what is often cited in the literature, but consistent with prior results of mechanical modeling, we find that thicker axial lithosphere tends to reduce fault spacing and heave. In addition, fault spacing and heave are predicted to increase with decreasing rates of off-axis lithospheric thickening. The combination of low M, particularly when M approaches 0.5, as well as a reduced rate of off-axis lithospheric thickening produces long-lived, large-offset faults, similar to oceanic core complexes. Such long-lived faults produce a highly asymmetric axial thermal structure, with thinner lithosphere on the side with the active fault. This across-axis variation in thermal structure may tend to stabilize the active fault for longer periods of time and could concentrate hydrothermal circulation in the footwall of oceanic core complexes.

Description: Funding for this research was provided by NSF grants OCE-0327018 (M.D.B.), OCE-0548672 (M.D.B.), OCE- 0327051 (G.I.), and OCE-03-51234 (G.I.).

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 1985

Description: From October, 1982 to October, 1983 a current meter mooring reaching from the bottom into the thermocline was deployed for the first time in the Gulf Stream at 68°W. The temperatures, pressures, and velocities at the uppermost instrument indicate the Gulf Stream moved back and forth across the mooring site, so that the entire Stream was sampled in time; hence the data may be used to examine horizontal as well as vertical structure of the Stream. The two key points to the success of the analysis are: 1)the well-defined relationship between temperature and cross-stream distance in the thermocline, enabling the use of the former as a horizontal coordinate; and 2)a daily-changing definition of Gulf Stream flow direction based on the shear between the thermocline and 2000 m depth. Time-series of daily-rotated velocities may be used to calculate empirical orthogonal functions for the long- and cross-stream vertical structures, which are decoupled and are respectively baroclinic and barotropic. Using the inferred horizontal coordinate one can estimate mass, momentum and kinetic energy fluxes for four individual events when the entire Stream swept by the mooring. The results agree well with historical data. Bryden's (1980) method has been used to calculate vertical velocities from the temperature equation; the resulting time-series of w are visually coherent throughout the water column and their vertical amplitude structure is reminiscent of that for a two-layer system. The rms vertical velocities are large (0(.05 cm/s)), and these as well as other estimates have been used to explore the validity of the quasi-geostrophic approximation at the mooring site. The Rossby number for the thermocline flow is about 0.3, and for the deep flow is ≤ 0.1. The entire data set may also be used to construct a horizontal and vertical profile of velocity in the Gulf Stream, from which a cross-section of the mean potential vorticity can be produced. The latter shares many common feature with cross-sections from past work for a nearby site, as well as analogous data from a three-layer numerical model, thus suggesting that they are robust features of Gulf Stream-like currents. These features are, in particular, a strong jump from low to high values crossing the Stream from south to north; and a change in the sign of the potential vorticity gradient on isothermal surfaces for T 〉 12°C. To complement the analysis of the observational data, a set of diagnostic calculations has been performed on an eddy-resolving qeneral circulation model, to provide a complete picture of the kinetic energy budgets of the free jet and its environs. It is found that the downstream convergence of kinetic energy in the decelerating jet is balanced primarily by an ageostrophic flow against the pressure gradient, which in turn implies some conversion of kinetic to available potential energy in the region. Energetic analysis of the observations as well as the numerical data suggests barotropic and baroclinic instabilities may be equally important to the kinetic energy budgets in the Stream. Because there is but one mooring, the dynamics governing the fluctuations remain elusive. Nonetheless, a kinematic framework is proposed, which is consistent with the data and accounts for a variety of unusual features that arise in the original analysis (for example, distinct asymmetries in the four Gulf Stream crossings, and the rather large vertical velocities). It is sugqested that the data we are now capable of collecting is proffering fundamentally new attributes of the Gulf Stream, which must be included and accounted for in future theoretical work.

Description: Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 113 (2008): G02026, doi:10.1029/2007JG000470.

Description: Permafrost is a defining characteristic of the Arctic environment. However, climate warming is thawing permafrost in many areas leading to failures in soil structure called thermokarst. An extensive survey of a 600 km2 area in and around the Toolik Lake Natural Research Area (TLNRA) revealed at least 34 thermokarst features, two thirds of which were new since ∼1980 when a high resolution aerial survey of the area was done. Most of these thermokarst features were associated with headwater streams or lakes. We have measured significantly increased sediment and nutrient loading from thermokarst features to streams in two well-studied locations near the TLNRA. One small thermokarst gully that formed in 2003 on the Toolik River in a 0.9 km2 subcatchment delivered more sediment to the river than is normally delivered in 18 years from 132 km2 in the adjacent upper Kuparuk River basin (a long-term monitoring reference site). Ammonium, nitrate, and phosphate concentrations downstream from a thermokarst feature on Imnavait Creek increased significantly compared to upstream reference concentrations and the increased concentrations persisted over the period of sampling (1999–2005). The downstream concentrations were similar to those we have used in a long-term experimental manipulation of the Kuparuk River and that have significantly altered the structure and function of that river. A subsampling of other thermokarst features from the extensive regional survey showed that concentrations of ammonium, nitrate, and phosphate were always higher downstream of the thermokarst features. Our previous research has shown that even minor increases in nutrient loading stimulate primary and secondary production. However, increased sediment loading could interfere with benthic communities and change the responses to increased nutrient delivery. Although the terrestrial area impacted by thermokarsts is limited, the aquatic habitat altered by these failures can be extensive. If warming in the Arctic foothills accelerates thermokarst formation, there may be substantial and wide-spread impacts on arctic stream ecosystems that are currently poorly understood.

Description: The results presented in this report are based upon work supported by the U.S. National Science Foundation under grants to the Arctic Hyporheic project (OPP- 0327440) and the Arctic Long-Term Ecological Research Program (DEB- 9810222).

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

Description: Data from the COBOLT experiment, which investigated the first 12 km off Long Island's south shore, are analyzed and discussed. Moored current meter records indicate that the nearshore flow field is strongly polarized in the alongshore direction and its fluctuations are well correlated with local meteorological forcing. Complex empirical orthogonal function analysis suggests that subtidal velocity fluctuations are barotropic in nature and are strongly influenced by bottom friction. Wind-related inertial currents were observed within the coastal boundary layer (CBL) under favorable meteorological and hydrographical conditions. The magnitude of these oscillations increases with distance from shore, and they display a very clear 180° phase difference between surface and bottom layers. Nearshore inertial oscillations of both velocity and salinity records appear to lead those further seaward, suggesting local generation and subsequent radiation away from the coast. The response of the coastal zone to impulsive wind forcing is discussed using simple slab and two-layer models, and the behavior of the nearshore current field examined. The major features of the observed inertial motions are in good qualitative agreement with model predictions. It is found that, in a homogeneous domain, the coastal boundary condition effectively prohibits inertial currents over the entire coastal zone. In the presence of stratification the offshore extent of this prohibition is greatly reduced and significant inertial currents may occur within one or two internal deformation radii of the coast. The "coastal effect", in the form of surface and interfacial waves which propagate away from the coast, modifies the "pure" inertial response as it would exist far from shore. The kinematics of this process is such that a 180° phase difference between currents in the two layers is characteristic of the entire coastal zone even before the internal wave has had time to traverse the CBL. It is also suggested that, for positions seaward of several internal deformation radii, interference between the surface and internal components of the coastal response will cause maximum inertial amplitudes to occur for t 〉 x/c2, where c2 is the phase speed of the internal disturbance. The hydrographic structure of the CBL is observed to undergo frequent homogenization. These events are related to both advective and mixing processes. Horizontal and vertical exchange coefficients are estimated from the data, and subsequently used in a diffusive model which accurately reproduces the observed mean density distribution in the nearshore zone. Dynamic balance calculations are performed which indicate that the subtidal cross-shore momentum balance is very nearly geostrophic. The calculations also suggest that the longshore balance may be reasonably represented by a steady, linear equation of motion which includes surface and bottom stresses. Evidence is presented which shows that variations in the longshore wind-stress component are primarily responsible for the energetic fluctuations in the sea surface slope along Long Island. Depth-averaged velocities characteristically show net offshore transport in the study area, and often display dramatic longshore current reversals with distance from shore. These observations are interpreted in terms of a steady circulation model which includes realistic nearshore topography. Model results suggest that longshore current reversals within the CBL may be limited to the eastern end of Long Island, and that this unusual flow pattern is a consequence of flow convergence related to the presence of Long Island Sound.

Description: This work was supported by the Department of Energy through contract no. DE-AC02-EVI0005 entitled Coastal-Shelf Transport and Diffusion.

Description: Author Posting. © American Geophysical Union, 2010. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 115 (2010): C10018, doi:10.1029/2009JC005660.

Description: Variations in the Arctic central Canada Basin mixed layer properties are documented based on a subset of nearly 6500 temperature and salinity profiles acquired by Ice-Tethered Profilers during the period summer 2004 to summer 2009 and analyzed in conjunction with sea ice observations from ice mass balance buoys and atmosphere-ocean heat flux estimates. The July–August mean mixed layer depth based on the Ice-Tethered Profiler data averaged 16 m (an overestimate due to the Ice-Tethered Profiler sampling characteristics and present analysis procedures), while the average winter mixed layer depth was only 24 m, with individual observations rarely exceeding 40 m. Guidance interpreting the observations is provided by a 1-D ocean mixed layer model. The analysis focuses attention on the very strong density stratification at the base of the mixed layer in the Canada Basin that greatly impedes surface layer deepening and thus limits the flux of deep ocean heat to the surface that could influence sea ice growth/decay. The observations additionally suggest that efficient lateral mixed layer restratification processes are active in the Arctic, also impeding mixed layer deepening.

Description: Support for the ITP program and this study was provided by the U. S. National Science Foundation and the Woods Hole Oceanographic Institution. Support for the IMB program came from the National Science Foundation and the National Oceanographic and Atmospheric Administration.

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 2010

Description: Mid-ocean ridge basalts (MORBs) exhibit a wide range of CO2 concentrations, reflecting saturation to supersaturation (and rarely undersaturation) relative to their emplacement depths. In this study, we explore the mechanisms of CO2 degassing and the implications this has for estimating the advance rates and durations of seafloor eruptions. We present dissolved volatile concentrations (mainly of CO2 and H2O) and vesicle size distributions (VSDs) for a unique suite of MORB glasses collected at the East Pacific Rise, ~9° 50′ N. These MORB glasses were collected at ~200 m intervals along an across-axis track over a single flow pathway within the recently emplaced 2005-06 eruption boundaries; systematic sample collection provides one of the first opportunities to characterize intra-flow geochemical and physical evolution during a single eruption at a fast-spreading ridge. Compared to measurements of MORB volatiles globally, dissolved H2O concentrations are relatively uniform (0.10 - 0.16 weight percent), whereas dissolved CO2 contents exhibit a range of concentrations (154 - 278 ppm) and decrease with distance from the EPR axis (i.e., eruptive vent). Ion microprobe analyses of dissolved volatiles within the MORB glasses suggest that the magma erupted supersaturated (pressure equilibrium with 920 - 1224 mbsf) and in near-equilibrium with the melt lens of the axial magma chamber (~1250 - 1500 mbsf), and degassed to near equilibrium (299 - 447 mbsf) with seafloor depths over the length of the flow. The decrease in CO2 concentrations spans nearly the full range of dissolved CO2 contents observed at the EPR and shows that the varying degrees of volatile saturation that have been observed in other MORB sample suites may be explained by degassing during emplacement. Vesicularity (0.1 - 1.2%) increases with decreasing dissolved CO2 concentrations. We use vesicle size distributions (VSDs)—vesicle sizes and number densities—to quantify the physical evolution of the CO2 degassing process. VSDs suggest that diffusion of CO2 into preexisting vesicles, and not nucleation of new vesicles, is the dominant mechanism of increasing CO2 in the vapor phase. We also use VSDs, along with estimates of vesicle growth rates, to constrain emplacement time of the 2005-06 eruption to 〈~24 hours and to resolve variations in advance rate with downflow distance.

Description: Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 9 (2008): Q08001, doi:10.1029/2007GC001922.

Description: Multichannel seismic and bathymetric data from the Juan de Fuca Ridge (JDFR) provide constraints on axial and ridge flank structure for the past 4–8 Ma within three spreading corridors crossing Cleft, Northern Symmetric, and Endeavour segments. Along-axis data reveal south-to-north gradients in seafloor relief and presence and depth of the crustal magma lens, which indicate a warmer axial regime to the south, both on a regional scale and within individual segments. For young crust, cross-axis lines reveal differences between segments in Moho two-way traveltimes of 200–300 ms which indicate 0.5–1 km thicker crust at Endeavour and Cleft compared to Northern Symmetric. Moho traveltime anomalies extend beyond the 5–15 km wide axial high and coincide with distinct plateaus, 32 and 40 km wide and 200–400 m high, found at both segments. On older crust, Moho traveltimes are similar for all three segments (∼2100 ± 100 ms), indicating little difference in average crustal production prior to ∼0.6 and 0.7 Ma. The presence of broad axis-centered bathymetric plateau with thickened crust at Cleft and Endeavour segments is attributed to recent initiation of ridge axis-centered melt anomalies associated with the Cobb hot spot and the Heckle melt anomaly. Increased melt supply at Cleft segment upon initiation of Axial Volcano and southward propagation of Endeavour segment during the Brunhes point to rapid southward directed along-axis channeling of melt anomalies linked to these hot spots. Preferential southward flow of the Cobb and Heckle melt anomalies and the regional-scale south-to-north gradients in ridge structure along the JDFR may reflect influence of the northwesterly absolute motion of the ridge axis on subaxial melt distribution.

Description: This work was supported by U.S. National Science Foundation grants OCE00-02488 to S.M.C., OCE06-48303 to S.M.C. and M.R.N., OCE-0648923 to J.P.C., and OCE00-02600 to G.M.K. and A.J.H.

Description: Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 113 (2008): C02018, doi:10.1029/2007JC004429.

Description: Radioisotope evaluation of a cold-core, anticyclonic eddy surveyed in September 2004 on the Chukchi Sea continental slope was used to determine its age since formation over the shelf environment. Because the eddy can be shown to have been generated near the shelf break, initial conditions for several age-dependent tracers could be relatively well constrained. A combination of 228Ra/226Ra, excess 224Ra, and 228Th/228Ra suggested an age on the order of months. This age is consistent with the presence of elevated concentrations of nutrients, organic carbon, suspended particles, and shelf-derived neritic zooplankton within the eddy compared to ambient offshore water in the Canada Basin but comparable to values measured in the Chukchi shelf and shelf-break environment. Hence this feature, at the edge of the deep basin, was poised to deliver biogeochemically significant shelf material to the central Arctic Ocean.

Description: This work was supported by National Science Foundation Polar Programs grants OPP-662690 and OPP-66040N to the University of Miami (DK), and Office of Naval Research grant N00014-02-1-0317 (RP).

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 2010

Description: In a stratified rotating fluid, frictionally driven circulations couple with the buoyancy field over sloping topography. Analytical and numerical methods are used to quantify the impact of this coupling on the vertical circulation, spindown of geostrophic flows, and the formation of a shelfbreak jet. Over a stratified slope, linear spindown of a geostrophic along-isobath flow induces cross-isobath Ekman flows. Ekman advection of buoyancy weakens the vertical circulation and slows spindown. Upslope (downslope) Ekman flows tend to inject (remove) potential vorticity into (from) the ocean. Momentum advection and nonlinear buoyancy advection are examined in setting asymmetries in the vertical circulation and the vertical relative vorticity field. During nonlinear homogeneous spindown over a flat bottom, momentum advection weakens Ekman pumping and strengthens Ekman suction, while cyclonic vorticity decays faster than anticyclonic vorticity. During nonlinear stratified spindown over a slope, nonlinear advection of buoyancy enhances the asymmetry in Ekman pumping and suction, whereas anticyclonic vorticity can decay faster than cyclonic vorticity outside of the boundary layers. During the adjustment of a spatially uniform geostrophic current over a shelfbreak, coupling between the Ekman flow and the buoyancy field generates Ekman pumping near the shelfbreak, which leads to the formation of a jet. Scalings are presented for the upwelling strength, the length scale over which it occurs, and the timescale for jet formation. The results are applied to the Middle Atlantic Bight shelfbreak.

Description: Funding for my research and education was provided by MIT EAPS, the WHOI Academic Programs O ce and the MIT Presidential Fellowship. Financial assistance from the Houghton Fund is also 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 June 2010

Description: The Labrador Sea, as one of a few places of deep water formation, plays an important role in the Meridional Overturning Circulation. While the interior of the Labrador Sea, where the deepest convection takes place, is known to experience variability on time scales ranging from days to decades, little is known about the variability of the other components of the Labrador Sea circulation - the boundary current system and the eddies that connect it with the interior. Using various types of in situ data combined with the surface flux and satellite altimetry data products, I studied the variability of both the boundary current system and the eddies on different time scales as well as their influence on the post-convective restratification of the Labrador Sea interior. The analysis presented in the thesis supports the result of the previous theoretical studies that argue that lateral fluxes, driven by the boundary current/interior gradients, play an important role in the post-convective restratification of the Labrador Sea. I found that both components of the boundary current, the surface West Greenland Current and the subsurface Irminger Current, have a strong seasonal cycle. In the spring both the West Greenland and Irminger Currents are colder and fresher than in the fall. However, the West Greenland Current is faster and thicker in the spring while the Irminger Current is the fastest and thickest in the fall. My analysis suggests that the observed seasonal changes in the velocity are primarily due to the baroclinic component of the current while the barotropic component remains nearly unchanged. The Subpolar Gyre, and the Labrador Sea in particular, have experienced a decline in the circulation accompanied by the warming of the water column over the last decades. I found that a similar trend is seen in the West Greenland Current system which slowed down from 1992 to 2004, primarily due to a decrease in the barotropic flow. At the same time, the subsurface Irminger Current has become warmer, saltier, and lighter, something that is also reflected in the properties of the eddies. Two years exhibited pronounced anomalies: in 1997 and 2003 the velocity, temperature and salinity of the Irminger Current abruptly increase with respect to the overall trend. Finally, I discuss the impacts of the boundary current changes on the lateral fluxes that are responsible for the restratification of the Labrador Sea and the properties of the interior.

Description: The financial support for my research came from the Academic Programs Office and from the NSF grants OCE-0424492 and OCE-0137023.

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 2010

Description: The variability in the DWBC, its connection to the forcing in the northern North Atlantic and interaction with the Gulf Stream were explored from a combination of remote sensing and in-situ measurements in the western North Atlantic. Using satellite altimetry and Sea Surface Temperature (SST) we found evidence of the relation between changes in the Gulf Stream path and the variability in the temperature and velocity fields in the Slope Water. This relation was such that southward shifts of the main axis of the Gulf Stream were preceded by cold temperature anomalies and intensification of the southwestward flow. The analysis of 5.5 years of moored CTD and horizontal velocity data in the DWBC at 69°W recorded during the period 2002-2008, showed that the variability along the DWBC is linked to changes in the dense water formation regions. The evolution of potential vorticity (PV) at the mooring site, characterized by a transition from deep to upper Labrador Sea Water (LSW), was similar to that observed in the Labrador Sea 6 to 9 years earlier, and imply spreading rates for the LSW that varied over time from 1.5 to 2.5cm/s. The time dependence of the spreading rates was in good agreement with changes in the strength of the DWBC at the mooring site. The evolution of the DWBC transport was explored in more detail from a 5-element moored array, also at 69°W. The results, for the period of 2004-2008, were consistent with the single mooring analysis. The variability measured from the array showed that upper, intermediate and deep water mass layers expand and contract at each other’s expense, leading to alternating positive and negative PV anomalies at the upper-LSW, deep-LSW and Overflow Water (OW). Larger DWBC transports were associated with enhanced presence of recently ventilated upper-LSW and OW, rather than deep-LSW. The relative contribution of the different water masses to the observed circulation was investigated by inverting individual PV anomalies isolated from the observations. We found that changes in the depth-integrated circulation were mostly driven by changes in the OW.

Description: During the first two years, my time here was funded through a La Caixa Foundation fellowship as well as an MIT presidential fellowship. Financial support for my thesis work was provided by National Science Foundation grants OCE-0241354 and OCE-0726720. Support for traveling to meetings was also provided by the MIT Houghton Fund and the WHOI Academic Programs Office.

Description: Author Posting. © American Geophysical Union, 2007. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 112 (2007): G04S54, doi:10.1029/2006JG000353.

Description: Dramatic changes have been observed in the Arctic over the last century. Many of these involve the storage and cycling of fresh water. On land, precipitation and river discharge, lake abundance and size, glacier area and volume, soil moisture, and a variety of permafrost characteristics have changed. In the ocean, sea ice thickness and areal coverage have decreased and water mass circulation patterns have shifted, changing freshwater pathways and sea ice cover dynamics. Precipitation onto the ocean surface has also changed. Such changes are expected to continue, and perhaps accelerate, in the coming century, enhanced by complex feedbacks between the oceanic, atmospheric, and terrestrial freshwater systems. Change to the arctic freshwater system heralds changes for our global physical and ecological environment as well as human activities in the Arctic. In this paper we review observed changes in the arctic freshwater system over the last century in terrestrial, atmospheric, and oceanic systems.

Description: The authors gratefully acknowledge the National Science Foundation (NSF) for funding this synthesis work. This paper is principally the work of authors funded under the NSF-funded Freshwater Integration (FWI) study.

Description: Author Posting. © American Geophysical Union, 2004. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 109 (2004): C03051, doi:10.1029/2003JC001940.

Description: Arctic Ocean model simulations have revealed that the Arctic Ocean has a basin-wide oscillation with cyclonic and anticyclonic circulation anomalies (Arctic Ocean Oscillation (AOO)) that has a prominent decadal variability [Proshutinsky and Johnson, 1997]. This study explores how the simulated AOO affects the Arctic Ocean stratification and its relationship to the sea ice cover variations. The simulation uses the Princeton Ocean Model coupled to sea ice [Häkkinen and Mellor, 1992; Häkkinen, 1999]. The surface forcing is based on National Centers for Environmental Prediction/National Center for Atmospheric Research Reanalysis and its climatology, of which the latter is used to force the model spin-up phase. Our focus is to investigate the competition between ocean dynamics and ice formation/melt on the Arctic basin-wide freshwater balance. We find that changes in the Atlantic water inflow can explain almost all of the simulated freshwater anomalies in the main Arctic basin. The Atlantic water inflow anomalies are an essential part of AOO, which is the wind driven barotropic response to the Arctic Oscillation (AO). The baroclinic response to AO, such as Ekman pumping in the Beaufort Gyre, and ice melt/freeze anomalies in response to AO are less significant considering the whole Arctic freshwater balance.

Description: We gratefully acknowledge the support from National Science Foundation under Grant No OPP-0230184 (AP) and from NASA Headquarters (SH).

Description: Author Posting. © American Geophysical Union, 2004. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 109 (2004): C03042, doi:10.1029/2003JC002007.

Description: Sea level is a natural integral indicator of climate variability. It reflects changes in practically all dynamic and thermodynamic processes of terrestrial, oceanic, atmospheric, and cryospheric origin. The use of estimates of sea level rise as an indicator of climate change therefore incurs the difficulty that the inferred sea level change is the net result of many individual effects of environmental forcing. Since some of these effects may offset others, the cause of the sea level response to climate change remains somewhat uncertain. This paper is focused on an attempt to provide first-order answers to two questions, namely, what is the rate of sea level change in the Arctic Ocean, and furthermore, what is the role of each of the individual contributing factors to observed Arctic Ocean sea level change? In seeking answers to these questions we have discovered that during the period 1954–1989 the observed sea level over the Russian sector of the Arctic Ocean is rising at a rate of approximately 0.123 cm yr−1 and that after correction for the process of glacial isostatic adjustment this rate is approximately 0.185 cm yr−1. There are two major causes of this rise. The first is associated with the steric effect of ocean expansion. This effect is responsible for a contribution of approximately 0.064 cm yr−1 to the total rate of rise (35%). The second most important factor is related to the ongoing decrease of sea level atmospheric pressure over the Arctic Ocean, which contributes 0.056 cm yr−1, or approximately 30% of the net positive sea level trend. A third contribution to the sea level increase involves wind action and the increase of cyclonic winds over the Arctic Ocean, which leads to sea level rise at a rate of 0.018 cm yr−1 or approximately 10% of the total. The combined effect of the sea level rise due to an increase of river runoff and the sea level fall due to a negative trend in precipitation minus evaporation over the ocean is close to 0. For the Russian sector of the Arctic Ocean it therefore appears that approximately 25% of the trend of 0.185 cm yr−1, a contribution of 0.048 cm yr−1, may be due to the effect of increasing Arctic Ocean mass.

Description: This material is based upon work supported by the National Science Foundation under grant 0136432.

Description: Author Posting. © American Geophysical Union, 2004. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 109 (2004): B10310, doi:10.1029/2004JB003066.

Description: Multichannel seismic reflection data are used to infer crustal accretion processes along the intermediate spreading Galapagos Spreading Center. East of 92.5°W, we image a magma lens beneath the ridge axis that is relatively shallow (1.0–2.5 km below the seafloor) and narrow (∼0.5–1.5 km, cross-axis width). We also image a thin seismic layer 2A (0.24–0.42 km) that thickens away from the ridge axis by as much as 150%. West of 92.7°W, the magma lens is deeper (2.5–4.5 km) and wider (0.7–2.4 km), and layer 2A is thicker (0.36–0.66 km) and thickens off axis by 〈40%. The positive correlation between layer 2A thickness and magma lens depth supports the interpretation of layer 2A as the extrusive volcanic layer with thickness controlled by the pressure on the magma lens and its ability to push magma to the surface. Our findings also suggest that narrower magma lenses focus diking close the ridge axis such that lava flowing away from the ridge axis will blanket older flows and thicken the extrusive crust off axis. Flow of lava away from the ridge axis is probably promoted by the slope of the axial bathymetric high, which is largest east of 92.5°W. West of ∼94°W the “transitional” axial morphology lacks a prominent bathymetric high and layer 2A no longer thickens off axis. We detect no magma lens west of 94.7°W where a small axial valley appears. The above changes can be linked to the westward decrease in the magma and heat flux associated with the fading influence of the Galapagos hot spot on the Galapagos Spreading Center.

Description: This project was funded by NSF-OCE- 0002189.