Trace gas exchange at the air/water interface: Measurements of mass accommodation coefficients

A liquid jet of 90 m diameter and variable length has been utilized to determine absorption rates and, hence, mass accommodation coefficients , of atmospheric trace gases. The compounds investigated are HCl (0.01), HNO₃ (0.01), N₂O₅ (0.005), peroxyacetyl nitrate (>0.001), and HONO (0.005). It is concluded that the absorption of these trace gases by liquid atmospheric water is not significantly retarded by interfacial mass transport. The strengths and limitations of the liquid jet technique for measuring mass accommodation coefficients are explored.     

Re-analysis of the nuclear winter phenomenon

Summary An analysis of the report of the (U.S.) National Academy of Sciences (NAS) on atmospheric effects of a nuclear exchange leads to conclusions that differ from those of the NAS and of the earlier TTAPS and AMBIO studies. Any cooling of the earth's surface is likely to beshort-lived because of rapid removal of the smoke clouds originating from nuclear burst-initiated fires, andminor because of appreciable green-house effects due to several distinct physical causes. (One of these, neglected in prior analyses, is the infrared absorption from cirrus clouds produced directly by the nuclear bursts.) Taken together, these effects may even induce slight surface warming (nuclear summer) instead of cooling (nuclear winter). The consequences to atmospheric ozone are similarly ambiguous; depending on the detailed nuclear scenario, the net ozone content may increase-rather than decrease as argued by TTAPS. Experiments could settle some uncertainties.With 2 Figures     

Numerical modeling of the nocturnal urban boundary layer

A numerical case study with a second-order turbulence closure model is proposed to study the role of urban canopy layer (UCL) for the formation of the nocturnal urban boundary layer (UBL). The turbulent diffusion coefficient was determined from an algebraic stress model. The concept of urban building surface area density is proposed to represent the UCL. Calculated results were also compared with field observation data. The height of the elevated inversion above an urban center was simulated and found to be approximately twice the average building height. The turbulent kinetic energy k, energy dissipation rate , and turbulence intensities u ² and w ² increase rapidly at the upwind edge of the urban area. The Reynolds stress uw displayed a nearly uniform profile inside the UBL, and the vertical sensible heat flux w had a negative value at the inversion base height. This indicates that the downward transport of sensible heat from the inversion base may play an important role in the formation of the nocturnal UBL.     

A random-walk model for dispersion of heavy particles in turbulent air flow

A random-walk model is presented for calculating the dispersion of heavy particles in a turbulent air flow when only air turbulence statistics and the drag characteristics of the particle are known. Algebraic expressions for the modification of air velocity variance ² and Lagrangian autocorrelation tune-scale T _L,due to particle inertia effects, are derived. These expressions introduce only a very small computational overhead on the random-walk models for inertia-less particles of Wilson et al. (1983). Measurements of T _Land by Snyder and Lumley (1971) for four different particles are used to determine constants in the heavy-particle model. It is shown that the agreement between the model, for a single set of constants, and the dispersion measurements is good for the 47 m hollow glass, 87 m glass, and 47 m copper particles. The predictions for the 87 m corn pollen particles show less satisfactory agreement by underestimating dispersion measurements by 15% after 0.4s. Finally, some aspects of the model's application to spray dispersion in and above a crop canopy are considered.     

Analytical solutions for the Ekman layer

The PBL equation that governs the transition from the constant-stress surface layer to the geostrophic wind in a neutrally stratified atmosphere for which the eddy viscosityK(z) is assumed to vary smoothly from the surface-layer value U _*z (0.4,U _*=friction velocity,z=elevation) to the geostrophic asymptoteK _GU _*d forzd is solved through an expansion in fd/U _*1 (f=Coriolis parameter). The resulting solution is separated into Ekman's constant-K solution an inner component that reduces to the classical logarithmic form forzd and isO() relative to the Ekman component forzd. The approximationKU _*d is supported by the solution of Nee and Kovasznay's phenomenological transport equation forK(z), which yieldsK–U _*d exp(–z/d), where is an empirical constant for which observation implies, 1. The parametersA andB in Kazanskii and Monin's similarity relation forG/U _* (G=geostrophic velocity) are determined as functions of . The predicted values ofG/U _* and the turning angle are in agreement with the observed values for the Leipzig wind profile. The predicted value ofB based on the assumption of asymptotically constantK is 4.5, while that based on the Nee-Kovasznay model is 5.1; these compare with the observed value of 4.7 for the Leipzig profile. A thermal wind correction, an asymptotic solution for arbitraryK(z) and 1, and an exact (unrestricted ) solution forK(z)=U _*d[1–exp(–z/d)] are developed in appendices.     

