An experimental study of sulfur and NOx fluxes over grassland

Three independent sulfur sensors were used in a study of sulfur eddy fluxes to a field of wheat stubble and mixed grasses, conducted in Southern Ohio in September, 1979. Two of these sensors were modified commercial instruments; one operated with a prefilter to measure gaseous sulfur compounds and the other with a denuder system to provide submicron particulate sulfur data. The third sensor was a prototype system, used to measure total sulfur fluxes. The data obtained indicated that the deposition velocity for gaseous sulfur almost always exceeded that for particulate sulfur; average surface conductances were about 1.0 cm s^–1 for gaseous sulfur in the daytime and about 0.4 cm s^–1 for particulate sulfur. The data indicate that nighttime values were probably much lower. The total sulfur sensor provided support for these conclusions. The boundary-layer quantity ln(z ₀ /z _H )was found to be 2.75 ± 0.55, in close agreement with expectations and thus providing some assurance that the site was adequate for eddy flux studies. However, fluxes derived using a prototype NO_x sensor were widely scattered, partially as a consequence of sensor noise but also possibly because of the effects of nearby sources of natural nitrogen compounds.Formerly with Argonne National Laboratory.     

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.     

An Attempt to Partition Stomatal and Non-stomatal Ozone Deposition Parts on a Short Grassland

To evaluate the damaging effect of tropospheric ozone on vegetation, it is important to evaluate the stomatal uptake of ozone. Although the stomatal flux is a dominant pathway of ozone deposition onto vegetated surfaces, non-stomatal uptake mechanisms such as soil and cuticular deposition also play a vital role, especially when the leaf area index \({LAI}< 4\). In this study, we partitioned the canopy conductance into stomatal and non-stomatal components. To calculate the stomatal conductance of water vapour for sparse vegetation, we firstly partitioned the latent heat flux into effects of transpiration and evaporation using the Shuttleworth–Wallace (SW) model. We then derived the stomatal conductance of ozone using the Penman–Monteith (PM) theory based on the similarity to water vapour conductance. The non-stomatal conductance was calculated by subtracting the stomatal conductance from the canopy conductance derived from directly-measured fluxes. Our results show that for short vegetation (LAI \(=\) 0.25) dry deposition of ozone was dominated by the non-stomatal flux, which exceeded the stomatal flux even during the daytime. At night the stomatal uptake of ozone was found to be negligibly small. In the case of vegetation with \({LAI}\approx 1\), the daytime stomatal and non-stomatal fluxes were of the same order of magnitude. These results emphasize that non-stomatal processes must be considered even in the case of well-developed vegetation where cuticular uptake is comparable in magnitude with stomatal uptake, and especially in the case of vegetated surfaces with \({LAI}< 4\) where soil uptake also has a role in ozone deposition.     

Eddy correlation measurements of CO2, latent heat,and sensible heat fluxes over a crop surface

Eddy correlation equipment was used to measure mass and energy fluxes over a soybean crop. A rapid response CO₂ sensor, a drag anemometer, a Lyman-alpha hygrometer and a fine wire thermocouple were used to sense the fluctuating quantities.Diurnal fluxes of sensible heat, latent heat and CO₂ were calculated from these data. Energy budget closure was obtained by summing the sensible and latent heat fluxes determined by eddy correlation which balanced the sum of net radiation and soil heat flux. Peak daytime CO₂ fluxes were near 1.0 mg m^–2 (ground area) s^–1.The eddy correlation technique was also employed in this study to measure nocturnal CO₂ fluxes caused by respiration from plants, soil, and roots. These CO₂ fluxes ranged from - 0.1 to - 0.25 mg m^–2s^–1.From the data collected over mature soybeans, a relationship between CO₂ flux and photosynthetically active radiation (PAR) was developed. The crop did not appear to be light-saturated at PAR flux densities < 1800 Ei m^–2 s^–1. The light compensation point was found to be about 160 Ei m^–2 s^–1.Published as Paper No. 7402, Journal Series, Nebraska Agricultural Experiment Station. The work reported here was conducted under Nebraska Agricultural Experiment Station Project 27-003 and Regional Research Project 11–33.Post-doctoral Research Associate, Professor and Professor, respectively. Center for Agricultural Meteorology and Climatology, Institute of Agriculture and Natural Resources, University of Nebraska, Lincoln, NE 68583-0728.     

