NASA Cold Land Processes Experiment (CLPX 2002/03): ground-based and near-surface meteorological observations

A short-term meteorological database has been developed for the Cold Land Processes Experiment (CLPX). This database includes meteorological observations from stations designed and deployed exclusively for CLPXas well as observations available from other sources located in the small regional study area (SRSA) in north-central Colorado. The measured weather parameters include air temperature, relative humidity, wind speed and direction, barometric pressure, short- and long-wave radiation, leaf wetness, snow depth, snow water content, snow and surface temperatures, volumetric soil-moisture content, soil temperature, precipitation, water vapor flux, carbon dioxide flux, and soil heat flux. The CLPX weather stations include 10 main meteorological towers, 1 tower within each of the nine intensive study areas (ISA) and one near the local scale observation site (LSOS); and 36 simplified towers, with one tower at each of the four corners of each of the nine ISAs, which measured a reduced set of parameters. An eddy covariance system within the North Park mesocell study area (MSA) collected a variety of additional parameters beyond the 10 standard CLPX tower components. Additional meteorological observations come from a variety of existing networks maintained by the U.S. Forest Service, U.S. Geological Survey, Natural Resource Conservation Service, and the Institute of Arctic and Alpine Research. Temporal coverage varies from station to station, but it is most concentrated during the 2002/ 03 winter season. These data are useful in local meteorological energy balance research and for model development and testing. These data can be accessed through the National Snow and Ice Data Center Web site

Feasibility of Measuring Mean Vertical Motion for Estimating Advection

Numerous recent studies calculate horizontal and vertical advection terms for budget studies of net ecosystem exchange of carbon. One potential uncertainty in such studies is the estimate of mean vertical motion. This work addresses the reliability of vertical advection estimates by contrasting the vertical motion obtained from the standard practise of measuring the vertical velocity and applying a tilt correction, to the vertical motion calculated from measurements of the horizontal divergence of the flow using a network of towers. Results are compared for three different tilt correction methods. Estimates of mean vertical motion are sensitive to the choice of tilt correction method. The short-term mean (10 to 60 minutes) vertical motion based on the horizontal divergence is more realistic compared to the estimates derived from the standard practise. The divergence shows long-term mean (days to months) sinking motion at the site, apparently due to the surface roughness change. Because all the tilt correction methods rely on the assumption that the long-term mean vertical motion is zero for a given wind direction, they fail to reproduce the vertical motion based on the divergence

Comparison of Heat and Moisture Fluxes from a Modified Soil-plant-atmosphere Model with Observations from BOREAS

This study evaluates the prediction of heat and moisture fluxes from a new land surface scheme with eddy correlation data collected at the old aspen site during the Boreal Ecosystem-Atmosphere Study (BOREAS) in 1994. The model used in this study couples a multilayer vegetation model with a soil model. Inclusion of organic material in the upper soil layer is required to adequately simulate exchange between the soil and subcanopy air. Comparisons between the model and observations are discussed to reveal model misrepresentation of some aspects of the diurnal variation of subcanopy processes. Evapotranspiratio

Exploring the possible role of small scale terrain drag on stable boundary layers over land

This paper addresses the possible role of unresolved terrain drag, relative to the turbulent drag on the development of the stable atmospheric boundary layer over land. Adding a first-order estimate for terrain drag to the turbulent drag appears to provide drag that is similar to the enhanced turbulent drag obtained with the so-called long-tail mixing functions. These functions are currently used in many operational models for weather and climate, although they lack a clear physical basis. Consequently, a simple and practical quasi-empirical parameterization of terrain drag divergence for use in large-scale models is proposed and is tested in a column mode. As an outcome, the cross-isobaric mass flow (a measure for cyclone filling) with the new scheme, using realistic turbulent drag, appears to be equal to what is found with the unphysical long-tail scheme. At the same time, the new scheme produces a much more realistic less-deep boundary layer than is obtained by using the long-tail mixing function

The Cospectral Gap and Turbulent Flux Calculations

An alternative method to Fourier analysis is discussed for studying the scale dependence of variances and covariances in atmospheric boundary layer time series. Unlike Fourier decomposition, the scale dependence based on multiresolution decomposition depends on the scale of the fluctuations and not the periodicity. An example calculation is presented in detail. Multiresolution decomposition is applied to tower datasets to study the cospectral gap scale, which is the timescale that separates turbulent and mesoscale fluxes of heat, moisture, and momentum between the atmosphere and the surface. It is desirable to partition the flux because turbulent fluxes are related to the local wind shear and temperature stratification through similarity theory, while mesoscale fluxes are not. Use of the gap timescale to calculate the eddy correlation flux removes contamination by mesoscale motions, and therefore improves similarity relationships compared to the usual approach of using a constant averaging timescale. A simple model is developed to predict the gap scale. The goal here is to develop a practical formulation based on readily available variables rather than a theory for the transporting eddy scales. The gap scale increases with height, increases with instability, and decreases sharply with increasing stability. With strong stratification and weak winds, the gap scale is on the order of a few minutes or less. Implementation of the gap approach involves calculating an eddy correlation flux using the modeled gap timescale to define the turbulent fluctuations (e.g., w′ and T′). The turbulent fluxes (e.g., w′T′) are then averaged over 1 h to reduce random sampling errors

