## 1. Introduction

Measuring and modeling of the boundary layer is challenging, especially for the nighttime stable boundary layer (SBL). In particular, high-quality in situ measurements of profiles of mean and turbulent statistics of the nighttime SBL are logistically difficult using instrumented towers or instrumented research aircraft (Tjernström 1993). A suite of fast response turbulence sensors attached to a tethered system can operate effectively under these nighttime conditions. The Tethered Lifting System (TLS) was developed for the 1999 Cooperative Atmosphere–Surface Exchange Study (CASES-99) field campaign (Poulos et al. 2002) by the Cooperative Institute for Research in Environmental Sciences (CIRES) at the University of Colorado (Balsley et al. 1998, 2003, 2006) and modified for measurements in the open spaces of the Washington, D.C., area (Frehlich et al. 2006) to provide high-resolution in situ profiles of the SBL.

High-resolution in situ measurements are able to resolve unusual features in the temperature, velocity, and turbulence structure of the boundary layer, the residual layer, and the free atmosphere. An ideal example is given in Fig. 8 of Frehlich et al. (2006), where the height of a well-defined mixing layer coincides with a small layer of cooler air, which could only be resolved with the TLS in situ sensors. The high-resolution velocity and temperature data provide unique measurements of the turbulent structure of stable regions, which are difficult to mimic in the laboratory, especially the low-turbulence regions above the SBL. This offers a unique opportunity to verify new turbulent scaling laws for stably stratified flows (Smyth and Moum 2000a, b; Werne and Fritts 1999, 2001; Riley and de Bruyn Kops 2003; Lindborg 2006). TLS measurements of a stably stratified boundary layer are presented and compared with the new theoretical predictions. Critical instrumentation issues are identified for accurate measurements of small-scale turbulence in these low-turbulence regimes.

The development of eye-safe Doppler lidar (Henderson et al. 1991, 1993; Grund et al. 2001; Pearson and Collier 1999) with range resolution on the order of 20–70 m and new processing algorithms (Banakh and Smalikho 1997; Frehlich 1997, 2001a; Frehlich and Cornman 2002; Frehlich et al. 1994, 1997, 1998, 2006; Banta et al. 2003, 2006; Davies et al. 2004; Lothon et al. 2006) have produced high-quality measurements of the boundary layer. For solid-state lidars, the lidar pulse is well approximated by a Gaussian shape and the critical parameters can be determined from the monitor pulse (Frehlich et al. 1994), and the effects of the spatial filtering of the radial velocity by the pulse can then be corrected (Banakh and Smalikho 1997; Frehlich 1997; Frehlich et al. 1998, 2006; Davies et al. 2004; Smalikho et al. 2005) to produce reliable estimates of turbulent statistics. This correction is possible when the lidar pulse parameters are stable and when the mean Doppler radial velocity estimates are given by a convolution of the instantaneous radial velocity with an effective pulse range–weighting function (Banakh and Smalikho 1997; Frehlich 1997). Accurate corrections for the pulse smoothing from other Doppler lidar systems, such as CO_{2} Doppler lidars, would be possible if these conditions are met (see Post and Cupp 1990; Pearson and Rye 1992). In some cases, such as measurements in the convective boundary layer (Lothon et al. 2006), the effects of the spatial smoothing can be ignored if the range gate is smaller than the length scale of the turbulence. Measurements of hazardous turbulence conditions around airports are also feasible with careful processing of the lidar data (Chan 2006). However, for low-turbulence conditions, especially the stable boundary layer, careful corrections for the pulse smoothing and the radial velocity estimation error is required for accurate turbulence estimates (Frehlich and Cornman 2002; Frehlich et al. 2006). In addition, the statistical accuracy of the estimates is improved by scanning over a three-dimensional volume with well-chosen scanning geometries. Recent interest in the effects of atmospheric turbulence on wind energy generation requires more reliable measurements of turbulence profiles in the SBL (Fritts 2005; Banta et al. 2006; Kelley et al. 2006). Recently, profiles of turbulence for the transverse velocity component have been developed based on the variations of the radial velocity as a function of the arc distance along an azimuth scan (Frehlich et al. 2006). A new lidar processing algorithm is presented for profiles of the turbulence statistics of the longitudinal velocity component (radial velocity component in the radial direction), which provides simultaneous profiles of the key turbulence statistics of two orthogonal horizontal components over a three-dimensional volume with vertical resolution sufficient for stable boundary layer studies.

High spatial and temporal resolution TLS measurements of the small-scale turbulence statistics have also determined the spatial intermittency in the SBL, which is the essence of Kolmogorov’s refined similarity theory (Kolmogorov 1962) and also provides the required input for rigorous error analysis (i.e., the statistical accuracy of turbulent estimates that are one source of self-correlation in many formulations of similarity theory such as Monin–Obuhkov similarity theory; Klipp and Mahrt 2004; Baas et al. 2006). In addition, the separation of the atmospheric variables into a turbulent component and a mean or slowly variable component can be ill posed, especially for statistics such as variances and fluxes for challenging conditions such as stable boundary layers (Kaimal and Finnigan 1994; Mahrt 1998; Vickers and Mahrt 2003) and larger-scale forcing such as with a density current and a frontal passage (Piper and Lundquist 2004). The standard analysis technique for these situations is selecting a suitable averaging interval to define the mean and turbulent quantities for both in situ observations (Kaimal and Finnigan 1994; Mahrt 1998; Vickers and Mahrt 2003) and remotely sensed measurements such as Doppler lidar (Banta et al. 2003, 2006). The optimal choice of the averaging time is difficult to quantify when there is no obvious scale separation between the forcing mechanisms and turbulence. The effects of the larger mesoscale forcing on convective boundary layer turbulence has been identified from the Boreas campaign where no well-defined length scale was observed in the spatial spectra and hence there is no clear separation between small-scale turbulence and large-scale forcing (Lenschow and Sun 2007).

Mesoscale models use parameterizations of the boundary layer that are challenging for the SBL (Roth 2000; Zilitinkevich and Baklanov 2002; Sorbjan 2006; Esau and Zilitinkevich 2006; Mauritsen and Svensson 2007; Zilitinkevich et al. 2007). Large-eddy simulations (LESs) are being used to investigate NWP model parameterization in neutral and stably stratified flows (Esau and Zilitinkevich 2006); however, the numerical requirements for stable boundary layers are challenging (Saiki et al. 2000; Sullivan et al. 2003). Better measurements of the SBL is required to advance the LES and NWP modeling effort, especially a consistent measurement and definition of boundary layer height (Roth 2000; Seibert et al. 2000; Zilitinkevich and Baklanov 2002; Balsley et al. 2006; Zilitinkevich et al. 2007).

