Markov chain simulations of vertical dispersion in the neutral surface layer for surface and elevated releases

Vertical dispersion in the neutral surface layer is investigated using a Markov Chain simulation procedure. The conceptual basis of the procedure is discussed and computation procedures outlined. Wind and turbulence parameterizations appropriate to the neutral surface layer are considered with emphasis on the Lagrangian time scale. Computations for a surface release are compared with field data. Good agreement is found for the variation of surface concentration and cloud height to distances 500 m downwind of the source. The functional form of the vertical concentration profile is examined and an exponential with exponent 1.6 is found to give the best fit with simulations.For elevated releases, it is demonstrated that an initial dip of the mass mean height from the simulation can be normalized for various release heights using a non-dimensionalized downwind coordinate incorporating advective wind speed and wind shear. The vertical distribution standard deviation ( _z ), as employed in Gaussian models, shows a fair degree of independence with source height but close examination reveals an optimum source height for maximum _z at a given downwind distance,x. This source height increases with downwind distance. Also the simulations indicate that vertical wind shear is more important than vertical variation of Lagrangian time scale close to the source, with a reverse effect farther downwind.     

Turbulent transport from an arctic lead: A large-eddy simulation

The upward transfer of heat from ocean to atmosphere is examined for an Arctic lead, a break in the Arctic ice which allows contact between the cold atmosphere and the relatively warm ocean. We employ a large-eddy model to compute explicitly the three-dimensional turbulent response of the atmosphere to a lead of 200 m width. The surface heat flux creates a turbulent plume of individual quasi-random eddies, not a continuous updraft, which penetrate into the stable atmosphere and transport heat upward.Maximum updraft velocities and turbulence occur downwind of the lead rather than over the lead itself, because the development time of an individual thermal eddy is longer than its transit time across the lead. The affected vertical region, while shallow over the lead itself, grows to a height of 65m at 600 m downwind of the lead; beyond that, the depth of the turbulent region decreases as the eddies weaken. The maximum vertical turbulent heat flux occurs at the downwind edge of the lead, beyond which a relative maximum extends upward into the plume. Negative surface heat flux immediately downwind of the lead creates a growing stable layer, but above that internal boundary layer the turbulent heat flux is still positive. Updraft maxima are typically 28 cm/s, but compensating downdrafts result in time-averaged vertical velocities of less than 1 cm/s in the plume. Conditional sampling separates the updraft and downdraft contributions. Formulas for the horizontal eddy development distance and for the vertical plume penetration height are presented. The relative importance of mean and turbulent transport is compared for both vertical and horizontal heat transfer: turbulence dominates the vertical heat transport whereas mean advection dominates the horizontal transport, these offsetting transports producing a quasi-stationary state.     

A test of isotropy and Taylor's hypothesis in the atmospheric boundary layer

This study attempts to determine the scales of turbulence in a high Reynolds number shear flow near which transition to isotropy occurs and the scales for which Taylor's hypothesis is applicable. The flow studied was the wind near height x ₃ = 2 m above a flat land surface. Four hot-wire anemometers were mounted in a three-dimensional array with equal separations between 1.8 m and 2 cm in three different directions. Theoretical cross-spectra were computed from the observed spectra of downwind velocity fluctuations assuming isotropy and Taylor's hypothesis. Comparison between these and the observed cross-spectra revealed that the turbulence in the flow studied was consistent with both assumptions provided k _1x3&> 20, where k ₁ is the radian wavenumber; this was the lower bound to which no departure from isotropy could be detected by the experiment. For 4 k _1x3 20, the observations are consistent with symmetry of the turbulence about the downstream direction. That part of Taylor's hypothesis relating observed frequency at a stationary sensor to the downstream wavenumber component appears to be justified within experimental error for k _1x3& > 3.     

Size-Resolved Characterisation of Soluble Ions in the Particles in the Tropospheric Plume of Masaya Volcano, Nicaragua: Origins and Plume Processing

