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 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 Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution August 1987

Description: Several problems connected by the theme of thermal forcing are addressed herein. The main topic is the stratification and flow field resulting from imposing a specified heat flux on a fluid that is otherwise confined to a rigid insulating basin. In addition to the traditional eddy viscosity and diffusivity, turbulent processes are also included by a convective overturning adjustment at locations where the local density field is unstable. Two classes of problems are treated. The first is the large scale meridional pattern of a fluid in an annulus. The detailed treatment is carried out in two steps. In the beginning (chapter 2) it is assumed that the fluid is very diffusive, hence, to first approximation no flow field is present. It is found that the convective overturning adjustment changes the character of the stratification in all the regions that are cooled from the top, resulting in a temperature field that is nearly depth independent in the northernmost latitudes. The response to a seasonal cycle in the forcing, and the differences between averaging the results from the end of each season compared to driving the fluid by a mean forcing are analyzed. In particular, the resulting sea surface temperature is warmer in the former procedure. This observation is important in models where the heat flux is sensitive to the gradient of air to sea surface temperatures. The analysis of the problem continues in chapter 5 where the contribution of the flow field is included in the same configuration. The dimensionless parameter controlling the circulation is now the Rayleigh number, which is a measure of the relative importance of gravitational and viscous forces. The effects of the convective overturning adjustment is investigated at different Rayleigh numbers. It is shown that not only is the stratification now always stable, but also that the vigorous vertical mixing reduces the effective Rayleigh number; thereby the flow field is more moderate, the thermocline deepens, and the horizontal surface temperature gradients are weaker. The interior of the fluid is colder compared to cases without convective overturning, and, because the amount of heat in the system is assumed to be fixed, the surface temperature is warmer. The fluid is not only forced by a mean heat flux, or a seasonally varying one, but its behavior under permanent winter and summer conditions is also investigated. A steady state for the experiments where the net heat flux does not vanish is defined as that state where the flow field and temperature structure are not changing with time except for an almost uniform temperature decrease or increase everywhere. It is found that when winter conditions prevail the circulation is very strong, while it is rather weak for continuous summer forcing. In contrast to those results, if a yearly cycle is imposed, the circulation tends to reach a minimum in the winter time and a maximum in the summer. This suggests that, depending on the Rayleigh number, there is a phase leg of several months between the response of the ocean and the imposed forcing. Differences between the two averaging procedures mentioned before are also observed when the flow field is present, especially for large Rayleigh numbers. The circulation is found to be weaker and the sea surface temperature colder in the mean of the seasonal realizations compared to the steady state derived by the mean forcing. As an extension to the numerical results, an analytic model is presented in chapter 4 for a similar annular configuration. The assumed dynamics is a bit different, with a mixed layer on top of a potential vorticity conserving interior. It is demonstrated that the addition of the thermal wind balance to the conservation of potential vorticity in the axially symmetric problem leads to the result that typical fluid trajectories in the interior are straight lines pointing downward going north to south. The passage of information in the system is surprisingly in the opposite sense to the clockwise direction of the flow. A model for water mass formation by buoyancy loss in the absence of a flow field is introduced in chapter 3. The idea behind it is to use the turbulent mixing parameterization to generate chimney-like structures in open water, followed by along-isopycnal advection and diffusion. This model can be applied to many observations of mode water. In particular, in this work it is related to the chimneys observed by the MEDOC Group (1970), and the Levantine Intermediate Water in the Eastern Mediterranean Basin. An analytic prediction of the depth of the water mass is derived and depends on the forcing and initial stratification. It suggests that the depth of shallow mode water like the 18°C water or the Levantine Intermediate Water would not be very sensitive to reasonable changes in atmospheric forcing. Similar conclusions were also reached by Warren (1972) by assuming that the temperature in the thermocline decreases linearly with depth, and by approximating the energy balance in a water column by a Newtonian cooling law.

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

Description: Fits of an annual harmonic to depth-average along-shelf current time series longer than 200 days from 27 sites over the Middle Atlantic Bight (MAB) continental shelf have amplitudes of a few centimeters per second. These seasonal variations are forced by seasonal variations in the wind stress and the cross-shelf density gradient. The component of wind stress that drives the along-shelf flow over most of the MAB mid- and outer shelf is oriented northeast–southwest, perpendicular to the major axis of the seasonal variation in the wind stress. Consequently, there is not a significant seasonal variation in the wind-driven along-shelf flow, except over the southern MAB shelf and the inner shelf of New England where the wind stress components forcing the along-shelf flow are north–south and east–west, respectively. The seasonal variation in the residual along-shelf flow, after removing the wind-driven component, has an amplitude of a few centimeters per second with maximum southwestward flow in spring onshore of the 60-m isobath and autumn offshore of the 60-m isobath. The spring maximum onshore of the 60-m isobath is consistent with the maximum river discharges in spring enhancing cross-shelf salinity gradients. The autumn maximum offshore of the 60-m isobath and a steady phase increase with water depth offshore of Cape Cod are both consistent with the seasonal variation in the cross-shelf temperature gradient associated with the development and destruction of a near-bottom pool of cold water over the mid and outer shelf (“cold pool”) due to seasonal variations in surface heat flux and wind stress.

Description: This research was funded by the Ocean Sciences Division of the National Science Foundation under Grants OCE-820773, OCE-841292, and OCE- 848961.

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

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

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

Description: Author Posting. © American Meteorological Society, 2007. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 20 (2007): 2416-2433, doi:10.1175/JCLI4103.1.

Description: North Pacific decadal oceanic and atmospheric variability is examined from a 650-yr control integration of the Community Climate System Model version 2. The dominant pattern of winter sea surface temperature (SST) variability is similar to the observed “Pacific decadal oscillation,” with maximum amplitude along the Kuroshio Extension. SST anomalies in this region exhibit significant spectral peaks at approximately 16 and 40 yr. Lateral geostrophic heat flux divergence, caused by a meridional shift of the Kuroshio Extension forced by basin-scale wind stress curl anomalies 3–5 yr earlier, is responsible for the decadal SST variability; local surface heat flux and Ekman heat flux divergence act as a damping and positive feedback, respectively. A simple linear Rossby wave model is invoked to explicitly demonstrate the link between the wind stress curl forcing and decadal variability in the Kuroshio Extension. The Rossby wave model not only successfully reproduces the two decadal spectral peaks, but also illustrates that only the low-frequency (〉10-yr period) portion of the approximately white noise wind stress curl forcing is relevant. This model also demonstrates that the weak and insignificant decadal spectral peaks in the wind stress curl forcing are necessary for producing the corresponding strong and significant oceanic peaks in the Kuroshio Extension. The wind stress curl response to decadal SST anomalies in the Kuroshio Extension is similar in structure but opposite in sign and somewhat weaker than the wind stress curl forcing pattern. These results suggest that the simulated North Pacific decadal variability owes its existence to two-way ocean–atmosphere coupling.

Description: The first author gratefully acknowledges financial support from NOAA’s Office of Global Programs (grant to C. Deser) and the CCSM Project Office.

Description: Author Posting. © American Meteorological Society, 2007. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 20 (2007): 3785–3801, doi:10.1175/JCLI4234.1

Description: The influences of strong gradients in sea surface temperature on near-surface cross-front winds are explored in a series of idealized numerical modeling experiments. The atmospheric model is the Naval Research Laboratory Coupled Ocean–Atmosphere Mesoscale Prediction System (COAMPS) model, which is fully coupled to the Regional Ocean Modeling System (ROMS) ocean model. A series of idealized, two-dimensional model calculations is carried out in which the wind blows from the warm-to-cold side or the cold-to-warm side of an initially prescribed ocean front. The evolution of the near-surface winds, boundary layer, and thermal structure is described, and the balances in the momentum equation are diagnosed. The changes in surface winds across the front are consistent with previous models and observations, showing a strong positive correlation with the sea surface temperature and boundary layer thickness. The coupling arises mainly as a result of changes in the flux Richardson number across the front, and the strength of the coupling coefficient grows quadratically with the strength of the cross-front geostrophic wind. The acceleration of the winds over warm water results primarily from the rapid change in turbulent mixing and the resulting unbalanced Coriolis force in the vicinity of the front. Much of the loss/gain of momentum perpendicular to the front in the upper and lower boundary layer results from acceleration/deceleration of the flow parallel to the front via the Coriolis term. This mechanism is different from the previously suggested processes of downward mixing of momentum and adjustment to the horizontal pressure gradient, and is active for flows off the equator with sufficiently strong winds. Although the main focus of this work is on the midlatitude, strong wind regime, calculations at low latitudes and with weak winds show that the pressure gradient and turbulent mixing terms dominate the cross-front momentum budget, consistent with previous work.

Description: This work was supported by the Office of Naval Research Grant N00014-05-1-0300.

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

Description: Inverse methods are applied to historical hydrographic data to address two aspects of the general circulation of the Atlantic Ocean. The method allows conservation statements for mass and other properties, along with a variety of other constraints, to be combined in a dynamically consistent way to estimate the absolute velocity field and associated property transports. The method is first used to examine the exchange of mass and heat between the South Atlantic and the neighboring ocean basins. The Antarctic Circumpolar Current (ACC) carries a surplus of intermediate water into the South Atlantic through Drake Passage which is compensated by a surplus of deep and bottom water leaving the basin south of Africa. As a result, the ACC loses .25±.18x1015 W of heat in crossing the Atlantic. At 32°S the meridional flux of heat is .25±.19x1015 W equatorward, consistent in sign but smaller in magnitude than other recent estimates. This heat flux is carried primarily by a meridional overturning cell in which the export of 17 Sv of North Atlantic Deep Water (NADW) is balanced by an equatorward return flow equally split between the surface layers, and the intermediate and bottom water. No "leak" of warm Indian Ocean thermocline water is necessary to account for the equatorward heat flux across 32°S; in fact, a large transfer of warm water from the Indian Ocean to the Atlantic is found to be inconsistent with the present data set. Together these results demonstrate that the Atlantic as a whole acts to convert intermediate water to deep and bottom water, and thus that the global thermohaline cell associated with the formation and export of NADW is closed primarily by a "cold water path," in which deep water leaving the Atlantic ultimately returns as intermediate water entering the basin through Drake Passage. The second problem addressed concerns the circulation and property fluxes across 24°and 36°N in the subtropical North Atlantic. Conservation statements are considered for the nutrients as well as mass, and the nutrients are found to contribute significant information independent of temperature and salinity. Silicate is particularly effective in reducing the indeterminacy of circulation estimates based on mass conservation alone. In turn, the results demonstrate that accurate estimates of the chemical fluxes depend on relatively detailed knowledge of the circulation. The zonal-integral of the circulation consists of an overturning cell at both latitudes, with a net export of 19 Sv of NADW. This cell results in a poleward heat flux of 1.3±.2x1015 Wand an equatorward oxygen flux of 2900±180 kmol S-l across each latitude. The net flux of silicate is also equatorward: 138±38 kmol s-1 and 152±56 kmol s -1 across 36°and 24° N, respectively. However, in contrast to heat and oxygen, the overturning cell is not the only important mechanism responsible for the net silicate transport. A horizontal recirculation consisting of northward flow of silica-rich deep water in the eastern basin balanced by southward flow of low silica water in the western basin results in a significant silicate flux to the north. The net equatorward flux is thus smaller than indicated by the overturning cell alone. The net flux of nitrate across 36°N is n9±35 kmol 8- 1 to the north and is indistinguishable from zero at 24°N (-8±39 kmol 8-1 ), leading to a net divergence of nitrate between these two latitudes. Forcing the system to conserve nitrate leads to an unreasonable circulation. The dominant contribution to the nitrate flux at 36°N results from the correlation of strong northward flow and relatively high nitrate concentrations in the sub-surface waters of the Gulf Stream. The observed nitrate divergence between 24°and 36°N, and convergence north of 36°N, can be accounted for by a shallow cell in which the northward flow of inorganic nitrogen (nitrate) in the Gulf Stream is balanced by a southward flux of dissolved organic nitrogen in the recirculation gyre. Oxidation of the dissolved organic matter during its transit of the subtropical gyre supplies the required source of regenerated nitrate to the Gulf Stream and consumes oxygen, consistent with recent observations of oxygen utilization in the Sargasso Sea.

Description: This research was supported by NASA under contract NAG5-534 and NSF under contract OCE-8521685.

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

Description: Analysis of modern and historical observations demonstrates that the temperature of the intermediate-depth (150–900 m) Atlantic water (AW) of the Arctic Ocean has increased in recent decades. The AW warming has been uneven in time; a local 1°C maximum was observed in the mid-1990s, followed by an intervening minimum and an additional warming that culminated in 2007 with temperatures higher than in the 1990s by 0.24°C. Relative to climatology from all data prior to 1999, the most extreme 2007 temperature anomalies of up to 1°C and higher were observed in the Eurasian and Makarov Basins. The AW warming was associated with a substantial (up to 75–90 m) shoaling of the upper AW boundary in the central Arctic Ocean and weakening of the Eurasian Basin upper-ocean stratification. Taken together, these observations suggest that the changes in the Eurasian Basin facilitated greater upward transfer of AW heat to the ocean surface layer. Available limited observations and results from a 1D ocean column model support this surmised upward spread of AW heat through the Eurasian Basin halocline. Experiments with a 3D coupled ice–ocean model in turn suggest a loss of 28–35 cm of ice thickness after 50 yr in response to the 0.5 W m−2 increase in AW ocean heat flux suggested by the 1D model. This amount of thinning is comparable to the 29 cm of ice thickness loss due to local atmospheric thermodynamic forcing estimated from observations of fast-ice thickness decline. The implication is that AW warming helped precondition the polar ice cap for the extreme ice loss observed in recent years.

Description: This study was supported by JAMSTEC (IP and VI), NOAA (IP, VI, and ID), NSF (IP,VA,VI, ID, JT, andMS),NASA(IP andVI), BMBF (ID), and UK NERC (SB) grants.

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

Description: A freshwater plume often forms when a river or an estuary discharges water onto the continental shelf. Freshwater plumes are ubiquitous features of the coastal ocean and usually leave a striking signature in the coastal hydrography. The present study combines both hydrographic data and idealized numerical simulations to examine how ambient currents and winds influence the transport and mixing of plume waters. The first portion of the thesis considers the alongshore transport of freshwater using idealized numerical simulations. In the absence of any ambient current, the downstream coastal current only carries a fraction of the discharged fresh water; the remaining fraction recirculates in a continually growing "bulge" of fresh water in the vicinity of the river mouth. The fraction of fresh water transported in the coastal current is dependent on the source conditions at the river mouth. The presence of an ambient current augments the transport in the plume so that its freshwater transport matches the freshwater source. For any ambient current in the same direction as the geostrophic coastal current, the plume will evolve to a steady-state width. A key result is that an external forcing agent is required in order for the entire freshwater volume discharged by a river to be transported as a coastal current. The next section of the thesis addresses the wind-induced advection of a river plume, using hydrographic data collected in the western Gulf of Maine. The observations suggest that the plume's cross-shore structure varies markedly as a function of fluctuations in alongshore wind forcing. Consistent with Ekman dynamics, upwelling favorable winds spread the plume offshore, at times widening it to over 50 km in offshore extent, while downwelling favorable winds narrow the plume width to a few Rossby radii. Near-surface current meters show significant correlations between cross-shore currents and alongshore wind stress, consistent with Ekman theory. Estimates of the terms in the alongshore momentum equation calculated from moored current meter arrays also indicate an approximate Ekman balance within the plume. A significant correlation between alongshore currents and alongshore wind stress suggests that interfacial drag may be important. The final section of the thesis is an investigation of the advection and mixing of a surface-trapped river plume in the presence of an upwelling favorable wind stress, using a three-dimensional model in a simple, rectangular domain. Model simulations demonstrate that the plume thins and is advected offshore by the crossshore Ekman transport. The thinned plume is susceptible to significant mixing due to the vertically sheared horizontal currents. The first order plume response is explained by Ekman dynamics and a Richardson number mixing criterion. Under a sustained wind event, the plume evolves to a quasi-steady, uniform thickness. The rate of mixing slowly decreases for longer times as the stratification in the plume weakens, but mixing persists under a sustained upwelling wind until the plume is destroyed. Mixing is most intense at the seaward plume front due to an Ekman straining mechanism in which the advection of cross-shore salinity gradients balances vertical mixing. The mean mixing rate observed in the plume is consistent with the mixing power law suggested by previous studies of I-D mixing, in spite of the two-dimensional dynamics driving the mixing in the plume.

Description: This research was funded by a National Science Foundation graduate fellowship, and Gulf of Maine Regional Marine Research Program grants UM-S227 and UM-S276.

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 2011

