Atmospheric turbulent fluxes over snow

On March 26, 1971, eddy fluxes of momentum, sensible heat and water vapour were measured over Lake Mendota, Wisconsin, U.S.A., which was covered by an extensive snowfall. An evaporation rate of about 0.7mm day^–1 (2.2 mW cm^–2) was detected. Wind speeds were light and the atmosphere near the surface was highly stable. In these conditions, the average sensible heat transfer and Reynolds stress were -0.9 mW cm^–2 and 0.10 dyn cm^–2, respectively. Comparison with measured gradients of wind speed, temperature and humidity yield a drag coefficient of about 0.54 × 10^–3, and bulk transfer coefficients for sensible and latent heat of 0.41 × 10^–3 and 0.78 × 10^–3, respectively, applied to 10-m data. When corrected for the effect of atmospheric stability, these three coefficients become (in the same order) 1.2 × 10^–3, 0.9 × 10^–3 and 2.5 × 10^–3. The errors in these estimates are such that the drag coefficient is not significantly different from that corresponding to an aerodynamically smooth surface, while the heat coefficients are similar to those normally applied over liquid water surfaces.     

Some evaluations of drag and bulk transfer coefficients over water bodies of different sizes

Three recent experiments allow evaluation of the bulk transfer coefficients for momentum, water vapour and sensible heat over water bodies of different sizes. As part of a study of evaporation rates from a swamp, measurements of latent and sensible heat fluxes were made over Lake Wyangan in southern N.S.W., Australia. This lake is of several kilometers diameter. In a later experiment, Reynolds stress and sensible heat transfer were measured from a natural-gas platform standing in Bass Strait, south of mainland Australia. The most recent experiment involved the direct measurement of each of these turbulent fluxes from a fixed tower erected in Lake Michigan, U.S.A.Perhaps the most important of the results is the finding that drag coefficients measured over Bass Strait are not significantly different from those over Lake Michigan, despite the obvious differences in depth, fetch, and hence surface wave structure. At both locations, drag coefficients are found to increase slightly with increasing wind speed, while at low wind speeds they are not significantly different from those corresponding to aerodynamic smoothness.Near-neutral bulk transfer coefficients for sensible heat and for water vapour are found to be similar. An average value of about 1.4 × 10^–3 is obtained.It is emphasized that stability effects should be considered in any discussion of drag coefficients or bulk transfer coefficients. Large errors can result if near-neutrality is incorrectly assumed.     

Observations of the convective plume of a lake under cold-air advective conditions

Moderating effects of Lake Apopka, Florida on downwind surface temperatures were evaluated under cold-air advective conditions. Point temperature measurements north and south of the lake and data obtained from a thermal scanner flown at 1.6 km indicate that surface temperatures directly downwind may be higher than surrounding surface temperatures by as much as 5 °C under conditions of moderate winds (~4 m s^–1). No substantial temperature effects were observed with surface wind speed less than 1 m s^–1. Fluxes of sensible and latent heat from Lake Apopka were calculated from measurements of lake temperature, net radiation, relative humidity and air temperature above the lake. Bulk transfer coefficients and the Bowen ratio were calculated and found to be in agreement with reported data for non-advective conditions.IFAS Journal Series No. 1006.     

Convective heat transfer measurements of plants in a wind tunnel

Heat transfer was studied between intact leaves of various sizes and shapes in vivo under free and forced air conditions. Use of a wind tunnel and a microwave transmitter to heat the leaves facilitated measurements of convective, along with radiative and evaporative, heat losses from plant leaves. Knowledge of input energy, analysis of cooling curves, and established formulae, respectively, formed the basis of the steady-state, unsteady-state, and analytical methods for the determination of heat transfer coefficients.Typical values of steady-state free convection coefficients for Peperomia obtusifolia varied from 1.5 × 10^–4 to 1.9 × 10^–4 cal cm^–2 s^–1 C^–1 as the temperature difference was increased from 5.9 to 9.6°C, whereas the forced convection coefficient was found to be 4.2 × 10^–4 cal cm^–2 s^–1 C^–1 at 122 cm s^–1 wind velocity. For egg-plant, this value was about 9 × 10^–4 cal cm^–2 s^–1 C^–1 at 488 cm s^–1 wind velocity. Convection coefficients as determined under steady-state conditions are compared with those of the unsteady-state and with analytical values for a single leaf and leaves of three different plants. In general, experimental values were found to be higher than the analytical ones.     

