Inhibited upper ocean restratification in nonequilibrium swell conditions

Author Posting. © American Geophysical Union, 2013. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 40 (2013): 3672–3676, doi:10.1002/grl.50708.Diurnal restratification of the ocean surface boundary layer (OSBL) represents a competition between mixing of the OSBL and solar heating. Langmuir turbulence (LT) is a mixing process in the OSBL, driven by wind and surface waves, that transfers momentum, heat, and mass. Observations in nonequilibrium swell conditions reveal that the OSBL does not restratify despite low winds and strong solar radiation. Motivated by observations, we use large-eddy simulations of the wave-averaged Navier-Stokes equations to show that LT is capable of inhibiting diurnal restratification of the OSBL. Incoming heat is redistributed vertically by LT, forming a warmer OSBL with a nearly uniform temperature. The inhibition of restratification is not reproduced by two common Reynolds-averaged Navier-Stokes equation models, highlighting the importance of properly representing sea-state dependent LT dynamics in OSBL models.This work was supported by the U.S. National Science Foundation (Grant OCE-1130678).2014-01-3

Rapid mixed layer depening by the combination of Langmuir and shear instabilities : a case study

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): 2381-2400, doi:10.1175/2010JPO4403.1.Langmuir circulation (LC) is a turbulent upper-ocean process driven by wind and surface waves that contributes significantly to the transport of momentum, heat, and mass in the oceanic surface layer. The authors have previously performed a direct comparison of large-eddy simulations and observations of the upper-ocean response to a wind event with rapid mixed layer deepening. The evolution of simulated crosswind velocity variance and spatial scales, as well as mixed layer deepening, was only consistent with observations if LC effects are included in the model. Based on an analysis of these validated simulations, in this study the fundamental differences in mixing between purely shear-driven turbulence and turbulence with LC are identified. In the former case, turbulent kinetic energy (TKE) production due to shear instabilities is largest near the surface, gradually decreasing to zero near the base of the mixed layer. This stands in contrast to the LC case in which at middepth range TKE production can be dominated by Stokes drift shear. Furthermore, the Eulerian mean vertical shear peaks near the base of the mixed layer so that TKE production by mean shear flow is elevated there. LC transports horizontal momentum efficiently downward leading to an along-wind velocity jet below LC downwelling regions at the base of the mixed layer. Locally enhanced vertical shear instabilities as a result of this jet efficiently erode the thermocline. In turn, enhanced breaking internal waves inject cold deep water into the mixed layer, where LC currents transport temperature perturbation advectively. Thus, LC and locally generated shear instabilities work intimately together to facilitate strongly the mixed layer deepening process.This research was supported by the Office of Naval Research through Grants N00014-09-M-0112 (TK) and N00014-06-1-0178 (AP, JT). Author TK also received support from a Woods Hole Oceanographic Institution Cooperative Institute for Climate and Ocean Research Postdoctoral Scholarship

Significance of Langmuir circulation in upper ocean mixing : comparison of observations and simulations

Author Posting. © American Geophysical Union, 2009. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 36 (2009): L10603, doi:10.1029/2009GL037620.Representing upper ocean turbulence accurately in models remains a great challenge for improving weather and climate projections. Langmuir circulation (LC) is a turbulent process driven by wind and surface waves that plays a key role in transferring momentum, heat, and mass in the oceanic surface layer. We present a direct comparison between observations and large eddy simulations, based on the wave-averaged Navier-Stokes equation, of an LC growth event. The evolution of cross-wind velocity variance and spatial scales, as well as mixed layer deepening are only consistent with simulations if LC effects are included in the model. Our results offer a validation of the large eddy simulation approach to understanding LC dynamics, and demonstrate the importance of LC in ocean surface layer mixing.This research was supported by the Office of Naval Research through grants N00014-09-M-0112 (TK) and N00014-06-1-0178 (AP, JT). TK also received support from a Woods Hole Oceanographic Institution Cooperative Institute for Climate and Ocean Research Postdoctoral Scholarship

Langmuir Turbulence Parameterization in Tropical Cyclone Conditions

The Stokes drift of surface waves significantly modifies the upper-ocean turbulence because of the Craik–Leibovich vortex force (Langmuir turbulence). Under tropical cyclones the contribution of the surface waves varies significantly depending on complex wind and wave conditions. Therefore, turbulence closure models used in ocean models need to explicitly include the sea state–dependent impacts of the Langmuir turbulence. In this study, the K-profile parameterization (KPP) first-moment turbulence closure model is modified to include the explicit Langmuir turbulence effect, and its performance is tested against equivalent large-eddy simulation (LES) experiments under tropical cyclone conditions. First, the KPP model is retuned to reproduce LES results without Langmuir turbulence to eliminate implicit Langmuir turbulence effects included in the standard KPP model. Next, the Lagrangian currents are used in place of the Eulerian currents in the KPP equations that calculate the bulk Richardson number and the vertical turbulent momentum flux. Finally, an enhancement to the turbulent mixing is introduced as a function of the nondimensional turbulent Langmuir number. The retuned KPP, with the Lagrangian currents replacing the Eulerian currents and the turbulent mixing enhanced, significantly improves prediction of upper-ocean temperature and currents compared to the standard (unmodified) KPP model under tropical cyclones and shows improvements over the standard KPP at constant moderate winds (10 m s−1)