Description: Eastern oceanic boundary currents are subject to hydrodynamic instability, generate small scale features that are visible in satellite images and may radiate westward into the interior, where they can be modified by the large-scale circulations. This thesis studies the stability of an eastern boundary current with and without the large-scale flow influence in an idealized framework represented by barotropic quasi-geostrophic dynamics. The linear stability analysis of a meridional current with a continuous velocity profile shows that meridional eastern and western boundary currents support a limited number of radiating modes with long meridional and zonal wavelengths and small growth rates. However, the linearly stable, long radiating modes of an eastern boundary current can become nonlinearly unstable by resonating with short trapped unstable modes. This phenomenon is clearly demonstrated in the weakly nonlinear simulations. Results suggest that linearly stable longwave modes deserve more attention when the radiating instability of a meridional boundary current is considered. A large-scale flow affects the short trapped unstable mode and long radiating mode through different mechanisms. The large-scale flow modifies the structure of the boundary current to stabilize or destabilize the unstable modes, leading to a meridionally localized maximum in the perturbation kinetic energy field. The shortwave mode is accelerated or decelerated by the meridional velocity adjustment of the large-scale flow to have an elongated or a squeezed meridional structure, which is confirmed both in a linear WKB analysis and in nonlinear simulations. The squeezed or elongated unstable mode detunes the nonlinear resonance with the longwave modes, which then become less energetic. These two modes show different meridional structures in kinetic energy field because of the different mechanisms. In spite of the model simplicity, these results can potentially explain the formation of the zonal jets observed in altimeter data, and indicate the influence of the large-scale wind-driven circulation on eastern boundary upwelling systems in the real ocean. Studies with more realistic configurations remain future challenges.

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 January 1999

Description: Today, deep waters produced in the North Atlantic are exported through the western South Atlantic. Antarctic intermediate water AAJW also enters the Atlantic in this region. Circumpolar deep water (CDW) fills the depths below AAIW and above and below northern source waters. A depth transect of cores from 1567-3909 m water depth in the western South Atlantic are ideally located to monitor inter-ocean exchange of deep water, and variations in the relative strength of northern versus southern source water production. Last glacial maximum (LGM) Cd/Ca and δ13C data indicate a nutrient-depleted intermediate-depth water mass. In the mid-depth western South Atlantic, a simple conversion of LGM δ13C data suggests significantly less nutrient enrichment than LGM Cd/Ca ratios, but Cd/Ca and δ13C data can be reconciled when plotted in CdW/δ13C space. Paired LGM Cd/Ca and δ13C data from mid-depth cores suggest increasingly nutrient rich waters below 2000 m, but do not require an increase in Southern Ocean water contribution relative to today. Cd/Ca data suggest no glacial-interglacial change in the hydrography of the deepest waters ofthe region. To maintain relatively low Cd/Ca ratios low nutrients in the deepest western South Atlantic waters, and in CDW in general, during the LGM requires an increased supply ofnutrient-depleted glacial North Atlantic intermediate water (GNA1W) and/or nutrient-depleted glacial Subantarctic surface waters to CDW to balance reduced NADW contribution to CDW. LGM Cd/Ca and δ13C data suggest strong GNA1W influence in the western South Atlantic which in turn implies export of GNAIW from the Atlantic, and entrainment of GNA1W into the Antarctic Circumpolar current.

Description: Author Posting. © American Meteorological Society, 2011. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Atmospheric and Oceanic Technology 28 (2011): 1539–1553, doi:10.1175/JTECH-D-11-00001.1.

Description: Turbulent Reynolds stresses are now routinely estimated from acoustic Doppler current profiler (ADCP) measurements in estuaries and tidal channels using the variance method, yet biases due to surface gravity waves limit its use in the coastal ocean. Recent modifications to this method, including spatially filtering velocities to isolate the turbulence from wave velocities and fitting a cospectral model to the below-wave band cospectra, have been used to remove this bias. Individually, each modification performed well for the published test datasets, but a comparative analysis over the range of conditions in the coastal ocean has not yet been performed. This work uses ADCP velocity measurements from five previously published coastal ocean and estuarine datasets, which span a range of wave and current conditions as well as instrument configurations, to directly compare methods for estimating stresses in the presence of waves. The computed stresses from each were compared to bottom stress estimates from a quadratic drag law and, where available, estimates of wind stress. These comparisons, along with an analysis of the cospectra, indicated that spectral fitting performs well when the wave climate is wide-banded and/or multidirectional as well as when instrument noise is high. In contrast, spatial filtering performs better when waves are narrow-banded, low frequency, and when wave orbital velocities are strong relative to currents. However, as spatial filtering uses vertically separated velocity bins to remove the wave bias, spectral fitting is able to resolve stresses over a larger fraction of the water column.

Description: J. Rosman acknowledges funding from the National Science Foundation (OCE-1061108).

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 2012

Description: Observations from a three-year field program on the inner shelf south of Martha's Vineyard, MA and a numerical model are used to describe the effect of stratification on inner shelf circulation, transport, and sediment resuspension height. Thermal stratification above the bottom mixed layer is shown to cap the height to which sediment is resuspended. Stratification increases the transport driven by cross-shelf wind stresses, and this effect is larger in the response to offshore winds than onshore winds. However, a one-dimensional view of the dynamics is not sufficient to explain the relationship between circulation and stratification. An idealized, cross-shelf transect in a numerical model (ROMS) is used to isolate the effects of stratification, wind stress magnitude, surface heat flux, cross-shelf density gradient, and wind direction on the inner shelf response to the cross-shelf component of the wind stress. In well mixed and weakly stratified conditions, the cross-shelf density gradient can be used to predict the transport efficiency of the cross-shelf wind stress. In stratified conditions, the presence of an along-shelf wind stress component makes the inner shelf response to cross-shelf wind stress strongly asymmetric.

Description: This work was supported through National Science Foundation grant no. OCE-0548961, the WHOI Academic Programs Office, and the WHOI Coastal Ocean Institute.

Description: Author Posting. © American Meteorological Society, 2012. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 25 (2012): 5153–5172, doi:10.1175/JCLI-D-11-00463.1.

Description: Atlantic meridional overturning circulation (AMOC) variability is documented in the Community Climate System Model, version 4 (CCSM4) preindustrial control simulation that uses nominal 1° horizontal resolution in all its components. AMOC shows a broad spectrum of low-frequency variability covering the 50–200-yr range, contrasting sharply with the multidecadal variability seen in the T85 × 1 resolution CCSM3 present-day control simulation. Furthermore, the amplitude of variability is much reduced in CCSM4 compared to that of CCSM3. Similarities as well as differences in AMOC variability mechanisms between CCSM3 and CCSM4 are discussed. As in CCSM3, the CCSM4 AMOC variability is primarily driven by the positive density anomalies at the Labrador Sea (LS) deep-water formation site, peaking 2 yr prior to an AMOC maximum. All processes, including parameterized mesoscale and submesoscale eddies, play a role in the creation of salinity anomalies that dominate these density anomalies. High Nordic Sea densities do not necessarily lead to increased overflow transports because the overflow physics is governed by source and interior region density differences. Increased overflow transports do not lead to a higher AMOC either but instead appear to be a precursor to lower AMOC transports through enhanced stratification in LS. This has important implications for decadal prediction studies. The North Atlantic Oscillation (NAO) is significantly correlated with the positive boundary layer depth and density anomalies prior to an AMOC maximum. This suggests a role for NAO through setting the surface flux anomalies in LS and affecting the subpolar gyre circulation strength.

Description: The CCSM project is supported by NSF and the Office of Science (BER) of the U.S. Department of Energy. SGY and YOK were supported by the NOAA Climate Program Office under Climate Variability and Predictability Program Grants NA09OAR4310163 and NA10OAR4310202, respectively.

Description: 2013-02-01

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

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

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

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

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

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

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

Description: Meridional velocity, mass, and heat transport in the equatorial oceans are difficult to estimate because of the nonapplicability of the geostrophic balance. For this purpose a steady-state model is utilized in the equatorial Indian Ocean using NCEP wind stress and temperature and salinity data from the World Ocean Atlas 2005 (WOA05) and Argo. The results show a Somali Current flowing to the south during the winter monsoon carrying −11.5 ± 1.3 Sv (1 Sv ≡ 106 m3 s−1) and −12.3 ± 0.3 Sv from WOA05 and Argo, respectively. In the summer monsoon the Somali Current reverses to the north transporting 16.8 ± 1.2 Sv and 19.8 ± 0.6 Sv in the WOA05 and Argo results. Transitional periods are considered together and in consequence, there is not a clear Somali Current present in this period. Model results fit with in situ measurements made around the region, although Argo data results are quite more realistic than WOA05 data results.

Description: This study has been partly funded by the MOC Project (CTM 2008- 06438) and the Spanish contribution to the Argo network (AC2009 ACI2009-0998), financed by the Spanish Government and Feder.

Description: 2013-06-01

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

Description: Between 2002 and 2011 a single mooring was maintained in the core of the Pacific Water boundary current in the Alaskan Beaufort Sea near 152° W. Using velocity and hydrographic data from six year-long deployments during this time period, we examine the interannual variability of the current. It is found that the volume, heat, and freshwater transport have all decreased drastically over the decade, by more than 80%. The most striking changes have occurred during the summer months. Using a combination of weather station data, atmospheric reanalysis fields, and concurrent shipboard and mooring data from the Chukchi Sea, we investigate the physical drivers responsible for these changes. It is demonstrated that an increase in summertime easterly winds along the Beaufort slope is the primary reason for the drop in transport. The intensification of the local winds has in turn been driven by a strengthening of the summer Beaufort High in conjunction with a deepening of the summer Aleutian Low. Since the fluxes of mass, heat, and freshwater through Bering Strait have increased over the same time period, this raises the question as to the fate of the Pacific water during recent years and its impacts. We present evidence that more heat has been fluxed directly into the interior basin from Barrow Canyon rather than entering the Beaufort shelfbreak jet, and this is responsible for a significant portion of the increased ice melt in the Pacific sector of the Arctic Ocean.

Description: The majority of the data for this project was funded by grant # ARC-0856244 from the O ce of Polar Programs of the National Science Foundation. My time at WHOI was funded by the United States Navy, the National Science Foundation Graduate Research Fellowship Program and the WHOI Academic Programs O ffice.

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

Description: The equatorial Pacific heat flow low, a major oceanic geothermal anomaly centered on the equatorial sediment bulge, was investigated using deeply penetrating heat flow probes (6-11 meters penetration) within three detailed surveys (400 km2) and along over 10,000 km of continuous seismic profiles (CSP). Previous heat flow measurements in this region defined a broad region characterized by a heat flux well below 1 HFU. We report 98 new measurements collected during cruises PLEIADES 3 and KNORR 73-4 that verify the anomalous nature of the heat flux and also define non-linear temperature gradients (concave down). Temperature field disturbances due to perturbations of a purely conductive heat transport regime are incapable of suitably explaining either of these observations . A simple model incorporating heat transport by both conduction and fluid convection through the sediments fits the observations. A volume flux of (hydrothermal) fluid in the range of 10-6 to 10-5 cm3/sec/cm2 (0.1 liter/yr/cm2) is required. The sense of the flow for all measurements exhibiting non-linear gradients is upward out of the sediment column; no evidence for the recharging of the system was observed. Investigation of a well-defined boundary of the low zone at 4°N and 114°W showed a transition from low and variable heat flow to values compatible with thermal models that correlated with a change in the nature of the basement from rough to smooth. A few outcrops occur in the area of rough basement, but otherwise the region is well-sedimented (greater than 200 meters). Measurements within a detailed survey centered at this transition showed a dramatic increase in heat flow from 1.21 HFU to values greater than 3 HFU over a horizontal distance of 10km. A similar transition from non-linear to linear temperature gradients was not observed as nearly every measurement was non-linear. Heat flow measurements located in well-sedimented, outcrop-free areas (A environments) were associated with linear gradients and a heat flux greater than 1 HFU, however, several of these values were well below the theoretical heat flow for the appropriate age crust. Values measured in environments other than A exhibited variable heat flow and non-linear gradients. The average value of measurements located in A environments within the equatorial Pacific heat flow low was 1.37±0.27 HFU. The previously reported average was 0.92±0.48 HFU based on several measurements from L-DGO cruise VEMA 24-3. The average heat flow measured at a survey located outside the low heat flow zone on crust of 55 ±5 m.a. was 1.76 ±0.30 HFU which is in good agreement with the theoretical value of 1.60. The measurements in this survey were not located in A environments suggesting that crustal convection has ceased or is greatly attenuated within crust of this age. Error analysis of the geothermal data reduction using the convective/conductive heat transport model suggests that the volume flux parameter is sensitive to temperature measurement errors greater than a few millidegrees. Volume fluxes less than 10-7 cm/sec are difficult to distinguish from the purely conductive case assuming instrumental accuracies of 0.001°C. Resolution of the volume flux deteriorates as heat flow decreases and is poor for values less than 0.5 HFU. A detailed survey located within the low zone confirmed previous measurements of low heat flow, however, due to the low value of heat flow (about 0.5 HFU) the small-scale variability could not be clearly defined.

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

Description: This is the second part of a two-part investigation of a coupled wind and wave model that includes the enhanced form drag of breaking waves. The model is based on the wave energy balance and the conservation of air-side momentum and energy. In Part I, coupled nonlinear advance–delay differential equations were derived, which govern the wave height spectrum, the distribution of breaking waves, and vertical air side profiles of the turbulent stress and wind speed. Numeric solutions were determined for mature seas. Here, numeric solutions for a wide range of wind and wave conditions are obtained, including young, strongly forced wind waves. Furthermore, the “spatial sheltering effect” is introduced so that smaller waves in airflow separation regions of breaking longer waves cannot be forced by the wind. The solutions strongly depend on the wave height curvature spectrum at high wavenumbers (the “threshold saturation level”). As the threshold saturation level is reduced, the effect of breaking waves becomes stronger. For young strongly forced waves (laboratory conditions), breaking waves close to the spectral peak dominate the wind input and previous solutions of a model with only input to breaking waves are recovered. Model results of the normalized roughness length are generally consistent with previous laboratory and field measurements. For field conditions, the wind stress depends sensitively on the wave height spectrum. The spatial sheltering may modify the number of breaking shorter waves, in particular, for younger seas.

Description: This work was supported by the U.S. National Science Foundation (Grant OCE- 0526177) and the U.S. Office of Naval Research (Grant N00014–06–10729).

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

Description: Analyses of current time series longer than 200 days from 33 sites over the Middle Atlantic Bight continental shelf reveal a consistent mean circulation pattern. The mean depth-averaged flow is equatorward, alongshelf, and increases with increasing water depth from 3 cm s−1 at the 15-m isobath to 10 cm s−1 at the 100-m isobath. The mean cross-shelf circulation exhibits a consistent cross-shelf and vertical structure. The near-surface flow is typically offshore (positive, range −3 to 6 cm s−1). The interior flow is onshore and remarkably constant (−0.2 to −1.4 cm s−1). The near-bottom flow increases linearly with increasing water depth from −1 cm s−1 (onshore) in shallow water to 4 cm s−1 (offshore) at the 250-m isobath over the slope, with the direction reversal near the 50-m isobath. A steady, two-dimensional model (no along-isobath variations in the flow) reproduces the main features of the observed circulation pattern. The depth-averaged alongshelf flow is primarily driven by an alongshelf pressure gradient (sea surface slope of 3.7 × 10−8 increasing to the north) and an opposing mean wind stress that also drives the near-surface offshore flow. The alongshelf pressure gradient accounts for both the increase in the alongshelf flow with water depth and the geostrophic balance onshore flow in the interior. The increase in the near-bottom offshore flow with water depth is due to the change in the relative magnitude of the contributions from the geostrophic onshore flow that dominates in shallow water and the offshore flow driven by the bottom stress that dominates in deeper water.

Description: This research was funded by Ocean Sciences Division of the National Science Foundation under Grants OCE-820773, OCE-841292, and OCE-848961.

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

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

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

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

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

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

Description: Author Posting. © American Meteorological Society, 2012. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 25 (2012): 343–349, doi:10.1175/JCLI-D-11-00059.1.

Description: The Equatorial Undercurrent (EUC) is a major component of the tropical Pacific Ocean circulation. EUC velocity in most global climate models is sluggish relative to observations. Insufficient ocean resolution slows the EUC in the eastern Pacific where nonlinear terms should dominate the zonal momentum balance. A slow EUC in the east creates a bottleneck for the EUC to the west. However, this bottleneck does not impair other major components of the tropical circulation, including upwelling and poleward transport. In most models, upwelling velocity and poleward transport divergence fall within directly estimated uncertainties. Both of these transports play a critical role in a theory for how the tropical Pacific may change under increased radiative forcing, that is, the ocean dynamical thermostat mechanism. These findings suggest that, in the mean, global climate models may not underrepresent the role of equatorial ocean circulation, nor perhaps bias the balance between competing mechanisms for how the tropical Pacific might change in the future. Implications for model improvement under higher resolution are also discussed.

Description: KBK gratefully acknowledges the J. Lamar Worzel Assistant Scientist Fund. GCJ is supported by NOAA’s Office of Oceanic and Atmospheric Research. RM gratefully acknowledges the generous support and hospitality of the Divecha Centre for Climate Change and CAOS at IISc, Bangalore, and partial support by NASA PO grants.

Description: 2012-07-01

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

Description: Reciprocal acoustic transmissions made in a region just south of the Gulf Stream are analyzed to determine the structure and variability of temperature, current velocity, and vorticity fields at the northern extent of the southern recirculation gyre. For ten months (November, 1988 through August, 1989), a pentagonal array of tomographic transceivers was situated in a region centered at 38°N, 55°W as part of the eastern array of the SYNOP (SYNoptic Ocean Prediction) Experiment. The region of focus is one rich in mesoscale energy, with the influence of local Gulf Stream meandering and cold-core ring activity strikingly evident. Daily-averaged acoustic transmissions yielded travel times which were inverted to obtain estimates of range-averaged temperature and current velocity fields, and area-averaged relative vorticity fields. The acoustically determined estimates are consistent with nearby current meter measurements and satellite infrared imagery. The signature of cold-core rings is clearly evident in the sections. Spectral estimates of the fields are dominated by motions with periodicities ranging from 32-128 days. Second-order statistics, such as eddy kinetic energies, and heat and momentum fluxes, are also estimated. The integrating nature of the tomographic measurement has been exploited to shed some light on the radiation of eddy energy from the Gulf Stream. The Eliassen-Palm flux diagnostic has been applied to an investigation of wave radiation from the Gulf Stream. Results of the diagnosis suggest that the Gulf Stream itself is the source of wave energy radiating into the far field and found in the interior of the North Atlantic subtropical gyre.

