A Field Study of the Mean Pressure About a Windbreak

To provide additional field data for assessingwindbreak flow models, mean ground-level pressurehas been measured upstream and downstream from along porous fence (height H = 1.25 m, resistancecoefficient k _r = 2.4). Measurements were madeduring periods of near-neutral stability and near-normallyincident flow, with the fence standing on bare soil(roughness length, z ₀ 0.8 cm;H/z ₀ 160), or within a plant canopy. The mean pressure field,measured far from the ends of the fence, was foundto be quite insensitive to mean wind direction( , zero for perpendicular flow), for| | less than about 25°.In the absence of a canopy, during each measurementperiod the minimum pressure occurred at the closestsampling location to leeward of the windbreak, thepressure-gradient in most cases beingmaximally-adverse in the immediate lee, and decayingwith increasing downwind distance (x). On one day ofmeasurements, however, the pressure gradient over2 x/H 6 (H = windbreak height) resembled theleeward plateau identified by Wang and Taklein their numerical studies. Perhaps thisoccasional feature was only due to instrumenterror. Nevertheless a plateau of sorts wasindicated in similar measurements by Judd andPrendergast (with H = 1.92 m, z ₀ 1.2 cm;H/z ₀ 160, k _r 3). Therefore,existence of a leeward pressure plateau behind athin fence cannot be definitely ruled out.When the windbreak was placed in a canopy, minimumsurface pressure was displaced downwind. Thisagrees with the wind-tunnel study of Judd, Raupach and Finnigan,and is consistent with a simple simulation reported here.     

The Value Of The β Coefficient In The Relaxed Eddy Accumulation Method In Terms Of Fourth-Order Moments

The commonly measured value of in the relaxed eddy accumulationmethod of about 0.56is shown to arise from the non-Gaussiannature of turbulence. Fourth-orderGram–Charlier functions forthe two-dimensional probability distributionsof variation in the horizontal component of wind velocityand concentrations of water vapour, carbondioxide and methane with respect to thevertical component of wind velocity are used to examinethe value of .An analytical solution for ispresented in terms of fourth-order moments.Under mean conditions, this solution givesa value for of0.557. Variation of is shown to be controlledprimarily by the ratio of the mean ofc'w3 (where c'is relevant to the entity of interest andw' is vertical component of windvelocity) to the correlationcoefficient between the entity concentrationand vertical component of wind velocity.     

Characterization Of Atmospheric Turbulence By Dynamical Systems Techniques

We describe a dimensional analysis for different time series of the vertical component of wind velocity from half-hourly sonic anemometermeasurements. The goal is a characterization of the atmospheric turbulence fromthe point of view of Dynamical Systems Theory, based on the correlation dimension of the strange attractor. Our results suggest that the convective turbulence correlation dimension (values 6) is lowerthan the mechanical one (values 7–9).     

Synchronous Scan and Excitation-Emission Matrix Fluorescence Spectroscopy of Water-Soluble Organic Compounds in Atmospheric Aerosols

Three-dimensional excitation–emission matrix (EEM) fluorescence spectra of water-soluble organic compounds (WSOC) from aerosol samples were measured and compared with those reported in the literature for natural dissolved organic matter. The EEM profiles of the WSOC presented three characteristic excitation/emission (_Exc/_Em) peaks: 240/405 nm, 310/405 nm and 280/340 nm. The fluorescence intensities at _Exc/_Em240/405 nm and _Exc/_Em310/405 nm are located at wavelengths shorter than those reported for aquatic humic substances, indicating a smaller content of both aromatic structures and condensed unsaturated bond systems in the WSOC fraction. The EEM profiles of fractions obtained by the isolation procedure of the WSOC by the XAD resins showed that a fractionation has occurred and the XAD-8 eluate is highly representative of the total WSOC of collected aerosol. Synchronous scan spectra were more detailed than conventional fluorescence emission spectra, appearing more suitable for studying multicomponent samples such as the WSOC from atmospheric aerosols.     

THE EJECTION-SWEEP CHARACTER OF SCALAR FLUXES IN THE UNSTABLE SURFACE LAYER

