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

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

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

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

Interactions between downslope flows and a developing cold-air pool

A numerical model has been used to characterize the development of a region of enhanced cooling in an alpine valley with a width of order (Formula presented.) km, under decoupled stable conditions. The region of enhanced cooling develops largely as a region of relatively dry air which partitions the valley atmosphere dynamics into two volumes, with airflow partially trapped within the valley by a developing elevated inversion. Complex interactions between the region of enhanced cooling and the downslope flows are quantified. The cooling within the region of enhanced cooling and the elevated inversion is almost equally partitioned between radiative and dynamic effects. By the end of the simulation, the different valley atmospheric regions approach a state of thermal equilibrium with one another, though this cannot be said of the valley atmosphere and its external environment.Peer reviewe

Pollutant dispersion in a developing valley cold-air pool

Pollutants are trapped and accumulate within cold-air pools, thereby affecting air quality. A numerical model is used to quantify the role of cold-air-pooling processes in the dispersion of air pollution in a developing cold-air pool within an alpine valley under decoupled stable conditions. Results indicate that the negatively buoyant downslope flows transport and mix pollutants into the valley to depths that depend on the temperature deficit of the flow and the ambient temperature structure inside the valley. Along the slopes, pollutants are generally entrained above the cold-air pool and detrained within the cold-air pool, largely above the ground-based inversion layer. The ability of the cold-air pool to dilute pollutants is quantified. The analysis shows that the downslope flows fill the valley with air from above, which is then largely trapped within the cold-air pool, and that dilution depends on where the pollutants are emitted with respect to the positions of the top of the ground-based inversion layer and cold-air pool, and on the slope wind speeds. Over the lower part of the slopes, the cold-air-pool-averaged concentrations are proportional to the slope wind speeds where the pollutants are emitted, and diminish as the cold-air pool deepens. Pollutants emitted within the ground-based inversion layer are largely trapped there. Pollutants emitted farther up the slopes detrain within the cold-air pool above the ground-based inversion layer, although some fraction, increasing with distance from the top of the slopes, penetrates into the ground-based inversion layer.Peer reviewe

Energy- and flux-budget (EFB) turbulence closure model for the stably stratified flows. Part I: Steady-state, homogeneous regimes

We propose a new turbulence closure model based on the budget equations for the key second moments: turbulent kinetic and potential energies: TKE and TPE (comprising the turbulent total energy: TTE = TKE + TPE) and vertical turbulent fluxes of momentum and buoyancy (proportional to potential temperature). Besides the concept of TTE, we take into account the non-gradient correction to the traditional buoyancy flux formulation. The proposed model grants the existence of turbulence at any gradient Richardson number, Ri. Instead of its critical value separating - as usually assumed - the turbulent and the laminar regimes, it reveals a transition interval, 0.1< Ri <1, which separates two regimes of essentially different nature but both turbulent: strong turbulence at Ri<<1; and weak turbulence, capable of transporting momentum but much less efficient in transporting heat, at Ri>1. Predictions from this model are consistent with available data from atmospheric and lab experiments, direct numerical simulation (DNS) and large-eddy simulation (LES).Comment: 40 pages, 6 figures, Boundary-layer Meteorology, resubmitted, revised versio

Global Intermittency and Collapsing Turbulence in the Stratified Planetary Boundary Layer

Direct numerical simulation of the turbulent Ekman layer over a smooth wall is used to investigate bulk properties of a planetary boundary layer under stable stratification. Our simplified configuration depends on two non-dimensional parameters: a Richardson number characterizing the stratification and a Reynolds number characterizing the turbulence scale separation. This simplified configuration is sufficient to reproduce global intermittency, a turbulence collapse, and the decoupling of the surface from the outer region of the boundary layer. Global intermittency appears even in the absence of local perturbations at the surface; the only requirement is that large-scale structures several times wider than the boundary-layer height have enough space to develop. Analysis of the mean velocity, turbulence kinetic energy, and external intermittency is used to investigate the large-scale structures and corresponding differences between stably stratified Ekman flow and channel flow. Both configurations show a similar transition to the turbulence collapse, overshoot of turbulence kinetic energy, and spectral properties. Differences in the outer region resulting from the rotation of the system lead, however, to the generation of enstrophy in the non-turbulent patches of the Ekman flow. The coefficient of the stability correction function from Monin-Obukhov similarity theory is estimated as (Formula presented.) in agreement with atmospheric observations, theoretical considerations, and results from stably stratified channel flows. Our results demonstrate the applicability of this set-up to atmospheric problems despite the intermediate Reynolds number achieved in our simulations. © 2014 The Author(s)

Wake response to an ocean-feedback mechanism: Madeira Island case study

This discussion focused on the numerical study of a wake episode. The Weather Research and Forecasting model was used in a downscale mode. The current literature focuses the discussion on the adiabatic dynamics of atmospheric wakes. Changes in mountain height and consequently on its relation to the atmospheric inversion layer should explain the shift in wake regimes: from a 'strong-wake' to a 'weak-wake' scenario. Nevertheless, changes in SST variability can also induce similar regime shifts. Increase in evaporation, contributes to increase convection and thus to an uplift of the stratified atmospheric layer, above the critical height, with subsequent internal gravity wave activity.Comment: Under review proces