Model analysis of the influence of gas diffusivity in soil on CO and H2 uptake

CO and H₂ uptake by soil was studied as a diffusion process. A diffusion model was used to determine how the surface fluxes (net deposition velocities) were controlled by in‐situ microbial uptake rates and soil gas diffusivity calculated from the 3‐phase system (solid, liquid, gas) in the soil. Analytical solutions of the diffusion model assuming vertical uniformity of soil properties showed that physical properties such as air‐filled porosity and soil gas diffusivity were more important in the uptake process than in the emission process. To incorporate the distribution of in‐situ microbial uptake, we used a 2‐layer model incorporating "a microbiologically inactive layer and an active layer" as suggested from experimental results. By numerical simulation using the 2‐layer model, we estimated the effect of several factors on deposition velocities. The variations in soil gas diffusivity due to physical properties, i.e., soil moisture and air‐filled porosity, as well as to the depth of the inactive layer and in‐situ microbial uptake, were found to be important in controlling deposition velocities. This result shows that the diffusion process in soil is critically important for CO and H₂ uptake by soil, at least in soils with higher in‐situ uptake rates and/or with large variation in soil moisture. Similar uptake rates and the difference in deposition velocity between CO and H₂ may be attributable to differences in CO and H₂ molecular diffusivity. The inactive layer is resistant to diffusion and creates uptake limits in CO and H₂ by soil. The coupling of high temperature and a thick inactive layer, common in arid soils, markedly lowers net CO deposition velocity. The temperature for maximum uptake of CO changes with depth of the inactive layer.     

Observations on the use of analytical and numerical models for the description of transfer to porous surface vegetation such as lichen

Convective deposition of submicron-size aerosol to porous surface vegetation was studied by electrochemical simulation, under Reynolds and Schmidt similarity, to a rectangular array of closely-packed lichen and artificial wire roughness layers. Results, showing an approximate tenfold increase in deposition velocity over that of a flat plate placed at the same position, were compared with predictions made on the basis of various rough-surface transfer models, including those based on statistical eddy renewal, as well as with numerical solutions of the diffusion equation in statistically-renewed surface cavities. Most analytical models could be made to fit the observed data, at least for a limited range of flow velocities, but poorly known and poorly defined parameters limit their usefulness for predictive purposes; and their validity across a large variation in molecular diffusivity (or Schmidt number Sc) is generally not assured. Numerical models also depend on poorly substantiated physical assumptions but the effect of such assumptions on transfer can be calculated for a wider range of conditions than those permitting an analytical solution. This allows more direct feedback between model assumptions and calculated or observed transfer. Numerically calculated values for deposition velocity in air for Sc from 0.7 to 7000 and flow velocities from 0.2 to 5 m s^-1 are presented for different model assumptions, with values ranging from < 0.01 to > 1 cms^-1.     

Direct measurements of nitrogen oxides and ozone fluxes over grassland

Using the eddy correlation method, fluxes of nitric oxide, nitrogen dioxide, ozone, water, and sensible heat were measured at a site 20 km north of Denver, Colorado over mature crested wheat grass, 0.75 m high in late June and early July. During this period the weather was fair with no synoptic disturbances. In the early morning a well-mixed diluted urban pollution plume traversed the site, by late morning aged pollution had mixed downward into the local boundary layer, and by afternoon the air came from a relatively unpolluted area of the high plains. The mean trace gas concentrations reflect this repeated pattern of local air flow. The fluxes of the trace gases were influenced both by the variation of the means and by other factors including temperature and biological activity. Ozone fluxes were found to be always negative and proportional to the mean, with an average deposition velocity for this case of about 0.006 m s^-1. For the oxides of nitrogen this simple treatment was not appropriate. Both deposition and emission were observed, generally deposition predominated in the morning and emission in the afternoon with observed variations in the fluxes of NO_x=NO+NO₂ from –0.3 to +0.2 ppbv m s^-1.The National Center for Atmospheric Research is sponsored by the National Science Foundation     

Field measurement of small ozone fluxes to snow,wet bare soil,and lake water

