A Lagrangian Model to Predict the Modification of Near-Surface Scalar Mixing Ratios and Air–Water Exchange Fluxes in Offshore Flow

A model was developed to predict the modification with fetch in offshore flow of mixing ratio, air–water exchange flux, and near-surface vertical gradients in mixing ratio of a scalar due to air–water exchange. The model was developed for planning and interpretation of air–water exchange flux measurements in the coastal zone. The Lagrangian model applies a mass balance over the internal boundary layer (IBL) using the integral depth scale approach, previously applied to development of the nocturnal boundary layer overland. Surface fluxes and vertical profiles in the surface layer were calculated using the NOAA COARE bulk algorithm and gas transfer model (e.g., Blomquist et al. 2006, Geophys Res Lett 33:1–4). IBL height was assumed proportional to the square root of fetch, and estimates of the IBL growth rate coefficient, α, were obtained by three methods: (1) calibration of the model to a large dataset of air temperature and humidity modification over Lake Ontario in 1973, (2) atmospheric soundings from the 2004 New England Air Quality Study and (3) solution of a simplified diffusion equation and an estimate of eddy diffusivity from Monin–Obukhov similarity theory (MOST). Reasonable agreement was obtained between the calibrated and MOST values of α for stable, neutral, and unstable conditions, and estimates of α agreed with previously published parametrizations that were valid for the stable IBL only. The parametrization of α provides estimates of IBL height, and the model estimates modification of scalar mixing ratio, fluxes, and near-surface gradients, under conditions of coastal offshore flow (0–50 km) over a wide range in stability.     

Modeling microclimate environments: A verification study

A numerical model is developed for simulating microclimate of plants and bare soil. The model evaluates heat, mass, momentum, and radiative fluxes in the soil-plant-atmosphere system. Its vertical domain may extend throughout the whole Planetary Boundary Layer (PBL). The model requires, either, temporal meteorological data of solar radiation, wind speed, air temperature and humidity measured over the field, or, when applied to the whole PBL, initial values of the latter three at its top. Vegetation parameters (leaf area index, photometric properties, root distribution and density) as well as soil texture, hydraulic and photometric properties are considered. The model was verified with meteorological data taken from two different climatological regions, above a bare soil and two cotton fields.For all case studies, observed and calculated values of air (except for within-canopy) and soil temperatures, wind speed, net radiation, and soil-, latent-, and sensible heat fluxes, agreed well with measurements.     

Intra-field variability of scalar flux densities across a transition between a desert and an irrigated potato field

This paper reports on measurements of sensible and latent heat and CO₂ fluxes made over an irrigated potato field, growing next to a patch of desert. The study was conducted using two eddy correlation systems. One measurement system was located within the equilibrium boundary layer 800 m downwind from the edge of the potato field. The other measurement system was mobile and was placed at various downwind positions to probe the horizontal transition of vertical scalar fluxes. Latent (LE) and sensible (H) heat fluxes, measured at 4 m above the surface, exhibited marked variations with downwind distance over the field. Only after the fetch to height ratio exceeded 75 to 1 didLE andH become invariant with downwind distance. When latent and sensible heat fluxes were measured upwind of this threshold, significant advection of humidity-deficit occurred, causing a vertical flux divergence ofH andLE.The measured fluxes of momentum, heat, and moisture were compared with predictions from a second-order closure two-dimensional atmospheric boundary layer model. There is good agreement between measurements and model predictions. A soil-plant-atmosphere model was used to examine nonlinear feedbacks between humidity-deficits, stomatal conductance and evaporation. Data interpretation with this model revealed that the advection of hot dry air did not enhance surface evaporation rates near the upwind edge of the potato field, because of negative feedbacks among stomatal conductance, humidity-deficits, andLE. This finding is consistent with results from several recent studies.     

Wind profiles,momentum fluxes and roughness lengths at Cabauw revisited