A study of the seasonal migration of ITCZ and the quasi-periodic oscillations in a simple monsoon system using an energy balance model

Summary A zonally averaged global energy balance model with feedback mechanisms was constructed to simulate (i) the poleward limits of ITCZ over the continent and over the ocean and (ii) a simple monsoon system as a result of differential heating between the continent and the ocean. Three numerical experiments were performed with lower boundary as (1) global continent, (2) global ocean and (3) continent-ocean, with freezing latitudes near the poles. Over the continent, midlatitude deserts were found and the ITCZ migrates 25° north and south with seasons. Over a global swamp ocean results do not show migration of ITCZ with time but once the ocean currents are introduced the ITCZ migrates 5° north and south with seasons. It was found that the seasonal migration of ITCZ strongly depends on the meridional distribution of the surface temperature. It was also found that continent influences the location of the oceanic ITCZ. In the tropics northward progression of quasi-periodic oscillations called events are found during the pre- and post-monsoon periods with a period of 8 to 15 days. This result is consistent with the observed quasi-periodic oscillations in the tropical region. Northward propagation of the surface temperature perturbation appears to cause changes in the sensible heat flux which in turn causes perturbations in vertical velocity and latent heat flux fields.List of Symbols vertical average - ₀ zonal average - vertical mean of the zonal average - _0s zonal average at the surface - _0a zonal average at 500 mb level - latitude We now define the various symbols used in the model rate of atmospheric heating due to convective cloud formation (K/sec) - dp/dt (N/m²/sec) - density - potential temperature (K) - rate of rotation of the earth (rad/sec) - empirical constant - humidity mixing ratio - ^* saturated humidity mixing ratio - opacity of the atmosphere - ₁,₂ factors for downward and upward effective black body long wave radiation from the atmosphere - Stefan-Boltzmann constant - emissivity of the surface - _D subsurface temperature (K) - a specific volume - _0xs ,_0ys eastward and northward components of surface frictional stress - ^* vertical velocity at the top of the boundary layer (N/m²/sec) - P Thickness of the boundary layer (mb) - nondimensional function of pressure - P pressure - P _a pressure of the model atmosphere (N/m²) - P _s pressure at the surface (N/m²) - t time (sec) - U eastward wind speed (m/sec) - V northward wind speed (m/sec) - surface water availability - T absolute temperature (K) - heat addition due to water phase changes - g acceleration due to gravity (m²/sec) - a radius of the earth (m) - R gas constant for dry air (J/Kg/K) - C _p specific heat of air at constant pressure (J/Kg/K) - k R/C _p - L latent heat of condensation (J/Kg) - f coriolis parameter (rad/sec) - H _s H _0s ⁽¹⁾ +H _0s ⁽²⁾ +H _0s ⁽³⁾ +H _0s ⁽⁴⁾ +H _0s ⁽⁵⁾ (J/m²/Sec)=sum of the rates of vertical heat fluxes per unit surface area, directed toward the surface - H _a H _0a ⁽¹⁾ +H _0a ⁽²⁾ +H _0a ⁽³⁾ +H _0a ⁽⁴⁾ (J/m²/Sec)=sum of the rates of heat additions to the atmospheric column per unit horizontal area by all processes - H _0s ⁽¹⁾ ,H _0a ⁽¹⁾ heat flux due to short wave radiation - H _0s ⁽²⁾ ,H _0a ⁽²⁾ heat flux due to long wave radiation - H _0s ⁽³⁾ ,H _0a ⁽³⁾ heat flux due to small scale convection - H _0s ⁽⁴⁾ heat flux due to evaporation - H _0a ⁽⁴⁾ heat flux due to condensation - H _0s ⁽⁵⁾ heat flux due to subsurface conduction and convection - e ^* saturation vapor pressure - R solar constant (W/m²) - r _a albedo of the atmosphere - r _s albedo of the surface - b ₂ empirical constant (J/m²/sec) - c ₂ empirical constant (J/m²/sec) - e ₂ nondimensional empirical constant - f ₂ empirical constant (J/m₂/sec) - factor proportional to the conductive capacity of the surface medium - a _s constant used in Sellers model - b _s positive constant of proportionality used in the Sellers model (kg m²/J/sec²) - K _HT coefficient for eddy diffusivity of heat (m²/sec) - K _HE exchange coefficient for water vapor (m²/sec) - h depth of the water column (m) - z height (m) - V _0ws meridional component of surface current (m/sec) - n cloud amount - G _0,n long wave radiation form the atmosphere for cloud amount n (W/m²) - B ₀ long wave radiation from the surface (W/m²) - S _0,n short wave radiation from the atmosphere for cloud amount n (W/m²) - A _n albedo factor for a cloud amount n - R _f1 large scale rainfall (mm/day) - R _f2 small scale rainfall (mm/day) With 22 Figures     