Eddy-correlation measurements over snow, wet bare soil, and lake water indicate very small vertical ozone fluxes. Adjustments to the small vertical fluxes are needed to take into account the effect of mean Stefan flow associated with evaporation at the surface and the effects of correlation between density variations and vertical wind fluctuations. For snow, the residual resistance calculated for the surface is about 34 s cm^-1, indicating that the maximum deposition velocity is abut 0.03 cm s^-1. For cold bare soil well saturated with water, the surface resistance is about 10 s cm^-1 (maximum deposition velocity of about 0.1 cm s^-1). The highest resistances obtained are for transfer to the surface of Lake Michigan, yielding values near 90 s cm^-1 for resistance (0.01 cm s^-1 for deposition velocity).Work supported by the U.S. Dept. of Energy and the U.S. Environmental Protection Agency.     

Fluxes of gases and particles above a deciduous forest in wintertime

Eddy-correlation measurements of the vertical fluxes of ozone, carbon dioxide, fine particles with diameter near 0.1 m, and particulate sulfur, as well as of momentum, heat and water vapor, have been taken above a tall leafless deciduous forest in wintertime. During the experimental period of one week, ozone deposition velocities varied from about 0.1 cm s^–1 at night to more than 0.4 cm s^-1 during the daytime, with the largest variations associated primarily with changes in solar irradiation. Most of the ozone removal took place in the upper canopy. Carbon dioxide fluxes were directed upward due to respiration and exhibited a strong dependence on air temperature and solar heating. The fluxes were approximately zero at air temperatures less than 5 °C and approached 0.8 mg m^–2 s^–1 when temperatures exceeded 15 °C during the daytime. Fine-particle deposition rates were large at times, with deposition velocities near 0.8 cm s^–1 when turbulence levels were high, but fluxes directed upward were found above the canopy when the surface beneath was covered with snow. Diffusional processes seemed to dominate fine-particle transfer across quasilaminar layers and subsequent deposition to the upper canopy. Deposition velocities for particulate sulfur were highly variable and averaged to a value small in magnitude as compared to similar measurements taken previously over a pine forest in summer.     

Impact of waves on air-sea exchange of sensible heat and momentum

The impact of sea waves on sensible heat and momentum fluxes is described. The approach is based on the conservation of heat and momentum in the marine atmospheric surface layer. The experimental fact that the drag coefficient above the sea increases considerably with increasing wind speed, while the exchange coefficient for sensible heat (Stanton number) remains virtually independent of wind speed, is explained by a different balance of the turbulent and the wave-induced parts in the total fluxes of momentum and sensible heat.Organised motions induced by waves support the wave-induced stress which dominates the surface momentum flux. These organised motions do not contribute to the vertical flux of heat. The heat flux above waves is determined, in part, by the influence of waves upon the turbulence diffusivity.The turbulence diffusivity is altered by waves in an indirect way. The wave-induced stress dominates the surface flux and decays rapidly with height. Therefore the turbulent stress above waves is no longer constant with height. That changes the balance of the turbulent kinetic energy and of the dissipation rate and, hence the diffusivity.The dependence of the exchange coefficient for heat on wind speed is usually parameterized in terms of a constant Stanton number. However, an increase of the exchange coefficient with wind speed is not ruled out by field measurements and could be parametrized in terms of a constant temperature roughness length. Because of the large scatter, field data do not allow us to establish the actual dependence. The exchange coefficient for sensible heat, calculated from the model, is virtually independent of wind speed in the range of 3–10 ms^-1. For wind speeds above 10 ms^-1 an increase of 10% is obtained, which is smaller than that following from the constant roughness length parameterization.The investigation was in part supported by the Netherlands Geosciences Foundation (GOA) with financial aid from the Netherlands Organization for Scientific Research (NWO).     

均匀裸土地表土壤热通量计算方法对比及对能量闭合的影响