Description: This research was carried out under Office of Naval Research (ONR) University Research Initiative contract N00014-86-K-0751 and ONR contract N00014- 90-J-1481. Construction of the tomographic instruments was supported by grants and contracts with MIT: National Science Foundation grant OCE 85-12430 and by ONR. The field work was supported by ONR under contract N00014-85-G-0241 (Secretary of the Navy Professorship (C. Wunsch)).

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 1989

Description: The general theme of this thesis is the study of systematic mathematical techniques for determining the ocean circulation from classical hydrographic data. Two aspects of this theme are analyzed. The first is finding an efficient representation of hydrographic structure so as to make it most useful and informative. The second is application of inverse methods to the data to determine ocean circulation. Both subjects are examined in the North Atlantic Ocean. The efficient representation is examined in terms of empirical orthogonal functions (EOFs) among the variations in vertical hydrographic profiles. The data used are of a new set of high quality hydrography, all obtained in the early 1980s. Common EOFs are examined among temperature, salinity, oxygen, phosphate, silicate, and nitrate. The EOFs identify a fundamental simplicity in the spatial distributions of t hese properties. Although the volume of numbers involved in the raw data is large, the significant degrees of freedom are only six in space and two among the six properties; temperature and salinity are represented by one mode, while the nutrients by another. The modal structure reflects some underlying simplicity in ocean physics. EOFs form a quantitative basis from which models of the ocean's hydrographic structure can be constructed for various degrees of complexities. As for the second aspect, two applications of inverse methods are explored on small regional scales. The first problem addressed concerns the circulation inside a 12° square located in the eastern basin over the axis of the Mediterranean Water tongue. The study is based on an ocean model constructed by mapping the modes identified in the first half of the thesis over the entire North Atlantic Ocean. A combination of box model inverse and β-spiral method is used to determine the geostrophic reference level velocities. The circulation consists of an anticyclonic circulation near the surface, which is part of the eastern half of the wind-driven subtropical gyre. The flow at depth is weak, and is a cyclonic circulation around the core of the Mediterranean Water tongue. In the second inverse problem, we examine a decaying warm-core ring. Observations of a warm-core ring are used to formulate a model for diagnosing the physics of ring change over a two month period. About 30 hydrographic casts and acoustic doppler current measurements are used to generate estimates of an equivalent radially symmetric ring with radial contrasts of stratification, temperature, salinity, azimuthal velocity, angular momentum, and potential vorticity. A series of related models are inverted for the ring circulation and mixing coefficients. The circulation is insensitive to the model details, is well-resolved, and is a radial outflow and upwelling. Eddy coefficients are only partially resolved; determining the mixing with any degree of confidence appears to require a much more elaborate data set than the one available.

Description: This research was funded in part by the Office of Naval Research (Secretary of the Navy Chair) and the National Science Foundation under grant OCE 85-21685.

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: A convection experiment was done with a rotating rectangular tank as a model of oceanic meridional overturning circulation. Heat flux was fixed at one bottom end of the tank using an electrical heater. Temperature was fixed at the other end using a cooling plate. All other boundaries were insulated. The cross sections of temperature field were made at several locations. In equilibrium, the heat input to the fluid H was the same as the meridional heat flux (heat flux from the source to the sink), so it was possible to find a scaling law relating H to the temperature difference across the tank ΔT and rotation rate f. The experimental result suggests that the meridional heat transport in the experiment was mostly due to geostrophic flows with a minor correction caused by the bottom friction. If there was no friction, the scaling law from the experiment resembles the one verified in part in the numerical model by Bryan and Cox (1967). Flow visualization and temperature sections showed that there were meridional geostrophic currents that transported heat. When the typical values of the North Atlantic are introduced, the geostrophic scaling law predicts meridional heat flux comparable to that estimated in the North Atlantic when the vertical eddy diffusivity of heat is about 1cm2s-1. Naturally, this experiment is a only crude model of the oceanic convective circulation. We do not claim that the geostrophic scaling applies in detail to the oceans, however, it may have some important use in climate modeling. For example, almost all existing box models and two-dimensional numerical models of ocean circulation use a frictional scaling law for buoyancy transport. A box model with the geostrophic scaling law is shown to be more robust to a change in the boundary forcing so that it is less likely to have a thermohaline catastrophic transition under the present conditions. It is also shown that a restoring boundary condition for salinity introduces stability to a thermal mode circulation, unless the restoring time for salinity is several orders of magnitude larger than that for temperature.

Description: This study has been funded by NSF grant number OCE92-01464 and Korean Government Overseas Scholarship Grant.

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: Large-scale thermal forcing and freshwater fluxes play an essential role in setting temperature and salinity in the ocean. A number of recent estimates of the global oceanic freshwater balance as well as the global oceanic surface net heat flux are used to investigate the effects of heat- and freshwater forcing at the ocean surface. Such forcing induces changes in both density and density-compensated temperature and salinity changes (’spice’). The ratio of the relative contributions of haline and thermal forcing in the mixed layer is maintained by large-scale surface fluxes, leading to important consequences for mixing in the ocean interior. In a stratified ocean, mixing processes can be either along lines of constant density (isopycnal) or across those lines (diapycnal). The contribution of these processes to the total mixing rate in the ocean can be estimated from the large-scale forcing by evaluating the production of thermal variance, salinity variance and temperature-salinity covariance. Here, I use new estimates of surface fluxes to evaluate these terms and combine them to generate estimates of the production of density and spice variance under the assumption of a linear equation of state. As a consequence, it is possible to estimate the relative importance of isopycnal and diapycnal mixing in the ocean. While isopycnal and diapycnal processes occur on very different length scales, I find that the surface-driven production of density and spice variance requires an approximate equipartition between isopycnal and diapycnal mixing in the ocean interior. In addition, consideration of the full nonlinear equation of state reveals that surface fluxes require an apparent buoyancy gain (expansion) of the ocean, which allows an estimate of the amount of contraction on mixing due to cabbeling in the ocean interior.

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

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: Internal tides are internal waves of tidal period generated by tidal currents flowing over submarine topography. Tall ridges that are nominally two-dimensional (2-D) are sites of particularly strong generation. The subsequent dissipation of internal tides contributes to ocean mixing, thereby playing an important role in the circulation of the ocean. Strong internal tides can also evolve into internal wave solitons, which affect acoustic communication, offshore structures and submarine navigation. This thesis addresses the generation of internal tides by tall submarine ridges using a combined analytical and experimental approach. The first part of the thesis is an experimental investigation of a pre-existing Green function formulation for internal tide generation by a tall symmetric ridge in a uniform density stratification. A modal decomposition technique was developed to characterize the structure of the experimental wave fields generated by 2D model topographies in a specially configured wave tank. The theory accurately predicts the low mode structure of internal tides, and reasonably predicts the conversion rate of internal tides infinite tidal excursion regimes, for which the emergence of non-linearities was notable in the laboratory. In the second part of the thesis, the Green function method is advanced for asymmetric and multiple ridges in weakly non-uniform stratifications akin to realistic ocean situations. A preliminary investigation in uniform stratification with canonical asymmetric and double ridges reveals asymmetry in the internal tide that can be very sensitive to the geometric configuration. This approach is then used with realistic topography and stratification data to predict the internal tide generated by the ridges at Hawaii and at the Luzon Strait. Despite the assumption of two-dimensionality, there is remarkably good agreement between field data, simulations and the new theory for the magnitude, asymmetry and modal content of the internal tide at these sites. The final part of the thesis investigates the possibility of internal wave attractors in the valley of double-ridge configurations. A one-dimensional map is developed to identify the existence and stability of attractors as a function of the ridge geometry. The Green function method is further advanced to include a viscous correction to balance energy focusing and dissipation along an attracting orbit of internal wave rays, and very good agreement is obtained between experiment and theory, even in the presence of an attractor.

Description: My Ph. D. and the work in this thesis have been generously funded by the National Science Foundation under grants OCE 0645529 and OCE 04-25283 and the Office of Naval Research under grants N00014-08-0390, N00014-05-1-0573 and N00014-09-0282.

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 1987

Description: This thesis consists of two loosely related problems, both of which analyze some consequences of the failure of Sverdrup relation. In the first part, Chapters 2 and 3, the Sverdrup relation is invalidated because substantial flow is obtained at the bottom where topography exists. The eddies play an essential role in transfering momentum vertically from the surface, where the forcing is applied, to the bottom, which is otherwise unforced. If the topography has a structure in the longitudinal direction, then the inviscid theory predicts the occurence of strong jets in the interior of the model ocean. According to the structure of the topography these internal jets can occur in both vertically homeogenous and baroclinic oceans. If the topographic slope changes sign, then one kind of jets is observed both in stratified and in homogeneous oceans. This phenomenon is robust to moderate amounts of dissipation and is not disturbed by the occurrence of recirculating gyres within the basin. If the topographic slope is constant, then another kind of internal jets is observed, and it occurs in stratified models only. I was unable to observe this kind of jets in the presence of weak dissipation. The reason for this failure is twofold: on one hand friction, especially interfacial friction, tends to make the flow more barotropic (and we believe that indeed this is one of the processes that the eddies accomplish in a stratified fluid) and therefore the phenomena that rely strongly on baroclinicity are discouraged. On the other hand, reduction of the dissipation leads to the onset of a strong recirculating, inertial gyre which, although confined in space, affects the global properties of the flow. In the second part of the thesis (Chapters 4 and 5) I developed a simple model of the recirculating, inertial gyre. Again the dynamics of this feature are far from being in Sverdrup balance. In this case inertia is responsible for the failure of Sverdrup relation, together with the eddy field which provides a mean for transfering momentum vertically and laterally into regions away from where the forcing is applied. In this model there is no direct forcing in the recirculation region, and the input of momentum is confined to the boundary currents surrounding the gyre, for example the separated Gulf Stream. One of the results of the recirculation model is the prediction of its transport. It is shown that most of the transport is depth independent, i.e. it can be calculated without detailed knowledge of the density structure of the ocean. It is also shown that the "barotropic" part of the transport increases as the cube of the meridional extent of the gyre.

Description: The thesis work has been supported by a National Foundation grant from the Office of Atmospheric Sciences.

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

Description: Based on the classical Ekman layer theory, a simple analytical solution of the steady flow induced by a stationary hurricane in a homogenous ocean is discussed. The model consists of flow converging in an inward spiral in the deeper layer and diverging in the upper layer. The simple analytical model indicates that both the upwelling flux and the horizontal transport increase linearly with increasing radius of maximum winds. Furthermore, they both have a parabolic relationship with the maximum wind speed. The Coriolis parameter also affects the upwelling flux: the response to a hurricane is stronger at low latitudes than that at middle latitudes. Numerical solutions based on a regional version of an ocean general circulation model are similar to the primary results obtained through the analytical solution. Thus, the simplifications made in formulating the analytical solution are reasonable. Although the analytical solution in this paper is sought for a rather idealized ocean, it can help to make results from the more complicated numerical model understandable. These conceptual models provide a theoretical limit structure of the oceanic response to a moving hurricane over a stratified ocean.

Description: ZML was supported by CAS (kzcx2-yw-226, SQ200813), NSFC under Grants 40906009 and 40776008, and Dr. Xiaodong Shang through ‘‘100 Talents Program’’ of CAS.

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

Description: Recent work by S. Lentz et al. documents offshore transport in the inner shelf due to a wave-driven return flow associated with the Hasselmann wave stress (the Stokes–Coriolis force). This analysis is extended using observations from the central Oregon coast to identify the wave-driven return flow present and quantify the potential bias of wind-driven across-shelf exchange by unresolved wave-driven circulation. Using acoustic Doppler current profiler (ADCP) measurements at six stations, each in water depths of 13–15 m, observed depth-averaged, across-shelf velocities were generally correlated with theoretical estimates of the proposed return flow. During times of minimal wind forcing, across-shelf velocity profiles were vertically sheared, with stronger velocities near the top of the measured portion of the water column, and increased in magnitude with increasing significant wave height, consistent with circulation due to the Hasselmann wave stress. Yet velocity magnitudes and vertical shears were stronger than that predicted by linear wave theory, and more similar to the stratified “summer” velocity profiles described by S. Lentz et al. Additionally, substantial temporal and spatial variability of the wave-driven return flow was found, potentially due to changing wind and wave conditions as well as local bathymetric variability. Despite the wave-driven circulation found, subtracting estimates of the return flow from the observed across-shelf velocity had no significant effect on estimates of the across-shelf exchange due to along-shelf wind forcing at these water depths along the Oregon coast during summer.

Description: This work was performed with the Partnership for Interdisciplinary Studies of Coastal Oceans (PISCO), funded primarily by the Gordon and Betty Moore Foundation and David and Lucile Packard Foundation. SL acknowledges support from NSF Ocean Science Grant #OCE-0548961. AK acknowledges support from the WHOI Coastal Ocean Institute Fellowship.

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

Description: As a driving parameter is slowly altered, thermohaline ocean circulation models show either a smooth evolution of a mode of flow or an abrupt transition of temperature and salinity fields from one mode to another. An abrupt transition might occur at one value or over a range of the driving parameter. The latter has hysteresis because the mode in this range depends on the history of the driving parameter. Although assorted ocean circulation models exhibit abrupt transitions, such transitions have not been directly observed in the ocean. Therefore, laboratory experiments have been conducted to seek and observe actual (physical) abrupt thermohaline transitions. An experiment closely duplicating Stommel’s box model possessed abrupt transitions in temperature and salinity with distinct hysteresis. Two subsequent experiments with more latitude for internal circulation in the containers possessed abrupt transitions over a much smaller range of hysteresis. Therefore, a new experiment with even more latitude for internal circulation was designed and conducted. A large tank of constantly renewed freshwater at room temperature had a smaller cavity in the bottom heated from below with saltwater steadily pumped in. The cavity had either a salt mode, consisting of the cavity filled with heated salty water with an interface at the cavity top, or a temperature mode, in which the heat and saltwater were removed from the cavity by convection. There was no measurable hysteresis between the two modes. Possible reasons for such small hysteresis are discussed.

Description: Support is gratefully acknowledged from the Woods Hole Oceanographic Institution Climate Change Institute, the National Science Foundation, Physical Oceanography Section under Grant OCE-0081179, and the Paul M. Fye Chair of the Woods Hole Oceanographic Institution.

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

Description: Measurements collected in the York River estuary, Virginia, demonstrate the important impact that tidal and lateral asymmetries in turbulent mixing have on the tidally averaged residual circulation. A reduction in turbulent mixing during the ebb phase of the tide caused by tidal straining of the axial density gradient results in increased vertical velocity shear throughout the water column during the ebb tide. In the absence of significant lateral differences in turbulent mixing, the enhanced ebb-directed transport caused by tidal straining is balanced by a reduction in the net seaward-directed barotropic pressure gradient, resulting in laterally uniform two-layer residual flow. However, the channel–shoal morphology of many drowned river valley estuaries often leads to lateral gradients in turbulent mixing. Tidal straining may then lead to tidal asymmetries in turbulent mixing near the deeper channel while the neighboring shoals remain relatively well mixed. As a result, the largest lateral asymmetries in turbulent mixing occur at the end of the ebb tide when the channel is significantly more stratified than the shoals. The reduced friction at the end of ebb delays the onset of the flood tide, increasing the duration of ebb in the channel. Conversely, over the shoal regions where stratification is more inhibited by tidal mixing, there is greater friction and the transition from ebb to flood occurs more rapidly. The resulting residual circulation is seaward over the channel and landward over the shoal. The shoal–channel segregation of this barotropically induced estuarine residual flow is opposite to that typically associated with baroclinic estuarine circulation over channel–shoal bathymetry.

Description: Support for this research was provided by the National Science Foundation Division of Ocean Sciences Grant OCE- 9984941.

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

Description: The overturning and horizontal circulations of the Labrador Sea are deduced from a composite vertical section across the basin. The data come from the late-spring/early-summer occupations of the World Ocean Circulation Experiment (WOCE) AR7W line, during the years 1990–97. This time period was chosen because it corresponded to intense wintertime convection—the deepest and densest in the historical record—suggesting that the North Atlantic meridional overturning circulation (MOC) would be maximally impacted. The composite geostrophic velocity section was referenced using a mean lateral velocity profile from float data and then subsequently adjusted to balance mass. The analysis was done in depth space to determine the net sinking that results from convection and in density space to determine the diapycnal mass flux (i.e., the transformation of light water to Labrador Sea Water). The mean overturning cell is calculated to be 1 Sv (1 Sv ≡ 106 m3 s−1), as compared with a horizontal gyre of 18 Sv. The total water mass transformation is 2 Sv. These values are consistent with recent modeling results. The diagnosed heat flux of 37.6 TW is found to result predominantly from the horizontal circulation, both in depth space and density space. These results suggest that the North Atlantic MOC is not largely impacted by deep convection in the Labrador Sea.

Description: This work was funded by the National Science Foundation through Grants OCE-0450658 (RP) and OCE-024978 (MS).

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

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

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

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