Air mass modification over Lake Michigan

The modification of a relatively cold air mass over the warm water of Lake Michigan is studied by using a two-dimensional nonlinear mesoscale model. Considerable amounts of heat and water vapor are supplied from the water surface to the lower atmosphere by turbulent eddies. A convective mixed layer develops and grows toward the downwind region with stratocumulus clouds over the lake.The model simulates the warming and moistening of the mixed layer, the development of a boundary layer, the divergence and convergence of wind near the coastlines, and the turbulent fluxes.The model warming of the mixed layer across the lake was about 2.2 °K and the moistening of the mixed layer was about 0.8 g kg^–1, which are comparable to 2.7 °K and 0.8 g kg^–1 observed by Lenschow (1973). The convective boundary layer, which includes the cloud layer, subcloud layer, and superadiabatic layer near the water surface, is well simulated. The tilt of the inversion which coincides with the cloud top is also well reproduced. When a prescribed cooling rate is applied at the cloud top, stronger turbulence and a deeper cloud layer are generated. Without the cooling, the cloud is shallow and the shape of the cloud base is determined by surface conditions. The rise of the inversion is due to upward vertical motion, and deepening of the convective layer in the downwind region.     

An inverse method for determining wind stress from water-level fluctuations

The hydrodynamic equations governing the water-level response of a lake to wind stress are inverted to determine wind stress from water-level fluctuations. In order to obtain a unique solution, the wind-stress field is represented in terms of a finite number of spatially dependent basis functions with time-dependent coefficients. The discretized version of the inverse equation is solved by a least-squares procedure to obtain the coefficients, and thereby the stress. The method is tested for several ideal cases with Lake Erie topography. Real water-level data is then used to determine hourly values of vector wind stress over Lake Erie for the period 5 May–31 October, 1979. Results are compared with measurements of wind speed and direction from buoys deployed in the lake. Calculated stress direction agrees with observed wind direction for wind speeds > 7.5 m s⁻¹. Under neutral conditions, calculated drag coefficients increase with the wind speed from 1.53 × 10⁻³ for 7.5−10 m s⁻¹ winds to 2.04 × 10⁻³ for 15−17.5 m s⁻¹ winds. Drag coefficients are lower for stable conditions and higher for unstable conditions.     

Atmospheric Loss Processes of Dimethyl and Diethyl Carbonate

A combined study of the OH gas phase reaction and uptake on aqueous surfacesof two carbonates, dimethyl and diethyl carbonate has been carried out todetermine the atmospheric lifetimes of these compounds. Rate coefficients havebeen measured for gas phase reactions of OH radicals with dimethyl and diethylcarbonate. The experiments were carried out using pulsed laser photolysis– laser induced fluorescence over the temperature range 263–372K and the kinetic data were used to derive the following Arrhenius expressions(in units of cm³ molecule^–1 s^–1):for dimethyl carbonate, k₁ = (0.83±0.27)×10^–12 exp [–(247± 98)/T] and fordiethyl carbonate, k₂ = (0.46±0.15)×10^–12 exp [(503± 203)/T]. At 298 K, therate coefficients obtained (in units of 10^–12 cm³molecule^–1 s^–1) are: k₁ =(0.35± 0.04) and k₂ = (2.31± 0.29). The results arediscussed in terms of structure-activity relationships.The uptake coefficients of both carbonates on aqueous surfaces were measuredas a function of temperature and composition of the liquid phase, using thedroplet train technique coupled to a mass spectrometric detection. Dimethyland diethyl carbonate show very similar results. For both carbonates, themeasured uptake kinetics were found to be independent of the aqueous phasecomposition (pure water, NaOH solutions) but dependent on gas-liquid contacttime which characterises a surface saturation effect. The uptake coefficientvalues show a slight negative temperature dependence for both carbonates.These values vary from 1.4×10^–2 to0.6×10^–2 in the temperature range of 265–279 Kfor dimethyl carbonate, from 2.4×10^–2 to0.9×10^–2 in the temperature range of 270–279 Kfor diethyl carbonate. From the kinetic data, the following Henry's lawconstants were derived between 279 and 265 K: dimethyl carbonate,H₁ = 20–106 M atm^–1; and diethyl carbonate,H₂ = 30–98 M atm^–1. The reported data showthat the OH reaction is the major atmospheric loss process of these twocarbonates with lifetimes of 33 and 5 days, respectively, while the wetdeposition is a negligible process.     