利用"内蒙古微气象观测蒸发试验"的观测资料,对6种地表土壤热通量计算方法(Plate Cal法、TDEC法、谐波法、热传导对流法、振幅法和相位法)进行比较,检验了6种方法在不同干湿地表状况下的适用性,并研究了6种方法计算地表土壤热通量的差异以及对地表能量闭合度的影响。结果表明:一般情况下,Plate Cal法计算的2 cm土壤热通量与观测值最接近,计算结果的均方差为6.9 W/m2。在不同干湿地表状况下,干燥和降水条件下适合使用Plate Cal法,计算结果的均方差分别为14.0 W/m2和30.1 W/m2;湿润条件下适合使用谐波法,计算结果的均方差为21.4 W/m2。6种方法计算的地表土壤热通量存在明显差别,最大相差178.6 W/m2,不同方法计算地表土壤热通量的最大差值超过25 W/m2的时次占样本的96.3%。不同方法计算地表土壤热通量的差异对地表能量闭合度的大小有明显影响,但不影响近地层能量闭合度随湍流混合增强而增大的规律。     

The Summer Surface Energy Balance of the High Antarctic Plateau

The summertime surface energy balance (SEB) at Kohnen station, situated on the high Antarctic plateau (75°00′ S, 0°04′ E, 2892m above sea level) is presented for the period of 8 January to 9 February 2002. Shortwave and longwave radiation fluxes were measured directly; the former was corrected for problems associated with the cosine response of the instrument. Sensible and latent heat fluxes were calculated using the bulk method, and eddy-correlation measurements and the modified Bowen ratio method were used to verify these calculated fluxes. The calculated sub-surface heat flux was checked by comparing calculated to measured snow temperatures. Uncertainties in the measurements and energy-balance calculations are discussed. The general meteorological conditions were not extraordinary during the period of the experiment, with a mean 2-m air temperature of −27.5°C, specific humidity of 0.52×10⁻³kg kg⁻¹ and wind speed of 4.1ms⁻¹. The experiment covered the transition period from Antarctic summer (positive net radiation) to winter (negative net radiation), and as a result the period mean net radiation, sensible heat, latent heat and sub-surface heat fluxes were small with values of −1.1, 0.0, −1.0 and 0.7 Wm⁻², respectively. Daily mean net radiation peaked on cloudy days (16 Wm⁻²) and was negative on clear-sky days (minimum of −19 W m⁻²). Daily mean sensible heat flux ranged from −8 to +10 Wm⁻², latent heat flux from −4 to 0 Wm⁻² and sub-surface heat flux from −8 to +7 Wm⁻².     

Turbulent carbon exchange in very stable conditions

Turbulent fluxes obtained using the conventional eddy covariance approach result in erratic results with large time fluctuations in extremely stable conditions. This can limit efforts to estimate components of the nocturnal energy budget and respiratory CO₂ fluxes. Well-organized fluxes that show a clear dependence on turbulent intensity were obtained when multiresolution decomposition was used to estimate turbulent exchanges. CO₂, heat and water vapour fluxes were observed at a site in the eastern Amazon basin that had been cleared for agricultural purposes. Temporal scales of the carbon transfer were determined and shown to be similar to those of latent heat, but as much as three times larger than those of sensible heat. CO₂ eddy diffusivities at the temporal scales on which most of the vertical CO₂ exchange occurs are shown to be 50 times larger than the eddy diffusivity for heat. A process associated with the vertical scale of the scalar accumulation layer is suggested to explain these different scales and turbulent diffusivities of carbon and sensible heat transfer. For an appreciable range of turbulence intensities, the observed vertical turbulent carbon exchange is insufficient to account for the locally respired CO₂ estimated independently. Evidence that shallow drainage currents may account for this is given.     

Generation and atmospheric heat exchange of coastal polynyas in the Weddell Sea