Description: The steady states of two models of the double-gyre wind-driven ocean circulation are studied. The link between the steady state solutions of the models and their time-mean and low-frequency variability is explored to test the hypothesis that both stable and unstable fixed points influence shape the model's attractor in phase space. The steady state solutions of a barotropic double-gyre ocean model in which the wind-stress curl input of vorticity is balanced primarily by bottom friction are studied. The bifurcations away from a unique and stable steady state are mapped as a function of two nondimensional parameters, (δI,δS), which can be thought of as measuring respectively the relative importance of the nonlinear advection and bottom damping of relative vorticity to the advection of planetary vorticity. A highly inertial branch characterized by a circulation with transports far in excess of those predicted by Sverdrup balance is present over a wide range of parameters including regions of parameter space where other solutions give more realistic flows. For the range of parameters investigated, in the limit of large Reynolds number, δI,δS → ∞, the inertial branch is stable and appears to be unique. This branch is anti-symmetric with respect to the mid-basin latitude like the prescribed wind-stress curl. For intermediate values of δI,δS, additional pairs of mirror image non-symmetric equilibria come into existence. These additional equilibria have currents which redistribute relative vorticity across the line of zero wind-stress curl. This internal redist~ibution of vorticity prevents the solution from developing the large transports that are necessary for the anti-symmetric solution to achieve a global vorticity balance. Beyond some critical Reynolds number, the nonsymmetric solutions are unstable to time-dependent perturbations. Time-averaged solutions in' this parameter regime have transports comparable in magnitude to those of the non-symmetric steady state branch. Beyond a turning point, where the non-symmetric steady state solutions cease to exist, all the computed time-dependent model trajectories converge to the anti-symmetric inertial runaway solution. The internal compensation mechanism which acts through explicitly simulated eddies is itself dependent explicit dissipation parameter. Using the reduced-gravity quasigeostrophic model an investigation of the link between the steady state solutions and the model's low-frequency variability is conducted. If the wind-stress curl is kept anti-symmetric, successive pairs of non-symmetric equilibria come into existence via symmetry-breaking pitchfork bifurcations as the model's biharmonic viscosity is reduced. Succesive pairs of mirror image equilibria have an additional half meander in the jet. The distinct energy levels of the steady state solutiOris can be understood in part by there different inter-gyre fluxes of vorticity. Those solutions with weak inter-gyre fluxes of vorticity have large and energetic recirculation cells which remove excess vorticity through bottom friction. Those solutions with strong inter-gyre fluxes of vorticity have much smaller and ·less energetic recirculation cells. A significant fraction of the variance (30%) of the interface height anomaly can be accounted by four coherent structures which point away from the time-mean state and towards four steady state solutions in phase space. After removing the variance which projects onto the four modes, the remaining variance is reduced predominantly at low-frequencies, showing that these modes are linked to the low-frequency variability of the model. Furthermore, the time-averaged flow fields within distinct energy ranges show distinct patterns which are in turn similar to the distinct steady state solutions.

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 1998

Description: Double-diffusive processes are studied and parameterized, and their impacts on the oceanic thermohaline circulation are investigated by single-hemisphere numerical models and scaling analysis. Scaling analysis on the thermohaline circulation has been done under three types of surface boundary conditions. (a) Under "relaxation" conditions, there is a two-thirds power law dependence of the meridional overturning rate (and the poleward heat transport) on the diapycnal diffusivity. For any given external forcing, there is only one equilibrium state for the thermohaline circulation. (b) Under "flux" boundary conditions, there is a half power law dependence of the meridional overturning rate on the diapycnal diffusivity. Only one mode is possible for given external forcing. (c) Under "mixed" boundary conditions, multiple equilibria become possible. For given thermal forcing, the existence of multiple equilibria depends on the relative contributions of diapycnal diffusivity and the hydrologic forcing. Numerical experiments are implemented to test the above scaling arguments. Consistent results have been obtained under the above three types of boundary conditions. These provide a basis for understanding how the thermohaline circulation depends on the diapycnal diffusivity, which we know is influenced by the double-diffusive processes of "salt fingering" and "diffusive layering" in some parts of the ocean. In order to examine this issue, the double-diffusive processes are parameterized by diapycnal eddy diffusivities for heat and salt that are different and depend on the local density ratio, Rp= αTz/βSz. A background diffusivity is applied to represent turbulent mixing in the stratified environment. The implementation of this double-diffusive - parameterization in numerical models has significant impacts on the thermohaline circulation. (a) Under "relaxation" boundary conditions, the meridional overturning rate and the poleward heat transport are reduced, and water mass properties are also changed. Similar results are obtained under "flux" boundary conditions. (b) Under "mixed" boundary conditions, the critical freshwater flux for the existence of the thermal mode becomes smaller with the double-diffusive parameterization. The extent to which the thermohaline circulation is affected by double-diffusive processes depends on the magnitude of the freshwater forcing.

Description: This thesis is supported by a grant from the Ocean Sciences Division of the National Science Foundation, OCE94-155S9.

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

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

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

Description: 2012-11-01

Description: 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: Eighteen months of sea surface height data from the GEOSAT altimeter along collinear subtracks were analyzed for information on the circulation pattern in the Bering Sea. Seventy subtracks from both ascending and descending orbits, with as many as 35 repeat cycles along each subtrack, were analyzed. Orbit errors were removed from the height data using a least-squares fit to a cubic polynomial, weighted by the inverse of the height variance. Addition of the weights decreased contamination of residual height profiles by the large geoid signal. Composite maps of variability along each track revealed patterns of increased variability in the regions of the documented Bering slope current (BSC) and the proposed western boundary current (WBC); however, no evidence was found of the expected bifurcation of the BSC near the Siberian coast. Past observations of tides in the Bering Sea were reviewed along with a local tide model to detect tidal contributions to the mesoscale sea surface height variability. The tidal analysis suggested that residual tides contributed primarily to the longer wavelengths which were removed in the collinear processing. Examination of the Schwiderski tidal correction proved it to be a sensible correction, reducing the height variance by approximately 60%. Finally, using a Gaussian model for the BSC velocity profile, synthetic residual heights were generated and fit to the actual data to produce estimates of absolute surface geostrophic velocity and transport. Comparisons of mean flow, height fluctuations and seasonal trends across the BSC, the WBC and Bering Strait support the hypothesis that the BSC turns north at Cape Navarin into the WBC which, in turn, is capable of supplying a major part of the transport through the Bering Strait.

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

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

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

Description: This thesis focuses on ocean circulation and atmospheric forcing in the Atlantic Ocean at the Last Glacial Maximum (LGM, 18-21 thousand years before present). Relative to the pre-industrial climate, LGM atmospheric CO2 concentrations were about 90 ppm lower, ice sheets were much more extensive, and many regions experienced significantly colder temperatures. In this thesis a novel approach to dynamical reconstruction is applied to make estimates of LGM Atlantic Ocean state that are consistent with these proxy records and with known ocean dynamics. Ocean dynamics are described with the MIT General Circulation Model in an Atlantic configuration extending from 35°S to 75°N at 1° resolution. Six LGM proxy types are used to constrain the model: four compilations of near sea surface temperatures from the MARGO project, as well as benthic isotope records of δ18O and δ13C compiled byMarchal and Curry; 629 individual proxy records are used. To improve the fit of the model to the data, a least-squares fit is computed using an algorithm based on the model adjoint (the Lagrange multiplier methodology). The adjoint is used to compute improvements to uncertain initial and boundary conditions (the control variables). As compared to previous model-data syntheses of LGM ocean state, this thesis uses a significantly more realistic model of oceanic physics, and is the first to incorporate such a large number and diversity of proxy records. A major finding is that it is possible to find an ocean state that is consistent with all six LGM proxy compilations and with known ocean dynamics, given reasonable uncertainty estimates. Only relatively modest shifts from modern atmospheric forcing are required to fit the LGM data. The estimates presented herein successfully reproduce regional shifts in conditions at the LGM that have been inferred from proxy records, but which have not been captured in the best available LGM coupled model simulations. In addition, LGM benthic δ18O and δ13C records are shown to be consistent with a shallow but robust Atlantic meridional overturning cell, although other circulations cannot be excluded.

Description: Primary support was provided by a National Defense Science and Engineering Graduate Fellowship and two National Science Foundation awards: Award #OCE-0645936: “Beyond the Instrumental Record: the Case of Circulation at the Last Glacial Maximum” and Award #OCE-1060735: “Collaborative Research: Beyond the Instrumental Record - the Ocean Circulation at the Last Glacial Maximum and the de-Glacial Sequence”. Important secondary support came from the National Ocean Partnership Program and the National Aeronautics and Space Administration via the ECCO effort at MIT.

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

Description: Ocean modellers seek to understand the circulation of the oceans, or portions thereof, by developing models of the ocean they can solve. This tractability constraint forces ocean modellers to make choices. Naturally, they hope to make intelligent choices, but whenever a new model is being developed or an existing one extended, the issue of tractability lurks. The large-scale, basin-wide, circulation of the oceans can be divided into two components, classified by their driving force. The wind-driven circulation, whose flow occurs mainly above the thermocline, was first explained qualitatively by Stommel (1948) with a simple, elegant analytical model. The other component of the oceans' circulation, the density-driven, or thermohaline circulation, flows below the thermocline. Again, the first simple analytical model for the deep thermohaline flow was proposed by Stommel (1958) and developed by Stommel and Aarons (1959) whose basic ideas underlie even the most recent conceptual models of the large-scale circulation. The details of the thermohaline circulation and its interaction with the wind-driven circulation in a realistic ocean basin is a problem which is not tractable analytically. This has driven ocean modellers interested in this aspect of the oceans' circulation to numerical models: ocean circulation 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 1994

Description: This work investigates whether large-scale coherent vortex structures driven by wave-current interaction (Langmuir circulation) are responsible for maintaining the oceanic mixed layer. Langmuir circulations dominate the near-surface vertical transport of momentum and density when the characteristic scale for forcing (defined as the Craik-Leibovich instability parameter γCLS) is stronger than the characteristic scale for diffusive decay γdiff. Since the wave-current forcing is concentrated near the surface both terms depend on the cell geometry. Cells with long wavelengths penetrate more deeply into the water column. These cells grow more slowly than the fastest growing mode for most cases, but always dominate the solution in the absence of Coriolis forces. In the presence of Coriolis forces, the horizontal wavelength and thus the depth of penetration are limited. When a cell geometry is found such that γCLS » γdiff, the current profile produced by small-scale diffusion is unstable to Langmuir cells and the cells replace small-scale diffusion as the dominant vertical transport mechanism for momentum and density. The perturbation crosscell shear is predicted to scale as γCLS. Such a scaling is observed during two field experiments. The observed velocity profile during these experiments is more sheared than predicted by a model which implicitly assumes instantaneous mixing by large eddies, but less sheared than predicted by a model which assumes small-scale mixing by near-isotropic turbulence. The latter profile is unstable to Langmuir cells when waves are present. The inclusion of cells driven by wave-current interaction explains the failure of the mixed layer to restratify on two days with high waves and low wind. Wave-current interaction introduces a small but efficient source of energy for transporting density which goes as the surface stress times the Stokes drift.

Description: The Office of Naval Research supported me throughout graduate school, first as an ONR Graduate Fellow. and later as a research assistant under the Surface Waves Processes Program (ONR Grant N00014-90-J-1495).

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: Author Posting. © American Meteorological Society, 2012. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 42 (2012): 2283–2296, doi:10.1175/JPO-D-11-0227.1.

Description: The dynamic influence of thermohaline circulation on wind-driven circulation in the South China Sea (SCS) is studied using a simple reduced gravity model, in which the upwelling driven by mixing in the abyssal ocean is treated in terms of an upward pumping distributed at the base of the upper layer. Because of the strong upwelling of deep water, the cyclonic gyre in the northern SCS is weakened, but the anticyclonic gyre in the southern SCS is intensified in summer, while cyclonic gyres in both the southern and northern SCS are weakened in winter. For all seasons, the dynamic influence of thermohaline circulation on wind-driven circulation is larger in the northern SCS than in the southern SCS. Analysis suggests that the upwelling associated with the thermohaline circulation in the deep ocean plays a crucial role in regulating the wind-driven circulation in the upper ocean.

Description: G. Wang is supported by the National Science Foundation of China (NSFC Grants 41125019, 40725017, and 40976017).D.Chen is supported by grants from the Ministry of Science and Technology (2010DFA21012), the State Oceanic Administration (201105018), and the NSFC (91128204).

Description: 2013-06-01

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

Description: The water mass distribution in the southwestern Barents Sea, the thermohaline structure of the western Barents Sea Polar Front, and the formation of local water masses are described based on an analysis of historical hydrographic data and a recent process-oriented field experiment. This study concentrated on the frontal region between Bj0rn0ya and Hopen Island where Arctic water is found on the Spitzbergen Bank and Atlantic Water in the Bear Island Trough and Hopen Trench. Distributions of Atlantic, Arctic, and Polar Front waters are consistent with topographic control of Atlantic water circulation. Seasonal buoyancy forcing disrupts the topographic control in the surface layer, altering the frontal structure, and affecting local water mass formation. In the winter, the topographic control is firmly established and both sides of the front are vertically well-mixed. Winter cooling creates sea-ice over Spitzbergen Bank and convectively formed Modified Atlantic Water in the Bear Island Trough and Hopen Trench. In the summer, heating melts the sea-ice, producing a surface meltwater pool that can cross the polar front, disrupting topographic control and substantially increasing the vertical thermohaline gradients in the frontal region. The meltwater pool produces the largest geostrophic shear in the region.

Description: Support for this work was provided by a Department of Defense National Defense Science and Engineering Graduate Fellowship and Office of Naval Research grant N00014- 90-J-1359.

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

Description: An adiabatic, inertial, and quasigeostrophic model is used to discuss the interaction of surface Ekman transport with an island. The theory extends the recent work of Spall and Pedlosky to include an analytical and nonlinear model for the interaction. The presence of an island that interrupts a uniform Ekman layer transport raises interesting questions about the resulting circulation. The consequential upwelling around the island can lead to a local intake of fluid from the geostrophic region beneath the Ekman layer or to a more complex flow around the island in which the fluid entering the Ekman layer on one portion of the island's perimeter is replaced by a flow along the island's boundary from a downwelling region located elsewhere on the island. This becomes especially pertinent when the flow is quasigeostrophic and adiabatic. The oncoming geostrophic flow that balances the offshore Ekman flux is largely diverted around the island, and the Ekman flux is fed by a transfer of fluid from the western to the eastern side of the island. As opposed to the linear, dissipative model described earlier, this transfer takes place even in the absence of a topographic skirt around the island. The principal effect of topography in the inertial model is to introduce an asymmetry between the circulation on the northern and southern sides of the island. The quasigeostrophic model allows a simple solution to the model problem with topography and yet the resulting three-dimensional circulation is surprisingly complex with streamlines connecting each side of the island.

Description: This research was supported in part by NSF Grant OCE Grant 0925061.

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

Description: Ocean carbon uptake and storage simulated by the Community Earth System Model, version 1–Biogeochemistry [CESM1(BGC)], is described and compared to observations. Fully coupled and ocean-ice configurations are examined; both capture many aspects of the spatial structure and seasonality of surface carbon fields. Nearly ubiquitous negative biases in surface alkalinity result from the prescribed carbonate dissolution profile. The modeled sea–air CO2 fluxes match observationally based estimates over much of the ocean; significant deviations appear in the Southern Ocean. Surface ocean pCO2 is biased high in the subantarctic and low in the sea ice zone. Formation of the water masses dominating anthropogenic CO2 (Cant) uptake in the Southern Hemisphere is weak in the model, leading to significant negative biases in Cant and chlorofluorocarbon (CFC) storage at intermediate depths. Column inventories of Cant appear too high, by contrast, in the North Atlantic. In spite of the positive bias, this marks an improvement over prior versions of the model, which underestimated North Atlantic uptake. The change in behavior is attributable to a new parameterization of density-driven overflows. CESM1(BGC) provides a relatively robust representation of the ocean–carbon cycle response to climate variability. Statistical metrics of modeled interannual variability in sea–air CO2 fluxes compare reasonably well to observationally based estimates. The carbon cycle response to key modes of climate variability is basically similar in the coupled and forced ocean-ice models; however, the two differ in regional detail and in the strength of teleconnections.

Description: The CESM project is supported by the National Science Foundation and the Office of Science (BER) of the U.S. Department of Energy. SCD acknowledges support of Collaborative Research: Improved Regional and Decadal Predictions of the Carbon Cycle (NSFAGS- 1048827).

Description: 2014-03-15

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

Description: Aspects of the circulation and convection in the Irminger Sea are investigated using a variety of in-situ, satellite, and atmospheric reanalysis products. Westerly Greenland tip jet events are intense, small-scale wind phenomena located east of Cape Farewell, and are important to circulation and convection in the Irminger Sea. A climatology of such events was used to investigate their evolution and mechanism of generation. The air parcels constituting the tip jet are shown to have a continental origin, and to exhibit a characteristic deflection and acceleration around southern Greenland. The events are almost invariably accompanied both by a notable coherence of the lower-level tip jet with an overlying upper-level jet stream, and by a surface cyclone located in the lee (east) of Greenland. It is argued that the tip jet arises from the interplay of the synopticscale flow evolution and the perturbing effects of Greenland’s topography upon the flow. The IrmingerGyre is a narrow, cyclonic recirculation confined to the southwest Irminger Sea. While the gyre’s existence has been previously documented, relatively little is known about its specific features or variability. The mean strength of the gyre’s circulation between 1991 and 2007 was 6.8 ± 1.8 Sv. It intensified at a rate of 4.3 Sv per decade over the observed period despite declining atmospheric forcing. Examination of the temporal evolution of the LSW layer thickness across the Irminger Basin suggests that local convection formed LSW during the early 1990s within the Irminger Gyre. In contrast, LSW appeared outside of the gyre in the eastern part of the Irminger Sea with a time lag of 2-3 years, consistent with transit from a remote source in the Labrador Sea. In the winter of 2007-08 deep convection returned to both the Labrador and Irminger seas following years of shallow overturning. The transition to a convective state took place abruptly, without going through a preconditioning phase, which is contrary to general expectations. Changes in the hemispheric air temperature, tracks of storms, flux of freshwater to the Labrador Sea, and distribution of pack ice all conspired to enhance the air-sea heat flux, resulting in the deep overturning.

