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     

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.     

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     

The Askervein Hill Project: Mean wind variations at fixed heights above ground

This is one of a series of papers on the Askervein Hill Project. It presents results on the variations in mean wind speed at fixed heights (z) above the ground from linear arrays of anemometer posts and towers. Most of the data are for z = 10 m but some are for z = 3 m. Selected and directionally grouped data from the 55 Mean Flow runs are presented together with mean flow data from Askervein '83 Turbulence runs. Comparisons are made between the data and guideline estimates of fractional speed-up ratio at hilltop locations and between the data and MS3DJH/3 model predictions along the tower lines. There is good agreement in most cases.     

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.     

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

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

A radiation and energy balance study of mature forest and clear-cut sites

Components of the radiation and energy balances were measured over a clear-cut area and a mature, mixed forest during the summer of 1981 at the Petawawa National Forestry Institute, Chalk River, Ontario. The work concentrated on the clear-cut site which supported a canopy layer composed primarily of bracken fern and logging remnants.Forty days of radiation data were collected at the clear-cut site. After the first four weeks of measurements (the green season), most of the ferns quickly died, and their foliage changed appearance from a green to brownish colour (the brown season). The daily mean reflection coefficient of solar radiation determined over the green season was 0.20 and decreased to 0.13 for the brown season. The corresponding value for the forest was 0.13, based on a limited amount of data. The clear-cut site received 11% and 21% less net radiation than the forest on a 24-hr and daylight-hours basis, respectively, as a consequence of the higher reflection coefficient and larger daytime longwave radiation emission.A reversing temperature difference measurement system (RTDMS), incorporating ten-junction thermopiles was employed at each site in order to determine Bowen ratios () via differential psychrometry. Both systems performed well, especially the RTDMS over the forest which was capable of resolving very small differences of temperature, typically less than 0.2 °C over a height of 3 m. The mean hourly Bowen ratio, calculated from values from 0800 to 1600 hr, varied from 0.2 to 1.0 for the forest and from 0.4 to 0.8 for the clear-cut site in the green season.A significant canopy heat storage component of the energy balance, Q _S , was found at the clear-cut site. In the early morning, a portion of the available energy was used to heat the biomass materials and air within the canopy layer. The stored heat within the canopy was released later in the day, increasing the available energy total.The daily mean value of the Priestley-Taylor coefficient (Priestley and Taylor, 1972) for the green season at the clear-cut site was 1.14, and individual values tended to increase during wet surface conditions and decrease when the surface dried. The daylight mean value during dry canopy conditions at the forest was 1.05, and much higher values occurred when the canopy was wet. The enhancement of for the wet forest was a result of the evaporation of intercepted rain (which is not limited by stomatal resistance) and the concomitant transfer of sensible heat to the forest.     

Shear convection

The paper discusses convection in the presence of wind shear, a condition analysed previously by Zilitinkevich (1971). This region (between the forced and fully developed convective layers) was also considered by Betchov and Yaglom (1971). In the present paper the author endeavors to develop a consistent analysis from the basic hypothesis of a very weak interaction between the vertical convective motions and mechanical turbulence, employing a new similarity model of the turbulent regime. Additional experimental data are introduced. Unlike the notation used in the references quoted above, this regime is termed shear convection rather than free convection. The latter is traditionally regarded as synonymous with the terms pure or fully developed convection.     

Ozone production due to emissions from vegetation burning

Ozone has been observed in elevated concentrations by satellites over areas previously believed to be background. There is meteorological evidence, that these ozone plumes found over the Atlantic Ocean originate from vegetation fires on the African continent.In a previous study (DECAFE-88), we have investigated ozone and assumed precursor compounds over African tropical forest regions. Our measurements revealed large photosmog layers at altitudes from 1.5 to 4 km. Both chemical and meteorological evidence point to savanna fires up to several thousand km upwind as sources.Here we describe ozone mixing ratios observed over western Africa and compare ozone production ratios from different field measurement campaigns related to vegetation burning. We find that air masses containing photosmog ingredients require several days to develop their oxidation potential, similar to what is known from air polluted by emissions from fossil fuel burning. Finally, we estimate the global ozone production due to vegetation fires and conclude that this source is comparable in strength to the stratospheric input.     

Carbon dioxide and the δ18O record of late-Quaternary climatic change: a global model

A dynamical model for the late-Quaternary global variations of ¹⁸O, mean ocean surface tempeature , ice mass I, deep ocean temperature , and atmospheric carbon dioxide concentration , is constructed. This model consists of two diagnostic equations (for ¹⁸O and ), and three prognostic equations (for I, , and ) of a form studied extensively in previous articles. The carbon dioxide equation includes forcing by a representation of the Milankovitch earth-orbital radiation effects, and contains a basic instability that drives a free oscillation of period near 100,000 years. The system is constrained to conserve mass and energy, contain physically plausible feedbacks including a system time constant no greater than 10.000 years, and be robust (i. e., structurally stable in the presence of expected noise levels and uncertainties in values of coefficients). Within the limits of these constraints, coefficients are chosen such that (i) the solution gives a good fit to the observed SPECMAP ¹⁸O variations, and (ii) the ice mass variations are qualitatively similar to the ¹⁸O variations. The predicted long term variations of sea surface temperature and atmospheric carbon dioxide are in reasonably good agreement with the limited observational evidence available for these quantities, while the predicted variations of deep ocean temperature remain to be verified when paleoclimatic estimates of this quantity become available. The relative contributions of ice mass changes and surface water temperature changes to the variations of ¹⁸O at any time are given by the model.     

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

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

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     