Advanced data assimilation techniques are based on statistics of total observation error, which include the observation sampling error or “error of representativeness” resulting from the mismatch of the observation spatial average and the effective model spatial average (Frehlich 2001b, 2006; Frehlich and Sharman 2004). Improvements in next-generation models and data assimilation will require observations that more closely match the spatial average of the models and observations that provide estimates of the turbulence statistics that define the observation sampling error. New Doppler lidar scanning and processing algorithms were employed in Washington, D.C., to extract such spatially averaged mean and turbulence profiles of velocity (Frehlich et al. 2006). These results are extended to provide a more complete description of the boundary layer velocity statistics.

To address these diverse issues of atmospheric measurements, the paper is structured as follows: in section 2, we present an analysis of the instrumental requirements for accurate in situ measurements of stably stratified turbulence and comparisons of recent results to theoretical predictions; section 3 discusses a new Doppler lidar analysis technique for producing higher vertical resolution estimates of the radial velocity statistics; and section 4 focuses on the first simultaneous comparisons of Doppler lidar derived profiles of the three key turbulence statistics (i.e., energy dissipation rate *ε*, velocity variance *σ*^{2}, and outer scale of turbulence *L*_{0}) for the two horizontal velocity components (longitudinal and transverse) and the key instrumental and data processing issues for these measurements such as the impact of spatially filtering the radial velocity estimates to remove large-scale atmospheric processes. The paper concludes with a summary and discussion.

## 2. In situ TLS measurements

The TLS was developed for a variety of atmospheric measurements with a focus on measurements of the small-scale turbulence (Balsley et al. 1998, 2003, 2006; Muschinski et al. 2001). Temperature measurements are produced with a new low-frequency response solid-state temperature sensor and a high-frequency response fine cold-wire sensor. Low-frequency velocity measurements are made using a new two-axis sonic anemometer vaned into the wind and a high-frequency response fine hot-wire sensor. The temperature is calibrated to an absolute accuracy of better than 0.5 K and the accuracy of the linear calibration constant was better than 2%. The velocity is calibrated to an absolute accuracy of better than 0.5 m s^{−1} and the accuracy of the slope of the calibration curve is better than 5% (Frehlich et al. 2003).

*E*

_{1}(

*k*) at high wavenumber

*k*can be written as (Monin and Yaglom 1975, p. 354) where

*ν*is the kinematic viscosity,

*ε*is the energy dissipation rate, is the Kolmogorov microscale,

*C*

_{2}is the Kolmogorov constant, and the universal function is

*ϕ*

_{1}(

*x*) =

*C*

_{2}

*x*

^{−5/2}

*g*(

*βx*

^{2}) (see Frehlich et al. 2003; Azizyan et al. 1989), and where is the incomplete gamma function. The constant

*β*= 2.955 18 is based on the required normalization for the Kolmogorov constant

*C*

_{2}= 0.4977. The velocity spectrum

*E*

_{1}(

*k*) is only a function of the energy dissipation rate

*ε*and the kinematic viscosity of air

*ν.*This result relies on the universal scaling of high Reynolds number turbulence.

*ε*for homogeneous and isotropic turbulence is (Monin and Yaglom 1975, p. 56) which indicates that turbulent energy is dissipated by viscous forces at the high-wavenumber region (small spatial scales). We define the function

*P*(

*k*) as which defines the fraction of the energy dissipation rate produced from wavenumber less than

*k = y/η.*The functions

*g*(

*βη*

^{2}

*k*

^{2}) and

*P*(

*ηk*) are shown in Fig. 1. The median wavenumber for energy dissipation rate is

*k*≈ 0.2/

*η*and the 10% wavenumber

*k*

_{10}(10% on the energy dissipation rate is produced from wavenumber less than

*k*

_{10}) is

*k*

_{10}≈ 0.05/

*η*. Since

*η*≈ 0.005 m,

*k*

_{10}≈ 10 m

^{−1}, and the small scales define the energy dissipation rate and therefore high-resolution, in situ measurements are required to resolve these scales for a direct estimate of

*ε*.

*C*

^{2}

_{T}is the temperature structure constant; is the temperature inner scale;

*χ*is the thermal diffusivity of air; and describes the Hill bump (Hill 1978). This result also relies on the universal scaling of high Reynolds number turbulence.

*U*to scale temporal frequency

*f*to spatial frequency

*k*(i.e.,

*k*= 2

*πf/U*; Frehlich et al. 2003). Then, The best spectral estimate of

*ε*and

*C*

^{2}

_{T}are produced using the maximum likelihood algorithm with the temporal spectra

*S*(

_{u}*f*) and

*S*(

_{T}*f*), respectively, assuming Taylor’s frozen hypothesis is valid (Frehlich et al. 2003).

The TLS can be either operated at a constant altitude or in profiling mode by slowly ascending and descending at a rate of approximately 0.5–1.0 m s^{−1}, which provides similar resolution to the moving carriage of the Boulder Atmospheric Observatory (BAO) tower (Gossard and Frisch 1987). To better evaluate the high-wavenumber region of the temperature and velocity spectrum for low-turbulence conditions in stable regions, the hot-wire and cold-wire electronics and data acquisition system were modified and the sampling rate increased to 1 kHz to produce direct estimates of small-scale turbulence statistics such as *ε* in Eq. (5).

An example of along-stream velocity data *U*(*t*) from a low-turbulence region above the stable boundary layer (see Fig. 7, region A from 148- to 168-m altitude) is shown in Fig. 2 as well as the corresponding spectra and theoretical best-fit model with *ε* = 1.03 × 10^{−5} m^{2} s^{−3} using Eqs. (11) and (1) with only the good spectral data away from the obvious interference peaks. The error in the estimates of *ε* is approximately 10% based on the *χ*^{2} statistics of the useful spectral estimates in the frequency interval of 0.2–220 Hz. Note that a direct estimate of *ε* using Eq. (5) is not possible because of the interference, however, there is good agreement at the high frequencies above the interference indicating a direct measurement is possible with improved sensors. Also shown in Fig. 2 is the best-fit straight line and the ±10% deviation in the slope about the best fit to indicate a measure of the statistical accuracy of the gradient *S* = (∂*U*/∂*z*) = −0.0157 s^{−1}, where the altitude *z* is determined from a GPS sensor. The descend rate is 0.833 m s^{−1} and therefore the vertical and horizontal measurement region is 20 and 220 m, respectively. There is an obvious wave with a small amplitude of approximately 0.1 m s^{−1} and the standard deviation of the velocity about the best-fit line is *σ _{u}* = 0.051 m s

