Recurrence statistics of concentration fluctuations in plumes within a near-neutral atmospheric surface layer

This study examines the statistical properties of the concentration derivative, , for a dispersing plume in a near-neutrally stratified atmospheric surface layer. Towards this goal, the probability density function (pdf) of , and the conditional pdf of given a fixed concentration level, , have been measured. These pdfs are found to be modeled well by a generalizedq-Gaussian (gqG) distribution with intermittency exponent,q, equal to 0.3 and 3/4, respectively. These results highlight the strong intermittency effect (patchiness) of the small-scale concentration eddy structures in the plume. The distribution of time intervals between successive high peaks in the squared derivative process, x², is found to be well approximated by a power-law distribution, implying that occurrences of these high peaks are much more clustered than would be predicted by a Poisson or shot-noise process. The results are used to improve models for the joint pdf of and , and for the expected number of upcrossings per unit time interval of a fixed concentration level that have been proposed by Kristensenet al. (1989). The predictions of the improved models are in accord with observations, and suggest that the intercorrelation between and must be explicitly incorporated if good estimates of the upcrossing intensity are to be obtained.     

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.     

Similarity forms for ground-source surface-layer diffusion

Past work on analyzing ground-source diffusion data in terms of surface-layer similarity theory is reviewed; these analyses assume that _z /L orh/L is a function of u ^* x/L (where h = Q/ dy). It is argued that an alternative scaling, h ^*/L versus x/L, is nearly as universal in that it is very weakly influenced by surface roughness, except for a modest influence in the free convective case (h ^* = Q/u ^* dy); the advantage of this scaling is that it eliminates the need to reassess as vertical diffusion progresses. The Prairie Grass data set is adjusted for the difference in source and sampling heights, and is plotted with this scaling. Simple analytic equations are suggested that fit the resultant data plots for stable and unstable conditions, and suggestions are made towards practical application of these results.On assignment from the National Oceanic and Atmospheric Administration, U.S. Department of Commerce.     

Observed variations of σθ and σΦ in the nocturnal drainage flow in a deep valley

The variations of and in the drainage flow in the Brush Creek valley of western Colorado are investigated using data from Doppler acoustic sodars and instrumented towers. The data were obtained on two experimental nights during the 1984 ASCOT field study. There is good agreement between the variations derived from low-level observations of the sodars and those derived from the towers located throughout the valley. The observed hourly average and in the nocturnal drainage flow are about 20 ° to 25 ° and 5 °, respectively; these values are much larger than those generally observed over flat terrain during nighttime stable conditions. After sunrise (about 0600 MST), as the valley warms and the flow direction changes to up-valley, these parameters increase sharply to their peak values at about 0800 MST and then decrease to their normal daytime values after about two hours.In the drainage flow, the hourly average varies inversely with wind speed according to the relation u 0.7ms^-1. The vertical standard deviation is much less enhanced by complex terrain than the horizontal standard deviation. The observed values are predicted fairly well by the local similarity theory.Oak Ridge Associated Universities (ORAU) Summer Research Participant at ATDD in 1987 andOak Ridge Associated Universities (ORAU) Summer Research Participant at ATDD in 1987 and     

Turbulence mechanisms at an agricultural site

An extensive set of turbulence data from the 3- and 12-m heights taken over an agricultural site (Marsta, Sweden) are analyzed and compared with data from ideal sites.In unstable air, Monin-Obukhov similarity is found to be valid for the non-dimensional gradients of wind, _m , temperature, _h , and humidity, _e , for (only a few data), for _T /|T _*|,/ _E /|E _*| and for the non-dimensionalized inertial subrange spectra of temperature and humidity. Where comparison is possible, the unstable data also agree with those found in the Kansas study, with one remarkable exception, the inertial subrange constant of the temperature spectrum, ₁, being only 0.39, compared to the value 0.80 found at the Kansas site.On the stable side, most similarity predictions break down, with most of the data differing systematically from the corresponding Kansas results, the only exception being . The inertial subrange constants for temperature, ₁, and for humidity, ₁ are found to have the same values, 0.39 and 0.30, respectively, as they do on the unstable side. Remarkable similarity is found for the shape of the stable u- and - and e-spectra. In addition, this shape is found to be identical with that found in Kansas. The peak wavelength of the stable u-, and -spectra is found to be about four times larger than it is for the corresponding Kansas spectra. This is interpreted to be a result of the increased macro-roughness at the Marsta site as compared with that at the Kansas site. A possible explanation for the low ₁-value is discussed, suggesting that ₁ is not a universal constant, but instead dependent on the turbulent structure.     

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     

