doi:10.1016/j.physd.2006.02.004
Copyright © 2006 Elsevier Ltd All rights reserved.
An analysis of intermittency, scaling, and surface renewal in atmospheric surface layer turbulence
Gabriel Katula,
,
, Amilcare Porporatob, Daniela Cavac and Mario Siqueiraa
aNicholas School of the Environment and Earth Sciences, Box 90328, Duke University, Durham, NC 27708-0328, USA
bDepartment of Civil and Environmental Engineering, Box 90287 Hudson Hall, Durham, NC 27708-0287, USA
cCNR-Institute of Atmosphere Sciences and Climate Section of Lecce, Lecce, Italy
Received 22 May 2005;
revised 31 January 2006;
accepted 7 February 2006.
Communicated by C.K.R.T. Jones.
Available online 24 March 2006.
References and further reading may be available for this article. To view references and further reading you must
purchase this article.
Abstract
Turbulent velocity and scalar concentration time series were collected in the atmosphere above an ice sheet, a mesic grassland, and a temperate pine forest, thereby encompassing a wide range of roughness conditions encountered in nature. Intermittency and scaling properties of such series were then analyzed using Tsallis’s non-extensive thermostatistics. While theoretical links between the Tsallis’s non-extensive thermostatistics and Navier–Stokes turbulence remain questionable, the Tsallis distribution (a special interpretation of Student’s t-distribution) provides a unifying framework to investigate two inter-connected problems: similarity between scalars and velocity statistics within the inertial subrange and “contamination” of internal intermittency by “external” factors. In particular, we show that “internal” intermittency models, including the She–Leveque, Lognormal, and Log-stable, reproduce the observed Tsallis parameters well for velocities within the inertial subrange, despite the differences in surface roughness conditions, but fail to describe the fluctuations for the scalars (e.g., air temperature CO2 and water vapor). Such scalars appear more intermittent than velocity when the underlying surface is a large source or sink. The dissimilarity in statistics between velocity and scalars within the inertial subrange is shown to be strongly dependent on “external” intermittency. The genesis of “external” intermittency for scalars is linked to the classical Higbie surface renewal process and scalar source strength. Surface renewal leads to a ramp-like pattern in the scalar concentration (or temperature) time series with a gradual increase (rise-phase) associated with sweeping motion from the atmosphere onto the surface or into the canopy and a sharp drop associated with an ejection phase from the surface (or the canopy) back into the atmosphere. The duration of the rise-phase is on the order of the integral time scale, while the duration of the ejection phase is much shorter and is shown to impact the distributional tails at the small scales. Implications for “scalar turbulence” models are also discussed in the context of biosphere–atmosphere CO2 exchange.
Keywords: Antartica; Grassland; Pine forest; Atmosphere turbulence; Intermittency; Scalar transfer; Surface renewal; Tsallis statistics
Fig. 1. Comparison between measured (open circle) and modeled (dots) p(x) for separation distances r/zm=0.059, 0.12, 0.2, 0.4, 0.98, 1.96, 4.9 and for the longitudinal (u) (left) and vertical (w) (middle) velocities, and air temperature (T) differences at the grass site. Note that for each flow variable, p(x) is shifted upwards by four decades with increasing r/zm, with the lowest corresponding to r/zm=0.059. Also, for each r/zm, the 12 separate lines and data correspond to 12 runs collected under similar meteorological and mean wind conditions, each 1/2 h in duration.
Fig. 2a. Comparison between measured and modeled Tsallis statistics for the Antarctica velocity data using the models in Table 1 for ξn. The solid horizontal line is for K41.
Fig. 2b. The same as Fig. 2a, but for the grass site.
Fig. 3. Time series of the three velocity components, air temperature, CO2, and water vapor above a pine forest for near-neutral conditions. All variables are normalized to have zero mean and unit variance. The run is collected between 1800 and 1830 when photosynthetically active radiation (PAR)
0.
Fig. 4. The same as Fig. 2a, but for the Pine forest data in Fig. 3.
Fig. 5. Mechanism for the onset of ramp-like motion in the CO2 time series and its interpretation using surface renewal theory. The inset is a 200 s time series taken from Fig. 3 (starting at
) documenting the ramp-like pattern of the CO2 concentration. The mean residence time (or ramp duration) is 60 s.
Fig. 6a. A ramp-like series with a constant residence time of 60 s and a fractional Brownian motion (fBm) series with a Hurst exponent =1/3. The fBm series recovers K41 scaling (no internal or external intermittency).
Fig. 6b. Left: Composite time series constructed from linear combinations of ramp-like
and fBm
signals using
. Right: The computed K(r). The measured K(r) for CO2 is taken from Fig. 4 (circles) and is repeated for reference. A decreasing α corresponds to an increase in “external” intermittency produced by the surface renewal process.
Fig. 7. Comparison between measured and modeled p(x) for u and CO2 series in Fig. 3.
Table 1.
Models for the “internal” intermittency exponent ξn for
within the inertial subrange, where n is the order of the moment

Table 2.
Summary of the experimental setup at the three sites