^{−1}and the length scale of the velocity fluctuations is defined (e.g., Monin and Yaglom 1975, p. 198) as

*L*

_{u}=

*σ*

^{3}

_{u}/

*ε*≈ 12.8 m.

*T*is converted to potential temperature by where

*T*is temperature in C,

*P*

_{0}is the reference pressure at the surface, and

*P*is the pressure at the measurement height. An example of potential temperature data from the same low-turbulence region as in Fig. 2 is shown in Fig. 3 as well as the corresponding spectra and theoretical best-fit model with

*C*

^{2}

_{T}= 8.13 × 10

^{−6}

*K*

^{2}m

^{−2/3}. The gradient of potential temperature is ∂Θ/∂

*z*= 0.0073 Km

^{−1}and is clearly bounded by the ±10% deviations around the best-fit line. Detector noise and interference from the nearby sonic anemometer (vibrational noise induced in the high gain electronics) dominates the high-frequency region of the spectrum. The best-fit model includes the 1/

*f*noise from the sensor, which appears at frequencies above 30 Hz. For this example, there is sufficient uncontaminated spectral power at lower frequencies to extract an accurate estimate of

*C*

^{2}

_{T}with an error of approximately 10% based on the

*χ*

^{2}statistics of the spectral estimates in the frequency interval of 0.2–30 Hz. Other physical parameters are

*ν*= 2.606 × 10

^{−5}m

^{2}s

^{−1},

*η*= 0.0064 m,

*U*= 9.15 m s

^{−1}, and the average temperature is

*= 26.4*T

^{∘}C.

The high-rate turbulence measurements are essential for investigating the turbulence parameters predicted by direct numerical simulations (DNS) of stably stratified turbulence (Werne and Fritts 1999, 2001; Gibson-Wilde et al. 2000; Riley and de Bruyn Kops 2003; Hebert and de Bruyn Kops 2006) with a constant potential temperature gradient and constant velocity gradient. The critical parameter regime is defined by the Brunt–Väisälä frequency *N* = (*g*/Θ)(∂Θ/∂*z*)^{1/2} = 0.0151 s^{−1}, the Ozmidov scale *L*_{OZ} = (*ε/N ^{3}*)

^{1/2}= 1.73 m, the Ozmidov frequency corresponding to

*L*

_{OZ}(i.e.,

*f*

_{OZ}=

*U/L*

_{OZ}= 5.26 Hz), the gradient Richardson number

*R*=

_{i}*N*

^{2}/

*S*

^{2}= 0.92, the buoyancy Reynolds number

*R*=

_{b}*ε*/(

*νN*) = 1747, the Reynolds number

^{2}*R*= 24963, and the local Froude number

_{e}= σ_{u}L_{u}/ν*F*=

_{h}*σ*/(

_{u}*NL*) = 0.26. Since the buoyancy Reynolds number

_{u}*R*is larger than 1000, the calculations of energy dissipation rate

_{b}*ε*from the velocity spectrum assuming isotropy of the small scales is valid (Smyth and Moum 2000b).

The recent theories of stably stratified turbulence in low Froude number (*F _{h}* < 1) predict a −5/3 spectral slope for frequencies less than the Ozmidov frequency (Riley and de Bruyn Kops 2003; Lindborg 2006). This is clearly shown in Figs. 2 and 3 for both the velocity and the temperature field. Since the sensor moves a vertical distance of approximately 0.9 m for a horizontal sampling distance of 10 m, there is little contribution from the fluctuations in the vertical direction since the vertical turbulence length scale is approximately

*L*= 3.4 m (Lindborg 2006). Also note that the Richardson number of approximately 1 supports the recent interpretation of the Richardson number in stratified flows (Majda and Shefter 1998). These turbulence structures are most likely the manifestation of the decay phase of some earlier larger-scale instability (e.g., a Kelvin–Helmholtz instability). Identifying the exact nature of the instability is not feasible with the limited data.

_{υ}= L_{u}F_{h}*ε*in terms of

*C*

^{2}

_{T}for well-defined turbulent layers in the steady state, especially in the radar literature (Gossard and Frisch 1987; Bertin et al. 1997): where

*γ*varies from approximately 0.8 to 2. The estimates of

*ε*based on the measurements in Figs. 3 and 4 for this range of

*γ*vary from 1.46 to 5.76 × 10

^{−7}m

^{2}s

^{−3}, which is inconsistent with the measured value of

*ε*= 1.26 × 10

^{−5}m

^{2}s

^{−3}. This turbulent regime is known to have large scatter in predictions of

*ε*from

*C*

^{2}

_{T}and local gradients (Bertin at al. 1997).

## 3. Doppler lidar measurements

The Doppler lidar used in this study is a 2-*μ*m eye-safe WindTracer lidar manufactured by Coherent Technologies, Inc. (CTI; Henderson et al. 1991, 1993). The lidar parameters are a Gaussian pulse width Δ*r* = 66.0 m (full width at half maximum), a range-gate length Δ*p* = 72.0 m, the effective range resolution Δ*R* = 85.8 m (Banakh and Smalikho 1997), an azimuth scan rate of 2.5° s^{−1}, a pulse repetition frequency (PRF) of 500 Hz, a measurement time per radial velocity profile *T* = 0.2 s using 100 lidar pulses per profile, and the azimuth spacing Δ*ϕ* = 0.5°. The scanning geometry is similar to the operations at Washington, D.C. (Frehlich et al. 2006) and a top view is shown in Fig. 4 along with the local buildings and the location of the TLS in the open space between the residential buildings in Lafayette, Colorado. The azimuth angle *ϕ* varied from 40° to 85° to focus on the open space. The elevation angle varied from −0.8° to 24° to provide useful profiles in both nighttime and daytime conditions (see also Frehlich et al. 2006). The lidar was approximately 20 m above the flat open space region. Range-gate distances from 420 to 2508 m were used for all the lidar analysis and the average transverse dimensions of the lidar sensing volume for each radial velocity estimate is Δ*h* ≈ 1464 sin 0.5° = 12.8 m, which is much less than the effective range resolution Δ*R* = 85.8 m and the assumption of zero Δ*h* is valid (Frehlich and Cornman 2002). The field campaign covered a 1-week period from 7 to 13 June 2006 and TLS data was only collected from 2300 to 0600 LT because of Federal Aviation Administration (FAA) restrictions. High daytime temperatures existed for most of this period and therefore all the nighttime data was statically stable because of radiative cooling.

