Chalk-Ex—fate of CaCO3 particles in the mixed layer : evolution of patch optical properties

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 Journal of Geophysical Research 114 (2009): C07020, doi:10.1029/2008JC004902.The fate of particles in the mixed layer is of great relevance to the global carbon cycle as well as to the propagation of light in the sea. We conducted four manipulative field experiments called “Chalk-Ex” in which known quantities of uniform, calcium carbonate particles were injected into the surface mixed layer. Since the production term for these patches was known to high precision, the experimental design allowed us to focus on terms associated with particle loss. The mass of chalk in the patches was evaluated using the well-calibrated light-scattering properties of the chalk plus measurements from a variety of optical measurements and platforms. Patches were surveyed with a temporal resolution of hours over spatial scales of tens of kilometers. Our results demonstrated exponential loss of the chalk particles with time from the patches. There was little evidence for rapid sinking of the chalk. Instead, horizontal eddy diffusion appeared to be the major factor affecting the dispersion of the chalk to concentrations below the limits of detection. There was unequivocal evidence of subduction of the chalk along isopycnals and subsequent formation of thin layers. Shear dispersion is the most likely mechanism to explain these results. Calculations of horizontal eddy diffusivity were consistent with other mixed layer patch experiments. Our results provide insight into the importance of physics in the formation of subsurface particle maxima in the sea, as well as the importance of rapid coccolith production and critical patch size for maintenance of natural coccolithophore blooms in nature.We would like to thank the Office of Naval Research/Optical and Biological Oceanography Program for their support of Chalk-Ex with awards N000140110042 (WMB) and N00014-01-1-0141 (AJP). Additional funding for this work came from ONR (N00014-05-1- 0111) and NASA (NNG04Gl11G, NNX08AC27G, NNG04HZ25C) to W.M.B

Turbulent energy production and entrainment at a highly stratified estuarine front

Author Posting. © American Geophysical Union, 2004. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 109 (2004): C05004, doi:10.1029/2003JC002094.Rates of turbulent kinetic energy (TKE) production and buoyancy flux in the region immediately seaward (~1 km) of a highly stratified estuarine front at the mouth of the Fraser River (British Columbia, Canada) are calculated using a control volume approach. The calculations are based on field data obtained from shipboard instrumentation, specifically velocity data from a ship mounted acoustic Doppler current profiler (ADCP), and salinity data from a towed conductivity-temperature-depth (CTD) unit. The results allow for the calculation of vertical velocities in the water column, and the total vertical transport of salt and momentum. The vertical turbulent transport quantities (inline equation, inline equation) can then be estimated as the difference between the total transport and the advective transport. Estimated production is on the order of 10−3 m2 s−3, yielding a value of ɛ(νN2)−1 on the order of 104. This rate of TKE production is at the upper limit of reported values for ocean and coastal environments. Flux Richardson numbers in this highly energetic system generally range from 0.15 to 0.2, with most mixing occurring at gradient Richardson numbers slightly less than inline equation. These values compare favorably with other values in the literature that are associated with turbulence observations from regimes characterized by scales several orders of magnitude smaller than are present in the Fraser River.This work was performed as a part of D. MacDonald’s Ph.D. thesis, and was funded by Office of Naval Research grants N000-14-97-10134 and N000-14-97- 10566, National Science Foundation grant OCE-9906787, a National Science Foundation graduate fellowship, and support from the WHOI Academic Programs Office

Large-eddy simulation of the tidal-cycle variations of an estuarine boundary layer

Author Posting. © American Geophysical Union, 2010. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 115 (2010): C08003, doi:10.1029/2009JC005702.The estuarine boundary layer affected by a horizontal density gradient exhibits temporal evolution over a tidal cycle, in a manner similar to the diurnal cycle of the ocean surface mixed layer. A large eddy simulation (LES) model is developed to investigate the physics controlling the growth of the boundary layer during the flood tide and restratification during the ebb tide. Turbulent kinetic energy, momentum and salt fluxes, bottom stress, and energy dissipation rates calculated from the LES model all show a strong flood-ebb asymmetry. Analysis of the turbulent kinetic energy (TKE) budget shows a primary balance between shear production and dissipation in the well-mixed boundary layer over the tidal cycle. However, TKE transport term is found to be important across the edge of the boundary layer during the flood tide so turbulent energy generated in the bottom boundary layer can be transferred to the stratified pycnocline region. Tidal straining leads to a small and weakly convective region inside the boundary layer during the flood tide but the strain-induced buoyancy flux does not make a significant contribution to the turbulence generation. Additional LES runs are conducted by switching off the baroclinic pressure gradient term in the momentum equation and the tidal straining term in the salinity equation to show that the baroclinic pressure gradient is the main mechanism responsible for generating the flood-ebb mixing asymmetry.This work is supported by grants OCE-0451699 (M.L.), OCE-0452380 (U.P. and S.R.), and OCE-0451740 (W.R.G.) from the National Science Foundation

