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     

Breaking global temperature records after Mt. Pinatubo

The probability of setting global temperature records is reconsidered in light of cooling due to the Mt. Pinatubo eruption. The cooling resets temperature by moving it away from the top of its 100-year range. Depending on the statistical model for temperature, this reset can lead to a much lower probability for a record in the next few years. The exercise illustrates how record setting depends on the underlying model, the current record value, and - if there is serial correlation- the current temperature.     

On Modelling the One-Dimensional Atmospheric Boundary Layer

Several commonly used turbulence closure schemes forthe atmospheric boundarylayer (ABL) are applied to simulate neutral, nocturnal and diurnal cyclesituations in a one-dimensional ABL. Results obtained with the differentschemes, E-, E- and its modified versions, and twoversions ofthe q ² Level 2.5, are compared and discussed.     

A Comparative Analysis of Transpiration and Bare Soil Evaporation

Transpiration E^v and bare soil evaporation E^b processes are comparatively analysed assuming homogeneous and inhomogeneous areal distributions of volumetric soil moisture content . For a homogeneous areal distribution of we use a deterministic model, while for inhomogeneous distributions a statistical-deterministic diagnostic surface energy balance model is applied. The areal variations of are simulated by Monte-Carlo runs assuming normal distributions of .The numerical experiments are performed for loam. In the experiments we used different parameterizations for vegetation and bare soil surface resistances and strong atmospheric forcing. According to the results theE^v()-^Eb() differences are great, especially in dry conditions. In spite of this, the available energy flux curves of vegetation A^v() and bare soil A^b() surfaces differ much less than the E^v() and ^Eb() curves. The results suggest that E^v is much more non-linearly related to environmental conditions than ^Eb. Both E^v and ^Eb depend on the distribution of , the wetness regime and the parameterization used. With the parameterizations, ^Eb showed greater variations than E^v. These results are valid when there are no advective effects or mesoscale circulation patterns and the stratification is unstable.     

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.     

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.     

AN E – ε – ℓ TURBULENCE CLOSURE SCHEME FOR PLANETARY BOUNDARY-LAYER MODELS: THE NEUTRALLY STRATIFIED CASE

The standard E – model generates aplanetary boundary layerthat appears to be much too deep. The cause of theproblem is traced to the equation for the dissipationrate () of turbulent kinetic energy (E), specifically theparameterization of dissipation production anddestruction. In the context of atmosphericboundary-layer modelling, we argue that a part of thedissipation production should be modelled as the inputto the spectral cascade from the energy-containingpart of the spectrum, with a characteristic length , while the equilibrium imbalancebetween local production and destruction ofdissipation is modelled as proportional toE2/E, as in the standard model. Wepropose an E – – turbulence closurescheme, in which both the mixing length, m, and are prescribed. The importance ofthe equation is diminished, though itstill determines the dissipation rate in the Eequation.     

Equilibrium theory of the planetary boundary layer with an inversion lid

The maintenance of an elevated inversion in steady flow above a cold, rotating surface is shown to be possible for a certain range of the buoyancy number bfV _g, where b is the buoyant acceleration appropriate to the density deficiency of the fluid above the inversion, f is Coriolis parameter and V _gis geostrophic velocity (so that fV _gis also horizontal pressure gradient in kinematic units). The height of the inversion lid is determined by a balance of surface stress and buoyancy, in a way which may be deduced from laboratory experiments. With the aid of such empirical evidence a theory is constructed for the layer below the inversion lid. The cross-isobar angle of ground-level stress is found to increase with the buoyancy number, to a limiting value of 90, by which time the inversion descends to the ground. Under typical conditions, a temperature difference of order 10C is necessary to eliminate the possibility of an equilibrium, elevated inversion lid and reduce ground level wind stress to a vanishingly small value.Woods Hole Oceanographic Institution Contribution #3011On leave from the University of Waterloo, Ontario     

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.     

Gas and Particulate Products Distribution from the Photooxidation of α-Humulene in the Presence of NOx, Natural Atmospheric Air and Sunlight

The photooxidation of -humulene in the presence of NO_x, natural sunlight, and rural background air was investigated using a combination of gas chromatography-mass spectrometry (GC-MS) and high performance liquid chromatography (HPLC). Identification and quantification of gas and particulate reaction products were reported over the course of the reaction. The daytime photooxidation 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. On average, measured gas and particle phase products accounted for 44% of the reacted -humulene carbon. Measurements show that a number of reaction products with low vapor pressures (e.g. 3-seco--humulone aldehyde, 7-seco--humulone aldehyde, -humulal aldehyde, -humulene 3-oxide or -humulene 7-oxide, -humulaic/alic acid isomers, and 3-seco--14-hydroxyhumulone aldehyde) were found in the early stage 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.     

Similarity of atmospheric Reynolds shear stress and heat flux fluctuations over a rough surface

Horizontal u and vertical w velocity fluctuations have been measured together with temperature fluctuations in the atmospheric surface layer, at a small height above a wheat crop canopy. Marginal probability density functions are presented for both individual fluctuations u, w, and for the instantaneous Reynolds stress uw, and heat fluxes w and u. Probability density functions of the velocity fluctuations deviate less significantly from the Gaussian form than the probability density of temperature. There appears to be closer similarity between statistics of the instantaneous heat fluxes than between the momentum flux and either of the heat fluxes investigated. The mean momentum flux receives equal contributions from the events referred to as ejections and sweeps in laboratory boundary layers. Sweeps provide the largest contribution to the heat fluxes.     

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.     

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.     

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     