In the atmospheric surface layer, it is widely accepted that ejection andsweep eddy motions, typically associated with coherent structures, areresponsible for much of the land-surface evaporation, sensible heat, andmomentum fluxes. The present study analyzes the ejection-sweep propertiesusing velocity and scalar fluctuation measurements over tall natural grassand bare soil surfaces. It is shown that momentum ejections and sweeps occurat equal frequencies (D eject D sweep 0.29) irrespective of surfaceroughness length or atmospheric stability conditions. Also, their magnitudesare comparable to values reported from open channel velocity measurements (Dsweep 0.33; D eject : 0.30). The scalar D eject is constant andsimilar in magnitude to the momentum D eject( 0.29) over both surfacesand for a wide range of atmospheric stability conditions, in contrast to thescalar D sweep. The scalar sweep frequency is shown to depend on the scalarskewness for the dynamic convective and free convective sublayers, but isidentical to D eject for the dynamic sublayer. The threshold scalar skewnessat which the D sweep dependence occurs is 0.25, in agreement with theaccepted temperature skewness value at near-neutral conditions. In contrastto a previous surface-layer experiment, this investigation demonstrates thatthe third-order cumulant expansion method (CEM) reproduces the measuredrelative flux contribution of ejections and sweeps (S0) for momentumand scalars at both sites. Furthermore, a linkage between S0 and thescalar variance budget is derived via the third-order CEM in analogy tomomentum. It is shown that S0 can be related to the flux divergenceterm and that such a relationship can be estimated from surface-layersimilarity theory, and the three sublayer model of Kader and Yaglom andproposed similarity functions.     

Analysis of Astex-Stratocumulus Observational Data Using a Mass-Flux Approach

A mass-flux approach is applied to observational data obtained in a convective boundary layer topped with stratocumulus clouds. The observational data were obtained from aircraft measurements during the Atlantic Stratocumulus Transition Experiment (ASTEX). A conditional sampling method is used to calculate average updraft and downdraft values. The vertical fluxes calculated with the mass-flux approach are found to be proportional to the real (measured) fluxes, with a proportionality factor being about 0.6. This value is predicted by theory for two variables having a joint Gaussian distribution function; proportionality factor = 2^-1 0.637. The horizontal fractional entrainment and detrainment rates calculated from the data ( 1–2 × 10^-2 m^-1) are an order of magnitude higher than the rates obtained by large eddy simulations for cumulus convection ( 2–3 × 10^-3 m^-1) and two orders of magnitude higher than those used in modelling cumulus convection with a mass-flux scheme in an operational weather forecast model ( 3 × 10^-4 m^-1). A numerical mass-flux model for the thermodynamics was developed and showed that results are in good agreement when compared with measured profiles of the liquid water content.     

Final results of the CONDORS convective diffusion experiment

TheConvectiveDiffusionObserved byRemoteSensors (CONDORS) field experiment conducted at the Boulder Atmospheric Observatory used innovative techniques to obtain three-dimensional mappings of plume concentration fields, /Q, of oil fog detected by lidar and chaff detected by Doppler radar. It included extensive meteorological measurements and, in 1983, tracer gases measured at a single sampling arc. Final results from ten hours of elevated and surface release data are summarized here. Many intercomparisons were made. Oil fog /Q measured 40m above the arc are mostly in good agreement withSF ₆ values, except in a few instances with large spacial inhomogeneities over short distances. After a correction scheme was applied to compensate for the effect of its settling speed, chaff dy/Q agreed well with those of oil except in two cases of oil fog hot spots. Mass or frequency distribution vs. azimuth or elevation angle comparisons were made for chaff, oil, and wind, with mostly good agreements. Spacial standard deviations, _y and _z, of chaff and oil agree overall and are consistent at short range with velocity standard deviations _vand _w 0.6w^* (the convective scale velocity), as measured atz>100m. Surface release _y is enhanced up to 60% at smallx, consistent with the Prairie Grass measurements and with larger _v and reduced wind speed measured near the surface. Decreased _y at small dimensionless average times is also noted. Finally, convectively scaled dy, C _y, were plotted versus dimensionlessx andz for oil, chaff, and corrected chaff for each 30–60 min period. Aggregated CONDORSC _y fields compare well with laboratory tank and LES numerical simulations; surface-released oil fog compares expecially well with the tank experiments. However, large deviations from the norm occurred in individual averaging periods; these deviations correlated strongly with anomalies in measured distributions.On assignment to the US Environmental Protection Agency, Atmospheric Research and Exposure Assessment Laboratory, RTP, NC.     

