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     

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     

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.     

A study of the sensitivity of bare soil evaporation schemes to soil surface wetness,using the coupled soil moisture and surface temperature prediction model,BARESOIL

Summary The performance of evaporation schemes with and approach and their combination within resistance representation of evaporation from bare soil surface is discussed. For this purpose nine schemes, based on different functions of or , on the ratio of the volumetric soil moisture content and its saturated value are used.The quality of the chosen schemes has been evaluated using the results of time integration by the coupled soil moisture and surface temperature prediction model, BARESOIL, using in situ data. A sensitivity analysis was made using two sets of data derived from the volumetric soil moisture content of the top soil layer. One with values below the wilting point (0.17 m³m^–3) and the second with values above 0.20m³m^–3. Data sets were obtained at the experimental site Rimski anevi, Yugoslavia, from the bare surface of a chernozem soil.With 4 Figures     

Modification Of The Standard ∈-Equation For The Stable Abl Through Enforced Consistency With Monin–Obukhov Similarity Theory

A condition is derived for consistency of the standard-equation with Monin–Obukhov (MO) similarity theory of thestably-stratified surface layer. The condition is derivedby extending the procedure used to derive the analogous condition forneutral theory to stable stratification. It is shown that consistencywith MO theory requires a function of flux Richardson number, Ri_f, to be absorbed into either of two closure parameters, c ₁ or c ₂.Inconsistency, on the other hand, results if constant values of these are maintained for all Ri_f, as is done in standardapplication of the equation, and the large overpredictions ofturbulence found in such application to the one-dimensionalstable atmospheric boundary layer (1D-SBL) are traced to thisinconsistency. Guided by this, we formulate a MO-consistent-equation by absorbing the aforementioned function intoc ₁, and combine this with a Level-2.5 second-orderclosure model for vertical eddy viscosity and diffusivities.Numerical predictions of the 1D-SBL by the modified model converge to a quasi-steady state, rectifying the predictive failure of the standard -equation for the case.Quasi-steady predictions of non-dimensional variables agree stronglywith Nieuwstadt's theory. Qualitative accuracy of predictionsis inferred from comparisons to field data, large-eddy simulationresults and Rossby-number similarity relationships.     

ALPEX-Simulation

Summary In a project ALPEX-Simulation, sponsored by the Österreichischer Fond zur Förderung der wissenschaftlichen Forschung (FWF), all eight cases of ALPEX-SOP cyclones were numerically simulated with a fine mesh isentropic model of the atmosphere. These numerical simulations in six-hourly intervals allow a deeper insight into the synoptics and dynamics of the cyclogeneses in the Western Mediterranean, especially into the genesis of the two basic types of cyclones: the so-called Überströmungs-type and Vorderseiten-type. In the first phase of cyclogenesis of the Überströmungs-type, the blocking and flow splitting of the cold air due to the Alps and the canalization between the Alps and the Massif Central are important. Cold air flows cyclonically around the western part of the Alps, creating a vorticity maximum at the south western edge of the Alpine, bow and leads also to an enhanced PV. In connection with warm air in the Mediterranean, a strong baroclinic zone is generated. The interaction between the arriving PV maximum in the upper troposphere and the enhanced PV at the bottom leads to cyclogenesis in the Western Mediterranean. In the case of the Vorderseiten-type warm air advection dominates with the exception of a shallow layer of cold air in the inner Po-Valley, which is shielded by the Alpine ridge. A well-pronounced PV maximum builds up and couples with the PV maximum arriving at upper levels, even before the cold air, coming from the north-west, has surrounded the Alps. The cold air only intensifies the development by raising the baroclinity. Therefore, the Vorderseiten-cyclogenesis is an orographically modified cyclogenesis, in the course of which the cyclonic development is triggered by the Alps, whereas the Überströmungs-cyclogenesis is an orographically induced cyclogenesis i.e. a true lee cyclogenesis.With 14 FiguresDied in a tragic traffic accident on June 6, 1993.     

The behaviour of aerosol optical depth during clear summer days in Athens,Greece