Quality Control and Flux Sampling Problems for Tower and Aircraft Data

A series of automated tests is developed for tower and aircraft time series to identify instrumentation problems, flux sampling problems, and physically plausible but unusual situations. The automated procedures serve as a safety net for quality controlling data. A number of special flags are developed representing a variety of potential problems such as inconsistencies between different tower levels and the flux error due to fluctuations of aircraft height. The tests are implemented by specifying critical values for parameters representing each specific error. The critical values are developed empirically from experience of applying the tests to real turbulent time series. When these values are exceeded, the record is flagged for further inspection and comparison with the rest of the concurrent data. The inspection step is necessary to either verify an instrumentation problem or identify physically plausible behavior. The set of tests is applied to tower data from the Risø Air Sea Experiment and Microfronts95 and aircraft data from the Boreal Ecosystem–Atmosphere Study

Heat fluxes over weak SST heterogeneity

The spatial variability of turbulence and surface heat flux are examined for the case of small air-surface temperature difference and modest sea-surface temperature variability. As a result of nonlinearities in the bulk formula, the heterogeneity is predicted to shift the area-averaged heat flux toward more significant upward values compared to that computed from the usual bulk formula and the area-averaged air-surface temperature difference. This prediction is supported by a case study analysis of aircraft data collected over heterogeneity of the sea-surface temperature. Traditional approaches for computing the surface flux are found to be unreliable over surface heterogeneity where the scales of the transport vary horizontally. A method based on a multiresolution wavelet transform reveals the spatial variability of the flux for different scales. This information together with several additional scale-dependent indices are used to select the range of scales included in the flux computation. This approach provides more reliable estimates of the spatial variation of the flux, although uncertainties remain. Required improvements in observation strategies are discussed for quantitative evaluation of the bulk formula in conditions of surface heterogeneity

Revisiting the Local Scaling Hypothesis in Stably Stratified Atmospheric Boundary Layer Turbulence: an Integration of Field and Laboratory Measurements with Large-eddy Simulations

The `local scaling' hypothesis, first introduced by Nieuwstadt two decades ago, describes the turbulence structure of stable boundary layers in a very succinct way and is an integral part of numerous local closure-based numerical weather prediction models. However, the validity of this hypothesis under very stable conditions is a subject of on-going debate. In this work, we attempt to address this controversial issue by performing extensive analyses of turbulence data from several field campaigns, wind-tunnel experiments and large-eddy simulations. Wide range of stabilities, diverse field conditions and a comprehensive set of turbulence statistics make this study distinct

Flux-gradient relationship, self-correlation and intermittency in the stable boundary layer

The correlation between dimensionless shear ϕₘ and dimensionless height z/L, where L is the Obukhov length, for stable conditions is strongly influenced by self‐correlation for the present datasets. This effect is quite large for stronger stability but still significant for near‐neutral conditions. A conditional analysis of nocturnal stable boundary‐layer data, where ‘non‐turbulent’ parts of the record are removed, reduces the impact of nonstationarity and therefore reduces the scatter. The conditional analysis also reduces the relative importance of self‐correlation. Difficulties with estimating self‐correlation are also discussed.Keywords: Monin-Obukhov similarity, Nocturnal boundary layer, z-less similarity, CASES99Keywords: Monin-Obukhov similarity, Nocturnal boundary layer, z-less similarity, CASES9

On the Stratification of Turbulent Mixed Layers

The vertical distribution of density or temperature is studied in turbulent boundary layers stratified by either salt or temperature and driven in an annulus by a rotating screen. With rapidly growing boundary layers the buoyancy flux due to entrainment becomes sufficiently large that the boundary layer is no longer well mixed. Stratification in the boundary layer is largest near the entrainment interface, corresponding to extension of the density transition into the turbulent boundary layer. The stratification in the rest of the boundary layer is smaller and nearly independent of height. The mean boundary layer stratification and structure and integral properties of the density distribution are found to be strongly dependent on the ratio of the entrainment rate to the turbulent velocity scale