The influence of surface cover and climate on energy partitioning and evaporation in a subarctic wetland

Energy partitioning and evaporation were measured over three wetland surfaces in a subarctic coastal marsh during pre-growing and growing periods. These surfaces included an alder/willow woodland, a sedge marsh and a raised backshore sedge meadow. A combination model analysis was used to assess the relative importance of surface resistance and meteorological conditions on the magnitude of the Bowen ratio, , during the growing period.Overall, the three surfaces experienced important site-to-site and seasonal differences in and evaporation, Q _E. During the non-foliated period, Q _E was largest and was smallest for the open water marsh, while the dry backshore site experienced the smallest Q _E and largest . The non-foliated woodland assumed intermediate values of and Q _E. After the vegetation covers were established, the woodland assumed the smallest and largest Q _E flux. It was also found that at the marsh site increased with the presence of a vegetation cover.Wind direction was always an important factor in determining Q _E and at all sites. was substantially larger and Q _E was smaller for onshore winds (i.e., originating from James Bay) than for offshore winds. The combination model analysis showed that canopy resistance at all sites was largest during warm offshore winds, which were associated with large saturation deficits. However, the effect of increased canopy resistance on during offshore winds was offset by a large climatological resistance, resulting in small values and large Q _E. When winds originated from James Bay, canopy resistance was smaller than for offshore winds, but the climatological resistance also was much smaller, resulting in larger and small Q _E. The results have important implications for changes in land cover and climate on the regional water balance.     

Flux Footprints Within and Over Forest Canopies

The characteristics of turbulence within a forest arespatially heterogeneous and distinct from thoseassociated with the surface boundary layer. Consequently, the size and probability distribution offlux footprints emanating from sources below aforest canopy have the potential to differ from thoseobserved above forests.A Lagrangian random walk model was used toinvestigate this problem since no analytical solutionof the diffusion equation exists. Model calculations suggest that spatialcharacteristics of flux footprint probabilitydistributions under forest canopies are muchcontracted, compared to those evaluated in the surfaceboundary layer. The key factors affecting thestatistical spread of the flux footprint, and theposition of the peak of its probability distribution,are horizontal wind velocity and the standarddeviations of vertical and horizontal velocityfluctuations. Consequently, canopies, which attenuatemean horizontal wind speed, or atmospheric conditions,which enhance vertical velocity fluctuations, willcontract flux footprint distributions mostly near thefloor of a forest. It was also found that theprobability distributions of the flux footprint arenarrower when horizontal wind velocity fluctuationsare considered, instead of the simpler case that considers only vertical velocity fluctuations and meanhorizontal wind velocity.     

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.     

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.     

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.     

An observational study of heat fluxes and their relationships with net radiation

Experimental evidence indicates that the diurnal behaviour of the fluxes of heat into the ground and into the atmosphere versus the net flux of radiation can be modelled by closed curves, the hourly values folowing one another in either a clockwise or counter clockfashion. A general formulation to express the different heat fluxes as a function of net radiation is proposed. This relationship between the different heat fluxes and can be expressed as a sum of three terms: the first indicates a direct proportionality to , the second gives the deviation from linearity and depends on /t, and the third gives the value of the flux when = 0. The formulae are then expressed as a function of time and the ratios between the heat fluxes and are evaluated. A comparison with the approximations generally used shows that the latter may be considered as particular cases of the more general equations proposed here.     

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.     

Regular and seasonal behaviour of atmospheric radio noise field strength (ARNFS) over low latitude station calcutta

Summary The electromagnetic radiation of cloud discharge known as atmospheric radio noise field strength (ARNFS) shows a gradual fall from a frequency of 9 kHz to 81 kHz as studied over a period of two years at Calcutta, very close to Bay of Bengal. The main characteristic features of ARNFS at Calcutta are that-(i) ARNFS shows that midday median value is smaller than midnight median value in all months, (ii) level of daily minimum is higher in February and monsoon compared to other seasons, (iii) sunrise effect and sunset effect are well correlated with local sunrise and sunset times, (iv) the magnitude of sunrise fade and sunrise fade rate are maximum in April and lowest during winter period, (v) the magnitude of sunset fade is higher in premonsoon and postmonsoon while it is lowest in monsoon, (vi) number of occurrence of both sunrise effect and sunset effect is remark-ably smaller in monsoon. The positions of the sun and of atmospheric sources are jointly the causes of seasonal and diurnal variations. The missing of sunrise effect and sunset effect are due to local cloud activity and variation of electron density during geomagnetic storms.With 7 Figures