Evaluation of turbulent fluxes over Maitri,Antarctica

Turbulent fluxes have been evaluated for clear sunny days over the Indian Antarctic station, Maitri, using the basic meteorological data recorded at four levels of a 28 m tower. The data are supplemented with radiation data. The surface layer over Maitri remains thermally stratified during the hours of minimum solar insolation, the so-called nighttime period. The surface winds during this period are generally very strong resulting in high momentum fluxes. In particular, for high winds (>12 m s^–1), the temperature gradient is found to be less positive than for moderate winds (4 to 7 m s^–1). Solar insolation provided the daytime heating necessary for the diurnal variation of atmospheric stability, and hence, for the turbulent fluxes. Thus, on clear days daytime conditions are marked by upward transport of heat with reduced momentum flux, while stable nighttime conditions are marked by a downward heat flux with increased momentum fluxes.     

Surface-layer scintillation measurements of daytime sensible heat and momentum fluxes

Line-averaged measurements of the structure parameter of refractive index (C _n ² ) were made using a semiconductor laser diode scintillometer above two markedly different surfaces during hours of positive net radiation. The underlying vegetation comprised in the first instance a horizontally homogeneous, pasture sward well-supplied with water, and in the second experiment, a sparse thyme canopy in a semi-arid environment. Atmospheric stability ranged between near neutral and strongly unstable (–20). The temperature structure parameterC _T ² computed from the optical measurements over four decades from 0.001 to 2 K² m^–2/3 agreed to within 5% of those determined from temperature spectra in the inertial sub-range of frequencies. Spectra were obtained from a single fine thermocouple sensor positioned near the midway position of the 100m optical path and at the beam propagation height (1.5m).With the inclusion of cup anemometer measurements, rule-of-thumb assumptions about surface roughness, and Monin-Obukhov similarity theory, path-averaged optical scintillations allow calculation of surface fluxes of sensible heat and momentum via a simple iterative procedure. Excellent agreement was obtained between these fluxes and those measured directly by eddy correlation. For sensible heat, agreement was on average close to perfect over a measured range of 0 to 500 W m^–2 with a residual standard deviation of 30 W m^–2. Friction velocities agreed within 2% over the range 0–0.9 m s^–1 (residual standard deviation of 0.06 m s^–1). The results markedly increase the range of validation obtained in previous field experiments. The potential of this scintillation technique and its theoretical foundation are briefly discussed.     

Sulfur dioxide plume dispersion and subsequent absorption by coniferous forests: Field measurements and model calculations

A Forest SO₂ Absorption Model (ForSAM) was developed to simulate (1) SO₂ plume dispersion from an emission source, (2) subsequent SO₂ absorption by coniferous forests growing downwind from the source. There are three modules: (1) a buoyancy module, (2) a dispersion module, and (3) a foliar absorption module. These modules were used to calculate hourly abovecanopy SO₂ concentrations and in-canopy deposition velocities, as well as daily amounts of SO₂ absorbed by the forest canopy for downwind distances to 42 km. Model performance testing was done with meteorological data (including ambient SO₂ concentrations) collected at various locations downwind from a coal-burning power generator at Grand Lake in central New Brunswick, Canada. Annual SO₂ emissions from this facility amounted to about 30,000 tonnes. Calculated SO₂ concentrations were similar to those obtained in the field. Calculated SO₂ deposition velocities generally agreed with published values.Notation c air parcel cooling parameter (non-dimensional) - E foliar absorption quotient (non-dimensional) - f areal fraction of foliage free from water (non-dimensional) - f _w SO₂ content of air parcel - h height of the surface layer (m) - H height of the convective mixing layer (m) - H _stack stack height (m) - k time level - k _drag coefficient of drag on the air parcel (non-dimensional) - K _z eddy viscosity coefficient for SO₂ (m²·s^–1) - L Monin-Obukhov length scale (m) - L _A single-sided leaf area index (LAI) - n degree-of-sky cloudiness (non-dimensional) - N number of parcels released with every puff (non-dimensional) - PAR photosynthetically active radiation (W m^–2) - Q emission rate (kg s^–2) - r _b diffusive boundary-layer resistance (s m^–1) - r _c canopy resistance (s m^–1) - r _cuticle cuticular resistance (s m^–1) - r _m mesophyllic resistance (s m^–1) - r _s stomatal resistance (s m^–1) - r _exit smokestack exit radius (m) - R normally distributed random variable with mean of zero and variance of t (s) - u _* frictional velocity scale, (m s^–1) - v lateral wind vector (m s^–1) - v _d SO₂ dry deposition velocity (m s^–1) - VCD water vapour deficit (mb) - z _can mean tree height (m) - Z zenith position of the sun (deg) - environmental lapse rate (°C m^–1) - dry adiabatic lapse rate (0.00986°C m^–1) - von Kármán's constant (0.04) - _B vertical velocities initiated by buoyancy (m s^–1) - canopy extinction coefficient (non-dimensional) - ()a denotes ambient conditions - ()can denotes conditions at the top of the forest canopy - ()h denotes conditions at the top of the surface layer - ()H denotes conditions at the top of the mixed layer - ()s denotes conditions at the canopy surface - ()p denotes conditions of the air parcels     