Description: Financial support for this work was provided by National Science Foundation grant OCE-0450658.

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 2001

Description: As part of Phase III of the U.S. GLOBEC Georges Bank program, drogued drifters and dye tracer were released into the pycnocline on the southern flank of Georges Bank to measure advective and diffusive transport relative to the tidal mixing front in May 1999. Potential density measurements placed the tidal mixing front around the 50-55 m isobath on the southern flank. Drogued drifter movement relative to the front was on the order of the drifter's slip velocity and therefore did not support the existence of a mean, advective flow. No movement relative to the front of the dye patch center of mass also indicated a lack of advective flow. Diffusive transport did occur as the dye patch spread laterally both toward and away from the front much as would be predicted by the diffusion relationship of Okubo (1971), who summarized diffusion experiments in the surface ocean. The dye did not spread symmetrically, but was rather elongated along the isobaths. This can be attributed to vertical shear in the along-isobath current that was measured by the shipboard ADCP.

Description: The National Science Foundation (Grant #OCE9806498) sponsored the data collection described within. The United States Navy funded the author's participation in this thesis project.

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

Description: This thesis studies the role of cross-isopycnal mixing in general circulation dynamics, from both the theoretical and observational points of view. The first two chapters discuss some theoretical aspects of cross-isopycnal mixing in the oceans. In chapter one, an integral constraint relating the interior stratification and air-sea heat fluxes is derived, based on the condition that the total mass of water of given density is constant in a steady state ocean. Two simple models are then used to examine the way the numerically small mixing, together with air-sea fluxes, determines the average vertical density stratification of the oceans, and the deep buoyancy driven circulation. In chapter two, a more complete model of a deep flow driven by cross isopycnal diffusion is presented, motivated by the Mediterranean outflow into the North Atlantic. Mixing in this model is responsible for the determination of the detailed structure of the flow and density field, while in the models of the first chapter it was allowed to determine only the average vertical density stratification. In chapter three, a hydrographic data set from the Mediterranean sea is analyzed by inverse methods. The purpose is to examine the importance of mixing when trying to explain tracer distributions in the ocean. The time-mean circulation and the appropriate mixing coefficients are calculated from the hydrographic data. We conclude that the numerically small cross isopycnal mixing processes are crucial to the dynamics, yet difficult to parameterize and measure using available hydrographic data.

Description: NSF grants OCE-8521685 and OCE-8017791 supported me during my studies in the joint program.

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

Description: The various distributions of tracer associated with the Northern Recirculation Gyre of the Gulf Stream (NRG) are studied to try to obtain information about the flow. An advective-diffusive numerical model is implemented to aid in the investigation. The model is composed of a gyre adjacent to a boundary current in which a source of tracer is specified at the upstream edge of the current. This set up attempts to simulate the lateral transfer of properties from the Deep Western Boundary Current (DWBC) to the NRG in the region where the two flows are in close contact west of the Grand Banks. The results of the model are analyzed in some detail. Tracer is entrained into the gyre as a plume which extends from the boundary current and spirals across streamlines toward the gyre center. The maintenance of the spiral during spin-up and its relationship to the occurrence of homogenization at steady state is examined. An asymmetry in the spiral exists due to the ellipticity of the gyre, which also effects homogenization. The anomalous properties that are fluxed into the NRG include salt, oxygen, and freon. These particular tracers are independent from each other, the former two because they are characterized by different vertical profiles in the deep layer. This results in a decay of oxygen but not salt, due to the presence of vertical mixing as discussed by Hogg et al. (1986, Deep-Sea Research, 33, 1139-1165). Their analysis is expanded upon here. The effect of vertical mixing on the gyre/boundary current system is examined within the context of the numerical model. Results are applied to recently collected water sample data from the region which leads to an estimate of the lateral and vertical eddy diffusion coefficients and an estimate of the amount of oxygen in the NRG that has diffused from the DWBC. The accumulation of freon within the NRG is considered in addition to salt and oxygen. Appreciable levels of freon have been present in the ocean only since 1950, and the atmospheric source functions have been increasing steadily since then. A simple overflow model is presented of the manner in which freon may be stirred in the Norwegian-Greenland basin prior to overflowing and entering the DWBC. Once in the boundary current the concentrations are diluted by way of mixing with surrounding water. Two different schemes are considered in which the immediate surrounding water accumulates a substantial amount of freon as time progresses. These models suggest that the freon-11:freon-12 ratio may not be a conserved quantity for the water in the core of the DWBC. It is found that the level of freon in the NRG is barely above the existing background level.

Description: This work was supported by the office of Naval Research through contracts N00014-76-C-0197 and N00014-84-C-0134, NR 083-400; and N00014-82-C-0019 and N00014-85-C-0001, NR 083-004, and the National Science Foundation through grant 0CE82-14925.

Description: Author Posting. © American Meteorological Society, 2010. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 23 (2010): 3146–3154, doi:10.1175/2010JCLI3310.1.

Description: The Atlantic meridional overturning circulation (AMOC) simulated in various ocean-only and coupled atmosphere–ocean numerical models often varies in time because of either forced or internal variability. The path of the Gulf Stream (GS) is one diagnostic variable that seems to be sensitive to the amplitude of the AMOC, yet previous modeling studies show a diametrically opposed relationship between the two variables. In this note this issue is revisited, bringing together ocean observations and comparisons with the GFDL Climate Model version 2.1 (CM2.1), both of which suggest a more southerly (northerly) GS path when the AMOC is relatively strong (weak). Also shown are some examples of possible diagnostics to compare various models and observations on the relationship between shifts in GS path and changes in AMOC strength in future studies.

Description: We wish to acknowledge support (TJ) from WHOI’s Paul Fye Chair and NASA (NNXZX09AF35G) and to NOAA/OAR (RZ) for this work.

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

Description: Middepth, time-mean circulation in the western North Pacific Ocean (28°–45°N, 140°–165°E) is investigated using drift information from the profiling floats deployed in the Kuroshio Extension System Study (KESS) and the International Argo programs. A well-defined, cyclonic recirculation gyre (RG) is found to exist north of the Kuroshio Extension jet, confined zonally between the Japan Trench (145°E) and the Shatsky Rise (156°E), and bordered to the north by the subarctic boundary along 40°N. This northern RG, which is simulated favorably in the eddy-resolving OGCM for the Earth Simulator (OFES) hindcast run model, has a maximum volume transport at 26.4 Sv across 159°E and its presence persists on the interannual and longer time scales. An examination of the time-mean x-momentum balance from the OFES hindcast run output reveals that horizontal convergence of Reynolds stresses works to accelerate both the eastward-flowing Kuroshio Extension jet and a westward mean flow north of the meandering jet. The fact that the northern RG is eddy driven is further confirmed by examining the turbulent Sverdrup balance, in which convergent eddy potential vorticity fluxes are found to induce the cyclonic RG across the background potential vorticity gradient field. For the strength of the simulated northern RG, the authors find the eddy dissipation effect to be important as well.

Description: This study was supported by NSF through Grant OCE-0220680 (UH) and OCE-0220161 (WHOI).

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

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

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

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

Description: Labrador Sea Water (LSW), a dense water mass formed by convection in the subpolar North Atlantic, is an important constituent of the meridional overturning circulation. Understanding how the water mass enters the deep western boundary current (DWBC), one of the primary pathways by which it exits the subpolar gyre, can shed light on the continuity between climate conditions in the formation region and their downstream signal. Using the trajectories of (profiling) autonomous Lagrangian circulation explorer [(P)ALACE] floats, operating between 1996 and 2002, three processes are evaluated for their role in the entry of Labrador Sea Water in the DWBC: 1) LSW is formed directly in the DWBC, 2) eddies flux LSW laterally from the interior Labrador Sea to the DWBC, and 3) a horizontally divergent mean flow advects LSW from the interior to the DWBC. A comparison of the heat flux associated with each of these three mechanisms suggests that all three contribute to the transformation of the boundary current as it transits the Labrador Sea. The formation of LSW directly in the DWBC and the eddy heat flux between the interior Labrador Sea and the DWBC may play leading roles in setting the interannual variability of the exported water mass.

Description: We are also grateful to the NSF for their support of this research.

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

Description: A theory for the exchange between a rotating, buoyancy-forced marginal sea and an ocean is developed and tested numerically. Cooling over the marginal sea leads to sinking and sets up a two-layer exchange flow, with a warm surface layer entering from the ocean and a cool layer exiting at depth. The connecting strait is sufficiently narrow and shallow to cause the exchange flow to be hydraulically controlled. The incoming surface layer forms a baroclinically unstable boundary current that circles the marginal sea in a cyclonic sense and feeds heat to the interior by way of eddies. Consistent with the overall heat and volume balances for the marginal sea, there is a continuous family of hydraulically controlled states with critical flow at the most constricted section of the strait. Included in this family is a limiting “maximal-exchange” solution with two sections of hydraulic control in the strait and with fixed layer depths at the most constricted section. The state of exchange for a given forcing is predicted using a theory that assumes energy conservation over a certain path connecting the strait to the marginal sea or, in some cases, the ocean. Depending on the configuration of the exchange, long-wave information may be blocked from entering the strait from the marginal sea, from the open ocean, or both. The scenario that holds determines what is predicted and what needs to be input. Numerical tests of the prediction for the temperature difference and the state of exchange are carried out for straits with a pure contraction in width and for a constant width strait with a topographic sill. The comparison is reasonable in most cases, though the numerical model is not able to reproduce cases of multiple states predicted by the theory for certain forcing values. The analytical model is an alternative to the Price and Yang and Siddall et al. models of a marginal sea outflow.

Description: This work was supported by the National Science Foundation under Grants OCE-0525729 and OCE-0423975.

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

Description: This is the first of a two-part investigation of a coupled wind and wave model that includes the enhanced form drag of breaking waves. In Part I here the model is developed and applied to mature seas. Part II explores the solutions in a wide range of wind and wave conditions, including growing seas. Breaking and nonbreaking waves induce air-side fluxes of momentum and energy above the air–sea interface. By balancing air-side momentum and energy and by conserving wave energy, coupled nonlinear advance–delay differential equations are derived, which govern simultaneously the wave and wind field. The system of equations is closed by introducing a relation between the wave height spectrum and wave dissipation due to breaking. The wave dissipation is proportional to nonlinear wave interactions, if the wave curvature spectrum is below the “threshold saturation level.” Above this threshold the wave dissipation rapidly increases so that the wave height spectrum is limited. The coupled model is applied to mature wind-driven seas for which the wind forcing only occurs in the equilibrium range away from the spectral peak. Modeled wave height curvature spectra as functions of wavenumber k are consistent with observations and transition from k1/2 at low wavenumbers to k0 at high wavenumbers. Breaking waves affect only weakly the wave height spectrum. Furthermore, the wind input to waves is dominated by nonbreaking waves closer to the spectral peak. Shorter breaking waves, however, can support a significant fraction, which increases with wind speed, of the total air–sea momentum flux.

Description: This work was supported by the U.S. National Science Foundation (Grant OCE-0526177) and the U.S. Office of Naval Research (Grant N00014-06-10729).

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

Description: Techniques from dynamical systems theory have been applied to study horizontal stirring of fluid in the Philippine Archipelago. The authors’ analysis is based on velocity fields produced by two high-resolution (3 and 6 km) numerical models. Particular attention is paid to identifying robust surface flow patterns and associating them with dominant Lagrangian coherent structures (LCSs). A recurrent wind-driven dipole in the lee of the coastline is considered in detail. The associated LCSs form a template for stirring, exchange, and biological transport in and around the dipole. Chaotic advection is argued to provide a relevant framework for interpreting mesoscale horizontal stirring processes in an archipelago as a whole. Implications for the formation of filaments, the production of tracer variance, and the scale at which stirring leads to mixing are discussed in connection with an observed temperature record.

Description: This work was supported by the Office of Naval Research: Grant N00014-07-1-0590 to Woods Hole Oceanographic Institution; Grant ONR- 0601153N to Stevens Institute of Technology; Grant N00014-07-1-0417 to Institute of Marine and Coastal Sciences, Rutgers University; and Grant N00014-09-1- 0582 to Lamont-Doherty Earth Observatory of Columbia University.

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

Description: Sections of temperature, salinity, dissolved oxygen, and velocity were made crossing the Gulf Stream in late January 2006 to investigate the role of frontal processes in the formation of Eighteen Degree Water (EDW), the Subtropical Mode Water of the North Atlantic. The sections were nominally perpendicular to the stream and measured in a Lagrangian frame by following a floating spar buoy drifting in the Gulf Stream’s warm core. During the survey, EDW was isolated from the mixed layer by the stratified seasonal pycnocline, suggesting that EDW was not yet actively being formed at this time in the season and at the longitudes over which the survey was conducted (64°–70°W). However, in two of the sections, the seasonal pycnocline in the core of the Gulf Stream was broken by an intrusion of cold, fresh, weakly stratified water, nearly saturated in oxygen, that appears to have been subducted from the surface mixed layer north of the stream. The intrusion was identified in three of the sections in profiles with a nearly identical temperature–salinity relation. From the western-to-easternmost sections, where the intrusion was observed, the depth of the intrusion’s salinity minimum descended by 90 m in the 71 h it took to complete this part of the survey. This apparent subduction occurred primarily on the upstream side of a meander trough, where the cross-stream velocity was confluent and frontogenetic. Using a variant of the omega equation, the vertical velocity driven by the confluent flow was inferred and yielded downwelling in the vicinity of the intrusion spanning 10–40 m day−1, a range of values consistent with the intrusion’s observed descent, suggesting that frontal subduction was responsible for the formation of the intrusion. In the easternmost section located downstream of the meander trough, the flow was diffluent, driving an inferred vertical circulation that was of the opposite sense to that in the section upstream of the trough. In transiting the two sides of the trough, the intrusion was observed to move toward the center of the stream between the downwelling branches of the opposing vertical circulations, resulting in a downward Lagrangian mean vertical velocity and net subduction. Hydrographic evidence of the subduction of weakly stratified surface waters was seen in the southern flank of the Gulf Stream as well. The solution of the omega equation suggests that this subduction was associated with a relatively shallow vertical circulation confined to the upper 200 m of the water column in the proximity of the front marking the southern edge of the warm core.

Description: Support came from the National Science Foundation Grant OCE-0424865 (TJ and LT). LT was also supported by NSF Grant OCE- 0549699 and the Worzel Assistant Scientist Endowed Fund at WHOI.

Description: Author Posting. © American Meteorological Society, 2010. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 23 (2010): 4841–4855, doi:10.1175/2010JCLI3273.1.

Description: A 1-Myr-long time-dependent solution of a zonally averaged ocean–atmosphere model subject to Milankovitch forcing is examined to gain insight into long-term changes in the planetary-scale meridional moisture flux in the atmosphere. The model components are a one-dimensional (latitudinal) atmospheric energy balance model with an active hydrological cycle and an ocean circulation model representing four basins (Atlantic, Indian, Pacific, and Southern Oceans). This study finds that the inclusion of an active hydrological cycle does not significantly modify the responses of annual-mean air and ocean temperatures to Milankovitch forcing found in previous integrations with a fixed hydrological cycle. Likewise, the meridional overturning circulation of the North Atlantic Ocean is not significantly affected by hydrological changes. Rather, it mainly responds to precessionally driven variations of ocean temperature in subsurface layers (between 70- and 500-m depth) of this basin. On the other hand, annual and zonal means of evaporation rate and meridional flux of moisture in the atmosphere respond notably to obliquity-driven changes in the meridional gradient of annual-mean insolation. Thus, when obliquity is decreased (increased), the meridional moisture flux in the atmosphere is intensified (weakened). This hydrological response is consistent with deuterium excess records from polar ice cores, which are characterized by dominant obliquity cycles.

Description: A. A. thanks the Global Environmental and Climate Change Centre of McGill University for a Network Grant that made possible an enriching twoweek stay at WHOI during June 2007. O. M. acknowledges support from theU.S.National Science Foundation. Support from a Canadian NSERC Discovery Grant awarded to L.A.M. is gratefully acknowledged.

Description: Author Posting. © American Meteorological Society, 2009. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 22 (2009): 4066–4082, doi:10.1175/2009JCLI2629.1.

Description: Small-scale variation in wind stress due to ocean–atmosphere interaction within the atmospheric boundary layer alters the temporal and spatial scale of Ekman pumping driving the double-gyre circulation of the ocean. A high-resolution quasigeostrophic (QG) ocean model, coupled to a dynamic atmospheric mixed layer, is used to demonstrate that, despite the small spatial scale of the Ekman-pumping anomalies, this phenomenon significantly modifies the large-scale ocean circulation. The primary effect is to decrease the strength of the nonlinear component of the gyre circulation by approximately 30%–40%. This result is due to the highest transient Ekman-pumping anomalies destabilizing the flow in a dynamically sensitive region close to the western boundary current separation. The instability of the jet produces a flux of potential vorticity between the two gyres that acts to weaken both gyres.