The Compounding Effects of Tropical Deforestation and Greenhouse Warming on Climate

This study reports the first assessment of the compounding effects of land-use change and greenhouse gas warming effects on our understanding of projections of future climate. An AGCM simulation of the potential impacts of tropical deforestation and greenhouse warming on climate, employing a version of NCAR Community Climate Model (CCM1-Oz), is presented. The joint impacts of tropical deforestation and greenhouse warming are assessed by an experiment in which removal of tropical rainforests is imposed into a greenhouse-warmed climate. Results show that the joint climate changes over tropical rainforest regions comprise large reductions in surface evapotranspiration (by about –180 mm yr^–1) andprecipitation (by about –312 mm yr^–1) over the Amazon Basin, along with anincrease of surface temperature by +3.0 K. Over Southeast Asia, similar but weaker changes are found in this study. Precipitation is decreased by –172 mmyr^–1, together with the surface warming of 2.1 K. Over tropical Africa, changes in regional climate is much weaker and with some different features, such as the increase of precipitation by 25 mm yr^–1. Energy budgetanalyses demonstrates that the large increase of surface temperature in the joint experiment is not solely produced by the increase of CO₂concentration, but is a joint effect of the reduction of surface evaporation (due to deforestation) and the increase of downward atmospheric longwave radiation (due to the doubling of CO₂ concentration). Furthermore, impactsof tropical deforestation on the greenhouse-warmed climate are estimated by comparing a pair of tropical deforestation simulations. It is found that in CCM1-Oz, deforestation has very similar impacts on greenhouse-warmed regional climates as on current climates over tropical rainforest regions. The extra-tropical climatic response to tropical deforestation is identified in both sets of tropical deforestation experiments. Statistically significant responses are seen in the large-scale atmospheric circulation such as changes in the velocity potential and vertically integrated kinetic and potential energy fields. Wave propagation patterns are identified in the large-scale circulation anomalies, which provides a mechanism for interpreting the model responses in the extra-tropics. In addition, this study suggests that land-use change such as tropical deforestation may affect projections of future climate.     

A theory for the scalar roughness and the scalar transfer coefficients over snow and sea ice

Although the bulk aerodynamic transfer coefficients for sensible (C _H ) and latent (C _E ) heat over snow and sea ice surfaces are necessary for accurately modeling the surface energy budget, they have been measured rarely. This paper, therefore, presents a theoretical model that predicts neutral-stability values of C _H and C _E as functions of the wind speed and a surface roughness parameter. The crux of the model is establishing the interfacial sublayer profiles of the scalars, temperature and water vapor, over aerodynamically smooth and rough surfaces on the basis of a surface-renewal model in which turbulent eddies continually scour the surface, transferring scalar contaminants across the interface by molecular diffusion. Matching these interfacial sublayer profiles with the semi-logarithmic inertial sublayer profiles yields the roughness lengths for temperature and water vapor. When coupled with a model for the drag coefficient over snow and sea ice based on actual measurements, these roughness lengths lead to the transfer coefficients. C _E is always a few percent larger than CH. Both decrease monotonically with increasing wind speed for speeds above 1 m s^–1, and both increase at all wind speeds as the surface gets rougher. Both, nevertheless, are almost always between 1.0 × 10^–3 and 1.5 × 10^–3.     

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.     

One year measurements of aerosol optical properties over an urban coastal site: Effect on local direct radiative forcing

