ROLAND B.

STULL
Atmospheric Science Programme, Department of Geography
The University of British Columbia, Vancouver, Canada.

An Introduction to

Boundary
Layer
Meteorology

KLUWER ACADEMIC PUBLISHERS

DORDRECHT / BOSTON / LONDON
Libfary 01 Congress Cataloging in Publication Data

SIUI I. RolinG S .. 1950-

An ,n'"OOuct lon to b. "n alrY loy," • • ' "O"OlogV ' Rollno e. 5'ul'.
ú K= CI. -- E^ r ç ë ú ú ç ê I Å= Hllne U " o," 'V!
IncluUS l D l f f K d ? ? l ú D ? = Ona I na ...
ISIIN,'3: 918-9().2n·2769-5
,. Bou nOlry I IY", ' [Link]'
OCe90.4.B655194 ,999
551.!I--ocI9 aa- 15564
ISBN-I): 978·90-2n-2769·5
001: IO.1007/978·94..(X)9·3027_8
c· ISBN·I): 978·94..((19..3027·8 '"

Published by Kluwe. Academic Publishe.s.

P.O. 80x 17. 3300 AA Dordred1t. The NetheflandS.

Kluwer Academic Publishers incorporates

the publishing proglammes 01
D. Reidel. Martinus Nijhotl. Dr W. Junk and MTP Press.

Sold and distributed in the U.S.A. and Canada

by Kluwer Academic: Publishers.
tOt Philip Orive. Norwell, MA 02061. U.S.A.

In all othel countries. sold and distributed

by Kluwer Academic Publ ishers GrOl)P.
P.O. 80x 322, 3.300 AH Dordlecht. The Netherlands.

e ú í =publislled 1999
RCIlO"inted 1989. 199 1. 1993. 1994
Reprinted with cITata 1991
R""rintoo 1999. 2001. 2003

Printtd on acid-fru paptr

All Rig hts Reserved

© t 988 by Kh./We, ACademic Publishers
SoftOOvef reprint of the hardoover tst edition t988
No part 01 the material protected by this copyright nOlice may be reproduced or
uhlized In any lorm or by any means. electronic or mechanical
Includin g photocopying. recording or by any inlormation "Iorage arid
retrieval system , without written permissioo Irom the copyright owner
 C!oo[co[s
Preface XI

1 Mean Boundary Layer Characteristics

1.1 A bowlCjary-layer definition 2
1.2 Wind and flow 3
1.3 Turbulent transport 4
1.4 Taylor's hypothesis 5
l.5 Virtual potential temperature 7
1.6 Boundary layer depth and slructure 9
1.7 Micrometeorology 19
1.8 Significance of the boundary layer 21
1.9 General references 23
1.10 References for this chapter 25
1.11 Exercises 26

2 Some Mathematical and Conceptual Tools: Part 1. Statistics 29

2.1 The significance of turbulence and its spectrum 29
2.2 The spectral gap 33
2.3 Mean and turbulent parts 33
2.4 Some basic statistical methods 35
2.5 Turbulence kinetic energy 45
2.6 Kinematic flux 47
2.7 Eddy flux 51
2.8 Summation notation 57
2.9 Stress 63
2.10 Friction velocity 67
2.11 References 68
2.12 Exercises 70
vi _ l r k a ^ o v i ^ v b o ú l o l i l d v =

3 Application or the Governing Equations to Turbulent Flow 75

3.1 Methodology 76
3.2 Basic governing equations 76
3.3 Simplifications, approximations, and scaling arguments 80
3.4 Equations for mean variables in a turbulent flow 87
3.5 Summary of equations, with simplifications 93
3.6 Case studies 97
3.7 References 110
3.8 Exercises 111

4 Prognostic Equations for Turbulent Fluxes and Variances 115

4.1 Prognostic equations for the turbulent departures 115
4.2 Free convection scaling variables 117
4.3 Prognostic equations for variances 120
4.4 Prognostic equations for turbulent fluxes 134
4.5 References 147
4.6 Exercises 148

5 Turbulence Kinetic Energy, Stability, and Scaling 151

5.1 The TKE budget derivation 151
5.2 Contributions 10 the TKE budget 153
5.3 TKE budget contributions as a function of eddy size 166
5.4 Mean kinetic energy and its interaction with turbulence 168
5.5 Stability concepts 169
5.6 The Richardson number 175
5.7 The Obukhov length 180
5.8 Dimensionless gradients 183
5.9 Miscellaneous scaling parameters 184
5.10 Combined stability tables 186
5.11 References 187
5.12 Exercises 189

6 Turbulence Closure Techniques 197

6.1 The closure problem 197
6.2 Parameterization rules 200
6.3 Local closure - zero and half order 202
6.4 Local closure - first order 203
6.5 Local closure - one-and·a·half order 214
6.6 Local closure - second order 220
6.7 Local closure - third order 224
6.8 Nonlocal closure - transilient turbulence theory 225
6.9 Nonlocal closure - spectral diffusivity theory 240
6.10 References 242
6.11 Exercises 245
 COmENTS vii

7 Boundary Conditions and External Forcings 251

7.1 Effective surface turbulent flux 251
7.2 Heat budget at the surface 253
7.3 Radiation budget 256
7.4 Auxes at interfaces 261
7.5 Partitioning of flux into sensible and latent portions 272
7.6 Aux 10 and from the ground 282
7.7 References 289
7.8 Exercises 292

8 Some Mathematical and Conceptual Tools: Part 2. Time Series 295

8.1 Time and space series 295
8.2 Autocorrelation 296
8.3 Structure function 300
8.4 Discrete Fourier transfonn 303
8.5 Fast Fourier Transfonn 310
8.6 Energy spectrum 312
8.7 Spectral characteristics 318
8.8 Spectra of two variables 329
8.9 Periodogram 335
8.10 Nonlocal spectra 336
8.11 Spectral decomposition of the TKE equation 340
8.12 References 342
8.13 Exercises 344

9 Similarity Theory 347

9.1 An overview 347
9.2 Buckingham Pi dimensional analysis methods 350
9.3 Scaling variables 354
9.4 Stable boundary layer similarity relationship lislS 360
9.5 Neutral boundary layer similarity relationship lislS 364
9.6 Convective boundary layer similarity relationship IislS 368
9.7 The log wind profile 376
9.8 Rossby-number similarity and proftle matching 386
9.9 Spectral similarity 389
9.10 Similarity scaling domains 394
9.11 References 395
9.12 Exercises 399
viii BOUNDARY LAYER MmTEOROLOGY

10 Measurement and Simulation Techniques 405

10.1 Sensor and measurement categories 405
10.2 Sensor lists 407
10.3 Active remote sensor observations of morphology 410
10.4 Instrument platforms 413
10.5 Field experiments 417
10.6 Simulation methods 420
10.7 Analysis methods 427
10.8 References 434
10.9 Exercises 438

11 Convective Mixed Layer 441

11.1 The unstable surface layer 442
11.2 The mixed layer 450
11.3 Enrrainment zone 473
11.4 Enrrainment velocity and its parameterization 477
11.5 Subsidence and advection 483
11.6 References 487
11.7 Exercises 494

12 Stable Boundary Layer 499

12.1 Mean Characteristics 499
12.2 Processc:; 506
12.3 Evolution 516
12.4 Other Depth Models 518
12 .5 Low-level (nocturnal) jet 520
12.6 Buoyancy (gravity) waves 526
12.7 Terrain slope and drainage winds 534
12.8 References 538
12.9 Exercises 542

13 Boundary Layer Clouds 545

13.1 Thermodynamics 545
13.2 Radiation 555
13.3 Cloud enrrainment mechanisms 558
13.4 Fair-weather cumulus 562
13.5 Stratocumulus 570
13 .6 Fog 576
13.7 References 578
13.8 Exercises 583
 CONTENTS ix

14 Geographic Effects 587

14.1 Geographically generated local winds 587
14.2 Geographically modified flow 595
14.3 Urban heat island 609
14.4 References 613
14.5 Exercises 617

Appendices 619
A. Scaling variables and dimensionless groups 621
B. Notation 629
C. Useful constants, parameters and conversion factors 639
D. Derivation of virtual potential temperature 645

Subject Index 649

Errata section 667

Preface

Part of the excitement of boundary-layer meteorology is the challenge in :;tudying and

understanding turbulent flow - one of the unsolved problems of classical physics.
Additional excitement stems from the rich diversity of topics and research methods that we
collect under the umbrella of boundary-layer meteorology. That we live our lives within
the boundary layer makes it a subject that touches us, and allows us to touch it. I've tried
to capture some of the excitement, challenges and diversity within this book.
I wrote this book with a variety of goals in mind. First and foremost, this book is
designed as a textbook. Fundamental concepts and mathematics are presented prior to
their use, physical interpretations of equation tenns are given, sample data is shown,
examples are solved, and exercises are included. Second, the book is organized as a
reference, with tables of parameterizations, procedures, field experiments, useful
constants, and graphs of various phenomena in a variety of conditions. Third, the last
several chapters are presented as a literature review of the current ideas and methods
in boundary layer meteorology.
It is assumed that this book will be used at the beginning graduate level for students
with an undergraduate background in meteorology. However, a diversity in the
background of the readers is anticipated. Those with a strong mathematical background
can skip portions of Chapter 2 on statistics, and those with experience with time series can
skip Chapter 8. These two mathematics chapters were separated to offer the reader a
chance to apply the first dose of statistics to boundary layer applications before delving
back into more math. Some students might have had a course on geophysical turbulence
or statistical fluid mechanics, and .can skim through the fll'st 5 chapters to get to the
boundary-layer applications. By excluding a few chapters, instructors can easily fit the
remaining material into a one-semester course. With supplemental readings, the book can
serve as a two-semester sequence in atmospheric turbulence and boundary-layer
meteorology.
Notational diversity proved to be the greatest difficulty. Each subdiscipline appeared
to have its own set of notation, which often conflicted with the notation of other
subdisciplines. To use the published notation would have lead to confusion. I was
therefore forced to select a consistent set of notation to use throughout the book. In most
cases I've tried to retain the notation most frequently used in the literature, or to add
subscripts or make logical extensions to existing notation. In other cases, I had to depart
from previously published notation. Readers are referred to Appendix B for a
comprehensive list of notation.
xi
ú = _ l r k a ^ o v i ^ v b o ú l o l i l d v =

I certainly cannot claim to be an expert in all the myriad subdisciplines of boundary-

layer meteorology. yet I knew the book should be comprehensive to be useful. My
interest and enthusiasm in writing this book motivated many trips to the library. and
stimulated my analysis of many research papers to learn the underlying themes and
common tools used in the diverse areas of boundary layer meteorology. Unfortunately.
the limited space within this book necessitated some difficult decisions regarding the
amount and level of material to include. Hopefully I've presented sufficient background
to lay the building blocks upon which more advanced concepts can be built by other
instructors and researchers. Certain topics such as atmospheric dispersion and agricultural
micrometeorology are not covered here. because there are other excellent books on these
subjects.
Many colleagues and friends helped with this book and contributed significantly to its
final form and quality. Michelle Vandall deserves special recognition and thanks for
drawing most of the figures. designing the page and chapter headers. and for her overall
dedication to this project. Colleagues Steve Stage. Larry Mahrt. George Young. Jacq
Schols. Chandran Kaimal. Steve Silberberg. Beth Ebert. and Bruce Albrecht reviewed
various chapters and provided valuable suggestions. Information regarding field
experiments was provided by Anton Beljaars. Ad Driedonks. Jean-Claude Andre. Jean-
Paul Goutorbe. Anne Jochum. Steve Nelson. Bob Murphy. Ruwim Berkowicz. Peter
Hildebrand. Don Lenschow. and others. Eric Nelson proofread the manuscript. and
helped with the list of notation. Sam Zhang compiled the index. Some of the equations
were typeset by Camille Riner and Michelle Vlllldall (her name keeps reappearing). Four
years of students in my micrometeorology courses at Wisconsin graciously tolerated
various unfinished drafts of the book. and caught many mistakes. The patient editors of
Kluwer Academic Publishers (formerly D. Reidel) provided constant encouragement. and
are to be congratulated for their foresight and advice. The American Meteorological
Society is acknowledged for permission to reproduce figures 10.7. 12.10 and 12.19. To
all of these people and the many more to whom I apologize for not naming here. I thank
you.
To my wife. Linda. I give my gratitude for her devotion during this exciting and
exhausting episode of my life.

Roland B. Stull
Boundary Layer Research Team
Madison. Wisconsin
 Z 1000
(m)
1000

500

Range from
lidar (km)

(a)

z
(m)

500

250

o
Range from
lidar (km)

(b)

Frontispiece lidar images of the aerosol-laden boundary layer, obtained during

the FIFE field experiment in Kansas. (a) Convective mixed layer
observed at 1030 local time on 1 July 1987, when winds were
generally less than 2 mls. (b) Slightly-stable boundary layer with
shear-generated turbulence, observed at 530 local time on 7 July
1987. Winds ranged from 5 mls near the surface to 15 mls near the
top of the boundary layer. Photographs from the Univ. of Wisconsin
lidar are courtesy of E. Eloranta, Boundary Layer Research Team.
 Layer Characteristics
From our flfSt breath, we spend most of our lives near the earth's surface. We feel the
warmth of the daytime sun and the chill of the nighttime air. It is here where our crops
are grown, our dwellings are constructed, and much of our commerce takes place. We
grow familiar with our local breezes and microclimates, and we sense the contrasts when
we travel to other places.
Such near-earth characteristics, however, are not typical of what we observe in the rest
of the atmosphere. One reason for this difference is the dominating influence of the earth
on the lowest layers of air.
The earth's surface is a boundary on the domain of the atmosphere. Transport
processes at this boundary modify the lowest tOO to 3000 m of the atmosphere, creating
what is called the boundary layer (Fig 1.1). The remainder of the air in the
troposphere is loosely called the/ree atmosphere. The nature of the atmosphere as
perceived by most individuals is thus based on the rather peculiar characteristics found in a
relatively shallow portion of the air.

_ _ _ Tropopause --------r----_
Free Atmosphere

Fig. 1.1 The troposphere can be divided into two parts: aboundary
layer (shaded) near the surface and the free atmosphere
above it.
2 BOUNDARY LAYER METEOROLOOY

This chapter is meant to provide a descriptive overview of the boundary layer.

Emphasis is placed on mid-latitude boundary layers over land, because that is where much
of the world's population resides and where many boundary layer measurements have
been made. In this region the diurnal cycle and the passage of cyclones are the dominant
forcing mechanisms. Tropical and maritime boundary layers will also be briefly
reviewed. Theories and equations presented in later chapters will make more sense when
put into the context of the descriptive nature of the boundary layer as surveyed here.

1.1 A Boundary-Layer Definition

The troposphere extends from the ground up to an average altitude of 11 lan, but often
only the lowest couple kilometers are directly modified by the underlying surface. We can
define the boundary layer as that part of the troposphere that is directly
influenced by the presence of the earth's surface, and responds to surface
forcings with a timescale of about an hour or less. These forcings include
frictional drag, evaporation and transpiration, heat transfer, pollutant emission, and terrain
induced flow modification. The boundary layer thickness is quite variable in time and
space, ranging from hundreds of meters to a few kilometers.
An example of temperature variations in the lower troposphere is shown in Fig 1.2.
These time-histories were constructed from rawinsonde soundings made every several
hours near Lawton, Oklahoma. They show a diurnal variation of temperature near the
ground that is not evident at greater altitudes. Such diurnal variation is one of the key
characteristics of the boundary layer over land. The free atmosphere shows little diurnal
variation.

30
Fig. 1.2
Evolution of
temperatures U
.
0
measured near - 20
GI
the ground
(97.5 kPa) and
at a height of
roughly 1100 m
.
:J
i;
GI
Q. 10
85 (kPa)

E
above ground
(85 kPa) . ...
GI
Lawton,Oklahoma
1983
Based on Noon Noon Noon
rawinsonde 0
June 7 June 10
lauches from
[Link],OK. Time

This diurnal variation is not caused by direct forcing of solar radiation on the boundary
layer. Little solar radiation is absorbed in the boundary layer; most is transmitted to the
ground where typical absorptivities on the order of 90% result in absorption of much of
 MEAN BOUNDARY LAYER CHARACfERlSTICS 3

the solar energy. It is the ground that warms and cools in response to the radiation, which
in turn forces changes in the boundary layer via transport processes. Turbulence is one
of the important transport processes, and is sometimes also used to define the boundary
layer.
Indirectly, the whole troposphere can change in response to surface characteristics, but
this response is relatively slow outside of the boundary layer. Hence, our defmition of the
boundary layer includes a statement about one-hour time scales. This does not imply that
the boundary layer reaches an equilibrium in that time, just that alterations have at least
begun.
Two types of clouds are often included in boundary-layer studies. One is the fair-
weather cumulus cloud. It is so closely tied to thermals in the boundary layer that it
is difficult to study the dynamics of this cloud type without focusing on the triggering
boundary-layer mechanisms. The other type is the stratocumulus cloud. It fills the
upper portion of a well-mixed, humid boundary layer where cooler temperatures allow
condensation of water vapor. Fog, a stratocumulus cloud that touches the ground, is
also a boundary-layer phenomenon.
Thunderstorms, while not a surface forcing, can modify the boundary layer in a matter
of minutes by drawing up boundary-layer air into the cloud, or by laying down a carpet of
cold downdraft air. Although thunderstorms are rarely considered to be boundary layer
phenomena, their interaction with the boundary layer will be reviewed in this book.

1.2 Wind and Flow

Air flow, or wind. can be divided into three broad categories: mean wind,
turbulence, and waves (Fig 1.3). Each can exist separately, or in the presence of any
of the others. Each can exist in the boundary layer, where transport of quantities such as
moisture, heat, momentum. and pollutants is dominated in the horizontal by the mean
wind, and in the vertical by turbulence.

----
Fig. 1.3 (a>
Idealization of
(a) Mean wind
alone. (b) waves
alone. and (c)
h (b)
•t
turbulence alone .
In reality waves
or turbulence are
often super-
Kú= r
imposed on a (e)
mean wind. U is
the component
of wind in the
x·direction.
4
Wú g í ê j K = t
4 BOUNDARY LAYER METEOROLOGY

Mean wind is responsible for very rapid horizontal transport, or advection.

Horizontal winds on the order of 2 to 10 mls are common in the boundary layer. Friction
causes the mean wind speed to be slowest near the ground. Vertical mean winds are much
smaller, usually on the order of millimeters to centimeters per second.
Waves, which are frequently observed in the nighttime boundary layer, transport little
heat, humidity, and other scalars such as pollutants. They are, however, effective at
transporting momentum and energy. These waves can be generated locally by mean-wind
shears and by mean flow over obstacles. Waves can also propagate from some distant
source, such as a thunderstorm or an explosion.
The relatively high frequency of occurrence of turbulence near the ground is one of the
characteristics that makes the boundary layer different from the rest of the atmosphere.
Outside of the boundary layer, turbulence is primarily found in convective clouds, and
near the jet stream where strong wind shears can create clear air turbulence (CAT).
Sometimes atmospheric waves may enhance the wind shears in localized regions,
causing turbulence to form. Thus, wave phenomena can be associated with the turbulent
transport of heat and pollutants, although waves without turbulence would not be as
effective.
A common approach for studying either turbulence or waves is to split variables such
as temperature and wind into a mean part and a perturbation part. The mean part
represents the effects of the mean temperature and mean wind, while the penurbation part
can represent either the wave effect or the turbulence effect that is superimposed on the
mean wind.
As will be seen in later chapters of this book, such a splitting technique can be applied
to the equations of motion,. creating a number of new terms. Some of these terms,
consisting of products of perturbation variables, describe lIolllinear interactions
between variables and are associated with turbulence. These terms are usually neglected
when wave motions are of primary interest. Other tetms, containing only one penurbation
variable, describe linear motions that are associated with waves. These terms are
neglected when turbulence is emphasized.

1.3 Turbulent Transport

Turbulence, the gustiness superimposed on the mean wind, can be visualized as
consisting of irregular swirls of motion called eddies. Usually turbulence consists of
many different size eddies superimposed on each other. The relative strengths of these
different scale eddies define the turbulence spectrum.
Much of the boundary layer turbulence is generated by forcings from the ground. For
example, solar heatillg of the ground during sunny days causes thermals of warmer
air to rise. These themlals are just large eddies. Frictional drag on the air flowing over
the ground c:tuses willd shears to develop, which frequently become turbulent.
Obstacles like trees and buildings deflect the flow, causing turbulent wakes adjacent
to, and downwind of the obstacle.
The largest boundary layer eddies scale to (ie, have sizes roughly equal to) the depth
of the bound:try layer; that is, 100 to 3000 m in diameter. These are the most intense
 MEAN BOUNDARY LAYER CHARACfERISTICS 5

eddies because they are produced directly by the forcings discussed above. "Cats paws"
on lake surfaces and looping smoke plumes provide evidence of the larger eddies.
< Smaller size eddies are apparent in the swirls of leaves and in the wavy motions of the
grass. These eddies feed on the larger ones. The smallest eddies, on the order of a few
millimeters in size, are very weak because of the dissipating effects of molecular viscosity.
Turbulence is several orders of magnitude more effective at transporting quantities
than is molecular diffusivity. It is turbulence that allows the boundary layer to respond to
changing surface forcings. The frequent lack of turbulence above the boundary layer
means that the rest of the free atmosphere cannot respond to surface changes. Stated more
directly, the free atmosphere behaves as if there were no boundary to contend with, except
in sense of mean wind flowing over the boundary-layer-top height contours.

1.4 Taylor's Hypothesis

We often need information on the size of eddies and on the scales of motions in the
boundary layer. Unfortunately, it is difficult to create a snapshot picture of the boundary
layer. Instead of observing a large region of space at an instant in time, we find it easier to
make measurements at one point in space over a long time period. For example,
meteorological instruments mounted on a tower can give us a time record of the boundary
layer as it blows past our sensors.
In 1938, G. I. Taylor suggested that for some special cases, turbulence might be
considered to be frozen as it advects past a sensor. Thus, the wind speed could be used
to translate turbulence measurements as a function of time to their corresponding
measurements in space. We must keep in mind that turbulence is not really frozen.
Taylor's simplification is thus useful for only those cases where the turbulent eddies
evolve with a timescale longer than the time it takes the eddy to be advected past a sensor
(Powell and Elderkin, 1974).
Let U and Y represent the eastward-moving and northward-moving Cartesian wind
components, and let M represent the total wind magnitude (speed) given by M2 = U2 +
y2. If an eddy of diameter A. is advected at mean wind speed M, then the time period P for
it to pass by a stationary sensor is given by P = A. / M .
Suppose that some variable like temperature varies from one side of the eddy to the
other. Then the temperature measured at our sensors would vary with time as the eddy
advects past. For example, Fig l.4a shows an initial condition where an eddy of 1()() m in
diameter is beginning to advect past our sensor tower. The leading side of the eddy, at a
temperature of 10°C, is warmer than the trailing side, at only 5°C. This is our spatial
description of the eddy at an instant in time. Namely, the temperature gradient across the
eddy is i)T/i)xd =0.05 KIm, where xd is in a direction parallel to the mean wind.
At that initial instant, our sensor on the tower measures a temperature of 10°C. If the
wind speed were 10 mis, then the sensor would measure a temperature of 5°C ten seconds
later, assuming that the eddy has not changed as it advected by (Fig l.4b). The local
change of temperature with time measured at our sensor is i)T/i)t = -0.5 Kls. We see that:
6 BOUNDARY LAYER METEOROLOOY