Transport-Dissipation Analytical Solutions to the E-∈Turbulence Model and their Role in Predictions of the Neutral ABL

E- turbulence model predictions of the neutralatmospheric boundary layer (NABL) are reinvestigated to determine thecause for turbulence overpredictions found in previous applications. Analytical solutions to the coupled E and equations for the case of steady balance between transport and dissipation terms, the dominant balance just below the NABL top, are derived. It is found that analytical turbulence profiles laminarizeat a finite height only for values of closure parameter ratio c ₂ /_e equal toor slightly greater than one, with laminarization as z for greater . The point = 2 is additionally foundthat where analytical turbulent length scale (l) profilesmade a transition from ones ofdecreasing ( < 2) to increasing ( > 2)values with height. Numerically predicted profiles near the NABL topare consistent with analytical findings. The height-increasingvalues of l predicted throughout the NABL with standard values ofclosure parameters thus appear a consequence of 2.5(> 2), implied by these values (c ₂ = 1.92, _{= 1.3}, _{e = 1}). Comparison of numericalpredictions with DNS data shows that turbulence overpredictions obtained with standard-valued parameters are rectifiedby resetting and _e to 1.1 and 1.6, respectively, giving, with c ₂ = 1.92, 1.3, and laminarization of the NABL's cappingtransport-dissipation region at a finite height.     

The Effect of Atmospheric Stability on the Surface-Layer Characteristics in a Low-Wind Area of Tropical West Africa

This study details the observed effects ofatmospheric stability on characteristics of thesurface layer in a low wind speed (U = 1.5 m s-1)regime of tropical West Africa. Theaerodynamic roughness length, z0, anddisplacement height, d, obtained from profilewind-speed data at our bush land site (height 2 m)have values of 0.24 ± 0.10 m and 1.54 ± 0.04 mrespectively. In the unstable range (-2.5 < Ri < -0.1; Riis gradient Richardson number), thestandard deviation in wind speed fluctuations, u, increased from 0.57 ± 0.19 m s-1 toa maximum of 0.7 ± 0.2 m s-1 in near-neutralconditions, and in the stable range, the parameterdecreased rapidly to 0.41 ± 0.15 m s-1 at Ri 0.2.In the same stability range, the horizontal winddispersion, , decreased withincreasing stability from 19 ± 8 deg. to 13 ± 5 deg.The surface-layer integral quantity, u/u*, when plottedas a function of stability, is in agreement with theempirical results. The ratio ofsensible heat flux (estimated) to the net radiationranged between 0.1 and 0.2 at nighttime,increasing to about 0.5 during the daytime, and showeda strong dependency on season.     

A study of a fair weather boundary layer in TOGA-COARE: Parameterization of surface fluxes in large scale and regional models for light wind conditions

The study focuses on a way to parameterize the effect of subgrid scale convective motions on surface fluxes in large scale and regional models for the case of light surface winds. As previously proposed, these subgrid effects are assumed to scale with the convection intensity through the relationship: where is the mean velocity of the wind, U₀ the velocity of the mean wind, w_* the free convection velocity, and an empirical coefficient to be determined. Both observations and numerical simulation are presently used to determine the free convection coefficient .Large eddy simulation of a fair weather convective boundary layer case observed during TOGA-COARE is performed. Comparisons between observations and the simulation of surface properties and vertical profiles in the planetary boundary layer are presented. The simulated vertical turbulent fluxes of heat, moisture and buoyancy range well within estimates from aircraft measurements.The most important result is that the true free convection coefficient , directly estimated from simulation, leads to a value of 0.65, smaller than the ones estimated from temporal and spatial variances. Using observations and simulation, estimates of from temporal and spatial variances are obtained with similar values 0.8. From both theoretical derivations and numerical computations, it is shown that estimates of the true from variances are possible but only after applying a correction factor equal to 0.8. If this correction is not used, is overestimated by about 25%. The time and space sampling problem is also addressed in using numerical simulations.     

On some new statistical characteristics of the wind “Koshava”