We describe the results of an experiment focusing on wind speed and momentum fluxes in the atmospheric boundary layer up to 200 m. The measurements were conducted in 1996 at the Cabauw site in the Netherlands. Momentum fluxes are measured using the K-Gill Propeller Vane. Estimates of the roughness length are derived using various techniques from the wind speed and flux measurements, and the observed differences are explained by considering the source area of the meteorological parameters. A clear rough-to-smooth transition is found in the wind speed profiles at Cabauw. The internal boundary layer reaches the lowest k-vane (20 m) only in the south-west direction where the obstacle-free fetch is about 2 km. The internal boundary layer is also reflected in the roughness lengths derived from the wind speed profiles. The lower part of the profile (< 40 m) is not in equilibrium and no reliable roughness analysis can be given. The upper part of the profile can be linked to a large-scale roughness length. Roughness lengths derived from the horizontal wind speed variance and gustiness have large footprints and therefore represent a large-scale average roughness. The drag coefficient is more locally determined but still represents a large-scale roughness length when it is measured above the local internal boundary layer. The roughness length at inhomogeneous sites can therefore be determined best from drag coefficient measurements just above the local internal boundary layers directly, or indirectly from horizontal wind speed variance or gustiness. In addition, the momentum and heat fluxes along the tower are analysed and these show significant variation with height related to stability and possibly surface heterogeneity. It appears that the dimensionless wind speed gradients scale well with local fluxes for the variety of conditions considered, including the unstable cases.     

Mesoscale Vegetation-Breeze Circulations and Their Impact on Boundary Layer Structures at Night

The impact of well watered mesoscale wheat over mid-latitude arid areas on mesoscale boundary layer structures (MBLS) and climate has been investigated in the study .using a mesoscale biophysical, meteorological model (BM) developed in the current study. The BM is composed of six modules:mesoscale atmospheric module, soil module, vegetation module, snow-atmosphere interaction module, underlying surface meteorology module and subgrid scale flux parameterization module. The six modules constitute an interacting system by supplying boundary conditions to each other.The investigation indicates that a horizontal pressure gradient associated with mesoscale perturbations in temperature and humidity is created during the day, which results from more water transpired from the vegetation canopy (VC) and evaporated from underlying wet soil. Non-classical mesoscale circulations (called as vegetation-breeze) are forced by the pressure perturbations with wind speeds about 5 m / s, flowing from the VC to the adjacent     

Increased heat fluxes near a forest edge

Summary ?Observations of sensible and latent heat flux above forest downwind of a forest edge show these fluxes to be larger than the available energy over the forest. The enhancement averages to 56 W m⁻², or 16% of the net radiation, at fetches less than 400 m, equivalent to fetch to height ratios less than 15. The enhancement of turbulent energy fluxes is explained by advection and increases with the difference in temperature and humidity of the air over the upwind area as compared to the forest. The relatively high temperature and humidity of the upwind air are not caused by high surface heat fluxes, but are explained by the relatively low aerodynamic roughness of the upwind surface. Although the heat fluxes over forest are enhanced, the momentum fluxes are almost adjusted to the underlying forest. The different behaviour of heat and momentum fluxes is explained by absorption of momentum by pressure gradients near the forest edge. It is concluded that fetch requirements to obtain accurate surface fluxes from atmospheric observations need to be more stringent for scalar fluxes as compared to momentum fluxes. Received November 23, 2001; accepted May 13, 2002     

Spectra and Cospectra of Turbulence in an Internal Boundary Layer over a Heterogeneously Irrigated Cotton Field