We present results of direct aerosol radiative forcing over a French Mediterranean coastal zone based on one year of continuous observations of aerosol optical properties during 2005–2006. Monthly-mean aerosol optical depth at 440 nm ranged between 0.1 and 0.34, with high Angstrom coefficient (α > 1.2). The single scattering albedo (at 525 nm) estimated at the surface ranged between 0.7 and 0.8, indicating significant absorption. The presence of aerosols over the Mediterranean zone during summer decreases the shortwave radiation reaching the surface by as much as 26 ± 3.9 W m^− 2, and increases the top of the atmosphere reflected radiation by as much as 5.2 ± 1.0 W m^− 2. The shortwave atmospheric absorption translates to an atmospheric heating of 2.5 to 4.6 K day^− 1. Concerted efforts are needed for investigating the possible impact of the increase in heating rate on the maintenance of heat-waves frequently occurring over this coastal region during summer time.     

Transfer processes in vegetation by electrochemical analog

Electrochemical modelling may provide fast order-of-magnitude estimates of energy and mass transfer in crops and partly fill the gap between rigid, expensive field experiments and oversimplified mathematical modelling. The technique comprises breaking up a flowing electrolyte in a manner analogous — but not necessarily completely similar — to the way the atmospheric flow is broken up by crops and studying the flow of ions in the electrolyte.The effects of variations in free-flow velocity, plant spacing, row spacing and orientation on velocity field, transfer coefficients and eddy diffusivities have been studied in an electrochemical model under conditions relevant to well-ventilated crops with relatively simple canopy structure. A one-dimensional analytical model is developed for the transfer at solid surfaces, with direct proportionality between transfer coefficient and eddy diffusivity. The proportionality constant has values in the model of (3.5±0.5)×10^–3cm^–1 at the ground and (5+-2)×10^–4cm^–1 (for a cylindrical probe) near canopy top. The latter is not too different from the values proposed for foliage surfaces in real canopies, viz., approximately 1×10^–3 cm^–1 (Philip, 1964) and 7×10^–4 cm^–1 (Uchijima, 1966). may therefore become a useful parameter in scaling transfer coefficients from systems with different molecular diffusion properties.A tentative extrapolation of measured data to corresponding values in air is given. The agreement between predicted values and presently available field data is encouraging although more precise field data are required for a final judgment on the validity of the model.This paper is published with the permission of the International Institute of Heat and Mass Transfer; an abbreviated version of the paper is given in the Proceedings, Intl. Seminar on Heat and Mass Transfer in the Environment of Vegetation, Dubrovnik, Aug. 26–30, 1974.     

Sensitivity studies on the calculation of the radiation balance of urban surfaces: II. Longwave radiation

The model developed in Verseghy and Munro (1989) is extended to the calculation of the longwave radiation incident on building surfaces. When compared with field measurements, the average magnitude of error associated with model predictions is found to be 10 W m^–2. The effects of six simplifying assumptions are investigated. The neglect of horizon obstructions is found to lead to errors of up to 60 W m^–2; the assumption of wall temperatures equal to air temperatures results in errors of up to 35 W m^–2. The neglect of absorption and emission by air between pairs of walls causes errors of the same magnitude as those associated with the predictions of the rigorous model itself. Of the three remaining simplifying assumptions tested (the assumption of isotrophic sky radiation, the use of published values of emissivities instead of measured values, and the blackbody surface assumption), none results om errors >5 W m^–2. As in the shortware case, the errors are site-specific, but nevertheless indicate the care with which the use of simplifying assumptions must be approached.     

青藏高原总体输送系数的特征

利用中日亚洲季风机制合作研究计划设在西藏的 4个自动气象站（AWS）获得的5a多（199年7月～1998年12月）时次密集、观测连续的近地层梯度资料,以最小二乘法确定出相应站点各季节的地表粗糙度,并应用廓线-通量法计算了4站逐日的总体输送系数,分析了其随时间的变化特征。结果表明:青藏高原动量输送系数的多年平均值为3.53×10-3～4.99×10-3,热量输送系数为4.67×10-3～6.73×10-3,并且两种输送系数都存在明显的日变化和季节变化,部分站点还存在较明显的年际变化。另外,还讨论了总体输送系数与近地层大气层结稳定度、地表粗糙度以及地面风速等因子的关系,初步建立了可用常规气象站地面观测资料计算青藏高原总体输送系数的拟合公式。     