Diurnal variation of horizontal wind direction fluctuations in complex terrain at Geysers,Cal.

Horizontal diffusion in the surface layer is dependent on the standard deviation of wind direction fluctuations . Diurnal variation of this parameter in complex terrain was studied for the July 1979 Geysers, Cal., experiment using data from a network of 11 short meteorological towers in the 25 km² Anderson Creek watershed Valley side slopes are roughly 20 ° and maximum terrain difference is about 1 km.Values of for wind directions sampled for one hour at a height of 10 m are about 35 ° during the daytime. They slowly decrease to about 20 ° by 8 to 10 p.m. as stability increases but wind speeds are still relatively high. After 10 p.m. the drainage flow sets in at most stations, with speeds of 1 to 2 m s^-1, and average increases to about 30° during the period 11 p.m. to 6 a.m. In general, highest values of at night are associated with lowest values of wind speed and greatest static stability. This enhancement of by the terrain suggests that horizontal diffusion at night always conforms to that expected during nearly neutral stabilities. That is, Pasquill class D diffusion applies to the horizontal component all night in complex terrain.     

The stably stratified internal boundary layer for steady and diurnally varying offshore flow

A two-dimensional numerical mesoscale model is used to investigate the internal structure and growth of the stably stratified internal boundary layer (IBL) beneath warm, continental air flowing over a cooler sea. Two situations are studied — steady-state and diurnally varying offshore flow. In the steady-state case, vertical profiles of mean quantities and eddy diffusion coefficients (K) within the IBL show small, but significant, changes with increasing distance from the coast. The top of the IBL is well defined, with large vertical gradients within the layer and a maximum in the coast-normal wind component near the top. Well away from the coast, turbulence, identified by non-zero K, decreases to insignificant levels near the top of the IBL; the IBL itself is characterised by a critical value of the layer-flux Richardson number equal to 0.18. The overall behaviour of the mean profiles is similar to that found in the horizontally homogeneous stable boundary layer over land.A simple physical model is used to relate the depth of the layer h to several relevant physical parameters viz., x, the distance from the coast and U, the large-scale wind (both normal to the coastline) and g/, being the temperature difference between continental mixed-layer air and sea surface, is the mean potential temperature and g is the acceleration due to gravity. Excellent agreement with the numerical results is found, with h = 0.014x ^1/2 U (g/)^–1/2.In the diurnally varying case, the mean profiles within the IBL show only small differences from the steady-state case, although diurnal variations, particularly in the wind maximum, are evident within a few hundred kilometres of the coast. A mesoscale circulation normal to the coast, and superimposed upon the mean offshore flow, develops seawards of the coastline with maximum vertical velocities about sunset, of depth about 2 km and horizontal scale 500 km. The circulation is related to the advection, and subsequent decay, of daytime convective turbulence over the sea.     

The carbon isotopic ratio of atmospheric carbon dioxide at Tsukuba,Japan