During the Energy Balance Experiment, patch-to-patch irrigation generated gradients in soil moisture in a north-south oriented cotton field. An internal boundary layer (IBL) developed as a result of strong horizontal advection from relatively dry upstream patches to relatively wet downstream patches associated with the prevailing northerly winds. This generated large eddies of multiple sizes, which had significant influences on the structure of turbulence in the IBL. The power spectra and cospectra of wind speed, temperature, humidity, and energy fluxes measured at two heights within the IBL are presented and used to investigate the influence of the IBL on surface layer turbulence. The spectra and cospectra were greatly enhanced by external disturbances at low frequencies. The peak frequencies of these disturbances did not change with height. The spectra and cospectra typically converged and were parallel to the Kansas spectrum at high frequencies (in the inertial subrange). A clear gap in the spectra of horizontal wind velocity existed at intermediate frequencies when the surface layer was stable. The results indicate that large eddies that originated in the upstream convective boundary layer had considerable impacts on the spectra and cospectra of surface layer turbulence. The influence of these large eddies was greater (1) when the IBL was well-developed in the near surface layer than when the IBL did not exist, (2) at higher levels than at lower levels, and (3) when the atmospheric surface layer (ASL) was unstable than when the ASL was stable. The length scales of these large eddies were consistent with the dominant scales of surface heterogeneity at the experiment site.     

An Improved Three-Dimensional Simulation of the Diurnally Varying Street-Canyon Flow

The impact of diurnal variations of the heat fluxes from building and ground surfaces on the fluid flow and air temperature distribution in street canyons is numerically investigated using the PArallelized Large-eddy Simulation Model (PALM). Simulations are performed for a 3 by 5 array of buildings with canyon aspect ratio of one for two clear summer days that differ in atmospheric instability. A detailed building energy model with a three-dimensional raster-type geometry—Temperature of Urban Facets Indoor-Outdoor Building Energy Simulator (TUF-IOBES)—provides urban surface heat fluxes as thermal boundary conditions for PALM. In vertical cross-sections at the centre of the spanwise canyon the mechanical forcing and the horizontal streamwise thermal forcing at roof level outweigh the thermal forces from the heated surfaces inside the canyon in defining the general flow pattern throughout the day. This results in a dominant canyon vortex with a persistent speed, centered at a constant height. Compared to neutral simulations, non-uniform heating of the urban canyon surfaces significantly modifies the pressure field and turbulence statistics in street canyons. Strong horizontal pressure gradients were detected in streamwise and spanwise canyons throughout the day, and which motivate larger turbulent velocity fluctuations in the horizontal directions rather than in the vertical direction. Canyon-averaged turbulent kinetic energy in all non-neutral simulations exhibits a diurnal cycle following the insolation on the ground in both spanwise and streamwise canyons, and it is larger when the canopy bottom surface is paved with darker materials and the ground surface temperature is higher as a result. Compared to uniformly distributed thermal forcing on urban surfaces, the present analysis shows that realistic non-uniform thermal forcing can result in complex local airflow patterns, as evident, for example, from the location of the vortices in horizontal planes in the spanwise canyon. This study shows the importance of three-dimensional simulations with detailed thermal boundary conditions to explore the heat and mass transport in an urban area.     

Horizontal and Vertical Co2 Advection In A Sloping Forest

A system measuring the horizontal and vertical advection was devised and installed in a sloping forest at the Vielsalm site, Belgium. The measurements showed that under stable conditions a flow regime established below the canopy: air flowed horizontally along the slope and entrained the air above the canopy vertically. This movement occurs during stable nights characterised by strongly negative net radiation. It creates negative air concentration gradients in both the vertical and horizontal directions. The advection fluxes associated with these movements are opposite and of a similar order of magnitude. This implies that the horizontal advection cannot be ignored in the carbon budget equation at night. Unfortunately, the large variability of, and considerable uncertainty about, advection fluxes does not enable one to produce estimates of the source term from these equations. Advection measurement systems should be improved in order to enable such estimates to be made. Particular attention should be paid to the estimation of the vertical velocity above the canopy and to the vertical profiles of the horizontal velocity and horizontal CO₂ gradient below the canopy.     

A Vegetated Urban Canopy Model for Meteorological and Environmental Modelling