Estimating product yields of carbon-containing products from the atmospheric photooxidation of ambient air alkenes

The yields of products have been calculated for the reactions of hydroxyl radicals and ozone with 19 of the two-through-six carbon anthropogenic alkenes. Based on their rate of reaction, mechanisms of reactions and the ambient air distribution for these alkenes their seasonal ambient air yields have been estimated.Aldehydes predominate as products irrespective of season, with smaller yields of several ketones. Other minor products include carboxylic acids, carbon monoxide, carbon dioxide, and alkenes. About a two-fold increase is estimated in the yields of hot biradicals and their products from summer to winter.One sensitivity analysis was made by recomputing yields at a different OH radical to O₃ concentration than assumed most likely in the calculations discussed above. In addition, the sensitivity of product yields to an estimated range of seasonally averaged sunset-to-sunrise NO₃ radical concentrations was calculated. The effects of free radical reactions are discussed, but these are believed to make a relatively minor contribution within the NO _x -rich atmospheres that contain anthropogenic alkenes.The uncertainties in product yields associated with the range of NO₃ radical concentrations assumed present is relatively small for aldehydes, as is the decrease in yield of the one carbon hot biradical. Larger uncertainties occur for ketones. Significant decreases in yields occur for larger hot biradicals, especially the branched-chain hot radicals in the presence of NO₃ radicals.     

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.     

Calculation of velocity skewness in real and artificial plant canopies

Wind velocities within a plant canopy are much more strongly skewed than those of the air flow above. We have examined the governing Eulerian equations for the velocity products u _i, u _j u_k using data from a wind tunnel study with an artificial canopy consisting of an array of 5 cm lengths of monofilament fishing line, and from measurements in corn (Zea mays L).Simple parameterizations for pressure-velocity correlations, and for the quadruple velocity products allowed reasonably accurate calculations of the third moments using measured profiles of the mean velocity, variance and covariance fields. Comparisons of individual terms in the rate equations for ovu _i, u _j u _krevealed that diffusion (from above) and mean shear were most important in creating large skewness in the canopy. A drag term also contributed but was of lesser importance. These terms were balanced by return-to-isotropy and a turbulence interaction term. A quasi-Gaussian approximation considerably underestimated the magnitude of the fourth moments within the canopy.     

Effect of finite sampling on atmospheric spectra

The effect of a finite averaging time on variances is well known, but its effect on power spectra is less clearly understood. We present numerical solutions for the spectral distortion arising from sampling over a finite time interval T and show that the commonly used filter function (1 – sinc²f T), valid for variances, is a reasonable approximation for power spectra only when T 10 _m , where f is the cyclic frequency, and _m is the dominant time scale of the process. Our results exhibit an increasingly steeper low-frequency roll-off as T decreases relative to _m , indicating that the measured spectrum is subject to a greater suppression of the lower frequencies (f > 1/T) than predicted by (1 – sinc²f T). This suppression is, in a sense, compensated by an overestimation of spectral estimates in the frequency range f 1/T.     

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     

Atmospheric effects in the remote sensing of phytoplankton pigments

We investigate the accuracy with which relevant atmospheric parameters must be estimated to derive phytoplankton pigment concentrations (chlorophyll a plus phaeophytin a ) of a given accuracy from measurements of the ocean's apparent spectral radiance at satellite altitudes. The analysis is limited to an instrument having the characteristics of the Coastal Zone Color Scanner scheduled to orbit the Earth on NIMBUS-G. A phytoplankton pigment algorithm is developed which relates the pigment concentration (C) to the three ratios of upwelling radiance just beneath the sea surface which can be formed from the wavelengths () 440, 520 and 550 nm. The pigment algorithm explains from 94 to 98% of the variance in log₁₀ C over three orders of magnitude in pigment concentration. This is combined with solutions to the radiative transfer equation to simulate the ocean's apparent spectral radiance at satellite altitudes as a function of C and the optical properties of the aerosol, the optical depth of which is assumed to be proportioned to ^-n . A specific atmospheric correction algorithm, based on the assumption that the ocean is totally absorbing at 670 nm, is then applied to the simulated spectral radiance, from which the pigment concentration is derived. Comparison between the true and derived values of C show that: (1) n is considerably more important than the actual aerosol optical thickness; (2) for C 0299-1 0.2 g l^-1 acceptable concentrations can be determined as long as n is not overestimated; (3) as C increases, the accuracy with which n must be estimated, for a given relative accuracy in C, also increases; and (4) for C greater than about 0.5 g 1^-1, the radiance at 440 nm becomes essentially useless in determining C. The computations also suggest that if separate pigment algorithms are used for C 1gl^-1 and C 1 gl^-1, accuracies considerably better than ±± in log C can be obtained for C 1 g l^-1 with only a coarse estimate of n, while for C 10 gl^-1, this accuracy can be achieved only with very good estimates of n.Contribution No. 387 from the NOAA/ERL Pacific Environmental Laboratory.On leave from Department of Physics, University of Miami, Coral Gables, Florida.