Profiles of wind speed and direction are produced by a best fit to the radial velocity as a function of azimuth angle for a given fixed elevation angle and range-gate distance and combining all the estimates occupying a height interval of dimensions of 10 m (Browning and Wexler 1968; Frehlich et al. 2006). The turbulence component of the radial velocity is defined as the deviations of the radial velocity from the best-fit model for each range gate (see Fig. 3 of Frehlich et al. 2006). To investigate the effects of filtering the data to remove the large-scale processes, the turbulent velocity fluctuations are produced by subtracting a polynomial fit of order M as a function of azimuth for each range gate.

Unbiased estimates of the turbulence statistics of the velocity field can be produced from structure functions of the Doppler lidar radial velocities if 1) a universal description of the structure functions are known a priori (e.g., a von Kármán model), 2) the effects of the spatial averaging by the finite extent of the lidar pulse and processing algorithms are included in the analysis, and 3) the contribution from the estimation error in the radial velocities is removed, which is essential for the low-turbulence regimes (Frehlich and Cornman 2002; Frehlich et al. 2006). Two versions of the structure function method will be employed: the structure functions with the separation variable *s* along range (i.e., longitudinal velocity statistics) and with the separation variable *s* based on the azimuth angle separation (i.e., transverse velocity statistics). Both processing techniques correct for the effects of the spatial averaging by the lidar pulse using standard deconvolution algorithms (Banakh and Smalikho 1997; Frehlich et al. 1998, 2006; Frehlich and Cornman 2002; Smalikho et al. 2005).

*R*is the effective range resolution of the radial velocity estimates and

*θ*is the lidar elevation angle. Profiles of the transverse turbulent statistics are produced by processing all the data within a given altitude bin. An example of the average azimuth structure function, the corrected structure function, the best-fit von Kármán model, and the best-fit lidar structure function model is shown in Fig. 5 from the time period 0523–0540 UTC 7 June 2006 for raw data (no spatial filtering) and a fourth-order polynomial filter (

*M*= 4). Range gates with an altitude inside the 80–90-m height interval were processed for elevation angles up to 6°. Therefore, the largest vertical resolution is Δ

*z*

_{turb}= 15 m. The average correction to the structure functions from the radial velocity estimation error is 2

*σ*

^{2}

_{e}= 0.004 02 m

^{2}s

^{−2}and therefore

*σ*= 0.0449 m s

_{e}^{−1}is the average estimation error determined with the covariance technique (Frehlich 2001a; Frehlich et al. 2006) applied to the azimuth structure function for each range gate. The raw data (no filtering) indicates a large turbulent length scale of

*L*

_{0}

*= 748 m and the filtered data has a much lower length scale and lower standard deviation of*

_{υ}*σ*. However, the estimates of

_{υ}*ε*, the energy dissipation rate that defines the small-scale turbulence, are independent of the filtering. Similar results were produced for filter order from 1 to 4.

_{υ}Lidar-derived profiles of the turbulent statistics (i.e., *σ _{u}*,

*ε*, and

_{u}*L*

_{0}

*) for the radial velocity (i.e., longitudinal component) have been produced for a fixed lidar beam using the time-averaged radial velocity to define the mean and turbulent components (Frehlich et al. 1998). These profiles have a poor vertical resolution of ≈200 m, which is suitable for convective boundary layer studies but not for stable boundary layers research. The longitudinal and transverse turbulence parameters were estimated simultaneously using data from a zero elevation angle azimuth scan and using the definition of turbulence as the fluctuation about the best-fit model (Frehlich et al. 2006, their Fig. 4). Since the domain was statistically homogeneous, large lags in the radial direction were possible and isotropic statistics were observed. These results will now be extended to produce simultaneous profiles of the longitudinal and transverse turbulence parameters.*

_{u}*R*

_{min}and

*R*

_{max}are the minimum and maximum range-gate distances, respectively, for the given altitude region. To improve the vertical resolution for SBL measurements, the elevation angle, the number of range gates (

*R*

_{max}−

*R*

_{min}) and the range-resolution Δ

*R*for a given altitude region must be made as small as possible. Therefore, for new lidar-derived profiles of the longitudinal velocity statistics, we restrict the analysis domain to elevation angles from 0° to 4° and select

*R*

_{max}−

*R*

_{min}= 360 m for a maximum vertical resolution of Δ

*z*

_{turb}≈ 30 m. An example of the radial or longitudinal velocity structure functions

*D*(

_{u}*s*), the corrected structure functions, the best-fit von Kármán model, and the best-fit lidar structure function model for the same altitude and range region as in Fig. 5 is shown in Fig. 6. The correction for the estimation error of the radial velocity is based on the temporal covariance of the velocity differences between range-gate pairs (Frehlich and Cornman 2002) and

*σ*≈ 0.047 m s

_{e}^{−1}. The longitudinal turbulence parameters are approximately equal to the transverse parameters of Fig. 5, indicating approximately isotropic horizontal velocity statistics. The estimate of

*ε*= 1.32 10

_{u}^{−5}m

^{2}s

^{−3}for the spatially filtered data is noticeable lower than the estimate

*ε*= 1.72 10

_{u}^{−5}m

^{2}s

^{−3}from the raw data indicating more sensitivity to the large-scale features than the azimuth technique. However, the spatially filtered estimates are in good agreement for both techniques, which may provide more accurate information on isotropy, especially for

*ε*. A shorter pulse length and shorter processing range gate is required to resolve the turbulent structures. However, since the estimation error is inversely proportional to the pulse width, the correction for the estimation error then becomes more difficult in the low-turbulence regions because a larger fraction of the structure functions at small spacings is produced from the estimation error contribution (see Fig. 5 in the current paper and Frehlich and Cornman 2002).

## 4. Measurements of atmospheric profiles