The modification of bottom boundary layer turbulence and mixing by internal waves shoaling on a barrier reef

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): 2223–2241, doi:10.1175/2011JPO4344.1.Results are presented from an observational study of stratified, turbulent flow in the bottom boundary layer on the outer southeast Florida shelf. Measurements of momentum and heat fluxes were made using an array of acoustic Doppler velocimeters and fast-response temperature sensors in the bottom 3 m over a rough reef slope. Direct estimates of flux Richardson number Rf confirm previous laboratory, numerical, and observational work, which find mixing efficiency not to be a constant but rather to vary with Frt, Reb, and Rig. These results depart from previous observations in that the highest levels of mixing efficiency occur for Frt < 1, suggesting that efficient mixing can also happen in regions of buoyancy-controlled turbulence. Generally, the authors find that turbulence in the reef bottom boundary layer is highly variable in time and modified by near-bed flow, shear, and stratification driven by shoaling internal waves.Funding was provided by grants from the National Oceanic and Atmospheric Administration’s National Undersea Research Program, National Science Foundation Grants OCE-0622967 and OCE- 0824972 to SGM, and the Singapore Stanford Program. Kristen Davis was supported by a National Defense Science and Engineering Graduate Fellowship and an ARCS Foundation Fellowship

Stratified turbulence in the nearshore coastal ocean: Dynamics and evolution in the presence of internal bores

High-frequency measurements of stratified turbulence throughout the water column were collected over a 2 week period in the nearshore environment of southern Monterey Bay, CA, using a cabled observatory system and an underwater turbulence flux tower. The tower contained a vertical array of acoustic Doppler velocimeters and fast-response conductivity-temperature sensors, providing a nearly continuous data set of turbulent velocity and density fluctuations and a unique look into the stratified turbulence field. The evolution of various turbulence quantities and direct measurements of the vertical turbulent diffusivity is examined in the presence of nearshore internal bores, both in the near-bed region and in the stratified interior. We show that individual bores can drive substantial changes in local turbulence and mixing dynamics, with considerable differences between the leading and trailing edges of the bores. Using direct observations of the flux Richardson number, our measurements confirm previous observations that show the highest mixing efficiencies (Γ) occurring in regions of buoyancy-controlled turbulence. Parameterizations of the flux Richardson number as a function of the turbulence activity number are also presented. Finally, we demonstrate that the commonly used assumption of a constant mixing efficiency (Γ = 0.2) for calculating turbulent diffusivities leads to significant overestimates compared to diffusivity values calculated using the directly measured mixing efficiency. Implications of the results are discussed

Turbulence and Fossil Turbulence in Oceans and Lakes

Turbulence is defined as an eddy-like state of fluid motion where the inertial-vortex forces of the eddies are larger than any of the other forces that tend to damp the eddies out. Energy cascades of irrotational flows from large scales to small are non-turbulent, even if they supply energy to turbulence. Turbulent flows are rotational and cascade from small scales to large, with feedback. Viscous forces limit the smallest turbulent eddy size to the Kolmogorov scale. In stratified fluids, buoyancy forces limit large vertical overturns to the Ozmidov scale and convert the largest turbulent eddies into a unique class of saturated, non-propagating, internal waves, termed fossil-vorticity-turbulence. These waves have the same energy but different properties and spectral forms than the original turbulence patch. The Gibson (1980, 1986) theory of fossil turbulence applies universal similarity theories of turbulence and turbulent mixing to the vertical evolution of an isolated patch of turbulence in a stratified fluid as its growth is constrained and fossilized by buoyancy forces. These theories apply to the dynamics of atmospheric, astrophysical and cosmological turbulence.Comment: 31 pages, 11 figures, 2 tables, see http://www-acs.ucsd.edu/~ir118 Accepted for publication by the Chinese Journal of Oceanology and Limnolog

The Critical Richardson Number and Limits of Applicability of Local Similarity Theory in the Stable Boundary Layer

Measurements of atmospheric turbulence made over the Arctic pack ice during the Surface Heat Budget of the Arctic Ocean experiment (SHEBA) are used to determine the limits of applicability of Monin-Obukhov similarity theory (in the local scaling formulation) in the stable atmospheric boundary layer. Based on the spectral analysis of wind velocity and air temperature fluctuations, it is shown that, when both of the gradient Richardson number, Ri, and the flux Richardson number, Rf, exceed a 'critical value' of about 0.20 - 0.25, the inertial subrange associated with the Richardson-Kolmogorov cascade dies out and vertical turbulent fluxes become small. Some small-scale turbulence survives even in this supercritical regime, but this is non-Kolmogorov turbulence, and it decays rapidly with further increasing stability. Similarity theory is based on the turbulent fluxes in the high-frequency part of the spectra that are associated with energy-containing/flux-carrying eddies. Spectral densities in this high-frequency band diminish as the Richardson-Kolmogorov energy cascade weakens; therefore, the applicability of local Monin-Obukhov similarity theory in stable conditions is limited by the inequalities Ri < Ri_cr and Rf < Rf_cr. However, it is found that Rf_cr = 0.20 - 0.25 is a primary threshold for applicability. Applying this prerequisite shows that the data follow classical Monin-Obukhov local z-less predictions after the irrelevant cases (turbulence without the Richardson-Kolmogorov cascade) have been filtered out.Comment: Boundary-Layer Meteorology (Manuscript submitted: 16 February 2012; Accepted: 10 September 2012