Description: AH and WD were supported by an ARC Linkage International Grant (LX0668781). WD was also supported by NSF Grants OCE 0424227 and OCE 0550139. Funding for PB was provided by NSF Grants OCE 0344094 and OCE 0725796 and by the research grant from the Newton Trust of the University of Cambridge. SK was supported by U.S. DOE Grant DE-FG02–02ER63413 and NASA Grant NNG-06- AG66G-1.

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

Description: Motivated by observations of surface drifters in the Adriatic Sea, transport in a three-gyre system is studied with the aid of dynamical systems techniques. Particular attention is paid to the issue of intergyre transport. The velocity field is assumed to be two-dimensional and incompressible and composed of a steady three-gyre background flow on which a time-dependent perturbation is superimposed. Two systems of this type are considered: 1) an observationally motivated, analytically prescribed model consisting of a steady background on which a multiperiodic time-dependent perturbation is superimposed, and 2) an observationally based model of the Adriatic Sea consisting of the mean surface circulation derived from surface drifter trajectories on which a time-dependent altimetry-based perturbation velocity field is superimposed. It is shown that for a small perturbation to the steady three-gyre background, two of the gyres exchange no fluid with the third gyre. When the perturbation strength exceeds a certain threshold, transport between all three gyres occurs. This behavior is described theoretically, illustrated using the analytic model and shown to be consistent with the observationally based model of the Adriatic. The relevance of the work presented to more complicated multiple-gyre problems is discussed.

Description: This work was supported by the National Science Foundation Grants CMG0417425 and CMG0825547.

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

Description: Along the Taiwan Strait (〈100 m in depth) a northeastward flow persists in all seasons despite the annually averaged wind stress that is strongly southwestward. The forcing mechanism of this countercurrent is examined by using a simple ocean model. The results from a suite of experiments demonstrate that it is the Kuroshio that plays the deciding role for setting the flow direction along the Taiwan Strait. The momentum balance along the strait is mainly between the wind stress, friction, and pressure gradient. Since both wind stress and friction act against the northward flow, it is most likely the pressure gradient that forces the northward flow, as noted in some previous studies. What remains unknown is why there is a considerable pressure difference between the southern and northern strait. The Kuroshio flows along the east coast of Taiwan, and thus the western boundary current layer dynamics applies there. Integrating the momentum equation along Taiwan’s east coast shows that there must be a pressure difference between the southern and the northern tip of Taiwan to counter a considerable friction exerted by the mighty Kuroshio. This same pressure difference is also felt on the other side of the island where it forces the northward flow through Taiwan Strait. The model shows that the local wind stress acts to dampen this northward flow. This mechanism can be illustrated by an integral constraint for flow around an island.

Description: This study has been supported by National Science Foundation through Grant OCE-0351055.

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

Description: Efforts to understand the Arctic system have recently focused on the role in local and global circulation of waters from the Arctic shelf seas. In this study, steady-state exchanges between the Arctic shelves and the central basins are estimated using an inverse box model. The model accounts for data uncertainty in the estimates, and quantifies the solution uncertainty. Other features include resolution of the two-basin Arctic hydrographic structure two-way shelf-basin exchange in the surface mixed layer, the capacity for shelfbreak upwelling, and recognition that most inflows enter the Arctic via the shelves. Aggregate estimates of all fluxes across the Arctic boundary, with their uncertainties, are generated from flux estimates published between 1975 and 1997. From the aggregate estimates, mass-, heat-, and salt-conserving boundary flux estimates are derived, which imply a net flux of water from the shelves to the basins of 1.2±0.4 Sv. Due primarily to boundary flux data uncertainty, constraints of mass, heat, and salt conservation alone cannot determine how much shelf-basin exchange occurs via dense overflows, and how much via the surface mixed layer. Adding δ180 constraints, however, greatly reduces the uncertainty. Dense water flux from the shelves to the basins is necessary for maintaining steady state, but shelfbreak upwelling is not required. Proper representation of external sources feeding the shelves, rather than the basins, is important to obtain the full range of plausible steady solutions. Implications of the results for the study of Arctic change are discussed.

Description: This work was supported by National Science Foundation grant OPP-9422292 as part of the Arctic System Science ARCSS program, administered by the Office of Polar Programs.

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

Description: Trace metal cycling is one of many processes that influence ocean ecosystem dynamics. Cobalt, iron, and manganese are redox active trace metal micronutrients with oceanic distributions that are influenced by both biological and abiotic sources and sinks. Their open ocean concentrations range from picomolar to nanomolar, and their bioavailabilities can impact primary production. Understanding the biogeochemical cycling of these hybrid-type metals with an emphasis on cobalt was the focus of this thesis. This was accomplished by determining the dissolved distributions of these metals in oceanic regions that were characterized by different dominant biogeochemistries. A large subsurface plume of dissolved cobalt, iron, and manganese was found in the Eastern South Atlantic. The cause of this plume is a combination of reductive dissolution in coastal sediments, wind-driven upwelling, advection, biological uptake, and remineralization. Additional processes that are discussed as sources of metals to the regions studied during this thesis include isopycnal uplift within cold-core eddies (Hawaii), ice melt (McMurdo Sound, Antarctica), riverine input (Arctic Ocean), and winter mixing (McMurdo Sound). The biological influence on surface ocean distributions of cobalt was apparent by the observation of linear relationships between cobalt and phosphate in mid to low latitudes. The cobalt:phosphate ratios derived from these correlations changed over orders of magnitude, revealing dynamic variability in the utilization, demand, and sources of this micronutrient. Speciation studies suggest that there may be two classes of cobalt binding ligands, and that organic complexation plays an important role in preventing scavenging of cobalt in the ocean. These datasets provided a basis for comparing the biogeochemical cycles of cobalt, iron, and manganese in three oceanic regimes (Hawaii, South Atlantic, McMurdo Sound). The relative rates of scavenging for these metals show environmental variability: in the South Atlantic, cobalt, iron, and manganese were scavenged at very different rates, but in the Ross Sea, mixing and circulation over the shallow sea was fast, scavenging played a minor role, and the cycles of all three metals were coupled. Studying the distributions of these metals in biogeochemically distinct regions is a step toward a better understanding of their oceanic cycles.

Description: Funding for this research was provided by the the National Science Foundation Chemical Oceanography (Division of Ocean Sciences OCE-0452883, OCE-0752291, OCE-0928414, OCE-0732665, OCE-0440840, OCE-0327225), the Center for Microbial Research and Education, the WHOI Coastal Ocean Institute, and the WHOI Ocean Life Institute, WHOI Academic Programs Office, and a Fye Teaching Fellowship.

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

Description: The relationship between depth-averaged velocity and bottom stress for wind-driven flow in unstratified coastal waters is examined here. The adequacy of traditional linear and quadratic drag laws is addressed by comparison with a 2 1/2-D model. A 2 1/2-D model is one in which a simplified 1-D depth-resolving model (DRM) is used to provide an estimate of the relationship between the flow and bottom stress at each grid point of a depth-averaged model (DAM). Bottom stress information is passed from the DRM to the DAM in the form of drag tensor with two components: one which scales the flow and one which rotates it. This eliminates the problem of traditional drag laws requiring the flow and bottom stress to be collinear. In addition , the drag tensor field is updated periodically so that the relationship between the velocity and bottom stress can be time-dependent. However, simplifications in the 2 1/2-D model that render it computationally efficient also impose restrictions on the time-scale of resolvable processes. Basically, they must be much longer than the vertical diffusion time scale. Four progressively more complicated scenarios are investigated. The important factors governing the importance of bottom friction in each are found to be 1) non-dimensional surface Ekman depth, u.5/fh where u.s is the surface shear velocity, f is the Coriolis parameter and h is the water depth 2) the non-dimensional bottom roughness, zo/h where zo is the roughness length and 3) the angle between the wind stress and the shoreline. Each has significant influence on the drag law. The drag tensor magnitude, r, and the drag sensor angle, θ are functions of all three, while a drag tensor which scales with the square of the depth-averaged velocity has a magnitude, Cd, that only depends on zo/h. The choice of drag Jaw is found to significantly affect the response of a domain. Spin up times and phase relationships vary between models. In general, the 2 1/2-D model responds more quickly than either a constant r or constant Cd model. Steady-state responses are also affected. The two most significant results are that failure to account for θ in the drag law sometimes leads to substantial errors in estimating the sea surface height and to extremely poor resolution of cross-shore bottom stress. The latter implies that cross-shore near-bottom transport is essentially neglected by traditional DAMs.

Description: Financial support during my time in graduate school came from the Woods Hole Oceanographic Institution and grants from the National Science Foundation (OCE84-03249) and the Office of Naval Research (N00014-86-K-0061).

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: This thesis studies mixing and convection in a rectangular basin driven by a specified heat flux at the surface. A numerical model is constructed for this purpose. The main focus of the study is on the density and circulation structure resulting from the thermal forcing. In chapter two, a simple vertical one-dimensional model is developed to examine the mixing processes under a given surface heat flux. In order to simulate strong vertical mixing in the region where stratification is unstable, turbulent processes are modeled by a convective overturning parameterization of eddy viscosity and diffusivity. The results show that the density structure is strongly affected by the convective overturning adjustment as surface cooling prevails, and the resulting density field is nearly depth independent. In chapter three, a more complicated two-dimensional model is constructed to simulate mixing and circulation in a vertical rectangular basin with rigid boundaries. The aspect ratio of the basin ranges from 1 to 0.001 and Rayleigh number from 104 to 2 x 1012. It is found that the circulation pattern is dominated by these two important numbers. The roles of density overturning and density-momentum overturning mixing are further investigated. The results show that the convective overturning not only homogenizes the density field in the unstably stratified region but also contributes to increase the circulation. A crude scale analysis of the system shows that the characteristics of the density and momentum fields from the analysis agree well with the numerical results.

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 1993

Description: An extensive set of new high-quality hydrographic data is assembled in order to determine the mean circulation in the equatorial Pacific, and thus the pathways for cross-equatorial and cross-gyre exchange. Making up the core of the data set are two onetime transpacific zonal sections nominally at 10°N and 14°S. Supplementing these are repeat surveys of the equatorial currents along the 165°E meridian with direct shear measurements, and repeat surveys of the western boundary current at 8°N including direct velocity measurements. The repeat survey data are crucial for obtaining a good estimate of the mean conditions in the face of strong annual and interannual variability of the near-equatorial flow field. A comparison with historical XBT and hydrographic data shows that the interior thermocline transports in the one-time sections are fortuitously representative of the mean conditions. A detailed study of the water mass distribution along the sections is the basis for choosing reference levels for the thermal wind shear in an initial guess circulation field. Using an inverse model, the initial guess circulation is adjusted such that volume, heat and salt arc conserved in a set of subthermocline layers (δΘ 〉 26.7). Cross-isopycnal diffusion and advection are explicitly accounted for in the inverse model, and the diapycnal diffusivity is constrained to be positive, though its value is allowed to vary with depth and location. Net mass conservation constraints are applied to the enclosed volumes of the North Pacific and eastern Pacific, and essentially require that the Ekman divergence be equal to the geostrophic convergence. The Ekman fluxes as estimated from wind-stress climatologies are an important element of the mass budget, and yet are subject to large uncertainties. The model is therefore given the freedom to determine the Ekman fluxes within the range of error of the wind-stresses. The circulation of the coldest waters (Θ 〈 1.2°C) is dominated by the northward flow of Lower Circumpolar Water (LCPW) in a system of narrow western boundary currents. A net transport of 12.1 Sv of LCPW flows across 14°S, 9.6 Sv of which flows into the North Pacific across 10°N. The bulk of the LCPW flux across the equator appears to occur in the denser part of the western boundary current which follows topography directly across the equator. Dissipation in the boundary layer can thus modify the potential vorticity of the fluid and allow it to cross the equator. The circulation of the upper part of the LCPW is dominated by a strong westward jet at the equator which is supplied both by upwelling from below and the recirculation of modified LCPW from the North Pacific. At mid-depth (4.0 〉 Θ 〉 1.2°C) high silica and low oxygen concentrations mark the North Pacific Deep Water (NPDW) which is present in both the North and South Pacific Oceans. Across both 10°N and 14°S, a net of 11 Sv of NPDW flows southward, returning the northward mass flux associated with the LCPW. In contrast to the LCPW, narrow western boundary currents are not present in this layer, and it is not clear how the deep water flows across the equator. Strong zonal jets on and about the equator may be important in allowing mass to cross the equator by increasing the time available for the cross-equatorial diffusion of potential vorticity to act on a fluid parcel. At intermediate depths equatorward advection is suggested by the presence of intermediate water salinity minima formed in the subpolar latitudes: Antarctic Intermediate Water dominates the 4 to 8°C classes south of the equator, while North Pacific Intermediate Water occupies this range north of the equator. Determination of the mean circulation of the intermediate waters is, however, confounded by the large eddies that dominate the geostrophic transport stream function along the onetime zonal sections. The equatorial thermocline is occupied by waters of subtropical origin: the shallow salinity minimum waters and saline Central Water from both the North and South Pacific Ocean. The equator marks the location of a front between northern and southern subtropical gyre waters, except in the lower thermocline where water from the South Pacific subtropical gyre penetrates to about 4°N to feed the Northern Subsurface Countercurrent at 165°E. All of the equatorward flowing thermocline waters are entrained in the eastward equatorial currents which in turn feed the upwelling system in the eastern Pacific. The upwelled waters largely supply the South Equatorial Current in the eastern Pacific, accounting for its large transport compared to that predicted by Sverdrup dynamics. Northward flow across the equator of the upwelled waters in the thermocline or surface layer in the western Pacific is necessary to supply the Ekman flux into the North Pacific. The analysis indicates that the Pacific Ocean does not convert a large amount of abyssal water to thermocline water, as required by several theories of the global thermohaline circulation. In contrast to the Atlantic Ocean, the thermocline circulation in the Pacific appears decoupled from the abyssal overturning, with little upwelling of abyssal waters occurring in either the North Pacific or the equatorial Pacific. The leakage of Pacific water into the Indian Ocean is deduced to be essentially zero, though an error analysis allows a range of 0-8 x 106m3s-1.

Description: I was supported by the 1986 Caltex Graduate Women Scholarship, and a NASA Scholarship in Global Change Research.

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 1994

Description: Paleo-tracers such as carbon 13 and cadmium show that the deep Atlantic was enriched in nutrients during the Last Ice Age. The conventionally accepted interpretation of these higher nutrient levels is that a reduction of the rate of formation of nutrient-depleted Lower North Atlantic Deep Water (Lower NADW) allowed nutrient-rich Antarctic Bottom Water (AABW) to push further north during the Last Glacial Maximum (LGM) (Boyle and Keigwin, 1982; 1987; Duplessy et al., 1988). The evidence for this interpretation is re-examined in this work, with an emphasis on the quantitative analysis of the paleo-data. An end-member analysis of the δ13C data indicates a larger volume of AABW and a smaller volume of Lower NADW during the LGM. It is not yet possible, however, to quantify the extent of the volume differences between the modern and the glacial distributions, because the LGM δ13C end-members are poorly known. The second issue examined in this thesis deals with the interpretation of the water mass distribution, inferred from paleo-tracers, in terms of the oceanic circulation. Using a dynamical inverse model of the North Atlantic and a kinematic inverse model of the South Atlantic, it is shown that a tracer distribution corresponding to a significantly reduced volume of Lower NADW does not necessarily correspond to a reduced flux of NADW. Indeed, a circulation almost identical to a modem ocean reference circulation is consistent with the available LGM δ13C and δ18O data A flux of Lower NADW reduced by 50%, though not needed to explain the LGM tracer distribution, is also consistent with the data Thus, the paleo-tracers δ13C and δ18O do not suffice to quantify the flux of NADW in the glacial ocean. The modem ocean circulation is one of many possible circulations consistent with the available δ13C and δ18O data.

Description: This research was funded by the National Science Foundation under grant OCE-9205942.

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 October 1993

