Ramsey R. Harcourt

Enhanced Turbulence and Energy Dissipation at Ocean Fronts

Energy in surface ocean currents can dissipate into deep water via enhanced turbulence at the boundaries between water masses. The ocean surface boundary layer mediates air-sea exchange. In the classical paradigm and in current climate models, its turbulence is driven by atmospheric forcing. Observations at a 1-kilometer-wide front within the Kuroshio Current indicate that the rate of energy dissipation within the boundary layer is enhanced by one to two orders of magnitude, suggesting that the front, rather than the atmospheric forcing, supplied the energy for the turbulence. The data quantitatively support the hypothesis that winds aligned with the frontal velocity catalyzed a release of energy from the front to the turbulence. The resulting boundary layer is stratified in contrast to the classically well-mixed layer. These effects will be strongest at the intense fronts found in the Kuroshio Current, the Gulf Stream, and the Antarctic Circumpolar Current, all of which are key players in the climate system.

Large-Eddy Simulation of Langmuir Turbulence in Pure Wind Seas

Abstract The scaling of turbulent kinetic energy (TKE) and its vertical component (VKE) in the upper ocean boundary layer, forced by realistic wind stress and surface waves including the effects of Langmuir circulations, is investigated using large-eddy simulations (LESs). The interaction of waves and turbulence is modeled by the Craik–Leibovich vortex force. Horizontally uniform surface stress τ0 and Stokes drift profiles uS(z) are specified from the 10-m wind speed U10 and the surface wave age CP/U10, where CP is the spectral peak phase speed, using an empirical surface wave spectra and an associated wave age–dependent neutral drag coefficient CD. Wave-breaking effects are not otherwise included. Mixed layer depths HML vary from 30 to 120 m, with 0.6 ≤ CP/U10 ≤ 1.2 and 8 m s−1 < U10 < 70 m s−1, thereby addressing most possible oceanic conditions where TKE production is dominated by wind and wave forcing. The mixed layer–averaged “bulk” VKE 〈w2〉/u*2 is equally sensitive to the nondimensional Stokes e-fol...

Determining Vertical Water Velocities from Seaglider

Vertical velocities in the worlds oceans are typically small, less than 1 cm/s, posing a significant challenge to observation techniques. Seaglider, an autonomous profiling instrument, can be used to estimate vertical water velocity in the ocean to about half a centimeter per second. Using a Seaglider flight model and pressure observations, vertical water velocities are estimated along glider trajectories in the Labrador Sea before, during and after deep convection. Results indicate that vertical velocities in the stratified ocean agree with theoretical WKB-scaling of w, and in the turbulent mixed layer, scale with buoyancy and wind forcing. We estimate that accuracy is within 0.6 cm/s. Due to uncertainties in the flight model, velocities are poor near the surface and deep apogees, and during extended roll maneuvers. Some of this may be improved by using a dynamic flight model permitting acceleration, and by better constraining flight parameters through pilot choices during the mission.

The latmix summer campaign: Submesoscale stirring in the upper ocean

AbstractLateral stirring is a basic oceanographic phenomenon affecting the distribution of physical, chemical, and biological fields. Eddy stirring at scales on the order of 100 km (the mesoscale) is fairly well understood and explicitly represented in modern eddy-resolving numerical models of global ocean circulation. The same cannot be said for smaller-scale stirring processes. Here, the authors describe a major oceanographic field experiment aimed at observing and understanding the processes responsible for stirring at scales of 0.1–10 km. Stirring processes of varying intensity were studied in the Sargasso Sea eddy field approximately 250 km southeast of Cape Hatteras. Lateral variability of water-mass properties, the distribution of microscale turbulence, and the evolution of several patches of inert dye were studied with an array of shipboard, autonomous, and airborne instruments. Observations were made at two sites, characterized by weak and moderate background mesoscale straining, to contrast diff...

Fully Lagrangian Floats in Labrador Sea Deep Convection: Comparison of Numerical and Experimental Results

Measurements of deep convection from fully Lagrangian floats deployed in the Labrador Sea during February and March 1997 are compared with results from model drifters embedded in a large eddy simulation (LES) of the rapidly deepening mixed layer. The deep Lagrangian floats (DLFs) have a large vertical drag, and are designed to nearly match the density and compressibility of seawater. The high-resolution numerical simulation of deep convective turbulence uses initial conditions and surface forcing obtained from in situ oceanic and atmospheric observations made by the R/V Knorr. The response of model floats to the resolved large eddy fields of buoyancy and velocity is simulated for floats that are 5 g too buoyant, as well as for floats that are correctly ballasted. Mean profiles of potential temperature, Lagrangian rates of heating and acceleration, vertical turbulent kinetic energy (TKE), vertical heat flux, potential temperature variance, and float probability distribution functions (PDFs) are compared for actual and model floats. Horizontally homogeneous convection, as represented by the LES model, accounts for most of the first and second order statistics from float observations, except that observed temperature variance is several times larger than model variance. There are no correspondingly large differences in vertical TKE, heat flux, or mixed layer depth. The augmented temperature variance may be due to mixing across large-scale temperature and salinity gradients that are largely compensated in buoyancy. The rest of the DLF statistics agree well with the response of correctly ballasted model floats in the lowest 75% of the mixed layer, and are less consistent with results from buoyantly ballasted model floats. Other differences between observation and simulation in the mean profiles of heat flux, vertical TKE, and Lagrangian heating and vertical acceleration rates are confined to the upper quarter of the mixed layer. These differences are small contributions to layer-averaged quantities, but represent statistically significant profile features. Larger observed values of heat flux and vertical TKE in the upper quarter of the mixed layer are more consistent with model floats ballasted light. Float buoyancy, however, cannot fully account for the observed PDFs, temperature profiles, and Lagrangian rates of heating and acceleration. A test of Lagrangian self-consistency comparing vertical TKE and Lagrangian acceleration also shows that DLF measurements are not significantly affected by excess float buoyancy. These upper mixed layer features may instead be due to the interaction of wind-driven currents and baroclinicity.