(8) (b)
N Q J J N MMã ú =

--/

Fig. 1.4 illustration of Taylofs hypothesis. (a) An eddy that is 100 m in

diameter has a 5 0 C temperature difference across it. (b) The
same eddy 10 seconds later is blown downwind at a wind speed
of 10 mls.

dT
(1.4a)
at
which is an expression of Taylor's hypothesis for temperature in one dimension.
For any variable ú I = Taylor's hypothesis states that turbulence is frozen when d!;/dt =
O. But the total derivative is defined by: d!;/dt = d!;/dt + U d!;/dX + V d!;/dy + W!;/dZ .
Thus, the general form of Taylor's hypothesis is

-U Ç ú = _V Ç ú = _W Ç ú = ( lAb)
dX dy dZ

This hypothesis can also be stated in terms of a wavenumber, K, and frequency, f:

K = flM (lAc)

where K = 21t().., and f = 21t/'P, for wavelength A. and wave period lP (Wyngaard and
Clifford, 1977). The dimensions of K are radians per unit length, while f has dimensions
of radians per unit time.
To satisfy the requirements that the eddy have negligible change as it advects past a
sensor, Willis and Deardorff (1976) suggest that

(lAd)

where aM' the standard deviation (see chapter 2 for a review of statistics) of wind speed ,
 MEAN BOUNDARY LAYER CHARACfERISTICS 7