Summary The integral aerosol optical depths (k ) at the hour of 08:20 Local Standard Time (LST), are compared with those calculated previously at 11:20 and 14:20 LST, for clear days during summer in Athens over the period 1962–1988. The mean values at 08:20 LST were consistently lower than the values at 11:20 and 14:20 LST. The influence of the vertical wind profile on the values ofk was also investigated. A comparison was made of the wind profiles at 02:00 and 14:00 LST, for days in which the 11:20 and 14:20 LST values ofk were 0.200 andk 0.350, respectively. The corresponding bulk wind shear s was also found for the period 1980–1988. The most significant results occurred with the first category of days. The resultant wind velocities from the surface to the 900 hPa level, in each hour were higher by 2–4 m·s^–1 with respect to the corresponding values for the second category. At 02:00 LST the bulk wind shear showed a considerable difference (1.8) between the two categories of days in the surface to 700 hPa layer at 02:00 LST. Finally, the associated weather conditions that appear to initiate a period of low values ofk (k 0.200) at 11:20 and 14:20 LST were examined for the period 1980–1988. Fifteen such cases were identified and it was found that they all occurred after the passage of weak cold fronts.With 6 Figures     

The validation of various schemes for parameterizing evaporation from bare soil for use in meteorological models: A numerical study using in situ data

Parameterization of evaporation from a non-plant-covered surface is very important in the hierarchy strategy of modelling land surface processes. One of the representations frequently used in its computation is the resistance formulation. The performance of the evaporation schemes using the , , and their combination resistance approaches to parameterize evaporation from bare soil surfaces is discussed. For that purpose, the nine schemes, based on a different dependence of and on volumetric soil moisture content and its saturated value, are used.The tests of performances of the considered schemes are based on time integrations by the land surface module (BARESOIL) using observed data. The 23 data sets at a bare surface experimental site in Rimski anevi, Yugoslavia on chernozem soil, were used for the resistance algorithm evaluation. The quality of the schemes was compared with the observed values of the latent heat flux using several statistical parameters.     

Flux-profile Relationships for Wind Speed and Temperature in the Stable Atmospheric Boundary Layer

Wind and temperature profiles in the stable boundary layer were analyzed in the context of MoninObukhov similarity. The measurements were made on a 60-m tower in Kansas during October 1999 (CASES-99). Fluxprofile relationships, obtained from these measurements in their integral forms, were established for wind speed and temperature. Use of the integral forms eliminates the uncertainty and accuracy issues resulting from gradient computations. The corresponding stability functions, which were nearly the same for momentum and virtual sensible heat, were found to exhibit different features under weakly stable conditions compared to those under strongly stable conditions. The gradient stability functions were found to be linear, namely _m = 1+ 5.8 and _h = 1 + 5.4 up to a limit of the MoninObukhov stability parameter = 0.8; this is consistent with earlier findings. However, for stronger stabilities beyond a transition range, both functions were observed gradually to approach a constant, with a value of approximately 7. To link these two distinct regimes, a general but pliable functional form with only two parameters is proposed for the stability functions, covering the entire stability range from neutral to very stable conditions.     

Application of a “Big-Tree” model to regional climate modeling: a sensitivity study

Summary ¶In order to better understand land-atmosphere interactions and increase the predictability of climate models, it is important to investigate the role of forest representation in climate modeling. Corresponding to the big-leaf model commonly employed in land surface schemes to represent the effects of a forest, a so called big-tree model, which uses multi-layer vegetation to represent the vertical canopy heterogeneity, was introduced and incorporated into the National Center for Atmospheric Research (NCAR) regional climate model RegCM2, to make the vegetation model more physically based. Using this augmented RegCM2 and station data for China during 1991 Meiyu season, we performed 10 experiments to investigate the effects of the application of the big-tree model on the summer monsoon climate.With the big-tree model incorporated into the regional climate model, some climate characteristics, e.g. the 3-month-mean surface temperature, circulation, and precipitation, are significantly and systematically changed over the model domain, and the change of the characteristics differs depending on the area. Due to the better representation of the shading effect in the big-tree model, the temperature of the lower layer atmosphere above the plant canopy is increased, which further influences the 850hPa temperature. In addition, there are significant decreases in the mean latent heat fluxes (within 20–30W/m²) in the three areas of the model domain.The application of the big-tree model influences not only the simulated climate of the forested area, but also that of the whole model domain, and its impact is greater on the lower atmosphere than on the upper atmosphere. The simulated rainfall and surface temperature deviate from the originally simulated result and are (or seem to be) closer to the observations, which implies that an appropriate representation of the big-tree model may improve the simulation of the summer monsoon climate.We also find that the simulated climate is sensitive to some big-tree parameter values and schemes, such as the shape, height, zero-plane displacement height and mixing-length scheme. The simulated local/grid differences may be very large although the simulated areal-average differences may be much lower. The area-average differences in the monthly-mean surface temperature and heat fluxes can amount to 0.5°C and 4W/m², respectively, which correspond to maximum local/grid differences of 3.0°C and 40W/m² respectively. It seems that the simulated climate is most sensitive to the parameter of the zero-plane displacement among the parameters studied.     