A comparison of some aerodynamic resistance methods using measurements over cotton and grass from the 1991 California ozone deposition experiment

Measurements of dry, deposition velocities (V _d ) of O₃ (using the eddy correlation technique) made over a cotton field and senescent grass near Fresno California during July and August 1991 were used to test some dry deposition velocity models. Over the cotton field, the observed maximum daytimeV _d was about 0.8 cm s^–1 and the average nighttime value was about 0.2 cm s^–1. Over the grass, daytime values averaged about 0.2 cm s^–1 and nighttime values about 0.05 cm s^–1. Application of a site-specific model known as ADOM (Acid Deposition and Oxidant Model) over the cotton field generally overestimated the observations except for a few hours in the afternoon when the observations were underestimated The overestimation was attributed to inadequacies in the surface resistance formulation and the underestimation to uncertainties in the aerodynamic formulation. Unlike previous studies which focused on the role of surface resistance, we perform additional tests using a large variety of aerodynamic resistance formulae, in addition to those in ADOM, to determine their influence on the modelledV _d of O₃ over, cotton. Over grass, ADOM considerably overestimated the observations but showed improvement when other surface resistance formulations were applied.     

Atmospheric deposition of nitrogen, sulfur and chloride in Thessaloniki, Greece

In the present study, the wet and dry depositions of particulate NO₃⁻, SO₄²⁻, Cl⁻ and NH₄⁺ were measured using a wet/dry sampler as a surrogate surface. Gas phase compounds of nitrogen, sulfur and chloride (HNO₃, NH₃, SO₂ and HCl) were measured by an annular denuder system (ADS) equipped with a back up filter for the collection of particles with diameter ≤ 5 μm. Ambient concentrations of NO, NO₂ and SO₂ were also taken into consideration. Sampling was conducted at an urban site in the center of the city of Thessaloniki, northern Greece. The presence of the aerosol species was examined by cold/warm period and the possible compounds in dry deposits were also considered. Dry deposition fluxes were found to be well correlated with ambient particle concentrations in order to be used for the calculation of particle deposition velocity. Average particulate deposition velocities calculated were 0.36, 0.20, 0.20 and 0.10 cm s^− 1 for Cl⁻, NO₃⁻, SO₄²⁻ and NH₄⁺, respectively. Total dry deposition fluxes (gas and particles) were estimated at 3.24 kg ha^− 1 year^− 1 for chloride (HCl + p-Cl), 9.97 kg ha^− 1 year^− 1 for nitrogen oxidized (NO + NO₂ + HNO₃ + p-NO₃), 5.32 kg ha^− 1 year^− 1 for nitrogen reduced (NH₃ + p-NH₄) and 15.77 kg ha^− 1 year^− 1 for sulfur (SO₂ + p-SO₄). 70–90% total dry deposition was due to gaseous species deposition. The contribution of dry deposition to the total (wet + dry) was at the level of 60–70% for sulfur and nitrogen (oxidized and reduced), whereas dry chloride deposition contributed 35% to the total. The dry-to-wet deposition ratio of all the studied species was found to be significantly associated with the precipitation amount, with nitrogen species being better and higher correlated. Wet, dry and total depositions measured in Thessaloniki, were compared with other countries of Europe, US and Asia.     

Emissions of marine biogenic sulfur to the atmosphere of northern Europe