is a measure of the intensity of turbulence. Thus, Taylor's hypothesis should be

satisfactory when the turbulence intensity is small relative to the mean wind speed.

1.5 Virtual Potential Temperature

Buoyancy is one of the driving forces for turbulence in the BL. Thermals of warm
air rise because they are less dense than the surrounding air, and hence positively
buoyant. Virtual temperature is a popular variable for these studies because it is the
temperature that dry air must have to equal the density of moist air at the same pressure.
Thus, variations of virtual temperature can be studied in place of variations in density.
Water vapor is less dense than dry air; thus, moist unsaturated air is more buoyant
than dry air of the same temperature. The virtual temperature of unsaturated moist air is
therefore always greater than the absolute air temperature, T. Liquid water, however, is
more dense than dry air, making cloudy air heavier or less buoyant than the corresponding
cloud-free air. The suspension of cloud droplets in an air parcel is called liquid water
loading, and it always reduces the virtual temperature.
Virtual Dotential temperatures are analogous to potential temperatures in that
they remove the temperature variation caused by changes in pressure altitude of an air
parcel. Turbulence includes vertical movement of air, making a variable such as virtual
potential temperature not just attractive but almost necessary.

1.5.1 Definitions

For saturated (cloudy) air, the virtual potential temperature, av , is defined by:
(1.5 . la)