Lidar-derived profiles of wind speed and turbulence are produced by processing the structure functions using data contained in a selected 10-m height interval as shown in Figs. 5 and 6 (see also Frehlich et al. 2006). The effective vertical resolution is then approximately 35 m for the longitudinal method and 20 m for the azimuth method. An example of 1-s TLS data (0.83-m vertical resolution) and lidar-derived profiles using the azimuth structure function technique is shown in Fig. 7 for the time period 0523–0540 UTC 7 June 2006. The roughness sublayer is clearly visible as enhanced turbulence and enhanced relative humidity in the first 50 m. Above the roughness sublayer is a region of lower turbulence (lidar derived *ε* ≈ 10 ^{−5} m^{2} s^{−3}) up to the mixing height *H* ≈ 120 m defined as the large drop in small-scale turbulence *ε* at the maximum gradient in potential temperature (see Balsley et al. 2006). These conditions are challenging for remote sensing estimates of turbulence because the velocity estimates must have a very small estimation error. The lidar-derived profiles of *ε* from the raw data and with spatial filtering are approximately equal for most of the profile. However, the length scale *L*_{0} and standard deviation *σ* are sensitive to the spatial filtering. Note the large variations or “small-scale intermittency” of the high-resolution TLS estimates of *ε* and *C*^{2}_{T}, which is typical of fully developed turbulence and is dominated by the intrinsic atmospheric processes and not the estimation error of approximately 10%–15% (Frehlich et al. 2004). There are also large difference between the TLS and lidar-derived *ε* profiles, which reflects the large variability of the small-scale turbulence in the stable boundary layer, which has been observed previously with sodar and frequency-modulated continuous-wave (FMCW) radar (Gossard and Frisch 1987; Eaton et al. 1995).

A comparison of the turbulence statistics for the transverse velocity (i.e., the azimuth structure functions) and the longitudinal velocity (i.e., the radial velocity structure functions) is shown in Fig. 8 for the raw lidar velocity estimates and the spatially filtered velocity estimates. The estimates of *ε* agree in the height interval from 50 to 120 m, which indicates a region of approximately isotropic small-scale turbulence. This interval also has similar estimates of *σ* for the raw and filtered lidar data. The length scales have larger variability, indicating the difficulty in defining a valid length scale for these intermittent turbulent processes with larger-scale forcing or gravity waves. There are regions with approximately equal length scales for the same spatial filtering (e.g., around 80–120 m).

The lower-altitude region below 50 m is the most challenging since the turbulence outer-scale *L*_{0} is small and therefore the effects of the spatial filtering by the lidar pulse are larger. In this regime, a shorter lidar pulse is desirable as well as extending the theoretical turbulence model to include anisotropy of the velocity field (e.g., Lenschow and Kristensen 1988). Variations in surface altitude and surface roughness from the suburban location (see Fig. 4) can also introduce anisotropy in the velocity statistics.

## 5. Summary and discussion

New measurements of atmospheric profiles are produced using in situ sensors and Doppler lidar processing algorithms with a focus on studies of the stable boundary layer (SBL). The low-turbulence conditions in and above the SBL are challenging conditions for high-rate turbulence measurements because a low noise floor with no interference is required. These stably stratified regions are ideally suited for investigations of turbulence processes since they contain many regions of approximately constant gradient of velocity and temperature (Figs. 2 and 3), which are the subject of many theoretical studies. For many of the regions investigated in this preliminary study, the temperature gradients are statistically more reliable than the velocity gradients because they typically have weaker wave contributions (Figs. 2 and 3). The small-scale velocity and temperature turbulence consists of regions of stably stratified turbulence with a low Froude number and a large inertial range scaling ( *f* ^{−5/3}) both below and above the Ozmidov frequency *f*_{OZ} (Riley and de Bruyn Kops 2003; Lindborg 2006). The gradient Richardson number is on the order of unity, which supports the hypothesis that turbulence-scaling laws are valid when the Richardson number is larger than 0.25 (Majda and Shefter 1998; Strang and Fernando 2001; Balsley et al. 2003, 2006) and are important for NWP model parameterization in stable conditions (Sorbjan 2006; Esau and Zilitinkevich 2006; Mauritsen and Svensson 2007). Reliable estimates of the energy dissipation rate *ε* and the temperature structure constant *C*^{2}_{T} can be produced from the one-dimensional spectra assuming universal scaling of high Reynolds number turbulence and assuming Taylor’s frozen assumption [see Eqs. (1) and (7)]. However, improvements in the sensors and electronics to remove interference and reduce the detector noise are required to produce accurate direct estimates of the small-scale statistics like *ε* [see Fig. 1 and Eq. (5)]. These in situ sensors can provide high-quality turbulence and local mean quantities and statistical accuracy can be improved using multiple sensors and a slower ascent/descent rate. In addition, the statistical description of small-scale intermittency (Frehlich et al. 2004) is required for better error analysis to fully evaluate theoretical predictions (Riley and de Bruyn Kops 2003; Lindborg 2006) of stably stratified turbulence, which requires many key turbulence parameters (e.g., *N*, *F _{h}*,

*R*,

_{e}*R*,

_{b}*R*,

_{i}*L*

_{OZ},

*ε*, and

*C*

^{2}

_{T}). Since the stably stratified free atmosphere above the boundary layer may be predominantly populated with the decay stage of various instabilities (Kelvin–Helmholtz, Taylor–Green, propagating and ducted wave instabilities, etc), fundamental turbulence-scaling laws can be extracted from this type of data and the results may be valid for all similar dimensionless parameter regimes such as the upper troposphere and stratosphere (Lindborg 2006). There are still fundamental questions regarding the definition of local gradients since waves can modulate both the background fields and the small-scale turbulence, especially when convective instabilities are produced (Meillier et al. 2008).