Occurrence of roll circulations in a shallow boundary layer

Meteorological measurements taken at the Näsudden wind turbine site during slightly unstable conditions have been analyzed. The height of the convective boundary layer (CBL) was rather low, varying between 60 and 300 m. Turbulence statistics near the ground followed Monin-Obukhov similarity, whereas the remaining part of the boundary layer can be regarded as a near neutral upper layer. In 55% of the runs, horizontal roll vortices were found. Those were the most unstable runs, with -z _i/L > 5. Spectra and co-spectra are used to identify the structures. Three roll indicators were identified: (i) a low frequency peak in the spectrum of the lateral component at low level; (ii) a corresponding increase in the vertical component at mid-CBL; (iii) a positive covariance {ovvw} together with positive wind shear in the lateral direction (V/z) in the CBL. By applying these indicators, it is possible to show that horizontal roll circulations are likely to be a common phenomenon over the Baltic during late summer and early winter.     

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.     

Equatorial disturbances,monsoons, and 1982–1983 ENSO

Summary Interannual modes are described in terms of three-month running mean anomaly winds (u,v), outgoing longwave radiation (OLR), and sea surface temperature (T _* ). Normal atmospheric monsoon circulations are defined by long-term average winds (u _n,v _n) computed every month from January to December. Daily winds are grouped into three frequency bands, i.e., 30–60 day filtered winds (u _L,v _L); 7–20 day filtered winds (u _M,v _M); and 2–6 day filtered winds (u _S,v _S). Three-month running mean anomaly kinetic energy (signified asK _L , K _M , andK _S , respectively) is then introduced as a measure of interannual variation of equatorial disturbance activity. Interestingly, all of theseK _L , K _M , andK _S perturbations propagate slowly eastward with same phase speed (0.3 ms^–1) as ENSO modes. Associated with this eastward propagation is a positive (negative) correlation between interannual disturbance activity (K _L , K _M , K _S ) and interannualu (OLR) modes. Namely, (K _L , K _M , K _S ) becomes more pronounced than usual nearly simultaneously with the arrival of westerlyu and negativeOLR (above normal convection) perturbutions. In these disturbed areas with (K _L , K _M , K _S >0), upper ocean mixing tends to increase, resulting in decreased sea surface temperature, i.e.T _* 0. Thus, groups (not individual) of equatorial disturbances appear to play an important role in determiningT _* variations on interannual time scales. HighestT _* occurs about 3 months prior to the lowestOLR (convection) due primarily to radiational effects. This favors the eastward propagation of ENSO modes. The interannualT _* variations are also controlled by the prevailing monsoonal zonal windsu _n, as well as the zonal advection of sea surface temperature on interannual time scales. Over the central Pacific, all of the above mentioned physical processes contribute to the intensification of eastward propagating ENSO modes. Over the Indian Ocean, on the other hand, some of the physical processes become insignificant, or even compensated for by other processes. This results in less pronounced ENSO modes over the Indian Ocean.With 10 FiguresContribution No. 89-6, Department of Meteorology, University of Hawaii, Honolulu, Hawaii.     

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     

Equilibrium evaporation beneath a growing convective boundary layer

Expressions for the equilibrium surface Bowen ratio ( _s ) and equilibrium evaporation are derived for a growing convective boundary layer (CBL) in terms of the Bowen ratio at the top of the mixed layer _i and the entrainment parameter A _R . If A_R is put equal to zero, the solution for _s becomes-that previously obtained for the zero entrainment or closed box model. The Priestley-Taylor parameter is also calculated and plotted in terms ofA _R and _i . Realistic combinations of the atmospheric parameters give values of in the range 1.1 to 1.4.     

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.     

A new method of cluster analysis for numerical classification of climate

Summary In this paper, a new method of cluster analysis incorporating a statistical test, was proposed in order to decide whether to accept or reject the assumption of a numerical classification of climate according to a ² test. On the basis of normality of errors, the ² test has been derived from the variance of the Euclidean distance or from that of the reduced distance.This new method of cluster analysis consists of uniform clustering and central clustering. The former is adapted to a more detailed numerical classification of climate, while the latter is appropriate for the large-scale numerical classification of climate. Examples of calculation and application of this new method of cluster analysis are presented in order to illustrate calculations undertaken in the processes of climate classification.     