The perturbed structure of the neutral atmospheric boundary layer over irregular terrain

The first-order (linear) response of the planetary boundary layer is calculated for flow over periodic terrain, for variations in both surface roughness and terrain elevation. Calculations are made for horizontal wavenumbers varying from 10^–4m^–1 to 3 × 10^–3m^–1. A simple second-order closure model of the turbulence is used, and Coriolis and buoyancy forces are neglected. As expected, flow over a periodic terrain produces corresponding periodic structure in all meteorological fields above the surface. The periodic structure consists of two components. The first is very nearly evanescent with height, showing little vertical structure. It corresponds to the motion that would be observed were the atmosphere inviscid. The second component, introduced by turbulent viscosity, exhibits considerable vertical structure, with vertical wavelengths the order of 100 m, and thus could be responsible for the layering often seen on acoustic sounder observations of the atmospheric boundary layer.Wave Propagation Laboratory.Environmental Science Group.     

Solar absorption at the surface and in the atmosphere as determined by different observational methods

Summary The concept of effective cloud cover, elaborated on the basis of an assumption that changes in the net radiation at the top of the atmosphere are mainly caused by changing cloudiness, has been used to deduce solar surface radiation from satellite data. It has been shown that the method permits a calculation of solar surface absorption distributions that agree well with the results obtained by other authors and that the existing disagreement can be to a great extent ascribed to the differences in the data sets and analysis periods. The method allows use of early satellite measurements to get longer time series of the surface radiation budget. In this study, it has been applied to the Nimbus-7 ERB WFOV data for 1979–1986.The net solar flux at the TOA (top of the atmosphere) can be partitioned into absorption at the surface and within the atmosphere. The geographical distributions of all the three quantities as well as the zonal averages of the surface absorption for January and July have been described. Special objectives of the present study are to estimate the interannual standard deviation for the 8-year period and to analyse the shortwave cloud-radiative forcing distributions at the surface and especially within the atmosphere.The standard deviation of the TOA and the surface solar absorption shows a temporal asymmetry, being much larger in January than in July. Noticeable is the disappearance of the wintertime strong variability over the central Pacific in July. As can be expected, the strong variability areas coincide with the strong variability areas of the cloud amount, showing the values up to 27 Wm^–2 at the surface.According to our estimate, the shortwave cloud forcing at the surface is everywhere stronger than that at the TOA, so that the cloud forcing of the atmosphere is negative. This means that in the belt of 58.5° N–58.5° S a cloudy atmosphere absorbs more solar energy than a cloud-free atmosphere. Our mean annual value of the atmospheric cloud forcing for this belt is –11 Wm^–2 which is somewhat stronger than that obtained by other investigators. It must be stressed that this value is within the uncertainty limits.Shortwave cloud forcing of the atmosphere is the strongest in the lower latitude areas of heavy cloudiness above the continents and negligible in the midlatitudes in winter. This gives evidence that the value of the shortwave cloud forcing of the atmosphere is modified by a combination of cloud absorption and cloud albedo.With 4 Figures     

Simulation of lake ice and its effect on the late-Pleistocene evaporation rate of Lake Lahontan

A model of lake ice was coupled with a model of lake temperature and evaporation to assess the possible effect of ice cover on the late-Pleistocene evaporation rate of Lake Lahontan. The simulations were done using a data set based on proxy temperature indicators and features of the simulated late-Pleistocene atmospheric circulation over western North America. When a data set based on a mean-annual air temperature of 3° C (7° C colder than present) and reduced solar radiation from jet-stream induced cloud cover was used as input to the model, ice cover lasting 4 months was simulated. Simulated evaporation rates (490–527 mm a^–1) were 60% lower than the present-day evaporation rate (1300 mm a^–1) of Pyramid Lake. With this reduced rate of evaporation, water inputs similar to the 1983 historical maxima that occurred in the Lahontan basin would have been sufficient to maintain the 13.5 ka BP high stand of Lake Lahontan.     