Description: In this thesis, production of dense water that feeds the dense overflows across the Greenland-Scotland Ridge has been considered. A new circulation scheme is developed which is consistent with the water masses, currents and air-sea fluxes in the region, and with the important observation that the dense overflows show little or no seasonal or interannual variability. An inverse box model has been constructed that shows that the new circulation scheme is consistent with conservation statements for mass, heat and salt as well. According to the new circulation scheme the major buoyancy is lost in the North Atlantic Current, which enters the Norwegian Sea between Iceland and Scotland, and flows northward towards the Arctic Ocean and the Barents Sea. The transformation is due to a large net annual heat loss over the North Atlantic Current, combined with a long residence time (2-3 years) and a large surface area. After subduction, one branch of the North Atlantic Current enters the Arctic Ocean, is modified in hydrographic properties into those associated with the Denmark Strait Overflow Waters in the western North Atlantic, exits the Arctic Ocean in the western Fram Strait and flows with the East Greenland Current towards the Denmark Strait Another branch of the North Atlantic Current recirculates directly in the Fram Strait and flows towards the Denmark Strait with the East Greenland Current This branch will not sink to the bottom of the North Atlantic as it is less compressible than the Arctic branch. The third branch of the North Atlantic Current enters the Barents Sea, continues to lose buoyancy, and enters the Arctic Ocean at intermediate depth. This branch exits the Arctic Ocean in the western Fram Strait, circulates around the Greenland Sea, enters the Norwegian Sea, and flows towards the Frer¢-Shetland Channel. The traditional view holds that the major sources of the dense overflows are the Iceland and Greenland gyres, west of the North Atlantic Current. Aside from the finding that the new circulation scheme is more likely in terms of water mass properties, currents etc., one fundamental problem with the old scheme lies with supplying a substantial overflow. There are indications that the production of dense water in the gyres is sensitive to the highly variable surface conditions and that indeed the production tends to shut on and off. The reservoirs in the gyres are so small that they would be drained within a few years if they were to supply the overflows during a shutdown period. Production of dense water within the North Atlantic Current is less sensitive to surface conditions. The density in the gyres is gained at a temperature around freezing, whereas in the North Atlantic Current the density is gained well above freezing. Therefore a freshwater anomaly in the two domains will have different consequences for vertical · overturning: within the North Atlantic Current the freshening can be overcome by further cooling, whereas in the gyres freezing will occur and the vertical overturning will cease. The observed lack of a significant seasonal signal associated with the dense overflows is consistent with the new circulations scheme. The net annual cooling dominates the seasonal oscillation in the atmospheric heat loss for time scales comparable with the residence time of the Atlantic Water within the domain. Thus winter formation of dense water within the North Atlantic Current does not induce a seasonal signal in the transport field of the dense water.

Description: Funding for this work was partly provided by a NASA Global Change Fellowship.

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 July 1993

Description: In this thesis, the dynamic role of bottom topography in a β-plane channel is systematically studied in both linear homogeneous and stratified layer models in the presence of either wind stress (Chapters 2, 3, 4, and 6) or buoyancy forcing (Chapter 5). In these studies, the structure of the geostrophic contour plays a fundamental role, and the role of bottom topography is looked at from two different angles. It is shown that blocking all the geostrophic contours leads to two different physical processes in which bottom topographic form drag is generated (Chapters 2, 3 and 4) and enables geostrophic flow in a β-plane channel to support a net cross-channel volume transport (Chapters 5 and 6). It is demonstrated that by blocking all the geostrophic contours in the presence of a sufficiently high ridge, the dynamics of both source-sink and wind driven circulations in a β-plane is similar to that in a dosed basin. First, wind-driven circulation in the inviscid limit is discussed in a linear barotropic channel model in the presence of a bottom ridge. There is a critical height of the ridge, above which all geostrophic contours in the channel are blocked. In the subcritical case, the Sverdrupian balance does not apply and there is no solution in the inviscid limit. In the supercritical case, however, the Sverdrupian balance applies. The form drag is generated through two different physical processes: the through-channel recirculating flow and the Sverdrupian gyre flow. These processes are fundamentally different from the nonlinear Rossby wave drag generation. In this linear model, the presence of a supercritical high ridge is essential in the inviscid limit. With this form drag generation determined, an explicit form for the zonal transport in the channel is obtained, which shows what model parameters determine the through-channel transport. In addition, the model demonstrates that most of the potential vorticity dissipation occurs at the northern boundary where the ridge intersects. The result from the homogeneous channel model in Chapter 2 is then extended to a model whose geometry consists of a zonal channel and two partial meridional barriers along each boundary at the same longitude. Both the model transport and especially the model circulation are significantly affected by the presence of the two meridional barriers. The presence of the northern barrier always leads to a decrease in the transport. The presence of the southern barrier, however, increases the transport for a narrow ridge. The northern barrier only has a localized influence on the circulation pattern, while the southern barrier has a global influence in the channel. Then a multi-layer Q-G model is constructed by assuming that potential vorticity in all subsurface layers is homogenized. The circulation is made up of baroclinic and the barotropic part. The barotropic part is same as that in a corresponding barotropic model, and is solely determined by the wind stress, while the baroclinic part is not directly related to the wind stress. It is determined by the potential vorticity homogenization and lateral boundary conditions. The presence of the stratification does not affect the bottom topographic form drag generation. The interfacial form drag is generated by the stationary eddies. Corresponding to the circulation structure, the zonal through-channel transport associated with the barotropic circulation is determined by the wind stress and bottom topography. The other part associated with the baroclinic circulation, however, is not directly related to the wind stress and it is determined by the background stratification. Based upon the discussion on the geostrophic contour, a simple barotropic model of abyssal circulation in a circumpolar ocean basin is constructed. The presence of a supercritically high ridge is both necessary and sufficient for geostrophic flow in a β-plane channel to support a net cross-channel volume flux. In the presence of a sufficiently high ridge, the classical Stommel & Arons theory applies here, but with significant modifications. The major novelty is that a throughchannel recirculation is generated. Both its strength and direction depend critically upon the model parameters. Then, a schematic picture of the abyssal circulation in a rather idealized Southern Ocean is obtained. The most significant feature is the narrow current along the northern boundary of the circumpolar basin, which feeds the deep western boundary currents of the Indian Ocean and Pacific Ocean and connects all the oceanic basins in the Southern Ocean. Finally, the question of how the northward surface Ekman transport out of the circumpolar ocean is returned is discussed in a two-layer model with an infinitesimally thin surface Ekman layer on top of a homogeneous layer of water in a rather idealized Southern Ocean basin. First, the case with a single subtropical ocean basin is discussed. In the case with a sufficiently high ridge connecting the Antarctic and the meridional barrier, an explicit solution is found. The surface Ekman layer sucks water from the lower layer in the circumpolar basin. This same amount of water flows northward as the surface Ekman drift. It downwells in the subtropical gyre, and is carried to the western boundary layer. From there, the same amount of water flows southward as a western boundary current across the inter-gyre boundary between the circumpolar ocean and the subtropical gyre along the west coast to the southern boundary of the meridional barrier. Then, the same amount of water is carried southward and feeds the water loss to the surface Ekman layer due to the Ekman sucking in the interior circumpolar ocean. The case with multiple subtropical ocean basins such as the Southern Ocean is also discussed. It is demonstrated that the surface Ekman drift drives a strong inter-basin water mass exchange.

Description: This thesis was supported by National Science Foundation through grant OCE OCE90-17158. Part of the numerical simulation was performed at NCAR's supercomputer, and was supported by SCD/NCAR.

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

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

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

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