Summary Koshava is a gusty wind of moderate to strong intensity, blowing from the south-eastern direction, over the area of the Republic of Serbia. It is caused by the interaction between the synoptic circulation and the orography of the Carpathian and the Balkan mountains. The Koshava wind can damage buildings, factories and industrial plants or city infrastructure. Therefore it is important to estimate its gust and the gustiness factor on the basis of the measured data.This paper discusses a statistical analysis of wind data in the maximum influence area of the Koshava wind in the periods of maximum duration of Koshava. The focus of the paper is the examination of urban and suburban effects on Koshava wind and the correlation between the instantaneous maximum wind speed and the hourly mean wind speed. The best fitting with various empirical distributions is proposed.With 10 Figures     

Scavenging Efficiency of ‘Aerosol Carbon’ and Sulfate in Supercooled Clouds at Mt. Sonnblick (3106 m a.s.l., Austria)

Cloud water and interstitial aerosol samples collected at Mt. Sonnblick (SBO) were analyzed for sulfate and aerosol carbon to calculate in-cloud scavenging efficiencies. Scavenging efficiencies for sulfate (_SO) ranged from 0.52 to 0.99 with an average of 0.80. Aerosol carbon was scavenged less efficiently with an average value (_AC) of 0.45 and minimum and maximum values of 0.14 and 0.81, respectively. Both _SO and _AC showed a marked, but slightly different, dependence on the liquid water content (LWC) of the cloud. At low LWC, _SO increased with rising LWC until it reached a relatively constant value of 0.83 above an LWC of 0.3 g/m³. In the case of aerosol carbon, we obtained a more gradual increase of _AC up to an LWC of 0.5 g/m³. At higher LWCs, _{_} remained relatively constant at 0.60. As the differences between _SO and _A varied across the LWC range observed at SBO, we assume that part of the aerosol carbon was incorporated into the cloud droplets independently from sulfate. This hypothesis is supported by size classified aerosol measurements. The differences in the size distributions of sulfate and total carbon point to a partially external mixture. Thus, the different chemical nature and the differences in the size and mixing state of the aerosol particles are the most likely candidates for the differences in the scavenging behavior.     

Development of sea surface temperature,surface wind and divergence anomalies during a composite ENSO episode

Summary The interannual variations in sea surface temperature (SST) in the equatorial east Pacific, which are dominated by the El Niño phenomenon, are shown for the period 1870–1983. Since 1870 25 significant warm events have occurred. These events are classified as weak, moderate, strong and very strong, according to the normalized SST anomalies in the region 130° W–80° W, 0°–5° S.The spatial and temporal development of a composite El Niño/Southern Oscillation (ENSO) episode, based on 10 very strong or strong events, is presented in terms of SST, surface wind and divergence anomalies for the tropical Pacific (10° N–30° S). During its evolution the following phases are distinguished: Antecedent Conditions, Onset Phase, Peak Phase, Mature Phase and Dissipation Stage.Some aspects of ocean-atmosphere interaction associated with this evolution and, more specifically, the initiation of the composite event, are described. Seasonally varying feedback processes between SST, surface wind and convergence anomaly patterns in the western Pacific/Indonesian region suggest a possible mechanism for the initiation of typical ENSO events.

---

Zusammenfassung Die interannuären Variationen der Meeresoberflächentemperatur (SST) im äquatorialen Ostpazifik, die von dem El Niño-Phänomen dominiert werden, werden für die Periode 1870–1983 aufgezeigt. Seit 1870 traten 25 signifikante Ereignisse auf. Diese Ereignisse werden entsprechend den normierten SST-Anomalien in der Region 130° W–80° W, 0°–5° S als schwach, mittel, stark und sehr stark klassifiziert.Die räumliche und zeitliche Entwicklung einer Komposit-El Niño/Southern Oscillation (ENSO)-Episode, die auf 10 sehr starken bzw. starken Ereignissen basiert, wird anhand von SST-, Bodenwind- und Divergenzanomalien für den tropischen Pazifik (10° N–30° S) dargestellt. Während ihrer Entwicklung werden die folgenden Phasen unterschieden: Vorausgehende Bedingungen, Einsetzphase, Spitzenphase, Reifestadium und Auflösungsstadium.Einige Aspekte der Wechselbeziehungen Ozean—Atmosphäre werden im Zusammenhang mit der Entwicklung und insbesondere der Auslösung des Komposit-Ereignisses beschrieben. Jahreszeitlich variierende Rückkopplungsprozesse zwischen SST-, Bodenwind- und Konvergenzanomalien in der westpazifischen/indonesischen Region stellen einen möglichen Mechanismus für die Auslösung typischer ENSO-Ereignisse dar.