An urban canopy model is developed for use in mesoscale meteorological and environmental modelling. The urban geometry is composed of simple homogeneous buildings characterized by the canyon aspect ratio (h/w) as well as the canyon vegetation characterized by the leaf aspect ratio (σ _l ) and leaf area density profile. Five energy exchanging surfaces (roof, wall, road, leaf, soil) are considered in the model, and energy conservation relations are applied to each component. In addition, the temperature and specific humidity of canopy air are predicted without the assumption of thermal equilibrium. For radiative transfer within the canyon, multiple reflections for shortwave radiation and one reflection for longwave radiation are considered, while the shadowing and absorption of radiation due to the canyon vegetation are computed by using the transmissivity and the leaf area density profile function. The model is evaluated using field measurements in Vancouver, British Columbia and Marseille, France. Results show that the model quite well simulates the observations of surface temperatures, canopy air temperature and specific humidity, momentum flux, net radiation, and energy partitioning into turbulent fluxes and storage heat flux. Sensitivity tests show that the canyon vegetation has a large influence not only on surface temperatures but also on the partitioning of sensible and latent heat fluxes. In addition, the surface energy balance can be affected by soil moisture content and leaf area index as well as the fraction of vegetation. These results suggest that a proper parameterization of the canyon vegetation is prerequisite for urban modelling.     

Heat and moisture budgets from airborne measurements and high-resolution model simulations

High-resolution simulations with a mesoscale model are performed to estimate heat and moisture budgets of a well-mixed boundary layer. The model budgets are validated against energy budgets obtained from airborne measurements over heterogeneous terrain in Western Germany. Time rate of change, vertical divergence, and horizontal advection for an atmospheric column of air are estimated. Results show that the time trend of specific humidity exhibits some deficiencies, while the potential temperature trend is matched accurately. Furthermore, the simulated turbulent surface fluxes of sensible and latent heat are comparable to the measured fluxes, leading to similar values of the vertical divergence. The analysis of different horizontal model resolutions exhibits improved surface fluxes with increased resolution, a fact attributed to a reduced aggregation effect. Scale-interaction effects could be identified: while time trends and advection are strongly influenced by mesoscale forcing, the turbulent surface fluxes are mainly controlled by microscale processes.     

水稻不同生长期稻田能量收支、CO2通量模拟研究

利用国家自然科学基金重大项目“长江三角洲低层大气物理化学过程与生态系统的相互作用”野外观测资料，将1999年常熟稻田水稻四个不同生长期（插秧、拔节、抽穗和成熟）的太阳辐射、水汽压、气温和水平风速及降水作为简化生物圈（SiB2）模式的5个强迫变量，在合理设计模式参数的前提下，分别模拟了上述四个阶段的能量收支、CO₂通量。模拟结果表明:在水稻不同生长阶段的感热、潜热和下垫面热通量占净辐射的比例不同；CO₂通量的模拟结果与1996年国际稻田试验结果相似。该研究为利用自动气象站网的资料估计各站能量收支、CO₂通量、土壤水分平衡、以及土壤温度廓线提供了新证据。     

Acoustic tomography and conventional meteorological measurements over heterogeneous surfaces

Summary Conventional micro-meteorological measuring methods are not particularly suitable for the investigation of the energy exchange under heterogeneous surface conditions. To consider the influence of the different surface properties, area-covered and spatially averaged meteorological measurements in combination with highly resolved simulations are necessary. In this context, the method of acoustic travel time tomography is introduced to provide information about the horizontal temperature and wind field.To check the applicability of the tomographic method for this problem and to provide data for the model initialisation a field experiment was carried out.The tomographic system and conventional meteorological equipment were utilised inside the investigation area with an extension of several hundred meters (300×700m²), which was arranged over areas with different surface properties: grassland and bare soil.The results of the field experiment show that differences between the measuring systems (in-situ and remote sensing) exist and the heterogeneity of the underlying surface is visible in the near surface temperature and wind field at a scale, which can be resolved with highly resolved numerical models. Depending on the incoming solar radiation and the local advection regional distinctions in the air temperature and wind field as well as in the vertical sensible heat fluxes were observed.The investigations demonstrate that the sensitivity of the acoustic tomography is sufficient to verify gradients in the meteorological fields even when the horizontal differences in sensible heat flux are small.     