A Second-Moment Closure Model of Langmuir Turbulence

AbstractThe Reynolds stress equation is modified to include the Craik–Leibovich vortex force, arising from the interaction of the phase-averaged surface wave Stokes drift with upper-ocean turbulence. An algebraic second-moment closure of the Reynolds stress equation yields an algebraic Reynolds stress model (ARSM) that requires a component of the vertical momentum flux to be directed down the gradient of the Stokes drift, in addition to the conventional component down the gradient of the ensemble-averaged Eulerian velocity. For vertical and horizontal component fluctuations, the momentum flux must be closed using the form , where the coefficient is generally distinct from the eddy viscosity or eddy diffusivity . Rational expressions for the stability functions , , and are derived for use in second-moment closure models where the turbulent velocity and length scales are dynamically modeled by prognostic equations for and . The resulting second-moment closure (SMC) includes the significant effects of the vo...

Measurement of Vertical Kinetic Energy and Vertical Velocity Skewness in Oceanic Boundary Layers by Imperfectly Lagrangian Floats

Abstract The effects of upward buoyancy on the accuracy with which Lagrangian floats can measure the Eulerian mean variance 〈ww〉E and skewness SwE of vertical fluid velocity w in the wind-driven upper-ocean boundary layer is investigated using both simulated floats in large-eddy simulation (LES) models and two float datasets. Nearly neutrally buoyant floats are repeatedly advected by the turbulent velocities across the boundary layer. Their vertical position Z is determined from pressure measurements; their W variance 〈WW〉F and skewness SWE are determined from the time series of float W = dZ/dt. If the float buoyancy is small, then the simulated floats measure the Eulerian velocity accurately; that is, δW2 = 〈WW〉F − 〈ww〉E and δSW = SWF − SwE are small compared to 〈ww〉E and SwE respectively. If the floats are buoyant, and thus have an upward vertical velocity Wbias relative to the water, then δW2 and δSW can become significant. Buoyancy causes the floats to oversample both shallow depths and strong vertica...

Characterizing Thermohaline Intrusions in the North Pacific Subtropical Frontal Zone

Abstract A monthlong field survey in July 2007, focused on the North Pacific subtropical frontal zone (STFZ) near 30°N, 158°W, combined towed depth-cycling conductivity–temperature–depth (CTD) profiling with shipboard current observations. Measurements were used to investigate the distribution and structure of thermohaline intrusions. The study revealed that local extrema of vertical salinity profiles, often used as intrusion indicators, were only a subset of a wider class of distortions in thermohaline fields due to interleaving processes. A new method to investigate interleaving based on diapycnal spiciness curvature was used to describe an expanded class of laterally coherent intrusions. STFZ intrusions were characterized by their overall statistics and by a number of case studies. Thermohaline interleaving was particularly intense within 5 km of two partially compensated fronts, where intrusions with both positive and negative salinity anomalies were widespread. The vertical and cross-frontal scales o...

Three-Dimensional Structure and Temporal Evolution of Submesoscale Thermohaline Intrusions in the North Pacific Subtropical Frontal Zone

Abstract Four instances of persistent intrusive deformation of the North Pacific Subtropical Front were tagged individually by a Lagrangian float and tracked for several days. Each feature was mapped in three dimensions using repeat towed observations referenced to the float. Isohaline surface deformations in the frontal zone included sheetlike folds elongated in the alongfront direction and narrow tongues extending across the front. All deformations appeared as protrusions of relatively cold, and fresh, water across the front. No corresponding features of the opposite sign or isolated lenslike structures were observed. The sheets were O(10 m) thick, protruded about 10 km into the warm saline side of the front, and were coherent for 10–30 km along the front. Having about the same thickness and cross-frontal extent as the sheets, tongues extended less than 5 km along the front. All of the intrusions persisted as long as they were followed, several days to one week. Their structures evolved on both inertial...

Langmuir Turbulence under Hurricane Gustav (2008)

Extreme winds and complex wave fields drive upper-ocean turbulence in tropical cyclone conditions. MotivatedbyLagrangianfloatobservationsofbulkverticalvelocityvariance(VVV)underHurricaneGustav (2008), upper-ocean turbulence is investigated based on large-eddy simulation (LES) of the wave-averaged Navier–Stokesequations. To realistically capture wind- and wave-driven Langmuir turbulence (LT), the LES model imposes the Stokes drift vector from spectral wave simulations; both the LES and wave model are forced by the NOAA Hurricane Research Division (HRD) surface wind analysis product. Results strongly suggest that without LT effects simulated VVV underestimates the observed VVV. LT increases the VVV, indicating that it plays a significant role in upper-ocean turbulence dynamics. Consistent with observations, theLES predictsasuppressionofVVVnearthehurricaneeyeduetowind-wavemisalignment. However,this decrease is weaker and of shorter duration than that observed, potentially due to large-scale horizontal advection not present in the LES. Both observations and simulations are consistent with a highly variable upper ocean turbulence field beneath tropical cyclone cores. Bulk VVV, a TKE budget analysis, and anisotropy coefficient(ratioofhorizontaltoverticalvelocityvariances)profilesallindicatethatLTissuppressedtolevels closer to that of shear turbulence (ST) due to misaligned wind and wave fields. VVV approximately scales with the directional surface layer Langmuir number. Such a scaling provides guidance for the development of an upper-ocean boundary layer parameterization that explicitly depends on sea state.