A Pulsed Laser Photolysis-Resonance Fluorescence Kinetic Study of the Atmospheric Cl Atom-Initiated Oxidation of Propene and a Series of 3-Halopropenes at Room Temperature

Absolute rate coefficient measurements have been carried out for the reactions of Cl atoms with propene and a series of 3-halopropenes, at room temperature (298 ± 2) K using a newly constructed laser photolysis-resonance fluorescence (PLP-RF) system. The rate coefficients obtained (in units of cm³ molecule^–1 s^–1) are: propene (1.40± 0.24) ×10^–10, 3-fluoropropene (4.92 ± 0.42) ×10^–11, 3-chloropropene (7.47 ± 1.50) × 10^–11, 3-bromopropene (1.23± 0.14) ×10^–10 and 3-iodopropene (1.29± 0.15) ×10^–10. In order to test this new system, the reactions of Cl atoms with acetone and isoprene have also been studied and compared with data previously reported. The rate coefficients determined at room temperature for these last two reactions are (2.93 ± 0.20) ×10^–12 cm³ molecule^–1 s^{– 1} and (3.64± 0.20)×10^–10 cm³ molecule^–1 s^–1, respectively. The measured values were independent of pressure over the range 20–200 Torr. The influence of the different halogen atoms substituents on the reactivity of these alkenes with Cl atoms as well as the atmospheric implications of these measurements are studied and discussed for the first time in this work and compared with the reactivity with NO₃ and OH radicals.     

Stable isotope (O and C) and pollen trends in eastern Lake Erie,evidence for a locally-induced climatic reversal of Younger Dryas age in the Great Lakes basin

A cool period from about 11000 to 10 500 BP (11 to 10.5 ka) is recognized in pollen records from the southern Great Lakes area by the return of Picea and Abies dominance and by the persistence of herbs. The area of cooling appears centred on the Upper Great Lakes. A high-resolution record (ca. 9 mm/y) from a borehole in eastern Lake Erie reveals, in the same time interval, this pollen anomaly, isotope evidence of meltwater presence (a — 3 per mil shift in ¹⁸O and a +1.1 per mil shift in ¹³C), increased sand, and reduced detrital calcite content, all suggesting concurrent cooling of Lake Erie. The onset of cooling is mainly attributed to the effect of enhanced meltwater inflow on the relatively large upstream Main Lake Algonquin during the first eastward discharge of glacial Lake Agassiz. Termination of the cooling coincides with drainage of Lake Algonquin, and is attributed to loss of its cooling effectiveness associated with a substantial reduction in its surface area. It is hypothesized that the cold extra inflow effectively prolonged the seasonal presence of lake ice and the period of spring overturn in Lake Algonquin. The deep mixing would have greatly increased the thermal conductive capacity of this extensive lake, causing suppression of summer surface lakewater temperatures and reduction of onshore growing-degree days. Alternatively, a rapid flow of meltwater, buoyed on sediment-charged (denser) lakewater, may have kept the lake surface cold in summer. Other factors such as wind-shifted pollen deposition and possible effects from the Younger Dryas North Atlantic cooling could have contributed to the Great Lakes climatic reversal, but further studies are needed to resolve their relative significance.Contribution to Climo Locarno — Past and Present Climate Dynamics; Conference September 1990, Swiss Academy of Sciences — National Climate ProgramGeological Survey of Canada Contribution 58 890     

Calibration and use of a sailplane variometer to measure vertical velocity fluctuations

The calibration of a sailplane variometer to measure vertical velocity fluctuations in the atmospheric boundary layer is described. Its usefulness is demonstrated with typical results from a boundary-layer development study. The atmospheric calibrations gave the ratio of standard deviations of vertical velocity fluctuations measured by a standard tower-mounted turbulence instrument to the values measured by variometer as 2.5 m s^–1 V^–1.