---

---

With 9 Figures     

An improved calibration for tethered balloon wind measurements

A previously published technique for using tethered spherical balloons as anemometers for measuring light low-level winds has been further developed. Earlier data on the relationship between the aerodynamic drag coefficient and the Reynolds number of spherical rubber balloons were combined with a large number of new data and re-analysed; and the errors in the relationship were estimated. The results allowed a more accurate calculation of wind speed from the deflection of a tethered balloon from the vertical. When combined with a new technique for calculating the effects of the tether, this enabled light to moderate low-level winds at fixed heights up to 600 m or more to be measured with simple, cheap, and readily mobile equipment; and a slight modification of the technique allowed measurement of winds in and above fog. Wind speeds measured by the ballon technique showed reasonably good agreement with measurements by an anemometer carried beneath the balloon.Glossary of Symbols a, b, c Coefficients in the relationship between lnC _d and lnR - A Quantity under square root in solution for lnV whena0 - C _d Wind drag coefficient for balloon - C _dc Value ofC _d given by calibration curve of Table I - D Dynamic wind pressure force on balloon - F Buoyant free lift of balloon with load - Re Reynold's number of balloon (sphere) - R = Re/10⁵ - r Radius of sphere - T Tension in tether - V Wind speed - ₈₃() =(lnC _dc -lnC _d ) when 83° , or 0 for other - Error in lnC _d - Elevation of tether where attached to balloon - Elevation of balloon from ground tether point - Molecular viscosity of air - Ratio of circumference to diameter of circle - Density of air     

Estimates of diabatic wind speed profiles from near-surface weather observations

In this paper we analyse diabatic wind profiles observed at the 213 m meteorological tower at Cabauw, the Netherlands. It is shown that the wind speed profiles agree with the well-known similarity functions of the atmospheric surface layer, when we substitute an effective roughness length. For very unstable conditions, the agreement is good up to at least 200 m or z/L–7(z is height, L is Obukhov length scale). For stable conditions, the agreement is good up to z/L1. For stronger stability, a semi-empirical extension is given of the log-linear profile, which gives acceptable estimates up to ~ 100 m. A scheme is used for the derivation of the Obukhov length scale from single wind speed, total cloud cover and air temperature. With the latter scheme and the similarity functions, wind speed profiles can be estimated from near-surface weather data only. The results for wind speed depend on height and stability. Up to 80 m, the rms difference with observations is on average 1.1 m s^–1. At 200 m, 0.8 m s^–1 for very unstable conditions increasing to 2.1 m s^–1 for very stable conditions. The proposed methods simulate the diurnal variation of the 80 m wind speed very well. Also the simulated frequency distribution of the 80 m wind speed agrees well with the observed one. It is concluded that the proposed methods are applicable up to at least 100 m in generally level terrain.     

Momentum and heat transfers in the surface layer over a frozen sea

Wind and temperature profiles from a 10-m mast were measured over the frozen Baltic Sea. The Monin-Obukhov similarity theory applies well to runs carefully selected according to stationarity criteria. This provides relatively low-scattered results with the roughness length z ₀ 0.04 cm, the drag coefficient C _D 1.5 × 10^-3 and the Stanton number C _H 1.00 × 10^-3 under near-neutral conditions. The roughness length, however, behaves in a peculiar way under intense stratification conditions. The reasons proposed for this could lead to an extension of the theory. The mechanisms for momentum and heat transfers are also examined, revealing that there are two different regimes, the smooth and the rough, for the wall friction but not for the heat transfer. Further, the scatter of the data for momentum and heat transfer under the aerodynamically rough regime can be explained by the distinction between the type of thermal stability conditions.     

Fossil-pollen evidence for abrupt climate changes during the past 18 000 years in eastern North America