Measurements of DMS and other reduced sulfur compounds in surface waters have been carried out from a helicopter in the seas surrounding Scandinavia. Average summer time concentrations of DMS ranged from 70 to 150 ngS L^-1. Simultaneous measurements of biological and physical parameters revealed no correlation between DMS and phytoplankton species, species assemblages, total phytoplankton biomass, chlorophyll a, temperature, and salinity. The only exception was a correlation between DMS concentration, Chrysochromulina spp. belonging to the Prymnesiophyceae, and salinity over a narrow range of salinity in the Baltic Sea.The flux of reduced sulfur to the atmosphere in July in this region is estimated to be 120–170 gS m^-2 d^-1 from the Baltic, 240–810 in the Kattegat/Skagerrak, and 120–690 in the North Sea. Annual fluxes are roughly 100 times higher than these daily fluxes. On an annual basis, biogenic sulfur emissions from the coastal seas are negligible (<1%) compared to the anthropogenic emissions in northern Europe. However, during the summer months, the biogenic sulfur emissions from the seas surrounding the Scandinavian peninsula are estimated to be as high as 20–70% of the anthropogenic emissions in Scandinavia. This makes it of interest to incorporate the biogenic emissions in calculations of long-range transport and deposition of sulfur within the region.Other volatile sulfur species, mainly methyl mercaptan, contribute about 10% of the total flux of reduced sulfur. Estimated fluxes of CS₂ to the atmosphere ranged from 1 gS m^-2 d^-1 in the Baltic Sea to 6 gS m^-2 d^-1 in the North Sea. No emissions for H₂S or COS were detected.     

Application of an infrared carbon dioxide and humidity instrument to studies of turbulent transport

A calibration equation and some results of the field performance of an infrared instrument, which is designed to measure simultaneous fluctuations of atmospheric carbon dioxide and water vapor, are described. Field observations show that the instrument is suitable for simultaneous measurement of turbulent fluxes of carbon dioxide and water vapor in conjunction with a sonic anemometer. Measured values of carbon dioxide and water vapor fluxes show diurnal variations characterized by crop activity with respect to assimilation, respiration and evapotranspiration. Carbon dioxide is transferred downward during the daytime and upward at night, while latent heat and sensible heat are transferred in the opposite sense. The non-dimensional gradient of carbon dioxide is expressed in the following form under weak unstable conditions: _c = (1 – 16 _v )^-1/2. Here, _v is the Monin-Obukhov stability parameter including the humidity effect. This relation was originally proposed for temperature and humidity. Thus, the results indicate that the turbulent mechanisms of carbon dioxide fluctuations are similar to those of other scalar entities. This is strongly supported by the high correlation coefficient found between fluctuations of carbon dioxide and temperature or humidity in the air layer over crop fields.     

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     

Influence of Transport over a Mountain Ridge on the Chemical Composition of Marine Aerosols during the ACE-2 Hillcloud Experiment