To find out the secular and seasonal trends of the ¹³C value and CO₂ concentration in the surface air and the determination of the ¹³C in the atmospheric CO₂ collected at Tsukuba Science City was carried out during the period from July 1981 to October 1983. The monthly average of the ¹³C value of CO₂ in the surface air collected at 1400 LMT ranged from -7.52 to \s-8.45 with an average of -7.96±0.25 and the CO₂ concentration in the air varied from 334.5 l 1^-1 to 359 l 1^-1 with an average of 347.2±6.3 l 1^-1. The ¹³C value is high in summer and low in winter and is negatively correlated with the CO₂ concentration. In general, the relationship between the ¹³C and the CO₂ concentration is explainable by a simple mixing model of two different constant carbon isotopic species but the relationship does not always follow the model. The correlation between the ¹³C value and the CO₂ concentration is low during the plant growth season and high at other times. The observed negative deviation of the ¹³C value from the simple mixing model in the plant growth season is partly due to the isotopic fractionation process which takes place in the land biota.     

Airflow above changes in surface heat flux,temperature and roughness; an extension to include the stable case

A numerical model of airflow above changes in surface roughness and thermal conditions is extended to include cases with stable thermal stratification within the internal boundary-layer. The model uses a mixing-length approach with empirical forms for M and H.Results are presented for some basic cases and an attempt is then made to compare results given by the model with the experimental results of Rider, Philip and Bradley. Tolerable agreement is achieved. The importance of roughness change and thermal stability effects in the diffusion of heat and moisture near a leading edge is emphasised.Notation A Refers to Taylor (1970) - B Businger-Dyer constant (= 16.0) in forms for _M and _H - C Constant in form for in stable case - c _p Specific heat at constant pressure - E Scaled absolute humidity - g Acceleration due to gravity - H Upward vertical heat flux - H ₀, H ₁ Surface heat fluxes for x <0, x0 - H _E Upward latent heat flux - k Von Kármán's constant (= 0.4) - K _H K _W Eddy transfer coefficients for heat and water vapour - L Monin-Obukhov length - L _H Latent heat of evaporation for water - m Ratio of roughness lengths ( = z ₁/z ₀) - RPB Refers to Rider et al. (1964) - RL* Non-dimensional parameter (see Equations (9), (20a), (22a), (24a)) - R* Net radiation less ground heat flux (see Equations (15), (16)) - T Scaled temperature - T ₁ Downstream scaled surface temperature - u ₀ u ₁(x) Surface friction velocities for x <0, x0 - U, W Horizontal and vertical mean velocities - x, z Horizontal and vertical co-ordinates - Z _i Local roughness length - z ₀, z _i Roughness lengths for x < 0, x 0 - Temperature - ₀, ₁ Surface temperatures for x<0, x0 - _E Non-dimensional absolute humidity gradient - _H Non-dimensional temperature gradient of heat flux - _M Non-dimensional wind shear - = _M = _H = _E an assumption used in stable conditions - Air density - Absolute humidity - _w Density of water - Kinematic shear stress - Logarithmic height scale (= ln(z+z ₁)/z ₁)     

Effect of land management on ecosystem carbon fluxes at a subalpine grassland site in the Swiss Alps

Summary The influence of agricultural management on the CO₂ budget of a typical subalpine grassland was investigated at the Swiss CARBOMONT site at Rigi-Seebodenalp (1025m a.s.l.) in Central Switzerland. Eddy covariance flux measurements obtained during the first growing season from the mid of spring until the first snow fall (17 Mai to 25 September 2002) are reported. With respect to the 10-year average 1992–2001, we found that this growing season had started 10 days earlier than normal, but was close to average temperature with above-normal precipitation (100–255% depending on month). Using a footprint model we found that a simple approach using wind direction sectors was adequate to classify our CO₂ fluxes as being controlled by either meadow or pasture. Two significantly different light response curves could be determined: one for periods with external interventions (grass cutting, cattle grazing) and the other for periods without external interventions. Other than this, meadow and pasture were similar, with a net carbon gain of –128±17g Cm^–2 on the undisturbed meadow, and a net carbon loss of 79±17g Cm^–2 on the managed meadow, and 270±24g Cm^–2 on the pasture during 131 days of the growing season, respectively. The grass cut in June reduced the gross CO₂ uptake of the meadow by 50±2% until regrowth of the vegetation. Cattle grazing reduced gross uptake over the whole vegetation period (37±2%), but left respiration at a similar level as observed in the meadow.     

