skip to main content


Title: Turbulence in a small boreal lake: Consequences for air–water gas exchange
The hydrodynamics within small boreal lakes have rarely been studied, yet knowing whether turbulence at the air-water interface and in the water column scales with metrics developed elsewhere is essential for computing metabolism and fluxes of climate-forcing trace gases. We instrumented a humic, 4.7 ha, boreal lake with 2 meteorological stations, 3 thermistor arrays, an infra-red (IR) camera to quantify surface divergence, obtained turbulence as dissipation rate of turbulent kinetic energy (ε) using an acoustic Doppler velocimeter and a temperature-gradient microstructure profiler, and conducted chamber measurements for short periods to obtain fluxes and gas transfer velocities (k). Near-surface ε varied from 10-8 m2 s-3 to 10-6 m2 s-3 for the 0 to 4 m s-1 winds and followed predictions from Monin-Obukhov similarity theory. The coefficient of eddy diffusivity in the mixed layer was up to 10-3 m2 s-1 on the windiest afternoons, an order of magnitude less other afternoons, and near molecular at deeper depths. The upper thermocline upwelled when Lake numbers (LN) dropped below 4 facilitating vertical and horizontal exchange. k computed from a surface renewal model using ε agreed with values from chambers and surface divergence and increased linearly with wind speed. Diurnal thermoclines formed on sunny days when winds were < 3 m s-1, a condition that can lead to elevated near-surface ε and k. Results extend scaling approaches developed in the laboratory and for larger water bodies, illustrate turbulence and k are greater than expected in small wind-sheltered lakes, and provide new equations to quantify fluxes.  more » « less
Award ID(s):
1737411 1753856
NSF-PAR ID:
10211535
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Limnology and Oceanography
Volume:
9999
ISSN:
0024-3590
Page Range / eLocation ID:
1-28
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Measurements of turbulence, as rate of dissipation of turbulent kinetic energy (ε), adjacent to the air-water interface are rare but essential for understanding of gas transfer velocities (k) used to compute fluxes of greenhouse gases. Variability in ε is expected over diel cycles of stratification and mixing. Monin-Obukhov similarity theory (MOST) predicts an enhancement in ε during heating (buoyancy flux, β+) relative to that for shear (u*w 3/κz where u*w is water friction velocity, κ is von Karman constant, z is depth). To verify and expand predictions, we quantified ε in the upper 0.25 m and below from profiles of temperature-gradient microstructure in combination with time series meteorology and temperature in a tropical reservoir for winds <4 m s−1. Maximum likelihood estimates of near-surface ε during heating were independent of wind speed and high, ∼5 × 10−6 m2 s−3, up to three orders of magnitude higher than predictions from u*w 3/κz, increased with heating, and were ∼10 times higher than during cooling. k, estimated using near-surface ε, was ∼10 cm hr−1, validated with k obtained from chamber measurements, and 2–5 times higher than computed from wind-based models. The flux Richardson number (Rf) varied from ∼0.4 to ∼0.001 with a median value of 0.04 in the upper 0.25 m, less than the critical value of 0.2. We extend MOST by incorporating the variability in Rf when scaling the influence of β+ relative to u*w 3/κz in estimates of ε, and by extension, k, obtained from time series meteorological and temperature data. 
    more » « less
  2. Abstract

    Seasonally flooded forests along tropical rivers cover extensive areas, yet the processes driving air‐water exchanges of radiatively active gases are uncertain. To quantify the controls on gas transfer velocities, we combined measurements of water‐column temperature, meteorology in the forest and adjacent open water, turbulence with an acoustic Doppler velocimeter, gas concentrations, and fluxes with floating chambers. Under cooling, measured turbulence, quantified as the rate of dissipation of turbulent kinetic energy (ε), was similar to buoyancy flux computed from the surface energy budget, indicating convection dominated turbulence production. Under heating, turbulence was suppressed unless winds in the adjacent open water exceeded 1 m/s. Gas transfer velocities obtained from chamber measurements ranged from 1 to 5 cm/hr and were similar to or slightly less than predicted using a turbulence‐based surface renewal model computed with measuredεandεpredicted from wind and cooling.

     
    more » « less
  3. Abstract Daytime atmospheric boundary layer (ABL) dynamics—including potential temperature budgets, water vapour budgets, and entrainment rates—are presented from in situ flight data taken on six afternoons near Fresno in the San Joaquin Valley (SJV) of California during July/August 2016. The flights took place as a part of the California Baseline Ozone Transport Study aimed at investigating transport pathways of air entering the Central Valley from offshore and mixing down to the surface. Midday entrainment velocity estimates ranged from 0.8 to 5.4 cm s −1 and were derived from a combination of continuously determined ABL heights during each flight and model-derived subsidence rates, which averaged -2.0 cm s −1 in the flight region. A strong correlation was found between entrainment velocity (normalized by the convective velocity scale) and an inverse bulk ABL Richardson number, suggesting that wind shear at the ABL top plays a significant role in driving entrainment. Similarly, we found a strong correlation between the entrainment efficiency (the ratio of entrainment to surface heat fluxes with an average of 0.23 ± 0.15) and the wind speed at the ABL top. We explore the synoptic conditions that generate higher winds near the ABL top and propose that warm anomalies in the southern Sierra Nevada mountains promote increased entrainment. Additionally, a method is outlined to estimate turbulence kinetic energy, convective velocity scale ( w * ), and the surface sensible heat flux in the ABL from a slow, airborne wind measurement system using mixed-layer similarity theory. 
    more » « less
  4. Abstract Microstructure observations in the Pacific cold tongue reveal that turbulence often penetrates into the thermocline, producing hundreds of watts per square meter of downward heat transport during nighttime and early morning. However, virtually all observations of this deep-cycle turbulence (DCT) are from 0°, 140°W. Here, a hierarchy of ocean process simulations, including submesoscale-permitting regional models and turbulence-permitting large-eddy simulations (LES) embedded in a regional model, provide insight into mixing and DCT at and beyond 0°, 140°W. A regional hindcast quantifies the spatiotemporal variability of subsurface turbulent heat fluxes throughout the cold tongue from 1999 to 2016. Mean subsurface turbulent fluxes are strongest (∼100 W m −2 ) within 2° of the equator, slightly (∼10 W m −2 ) stronger in the northern than Southern Hemisphere throughout the cold tongue, and correlated with surface heat fluxes ( r 2 = 0.7). The seasonal cycle of the subsurface heat flux, which does not covary with the surface heat flux, ranges from 150 W m −2 near the equator to 30 and 10 W m −2 at 4°N and 4°S, respectively. Aseasonal variability of the subsurface heat flux is logarithmically distributed, covaries spatially with the time-mean flux, and is highlighted in 34-day LES of boreal autumn at 0° and 3°N, 140°W. Intense DCT occurs frequently above the undercurrent at 0° and intermittently at 3°N. Daily mean heat fluxes scale with the bulk vertical shear and the wind stress, which together explain ∼90% of the daily variance across both LES. Observational validation of the scaling at 0°, 140°W is encouraging, but observations beyond 0°, 140°W are needed to facilitate refinement of mixing parameterization in ocean models. Significance Statement This work is a fundamental contribution to a broad community effort to improve global long-range weather and climate forecast models used for seasonal to longer-term prediction. Much of the predictability on seasonal time scales is derived from the slow evolution of the upper eastern equatorial Pacific Ocean as it varies between El Niño and La Niña conditions. This study presents state-of-the-art high-resolution regional numerical simulations of ocean turbulence and mixing in the eastern equatorial Pacific. The results inform future planning for field work as well as future efforts to refine the representation of ocean mixing in global forecast models. 
    more » « less
  5. null (Ed.)
    Lakes and reservoirs contribute to regional carbon budgets via significant emissions of climate forcing trace gases. Here, for improved modelling, we use 8 years of floating chamber measurements from three small, shallow subarctic lakes (2010–2017, n = 1306) to separate the contribution of physical and biogeochemical processes to the turbulence-driven, diffusion-limited flux of methane (CH4) on daily to multi-year timescales. Correlative data include 9 years of surface water concentration measurements (2009–2017, n = 606), total water column storage (2009–2017, n = 1593) and in situ meteorological observations. We used the latter to compute near surface turbulence based on similarity scaling and then applied the surface renewal model to compute gas transfer velocities. Chamber fluxes averaged 6.9 ± 0.3 mg CH4 m−2 d−1 and gas transfer velocities (k600) averaged 4.0 ± 0.1 cm h−1. Chamber derived gas transfer velocities tracked the power-law wind speed relation of the model. Coefficients for the model and dissipation rates depended on shear production of turbulence, atmospheric stability, and exposure to wind. Fluxes increased with wind speed until daily average values exceeded 6.5 m s−1, at which point emissions were suppressed due to rapid water column degassing reducing the water–air concentration gradient. Arrhenius-type temperature functions of the CH4 flux (Ea‘ = 0.90 ± 0.14 eV) were robust (R2 ≥ 0.93, p < 0.01) and also applied to the surface CH4 concentration (Ea‘ = 0.88 ± 0.09 eV). These results imply that emissions were strongly coupled to production and supply to the water column. Spectral analysis indicated that on timescales shorter than a month, emissions were driven by wind shear whereas on longer timescales variations in water temperature governed the flux. Long-term monitoring efforts are essential to identify distinct functional relations that govern flux variability on timescales of weather and climate change. 
    more » « less