We present the first application of a multi-stage impactor to study volcanic particle emissions to the troposphere from Masaya volcano, Nicaragua. Concentrations of soluble SO₄ ^2–,Cl^–, F^–, NO₃ ^–, K⁺, Na⁺,NH₄ ⁺, Ca²⁺ and Mg²⁺ were determined in 11 size bins from 0.07 m to >25.5 m. The near-source size distributions showed major modes at 0.5m (SO₄ ^2–, H⁺,NH₄ ⁺); 0.2 m and 5.0 m (Cl^–) and 2.0–5.0 m(F^–). K⁺ and Na⁺ mirrored the SO₄ ^2– size-resolvedconcentrations closely, suggesting that these were transported primarily asK₂SO₄ and Na₂SO₄ in acidic solution, while Mg²⁺ andCa²⁺ presented modes in both <1 m and >1 m particles. Changes in relative humidity were studied by comparing daytime (transparent plume) and night-time (condensed plume) results. Enhanced particle growth rates were observed in the night-time plume as well as preferential scavenging of soluble gases, such as HCl, by condensed water. Neutralisation of the acidic aerosol by background ammonia was observed at the crater rim and to a greater extent approximately 15 km downwind of the active crater. We report measurements of re-suspended near-source volcanic dust, which may form a component of the plume downwind. Elevated levels ofSO₄ ^2–, Cl^–, F^–,H⁺, Na⁺, K⁺ and Mg²⁺ were observed around the 10 m particle diameter in this dust. The volcanic SO₄ ^2– flux leaving the craterwas 0.07 kg s^–1.     

THE WINDFIELD AND RAINFALL DISTRIBUTION INDUCED WITHIN A SMALL VALLEY: FIELD OBSERVATIONS AND 2-D NUMERICAL MODELLING

The distribution of rainfall over small-scale topography of uniformroughness was studied both numerically and from field observations. Field data haveshown that small-scale topographic inhomogeneities (hills and valleys of afew tens of metres in height) substantially influence the rainfalldistribution over the inhomogeneity itself and also some distance downwind from it.According to comprehensive rainfall measurements carried out in a smallwatershed, windward slopes received less conventionally-measured rain than lee-sides.Moreover, on the windward slope rainfall decreased uphill, to apronounced minimum near the crest. In addition, directional raingages,situated at the bottom of both slopes revealed a pattern suggesting a reversedsurface flow on the lower slopes. This pattern persisted through a variety ofbackground conditions.To explain the observed data, a numerical simulation of both the windfieldand of drop trajectories was carried out. The windfield was generated insimulations using the Colorado State University RAMS model with a very highresolution (5 m in the horizontal). Drop trajectories were calculated usingan equation for the motion of drops in a flow field that incorporates currentexperimental results.A strong effect of small-scale topography on drop trajectories wasfound. The simulation results are encouraging and help to explain manycharacteristic features of precipitation distributions observed in fieldexperiments.     

Experimental investigation of atmospheric dispersion over the Swiss Plain — Experiment “SIESTA”

Between November 15 and 30, 1985, an international mesoscale transport experiment was performed on the Swiss Plain. Seven meteorological groups from Denmark, Germany, Italy and Switzerland measured diffusion properties of near neutral planetary boundary layers during six completely overcast days: four Bise (north-east wind) situations, one transitional situation and one west-wind situation. Diffusion was measured using SF₆ as tracer, which was released at the meteotower of the Gösgen nuclear power plant at 6 m above ground level. One hundred automatic samplers plus a mobile gas chromatograph were used to measure the concentration fields at distances up to 90 km downwind. Meteorological parameters were observed using radar-tracked constant-level balloons, tethered balloon soundings, sonic anemometers, an acoustic sounder and several meteorological ground stations, including a 110 m mast. All data were collected on a magnetic tape with free access to interested persons.The aim of the experiments was to obtain knowledge about the general nature of the turbulence, advection and atmospheric dispersion during neutral weak-wind situations over complex terrain. The collected data set will be useful for testing mesoscale transport and diffusion models. The results clearly show the channelling effect of the Jura mountains and the hilly prealpine region. An interesting result is that the SF₆ plumes showed intensive horizontal spread but only limited diffusion in the vertical direction.     