Impact of Sea-State-Dependent Langmuir Turbulence on the Ocean Response to a Tropical Cyclone

Tropical cyclones are fueled by the air–sea heat flux, which is reduced when the ocean surface cools due to mixed layer deepening and upwelling. Wave-driven Langmuir turbulence can significantly modify these processes. This study investigates the impact of sea-state-dependent Langmuir turbulence on the three-dimensional ocean response to a tropical cyclone in coupled wave–ocean simulations. The Stokes drift is computed from the simulated wave spectrum using the WAVEWATCH III wave model and passed to the three-dimensional Princeton Ocean Model. The Langmuir turbulence impact is included in the vertical mixing of the ocean model by adding the Stokes drift to the shear of the vertical mean current and by including Langmuir turbulence enhancements to the K-profile parameterization (KPP) scheme. Results are assessed by comparing simulations with explicit (sea-state dependent) and implicit (independent of sea state) Langmuir turbulence parameterizations, as well as with turbulence driven by shear alone. The results demonstrate that the sea-state-dependent Langmuir turbulence parameterization significantly modifies the three-dimensional ocean response to a tropical cyclone. This is due to the reduction of upwelling and horizontal advection where the near-surface currents are reduced by Langmuir turbulence. The implicit scheme not only misses the impact of sea-state dependence on the surface cooling, but it also misrepresents the impact of the Langmuir turbulence on the Eulerian advection. This suggests that explicitly resolving the sea-state-dependent Langmuir turbulence will lead to increased accuracy in predicting the ocean response in coupled tropical cyclone–ocean models

The effect of breaking waves on a coupled model of wind and ocean surface waves. Part I : mature seas

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.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.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)

A Novel Approach to Flow Estimation in Tidal Rivers

Reliable estimation of river discharge to the ocean from large tidal rivers is vital for water resources management and climate analyses. Due to the difficulties inherent in measuring tidal-river discharge, flow records are often limited in length and/or quality and tidal records often predate discharge records. Tidal theory indicates that tides and river discharge interact through quadratic bed friction, which diminishes and distorts the tidal wave as discharge increases. We use this phenomenon to develop a method of estimating river discharge for time periods with tidal data but no flow record. Employing sequential 32 day harmonic analyses of tidal properties, we calibrate San Francisco (SF), CA tide data to the Sacramento River delta outflow index from 1930 to 1990, and use the resulting relationship to hindcast river flow from 1858 to 1929. The M2 admittance (a ratio of the observed M2 tidal constituent to its astronomical forcing) best reproduces high flows, while low-flow periods are better represented by amplitude ratios based on higher harmonics (e.g.). Results show that the annual inflow to SF Bay is now 30% less than before 1900 and confirm that the flood of January 1862 was the largest since 1858

The effect of breaking waves on a coupled model of wind and ocean surface waves. Part II : growing seas

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.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.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)

Langmuir Turbulence under Hurricane Gustav (2008)

Extreme winds and complex wave fields drive upper-ocean turbulence in tropical cyclone conditions. Motivated by Lagrangian float observations of bulk vertical velocity variance (VVV) under Hurricane Gustav (2008), upper-ocean turbulence is investigated based on large-eddy simulation (LES) of the wave-averaged Navier–Stokes equations. 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, the LES predicts a suppression of VVV near the hurricane eye due to wind-wave misalignment. 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 (ratio of horizontal to vertical velocity variances) profiles all indicate that LT is suppressed to levels 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

The physical oceanography of the transport of floating marine debris

Marine plastic debris floating on the ocean surface is a major environmental problem. However, its distribution in the ocean is poorly mapped, and most of the plastic waste estimated to have entered the ocean from land is unaccounted for. Better understanding of how plastic debris is transported from coastal and marine sources is crucial to quantify and close the global inventory of marine plastics, which in turn represents critical information for mitigation or policy strategies. At the same time, plastic is a unique tracer that provides an opportunity to learn more about the physics and dynamics of our ocean across multiple scales, from the Ekman convergence in basin-scale gyres to individual waves in the surfzone. In this review, we comprehensively discuss what is known about the different processes that govern the transport of floating marine plastic debris in both the open ocean and the coastal zones, based on the published literature and referring to insights from neighbouring fields such as oil spill dispersion, marine safety recovery, plankton connectivity, and others. We discuss how measurements of marine plastics (both in situ and in the laboratory), remote sensing, and numerical simulations can elucidate these processes and their interactions across spatio-temporal scales