Turbulence Reynolds number and the turbulent kinetic energy balance in the atmospheric surface layer

The relation between the turbulence Reynolds numberR and a Reynolds numberz* based on the friction velocity and height from the ground is established using direct measurements of the r.m.s. longitudinal velocity and turbulent energy dissipation in the atmospheric surface layer. Measurements of the relative magnitude of components of the turbulent kinetic energy budget in the stability range 0 >z/L 0.4 indicate that local balance between production and dissipation is maintained. Approximate expressions, in terms of readily measured micrometeorological quantities, are proposed for the Taylor microscale and the Kolmogorov length scale .     

Experiments On Buoyant Plume Dispersion In A Laboratory Convection Tank

Plume dispersion in the convective boundary layer (CBL) is investigated experimentally in a laboratory convection tank. The focusis on highly-buoyant plumes that loft near or become trapped in the CBL capping inversion and resistdownward mixing. Such plumes are defined by dimensionless buoyancy fluxes F_* 0.1, where F_* = F_b/(U w_* ² z_i), F_b is the stack buoyancy flux,U is the mean wind speed, w_* is the convective velocity scale, and z_i is the CBL depth. The aim is to obtain statistically-reliable mean (C) and root-mean-square (rms, _c) concentration fields as a function of F_* and the dimensionless distance X = w_*x/(U z_i), where x is the distance downstream of the source.The experiments reveal the following mainresults: (1) For 3 X 4and F_* 0.1, the crosswind-integrated concentration (CWIC) fields exhibit distinctly uniform profiles below z_i with a CWIC maximum aloft, in contrast to the nonuniform profiles obtained earlier by Willis and Deardorff. (2) The lateral dispersion (_y) variation with X is consistent with Taylor's theory for _* 0.1 and a buoyancy-enhanced dispersion, _y/z_i F_* ^1/3X^2/3, forF_* = 0.2 and 0.4. (3) The entrapment, the plume fraction above z_i, has a mean (E) that follows a systematic variationwith X and F_*, and a variability (_e/E) that is broad ( 0.3 to 2) near the source but subsides to 0.25 far downstream. (4) Vertical profiles of the concentration fluctuation intensity (c/C) are uniform for z < z_i and X > 1.5, but exhibit significant increases: (a) at the surface and close to the source (X 1.5), and(b) in the entrainment zone. (5) The cumulative distribution functions (CDFs) of the scaled concentration fluctuations (c/_c) separate into mixed-layer and entrainment-layer CDFs for X 2, with the mixed-layer group collapsing to a single distribution independent of z.These are the first experiments to obtain all components of the lateral and vertical dispersion parameters (rms meander, relative dispersion, total dispersion) for continuous buoyant releases in a convection tank. They also are the first tank experiments to demonstrate agreement with field observations of: (1) the scaled ground-level concentration along the plume centreline, and (2) the dimensionless lateral dispersion _y/z_i of buoyant plumes.     

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.     

Charakteristische Zirkulationstypen in mittleren Breiten der nördlichen Hemisphäre

Zusammenfassung In der vorliegenden Arbeit werden neue Klassifikationsprinzipien für Großwetterlagen entwickelt. Bisher wurde bei Wetterlagenklassifikationen das Druckfeld zugrunde gelegt, wobei quasistationären Druckzentren eine nicht berechtigte Vorrangstellung eingeräumt wurde. In der hier versuchten Klassifikation wird vom Strömungsfeld ausgegangen, das in elementare Formen zerlegt wird. Eine zu diesem Zweck durchgeführte statistische Untersuchung ergab, daß alle im Strömungsfeld auftretenden Zirkulationstypen auf drei Grundformen zurückgeführt werden können. Diese Grundformen sind:Driften, Wellen undWirbel.Die Untersuchung ergab im einzelnen, daß in mittleren Breiten der nördlichen Hemisphäre bei 49% aller untersuchten Fälle Driften, bei 23% Wellen und bei 28% Wirbel auftraten.In der hier durchgeführten Klassifikation wird das Druckfeld durch das Strömungsfeld und der Begriff Großwetterlage durch den umfassenderen Begriff des Zirkulationstyps ersetzt. Damit wird der unberechtigte Vorrang der Druckformen bei der Wetterlagenklassifikation aufgegeben. Die Klassifizierung der Zirkulationstypen ergibt sich schließlich durch Kombination der drei Zirkulationselemente: Drift, Welle und Wirbel.

---