Aerosol chemical composition and trace gas measurements were made at twolocations on the northeastern peninsula of Tenerife during the ACE-2HILLCLOUD experiment, between 28 June and 23 July 1997. Measurementswere made of coarse (#gt;2.5 m aerodynamic diameter) and fine (#lt; 2.5m) aerosol Cl^–, NO₃ ^–,SO₄ ^2–, non-sea saltSO₄ ^2– (NSSS),CH₃SO₃ ^– (MSA) andNH₄ ⁺, and gas phase dimethylsulphide (DMS), HCl,HNO₃, SO₂, CH₃COOH, HCOOH andNH₃. Size distributions were measured using a cascadeimpactor. Results show that in marine air masses NSSS and MSA wereformed via DMS oxidation, with additional NSSS present in air massescontaining a continental component. Using a Eulerian box model approachfor aerosols transported between upwind and downwind sites, a mean NSSSproduction rate of 4.36 × 10^–4 gm^–3 s^–1 was calculated for daytimeclear sky periods (highest insolation), with values for cloudy periodsduring daytime and nighttime of 3.55 × 10^–4 and2.40 × 10^–4 g m^–3s^–1, respectively. The corresponding rates for MSA were6.23 × 10^–6, 8.49 × 10^–6and 6.95 × 10^–6 g m^–3s^–1, respectively. Molar concentration ratios forMSA/NSSS were 8.7% (1.8–18.2%) and 1.9%(1.3–3.5%) in clean and polluted air masses, respectively.Reactions occurring within clouds appeared to have a greater influenceon rates of MSA production, than of NSSS, while conversely daytime gasphase reactions were more important for NSSS. For MSA, nighttimein-cloud oxidation rates exceeded rates of daytime gas phase productionvia OH oxidation of DMS. NSSS, MSA and ammonium had trimodal sizedistributions, with modes at 0.3, 4.0 and >10.0 m (NSSS andNH₄ ⁺), and 0.3, 1.5 and 4.0 m (MSA). Nosignificant production of other aerosol species was observed, with theexception of ammonium, which was formed at variable rates dependent onneutralisation of the aerosol with ammonia released from spatiallynon-uniform surface sources. Seasalt components were mainly present incoarse particles, although sub-micrometre chloride was also measured.Losses by deposition exceeded calculated expectations for all species,and were highest for the seasalt fraction and nitrate.     

Atmospheric Dry Deposition to Marble and Red Stone

This paper presents dry deposition flux and deposition velocity of atmospheric particles on white marble and red stone at Dayalbagh, a suburban site of semi arid region, which is 10 km away from the industrial sector of the Agra city where due to agricultural practices vegetation predominates. The wind speed at Agra is mostly in the range of 1–2 m s^–1. The atmospheric calm conditions at Agra in summer, monsoon, and winter seasons are 47%, 35%, and 76%, respectively. Industrial areas of the city are away from Dayalbagh and are located in the NE, E, SE, and SW sectors. The main industrial activities, which are in operation in Agra city and its outskirts, are foundry and forging industry. The other industrial activities in Agra are rubber processing, lime oxidation and pulverization, chemicals, engineering and brick refractory kilns. Dry deposition samples were collected on dry days on white marble and red stone (0.224 m × 0.224 m × 0.02 m) using surface washing method. Both slabs were fixed to an iron stand (1.5 m height) at an angle of about 80 from the horizontal and exposed for 24 h on the roof of the faculty building. The order of deposition flux on white marble is NH₄⁺ > Mg²⁺ > Ca²⁺ > Na⁺ > Cl^– > K⁺ > NO₃^– > SO₄^2– > F^– and that on red stone is NH₄⁺ > Mg²⁺ > Ca²⁺ > SO₄^2– > Na⁺ > NO₃^– > K⁺ > F^– > Cl^–. Average dry deposition flux of major ions varies from 3.4 to 128.5 M m^–2 d^–1. The sum of major cations on white marble and red stone are 516.4 and 450.4 eq m^–2 d^–1, respectively while sum of major anions are 425.3 and 400.4 eq m^–2 d^–1 on white marble and red stone, respectively. Higher deposition of all ions was observed when wind blows from NE as most of the Agra Iron foundries and Ferozabad glass industries lie in this direction. The mean values of dry deposition velocity of ions vary between 0.22 cm s^–1 to 1.49 cm s^–1. Deposition velocity for all ions is higher on white marble than red stone inspite of rougher surface of red stone as compared to white marble. This could be due to the chemical nature of white marble, which is made of dolomite and hence adds significant amount of ions by dissolution during washing. Seasonally the deposition velocity was highest in winter.     

Dispersion of spores released from an elevated line source within a wheat canopy

Turbulent dispersion of spores was studied near a source located inside a wheat canopy. Two colors of Lycopodium spores were released simultaneously at a steady rate from line sources at two heights (0.4–0.5 m and 0.7–0.8 m) in a 0.8 to 1.0 m tall crop. The number of spores of each color released was estimated by weighing the sources before and after each release. Aerial spore concentrations were measured at 2 and 4 m downwind of the sources using rotorods placed at four heights above the canopy and small suction traps at two heights inside the canopy. Concentrations near the ground were estimated from deposits on sticky glass microscope slides placed on the ground. Experiments were conducted on six different days. Friction velocities ranged from about 0.3 to 0.5 m s^–1. The average horizontal fluxes of spores were calculated as the product of the observed concentrations and average wind speeds. At a distance of 2 m downwind from the sources, more than 16 to 44% of the flux of spores released from the lower source and more than 41 to 50% of the flux of spores released from the upper source were estimated to be above the canopy. These fluxes were compared with fluxes calculated using both a K-theory model and a random-flight-fluid-parcel-trajectory simulation model. The fluxes predicted by the models were generally considerably less than the values determined experimentally.     

Characterization and Temporal Variation of Asian Dust Aerosol from a Site in the Northern Chinese Deserts

Time series for total aerosol mass and the concentrations of nineteen elements, and the mass particle-size distributions of the elements, were determined for samples collected from a site in the northern deserts of China -Zhenbeitai, one of ground sites of the Asia-Pacific Regional Aerosol Characterization Experiment (ACE-Asia). Nine dust storm (DS) events were observed during the spring of 2001, lasting a combined total of twenty-six days. Peak mass loadings greatly exceeded the average (260 g m^–3), and higher than average dust-element concentrations and wind speeds were also observed during the three-month study. Material balance calculations showed that 82% mass of the total aerosol particles could be ascribed to Asian dust, of which Al, Ca, Fe, K, Mn, Si and Ti accounted for 7%, 6%, 4%, 2%, 0.1%, 32% and 1% by weight, respectively. Modelcalculated dry deposition velocities for 7 dust-elements during dust storm periods averaged 17 cm s^–1 and ranged from 14 to 21 cm s^–1. The estimated dry depositiona of Asian dust for the spring of 2001 was 189 g m^–2, of which 85% was due to dust storms. Factor analysis indicated that 89% of the dust loading during this period was due to remote or regional transport; 11% to local or background dust.     

The prediction of dust erosion by wind: An interactive model

We have devised a partial differential equation for the prediction of dust concentration in a thin layer near the ground. In this equation, erosion (detachment), transport, deposition and source are parameterised in terms of known quantities. The interaction between a wind prediction model in the boundary layer and this equation affects the evolution of the dust concentration at the top of the surface layer. Numerical integrations are carried out for various values of source strength, ambient wind and particle size. Comparison with available data shows that the results appear very reasonable and that the model should be subjected to further development and testing.Notation (x, y, z, t) space co-ordinates and time (cm,t) - u, v components of horizontal wind speed (cm s^–1) - u _g, v_g components of the geostrophic wind (cm s^–1) - V=(u²+v²)^1/2 (cm s^–1) - (û v)= 1/(h – k) _k ^h(u, v)dz(cm s^–1) - V _* friction velocity (cm s^–1) - z ₀ roughness length (cm) - k ₁ von Karman constant =0.4 - V _d deposition velocity (cm s^–1) - V _g gravitational settling velocity (cm s^–1) - h height of inversion (cm) - k height of surface layer (cm) - potential temperature (°K) - _gr potential temperature at ground (°K) - _K potential temperature at top of surface layer (°K) - P pressure (mb) - P ₀ sfc pressure (mb) - C _p/C_v - (t)= /z lapse rate of potential temperature (°K cm^–1) - A(z) variation of wind with height in transition layer - B(z) variation of wind with height in transition layer - Cd drag coefficient - C _HO transfer coefficient for sensible heat - C dust concentration (g m^–3) - C _K dust concentration at top of surface layer (g m^–3) - D(z) variation with height of dust concentration - u, v, w turbulent fluctuations of the three velocity components (cm s^–1) - A ₁ constant coefficient of proportionality for heat flux =0.2 - Ri Richardson number - g gravitational acceleration =980 cm s^–2 - Re Reynolds number = - D _s thickness of laminar sub-layer (cm) - v molecular kinematic viscosity of air - coefficient of proportionality in source term - dummy variable - t time step (sec) - n time index in numerical equations On sabbatical leave at University of Aberdeen, Department of Engineering, September 1989–February 1990.     

A numerical modeling study of the marine boundary layer over the Gulf Stream during cold air advection

A two-dimensional mesoscale model has been developed to simulate the air flow over the Gulf Stream area where typically large gradients in surface temperature exist in the winter. Numerical simulations show that the magnitude and the maximum height of the mesoscale circulation that develops downwind of the Gulf Stream depends on both the initial geostrophic wind and the large-scale moisture. As expected, a highly convective Planetary Boundary Layer (PBL) develops over this area and it was found that the Gulf Stream plays an important role in generating the strong upward heat fluxes causing a farther seaward penetration as cold air advection takes place. Numerical results agree well with the observed surface fluxes of momentum and heat and the mesoscale variation of vertical velocities obtained using Doppler Radars for a typical cold air outbreak. Precipitation pattern predicted by the numerical model is also in agreement with the observations during the Genesis of Atlantic Lows Experiment (GALE).List of Symbols u east-west velocity [m s^–1] - v north-south velocity [m s^–1] - vertical velocity in coordinate [m s^–1] - w vertical velocity inz coordinate [m s^–1] - gq potential temperature [K] - q moisture [kg kg^–1] - scaled pressure [J kg^–1 K^–1] - U _g the east-south component of geostrophic wind [m s^–1] - V _g the north-south component of geostrophic wind [m s^–1] - vertical coordinate following terrain - x east-west spatial coordinate [m] - y north-south spatial coordinate [m] - z vertical spatial coordinate [m] - t time coordinate [s] - g gravity [m² s^–1] - E terrain height [m] - H total height considered in the model [m] - q _s saturated moisture [kg kg^–1] - p pressure [mb] - p ₀₀ reference pressure [mb] - P precipitation [kg m^–2] - vertical lapse rate for potential temperature [K km^–1] - L latent heat of condensation [J kg^–1] - C _p specific heat at constant pressure [J kg^–1 K^–1] - R gas constant for dry air [J kg^–1 K^–1] - R _v gas constant for water vapor [J kg^–1 K^–1] - f Coriolis parameter (2 sin ) [s^–1] - angular velocity of the earth [s^–1] - latitude [^o] - K _H horizontal eddy exchange coefficient [m² s^–1] - t integration time interval [s] - x grid interval distance inx coordinate [m] - y grid interval distance iny coordinate [m] - adjustable coefficient inK _H - subgrid momentum flux [m² s^–2] - subgrid potential temperature flux [m K s^–1] - subgrid moisture flux [m kg kg^–1 s^–1] - u _* friction velocity [m s^–1] - _* subgrid flux temperature [K] - q _* subgrid flux moisture [kg kg^–1] - w _* subgrid convective velocity [m s^–1] - z ₀ surface roughness [m] - L Monin stability length [m] - _s surface potential temperature [K] - k von Karman's constant (0.4) - v air kinematic viscosity coefficient [m² s^–1] - K _M subgrid vertical eddy exchange coefficient for momentum [m² s^–1] - K subgrid vertical eddy exchange coefficient for heat [m² s^–1] - K _q subgrid vertical eddy exchange coefficient for moisture [m² s^–1] - z _i the height of PBL [m] - h _s the height of surface layer [m]     

Seasonal Variation in Amount and Composition of Monoterpenes Emitted by Young Pinus pinea Trees – Implications for Emission Modeling

Current inventories of terpenes released from vegetation consider only the short-term influences of light and temperature on emissions to simulate temporal variation during the year. We studied whole canopy emissions from young Pinus pinea during a 15-month enclosure in greenhouse chambers and examined data for other long-term influences. Mean daytime emission rates strongly increased during spring, reached an annual maximum of 200 pmol m^–2 total needle area s^–1 (1.1 g g^–1 leaf dry weight h^–1) between mid June and mid August, strongly declined in fall and reached an annual minimum of 1 pmol m^–2 s^–1 (0.006 g g^–1 h^–1) between January and February. Normalization to standard temperature and light conditions did not change the annual time course of emissions, but reduced summer to winter ratio from a factor of 200 to about 45. Seasonal variation was characterized also by changes in terpene composition: among the six main compounds, three (t--ocimene, linalool, 1.8-cineol) were exclusively emitted during sunlit hours in the main vegetation period, whereas the other (limonene, -pinene, myrcene) were emitted day and night and throughout the seasons. The results suggest that different terpene sources in P. pinea foliage exist and that a great part of the annual emission course observed here results from seasonal influences on these sources. A global model to simulate plant emissions is proposed, which accounts for seasonal influences on emissions in addition to the short-term effects of temperature and light. The model is tested on field data and discussed for its general application.     

Water loss and canopy resistance of a young Sitka spruce plantation

Transpiration of a 7 m-high Stika spruce forest was investigated using measurments of net radiation, sensible heat and ground heat fluxes in an energy balance to give latent heat flux, and hence canopy resistance from the Penman-Monteith equation. Sensible heat flux was measured by the eddy-correlation method using a Fluxatron circuit.During six consecutive days of measurement in July/August, canopy resistance typically followed a decreasing trend from high values (150 s m-) at dawn to around 40 s m^–1 at midday and then returning steadily to > 100 s m^–1 at sunset. Transpiration was 3 mm day^–1 on average over the period studied and changes in the rate within the day were significantly correlated with changes in net radiation.Comparisons are drawn with published data from other forest sites and the conclusion is reached that it is imprudent to generalise about transpiration rates and canopy resistances of different species at different sites from results gathered at one or two places.