Tobias Kukulka

Distribution of surface plastic debris in the eastern Pacific Ocean from an 11-Year data set

We present an extensive survey of floating plastic debris in the eastern North and South Pacific Oceans from more than 2500 plankton net tows conducted between 2001 and 2012. From these data we defined an accumulation zone (25 to 41 °N, 130 to 180 °W) in the North Pacific subtropical gyre that closely corresponds to centers of accumulation resulting from the convergence of ocean surface currents predicted by several oceanographic numerical models. Maximum plastic concentrations from individual surface net tows exceeded 10(6) pieces km(-2), with concentrations decreasing with increasing distance from the predicted center of accumulation. Outside the North Pacific subtropical gyre the median plastic concentration was 0 pieces km(-2). We were unable to detect a robust temporal trend in the data set, perhaps because of confounded spatial and temporal variability. Large spatiotemporal variability in plastic concentration causes order of magnitude differences in summary statistics calculated over short time periods or in limited geographic areas. Utilizing all available plankton net data collected in the eastern Pacific Ocean (17.4 °S to 61.0 °N; 85.0 to 180.0 °W) since 1999, we estimated a minimum of 21,290 t of floating microplastic.

Observations of Turbulence in the Ocean Surface Boundary Layer: Energetics and Transport

Observations of turbulent kinetic energy (TKE) dynamics in the ocean surface boundary layer are presented here and compared with results from previous observational, numerical, and analytic studies. As in previous studies, the dissipation rate of TKE is found to be higher in the wavy ocean surface boundary layer than it would be in a flow past a rigid boundary with similar stress and buoyancy forcing. Estimates of the terms in the turbulent kinetic energy equation indicate that, unlike in a flow past a rigid boundary, the dissipation rates cannot be balanced by local production terms, suggesting that the transport of TKE is important in the ocean surface boundary layer. A simple analytic model containing parameterizations of production, dissipation, and transport reproduces key features of the vertical profile of TKE, including enhancement near the surface. The effective turbulent diffusion coefficient for heat is larger than would be expected in a rigid-boundary boundary layer. This diffusion coefficient is predicted reasonably well by a model that contains the effects of shear production, buoyancy forcing, and transport of TKE (thought to be related to wave breaking). Neglect of buoyancy forcing or wave breaking in the parameterization results in poor predictions of turbulent diffusivity. Langmuir turbulence was detected concurrently with a fraction of the turbulence quantities reported here, but these times did not stand out as having significant differences from observations when Langmuir turbulence was not detected.

Impacts of Columbia River discharge on salmonid habitat: 2. Changes in shallow‐water habitat

[1] This is the second part of an investigation that analyzes human alteration of shallow-water habitat (SWH) available to juvenile salmonids in the tidal Lower Columbia River. Part 2 develops a one-dimensional, subtidal river stage model that explains ∼90% of the stage variance in the tidal river. This model and the tidal model developed in part 1 [Kukulka and Jay, 2003] uncouple the nonlinear interaction of river tides and river stage by referring both to external forcing by river discharge, ocean tides, and atmospheric pressure. Applying the two models, daily high-water levels were predicted for a reach from rkm-50 to rkm-90 during 1974 to 1998, the period of contemporary management. Predicted water levels were related to the bathymetry and topography to determine the changes in shallow-water habitat area (SWHA) caused by flood control dikes and altered flow management. Model results suggest that diking and a >40% reduction of peak flows have reduced SWHA by ∼62% during the crucial spring freshet period during which juvenile salmon use of SWHA is maximal. Taken individually, diking and flow cycle alteration reduced spring freshet SWHA by 52% and 29%, respectively. SWHA has been both displaced to lower elevations and modified in its character because tidal range has increased. Our models of these processes are economical for the very long simulations (seasons to centuries) needed to understand historic changes and climate impacts on SWH. Through analysis of the nonlinear processes controlling surface elevation in a tidal river, we have identified some of the mechanisms that link freshwater discharge to SWH and salmonid survival.