Lidar-derived profiles of wind speed, direction, and turbulence were produced as a spatial average over a suburban location using two methods: the azimuth structure function and the radial structure function. The azimuth technique senses the velocity component transverse to the structure function separation variable, has the highest vertical resolution, and is most sensitive to the smaller scales of turbulence since the lidar sensing volume is much smaller in the direction transverse to the lidar beam than the radial or longitudinal direction (see Figs. 5 and 6). Turbulent total kinetic energy (TKE) is defined as the variance of all three of the turbulent velocity components, which usually employs a time average or space average to separate the turbulent fluctuations from the background mean quantities and has been related to the longitudinal velocity variance *σ*^{2}_{u} by Banta et al. (2006). Spatially average lidar-derived turbulence profiles were determined with the raw lidar radial velocity estimates and with a spatial filter produced from a higher-order polynomial fit to the radial velocity as a function of the lidar beam azimuth for each range gate and fixed elevation angle. The estimates of energy dissipation rate *ε* from the azimuth structure function method were similar over the altitude region of well-defined turbulence for both the raw data and the spatially filtered data (see Figs. 5 and 7). This implies that the exact shape of the turbulence model used in the best fit has little effect on the estimates of *ε*. More high-quality verification data that samples the same space–time region is required to verify this hypothesis. The radial structure function method is not as robust for these conditions because the effects of the smoothing of the velocity field by the lidar pulse in the radial direction is more extreme as shown in Fig. 6. A shorter pulse is required to better resolve the turbulent eddies in the inertial range. However, even with the large corrections for the pulse smoothing, the estimates of *ε* from the radial structure function method are in good agreement with the azimuth structure function method in the upper region of the SBL (i.e., from 50 to 120 m, for both raw data and filtered data; Fig. 8). The estimates of velocity standard deviation *σ* (related to TKE) and the turbulence length scale *L*_{0} are not well defined for the higher altitudes above the SBL because of larger-scale atmospheric processes. A spatial filter that removes the larger-scale radial velocity components does produce well-behaved profiles but the results are a function of the amount of the spatial filtering. This is a fundamental issue, which is critical for in situ tower measurements (Vickers and Mahrt 2003) and now lidar-derived turbulence algorithms (Banta et al. 2006), especially for the SBL. Another challenge is the characterization of the low-turbulence regions above the SBL where accurate corrections for the estimation error of the radial velocities is more important (see Figs. 5 and 6 in the current paper; Frehlich 2001a; Frehlich and Cornman 2002; Frehlich et al. 2006). Improvements in lidar design and processing algorithms should resolve the complex features of the roughness sublayer and the low-turbulence area of the stably stratified free atmosphere (see Figs. 7 and 8).

The high-resolution TLS data displays large vertical variability in the small-scale turbulence statistics *ε* and *C*^{2}_{T} compared with the spatially averaged lidar-derived profiles. However, the location of the mixing height is still clearly observed as a drop in *ε* at approximately 120 m. The lidar-derived profiles of *ε* also indicate a noticeable drop in *ε* at the same altitude. More research is required to provide the most robust statistical description of the atmospheric turbulence profiles for the SBL. For NWP applications, the spatially averaged lidar profiles are more representative of the inherent spatial average of the NWP model and therefore provide smaller total observation error than rawinsondes and in situ data, which is a critical requirement for advanced data assimilation (Frehlich 2001b, 2006; Frehlich and Sharman 2004). This may be most critical for NWP model parameterization, which are essentially based on gradients of the model variables or Richardson number (Sorbjan 2006; Mauritsen and Svensson 2007) and better parameterizations may require data with a similar spatial average as the NWP model. Recent and future advances in high-resolution remote measurements of temperature with Raman lidar (Behrendt et al. 2004) may produce spatially averaged profiles of mean and turbulent statistics of the temperature field using similar algorithms as for the velocity turbulence profiles from Doppler lidar. Will the Richardson number profile from in situ data agree with a spatially averaged profile from a Doppler lidar and a Raman temperature lidar? These fundamental questions can be answered with improvements in both in situ and remote sensing techniques.

This work was supported by the National Science Foundation (NSF) under Grants ATM-0522004 and ATM-0335205 and by the Army Research Office (ARO) under Grant W911NF-06-1-0256. We would like to acknowledge the support of Walter Bach, program manager. The authors acknowledge Ben Balsley for useful discussions and field campaign support. The valuable feedback of Jim Riley is also acknowledged. Mark Vercauteren, Steve Hannon, and Gerard Pelk of Coherent Technologies, Inc., provided the Doppler lidar data.

## REFERENCES

Azizyan, G. V., , Koprov B. M. , , and Kalugin V. V. , 1989: Spectra of turbulent velocity fluctuations in the viscous range.

,*Izv. Atmos. Oceanic Phys.***25****,**935–938.Baas, P., , Steeneveld G. J. , , van de Wiel B. J. H. , , and Holtslag A. A. M. , 2006: Exploring self-correlation in flux–gradient relationships for stably stratified conditions.

,*J. Atmos. Sci.***63****,**3045–3054.Balsley, B. B., , Jensen M. L. , , and Frehlich R. G. , 1998: The use of state-of-the-art kites for profiling the lower atmosphere.

,*Bound.-Layer Meteor.***87****,**1–25.Balsley, B. B., , Frehlich R. G. , , Jensen M. L. , , Meillier Y. , , and Muschinski A. , 2003: Extreme gradients in the nighttime boundary layer: Structure, evolution, and potential causes.

,*J. Atmos. Sci.***60****,**2496–2508.Balsley, B. B., , Frehlich R. G. , , Jensen M. L. , , and Meillier Y. , 2006: High-resolution in situ profiling through the stable boundary layer: Examination of the SBL top in terms of minimum shear, maximum stratification, and turbulence decrease.

,*J. Atmos. Sci.***63****,**1291–1307.Banakh, V. A., , and Smalikho I. N. , 1997: Estimation of the turbulence energy dissipation rate from pulsed Doppler lidar data.

,*Atmos. Oceanic Opt.***10****,**957–965.Banta, R. M., , Pichugina Y. L. , , and Newsom R. K. , 2003: Relationship between low-level jet properties and turbulence kinetic energy in the nocturnal stable boundary layer.

,*J. Atmos. Sci.***60****,**2549–2555.Banta, R. M., , Pichugina Y. L. , , and Brewer W. A. , 2006: Turbulent velocity-variance in the stable boundary layer generated by a nocturnal low-level jet.

,*J. Atmos. Sci.***63****,**2700–2719.Behrendt, A., , Nakamura T. , , and Tsuda T. , 2004: Combined temperature lidar for measurements in the troposphere, stratosphere, and mesosphere.

,*Appl. Opt.***43****,**2930–2939.Bertin, F., , Barat J. , , and Wilson R. , 1997: Energy dissipation rates, eddy diffusivity, and the Prandtl number: An in situ experimental approach and its consequences on radar estimate of turbulent parameters.

,*Radio Sci.***32****,**791–804.Browning, K. A., , and Wexler R. , 1968: The determination of kinematic properties of a wind field using Doppler radar.

