skip to main content

Title: Scalewise Return to Isotropy in Stratified Boundary Layer Flows

Anisotropic turbulence is ubiquitous in atmospheric and oceanic boundary layers due to differences in energy injection mechanisms. Unlike mechanical production that injects energy in the streamwise velocity component, buoyancy affects only the vertical velocity component. This anisotropy in energy sources, quantified by the flux Richardson numberRif, is compensated by a “return to isotropy” (RTI) tendency of turbulent flows. Describing RTI in Reynolds‐averaged models and across scales continues to be a challenge in stratified turbulent flows. Using phenomenological models for spectral energy transfers, the necessary conditions for which the widely‐used Rotta model captures RTI across variousRifand eddy sizes are discussed for the first time. This work unravels adjustments to the Rotta constant, withRifand scale, necessary to obtain consistency between RTI models and the measured properties of the atmospheric surface layer for planar‐homogeneous and stationary flows in the absence of subsidence. A range ofRifand eddy sizes where the usage of a conventional Rotta model is prohibited is also found. Those adjustments lay the groundwork for new closure schemes.

more » « less
Award ID(s):
1644382 1754893
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Atmospheres
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Here we have determined the nature of turbulent flow associated with oceanic nonbreaking waves, which are on average much more prevalent than breaking waves in most wind conditions. We found this flow to be characterized by a low turbulence microscale Reynolds number of30 < Reλ < 100. We observed that the turbulent kinetic energy dissipation rate associated with nonbreaking wavesϵ, ranged to3 · 10−4 W/kg for a wave amplitude 50 cm. Theϵ, under nonbreaking waves, was consistent with;Sijis the large‐scale (energy‐containing scales) wave‐induced mean flow stress tensor. The turbulent Reynolds stress associated with nonbreaking waves was consistent with experimental data when parameterized by an amplitude independent constant turbulent eddy viscosity, 10 times larger than the molecular value. Given that nonbreaking waves typically cover a much larger fraction of the ocean surface (90–100%) than breaking waves, this result shows that their contribution to wave dissipation can be significant.

    more » « less
  2. Abstract

    Wetlands play an important role in regulating the atmospheric carbon dioxide (CO2) concentrations and thus affecting the climate. However, there is still lack of quantitative evaluation of such a role across different wetland types, especially at the global scale. Here, we conducted a meta‐analysis to compare ecosystemCO2fluxes among various types of wetlands using a global database compiled from the literature. This database consists of 143 site‐years of eddy covariance data from 22 inland wetland and 21 coastal wetland sites across the globe. Coastal wetlands had higher annual gross primary productivity (GPP), ecosystem respiration (Re), and net ecosystem productivity (NEP) than inland wetlands. On a per unit area basis, coastal wetlands provided largeCO2sinks, while inland wetlands provided smallCO2sinks or were nearlyCO2neutral. The annualCO2sink strength was 93.15 and 208.37 g C m−2for inland and coastal wetlands, respectively. AnnualCO2fluxes were mainly regulated by mean annual temperature (MAT) and mean annual precipitation (MAP). For coastal and inland wetlands combined,MATandMAPexplained 71%, 54%, and 57% of the variations inGPP,Re, andNEP, respectively. TheCO2fluxes of wetlands were also related to leaf area index (LAI). TheCO2fluxes also varied with water table depth (WTD), although the effects ofWTDwere not statistically significant.NEPwas jointly determined byGPPandRefor both inland and coastal wetlands. However, theNEP/ReandNEP/GPPratios exhibited little variability for inland wetlands and decreased for coastal wetlands with increasing latitude. The contrasting ofCO2fluxes between inland and coastal wetlands globally can improve our understanding of the roles of wetlands in the global C cycle. Our results also have implications for informing wetland management and climate change policymaking, for example, the efforts being made by international organizations and enterprises to restore coastal wetlands for enhancing blue carbon sinks.

    more » « less
  3. Abstract

    Accurate estimation of terrestrial gross primary productivity (GPP) remains a challenge despite its importance in the global carbon cycle. Chlorophyll fluorescence (ChlF) has been recently adopted to understand photosynthesis and its response to the environment, particularly with remote sensing data. However, it remains unclear how ChlF and photosynthesis are linked at different spatial scales across the growing season. We examined seasonal relationships between ChlF and photosynthesis at the leaf, canopy, and ecosystem scales and explored how leaf‐level ChlF was linked with canopy‐scale solar‐induced chlorophyll fluorescence (SIF) in a temperate deciduous forest at Harvard Forest, Massachusetts,USA. Our results show that ChlF captured the seasonal variations of photosynthesis with significant linear relationships between ChlF and photosynthesis across the growing season over different spatial scales (R= 0.73, 0.77, and 0.86 at leaf, canopy, and satellite scales, respectively;P < 0.0001). We developed a model to estimateGPPfrom the tower‐based measurement ofSIFand leaf‐level ChlF parameters. The estimation ofGPPfrom this model agreed well with flux tower observations ofGPP(R= 0.68;P < 0.0001), demonstrating the potential ofSIFfor modelingGPP. At the leaf scale, we found that leafFq/Fm, the fraction of absorbed photons that are used for photochemistry for a light‐adapted measurement from a pulse amplitude modulation fluorometer, was the best leaf fluorescence parameter to correlate with canopySIFyield (SIF/APAR,R= 0.79;P < 0.0001). We also found that canopySIFandSIF‐derivedGPP(GPPSIF) were strongly correlated to leaf‐level biochemistry and canopy structure, including chlorophyll content (R= 0.65 for canopyGPPSIFand chlorophyll content;P < 0.0001), leaf area index (LAI) (R= 0.35 for canopyGPPSIFandLAI;P < 0.0001), and normalized difference vegetation index (NDVI) (R= 0.36 for canopyGPPSIFandNDVI;P < 0.0001). Our results suggest that ChlF can be a powerful tool to track photosynthetic rates at leaf, canopy, and ecosystem scales.

    more » « less
  4. Abstract

    The flux Richardson numberRf, also called the mixing efficiency of stratified turbulence, is important in determining geophysical flow phenomena such as ocean circulation and air‐sea transports. MeasuringRfin the field is usually difficult, thus parameterization ofRfbased on readily observed properties is essential. Here, estimates ofRfin a strongly turbulent, sediment‐stratified estuarine flow are obtained from measurements of covariance‐derived turbulent buoyancy fluxes (B) and spectrally fitted values of the dissipation rate of turbulent kinetic energy (ε). We test scalings forRfin terms of the buoyancy Reynolds number (Reb), the gradient Richardson number (Ri), and turbulent Froude number (Frt). Neither theReb‐based nor theRi‐based scheme is able to describe the observed variations inRf, but theFrt‐based parameterization works well. These findings support further use of theFrt‐ based parameterization in turbulent oceanic and estuarine environments.

    more » « less
  5. Abstract

    Though the Kelvin‐Helmholtz instability (KHI) has been extensively observed in the mesosphere, where breaking gravity waves produce the conditions required for instability, little has been done to describe quantitatively this phenomenon in detail in the mesopause and lower thermosphere, which are associated with the long‐lived shears at the base of this statically stable region. Using trimethylaluminum (TMA) released from two sounding rockets launched on 26 January 2018, from Poker Flat Research Range in Alaska, the KHI was observed in great detail above 100 km. Two sets of rocket measurements, made 30 min apart, show strong winds (predominantly meridional and up to 150 ms−1) and large total shears (90 ms−1 km−1). The geomagnetic activity was low in the hours before the launches, confirming that the enhanced shears that triggered the KHI are not a result of the E‐region auroral jets. The four‐dimensional (three‐dimensional plus time) estimation of KHI billow features resulted in a wavelength, eddy diameter, and vertical length scale of 9.8, 5.2, and 3.8 km, respectively, centered at 102‐km altitude. The vertical and horizontal root‐mean‐square velocities measured 29.2 and 42.5 ms−1, respectively. Although the wind structure persisted, the KHI structure changed significantly with time over the interval separating the two launches, being present only in the first launch. The rapid dispersal of the TMA cloud in the instability region was evidence of enhanced turbulent mixing. The analysis of the Reynolds and Froude numbers (Re = 7.2 × 103andFr = 0.29, respectively) illustrates the presence of turbulence and weak stratification of the flow.

    more » « less