On the similarity in atmospheric fluctuations of carbon dioxide,water vapor and temperature over vegetated fields

This paper describes the similarity between atmospheric fluctuations of carbon dioxide, water vapor and temperature using data which cover a wide range of instability (0.02 < < 10). The is the Monin-Obukhov stability parameter including the humidity effect.The spectral analysis shows that the coherency between fluctuations of carbon dioxide and water vapor or temperature is very close to unity, and the phase difference is basically out of phase for whole frequency ranges analyzed. The stability dependence of the normalized standard deviation of carbon dioxide is very similar to those of water vapor and temperature. The normalized standard deviation is about 2.5 under near neutral conditions, and it decreases with increasing instability following the -1/3; power law as (-)^-1/3. The skewness factors of carbon dioxide, water vapor and temperature show a systematic departure with increasing instabilities for 0.02 < s- < 1, and level off at high instabilities for 1 < -\s < 10. The stability dependence of the flatness factors is not so clear as that noted in the standrard deviation and skewness factors. Dissipation rates of carbon dioxide, water vapor and temperature variance are well related to the spectral peak wavelength. This seems to be real since the local production and local dissipation rates are the main terms, almost balancing one another in the variance budget equations for scalar entities.     

On the representation of frequency spectra in meteorology

The common representation of frequency spectra in meteorology and climatology is discussed. It is pointed out that this representation is misleading since spectral peaks and spectral gaps are obtained even when the spectrum density is monotonously decreasing in the whole frequency range. A plea is made for using the spectrum distribution function, F() which gives an unambiguous picture of the distribution of variance with frequency.     

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.     

Effects of water vapour on the structure parameter of the refractive index for near-infrared radiation

The refractive index of air (n) mainly depends on temperature andwater vapour content. For near-infrared radiation, temperature isthe main determining factor. To determinethe structure parameter of temperature ( ) fromthe structure parameter of the refractive index( ), the influence ofwater vapour content on n needs to be taken into account as acorrection. Three levels of approximation are presented. The firstinvolves the standard deviations of T and q (_T and_q) as well as the correlation coefficient betweenT and q (R_Tq). The second approximation involvesR_Tq and the Bowen ratio (), and the last uses only the Bowen ratio.The latter is the classical Bowen ratio correction.Evaluation of the validity of the assumptions used in the derivationreveals that a large error may be introduced (for situations withR_Tq < 0.8) by replacing _T/_qwith the Bowen ratio.Consequently, estimating from using the Bowen ratiocan contain errors of 5–40% for || < 1. Adding R_Tqas extra information reduces the errors only marginally. For || > 1,the correction is small, and all three approximations give errors of lessthan 1% in . When is used to compute the sensible heat flux,the influence of the quality of the correction forwater vapour fluctuations on the measured energy balance is small: forsmall ||, the correction is large, but the absolute value ofthe heat flux is small, whereas for large ||, the correctionis insignificant.     

Measurements of solar radiation and estimation of optical depth in the high Andes of Peru

Summary During an expedition to the high Andes of Southern Peru in June–July 1977, measurements of direct solar radiation in four spectral bands (0.270–0.530–0.630–0.695–2.900 ) were conducted at six sites in elevations ranging from sea level to 5645 m. These measurements were evaluated in Langley plots to determine total optical depths () and irradiances at the top of the atmosphere. In addition, water vapor optical depths (_wv) were calculated from the mean radiosounding over Lima during the expedition, and Rayleigh (_ray) and ozone (_oz) optical depths were obtained from published tabulations. Subtracting _ray, _oz, and _wv from yielded estimates of aerosol optical depth _aer. The components _ray and _oz decrease from the shorter towards the longer wavelength bands and from the lower towards the higher elevation sites; _aer also decreases towards the higher elevations. Particularly pronounced is the decrease of _aer and from the lowlands of the Pacific coast to the highlands of the interior, reflecting the effect of a persistent lower-tropospheric inversion and the contrast from the marine boundary layer to the clear atmosphere of the high Andes.With 4 Figures     