Summary New principles of classification for large-scale weather situations are outlined in this paper. Hitherto the pressure-field has been taken as a basis for such classifications by conceding a precedence of an unjustified position to semi-permanent centres of pressure. The new classification starts from the field of large-scale motions, which is dissected in elementary models. A statistical test yielded the possibility to reduce all types of atmospheric circulations in the following three elementary models:drifts, waves andeddies.In detail it was found out, that drifts occur in 49%, waves in 23% and eddies in 28% of all cases investigated.In the new classification the term pressure-field is substituted by field of motion and the expression large-scale weather situation by the more comprehensive conception type of circulation. By that the unjustified priority of pressure-centers in classifying weather situations is abolished. At last the classification of the types of circulation follows from a combination of the three elementary models: drift, wave and eddy.

---

résumé La présente étude développe de nouveaux principes de classification des situations météorologiques. Alors que jusqu'ici on s'est fondé sur le champ de pression ce qui conduisait à attribuer aux centres d'action quasi stationnaires un rôle trop important, l'auteur part ici du champ de mouvement décomposé en formes élémentaires. Un examen statistique lui a montré que tous les types de circulation peuvent se ramener à trois formes fondamentales:courants, ondulations ettourbillons.Aux latitudes moyennes de l'hémisphère Nord les courants représentent le 49%, les ondulations le 23% et les tourbillons le 28%.Au champ de pression se substitue donc le champ de courant, et les situations météorologiques se groupent en types de circulation ce qui supprime le rôle prépondérant des formes isobariques. Le classement final des types de circulation résulte de la combinaison des trois types mentionnés:courants, ondulations ettourbillons.

---

---

Mit 15 Textabbildungen     

Arctic hazes in summer over Greenland and the North American Arctic: II. Nature and concentrations of accumulation-mode and giant particles

Airborne measurements made during August 1985 over Greenland and its environs show that both accumulation-mode (0.1 m D2.0 m) and giant (D2 m) particles were present in relatively high concentrations in arctic haze layers and that the accumulation-mode particles dominated light scattering. Particles with diameters (D) between 1 and 4 m consisted predominately of mixed materials, small and dense inclusions, and probably organic compounds containing sulfur. Many of the particles from 0.1 to 1 m in diameter were also of mixed composition, with sulfuric acid, ammonium sulfate and organics probably the dominant constituents.     

Bioaccumulation of Cobalt in Parmelia sulcata

Six locations across mainland Portugal were selected for exposing Parmelia sulcata, for a one-year period (8 months for one site), with simultaneous measurement of total (dry + wet) deposition (one-month periods). The exposed lichens and the total (dry + wet) deposition were analysed for cobalt contents by INAA (instrumental neutron activation analysis) and ICP-MS (inductively coupled plasma mass spectroscopy), respectively. The designated wet deposition was evaluated through the collected water volume; the designated dry deposition was assessed after the (dried) residual mass of the wet deposition. An excellent agreement between Co contents in exposed lichens and the cumulative (1) Co contents in the dry deposition, (2) dry deposition, and (3) wet deposition has been found for the locations with alternate drought and precipitation months, high dry deposition, and high Co contents in the latter. Continuous rainfall was found to hinder the Co accumulation in the lichen due to its release from the lichen and/or lower Co contents in the dry deposition. At three locations, P. sulcata Co contents, after subtraction of the background (before exposure), equalled or exceeded the Co contents in the cumulative dry deposition at the end of the exposure time. The optimal exposure period for this species likely depends on the exposure conditions.     

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.     

Gas and Particle Products Distribution from the Reaction of β-Caryophyllene with Ozone

Gas and particulate reaction products from the ozonolysis of -caryophyllene (I) in the presence of atmospheric air were investigated using a combination of gas chromatography-mass spectrometry (GC-MS) and high performance liquid chromatography (HPLC). A Scanning Mobility Particle Sizer system (3936, TSI) and a Condensation Particle Counter (3025A, TSI) were used to study secondary organic aerosol formation. The nighttime oxidation was carried out in a large outdoor smog chamber (190 m³). A wide range of ring retaining and ring opening products in the gas and particle phase are reported over the course of the reaction. On average, measured gas and particle phase products accounted for 64% of the reacted -caryophyllene (I) carbon. Measurements show that a number of reaction products with low vapor pressure (e.g., -caryophyllone aldehyde (IV), -norcaryophyllone aldehyde (V), -caryophyllonic acid (VIII), -14-hydroxycaryophyllonic acid (XIV)) were found in the sample taken during the first 20 min of the reaction and may play an important role in the early formation of secondary organic aerosol. A detailed mechanism is proposed to account for most products observed in this investigation.     

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.