Rapid Mixed Layer Deepening by the Combination of Langmuir and Shear Instabilities: A Case Study

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 simulationsandobservationsof the upper-ocean response to a wind event with rapid mixed layer deepening. The evolution of simulated crosswind velocity varianceandspatialscales,aswellasmixedlayerdeepening,wasonlyconsistentwithobservationsifLCeffects 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.

A Model of the Air–Sea Momentum Flux and Breaking-Wave Distribution for Strongly Forced Wind Waves

Abstract Under high-wind conditions, breaking surface waves likely play an important role in the air–sea momentum flux. A coupled wind–wave model is developed based on the assumption that in the equilibrium range of surface wave spectra the wind stress is dominated by the form drag of breaking waves. By conserving both momentum and energy in the air and also imposing the wave energy balance, coupled equations are derived governing the turbulent stress, wind speed, and the breaking-wave distribution (total breaking crest length per unit surface area as a function of wavenumber). It is assumed that smaller-scale breaking waves are sheltered from wind forcing if they are in airflow separation regions of longer breaking waves (spatial sheltering effect). Without this spatial sheltering, exact analytic solutions are obtained; with spatial sheltering asymptotic solutions for small- and large-scale breakers are derived. In both cases, the breaking-wave distribution approaches a constant value for large wavenumbe...

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.

Langmuir Turbulence Parameterization in Tropical Cyclone Conditions

AbstractThe 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 a...

Passive buoyant tracers in the ocean surface boundary layer: 2. Observations and simulations of microplastic marine debris

This paper is the second of a two-part series that investigates passive buoyant tracers in the ocean surface boundary layer (OSBL). The first part examines the influence of equilibrium wind-waves on vertical tracer distributions, based on large eddy simulations (LESs) of the wave-averaged Navier-Stokes equation. Motivated by observations of buoyant microplastic marine debris (MPMD), this study applies the LES model and the parametric one-dimensional column model from part one to examine the vertical distributions of MPMD. MPMD is widely distributed in vast regions of the subtropical gyres and has emerged as a major open ocean pollutant whose distribution is subject to upper ocean turbulence. The models capture shear-driven turbulence, Langmuir turbulence (LT), and enhanced turbulent kinetic energy input due to breaking waves (BWs). Model results are only consistent with observations of MPMD profiles and the relationship between surface concentrations and wind speed if LT effects are included. Neither BW nor shear-driven turbulence is capable of deeply submerging MPMD, suggesting that the observed vertical MPMD distributions are a characteristic signature of wave-driven LT. Thus, this study demonstrates that LT substantially increases turbulent transport in the OSBL, resulting in deep submergence of buoyant tracers. The parametric model is applied to 11 years of observations in the North Atlantic and North Pacific subtropical gyres to show that surface measurements substantially underestimate MPMD concentrations by a factor of 3–13.

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

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

Passive buoyant tracers in the ocean surface boundary layer: 1. Influence of equilibrium wind‐waves on vertical distributions

This paper is the first of a two part series that investigates passive buoyant tracers in the ocean surface boundary layer. The first part examines the influence of equilibrium wind-waves on vertical tracer distributions, based on large eddy simulations (LES) of the wave-averaged Navier-Stokes equation. The second part applies the model to investigate observations of buoyant microplastic marine debris, which has emerged as a major ocean pollutant. The LES model captures both Langmuir turbulence (LT) and enhanced turbulent kinetic energy input due to breaking waves (BW) by imposing equilibrium wind-wave statistics for a range of wind and wave conditions. Concentration profiles of LES agree well with analytic solutions obtained for an eddy diffusivity profile that is constant near the surface and transitions into the K-Profile Parameterization (KPP) profile shape at greater depth. For a range of wind and wave conditions, the eddy diffusivity normalized by the product of water-side friction velocity and mixed layer depth, h, mainly depends on a single nondimensional parameter, the peak wavelength (which is related to Stokes drift decay depth) normalized by h. For smaller wave ages, BW critically enhances near-surface mixing, while LT effects are relatively small. For greater wave ages, both BW and LT contribute to elevated near-surface mixing, and LT significantly increases turbulent transport at greater depth. We identify a range of realistic wind and wave conditions for which only Langmuir (and not BW or shear driven) turbulence is capable of deeply submerging buoyant tracers.