The forcing mechanisms for Antarctic coastal polynyas and the thermodynamic effects of existing polynyas are studied by means of an air-sea-ice interaction experiment in the Weddell Sea in October and November 1986.Coastal polynyas develop in close relationship to the ice motion and form most rapidly with offshore ice motion. Narrow polynyas occur frequently on the lee side of headlands and with strong curvature of the coastline. From the momentum balance of drifting sea ice, a forcing diagram is constructed, which relates ice motion to the surface-layer wind vector v _z and to the geostrophic ocean current vector c _g . In agreement with the data, wind forcing dominates when the wind speed at a height of 3 m exceeds the geostrophic current velocity by a factor of at least 33. This condition within the ocean regime of the Antarctic coastal current usually is fulfilled for wind speeds above 5 m/s at a height of 3 m.Based on a nonlinear parameter estimation technique, optimum parameters for free ice drift are calculated. Including a drift dependent geostrophic current in the ice/water drag yields a maximum of explained variance (91%) of ice velocity.The turbulent heat exchange between sea ice and polynya surfaces is derived from surface-layer wind and temperature data, from temperature changes of the air mass along its trajectory and from an application of the resistance laws for the atmospheric PBL. The turbulent heat flux averaged over all randomly distributed observations in coastal polynyas is 143 W/m². This value is significantly different over pack ice and shelf ice surfaces, where downward fluxes prevail. The large variances of turbulent fluxes can be explained by variable wind speeds and air temperatures. The heat fluxes are also affected by cloud feedback processes and vary in time due to the formation of new ice at the polynya surface.Maximum turbulent fluxes of more than 400 W/m² result from strong winds and low air temperatures. The heat exchange is similarly intense in a narrow zone close to the ice front, when under weak wind conditions, a local circulation develops and cold air associated with strong surface inversions over the shelf ice is heated above the open water.     

A Case Study Of An On-Ice Air Flow Over The Arctic Marginal Sea-Ice Zone

A case study of warm air advection over the Arctic marginalsea-ice zone is presented, based on aircraft observations with direct flux measurements carriedout in early spring, 1998. A shallow atmospheric boundary layer (ABL) was observed, which wasgradually cooling with distance downwind of the ice edge. This process was mainly connected with astrong stable stratification and downward turbulent heat fluxes of about 10–20 W m^-2, but wasalso due to radiative cooling. Two mesoscale models, one hydrostatic and the other non-hydrostatic,having different turbulence closures, were applied. Despite these fundamental differences betweenthe models, the results of both agreed well with the observed data. Various closure assumptions had amore crucial influence on the results than the differences between the models.Such an assumption was, for example,the parameterization of the surface roughness for momentum (z₀) and heat (z_T). This stronglyaffected the wind and temperature fields not only close to the surface but also within and abovethe temperature inversion layer. The best results were achieved using a formulation for z₀ that took intoaccount the form drag effect of sea-ice ridges together withz_T = 0.1z₀. The stability within theelevated inversion strongly depended on the minimum eddy diffusivity K_min. A simple ad hocparameterization seems applicable, where K_min is calculated as 0.005 timesthe neutral eddy diffusivity. Although the longwave radiative cooling was largest within the ABL, theapplication of a radiation scheme was less important there than above the ABL. This was related to theinteraction of the turbulent and radiative fluxes. To reproduce the strong inversion, it wasnecessary to use vertical and horizontal resolutions higher than those applied in most regional andlarge-scale atmospheric models.     

Vertical ozone fluxes and related deposition parameters over agricultural and forested landscapes

The spatial variability and temporal behavior of the vertical flux of ozone have been investigated from turbulence measurements collected on aircraft flight legs in the daytime period during two consecutive summer experimental field programs. The data were obtained during horizontal flight legs conducted over agricultural crops and forested land in three different regions of the eastern United States.Results from individual experimental cases and statistics derived from all cases in each region are presented. Ozone flux generally exhibited a significant height dependency. The strongest negative (downward) fluxes in the lowest-level flight legs were primarily attributed to the uptake of ozone by the surface and vegetative cover. Fluxes were near-zero in the middle of the convective boundary layer (CBL) in the afternoon period. As ozone flux was proportional to concentration, slightly stronger fluxes were found in low-level urban plume segments where ozone concentrations were 10–20 ppb higher than in the surrounding area. The derived deposition velocity showed no such bias as a function of position across the urban plume. Ozone flux differences were not apparent between the more heavily forested sections and the primarily agricultural cropland areas in these regions. During the afternoon period when no clear temporal trend was evident, means from values obtained below 0.15Z _i (Z _i being the CBL height) were -0.421 and -0.431 ppb m^-2 s^-1 for ozone flux and 0.81 and 0.82 cm s^-1 for the derived mean deposition velocity in the southeastern Pennsylvania and central Ohio areas, respectively. These experimental results for ozone provide support to a dry deposition parameterization module which computes grid-area averaged deposition velocities for use in regional-scale models.On assignment from the National Oceanic and Atmospheric Administration, U.S. Department of Commerce.     