The Atmospheric Boundary Layer In An Arctic Wintertime On-Ice Air Flow

A warm on-ice air flow from the open water over the Arctic sea ice in the Fram Straitwas, for the first time, systematically measured on 12 March 1998 by aircraft in thelowest 3 km over a 300-km long distance. The air mass modification and the processesinvolved are discussed.Over the water, air temperature was lower than water temperature so that a convectiveboundary layer (CBL) was present as initial condition. As soon as the CBL passed theice edge, a shallow stable internal boundary layer (IBL) was formed. In the residual CBL, turbulence and pre-existing convective clouds dissolved within about 20 km. Within about the same distance, due to the transition from unstable to stable stratification, the influence of surface friction increased in the IBL and decreased above the IBL with consequent generation of a low-level jet at IBL top. The IBL was strongly stratified with respect to both temperature and wind. The wind shear was around 0.1 s^-1 so that the Richardson number in the IBL was subcritical and turbulence was generated. The IBL top grew to about 145 m over 230 km distance. The growth of the IBL was not monotonic and was influenced by (a) inhomogeneous ice surface temperatures causedby both different ice thickness and changes in the cloud conditions, and (b) leads in theice deck. At the front side of the on-ice flow, the air mass boundary between the warmair and the cold Arctic air was sharp (12 K over 10 km) at low levels and tilted withheight. Observations suggest that the stratified IBL was lifted as a slab on top of thecold air.     

Numerical simulation of internal boundary-layer development and comparison with atmospheric data

A finite-volume numerical model is employed to investigate the adaptation of the atmospheric boundary layer to a change in the underlying surface roughness, such as that existing in the transition from land to the free surface of a water body. Numerical results are validated by comparison with neutral stratification atmospheric data and compared with the internal boundary-layer (IBL) heights computed using a number of existing empirical formulae. The numerical analysis allows an extension of the fetch range in which the existing formulae, calibrated only by comparison with short fetch data, may be applied. An argument is offered that the spatial variability of the water surface roughness should be also taken into account for the IBL development over the surface of a water body.     

Conditional Analysis of an Internal Boundary Layer

Data from a convective internal boundary layer (IBL) are analyzed by focusing on the instantaneousstructure of the top of the IBL instead of the time-average structure.A conditional averaging technique is developed todiscriminate between air from above the IBL and air from below the IBL , which alternately invade some instrument levels due tosubstantial variation of the top of the convective IBL.Sensitivity to the conditional sampling criteria is examined.Inside the IBL , buoyant and mechanicalproduction and dissipation dominate the turbulent kinetic energy budget.The horizontal advection and turbulent transport terms are smaller, but not negligible. The inferred pressure correlation term is negligible.Above the IBL , buoyant production and dissipation, although weak,dominate the turbulent kinetic energy budget. Shear generation andturbulent transport are smaller but significant.     

The structure of the stably stratified internal boundary layer in offshore flow over the sea