,*J. Appl. Meteor.***7****,**105–113.Chan, P. W., 2006: Generation of eddy dissipation rate map at the Hong Kong International Airport based on Doppler lidar data. Preprints,

*12th Conf. on Aviation, Range, and Aerospace Meteorology,*Atlanta, GA, Amer. Meteor. Soc., P5.6. [Available online at http://ams.confex.com/ams/pdfpapers/100309.pdf.].Davies, F., , Collier C. G. , , Pearson G. N. , , and Bozier K. E. , 2004: Doppler lidar measurements of turbulent structure function over an urban area.

,*J. Atmos. Oceanic Technol.***21****,**753–761.Eaton, F. D., , McLaughlin S. A. , , and Hines J. R. , 1995: A new frequency-modulated continuous wave radar for studying planetary boundary layer morphology.

,*Radio Sci.***30****,**75–88.Esau, I. N., , and Zilitinkevich S. S. , 2006: Universal dependencies between turbulent and mean flow parameters in stably and neutrally stratified planetary boundary layers.

,*Nonlinear Processes Geophys.***13****,**135–144.Frehlich, R., 1997: Effects of wind turbulence on coherent Doppler lidar performance.

,*J. Atmos. Oceanic Technol.***14****,**54–75.Frehlich, R., 2001a: Estimation of velocity error for Doppler lidar measurements.

,*J. Atmos. Oceanic Technol.***18****,**1628–1639.Frehlich, R., 2001b: Errors for space-based Doppler lidar wind measurements: Definition, performance, and verification.

,*J. Atmos. Oceanic Technol.***18****,**1749–1772.Frehlich, R., 2006: Adaptive data assimilation including the effect of spatial variations in observation error.

,*Quart. J. Roy. Meteor. Soc.***132****,**1225–1257.Frehlich, R., , and Cornman L. , 2002: Estimating spatial velocity statistics with coherent Doppler lidar.

,*J. Atmos. Oceanic Technol.***19****,**355–366.Frehlich, R., , and Sharman R. , 2004: Estimates of turbulence from numerical weather prediction model output with applications to turbulence diagnosis and data assimilation.

,*Mon. Wea. Rev.***132****,**2308–2324.Frehlich, R., , Hannon S. , , and Henderson S. , 1994: Performance of a 2-

*μ*m coherent Doppler lidar for wind measurements.,*J. Atmos. Oceanic Technol.***11****,**1517–1528.Frehlich, R., , Hannon S. , , and Henderson S. , 1997: Coherent Doppler lidar measurements of winds in the weak signal regime.

,*Appl. Opt.***36****,**3491–3499.Frehlich, R., , Hannon S. , , and Henderson S. , 1998: Coherent Doppler lidar measurements of wind field statistics.

,*Bound.-Layer Meteor.***86****,**233–256.Frehlich, R., , Meillier Y. , , Jensen M. L. , , and Balsley B. , 2003: Turbulence measurements with the CIRES tethered lifting system during CASES-99: Calibration and spectral analysis of temperature and velocity.

,*J. Atmos. Sci.***60****,**2487–2495.Frehlich, R., , Meillier Y. , , Jensen M. L. , , and Balsley B. , 2004: A statistical description of small-scale turbulence in the low-level nocturnal jet.

,*J. Atmos. Sci.***61****,**1079–1085.Frehlich, R., , Meillier Y. , , Jensen M. L. , , Balsley B. , , and Sharman R. , 2006: Measurements of boundary layer profiles in an urban environment.

,*J. Appl. Meteor. Climatol.***45****,**821–837.Fritts, D. C., 2005: Conditions contributing to adverse loading of wind turbines in the nocturnal boundary layer: Final report. National Renewable Energy Laboratory Tech. Rep. NREL/SR-500–37809, 44 pp.

Gibson-Wilde, D., , Werne J. , , Fritts D. , , and Hill R. , 2000: Direct numerical simulation of VHF radar measurements of turbulence in the mesosphere.

,*Radio Sci.***35****,**783–798.Gossard, E. E., , and Frisch A. S. , 1987: Relationship of the variances of temperature and velocity to atmospheric static stability—Application to radar and acoustic sounding.

,*J. Climate Appl. Meteor.***26****,**1021–1036.Grund, C. J., , Banta R. M. , , George J. L. , , Howell J. N. , , Post M. J. , , Richter R. A. , , and Weickmann A. M. , 2001: High-resolution Doppler lidar for boundary layer and cloud research.

,*J. Atmos. Oceanic Technol.***18****,**376–393.Hebert, D. A., , and de Bruyn Kops S. M. , 2006: Relationship between vertical shear rate and kinetic energy dissipation rate in stably stratified flows.

,*Geophys. Res. Lett.***33****.**L06602, doi:10.1029/2005GL025071.Henderson, S. W., , Hale C. P. , , Magee J. R. , , Kavaya M. J. , , and Huffaker A. V. , 1991: Eye-safe coherent laser radar system at 2.1

*μ*m using Tm,Ho:YAG lasers.,*Opt. Lett.***16****,**773–775.Henderson, S. W., , Suni P. J. M. , , Hale C. P. , , Hannon S. M. , , Magee J. R. , , Bruns D. L. , , and Yuen E. H. , 1993: Coherent laser radar at 2

*μ*m using solid-state lasers.,*IEEE Trans. Geosci. Remote Sens.***31****,**4–15.Hill, R. J., 1978: Models of the scalar spectrum for turbulent advection.

,*J. Fluid Mech.***88****,**541–562.Hill, R. J., 1996: Corrections to Taylor’s frozen turbulence approximation.

,*Atmos. Res.***40****,**153–175.Kaimal, J. C., , and Finnigan J. J. , 1994:

*Atmospheric Boundary Layer Flows: Their Structure and Measurement*. Oxford University Press, 302 pp.Kelley, N. D., , Jonkman B. J. , , and Scott G. N. , 2006: The Great Plains turbulence environment: Its origins, impact, and simulation. National Renewable Energy Laboratory Tech. Rep. NREL/CP-500–40176, 21 pp.

Klipp, C. L., , and Mahrt L. , 2004: Flux-gradient relationship, self-correlation and intermittency in the stable boundary layer.

,*Quart. J. Roy. Meteor. Soc.***130****,**2087–2103.Kolmogorov, A. N., 1962: A refinement of previous hypotheses concerning the local structure of turbulence in a viscous incompressible fluid at high Reynolds number.

,*J. Fluid Mech.***13****,**82–85.Lenschow, D. H., , and Kristensen L. , 1988: Applications of dual aircraft formation flights.

,*J. Atmos. Oceanic Technol.***5****,**715–726.Lenschow, D. H., , and Sun J. , 2007: The spectral composition of fluxes and variances over land and sea out to the mesoscale.

,*Bound.-Layer Meteor.***125****,**63–84. doi:10.1007/s10546-007-9191-8.Lindborg, E., 2006: The energy cascade in a strongly stratified fluid.

,*J. Fluid Mech.***550****,**207–242.Lothon, M., , Lenschow D. H. , , and Mayor S. D. , 2006: Coherence and scale of vertical velocity in the convective boundary layer from a Doppler lidar.

,*Bound.-Layer Meteor.***121****,**521–536.Mahrt, L., 1998: Flux sampling errors for aircraft and towers.

,*J. Atmos. Oceanic Technol.***15****,**416–429.Majda, A. J., , and Shefter M. G. , 1998: Elementary stratified flows with instability at large Richardson number.

,*J. Fluid Mech.***376****,**319–350.Mauritsen, T., , and Svensson G. , 2007: Observations of stably stratified shear-driven atmospheric turbulence at low and high Richardson numbers.

,*J. Atmos. Sci.***64****,**645–655.Meillier, Y. P., , Frehlich R. G. , , Jones R. M. , , and Balsley B. B. , 2008: Modulation of small-scale turbulence by ducted gravity waves in the nocturnal boundary layer.

,*J. Atmos. Sci.***65****,**1414–1427.Monin, A. S., , and Yaglom A. M. , 1975:

*Statistical Fluid Mechanics: Mechanics of Turbulence*. Vol. 2. MIT Press, 874 pp.Muschinski, A., , Frehlich R. , , Jensen M. L. , , Hugo R. , , Hoff A. , , Eaton F. , , and Balsley B. , 2001: Fine-scale measurements of turbulence in the lower troposphere: An intercomparison between a kite- and balloon-borne, and a helicopter-borne measurement system.

,*Bound.-Layer Meteor.***98****,**219–250.Pearson, G. N., , and Rye B. J. , 1992: Frequency fidelity of a compact CO2 Doppler lidar transmitter.

,*Appl. Opt.***31****,**6475–6484.Pearson, G. N., , and Collier C. G. , 1999: A pulsed coherent CO2 lidar for boundary-layer meteorology.

,*Quart. J. Roy. Meteor. Soc.***125****,**2703–2721.Piper, M., , and Lundquist J. K. , 2004: Surface layer turbulence measurements during a frontal passage.

,*J. Atmos. Sci.***61****,**1768–1780.Post, M. J., , and Cupp R. E. , 1990: Optimizing a pulsed Doppler lidar.

,*Appl. Opt.***29****,**4145–4158.Poulos, G. S., and Coauthors, 2002: CASES-99: A comprehensive investigation of the stable nocturnal boundary layer.

,*Bull. Amer. Meteor. Soc.***83****,**555–581.Riley, J. J., , and de Bruyn Kops S. M. , 2003: Dynamics of turbulence strongly influenced by buoyancy.

,*Phys. Fluids***15****,**2047–2059.Roth, M., 2000: Review of atmospheric turbulence over cities.

,*Quart. J. Roy. Meteor. Soc.***126****,**941–990.Saiki, E. M., , Moeng C-H. , , and Sullivan P. , 2000: Large-eddy simulation of the stably stratified planetary boundary layer.

,*Bound.-Layer Meteor.***95****,**1–30.Seibert, P., , Beyrich F. , , Gryning S. E. , , Joffre S. , , Rasmussen A. , , and Tercier P. , 2000: Review and intercomparison of operational methods for the determination of the mixing height.

,*Atmos. Environ.***34****,**1001–1027.Smalikho, I., , Köpp F. , , and Rahm S. , 2005: Measurement of atmospheric turbulence by 2-

*μ*m Doppler lidar.,*J. Atmos. Oceanic Technol.***22****,**1733–1747.Smyth, W. D., , and Moum J. N. , 2000a: Length scales of turbulence in stably stratified mixing layers.

,*Phys. Fluids***12****,**1327–1342.Smyth, W. D., , and Moum J. N. , 2000b: Anisotropy of turbulence in stably stratified mixing layers.

,*Phys. Fluids***12****,**1343–1362.Sorbjan, Z., 2006: Local structure of turbulence in stably stratified boundary layers.

,*J. Atmos. Sci.***63****,**1526–1537.Strang, E. J., , and Fernando H. J. S. , 2001: Vertical mixing and transports through a stratified shear layer.

,*J. Phys. Oceanogr.***31****,**2026–2048.Sullivan, P. P., , Horst T. W. , , Lenschow D. H. , , Moeng C-H. , , and Weil J. C. , 2003: Structure of subfilter-scale fluxes in the atmospheric surface layer with application to large-eddy simulation modeling.

,*Bound.-Layer Meteor.***482****,**101–139.Taylor, G. I., 1938: The spectrum of turbulence.

,*Proc. Roy. Soc. London***A132****,**476–490.Tjernström, M., 1993: Turbulence length scales in stably stratified free shear flow analyzed from slant aircraft profiles.

,*J. Appl. Meteor.***32****,**948–963.Vickers, D., , and Mahrt L. , 2003: The co-spectral gap and turbulent flux calculations.

,*J. Atmos. Oceanic Technol.***20****,**660–672.Werne, J., , and Fritts D. C. , 1999: Stratified shear turbulence: Evolution and statistics.

,*Geophys. Res. Lett.***26****,**439–442.Werne, J., , and Fritts D. C. , 2001: Anisotropy in a stratified shear layer.

,*Phys. Chem. Earth***26B****,**263–268.Wyngaard, J. C., , and Clifford S. F. , 1977: Taylor’s hypothesis and high-frequency turbulence spectra.

,*J. Atmos. Sci.***34****,**922–929.Zilitinkevich, S., , and Baklanov A. , 2002: Calculation of the height of the stable boundary layer in practical applications.

,*Bound.-Layer Meteor.***105****,**389–409.Zilitinkevich, S., , Esau I. , , and Baklanov A. , 2007: Further comments on the equilibrium height of neutral and stable planetary boundary layers.

,*Quart. J. Roy. Meteor. Soc.***133****,**265–271.