Deposition of atmospheric particulate PCBs in suburban site of Turkey

Dry deposition and air concentration samples were collected from July 2004 to May 2005 at a suburban site in Turkey. A water surface sampler (WSS) was used to measure directly the dry deposition flux of particulate polychlorinated biphenyls (PCBs) while a high volume air sampler (HVAS) was employed to collect air samples. Particulate PCB concentrations accounted for 15% of total PCBs (gas + particle phase) at the site. The overall particulate phase PCB flux ranged from 2 to 160 ng m^− 2 d^− 1 with an average of 46.3 ± 40.6 ng m^− 2 d^− 1. Forty one PCB congeners were targeted in the samples while twenty one congeners were found to be higher than detection limits in deposition samples. Fluxes for homolog groups ranged between 0.9 (7-CBs) and 21.0 (3-CBs) ng m^− 2 d^− 1. Measured dry deposition fluxes were lower than the ones usually reported for urban sites. Average PCB dry deposition velocity, calculated using flux values and concurrently measured atmospheric concentrations, was 1.26 ± 1.86 cm s^− 1 depended on size distribution of particles, atmospheric PCB concentrations and meteorological conditions.     

Sensitivity of PBL model predictions to model design and uncertainties in environmental inputs

A simple time-dependent one-dimensional model of the planetary boundary layer (PBL) is described and used to examine the degree to which model design decisions affect model output variables. The model's sensitivity to changes in the environmental conditions is also explored. Averages of the surface fluxes, near-ground wind speeds and other PBL properties from 48 h simulations are compared to control runs. The model-calculated surface fluxes are most sensitive, in decreasing order of importance, to the vertical grid spacing, the form of closure between the surface temperature and the atmosphere, the use of vertical diffusivity smoothing, the choice of maximum time step and choice of turbulence closure scheme. These fluxes are relatively insensitive to mixing-length scaling or choice of implicit time step weighting factor. Sensitivity to changes in soil type exceeds any of the design criteria tested. The modeled fluxes are moderately sensitive to small variations in the horizontal pressure gradient, to unsteadiness in the geostrophic wind and to variations in surface roughness. They are relatively insensitive to uncertainties in local vertical velocities and small (25%) variations applied separately to soil thermal diffusivity or heat capacity. The sensitivity of the average PBL depth (Z _i ) to model and environmental changes are similar to those of surface fluxes except thatZ _i is more sensitive to changes in mixing length, albedo and imposed vertical velocity then are the surface fluxes.     

Relationship between turbulent flux and variance in the urban canopy

The relationship between fluxes and variances has been widely studied in the frame of the similarity theory. However, the dynamics of urban morphology should be highlighted in the urban canopy. Therefore, the relationship between fluxes and variances of turbulent quantities (wind velocity, temperature, carbon dioxide, and water vapor concentration) without stability parameters has been investigated with data collected from the Beijing 325-m meteorological tower in the urban canopy. Though the similarity theory is available for the statistics of wind velocity and temperature, urban land use plays an important role in the linear relationship between momentum fluxes and velocity variances, while the linear relationship between sensible heat fluxes and temperature variances is influenced at large by stratification. As for carbon dioxide and water vapor, the linear relationship between their fluxes and variances is a better way compared to the similarity theory to describe their motion, though the plots are a little scattered owing to their fugitive source.     

Parameters of Land-Surface Processes for Gobi in North-West China

Using data from the intensive observation period (May–June 2000) of the NWC-ALIEX (The Field Experiment on Interaction between Land and Atmosphere in Arid Region of North-west China), the characteristics and mechanisms of some key land-surface process over Gobi in a typical arid region of north-west China are analyzed and several parameters of land-surface processes are calculated. The weighted mean of the surface albedo over a typical arid region of the Dunhuang Gobi is calculated using the relative reflection as a weighting factor, and its value is 0.255 ± 0.021. After removing the influence of precipitation, the mean soil heat capacity over a typical arid region of the Dunhuang Gobi is 1.12 × 10⁶ J m^-3 K^-1,which is smaller than that observed in the Heihe (China) River basin Field Experiment (HEIFE). The mean soil heat diffusivity and conductivity are about half of those observed in HEIFE.     