Observations obtained mainly from a research aircraft are presented of the mean and turbulent structure of the stably stratified internal boundary layer (IBL) over the sea formed by warm air advection from land to sea. The potential temperature and humidity fields reveal the vertical extent of the IBL, for fetches out to several hundred of kilometres, geostrophic winds of 20–25 m s^–1, and potential temperature differences between undisturbed continental air and the sea surface of 7 to 17 K. The dependence of IBL depth on these external parameters is discussed in the context of the numerical results of Garratt (1987), and some discrepancies are noted.Wind observations show the development of a low-level wind maximum (wind component normal to the coast) and rotation of the wind to smaller cross-isobar flow angles. Potential temperature () profiles within the IBL reveal quite a different structure to that found in the nocturnal boundary layer (NBL) over land. Over the sea, profiles have large positive curvature with vertical gradients increasing monotonically with height; this reflects the dominance of turbulent cooling within the layer. The behaviour is consistent with known behaviour in the NBL over land where curvature becomes negative (vertical gradients of decreasing with height) as radiative cooling becomes dominant.Turbulent properties are discussed in terms of non-dimensional quantities, normalised by the surface friction velocity, as functions of normalised height using the IBL depth. Vertical profiles of these and the normalised wavelength of the spectral maximum agree well with known results for the stable boundary layer over land (Caughey et al., 1979).     

Simulation of Meteorological Fields Within and Above Urban and Rural Canopies with a Mesoscale Model

Accurate simulation of air quality at neighbourhood scales (on order of 1-km horizontal grid spacing) requires detailed meteorological fields inside the roughness sub-layer (RSL). Since the assumptions of the roughness approach, used by most of the mesoscale models, are unsatisfactory at this scale, a detailed urban and rural canopy parameterisation, called DA-SM2-U, is developed inside the Penn State/NCAR Mesoscale Model (MM5) to simulate the meteorological fields within and above the urban and rural canopies. DA-SM2-U uses the drag-force approach to represent the dynamic and turbulent effects of the buildings and vegetation, and a modified version of the soil model SM2-U, called SM2-U(3D), to represent the thermodynamic effects of the canopy elements. The turbulence length scale is also modified inside the canopies. SM2-U(3D) assesses the sensible and latent heat fluxes from rural and urban surfaces in each of the computational layers inside the canopies by considering the shadowing effect, the radiative trapping by the street canyons, and the storage heat flux by the artificial surfaces. DA-SM2-U is tested during one simulated day above the city of Philadelphia, U.S.A. It is shown that DA-SM2-U is capable of simulating the important features observed in the urban and rural RSL, as seen in the vertical profiles of the shear stress, turbulent kinetic energy budget components, eddy diffusivity, potential air temperature, and specific humidity. Within the canopies, DA-SM2-U simulates the decrease of the wind speed inside the dense canopies, the skirting of the flow around the canopy blocks, warmer air inside the vegetation canopy than above open areas during the night and conversely during the day, and constantly warmer air inside the urban canopy. The comparison with measurements shows that the surface air temperature above rural and urban areas is better simulated by DA-SM2-U than by the `standard version' of MM5.     

Influence of leaf water potential on diurnal changes in CO2 and water vapour fluxes

Mass and energy fluxes between the atmosphere and vegetation are driven by meteorological variables, and controlled by plant water status, which may change more markedly diurnally than soil water. We tested the hypothesis that integration of dynamic changes in leaf water potential may improve the simulation of CO₂ and water fluxes over a wheat canopy. Simulation of leaf water potential was integrated into a comprehensive model (the ChinaAgrosys) of heat, water and CO₂ fluxes and crop growth. Photosynthesis from individual leaves was integrated to the canopy by taking into consideration the attenuation of radiation when penetrating the canopy. Transpiration was calculated with the Shuttleworth-Wallace model in which canopy resistance was taken as a link between energy balance and physiological regulation. A revised version of the Ball-Woodrow-Berry stomatal model was applied to produce a new canopy resistance model, which was validated against measured CO₂ and water vapour fluxes over winter wheat fields in Yucheng (36°57′ N, 116°36′ E, 28 m above sea level) in the North China Plain during 1997, 2001 and 2004. Leaf water potential played an important role in causing stomatal conductance to fall at midday, which caused diurnal changes in photosynthesis and transpiration. Changes in soil water potential were less important. Inclusion of the dynamics of leaf water potential can improve the precision of the simulation of CO₂ and water vapour fluxes, especially in the afternoon under water stress conditions.