TMP 31 AB
TMP 31 AB
BOUNDARY-LAYER TURBULENCE
GERMAN POVEDA-JARAMILLO
Postgrado en Aprovechamiento de Recursos Hidrdulicos, Universidad Nacional de Colombia, Seccional
de Medellin, Apartado Adreo 1027, Medellln, Colombia
and
C A R L O S E. P U E N T E
Hydrologic Science, Department of Land, Air and Water Resources, University of California, Davis,
CA 95616, U.S.A.
Abstract. The possible chaotic nature of the turbulence of the atmospheric boundary layer in and
above a decidious forest is investigated. In particular, this work considers high resolution temperature
and three-dimensional wind speed measurements, gathered at six alternative elevations at Camp
Borden, Ontario, Canada (Shaw et al., 1988). The goal is to determine whether these time series may
be described (individually) by sets of deterministic nonlinear differential equations, such that: (i) the
data's intrinsic (and seemingly random) irregularities are captured by suitable low-dimensional fractal
sets (strange attractors), and (ii) the equation's lack of knowledge of initial conditions translates into
unpredictable behavior (chaos). Analysis indicates that indeed all series exhibit chaotic behavior, with
strange attractors whose (correlation) dimensions range from 4 to 7. These results support the existence
of a low-dimensional chaotic attractor in the lower atmosphere.
1. Introduction
Turbulence has become one of the most elusive problems in physics during the
last century. Several approaches, which include those by Von Karman, Kolmogo-
rov and Landau, have failed to simultaneously describe, explain and predict turbu-
lent motion (Stanigid, 1985; Ruelle, 1991). New approaches to the problem have
recently emerged due to the advent of chaos theory, in particular through the
characterization of "strange attractors" (Ruelle and Takens, 1971). The idea of
chaos has been lurking since Poincar6's work on the stability of nonlinear differ-
ential equations. Recently, developments in computer technology and mathematics
have allowed an improved understanding of nonlinear chaotic phenomena, leading
to better visualization of its solutions and the identification of strange attractors.
The first of such attractors was observed in a set of nonlinear deterministic equa-
tions (a truncated form of the Navier-Stokes equations) in an attempt to model
fluid thermal convection in a gradient-perturbed atmosphere (Lorenz, 1963). This
work made clear the intrinsic differences that exist between linear and nonlinear
systems, and showed that complex and seemingly random behavior may be ob-
tained from nonlinear deterministic equations.
Nonlinear differential equations may be studied following the evolution of arbi-
1 . 5 L q , i i , , , i , , i , , , i , , , i , , , i , ' ~ i , , , i , , , i , , ,
1,0
0.5
(.3
o -0.5
-1.0
-1.5 ~
--2.0 , , , I , , , I , , I , , , I p , , I , , , J , ~ ~ I , , , I , , ~ I , , i
2.0 -1.8 -1.6 -1.4 -1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0
LOG (~)
Fig. 1. Powerspectrum of vertical wind velocityat site 6, and w-5/3 line.
trary initial conditions. For fixed times, the resulting coordinates define a point in
the phase-space of the system. As time passes, they delineate a trajectory (orbit)
in phase-space. The simplest of all behaviors results when all trajectories converge
to a single fixed point. In this case, all trajectories are "attracted" to the single
steady-state, irrespective of initial conditions. When the equations under study
are contractile, i.e., when volumes of arbitrary initial sets decrease as time passes,
more sophisticated behaviors may be observed (Lichtenberg and Lieberman,
1983). For these systems, termed dissipative, trajectories in phase-space may
converge to infinite (non-periodic) but bounded sets called "strange attractors"
(Berg6 e t a l . , 1984). These objects derive their names from their fractal nature
(Mandelbrot, 1982): they possess a built-in repetitive.structure that makes them
appear similar when viewed at different scales. When strange attractors exist,
extreme sensitivity to initial conditions is typically observed (Moon, 1987). Any
two trajectories which start close to one another in phase space, move exponen-
tially away from each other (for small times on the average). This extreme sensitiv-
ity to initial conditions defines chaotic motion.
Many types of turbulent flows have been shown to possess chaotic behavior,
with a strange attractor occurring on a low-dimensional phase-space (with few
coordinates). Among them there are: (i) Couette-Taylor flow (Brandstater et al.,
S T R A N G E A T T R A C T O R S IN A T M O S P H E R I C B O U N D A R Y - L A Y E R TURBULENCE 177
-2
-.3
-4
-5
"U
-7
-8
Fig. 2(a).
1983; Mullin and Price, 1989), (ii) Rayleigti-Benard convection (Libchaber, 1987),
(iii) fully developed turbulent flows in pipes (Sieber, 1987), (iv) meteorological
signals (Romanelli et al., 1988), (v) wind velocities (Tsonis and Elsner, 1989; Pool,
1989), and (vi) rainfall (Rodriguez-Iturbe et al., 1989), among others. All of these
strange attractors, living in Euclidean spaces of dimensions less than 10, have been
discovered using a powerful theorem proved by Takens (1981). The basis for the
procedure is the fact that a single signal (along any given coordinate) contains
information of all the variables involved in a given phenomenon liable to be
described by a set of coupled nonlinear differential equations. The idea is to
reconstruct the underlying attractor by means of successive shifts (delays) of a
single observable (signal), such that the obtained attractor preserves the topolog-
ical characteristics of the underlying process. At this stage, quantification of the
fundamental dynamics may be achieved.
The purpose of this work is to study the existence of strange attractors in the
turbulence field of the atmospheric boundary layer. For this purpose, the present
work reports on an analysis performed on high resolution time series of wind
velocities and temperature fluctuations, as measured in and above a deciduous
forest in Ontario, Canada (Shaw et al., 1988).
178 GERM,A~N P O V E D A - J A R A M I L L O AND CARLOS E. PUENTE
-2
--4
C" - 5
0
9
0
m
--6
-7
-8 -- -, ~ J r --
-g i i i i i i i ~ J [ r i i r p ~ t i i [ r i i I : i I t I I E I i I I i i i , I I , i i i i r , i
-4 -3 -2 -1 0
LOG r
Fig. 2(b).
2.1. C H O O S I N G T H E T I M E D E L A Y 7
Accurate definition of the delay z is crucial for the proper reconstruction of strange
attractors. This work uses the following optimality criteria due to Albano et al.
(1988, 1991). The time delay z is defined according to
S T R A N G E A T T R A C T O R S IN A T M O S P H E R I C B O U N D A R Y - L A Y E R TURBULENCE 179
-2
-5
-4
~, -5
,s
(_b
o
-6
-8
-9
-4 -5 -2 -1 0 1
LOG r
Fig. 2(c).
Wl
'r - (1)
N-I'
where wI is either the window length of the autocorrelation function of the data
(the first lag for which the autocorrelation first passes through zero) or the first
local minimum of the autocorrelation function, and N is the plausible embedding
dimension (the number of coordinates) being considered.
Since the motion of a chaotic signal is non-periodic, its power spectrum (the
Fourier transform of its autocorrelation function) does not exhibit defined peaks.
Instead, it typically yields a broad-band function, which in general shows an
exponential decay with frequency (Berg6 et al., 1984; Mayer-Kress, 1985; Sigeti
and Horsthemke, 1987), as opposed to an algebraic decay as observed for stochas-
tic processes (Greenside et al., 1982).
180 GERM~.N POVEDA-JARAMILLO AND CARLOS E, PUENTE
-6I
1
iI
--91 i I ~ i i i i i ~ I i ~ i p i F i i ] i I p I I i I i ~ I L i i I I I i I I I i i I r i i I I i--
a -3 -2 -1 0
LOG r
Fig. 2(d).
In s pq
1 i=1
Dq - lim - - (2)
q- 1,~ lne
STRANQE ATTRACTORS IN ATMOSPHERIC BOUNDARY-LAYER TURBULENCE 181
-2
-3
qp
(.9
0
-7
-8
m
-9L
-3 -2 -1 0 !
-4
LOG r
Fig. 2(e).
Notice that different q's weigh the attractor's concentration of points with varying
degrees and that Do simply yields the fractal dimension of the attractor.
In addition to the fractal dimension, the information dimension (D1) and the
correlation dimension (De) are obtained in the limit of Equation (2) as q tends to
1 and 2, respectively. D1 may be expressed as (Farmer et al., 1983):
m(~)
In Z Pi lnpi
D1 = lim i=1 (3)
~--,0 In 9
-2
~- -5
(D
L)
C)
12
-7-
/
-8- / _J ._Jj i -
-9 t ~ i i i i i i i ] T i i i i i i i i I i i I i i q ~ t r I i i i i i I r i r I i i i r i i E I i
-4 -3 -2 -1 0
LOG r
Fig. 2(f).
F i g . 2. Two-point correlation functions for selected time series. (a) w at site 7; (b) v at site 6; (c) w
at site 5; (d) v at site 4; (e) T at site 3; (f) u at site 2.
where 0 is the Heaviside step function, and M is the number of points in the
attractor. C(r) counts the number of pairs on the attractor ()(i,)?j), whose
distance 1 2 i - J(jl is less than r. Intuitively, C(r) accounts for the geometrical
correlation that points have within the attractor.
Notice that Equation (4) implies that for small values of r, the two-point corre-
lation function C(r) behaves like a power law:
C(r) ~ r D2 . (6)
14
12 0
10 D
o D
4
I
[]
0 I , r , l t , , I 1El , I
0 2 6 8 10 12 14
N
Fig. 3(a).
L.
L
12j- o
10
O X ~
0 r I i i I 1 ~ 1 I ~11 rll ,
0 2 6 8 10 12 14
N
Fig. 3(b)
184 GERMAN POVEDA-JARAM1LLO AND CARLOS E, PUENTE
q
i
12
r i
10
I
[]
8
L
o
L
(/3
iI"
i-
0 2 4 6 8 !0 12
N
Zig. 3(~).
12 D
[]
[]
10
0
o [3
[]
[]
4 []
0 2 4 6 8 10 12 14
N
Fig. 3(d).
S T R A N G E A T F R A C T O R S IN A T M O S P H E R I C B O U N D A R Y - L A Y E R TURBULENCE 185
14
12 []
10 []
[]
[]
&
0 []
[]
[]
0 2 4 6 8 10 12
N
Fig. 3(e).
!4
12 L- []
[3
10 O
[]
|
o_
o E1
[]
[]
4
i
0 2 4 6 8 10 12
N
Fig. 3(f).
Fig. 3. Slope D2 vs. embeddingdimensionN for selected series. Squares represent a purely random
signal. (a) w at site 7; (b) v at site 6; (e) w at site 5; (d) v at site 4; (e) T at site 3; (f) u at site 2.
2.5 ' ' ' I ' ' ' ' ' '
2.5 '
2.0
r'''l ~
I
i
i
I
''t
t
2511......1
#
2.0 ~i u
2.01 ~1
i ,
vI
1,5 u J" 1.5 1.5 w I
.)
,' s
1.0 * /)' 1.0
9N \
\\
\ \,~\
The desired computational savings are realized by considering some - but not all -
of the distances greater than ro. Theiler's method will be used in the experimental
part of this work.
As mentioned earlier, data requirements represent another constraint in the
estimation of the correlation dimension, D2. Several investigators have proposed
bounds on the actual number of data points needed. Among others, Smith (1988)
suggests a safe size of the order of 42 [D~ data points, where [Do] is the integer
part of the fractal dimension of the set; Ruelle (1990) proposes a size of 10 D2/2
points; while Essex and Nurenberg (1991) argue that 2D2(D2 + 1) points are needed
in the reconstructed attractor. These bounds will be checked later with meteoro-
logical data.
2.2.4. L y a p u n o v e x p o n e n t s
The Lyapunov exponents contain information about the speed at which trajectories
exponentially diverge or converge within the attractor. They are related to the
expanding or contracting directions in phase-space (Wolf et al., 1986). A chaotic
system is characterized by an exponential divergence of nearby initial conditions,
and therefore possesses at least one positive Lyapunov exponent. The magnitude
of the largest positive exponent (if existent) reflects the time scale on which the
188 GERMAN POVEDA-JARAMILLO AND CARLOS E. PUENTE
o.o ~ I A
-1.5 it
0.0
(.2
o
-I.5
c
o
z=] .Sin
ID
0.0 I*
23
4-F
(.9
-1.5
1.5
03
k_
3
-4J
-1.5
~D 1.5 F f=10.Sm ~1
L.
00
o_
E
0o
1.5 ~z
o.o
Time (s)
Fig. 5. Temperature fluctuation series for sites 1, 2, 3, 5, 6, 7, after Gao et al. (1989).
1 r Pi(t) 3
Ai = lira logz] |. (7)
~-~ t LP~(O)J
Several methods have been proposed to estimate these exponents from experi-
mental time series. Among them, Wolf et al. (1985) presented the first procedure
which follows the dynamic separation between nearby pairs of points on the
attractor. Landa and Chetverikov (1988) have suggested improvements to this
algorithm to shorten its computational time and decrease data requirements. Such
STRANGE ATFRACTORS IN ATMOSPHERIC BOUNDARY-LAYER TURBULENCE 189
0.220 ! I
0.207
0.i 93 - -
L
0.!53,
0.140 I I I
0 6 12 18 24 30
Time (Min)
Fig. 6(a).
improvements will be used together with the procedure of Wolf et al. (1985) to
estimate the largest Lyapunov exponent for the meteorological time series.
Basically, the method employed calculates the average divergence of initially
close orbits on the attractor, by considering many points along a trajectory. The
algorithm starts by selecting both a reference trajectory on the attractor and a
close point on a nearby trajectory. Then, it measures the rate of divergence as
both trajectories evolve. As the distance between the two trajectories becomes
too large, the algorithm looks for a new nearby trajectory and defines a new
small initial separation. Then, averaging over different regions of the attractor is
performed, and the largest Lyapunov exponent is estimated through the equation:
N
1
"~1 = - - E log2 P~(tk) (8)
t N -- to k=l Pi(tk-1)
3. Experimental Analysis
The mathematical tools outlined in the preceding sections were used to test the
existence of strange attractors (chaotic dynamics) in the turbulence of the atmo-
spheric boundary layer. Simultaneous measurements of wind velocity and tempera-
ture fluctuations at six locations constituted the data set of 24 meteorological time
190 GERM~dq P O V E D A - J A R A M I L L O AND CARLOS E, P U E N T E
0.180 I I I
0.163
i
i
I
0.147 I
I
~-0.130 ~
0.113
i
0.097 -
/
0.080 I I ~ J r
0 6 12 18 24 30
Time (Min)
Fig. 6(b).
series. The original data were measured at a deciduous forest experimental site in
Ontario, Canada, in an attempt to quantify the turbulent patterns observed in and
above the forest (Shaw et al., 1988, 1989; Gao et al., 1989, 1992). Such data
included: (i) longitudinal wind velocity (u), (ii) transversal wind velocity (v), (iii)
vertical wind velocity (w), and (iv) temperature fluctuations (T).
The set of meteorological time series was gathered using arrays of 3-D ane-
mometers located on two towers separated by a distance of 25 m in which 7
ultrasonic anemometers/thermometers were installed. The seven recording sites
were distributed in the following way: sites 1, 2, 3, and 4 at tower number one at
heights of 5.9, 10.5, 15.4 and 17.6 m above the ground, and sites 5, 6, and 7 at
tower number two at 17.9, 34.2, and 43.1 m above the ground. See Shaw et aL
(1988).
The simultaneous series considered in this work each consisted of 18,000 obser-
vations measured every 0.1s (10Hz) for a total length of 30min. They were
obtained on October 2, 1986 from 2 to 2:30 PM. Skies were cloudless. The Monin-
Obukhov length was computed at the forest height of 18 m, and its value of - 138
corroborated the presence of unstable atmospheric conditions.
S T R A N G E A T T R A C T O R S IN A T M O S P H E R I C B O U N D A R Y - L A Y E R T U R B U L E N C E 191
0.46 I I I I
0.44
0.42
~<- 0.40
0.38
0.36 -
0.34 I I I I
0 6 12 18 24 50
Time (Min)
Fig. 6(c).
We shall use the 24 (4 x 6) time series recorded at the six highest sites, i.e., 2
through 7. For further information on the time series used, see Poveda (1989).
0.SO01
i
i
0.27.5 "-
I
o.247
i
I
.<- 0.220 I
ii i
/
I
,J.ls5 "-
I
m
o.167L
i /
/
0.!40[_ I I i L
0 6 12 18 24 5O
Time (Min)
Fig. 6(d).
Fig. 6~ Time evolution of largest Lyapunov exponent for selected series. (a) u at site 6; (b) T at site
6; (c) u at site 5; (d) v at site 5.
space ranged from 2 through 12. For all cases, convergence of the slopes of
log C(r) vs. log r (found via linear regression) was obtained. This could be seen
in Figures 2a-f, which in order show the two-point correlation functions for the
series w7, v6, w5, v4, T3, and u2, as the embedding dimension N is varied. Figures
3a-f present the embedding dimension, N, vs. correlation dimension, D2, for the
time series considered in Figures 2. Also included in Figures 3, as squares, are
the correlation dimensions that would be obtained from purely stochastic pro-
cesses. Notice that the computed dimensions saturate as N increases. This fact
indicates the existence of low-dimensional attractors, and deterministic chaos,
whose confirmation requires positive Lyapunov exponents.
Table I includes the correlation dimensions found for the different phenomena.
The orders of magnitude obtained in this work agree well with those previously
found by Tsonis and Elsner (1989) and Pool (1989) for wind velocities at other
locations. Although the very stringent data requirements set forward by Smith
(1988) are not valid for any of the time series analyzed, both conditions suggested
by Ruetle (1990), and by Essex and Nurenberg (1991) are satisfied by all data
sets. Observe that Kolmogorov's - 5 / 3 law implies for a stochastic process with
power law spectra a correlation dimension of 3. Because all correlation dimensions
STRANGE ATrRACTORS IN A T M O S P H E R I C B O U N D A R Y - L A Y E R TURBULENCE 193
TABLEI
Correlationdimensionsforsites2to7
Site 2 3 4 5 6 7
phenomenon
u 4.4 3.5 5.0 4.1 5.7 6.4
v 4.0 4.4 3.7 4.9 5.7 6.2
w 4.3 4.1 3.6 4.1 6.7 6.8
T 4.5 5.0 5.3 5.8 6.7 5.5
in Table I are above 3, there is indeed no evidence that such data may be
represented by such a stochastic process.
Notice that although sites 4 and 5 are located at about the same elevation on
separate towers, their correlation dimensions were not equal. This observation,
based on a single data set, does not contradict the good agreement in Reynolds
stresses that Shaw et al. (1988) reported over a variety of stability conditions at
those elevations. Even though there are no obvious trends on all the obtained
dimensions, it is worth noticing that correlation dimensions for velocities do in-
crease for elevations at and above the forest canopy (sites 4, 6, 7, or, 5, 6, 7),
indicating qualitatively more turbulent behavior for higher elevations. Inside the
forest at the first tower (2, 3, 4), there are no consistent patterns. As seen in
Figure 4a, there is non-monotonic behavior in correlation dimensions for velocities
in the horizontal plane, while they decrease with height for vertical velocities.
This overall picture resembles the patterns observed by Shaw et al. (1988) regard-
ing horizontal correlations (uw) vs. height, with higher negative correlations above
the forest, and variable negative correlations inside the forest, see Figure 4b. Also,
mirror images of the observed profiles on correlation dimensions exhibit the same
qualitative trends of relative intensities of turbulence (a coefficient of variation
defined by dividing the velocity standard deviations by the cup wind speed), as
reported by Shaw et al. (1988), Figure 4c. Whether or not these structures could
be understood jointly remains an open question. Certainly, a chaotic analysis of
such surrogates of turbulence (and others like Reynolds stresses, heat fluxes, etc.)
may provide additional information, specially if alternative stability conditions for
different foliage coverages are considered.
Figure 5 shows the observed ramp structure of the temperature measurements
at sites 1, 2, 3, 5, 6, and 7, as reported by Gao et al. (1989). The ramps are
more clearly defined for locations inside the forest (below 18 m). The relative
unpredictability of temperature fluctuations for different elevations is nicely cap-
tured by the correlation dimensions in the last row of Table I. Observe that as
one moves upward, the correlation dimensions increase (except for the highest
site), indicating more erratic behavior and the need for higher order sets of
differential equations. Although the pattern is broken for the last site, notice that
the data complexity captured by the eye is entirely consistent with the obtained
194 GERMAN POVEDA-JARAMILLO AND CARLOS E. PUENTE
TABLE II
Lyapunovexponents at sites 5, 6, and 7
Site 5 6 7
phenomenon
u 0.370 0.205 0.195
v 0.250 0.250 0.150
w 0.185 0.273 0.209
T 0.211 0.156 0.155
Acknowledgments
.,<
Field data used in this work were from an experiment funded by Environment
Canada under the direction of Drs. G,. den H artog and H. H. Neumann of the
Atmospheric Environment Service. Th~ authors.thank Dr. Roger Shaw for provid-
ing the data, Dr. James Theiler for sharing his computer algorithm, Dr. Alfonso
Albano for sending his unpublished manuscripts, and Mr. Miguel Mejia for his
help on the final graphs of this work. This research was partially funded by
California Water Resources Center W-731 and by Kearney Foundation (Project
86-7). The suggestions of anonymous reviewers are gratefully acknowledged.
References
Abarbanel, H. D. I., Brown, R., and Kadtke, J. B.: 1990, 'Prediction in Chaotic Nonlinear Systems-
Method for Time Series with Broadband Fourier Spectra', Phys. Rev. A. 41, 1782-1807.
Albano, A. M., Muench, J., and Schwartz, C.: 1988, 'Singular-Value Decomposition and the
Grassberger-Procaccia Algorithm', Phys. Rev. A. 38, 3017-3025.
Albano, A. M., Passamante, A., and Farrell, M. E.: 1991, 'Using Higher-Order Correlations to Define
an Embedding Window', Physica D 54, 85-97.
Bath, M.: 1974, Spectral Analysis in Geophysics, Elsevier, Amsterdam.
Berg6, P., Pomeau Y., and Vidal, C.: 1984, Order Within Chaos, John Wiley, New York.
Brandstater, A., Swift, J., Sweeney, H. L., and Wolf, A.: 1983, 'Low Dimensional Chaos in a
Hydrodynamic System', Phys. Rev. Lett. 51, 1442-1444.
Crutchfield, J. P. and McNamara, B. S.: 1987, 'Equations of Motion from Data Series', Complex Sys.
1, 417-452.
Essex, C., Lookman, T., and Nerenberg, M. A. H.: 1987, 'The Climate Attractor Over Short Time
Scales', Nature 326, 64-67.
Essex, C. and Nerenberg, M. A. H.: 1991, 'Comment on "Deterministic Chaos: The Science and the
Fiction" by D. Ruelle', Proc. R. Soc. Lond. A, 435, 287-292.
Farmer, J. D., Ott, E., and Yorke, J. A.: 1983, 'The Dimension of Chaotic Attractors', Physica D 7,
153-180.
Fraedrich, K.: 1986, 'Estimating the Dimensions of Weather and Climatic Attractors', J. Atmos. Sci.
43, 419-432.
Gao, W., Shaw, R. H., and Paw U, K. T.: 1989, 'Observation of Organized Structure in Turbulent
Flow Within and Above a Forest Canopy', Boundary-Layer Meteorol. 47, 349-377.
Gao, W., Shaw, R. H., and Paw U, K. T.: 1992, 'Conditional Analysis of Temperature and Humidity
196 GER/vLX.N P O V E D A - J A R A M I L L O AND CARLOS E. PUENTE
Microfronts and Ejection/Sweep Motions Within and Above a Deciduous Forest', Boundary-Layer
MeteoroL 59, 35-57.
Grassberger, P. and Procaccia, I.: 1983, 'Measuring the Strangeness of Strange Attractors'. Physica
D 9, 189-208,
Grassberger, P. and Procaccia, I.: 1984, ~Dimensions and Entropies of Strange Attractors from a
Fluctuating Dynamics Approach', Physica D 13, 34-54.
Grassberger, P.: 1986, 'Estimating the Fractal Dimensions and Entropies of Strange Attractors', in A.
V. H01den (ed.), Chaos, Princeton University Press.
Grassberger, P.: 1986, 'Do Climatic Attractors Exist?', Nature 323, 609-612.
Greenside, H. S., Wolf I., Swift J. and Pignataro, T.: 1982. qmpracticahty of a Box-Counting Algo-
rithm for Calculating the Dimensionality of Strange Attractors', Phys. Rev. A. 25, 3453-3446.
Hentschel, H. G. E. and Procaccia, I.: 1983, 'The Infinite Number of Generalized Dimensions of
Fractals and Strange Attractors', Physica D 8, 435-444.
Landa, P. S., and Chetverikov, V. I.: 1988, 'On the Evaluation of the Maximum Lyapunov Exponent
from a Single Experimental Time Series', Soviet Phys. Tech. Phys. 33,236-268.
Libchaber, A.: 1987, 'From Chaos to Turbulence in Benard Convection', Proc. Royal Soc. London.
A 413, 63-69.
Lichtenberg, A. J. and Lieberman, M. A.: 1983, Regular and Stochastic Motion, Springer-Verlag, New
York: 63-69, 1987.
Lorenz, E. N.: 1963, 'Deterministic nonperiodic Flow;, J. Atmos. Sci. 20, 130-141.
Lorenz, E. N.: 1991, 'Dimension of Weather and Climatic Attractors', Nature 353, 241-244.
MacCready, J.: 1962, ~Turbulence Measurements by Sailplane', J. Geophys. Res. 67, 1041-1050.
Mandelbrot, B. B.: 1982, The Fractal Geometry of Nature, Freeman and Co., San Francisco.
Mayer-Kress, Ed.: 1985, Dimensions and Entropies in Chaotic Systems, Springer Verlag, Berlin.
Moller, M., Lange, W., Mitschke, F., and Abraham, N. B.: 1989 'Errors from Digitizing and Noise
in Estimating Attractors Dimensions', Physics LettersA 138, 176-182.
MullJn, T., and Price, T, J.: 1989 'An Experimental Observation of Chaos Arising from the Interaction
of Steady and qglme-Dependent Flows', Nature 340, 294-296.
Nicolis, C. and Nicolis, G.: 1984, 'Is There a Climatic Attractor?', Nature 311, 529-532.
Osborne, A, R. and Provenzate, A.: 1989, 'Finite Correlation Dimension for Stochastic Systems with
Power-Law Spectra', Physica D 35, 357-381.
Pincus, S. M.: 1991, 'Approximate Entropy as a Measure of System Complexity', Proc. Nat. Acad.
Sci. USA, 88, 2297-2301.
Pool, R~: 1989, 'Is Something Strange about the Weather?', Science 243, 1290-1293.
Poveda, G.: 1989, ~ for Chaos in Meteorological Time Series', M. S. Thesis, University of
California, Davis.
Press, W. H., Flannery, B. P, Teukolsky, S. A., and Vetterling, W. T.: 1986, Numerical Recipes: The
Art of Scientific Computing, Cambridge University Press, Cambridge.
Rodriguez-Iturbe, I., Febres de Power, B., Sharifi, M. B., and Georgakakos, K. P.: 1989, 'Chaos in
Rainfall', Water Res. Res. 25, 1667-1675.
Romanelli, L., Figliola, M. A., and Hirsch, F. A.: 1988, 'Deterministic Chaos and Natural Phenomena',
J. Star. Phys. 53, 991-994.
Ruelle, D., and Takens, F.: 1971, 'On the Nature of Turbulence', Comm. Theo. Phys. 20, 167-192.
Ruelle, D.: 1990, 'Deterministic Chaos: The Science and the Fiction', Proc. R. Soc. Lond. A., 427,
241-248.
Ruelle, D.: i991, Chance and Chaos, Princeton University Press.
Shaw, R. H., den Hartog, G., and Newmann, H. H.: 1988, 'Influence of Foliar Density and Thermal
Instability on Profiles of Reynolds Stress and Turbulence Intensity in a Deciduous Forest', Boun-
dary-Layer Meteorol. 45, 391-409.
Shaw, R. H., Paw U, K. T., and Gao, W.: 1989, 'Detection of Temperature Ramps and Flow
Structures at a Deciduous Forest Site', Agric. Forest Meteorol. 47, 123-138.
Sieber, M. : 1987, 'Experiments on the Attractor-Dimension for Turbulent Pipe Flow', Phys. Lett. A
122, 467-469.
Sigeti, D., and Horsthemke, W.: 1987, 'High-Power Spectra for Systems Subject to Noise', Phys. Rev.
A 35, 2276-2282.
S T R A N G E A T T R A C T O R S IN A T M O S P H E R I C B O U N D A R Y - L A Y E R T U R B U L E N C E 197
Smith, L. A.: 1988, 'Intrinsic Limits on Dimension Calculations', Phys. Lett. A 133(6), 283-288.
Stanigid, M. M.: 1985, The Mathematical Theory o f Turbulence, Springer-Verlag, New York.
Sugihara, G., and May, R. M.: 1990, 'Nonlinear Forecasting as a Way of Distinguishing Chaos Form
Measurement Error in Time Series', Nature 344, 734-741.
Takens, F.: 1981, 'Detecting Strange Attractors in Turbulence', In Lecture Notes in Mathematics 898,
Springer Verlag, New York.
Theiler, J.: 1986, 'Spurious Dimension from Correlation Algorithms Applied to'Limited Time-Series
Data', Phys. Rev. A 34, 2427-2432.
Theiler, J.: 1987, 'Efficient Algorithm for Estimating the Correlation Dimension from a Set of Discrete
Points', Phys. Rev. A 36, 4456-4462.
Tsonis, A. A. and Eisner, J. B.: 1989, 'Chaos, Strange Attractors and Weather', Bull. Amer. MeteoroI.
Soc. 70, 14-23.
Wolf, A., Swift, J. B., Swinney, H., and Vastano, J. A.: 1985, 'Determining Lyapunov Exponents
from a Time Series', Physica D 16, 286-317.
Wolf, A.: 1986, 'Quantifying Chaos with Lyapunov Exponents', in V. Holden (ed.), Chaos, Holden,
Manchester University Press.