Free convection scaling over the ocean derived from buoy measurements during gate

The structure of atmospheric turbulence in the surface layer over the open ocean is examined under conditions of local free convection. The raw data consist of profile and fluctuation measurements of wind and temperature as obtained from a meteorological buoy. For near neutral conditions and for waves running approximately along the wind direction, wave-induced wind fluctuations can be described by a simplified linear theory based on Miles (1957). In this case, the spectrum of wind velocity is given as the sum of two parts; for the turbulent part, the parameterization as obtained by Kaimal et al. (1972) applies, while the wave-induced part is parameterized using a simplification of Miles' linear theory. For cases of local free convection, the measurements of the vertical component of the wind velocity are well described by similarity theory; as expected, _w /(-uw)^1/2 is proportional to (- z/L)^1/3. In order to scale the longitudinal wind velocity component, it seems to be reasonable to extend the list of relevant parameters by the height of the mixed layer z _i. We obtain _u /(- uw)^1/2 (z/z _i)^1/3(- z/L)^1/3 with only a poor correlation coefficient of r = 0.6. Overall, the results of local free convection scaling obtained from direct measurements show good agreement with those obtained from profile measurements. A comparison between direct and indirect determination of turbulent fluxes of momentum shows an unexplained difference of about 20%. This discrepancy is mainly due to a gap in the uw-cospectrum at the swell frequency.     

A review of flux-profile relationships

Flux-profile relationships in the constant flux layer are reviewed. The preferred relationships are found to be those of Dyer and Hicks (1970), namely, _H = _W =(1–16(z/L))^–1/2, _M =(1–16(z/L))^–1/4 for the unstable region, and _H = _W = _M = 1+5(z/L) for the stable region.The carefully determined results of Businger et al. (1971) remain a difficulty which calls for considerable clarification.     

Near-field dispersion from instantaneous sources in the surface layer

This paper considers the near-field dispersion of an ensemble of tracer particles released instantaneously from an elevated source into an adiabatic surface layer. By modelling the Lagrangian vertical velocity as a Markov process which obeys the Langevin equation, we show analytically that the mean vertical drift velocity w(t) is w()=bu _*(1–e ^–(1+)), where is time since release (nondimensionalized with the Lagrangian time scale at the source), b Batchelor's constant, and u _*, the friction velocity. Hence, the mean height and mean depth of the ensemble are calculated. Although the derivation is formally valid only when 1, the predictions for w, mean height and mean depth are consistent in the downstream limit ( 1) with surface-layer Lagrangian similarity theory and with the diffusion equation. By comparing the analytical predictions with numerical, randomflight solutions of the Langevin equation, the analytical predictions are shown to be good approximations at all times, both near-field and far-field.     

Cold surges along the front range of the Rocky mountains: Development of a classification scheme

Summary A classification scheme for cold surges occurring along the Front Range of the Rocky Mountains has been developed using NMC surface and upper air data. A separation is seen in the data between cases with parent anticyclones that moved from Canada south into the United States (Northern type) and those which moved from the Eastern Pacific Ocean east into the central US (Western type). The Northern surges are further broken down into three types, A, B and B, when the degree of southward motion of the anticyclone and the leading edge of the cold air is considered. An anticyclone movement south of 40°N distinguishes type A, while anticyclone tracks remaining north of 40°N characterize types B and B. The leading edge of the surge (the surge line) associated with the type A anticyclone travels south of 25N, while the type B anticyclones produce two distinct surge line movements, one remaining north of 25°N and the other moving south of 25°N. It is not possible to break down the Western surges into distinct types. As a rule, the anticyclones associated with Western surges remain between 50°N and 30°N, with the surge line moving south to 18°N on average.Though the classification scheme is developed in a relatively subjective manner, the statistical robustness of the scheme is verified through application of the nonparametric Mann-Whitney test. Results of the test clearly provide support for the division of the surges that was obtained.With 16 Figures