Accuracy of the relaxed eddy-accumulation technique,evaluated using CO2 flux measurements

A system capable of measuring the fluxes of trace gases was developed. It is based on a simpler version of the eddy-accumulation technique (EA), known as the relaxed eddy-accumulation technique (REA). It accumulates air samples associated with updrafts and downdrafts at a constant flow rate in two containers for later analysis of the trace gas mean concentration. The flux integration is based on the durations of updraft and downdraft events, rather than on the vertical wind velocity (W) as is the case for EA and eddy-correlation (EC) techniques. The flux, calculated by the REA technique, is equal to the difference in the mean concentration of the trace gas of interest between the upward and downward moving eddies, multiplied by the standard deviation of the vertical wind velocity and an empirical coefficient. CO₂ fluxes measured for 162 half-hour periods over a soybean field by both EC and REA techniques showed excellent agreement (coefficient of determination,R ²=0.92). The slope (0.985) and the intercept (–0.042 mg m^–2 s^–1) were not significantly different from 1 and 0, respectively, at the 5% level; and the standard error of estimate was 0.074 mg m^–2 s^–1. It is also shown that the empirical coefficient can be calculated from either latent or sensible heat fluxes. A model describing the effect on this empirical coefficient of not sampling aroundW equal to zero is proposed.Centre for Land and Biological Resources Research Contribution No. 92-212.     

Deposition of trace substances via cloud interception on a coniferous forest at Kleiner Feldberg

A resistance model to calculate the deposition of cloud droplets on a coniferous forest and some improved parameterizations of the indispensable input parameters are described. The deposition model is adapted to the coniferous forest at the Kleiner Feldberg site and verified by the data of a drip water monitoring station below the forest canopy. The measurements of liqud water content, wind speed and trace substance compounds in cloud water of the Ground-based Cloud Experiment (GCE) at Kleiner Feldberg in 1990 are used to calculate the cloud water deposition fluxes and the deposition of trace substances via cloud water interception. The calculated deposition of trace substances via cloud water interceptions is three to six times higher than via rain during the experiment. On a long term data basis the yearly amount of cloud water deposition is 180 mm year^–1 at Kleiner Feldberg site (840 m a.s.l.) while the precipitation amount is 1030 mm year^–1. Due to higher trace substance concentrations in cloud water compared to rain the ionic deposition via cloud water interception and via precipitation were assessed to be of comparable magnitude.     

Flux Corrections Revisited

In a horizontally homogeneous atmosphere, the calculation of thevertical fluxes of mass, momentum, and heat from point measurementsare to lowest order calculated as the standard eddy covariance.Additional terms primarily related to the mean vertical velocity areadded to improve the accuracy of the calculation. A review ofprevious work is undertaken. The fluxes are reformulated beginningwith the equations of conservation of mass, momentum, and energy. Aformula for the mean vertical velocity is found that is accurate tothe ideal gas approximation. Relations are developed between fluxesat different heights. In addition to the relations for fluxes ofmass, sensible heat, and latent heat, the fluxes of momentum andturbulent kinetic energy are also considered.     

Stationarity of Turbulence in Light Winds during the Maritime Continent Thunderstorm Experiment

Non-stationary turbulence can invalidate eddy flux calculations. Two-hour longrecords of wind velocity, temperature and humidity are classified as stationaryor non-stationary based on the behaviour of the flux as a function of Reynoldsaveraging period; a number of indicators of stationarity are investigated. Thetwo-hour Maritime Continent Thunderstorm Experiment wind datasets are notcompletely stationary, as indicated by the lack of a spectral gap, but can beclassified as approximately stationary, when the mean wind speed is greaterthan 3.8 m s^-1, or the standard deviation of true wind direction is 10°, or the ratio of horizontal wind variance to wind speed is <0.25. In the stationary case the calculated friction velocity exhibits a 7% decrease on average when the Reynolds averaging period is doubled, while data classified as non-stationaryexhibit an increase of 32%. There is little non-stationary behaviour in the kinematicheat fluxes, and is independent of the non-stationarity of the friction velocity. Thekinematic heat fluxes show small decreases (around 3%) when the Reynolds averagingperiod is doubled.