Lidar measurements of plume statistics

Surface-layer aerosol diffusion experiments have been conducted using artificial smoke plume releases at ground level over flat and homogeneously vegetated terrain at the Meppen proving grounds in the Federal Republic of Germany (1989). At fixed downwind locations in the range out to 800 m from the source, instantaneous crosswind plume profiles were detected repetitively at high spatial (1.5 m) and temporal (3 sec) intervals by use of a mini LIDAR system. The experiments were accompanied by measurement of the surface-layer mean wind and turbulence quantities by sonic anemometers. On the basis of measured crosswind concentration profiles, the following statistics were obtained: 1) Mean profile, 2) Root mean square profile, 3) Fluctuation intensities, and 4) Intermittency factors. Furthermore, some experimentally determined probability density functions (pdf's) of the fluctuations are presented. All the measured statistics are referred to a fixed and a moving frame of reference, the latter being defined as a frame of reference from which the (low frequency) plume meander is removed. Finally, the measured statistics are compared with statistics on concentration fluctuations obtained with a simple puff diffusion model (RIMPUFF) developed at Risø.     

An approximate analytical solution to the diffusion equation for short-range dispersion from a continuous ground-level source

An easily-evaluated expression for the dimensionless concentration profile (z/z ₀,/z ₀, z ₀/L) = = cu _*/kQ (or z ₀cu_*/kQ) downwind of a continuous ground-level area (or line) source in the stable surface layer is obtained by integrating the diffusion equation using the Shwetz approximation method (c = concentration, Q = source strength, k = von Kárman's constant). The analytical solution compares closely with concentration profiles obtained using a trajectory-simulation model over a useful range of heights, the important discrepancies occurring at the upper edge of the plume. The analytical solution is used to generate predictions of ground-level concentration for the Project Prairie Grass experiments; good agreement with the observations is obtained at all downwind distances (50 to 800 m).     

Flux Footprint Simulation Downwind of a Forest Edge

Surface fluxes, originating from forest patches, are commonly calculated from atmospheric flux measurements at some height above that patch using a correction for flux arising from upwind surfaces. Footprint models have been developed to calculate such a correction. These models commonly assume homogeneous turbulence, resulting in a simulated atmospheric flux equal to the average surface flux in the footprint area. However, atmospheric scalar fluxes downwind of a forest edge have been observed to exceed surface fluxes in the footprint area. Variations in atmospheric turbulence downwind of the forest edge, as simulated with an E – model, can explain enhanced atmospheric scalar fluxes. This E – model is used to calculate the footprint of atmospheric measurements downwind of a forest edge. Atmospheric fluxes appear mainly enhanced as a result of a stronger sensitivity to fluxes from the upwind surface. A sensitivity analysis shows that the fetch over forest, necessary to reach equilibrium between atmospheric fluxes and surface fluxes, tends to be longer for scalar fluxes as compared to momentum fluxes. With increasing forest density, atmospheric fluxes deviate even more strongly from surface fluxes, but over shorter fetches. It is concluded that scalar fluxes over forests are commonly affected by inhomogeneous turbulence over large fetches downwind of an edge. It is recommended to take horizontal variations in turbulence into account when the footprint is calculated for atmospheric flux measurements downwind of a forest edge. The spatially integrated footprint is recommended to describe the ratio between the atmospheric flux and the average surface flux in the footprint.     

The Effect of Wind-Turbine Wakes on Summertime US Midwest Atmospheric Wind Profiles as Observed with Ground-Based Doppler Lidar

We examine the influence of a modern multi-megawatt wind turbine on wind and turbulence profiles three rotor diameters ( $D$ D ) downwind of the turbine. Light detection and ranging (lidar) wind-profile observations were collected during summer 2011 in an operating wind farm in central Iowa at 20-m vertical intervals from 40 to 220 m above the surface. After a calibration period during which two lidars were operated next to each other, one lidar was located approximately $2D$ 2 D directly south of a wind turbine; the other lidar was moved approximately $3D$ 3 D north of the same wind turbine. Data from the two lidars during southerly flow conditions enabled the simultaneous capture of inflow and wake conditions. The inflow wind and turbulence profiles exhibit strong variability with atmospheric stability: daytime profiles are well-mixed with little shear and strong turbulence, while nighttime profiles exhibit minimal turbulence and considerable shear across the rotor disk region and above. Consistent with the observations available from other studies and with wind-tunnel and large-eddy simulation studies, measurable reductions in wake wind-speeds occur at heights spanning the wind turbine rotor (43–117 m), and turbulent quantities increase in the wake. In generalizing these results as a function of inflow wind speed, we find the wind-speed deficit in the wake is largest at hub height or just above, and the maximum deficit occurs when wind speeds are below the rated speed for the turbine. Similarly, the maximum enhancement of turbulence kinetic energy and turbulence intensity occurs at hub height, although observations at the top of the rotor disk do not allow assessment of turbulence in that region. The wind shear below turbine hub height (quantified here with the power-law coefficient) is found to be a useful parameter to identify whether a downwind lidar observes turbine wake or free-flow conditions. These field observations provide data for validating turbine-wake models and wind-tunnel observations, and for guiding assessments of the impacts of wakes on surface turbulent fluxes or surface temperatures downwind of turbines.     

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]     

A Spectrum-Independent Procedure for Correcting Eddy Fluxes Measured with Separated Sensors

We investigate flux underestimates in eddy correlation measurements that are caused by horizontal separation of the sensors. A common eddy correlation setup consists of a sonic anemometer and a humidity sensor which, because of its bulk, must be placed some distance away from the sonic path, leading to a flux loss (of latent heat). Utilizing an additional fast temperature sensor placed near the humidity sensor, we develop a procedure for correcting for this loss. The procedure simultaneously corrects the sensible heat flux for the difference between true temperature and sonic temperature. Our correction procedure, which does not depend on the shape of the cospectrum, is then compared to the widely-used procedure following Moore (1986), which assumes a cospectral model ('Kansas Model). Both correction methods are applied to data collected within the internal boundary layer over a rice paddy, downwind of arid land. Under conditions of good fetch, they were found to agree well. Under poor fetch conditions, the model-based correction tended to be too small, while the spectrum-independent combined correction was robust. The latter is thus recommended for situations where the cospectral shape can be expected to deviate from the 'Kansas shape.     

Gaseous and Particulate Oxidation Products Analysis of a Mixture of α-pinene + β-pinene/O3/Air in the Absence of Light and α-pinene + β-pinene/NOx/Air in the Presence of Natural Sunlight

The gas and particle phase reaction products of a mixture of the atmospherically important terpenes -pinene and -pinene with the atmospheric oxidants O₃ and OH/NO_x were investigated using both gas chromatography-mass spectrometry (GC-MS) and high performance liquid chromatography (HPLC) for identification and quantification of reaction products. The nighttime oxidation of a mixture of -pinene and -pinene in the presence of O₃/air, and the daytime oxidation of a mixture of -pinene + -pinene with NO_x air in the presence of natural sunlight were carried out in the University of North Carolina's large outdoor smog chamber (190 m³) located in Chatham County, North Carolina. Mass balances for gaseous and aerosol reaction products are reported over the course of the reaction. More than twenty-nine products were identified and/or quantified in this study. On average, measured gas and particle phase products accounted for 74 to 80% of the reacted -pinen/-pinene mixture carbon. Measurements show that a number of reaction products were found in both O₃ and NO_x system [pinonaldehyde, pinic acid, pinonic acid, pinalic-3-acid, 4-hydroxypinalic-3-acid, 4-oxonopinone, 1-hydroxy-nopinone, 3-hydroxy-nopinone, and nopinone]. Pinonic acid, pinic acid, pinalic-3-acid, 4-hydroxypinalic-3-acid, and 10-hydroxypinonic acid were observed in the early stage in the aerosol phase and may play an important role in the early formation of secondary aerosols.     

Numerical simulation of particle trajectories in inhomogeneous turbulence,III: Comparison of predictions with experimental data for the atmospheric surface layer

It is shown that predictions of a numerical trajectory-simulation method agree closely with the Project Prairie Grass observations of the concentrations 100 m downwind of a continuous point source of sulphur dioxide if the height (z) dependence of the Lagrangian length scale Λ _L is chosen as: whereL is the Monin-Obukhov length. The value of 0.5 for Λ _L /z in neutral conditions is consistent with the findings of Reid (1979) for the Porton experiment, and is also shown to be the best choice for simulation of an experiment in which concentration profiles were measured a short distance (< 40 m) downwind of an elevated point source of glass beads (40 μn diameter). $$\begin{gathered} \Lambda _L = 0.5z\left( {1 - 6\frac{z}{L}} \right)^{{1 \mathord{\left/ {\vphantom {1 4}} \right. \kern-\nulldelimiterspace} 4}} L< 0 \hfill \\ \Lambda _L = 0.5z/\left( {1 + 5\frac{z}{L}} \right)L > 0 \hfill \\ \end{gathered} $$      

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 numerical investigation of wind speed effects on lake-effect storms

Observations of lake-effect storms that occur over the Great Lakes region during late autumn and winter indicate a high sensitivity to ambient wind speed and direction. In this paper, a two-dimensional version of the Penn State University/National Center for Atmospheric Research (PSU/NCAR) model is used to investigate the wind speed effects on lake-effect snowstorms that occur over the Great Lakes region.Theoretical initial conditions for stability, relative humidity, wind velocity, and lake/land temperature distribution are specified. Nine different experiments are performed using wind speeds ofU=0, 2, 4,..., 16 m s^–1. The perturbation wind, temperature, and moisture fields for each experiment after 36 h of simulation are compared.It is determined that moderate (4–6 m s^–1) wind speeds result in maximum precipitation (snowfall) on the lee shore of the model lake. Weak wind speeds (0U<4 m s^–1) yield significantly higher snowfall amounts over the lake along with a spatially concentrated and intense response. Strong wind speeds (6<U16 m s^–1), yield very little, if any, significant snowfall, although significant increases in cloudiness, temperature, and perturbation wind speed occur hundreds of kilometers downwind from the lake.     

Fine-resolution regional climate model simulations of the impact of climate change on tropical cyclones near Australia

Fine-resolution regional climate simulations of tropical cyclones (TCs) are performed over the eastern Australian region. The horizontal resolution (30 km) is fine enough that a good climatological simulation of observed tropical cyclone formation is obtained using the observed tropical cyclone lower wind speed threshold (17 m s^–1). This simulation is performed without the insertion of artificial vortices (bogussing). The simulated occurrence of cyclones, measured in numbers of days of cyclone activity, is slightly greater than observed. While the model-simulated distribution of central pressures resembles that observed, simulated wind speeds are generally rather lower, due to weaker than observed pressure gradients close to the centres of the simulated storms. Simulations of the effect of climate change are performed. Under enhanced greenhouse conditions, simulated numbers of TCs do not change very much compared with those simulated for the current climate, nor do regions of occurrence. There is a 56% increase in the number of simulated storms with maximum winds greater than 30 m s^–1 (alternatively, a 26% increase in the number of storms with central pressures less than 970 hPa). In addition, there is an increase in the number of intense storms simulated south of 30°S. This increase in simulated maximum storm intensity is consistent with previous studies of the impact of climate change on tropical cyclone wind speeds.     

Global radiation climate changes in Israel

A detailed study of the 26-year series of global radiation K, measurements at Bet Dagan, the Israel Meteorological Service's pyranometer station in the central coastal plain of Israel, confirmed earlier findings of a significant reduction in insolation which were based on a small sample of this data set (Stanhill and Moreshet, 1992). Between 1956 and 1987 the annual reduction averaged 45.2 ± 4.3 MJ m^–2, equivalent to –0.63% yr^–1. Relatively the reduction was greater in midwinter (–0.91% yr^–1), than midsummer (–0.56% yr^–1), and under average (–0.63% yr^–1), than cloudless (–0.48% yr^–1) sky conditions. No changes were found in the degree of cloud cover observed at Bet Dagan.The annual decrease of K at Bet Dagan was highly correlated (r = –0.78) with the increase in the number of motor vehicles using the major roads passing within 1 km of the site: each additional vehicle passing was associated with a 21.5 J m^–2 decrease in K The causal nature of this correlation was confirmed by the difference of 18% found in daily values of K measured at Bet Dagan under traffic-free and extremely congested road conditions.The reduction in the K at Bet Dagan could not, however, be attributed exclusively to the increase in motor traffic in the immediate vicinity of the site, as no significant difference was found in values measured at a relatively traffic-free site 2 km downwind of the pyranometer station. The effect of aerosol pollutants originating in Tel Aviv - the major urban and industrial connurbation upwind of Bet Dagan - was confirmed by the changes recorded in the relative size of the direct and diffuse components of K measured at this site.The importance of pollution from Tel Aviv would also explain the absence of any significant changes in the annual values of K measured at Jerusalem, a smaller and less industrialized urban center 46 km downwind of Bet Dagan, or at Qidron, an uninhabited, isolated site on the NW coast of the Dead Sea, 25 km further downwind.Contribution from the Agricultural Research Organization, The Volcani Center, Bet Dagan, Israel. No. 3074-E, 1990 series.     

Mass transfer to bean leaves

A vapour of radio-lead (²¹²Pb) has been used to measure the Sherwood number, Sh, of model leaves at various angles of incidence,, to the airstream in a wind tunnel. The results for=0 are compared with Pohlhausen's formula and the results for 0, with Powell's experiments. The local values of Sh on the upwind and downwind sides of discs have been obtained. For leaves in the canopy, Sh was found to be about 25% greater than would be predicted by applying Pohlhausen's equation without correction for orientation.     

Three-dimensional simulation of the Denver 11 July 1988 microburst-producing storm

Summary Near Denver on July 11 1988, a moderate-reflectivity thunderstorm produced a microburst of unusual intensity during the test operation of the Terminal Doppler Weather Radar (TDWR) system. This microburst-producing storm, unlike those previously investigated, contained a complex updraft structure and an intense microburst downstream from the main precipitation shaft. Of special significance is that several successive airliners inadvertently encountered the microburst during final approach to Stapleton Airport. This paper focuses on the numerical investigation of this storm-via simulation with the Terminal Area Simulation System (TASS). Evolution and structure of the storm, including hazard indices based onF-factor, are presented and compared with observed data, including that measured by Doppler radar and aircraft flight data recorders.Model results show that multiple low-to moderate-reflectivity microbursts formed downshear of the main precipitation shaft. The most intense of the modeled microbursts contained a velocity differential exceeding 40m/s, strong downdrafts, and hazardous windshear withF-factors approaching 0.2. Outflow from this microburst was roughly symmetric at ground level during peak intensity, but became more asymmetric with time and altitude above the ground. The intense microburst, not unlike the real microburst that was encountered by the aircraft, originated from light precipitation that was swept downwind of the primary source. In contrast to the intense microburst, downdrafts and outflows associated with the simulated storm's main precipitation shaft were less vigorous by comparison, even though accompanied by higher magnitudes of radar reflectivity. Additional simulations with the axisymmetric version of TASS configured from the 3-D simulation, suggest that the intense microburst was driven by cooling primarily from sublimating snow.With 20 Figures