A quantitative measure of the rate at which fossil-pollen abundances changed over the last 18 000 years at 18 sites spread across eastern North America distinguishes local from regionally synchronous changes. Abrupt regional changes occurred at most sites in late-glacial time (at 13700, 12 300, and 10000 radiocarbon yr BP) and during the last 1000 years. The record of abrupt late-glacial vegetation changes in eastern North America correlates well with abrupt global changes in ice-sheet volume, mountain snow-lines, North Atlantic deep-water production, atmospheric CO₂, and atmospheric dust, although the palynological signal varies from site to site. Changes in vegetation during most of the Holocene, although locally significant, were not regionally synchronous. The analysis reveals non-alpine evidence for Neoglacial/Little Ice Age climate change during the last 1000 years, which was the only time during the Holocene when climate change was of sufficient magnitude to cause a synchronous vegetational response throughout the subcontinent. During the two millennia preceding this widespread synchronous change, the rate of change at all sites was low and the average rate of change was the lowest of the Holocene.Contribution to Clima Locarno Past and Present Climate Dynamics; Conference September 1990, Swiss Academy of Sciences — National Climate Program     

Bistability of the thermohaline circulation identified through comprehensive 2-parameter sweeps of an efficient climate model

The effect of changes in zonal and meridional atmospheric moisture transports on Atlantic overturning is investigated. Zonal transports are considered in terms of net moisture export from the Atlantic sector. Meridional transports are related to the vigour of the global hydrological cycle. The equilibrium thermohaline circulation (THC) simulated with an efficient climate model is strongly dependent on two key parameters that control these transports: an anomaly in the specified Atlantic–Pacific moisture flux (F_a) and atmospheric moisture diffusivity (K_q). In a large ensemble of spinup experiments, the values of F_a and K_q are varied by small increments across wide ranges, to identify sharp transitions of equilibrium THC strength in a 2-parameter space (between Conveyor On and Off states). Final states from this ensemble of simulations are then used as the initial states for further such ensembles. Large differences in THC strength between ensembles, for identical combinations of F_a and K_q, reveal the co-existence of two stable THC states (Conveyor On and Off)—i.e. a bistable regime. In further sensitivity experiments, the model is forced with small, temporary freshwater perturbations to the mid-latitude North Atlantic, to establish the minimum perturbation necessary for irreversible THC collapse in this bistable regime. A threshold is identified in terms of the forcing duration required. The model THC, in a Conveyor On state, irreversibly collapses to a Conveyor Off state under additional freshwater forcing of just 0.1 Sv applied for around 100 years. The irreversible collapse is primarily due to a positive feedback associated with suppressed convection and reduced surface heat loss in the sinking region. Increased atmosphere-to-ocean freshwater flux, under a collapsed Conveyor, plays a secondary role.     

Analytical investigations of horizontal meso-scale momentum equations with Newtonian nudging

Summary This study presents an analytical investigation of the local behaviour of the solution to a mesoscale model with Newtonian nudging when observed winds are time varying. The analysis examines each Fourier component of the time series of observed winds. Unlike the case with a constant observed wind, the nudged wind vector does not asymptotically approach the observed wind. In response to sinusoidal oscillation of the observed wind, the nudged wind vector is always on a half circle connecting the vector ends of the observed and un-nudged modelled winds. When nudging parameter 0, the nudged wind vector approaches the un-nudged wind; when , the nudged wind vector approaches the observed wind. For commonly used values of nudging parameter , the modelled wind field always carries errors.A target nudging scheme is devised in this study in order to ensure the model result is identical to observed winds with sinusoidal oscillation. Investigation shows that such a target wind exists for a finite value of , and the magnitude of the target-nudging term is about the same as that of a normal nudging term if f, wheref is the Coriolis parameter and is the frequency of the wind oscillation.With 7 Figures     

Variangular wind spirals

The term variangular is introduced to emphasize a significant difference between the present and certain earlier solutions to the problem of organized airmotion within the planetary boundary layer. The latter belong to the family of equiangular wind spirals and have the characteristic that the angle () formed by the vectors of shearing stress and geostrophic departure is invariant with height; it is shown that in this spiral-family, parabolic height-dependency of the effective (eddy) diffusivity (K) alone is permitted, including the asymptotic case of constant K; the famous Ekman spiral as well as the Rossby spiral are two prominent members of the family of equiangular wind spirals. The new variangular theory, as the name implies, permits variation of with height (z) and produces more versatile profiles of wind and stress due to less restraint in K (z). As an example of comparison with observed data, monthly mean wind profiles obtained at Plateau Station, Antarctica, are selected since they exhibit a noteworthy degree of variangularity, in relatively satisfactory agreement with properties of the new theoretical model for wind spirals.National Research Council Visiting Scientist Research Associate, Regional Environments Division, Earth Sciences Laboratory.