Determination of similarity functions of the resistance laws for the planetary boundary layer using surface-layer similarity functions

Functional forms of the universal similarity functions A, B (for wind components parallel and normal to the surface stress), and C (for potential temperature difference) are determined based on the generalized theory of the resistance laws for the Planetary Boundary Layer (PBL). The similarity-profile functions for the surface layer are matched with the velocity and temperature-defect profiles that are assumed to have shapes modified by certain powers of nondimensional height z/h, where h is the PBL height. The powers of the outer-layer profile functions are determined, so that the functions become negligible in the surface layer. To close the temperature defect law, an assumption that the temperature gradient across the top of the PBL is continuous with the stratification of the overlying atmosphere is used. The result of this assumption is that nondimensional momentum and temperature profiles in the PBL can be described in terms of four basic ratios: (1) roughness ratio = /h (2) scale-height ratio =|f|h/u_*, (3) ambient stratification parameter =h/_*, and (4) stability parameter =h/L, where L is the Monin-Obukhov length, z₀ is the surface roughness, is the upper-air stratification, u _* is the friction velocity, and _* is the temperature scale at the surface. For stable conditions, the scale-height ratio can be related to the atmospheric stability and the upperair stratification, and the generalized similarity and Rossby number similarity theories become identical. Under appropriate boundary conditions, function A is explicitly dependent on the stability parameter , while B is a function of scale-height ratio , which in turn depends on the stability. Function C is shown to be dependent on the stability and the upper-air stratification, due to the closure assumption used for the temperature profile.The suggested functional forms are compared with other empirical approximations by several authors. The general framework used to determine the functional forms needs to be tested against good boundary-layer measurements.     

A multifractal analysis of lidar measured water vapour

Accurate and fast-response measurements of space-time observations of specific humidity were made above a drying land surface at the University of California at Davis, using the Los Alamos water Raman-lidar. In an attempt to quantify the space-time intermittency features of turbulent flows in the lower atmosphere, a multifractal analysis of these water vapour measurements was performed. The structure of the specific humidity, (x, t), was analyzed quantifying a scalar gradient measure both in time and space, for all possible one-dimensional cuts, i.e. and . The results confirm the multifractal nature of this scalar gradient measure (a type of scalar dissipation rate) and show that humidity measurements at fixed times (x) are more intermittent (e.g. have less entropy dimension) than those at fixed locations in space (t). Similar multifractal behaviour of the spatial data, with and without a transformation from the observed wind velocities, supports the validity of Taylor's hypothesis for the studied fields.     

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.     

MCC and gull flight behavior

Wintertime observations of mesoscale cellular convection (MCC) over the East China Sea have resulted in criteria that have a remarkable similarity to those reported by Woodcock (1975) in the study of thermals and gull flight behavior. It has been determined that the surface wind speed (V) and the air-sea temperature difference (T) prescribe unique and compatible conditons for both the occurrence of MCC and soaring by sea gulls. Specifically, the onset of MCC when V is between 5 and 9 m s^–1 is inversely proportional to T in the range 5 to 7 °C. Elsewhere, the onset of MCC occurs under conditions of direct proportionality between V and T. Necessary conditions for the occurrence of MCC due to heating from below are T 5 °C and V 5 m s^–1. The boundaries of the convective regime for MCC are discussed and interpreted in accordance with the regime for sea-gull soaring and similarity concepts.     

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.     

Radiation balance in a tropical city — Delhi (India)

The impact of urbanization is assessed by comparing values of the radiation parameters at an urban location with those of a rural site. Urban Delhi was divided according to land-use and the effects of urbanization was studied on incoming short-wave (K), albedo, incoming longwave (L, outgoing longwave (L), and net radiation (Q ^*), were individually studied at four representative sites (Rural, Commercial, Residential and Industrial). MaximumK was observed in the rural and commercial areas whereas highL was observed in the commercial and industrial locations. High depletion ofK of the order of 13% was observed for the industrial location in the winter season. An increase ofL in, the industrial location is of the same order as that of the commercial location, i.e. 20%. The residential location shows quite moderate (4.6%) depletion ofK in comparison with other sites.Nomenclature U Urban - R Rural - K Incoming Short-wave radiation - L Incoming Long-wave radiation - L Outgoing Long-wave radiation - Q ^* Net Radiation - Albedo - K ^* Net Short-wave radiation - L ^* Net Long-wave radiation