where rsat is the water-vapor saturation mixing ratio of the air parcel, and rL is the liquid-
water mixing ratio. In (1.S.1a) the potential temperatures are in units of K, and the
mixing ratios are in units of gig. For unsaturated air with mixing ratio r, the virtual
potential temperature is:

av = a·( 1 + 0.61·r) (1.S .1b)

A derivation of the virtual temperature is given in Appendix D.

As usual , the potential temperature, a, is defined as

(
Po Jo.286
a = T P (I.5.1c)

where P is air pressure and Po is a reference pressure. Usually, Po is set to 100 kPa
8 BOUNDARY LAYER METEOROLOGY

(1000 mb), but sometimes for boundary layer work the surface pressure is used instead.
To flfSt order, we can approximate the potential temperature by

([Link])

where z is the height above the 100 kPa (1000 mb) level, although sometimes height
above ground level (agl) is used instead. The quantity g/Cp = 0.0098 KIm is just the
negative of the dry adiabatic lapse rate (9.8 °C per kilometer), where g is the gravitational
acceleration and C p is the specific heat at constant pressure for air. Sometimes the
quantity Cp·S is called the dry static energy.

Fig. 1.5
Example of the
difference between
mean é ç í É å í á ú =
temperature. a .
and mean virtual
potential tempera-
ture. 9v • Qiven
observations of
mixing ratio. T.
e Elv
and ~ Ä ë ç ä ì í É ú = 90
temperature T .
Dew point.
shown.
t. is also 100 o!:-,-,...L..L+.l..4- OR? Dl ú OS MK K i J K ç K K K K K g ONi K K K l J DJ L F OU? ? l ú OV J DJ l ? ? ? ? P Mú l ú P K K i Dl ú P J DOM=
(g/kg) Temperature (K)

An example of the difference between potential temperature and virtual potential

temperature is shown in Fig 1.5 for a case of moist unsaturated air. The difference is
small, but not negligible. Only when the air is very dry can we neglect the difference.

1.5.2 Example

Problem. Given a temperature of 25°C and a mixing ratio of 20 g/kg measured at a

pressure of 90 kPa (900 mb), find the virtual potential temperature.

Solution. First, we must find the potential temperature:

S = T (PJP) 0.286 = 298.16.(100/90) 0.286 = 307.28 K

The air is unsaturated. allowing us to find the virtual potential temperature from:

Sv = S·( 1 + 0.61 ·r) = 307.28·[ 1 + 0.61·(0.020)] = 311.03 K

 MEAN BOUNDARY LAYER CHARACfERISTICS 9

Discussion. Even though the virtual potential temperature is only about 4 K warmer
than the potential temperature, this difference is on the same order as the difference
between the warm air rising in thermals and the surrounding environment. Thus, neglect
of the humidity in buoyancy calculations could lead to erroneous conclusions regarding
convection and turbulence.

1.6 Boundary Layer Depth and Structure

Over oceans, the boundary layer depth varies relatively slowly in space and time. The
sea surface temperature changes little over a diurnal cycle because of the tremendous
mixing within the top of the ocean. Also, water has a large heat capacity, meaning that it
can absorb large amounts of heat from the sun with relatively little temperature change.
Thus, a slowly varying sea surface temperature means a slowly varying forcing into the
bottom of the boundary layer.
Most changes in boundary layer depth over oceans are caused by synoptic and
mesoscale processes of vertical motion and advection of different air masses over the sea
surface. An air mass with a temperature different than that of the ocean will undergo a
modification as its temperature equilibrates with that of the sea surface. Once equilibrium
is reached, the resulting boundary layer depth might vary by only 10% over a horizontal
distance of 1000 km. Exceptions to this gentle variation can occur near the borders
between two ocean currents of different temperatures (Stage and Weller, 1976).

Subelden.,.

11

Fig. 1.6 Schematic of synoptic - scale variation of boundary layer

depth between centers of surface high (H) and low (L)
pressure. The dotted line shows the maximum height
reached by surface modified air during a one-hour period.
The solid line encloses the shaded region, which is most
studied by boudary-Iayer meteonogists.
10 BOUNDARY LAYER METEOROLOGY

Over both land and oceans, the general nature of the boundary layer is to be thinner in
high-pressure regions than in low-pressure regions (Fig 1.6). The subsidence and low-
level horizontal divergence associated with synoptic high pressure moves boundary layer
air out of the high towards lower pressure regions. The shallower depths are often
associated with cloud-free regions. If clouds are present, they are often fair-weather
cumulus or stratocumulus clouds.
In low pressure regions the upward motions carry boundary-layer air away from the
ground to large altitudes throughout the troposphere. It is difficult to define a boundary-
layer top for these situations. Cloud base is often used as an arbitrary cut-off for
boundary layer studies in these cases. Thus, the region studied by boundary layer
meteorologists may actually be thinner in low-pressure regions than in high pressure ones
(see Fig 1.6).
Over land surfaces in high pressure regions the boundary layer has a well defined
structure that evolves with the diurnal cycle (Fig 1.7). The three major components of this
structure are the mixed layer, the residual layer, and the stable boundary layer.
When clouds are present in the mixed layer, it is further subdivided into a cloud layer
'and a subcloud layer.
The sUrface layer is the region at the bottom of the boundary layer where turbulent
fluxes and stress vary by less than 10% of their magnitude. Thus, the bottom 10% of the
boundary layer is called the surface layer, regardless of whether it is part of a mixed layer
or stable boundary layer. Finally, a thin layer called a microlayer or interfacial layer
has been identified in the lowest few centimeters of air, where molecular transport
dominates over turbulent transport.
The following shorthand notatiol" is often used for the various parts of the boundary
layer. For the sake of completeness, some additional terms are listed here that will not be
discussed until later:
BL Boundary layer (also known as the planetary boundary layer, PBL, or the
atmospheric boundary layer, ABL)
CL Cloud layer
FA Free atmosphere
mL Internal boundary layer
ML Mixed layer (also known as the convective boundary layer, CBL)
RL Residual layer
SBL Stable boundary layer (also known as the nocturnal boundary layer, NBL)
SCL Subcloud layer
SL Surface layer (the bottom 10% of the boundary layer)

The tops of four of these layers are given the following symbols:
h Top of the stable boundary layer (often defined as the top of the NBL)
Zj Top of the mixed layer (often defmed as the average base of the overlying stable
layer)
1y Top of the residual layer (often defined as the average base of the overlying
stable layer)
zb Top of the subcloud layer (this is the height of cloud base; usually near the
lifting condensation level, LCL)
 2000

Free Atmosphere

Capping Inversion

g
ú = Residual layar
CO ú=
ú =
::t:

i
ú=
-<
ú=
o
ú=
g
Local Time

Fig. 1.7 The boundary layer in high pressure regions over land consists of three major parts: a very
turbulent mixed layer; aless-turbulent residual layer containing former mixed-layer air; and a
nocturnal stable boundary layer of sporadic turbulence. The mixed layer can be subdivided into a en
cloud layer and a subcloud layer. Time markers indicated by Sl-SS will be used in Fig. 1.12.
I
::l
Q

--
12 BOUNDARY LAYER METEOROLOGY

1.6.1 Mixed Layer

The turbulence in the mixed layer is usually convectively driven, although a

nearly well-mixed layer can form in regions of strong winds. Convective sources include
heat transfer from a warm ground surface, and radiative cooling from the top of the cloud
layer. The fIrst situation creates thermals of warm air rising from the ground, while the
second creates thermals of cool air sinking from cloud top. Both can occur
simultaneously, particularly when a cool stratocumulus topped mixed layer is being
advected over warmer ground.
Even when convection is the dominant mechanism, there is usually wind shear
across the top of the ML that contributes to the turbulence generation. This free-shear
situation is more akin to CAT, and is thought to be associated with the formation and
breakdown of waves in the air known as Kelvin-Helmholtz waves.
On initially cloud-free days, however, ML growth is tied to solar heating of the
ground. Starting about a half hour after sunrise, a turbulent ML begins to grow in depth.
This ML is characterized by intense mixing in a statically unstable situation where thermals
of warm air rise from the ground (Fig 1.8). The ML reaches its maximum depth in late
afternoon. It grows by entraining, or mixing down into it, the less turbulent air from
above.

Fig. 1.8 Idealization of thermals in a mixed layer. Smoke plumes

loop up and down in the mixed layer eventually becoming
uniformly distributed.

The resulting turbulence tends to mix heat, moisture, and momentum uniformly in the
vertical. Pollutants emitted from smoke stacks exhibit a characteristic looping as those
portions of the effluent emitted into warm thermals begin to rise (Fig l.8). The resulting
profIles of virtual potential temperature, mixing ratio, pollutant concentration, and wind
speed frequently are as sketched in Figure 1.9.
Virtual potential temperature profiles are nearly adiabatic in the middle portion of the
ML. In the surface layer one often finds a superadiabatic layer adjacent to the ground.
 MEAN BOUNDARY LAYER CHARACTERISTICS 13

Free
Atmosphere
...... i" Entrainment zone
...... (.

\ú=
úJ = \\-8ú =
0;

······l
yú=

M T
L -_ _ ú F =WWWúWW?=
Fig. 1.9 Typical daytime ruofiles of l]'1ean virtual potential temperature
9v • wind speed M (where M'=11+7). water vapor mixing
ratio T, and pollutant concentration c.

A stable layer at the top of the ML acts as a lid to the rising thermals, thus restraining the
domain of turbulence. It is called the entrainment zone because entrainment into the
ML occurs there.
At times this capping stable layer is strong enough to be classified as a temperature
inversion; that is, the absolute temperature increases with height. In fact, it is frequently
called an inversion layer regardless of the magnitude of the stability. The most
common symbol for ML depth is zi' which represents the average height of the inversion
base.
Wind speeds are subgeostrophic throughout the ML, with wind directions
crossing the isobars at a small angle towards low pressure. The middle portion of
the ML frequently has nearly constant wind speed and direction. Wind speeds decrease
towards zero near the ground. resulting in a wind speed profile that is nearly
logarithmic with height in the surface layer. Wind directions cross the isobars at
increasingly large angles as the ground is approached, with 45 degree angles not
uncommon near the surface.
Mixing ratios tend to decrease with height, even within the center portion of the ML.
This reflects the evaporation of soil and plant moisture from below, and the entrainment of
drier air from above. The moisture decrease across the top of the ML is very pronounced,
and is often used together with potential temperature profiles to identify the ML top from
rawinsonde soundings.
Most pollutant sources are near the earth's surface. Thus, pollutant concentrations can
build up in the ML while FA concentrations remain relatively low . Pollutants are
transported by eddies such as thermals; therefore, the inability of thermals to penetrate
very far into the stable layer means that the stable layer acts as a lid to the pollutants too.
Trapping of pollutants below such an "inversion layer" is common in high- pressure
regions. and sometimes leads to pollution alerts in large communities.
14 BOUNDARY LAYER METEOROLOOY

As the tops of the highest thennals reach greater and greater depths during the course
of the day, the highest thennals might reach their lifting condensation level, LCL, if
sufficient moisture is present. The resulting fair-weather clouds are often targets for
soaring birds and glider pilots, who seek the updraft of the thennals.
High or middle overcast can reduce the insolation at ground level. This, in turn,
reduces the intensity of thennals. On these days the ML may exhibit slower growth, and
may even become non turbulent or neutrally-stratified if the clouds are thick enough.

1.6.2 Residual Layer

About a half hour before sunset the thennals cease to fonn (in the absence of cold air
advection), allowing turbulence to decay in the fonnerly well-mixed layer. The resulting
layer of air is sometimes called the residual layer because its initial mean state variables
and concentration variables are the same as those of the recently-decayed mixed layer.
For example, in the absence of advection, passive tracers dispersed into the daytime
mixed layer will remain aloft in the RL during the night. The RL is neutrally stratified,
resulting in turbulence that is nearly of equal intensity in all directions. As a result, smoke
plumes emitted into the RL tend to disperse at equal rates in the vertical and lateral
directions, creating a cone-shaped plume. Figure 1.10 shows a sketch of coning.

Fig . 1.10 The static stability decreases with height in the nocturnal boundary
layer. gradually blending into the neutrally-stratified residual layer
aloft, as indicated by the isentropic surfaces sketched on the left.
Smoke emissions into the stable air fan out in the horizontal with
little vertical dispersion other than wavelike OScillations. Smoke
emissions in the neutral residual-layer air spread with an almost
equal rate in the vertical and horizontal, allowing the smoke plume
to assume a cone-like shape.
 MEAN BOUNDARY LAYER CHARACfERISTICS 15

Nonpassive pollutants may react with other constituents during the night to create
compounds that were not originally emitted from the ground. Sometimes gaseous
chemicals may react to form aerosols or particulates which can precipitate out. The RL
often exists for a while in the mornings before being entrained into the new ML. During
this time solar radiation may trigger photochemical reactions among the constituents in the
RL.
Moisture often behaves as a passive tracer. Each day, more moisture may be
evaporated into the ML and will be retained in the RL. During succeeding days, the re-
entrainment of the moist air into the ML might allow cloud formation to occur where it
otherwise might not.
Variables such as virtual potential temperature usually decrease slowly during the night
because of radiation divergence. This cooling rate is on the order of 1 °C/d. The cooling
rate is more-or-less uniform throughout the depth of the RL, thus allowing the RL virtual
potential temperature profile to remain nearly adiabatic. When the top of the next day's
ML reaches the base of the RL, the ML growth becomes very rapid.
The RL does not have direct contact with the ground. During the night, the nocturnal
stable layer gradually increases in thickness by modifying the bottom of the RL. Thus,
the remainder of the RL is not affected by turbulent transport of surface-related properties
and hence does not really fall within our definition of a boundary layer. Nevertheless, we
will include the RL in our studies as an exception to the rule.

1.6.3 Stable Boundary Layer

As the night progresses, the bottom portion of the residual layer is transformed by its
contact with the ground into a stable boundary layer. This is characterized by statically
stable air with weaker, sporadic turbulence. Although the wind at ground level frequently
becomes lighter or calm at night, the winds aloft may accelerate to supergeostrophic
speeds in a phenomenon that is called the low-level jet or nocturnal jet.
The statically stable air tends to suppress turbulence, while the developing nocturnal
jet enhances wind shears that tend to generate turbulence. As a result, turbulence
sometimes occurs in relatively short bursts that can cause mixing throughout the SBL.
During the non turbulent periods, the flow becomes essentially decoupled from the
surface.
As opposed to the daytime ML which has a clearly defined top, the SBL has a poorly-
defined top that smoothly blends into the RL above (Fig 1.10 and 1.11). The top of the
ML is defined as the base of the stable layer, while the SBL top is defined as the top of the
stable layer or the height where turbulence intensity is a small fraction of its surface value.
Pollutants emitted into the stable layer disperse relatively little in the vertical. They
disperse more rapidly, or "fan out", in the horizontal. This behavior is called/anning,
and is sketched as the bottom smoke plume in Fig 1.10. Sometimes at night when winds
are lighter, the effluent meanders left and right as it drifts downwind.
Winds exhibit a very complex behavior at night. Just above ground level the wind
speed often becomes light or even calm. At altitudes on the order of 200 m above ground,
the wind may reach 10-30 mls in the nocturnal jet. Another few hundred meters above
16 BOUNDARY LAYER METEOROLOGY

z z
Free atmosphere
---------- >._------ Former entrainment zone
----------- ú KJ J J J J J J
Residual Layer

Stable Boundary layer

=-----.... 9 y '"""'---'---.M
Fig. 1.11 Mean virtual potential temperature, 6v , and wind speed,
M ,profiles for an idealized stable boundary layer in a
high-pressure region .

that, the wind speed is smaller and closer to its geostrophic value. The strong shears
below the jet are accompanied by a rapid change in wind direction, where the lower level
winds are directed across the isobars towards low pressure.
Touching the ground, however, is a thin (order of a few meters) layer of katabatic
or drainage winds. These winds are caused by the colder air, adjacent to the ground,
flowing downhill under the influence of gravity. Wind speeds of 1 mls at a height of I m
are possible. This cold air collects in the valleys and depressions and stagnates there.
Unfortunately, many weather stations are located in or near Valleys, where the observed
surface Winds bear little relationship to the synoptic-scale forcings at night.
Wave motions are a frequent occurrence in the SBL. The strongly stable NBL not
only suppons gravity waves, but it can trap many of the higher-frequency waves near the
ground. Vertical wave displacements of 100 m have been observed, although the
associated wind and temperature oscillations are relatively small and difficult to observe
without sensitive instruments.
SBLs can also form during the day, as long as the underlying surface is colder than
the air. These situations often occur during warm-air advection over a colder surface,
such as after a warm frontal passage or near shorelines.

1.6.4 Virtual Potential Temperature Evolution

Given the virtual potential temperature profiles from the previous subsections, it is
useful to integrate these profiles into our concept of how the boundary layer evolves. If
rawinsonde soundings were made at the times indicated by flags S I through S6 in Fig
1.7, then Fig 1.12 shows the resulting virtual potential temperature proflle evolution.
We see from these soundings that knowledge of the vinual potential temperature
prome is usually sufficient to identify the parts of the boundary layer. The structure of the
BL is clearly evident.
 MEAN BOUNDARY LA YER CHARACTERISTICS 17

II ú = II

FA FA

Fig. 1.12
Pro files of mean RL
ML
virtual potentigl
tempe rature. ev,
showing the SBL
boundary-layer
evo lution during e. e.
a d ui rnal cycle
startin g at abo ut ú = ú = M
1600 local tim e.
$1-$6 iden tify FA
each so und ing FA
with an associated
CL
launch time
indicated in RL
Fig. U. ML
SC L
ML

e. e.

Stated another way, knowledge of the virtual potential temperature lapse rate is usually
sufficient for determining the static stability. An exception to this rule is evident by
comparing the lapse rate in the middle of the RL with that in the middle of the ML. Both
are adiabatic; yet, the ML corresponds to statically unstable air while the RL contains
statically neutral air.
One way around this apparent paradox for the classification of adiabatic layers is to
note the lapse rate of the air immediately below the adiabatic layer. If the lower air is
superadiabaric, then both that superadiabatic layer and the overlying adiabatic layer are
statically unstable. Otherwise, the adiabatic layer is statically neutral. A more precise
definition of static stability is presented later.
It is obvious that as the virtual potential temperature profile evolves with time, so must
the behavior of smoke plumes. For example, smoke emitted into the top of the NBL or
into the RL rarely is dispersed down to the ground during the night because of the limited
turbulence . These smoke plumes can be advected hundreds of kilometers downwind from
their sources during the night.
Smoke plumes in the RL may disperse to the point where the bottom of the plume hits
the top of the NBL. The strong static stability and frequent reduction in turbulence
reduces the downward mixing into the NBL. The top of the smoke plume sometimes can
con tinue to rise into the neutral air. This is called loftillg (see Fig 1.13).
18 BOUNDARYLAYERME1EOROLOOY

Sunset
Fig. 1.13 Lofting of a smoke plume occurs when the top of the plume
grows upward into a neutral layer of air while the bottom is
stopped by a stable layer.

After sunrise a new ML begins to grow, eventually reaching the height of the elevated
smoke plume from the previous night. At this time, the elevated pollutants are mixed
down to the ground by ML entrainment and turbulence in a process that is called
fumigation . A sketch of this process is shown in Fig 1.14. An analogous process is

Sunset Sunrise F1 F2
Time
Fig. 1.14 Sketch ofthe fumigation process, where a growing mixed
tayer mixes elevated smoke plumes down to the ground.
Smoke plume 1 is fumigated at time F1 . while plume 2 is
fumigated at time F2.
 MEAN BOUNDARY LAYER CHARACfERISTICS 19

often observed near shorelines, where elevated smoke plumes in stable or neutral air
upstream of the shoreline are continuously fumigated downstream of the shoreline after
advecting over a warmer bottom boundary that supports ML growth.

1.7 Micrometeorology
Compared to the other scales of meteorological motions, turbulence is on the small
end. Figure 1.15 shows a classification scheme for meteorological phenomena as a
function of their time and space scales. Phenomena such as turbulence with space scales
smaller than about 3 km and with time scales shorter than about I h are classified as
microscale. Micrometeorology is the study of such small-scale phenomena.
It is evident that the study of the boundary layer involves the study of microscale
processes. For this reason, boundary layer meteorology and micrometeorology are
virtually synonymous. Since many of the early micrometeorological measurements were
made with sensors on short stands and towers, micrometeorology was often associated
with surface-layer phenomena.
Regardless of what you call it, the small-scale phenomena being studied here are so
transient in nature that the detenninistic description and forecasting of each individual eddy
is virtually impossible. As a result, micrometeorologists have developed three primary
avenues for exploring their subject:
• stochastic methods
• similarity theory
• phenomenological classifications.
Stochastic methods deal with the average statistical effects of the eddies.
Similarity theory involves the apparent common-behavior exhibited by many
empirically-observed phenomena, when properly scaled. In the phenomenological
methods, the largest size structures such as thermals are classified and sometimes
approached in a partially detenninistic manner.
Micrometeorology has always relied heavily on field experiments to learn more
about the boundary layer. Unfortunately, the large variety of scales involved and the
tremendous variability in the vertical require a large array of sensors including airborne
platforms and remote sensors. The relatively large costs have limited the scope of many
field experiments. Only a few general-purpose, large-scale boundary layer experiments
have been conducted.
Alternative studies have used numerical and laboratory simulations. Much of
the turbulence work has been performed in laboratory tanks, usually using liquids such as
water as the working medium. Although there have been many successful laboratory
studies of small-scale turbulence, there have been only a few simulations of larger
phenomena such as thermals. Wind tunnel studies have been used to observe the flow of
neutral boundary layers over complex terrain and buildings, although the difficulty of
stratifying the air has meant that typical daytime and nighttime boundary layers could not
be adequately simulated.
 N
0

ttl
TIME SCALE g
1 month 1 day 1 hour 1 minute 1 second s>-
IHurricanes Meso a ;:0
Fronts Scale -<
200 km Low-level jet
Cloud cluster & MCC!
T ú=
j Éë ç ú = , .!!! -<
ú l Ö ê ~ m Ü á ` =disturbances Scale Iú= m
ualilines I ;:0
1-20km <II
! 0
Vhunderstorms Mesoy <II
W I
...J Urban eHects Scale
« f å í É ê å ~ ä V ú ó =waves
u 1---2 km
r/) Boundary layers ;:0
...J Cumulus Å ä ç ì Ç ú =
Mieron
Scale
ú=
0
« Short gravity waves , r
ú = 1-200 m
«
D-
Dust devilsl Micro ú=
'8
r/)

...J
«
1--20 m
I Thermals
Wakes
Scale

Microy
n .!!!
II II -<
I- ú=
z Surface-layer Scale
<II
oN f--2m Plumes ú=
0: ú=
o
J: -02m Mechanical
turbulence Micro S
Scale
1--2 em Isotropic
turbulence
ú Oã ã = 1
Fig. 1.15 Typical time and space orders-of-magnitude for micro and mesoscales. (After Orlanski, 1975.)
 MEAN BOUNDARY LAYER CHARACfERISTICS 21

Numerical simulation using digital computers has been very popular since the 1960's,
with many discoveries that were verified years later with field studies. Most of these
computer simulations employ the stochastic method of modeling fluid flow.
Unfortunately, a difficulty known as the closure problem has meant that each of these
models has at best been able to only approximate the governing equations, with
uncertainties introduced via a necessary parameterization of the unknowns.
That same closure problem has limited the avenues for theoretical studies
involving analytical solutions. For the most part, only highly simplified approximations
to the boundary Jayer have been amenable to direct solution.
If anything, these difficulties have stimulated, rather than stifled, the work of
micrometeorologists. There is an underlying assumption in meteorology that
subsynoptic-scale phenomena such as turbulence might be responsible, in part, for the
difficulty in making quality weather forecasts beyond a few days. Thus, part of the effort
in boundary layer meteorology involves the search for adequate turbulence
parameterization schemes for larger-scale numerical forecast models.
Additional motivation has come from concern over our environment. Every species of
animal and plant modifies its environment; the human species, however, is in a position to
recognize the consequences of its pollution and take appropriate action. Since most of the
anthropogenic effluents are emitted from near-surface sources, the resulting dispersion of
the pollutants is tied to boundary layer processes. As a result, air-pollution
meteorology is an applied form of micrometeorology.
Other applications include agricultural meteorology, where airborne transport of
chemicals necessary to plant life is governed by turbulence. Nocturnal processes such as
frost formation warrant improved study and forecast methods for crop protection. Fog
and low stratocumulus, which inhibit aviation operations, are essentially boundary layer
phenomena. Wind-generated power, a popular energy source for centuries, has had a
recent increase in interest as wind turbines have been designed to extract energy more
efficiently from the boundary layer wind. Other structures such as bridges and buildings
must be designed to withstand wind gusts appropriate to their sites.

1.8 Significance of the Boundary Layer

The role of the boundary layer on our lives is put into perspective when we compare
the characteristics of the boundary layer and free atmosphere (Table I-I). A taste of the
importance of the BL is given in the following summary:
• People spend most of their lives in the BL.
• Daily weather forecasts of dew, frost, and maximum and minimum temperatures are
really BL forecasts.
• Pollution is trapped in the BL.
• Fog occurs within the BL.
• Some aviation, shipping, and other commerce activities conducted within it.
• Air masses are really boundary layers in different parts of the globe that have
equilibrated with their underlying surface. Baroclinicity is generated this way.
22 BOUNDARY LAYER METEOROLOGY

Table 1·1, Comparison of boundary layer and free atmosphere characteristics.

Property Boundary Layer Free Atmosphere

Turbulence • Almost continuously turbulent • Turbulence in convective clouds,

over its whole depth. and sporadic CAT in thin layers
of large horizontal extent.
Friction • Strong drag against the
earth's surtace. Large • Small viscous dissipation.
energy dissipation.

Dispersion • Rapid turbulent mixing in • Small molecular diffusion. Otten rapid

the vertical and horizontal. horizontal transport by mean wind.

Winds • Near logarithmic wind speed • Winds nearly geostrophic.

profile in the surtace layer.
Subgeostrophic, cross-
isobaric flow common.

Vertical Transport • Turbulence dominates. • Mean wind and cumulus-scale

dominate

Thickness • Varies between 100 m to 3 km Less variable. 8-18 km.

in time and space. Slow time variations.
Diurnal oscillations over land.

• The primary energy source for the whole atmosphere is solar radiation, which for the
most part is absorbed at the ground and transmitted to the rest of the atmosphere by
BL processes. About 90% of the net radiation absorbed by oceans causes
evaporation, amounting to the evaporation of about I m of water per year over all
the earth's ocean area. The latent heat stored in water vapor accounts for 80% of
the fuel that drives atmospheric motions.
• Crops are grown in the BL. Pollen distributed by boundary layer circulations.
• Cloud nuclei are stirred into the air from the surface by BL processes.
• Virtually all water vapor that reaches the FA is first transported through the BL by
turbulent and advective processes.
• Thunderstorm and hurricane evolution are tied to the inflow of moist BL air.
• Turbulent transport of momentum down through the BL to the surface is the most
important momentum sink for the atmosphere.
• About 50% of the atmosphere's kinetic energy is dissipated in the BL.
• Turbulence and gustiness affects architecture in the design of structures.
• Warm and cold fronts separate boundary layers of different temperature.
• Wind turbines extract energy from the BL winds.
 MEAN BOUNDARY LAYER CHARACfERlSTICS 23

• Wind stress on the sea surface is the primary energy source for ocean currents .
• Turbulent transport and advection in the BL move water and oxygen to and from
immobile life forms like plants.

Obviously the list could go on, but the main point is BL processes affect our lives
directly, and indirectly via its influence on the rest of the weather. In this book we
examine some of the processes that occur in the BL, develop some schemes for coping
with turbulence, and show how they can be applied to benefit mankind.
For additional infonnation on boundary layers, turbulence and micrometeorology, the
General References section below lists books and other secondary sources. A reference
section at the end of each chapter lists the specific articles sited in the chapter.
A summary table of frequently-used scaling variables and dimensionless groups is
given in Appendix A. Appendix B lists notation, including abbreviations, symbols and
acronyms . Appendix C gives values of frequently-used parameters and constants .
Additional appendices are specifically referenced in the text.

1.9 General References

Turbulence & Waves:

Batchelor, G.K., 1953: The theory of homogeneous turbulence. Cambridge Science
Classics Series. Cambridge Univ. Press. 197pp.
Frost, W. and T.H. Moulden, 1977: Handbook of Turbulence, Vol.I, Fundamentals and
Applications. Plenum Press, NY. 498 pp.
Gossard, E.E. and W.H. Hooke, 1975: Waves in the Atmosphere, Atmospheric
lnfrasound and Gravity Waves - their Generation and Propagation. Elsevier
Scientific Pub!. Co., NY. 456pp.
Hinze, J.O., 1975 : Turbulence (2nd ed). McGraw-Hill Series in Mechanical
Engineering. McGraw-Hill Book Co., NY 790pp.
Lighthill, J., 1978: Waves in Fluids. Cambridge University Press. England. 504pp.
Lumley, J.L. and H.A. Panofsky, 1964: The Structure of Atmospheric Turbulence.
Monographs and Texts in Physics and Astronomy. Vol XII. Interscience Pub\., John
Wiley & Sons, NY 239pp.
Monin, A.S. and A.M. Yaglom, 1973: Statistical Fluid Mechanics, Vols 1 & 2. Edited
by John Lumley. The MIT Press, Cambridge, MA . 769pp.
Panofsky, H.A., and J.A. Dutton, 1984: Atmospheric Turbulence, Models and Methods
for Engineering Applications. Wiley-Interscience, John Wiley & Sons, NY. 397pp.
Scorer, R.S., 1978: Environmental Aerodynamics. Ellis Horwood, Halsted Press, John
Wiley & Sons. London. 488pp.
Stanisic, M.M., 1985: The Mathematical Theory of Turbulence. Universitext, Springer-
Verlag, NY. 429pp.
Tennekes, H. and J.L. Lumley, 1982: A First Course in Turbulence (2nd Ed) . The MIT
Press, Cambridge, MA 300pp.
Townsend, A.A., 1976: The Structure of Turbulent Shear Flow (2nd Ed). Cambridge
24 BOUNDARY LAYER METEOROLOGY

University Press, Cambridge, England. 429pp.

Tritton, D.J., 1977: Physical Fluid Dynamics, Van Nostrand Reinhold. NY. 362pp.
Turner, J.S., 1973: Buoyancy Effects in Fluids, Cambridge Univ. Press. 367pp.
Van Dyke, M., 1982: An Album of Fluid Motion. The Parabolic Press, Stanford. 176 pp.

Boundary Layers:
Bhumralkar, e.M., 1975: A Survey of Parameterization Techniques for the Planetary
Boundary Layer in Atmospheric Circulation Models. Report R-1653-ARPA July
1975. ARPA order no. 189-1. 6PlO Information Processing Techniques Office,
Rand Corp, Santa Monica, CA 90406. 84pp.
Boundary Layer Meteorology, a journal published by Reidel, is devoted to boundary
layer topics.
Coantic, M.F., 1978: An Introduction to Turbulence in Geophysics, and Air-Sea
Interactions. NATP-AGARD. Available from NTIS, Springfield, VA 22161.
242pp.
Kraus, E.B., 1972: Atmosphere-Ocean Interaction. Clarendon Press, Oxford, England.
271pp.
Nieuwstadt, F.T.M. and H. van Dop, 1982: Atmospheric Turbulence and Air Pollution
Modeling. Kluwer Academic Publishers, Dordrecht, The Netherlands. 358pp.
Plate, EJ., 1971: Aerodynamic Characteristics of Atmospheric Boundary Layers. AEC
Critical Review Series, US Atomic Energy Commission, Office of Information
Services. Available as TID-25465 from NTIS, Springfield, V A 2215l. 190pp.
Roll, H. U., 1965: Physics of the Marine Atmosphere, Academic Press, NY. 426pp.
Schlichting, H., 1968: Boundary Layer Theory, 6th Ed. McGraw-Hili Series in
Mechanical Engineering, McGraw-HiD Book Co., NY 747pp.
Sorbjan, Z., 1988: Structure of the Atmospheric Boundary Layer. Prentice-Hall, NY.
300pp.
Stun, R.B., 1986: Boundary Layer Basics, A Survey of Boundary Layer Meteorology.
(available from the author), 5lpp. Also published as "Atmospheric Boundary Layer",
The Encyclopedia of Physical Science and Technology (edited by R.A. Meyers).
Academic Press, Inc., NY.
Wyngaard, J.e., 1980: Workshop on the Planetary Boundary Layer. Am. Meteor. Soc.,
45 Beacon St., Boston, MA 02108. 322pp.

Boundary Layer Experimental Techniques:

Lenschow, D.H. (Ed.), 1986: Probing the Atmospheric Boundary Layer. Am. Meteor.
Soc., 45 Beacon St., Boston, MA 02108. 269pp.
Vinnichenko, N.K., N.Z. Pinus, S.M. Shmeter, and G.N. Shur, 1980: Turbulence in
the Free Atmosphere (2nd Ed.). Consultants Bureau, Plenum Pub!., NY. 310 pp.

Micrometeorolof;Y:
Arya, S.P., 1988: Introduction to Micrometeorology. International Geophysical Series,
Vol 42. Academic Press, Inc., NY. 307pp.
Brutsaert, W., 1982: Evaporation into the Atmosphere, Theory, History, and
Applications. Kluwer Academic Publishers, Dordrecht, The Netherlands. 258pp.
 MEAN BOUNDARY LAYER CHARACfERISTICS 25

Geiger, R., 1965: The Climate Near the Ground. Harvard Univ . Press, Cambridge,
MA. 611pp.
Haugen, D.A ., (Ed .), 1973 : Workshop on Micrometeorology. Amer. Meteor. Soc., 45
Beacon St., Boston, MA 02108. 392pp.
Oke, T.R., 1988. Boundary Layer Climates, 2nd Ed .. Halsted Press, NY 435pp.
Sutton , O.G., 1953: Micrometeoro[ogy. McGraw-Hill (1977 reprint by Kreiger Pub\.
Co., Inc ., 645 NY Ave., Huntington, NY 11743). 333pp.

Chaos:
Bai-Lin, H., 1984: Chaos. World Scientific, Singapore. 386pp.
Cvitanovic, P., 1984: Universality in Chaos. Adam Hilger Ltd., Bristol, England.
513pp.
Gleick, J., 1987: Chaos, Making a New Science. Viking Penguin, Inc., NY . 352 pp.
Mandelbrot, B., 1977: The Fractal Geometry of Nature. Freeman, NY. 285pp.
Peitgen, H.-O. and P.H. Richter, 1986: The Beauty of Fractals. Springer-Verlag,
Berlin. 210pp.
Prigogine, I., and I. Stengers, 1984: Order Out of Chaos, Man's New Dialogue with
Nature. Bantam, NY. 349pp.

1.10 References for this Chapter

Orlanski, I., 1975: A rational subdivision of scales for atmospheric processes. Bull.
Amer. Meteor. Soc ., 56, 529-530.
Powell, D.C. and C.E. Elderkin, 1974: An investigation of the application of Taylor's
hypothesis to atmospheric boundary layer turbulence. 1 . Atmos. Sci ., 31, 990-
1002.
Stage, S.A . and R.A. Weller, 1986 : The frontal air-sea interaction experiment
(FASINEX); part II: experimental plan. Bull. Amer. Meteor. Soc., 67, 16-20.
Taylor, GJ., 1938: The spectrum of turbulence. Proc. R. Soc., A164, 476-490.
Willis, G.E. and l .W. Deardorff, 1976: On the use of Taylor's translation hypothesis for
diffusion in the mixed layer. Quart. 1. Roy. Meteor. Soc., 102, 817-822.
Wyngaard, J.e. and S.F. Clifford , 1977 : Taylor's hypothesis and high-frequency
turbulence spectra. 1. Atmos. Sci ., 34, 922-929.
26 BOUNDARY LAYER MElEOROLOGY

1.11 Exercises
1) a) Resketch figure 1.7 without any shading.
b) Next, shade on your figure those regions where the lapse rate is superadiabatic. (A
color like red is appropriate.)
c) In a different color (like dark blue), shade those regions that have a strongly stable
lapse rate.
d) With a different color (like light blue), shade those regions that have a weakly stable
lapse rate.
e) Leave the nearly adiabatic lapse rate regions unshaded. (Hint, look at the vinual
potential temperature graphs shown after figure 1.7.)
2) From your answer above, we see that a large portion of the mixed layer and the residual
layer have nearly adiabatic lapse rates. Nevenheless, there is a difference between the
two. Comment on this difference. How would smoke disperse in these regions?
3) If a very tall smoke stack could inject a smoke plume into the free atmosphere, then
what characteristic dispersion pattern would you expect? Options include looping,
fanning, coning, lofting, fumigation, and trapping within a layer.
4) Given the boundary layer evolution shown in Fig 1.7 and 1.12, sketch the modified
evolution that would occur if the sun was turned off at time SS. Also comment on the
processes that could effectively tum off the effects of the sun on the boundary layer.
S) If you observed the following measurement of wind speed on a strip chan, then what
diameter (in meters) are the major eddies? (Hint: First find the typical time duration of
an eddy, and then employ Taylor's Hypothesis.)

Wind 10
speed 5
( mls)
o 123456
Time (min.)

6) Suppose that the boundary layer can be approximated by a I km constant thickness

layer over the whole continental United States. Above that, assume that a constant
thickness (10 km) free atmosphere exists. Assume that each of these layers is well
mixed for simplicity. (You may need to consult synoptics textbooks or study a
number of weather maps to answer the following questions.)
a) On the average in spring, what fraction of the area of the continental United States is
covered by cyclones (low pressure regions)?
b) What is the average venical velocity magnitude out of the top of the boundary layer
in those cyclonic regions?
c) What is the average subsidence velocity acting on the top of the boundary layer in
anticyclonic (high pressure) regions? Remember that the boundary layer is of constant
 MEAN BOUNDARY LAYER CHARACfERISTICS 27

thickness. so the total inflow of air must equal the outflow. (Nonnally. subsidence
would not be able to cause air to penetrate down into the top of the boundary layer. but
for simplicity we will assume that the turbulent entrainment rate of air into the top of
the boundary layer is exactly equal to the subsidence velocity. Thus. all of the
subsiding air would enter the boundary layer.)
d) What are the average residence times. P • of air within the boundary layer and the
free atmosphere? Residence time is defmed as P=V/F. where F = volumetric flow rate
(volume/time) and V = volume. Neglect horizontal advection out of the United States.
e) Compare these residence times with the typical Rossby wave time period (that is.
the average time period between cyclone passage over a fixed point). Comment on the
significance of these times.
7) Would there be a boundary layer on a planet that had an atmosphere. but that did not
experience a diurnal variation of net radiation at the ground?
8) Define the boundary layer in your own tenns.
9) Add five items to the list in section 1.8 regarding the significance of the boundary layer.
Hint. look for things close to home. school. or work.
10) a) Given air at a pressure height of 90 kPa (900 mb) with a temperature of 30 °C and a
mixing ratio of 20 g/kg. fmd the virtual potential temperature.
b) Given saturated air at 85 kPa (850 mb) with a temperature of 20 °C and a total
water mixing ratio (i.e.• sum of vapor and liquid mixing ratios) of 20 g/kg find the
virtual potential temperature. (Hint. you might want to employ a thennodynamic
diagram.)
11) Given the following virtual potential temperature sounding. identify each layer
(example. ML. RL. SBL. FA). Estimate what time of day that sounding was made.
and then sketch the virtual potential temperature profile that you might expect four
hours later.