Description: Efforts to monitor the ocean for signs of climate change are hampered by ever-present noise, in the form of stochastic ocean variability, and detailed knowledge of the character of this noise is necessary for estimating the significance of apparent trends. Typically, uncertainty estimates are made by a variety of ad hoc methods, often based on numerical model results or the variability of the data set being analyzed. We provide a systematic approach based on the four-dimensional frequency-wavenumber spectrum of low-frequency ocean variability. This thesis presents an empirical model of the spectrum of ocean variability for periods between about 20 days and 15 years and wavelengths of about 200{10,000 km, and describes applications to ocean circulation trend detection, observing system design, and satellite data processing. The horizontal wavenumber-frequency part of the model spectrum is based on satellite altimetry, current meter data, moored temperature records, and shipboard ADCP data. The spectrum is dominated by motions along a "nondispersive line". The observations considered are consistent with a universal ω-2 power law at the high end of the frequency range, but inconsistent with a universal wavenumber power law. The model spectrum is globally varying and accounts for changes in dominant phase speed, period, and wavelength with location. The vertical structure of the model spectrum is based on numerical model results, current meter data, and theoretical considerations. We find that the vertical structure of kinetic energy is surface intensified relative to the simplest theoretical predictions. We present a theory for the interaction of linear Rossby waves with rough topography; rough topography can explain both the observed phase speeds and vertical structure of variability. The improved description of low-frequency ocean variability presented here will serve as a useful tool for future oceanographic studies.

Description: This research was supported by NASA under grants NNG06GC28G and NNX08AR33G.

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

Description: Moored time series from the Coastal Ocean Dynamics Experiment (CODE), Shelf Mixed Layer Experiment (SMILE), Sediment Transport Events over the Shelf and Slope (STRESS) study, and Northern California Coastal Circulation Study (NCCCS) are used to study subtidal cross-shelf circulation over the northern California shelf. The northern California shelf, like much of the United States Pacific coast, is subject to strong wind forcing which exhibits characteristic seasonality. In winter and early spring, it is distinguished by poleward and equatorward fluctuations on time scales of days and by weak monthly means. In summer, it is distinguished by periods of equatorward stress lasting several weeks and by relatively strong monthly means. The intensive winter and spring SMILE and STRESS and summer CODE-2 field programs permit the examination of cross-shelf circulation under both types of wind forcing conditions at a mid-shelf site (~90 m) 6 km from the northern California coast. The primary thesis goal is to examine the applicability of a two-dimensional conceptual model of wind-forced cross-shelf circulation. In this conceptual model, surface and bottom cross-shelf flows are forced by along-shelf wind stress and bottom stress, and interior cross-shelf flow compensates such that the depth-averaged flow is zero. A secondary thesis goal is to use the seasonal coverage of available field programs to gain insight into seasonal variability of cross-shelf circulation on the northern California shelf. To accomplish these goals, the observed subtidal cross-shelf circulation is examined in the context of the winter and spring heat and salt balances, an analytic model of wind-forced cross-shelf circulation, and the spatial scales of subtidal velocity. Mean and fluctuating heat and salt balances estimated between December, 1988 and May, 1989 demonstrate the importance of cross-shelf fluxes and their general consistency with the simple conceptual model. Mean fluxes are consistent with the weak mean equatorward wind stress observed during SMILE. The dominant terms in the fluctuating balances are the cross-shelf fluxes and local changes in heat and salt content. These are well correlated with each other and with the local along-shelf wind stress. The along-shelf heat flux divergence is of secondary importance to the fluctuating heat balance. It is uncorrelated with the along-shelf wind stress, and occurrences when it is strong are interpreted as effects of mesoscale features. To examine the applicability of the wind-forced conceptual model in more detail, a simple analytic model incorporating the assumptions of the conceptual model and observed local wind forcing is compared quantitatively to estimates of surface mixed layer, interior, and bottom mixed layer cross-shelf transport for winter SMILE and STRESS and summer CODE-2 observations. This comparison suggests the model is more suited to the transient wind forcing observed during SMILE and STRESS than to the steady wind forcing observed during CODE-2. For 2-3 day wind events between December, 1988 and February, 1989, the model is well correlated with observed depthdependent (total minus depth-averaged) transports throughout the water column and with total surface mixed layer transports. For 2-3 week wind events between April and July, 1982, the model does not work nearly as well below the surface mixed layer. In the absence of other processes, the locally wind-forced model implies that the wind stress sets the horizontal scales of subtidal velocity. Correlation scales estimated for subtidal along-shelf velocity over the northern California shelf are for all field programs longer than the maximum mooring separation (60 km) and are similar to those of the wind stress. However, along-shelf correlation scales of cross-shelf velocity are shorter than minimum mooring separations for CODE. SMILE and NCCCS time series do resolve along-shelf correlation scales for near surface cross-shelf velocity. During this time, along-shelf correlation scales for near surface cross-shelf velocity vary on a monthly time scale. They are generally long (30 km or more) when correlation with wind stress is high and short (15 km or less) when correlation with wind stress is low. On at least one occasion, short along-shelf correlation scales coincide with the intrusion of an offshore mesoscale feature onto the shelf. Results of the three studies show the two-dimensional model offers some insight into the observed subtidal cross-shelf circulation, particularly in winter. During this time, the heat balance, analytical transport model, and correlation scales all provide evidence that the winter wind-forced circulation is quasi-two-dimensional. Threedimensional variability on the shelf, though important on occasion, does not appear to be wind-driven and may result from the influence of offshore mesoscale features. A quite different story emerges for summer when the simple conceptual model of crossshelf circulation fails to describe adequately subsurface cross-shelf flow. Two useful areas of further investigation may be the non-linear response of cross-shelf velocity to wind forcing and its response to other processes such as remotely generated mesoscale features.

Description: In my first year at MIT, Carl Wunsch supported me through NSF grant OCE 88- 23043. At WHOI, I first started work on SMILE with Bob Beardsley under NSF grant OCE 87-16937 and continued working on it with Steve Lentz under NSF grant OCE 91-15713.

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 1995

Description: Data from fifteen globally distributed, modern, high resolution, hydrographic oceanic transects are combined in an inverse calculation using large scale box models. The models provide estimates of the global meridional heat and freshwater budgets and are used to examine the sensitivity of the global circulation, both inter and intra-basin exchange rates, to a variety of external constraints provided by estimates of Ekman, boundary current and throughflow transports. A solution is found which is consistent with both the model physics and the global data set, despite a twenty five year time span and a lack of seasonal consistency among the data. The overall pattern of the global circulation suggested by the models is similar to that proposed in previously published local studies and regional reviews. However, significant qualitative and quantitative differences exist. These differences are due both to the model definition and to the global nature of the data set. The picture of the global circulation which emerges from the models IS a complex, turbulent flow. When integrated across ocean basins not one, but two major cells emerge. The first connects an Atlantic overturning cell (estimated at 18± 4x 109 kg s- 1) to the Southern Ocean where the Antarctic Circumpolar Current carries lower deep waters to the Indian and Pacific basins where they are converted to upper deep and intermediate waters before returning to the Atlantic. The second cell connects the Pacific and Indian Basins to the north and south of Australia. In t his cell deep waters pass into the Pacific and return within the Indian Basin as intermediate waters after passing through the Indonesian Passages. The two cells are found to be independent of one another, i.e. within the models, the Indonesian Passages do not represent a significant element in a net global circulation. While there is ample evidence of westward flow around the southern tip of South Africa which would support a "warm" water path scenario, the variability of flow in this region, rich with eddies makes hydrography a poor estimator of the relative strengths of the controversial "warm" and "cold" water paths. All existing estimates of Indonesian Passage throughflow, including the smallest (O x 106 m3 s-1) and the largest (20 x 106 m3 s-1), are consistent with the model constraints. When the Pacific- Indian throughflow is not constrained, the model produces an estimate of 11 ± 14x 109 kg s-1. The model heat flux estimates are both significantly different from zero and quite robust to changes in initial assumptions, with the exception of the choice of wind field. Although in this work it was not possible to compute freshwater fluxes which were significantly different from zero, future inclusion of salinity anomaly constraints along with terms describing vertical diffusion may yet make it possible to compute significant freshwater :flux estimates from hydrography.

Description: This research was partially funded by a NASA Global Change Fellowship and was also supported by NASA under contract NAGW-1048 and NSF under contract OCE-9205942.

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: Inverse modeling activities in oceanography have recently been intensified, aided by the oncoming observational data stream of WOCE and the advance of computer power. However, interpretations of inverse model results from climatological hydrographic data are far from simple. This thesis examines the behavior of an inverse model in the WOCE CME (Community Modeling Effort) results where the physics and the parameter values are known. The ultimate hypotheses to be tested are whether the inferred circulations from a climatological hydrographic data set (where limited time means and spatial smoothing are usually used) represent the climatological ocean general circulations, and what the inferred "diffusion" coefficients really are. The inverse model is first tested in a non-eddy resolving numerical GCM ocean. Numerical/scale analyses are used to test whether the inverse model properly represents the GCM ocean. Experiments show how biased answers could result from an incorrect model, and how a correct model must produce the right answers. When the inverse model is applied to the time-mean hydrographic data of an eddy-resolving GCM ocean in the fine grid resolution of the GCM, the estimated horizontal circulation is statistically consistent with the EGCM time means in both patterns and values. Although the flow patterns are similar, the uncertainties for the GCM time means and the inverse model estimates are different. The former are very large, such that the GCM time-mean circulation has no significance in the deep ocean. The latter are much smaller, and with them the estimated circulations are well defined. This is consistent with the concept that ocean motions are very energetic, while variations of tracers (temperature, salinity) are low frequency. The inverse model succeeded in extracting the ocean general circulation from the "climatological" hydrographic data. The estimated vertical velocities are also statistically indistinguishable from the GCM time means. However, significant differences between the estimated "diffusion" coefficients and the EGCM eddy diffusion coefficients are found at certain locations. These discrepancies are attributed to the differences in physics of the inverse model and the EGCM ocean. The "diffusion" coefficients from the inversion parameterize not only the eddy fluxes, but also (part of) the temporal variation and biharmonic terms which are not explicitly included in the inverse model. Given the essentially red spectrum of the ocean, it makes sense to look for smooth solutions. Aliasing due to subsampling on a coarse grid and the effects of spatial smoothing are addressed in the last part of this thesis. It is shown that this aliasing could be greatly reduced by spatial smoothing. The estimated horizontal circulation from the spatially smoothed time-mean EGCM hydrographic data with a coarse grid resolution (2.4° longitude by 2.0° latitude) is generally consistent with the spatially smoothed EGCM time means. Significant differences only occur at some grid points at great depths, where the GCM circulations are very weak. The conclusions of this study are different from some previous studies. These discrepancies are explained in the concluding chapter. Finally, it should be pointed out that the issue of properly representing a GCM ocean by an inverse model is not identical to the issue of represent ing the real ocean by the same inverse model, since the GCM ocean is not identical to the real ocean. Numerical calculations show that both the non-eddy resolving and the eddy-resolving GCM oceans used in this work are evolving towards a statistical equilibrium. In the real ocean, the importance of temporal variation terms in the property conservation equations should also be analyzed when a steady mverse model is applied to a limited time-mean (the climatological) data set.

Description: This research was carried out under National Science Foundation grant OCE- 90-04396.

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 2013

Description: This thesis explores the role that the circulation in the Gulf of Maine (GOM) plays in determining the distribution of dense aggregations of copepods. These aggregations are an important part of the marine ecosystem, especially for endangered North Atlantic right whales. Certain ocean processes may generate dense copepod aggregations, while others may destroy them; this thesis looks at how different characteristics of the GOM circulation fit into these two categories. The first part of the thesis investigates a hypothetical aggregation mechanism in which frontal circulation interacts with copepod behavior to generate a dense patch of copepods. The first two chapters of this thesis address this mechanism in the context of coastal river plumes and salinity fronts. One chapter describes the characteristics and variability of coastal freshwater and salinity fronts using a historical dataset and a realistic numerical model. The seasonal variability of freshwater is tied in part to seasonality in river discharge, while variability on shorter time scales in the frontal position is related to wind stress. Another chapter applies the hypothetical mechanism to idealized river plumes using a suite of numerical models. The structure of the plume and plume-relative circulation change the resulting copepod aggregation from what is expected from the hypothetical mechanism. The second part of the thesis discusses the GOM circulation and how it may eliminate copepod patches. The summertime mean surface circulation and eddy kinetic energy are computed from a Lagrangian drifter dataset and an adaptive technique that allows for higher spatial resolution while also keeping uncertainty low. Eddy diffusivity is also computed over different regions of the GOM in an attempt to quantify the spreading of a patch of copepods, and is found to be lower near the coast where right whales are often found feeding on copepod patches. In the next chapter, a numerical drifter dataset is used to understand how the results of the previous chapter depend upon the quantity of observations. It is found that the uncertainty in estimating eddy diffusivity is tightly coupled to the number of drifters in the calculation.

Description: This work was supported by the WHOI Coastal Ocean Institute Graduate Student Fellowship and Student Research Award, the WHOI Academic Programs O ce, the NOAA National Marine Fisheries Service Northeast Fisheries Science Center (NOAA Cooperative Agreement NA09OAR4320129), and the O ce of Naval Research Marine Mammals and Biology Program (Grant N00014-12-1-0208).

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

Description: This is the second part of a three-part paper on North American climate in phase 5 of the Coupled Model Intercomparison Project (CMIP5) that evaluates the twentieth-century simulations of intraseasonal to multidecadal variability and teleconnections with North American climate. Overall, the multimodel ensemble does reasonably well at reproducing observed variability in several aspects, but it does less well at capturing observed teleconnections, with implications for future projections examined in part three of this paper. In terms of intraseasonal variability, almost half of the models examined can reproduce observed variability in the eastern Pacific and most models capture the midsummer drought over Central America. The multimodel mean replicates the density of traveling tropical synoptic-scale disturbances but with large spread among the models. On the other hand, the coarse resolution of the models means that tropical cyclone frequencies are underpredicted in the Atlantic and eastern North Pacific. The frequency and mean amplitude of ENSO are generally well reproduced, although teleconnections with North American climate are widely varying among models and only a few models can reproduce the east and central Pacific types of ENSO and connections with U.S. winter temperatures. The models capture the spatial pattern of Pacific decadal oscillation (PDO) variability and its influence on continental temperature and West Coast precipitation but less well for the wintertime precipitation. The spatial representation of the Atlantic multidecadal oscillation (AMO) is reasonable, but the magnitude of SST anomalies and teleconnections are poorly reproduced. Multidecadal trends such as the warming hole over the central–southeastern United States and precipitation increases are not replicated by the models, suggesting that observed changes are linked to natural variability.

Description: The authors acknowledge the support of NOAA/Climate Program Office/Modeling, Analysis, Predictions and Projections (MAPP) program as part of the CMIP5 Task Force.

Description: 2014-06-01

Description: Author Posting. © American Meteorological Society, 2010. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 23 (2010): 6221–6233, doi:10.1175/2010JCLI3402.1.

Description: Enhanced decadal variability in sea surface temperature (SST) centered on the Kuroshio Extension (KE) has been found in the Community Climate System Model version 3 (CCSM3) as well as in other coupled climate models. This decadal peak has higher energy than is found in nature, almost twice as large in some cases. While previous analyses have concentrated on the mechanisms for such decadal variability in coupled models, an analysis of the causes of excessive SST response to changes in wind stress has been missing. Here, a detailed comparison of the relationships between interannual changes in SST and sea surface height (SSH) as a proxy for geostrophic surface currents in the region in both CCSM3 and observations, and how these relationships depend on the mean ocean circulation, temperature, and salinity, is made. We use observationally based climatological temperature and salinity fields as well as satellite-based SSH and SST fields for comparison. The primary cause for the excessive SST variability is the coincidence of the mean KE with the region of largest SST gradients in the model. In observations, these two regions are separated by almost 500 km. In addition, the too shallow surface oceanic mixed layer in March north of the KE in the subarctic Pacific contributes to the biases. These biases are not unique to CCSM3 and suggest that mean biases in current, temperature, and salinity structures in separated western boundary current regions can exert a large influence on the size of modeled decadal SST variability.

Description: Support for L.T. was provided by the NASA sponsored Ocean Surface Topography Science Team, under Contract 1267196 with the University of Washington, administered by the Jet Propulsion Laboratory. Support for Y.-O. K. comes from the NOAA Office of Global Programs (grant to C. Deser and Y.-O. Kwon) and the WHOI Heyman fellowship.

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

Description: A number of previous observational studies have found that the waters of the deep Pacific Ocean have an age, or elapsed time since contact with the surface, of 700–1000 yr. Numerical models suggest ages twice as old. Here, the authors present an inverse framework to determine the mean age and its upper and lower bounds given Global Ocean Data Analysis Project (GLODAP) radiocarbon observations, and they show that the potential range of ages increases with the number of constituents or sources that are included in the analysis. The inversion requires decomposing the World Ocean into source waters, which is obtained here using the total matrix intercomparison (TMI) method at up to 2° × 2° horizontal resolution with 11 113 surface sources. The authors find that the North Pacific at 2500-m depth can be no younger than 1100 yr old, which is older than some previous observational estimates. Accounting for the broadness of surface regions where waters originate leads to a reservoir-age correction of almost 100 yr smaller than would be estimated with a two or three water-mass decomposition and explains some of the discrepancy with previous observational studies. A best estimate of mean age is also presented using the mixing history along circulation pathways. Subject to the caveats that inference of the mixing history would benefit from further observations and that radiocarbon cannot rule out the presence of extremely old waters from exotic sources, the deep North Pacific waters are 1200–1500 yr old, which is more in line with existing numerical model results.

Description: GG is supported by the J. Lamar Worzel Assistant Scientist Fund and the Penzance Endowed Fund in Support of Assistant Scientists. PJH is supported by NSF Award 0960787.

Description: 2012-08-01

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

Description: Based on the Levitus atlas, we find that the application of the Montgomery streamfunction to the isopycnal surfaces induces an error which can not be ignored in some regions in the ocean. The error arises from the sloping effect of the specific volume anomaly along isopycnal surfaces. By including the major part of this effect, new streamfunctions, namely the pressure anomaly and main pressure streamfunctions, are suggested for the use in potential density coordinates. By using the newly proposed streamfunction and by including the variations of specific volume anomaly along isopycnal surfaces, the inverse model proposed by Hogg (1987) is modified for increasing accuracy and applied to the Brazil Basin to study the circulation, diffusion and water mass balances. The equations in the model, i.e. the dynamic equation, continuity equation, integrated vorticity equation, and conservation equations for heat, salt and oxygen (in which a consumption sink term is allowed), are written in centered finite difference form with lateral steps of 2 degree latitude and longitude and 8 levels in the vertical. This system of equations with constraints of positive diffusivities and oxygen consumption rates is solved by the inverse method. The results indicate that the circulation in the upper oceans is consistent with previous works, but that in the deep ocean is quite different. In the NADW region, we find a coincidence of the flows with the tongues of water properties. The diffusivities and diapycnal velocities seem stronger in the region near the equator than in the south, with reasonable values. Diffusion plays an important role in the water mass balance. Examples show that similar property fields may results from different processes.

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

Description: Six hydrographic basinwide sections, two in each of the three major ocean basins, are employed in a set of inverse calculations to determine the extent of exchange between the Pacific and Indian Oceans through the Indonesian Archipelago and the net global oceanic heat flux at 30°S. Using a model which combines the data for the South Pacific and South Indian Oceans, it is found that even the largest existing estimates of Indonesian Passage through flow (20 Sv) are consistent with the data. However, the available information cannot limit the extent of the exchange, i.e. both smaller and larger through flows produce physically reasonable circulation patterns. The seasonal and interannual variations which have been found by other investigators and which we are incapable of resolving, lead us to conclude that in the long term mean an estimate of ~10 Sv for the through flow is most reasonable. Globally, at 30°S, we find a net oceanic heat flux of -1.1 ± 1.7 PW, which is not significantly different from zero. It is dominated by a large (〉1 PW) southward heat flux in the Indian Ocean. Large equatorward (~0.8 PW) heat flux values in the South Atlantic Basin are not consistent with our data. We therefore conclude that although our data are consistent with some water following the warm water return path for NADW (Gordon 1986), the cold water path must play the dominant role in the maintenance of the global thermohaline cell associated with the formation process of NADW.

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

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

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

Description: 2014-09-01

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

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

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

Description: 2015-04-01

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

Description: A tomographic array is placed in a 2-layer, flat bottom, steady-wind driven quasi-geostrophic circulation model to investigate whether the analysis of acoustic travel time changes can detect large-scale barotropic oscillations. Time series of sea surface elevation and upper and lower layer meridional currents are generated for comparison against a series of acoustic travel times. The spectra of these time series exhibit a broad mesoscale peak near a period of 40 days. The spectrum of the acoustic travel time contains a significant peak due to a resonant barotropic oscillation with a period of 28.6 days which is not present in the spectra of the point measurements. In this numerical model, basin-scale tomographic measurements are a better method of sensing the large-scale resonant barotropic oscillations than are conventional point measurements because the tomographic system attenuates the "noise" from the mesoscale.

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

Description: The vertical dispersion of a tracer released on a density surface near 1500-m depth in the Antarctic Circumpolar Current west of Drake Passage indicates that the diapycnal diffusivity, averaged over 1 yr and over tens of thousands of square kilometers, is (1.3 ± 0.2) × 10−5 m2 s−1. Diapycnal diffusivity estimated from turbulent kinetic energy dissipation measurements about the area occupied by the tracer in austral summer 2010 was somewhat less, but still within a factor of 2, at (0.75 ± 0.07) × 10−5 m2 s−1. Turbulent diapycnal mixing of this intensity is characteristic of the midlatitude ocean interior, where the energy for mixing is believed to derive from internal wave breaking. Indeed, despite the frequent and intense atmospheric forcing experienced by the Southern Ocean, the amplitude of finescale velocity shear sampled about the tracer was similar to background amplitudes in the midlatitude ocean, with levels elevated to only 20%–50% above the Garrett–Munk reference spectrum. These results add to a long line of evidence that diapycnal mixing in the interior middepth ocean is weak and is likely too small to dictate the middepth meridional overturning circulation of the ocean.

Description: This material is based upon work supported by the National Science Foundation Grants OCE-0622825,OCE-0622670, OCE-0622630, and OCE-0623177.

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

Description: Climate change west of the Antarctic Peninsula is the most rapid of anywhere in the Southern Hemisphere, with associated changes in the rates and distributions of freshwater inputs to the ocean. Here, results from the first comprehensive survey of oxygen isotopes in seawater in this region are used to quantify spatial patterns of meteoric water (glacial discharge and precipitation) separately from sea ice melt. High levels of meteoric water are found close to the coast, due to orographic effects on precipitation and strong glacial discharge. Concentrations decrease offshore, driving significant southward geostrophic flows (up to ~30 cm s−1). These produce high meteoric water concentrations at the southern end of the sampling grid, where collapse of the Wilkins Ice Shelf may also have contributed. Sea ice melt concentrations are lower than meteoric water and patchier because of the mobile nature of the sea ice itself. Nonetheless, net sea ice production in the northern part of the sampling grid is inferred; combined with net sea ice melt in the south, this indicates an overall southward ice motion. The survey is contextualized temporally using a decade-long series of isotope data from a coastal Antarctic Peninsula site. This shows a temporal decline in meteoric water in the upper ocean, contrary to expectations based on increasing precipitation and accelerating deglaciation. This is driven by the increasing occurrence of deeper winter mixed layers and has potential implications for concentrations of trace metals supplied to the euphotic zone by glacial discharge. As the regional freshwater system evolves, the continuing isotope monitoring described here will elucidate the ongoing impacts on climate and the ecosystem.

Description: The Palmer LTER participants acknowledge Award 0823101 from the Organisms and Ecosystems program in NSF OPP

Description: 2013-09-01

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

Description: Estimates of the recent mean and time varying water mass transformation rates associated with North Atlantic surface-forced overturning are presented. The estimates are derived from heat and freshwater surface fluxes and sea surface temperature fields from six atmospheric reanalyses—the Japanese 25-yr Reanalysis (JRA), the NCEP–NCAR reanalysis (NCEP1), the NCEP–U.S. Department of Energy (DOE) reanalysis (NCEP2), the European Centre for Medium-Range Weather Forecasts (ECMWF) Interim Re-Analysis (ERA-I), the Climate Forecast System Reanalysis (CFSR), and the Modern-Era Reanalysis for Research and Applications (MERRA)—together with sea surface salinity fields from two globally gridded datasets (World Ocean Atlas and Met Office EN3 datasets). The resulting 12 estimates of the 1979–2007 mean surface-forced streamfunction all depict a subpolar cell, with maxima north of 45°N, near σ = 27.5 kg m−3, and a subtropical cell between 20° and 40°N, near σ = 26.1 kg m−3. The mean magnitude of the subpolar cell varies between 12 and 18 Sv (1 Sv ≡ 106 m3 s−1), consistent with estimates of the overturning circulation from subsurface observations. Analysis of the thermal and haline components of the surface density fluxes indicates that large differences in the inferred low-latitude circulation are largely a result of the biases in reanalysis net heat flux fields, which range in the global mean from −13 to 19 W m−2. The different estimates of temporal variability in the subpolar cell are well correlated with each other. This suggests that the uncertainty associated with the choice of reanalysis product does not critically limit the ability of the method to infer the variability in the subpolar overturning. In contrast, the different estimates of subtropical variability are poorly correlated with each other, and only a subset of them captures a significant fraction of the variability in independently estimated North Atlantic Subtropical Mode Water volume.

Description: JPG is funded by UK Natural Environment Research Council New Investigator Grant NE/I001654/1. Y-OK was supported by the U.S. National Science Foundation under Grant OCE-0424492. RJB is supported by a fellowship from the UK National Centre for Earth Observation.

Description: 2014-11-15

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

Description: Benthic foraminiferal δ13C, Cd/Ca, and Ba/Ca are important tools for reconstructing nutrient distributions, and thus ocean circulation, on glacial-interglacial timescales. However, each tracer has its own "artifacts" that can complicate paleoceanographic interpretations. It is therefore advantageous to measure multiple nutrient proxies with the aim of separating the various complicating effects. Zn/Ca is introduced as an important aid toward this goal. Benthic (Hoeglundina elegans) Cd/Ca ratios from the Bahama Banks indicate that the North Atlantic subtropical gyre was greatly depleted in nutrients during the last glacial maximum (LGM). A high-resolution Cd/Ca record from 965 m water depth suggests that Glacial North Atlantic Intermediate Water formation was strong during the LGM, weakened during the deglaciation, and strengthened again during the Younger Dryas cold period. Comparison of Cd/Ca and δ13C data reveals apparent short-term changes in carbon isotopic air-sea signatures. Benthic foraminiferal Zn/Ca could be a sensitive paleoceanographic tracer because deep water masses have characteristic Zn concentrations that increase about ten-fold from the deep North Atlantic to the deep North Pacific. A "core top calibration" shows that Zn/Ca is controlled by bottom water dissolved Zn concentration and, like Cd/Ca and BalCa, by bottom water saturation state with respect to calcite Since Zn/Ca responds to a different range of saturation states than Cd/Ca, the two may be used together to evaluate changes in deep water carbonate ion (CO32-) concentration. Zn/Ca and Cd/Ca ratios in the benthic foraminifer Cibicidoides wuellerstorfi exhibit large fluctuations over the past 100,000 years in a deep (3851 m) eastern equatorial Pacific sediment core. The data imply that bottom water CO32- concentrations were lowest during glacial Marine Isotope Stage 4 and highest during the last deglaciation. LGM CO32- concentrations appear to have been within a few μmol kg-1 of modern values. Deep North Atlantic Cd/Ca ratios imply much higher nutrient concentrations during the LGM. Although such data have usually been explained by a northward penetration of Southern Ocean Water (SOW), it has been suggested that they could result from increased preformed nutrient levels in the high-latitude North Atlantic or by increased aging of lower North Atlantic Deep Water (NADW). Glacial Zn/Ca data, however, require a substantially increased mixing with SOW and thus a reduction in NADW formation. Large changes in carbon isotopic air-sea exchange are invoked to reconcile benthic δ13C and trace metal data.

Description: This work was supported by a JOIlUSSAC Ocean Drilling Fellowship (subgrant JSG-CY 12-4), the R. H. Cole Ocean Ventures Fund, the Joint Program Education Office, and the National Science Foundation (grants OCE-9402804 and OCE-9503135 to W. Curry, and grant OCE-9633499 to D. Oppo).