Periodicity of annual precipitation in different climate regions of Croatia

Summary The periodicity of a 100-year series of annual precipitation over Croatia has been studied by means of power spectrum analysis at 3 stations representing the different climatic regions of Croatia. The annual precipitation variance spectra in the continental lowland (Osijek) and at the north East Adriatic coast (Crikvenica) can be fitted by Markov white noise continuum, but in the transitional region between the Dinaric Alps and the Pannonian lowland (Zagreb-Gri) a non-white noise continuum is necessary. Quasi-periodic oscillations appear in two spectra ranges: short (2.2 and 4.7 years) and medium (25.0 and 33.3 years). These results are compared with those of other authors for other parts of the Europe.With 2 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 Model to Calculate what a Remote Sensor ‘Sees’ of an Urban Surface

Whilst the measurement of radiation emissions from a surface is relatively straightforward, correct interpretation and proper utilization of the information requires that the surface area seen be known accurately. This becomes non-trivial when the target is an urban surface, due to its complex three-dimensional form and the different thermal, radiative and moisture properties of its myriad surface facets. The geometric structure creates shade patterns in combination with the solar beam and obscures portions of the surface from the sensor, depending on where it is pointing and its field-of-view (FOV). A model to calculate these surface-sensor-sun relations (SUM) is described. SUM is tested against field and scale model observations, and theoretical calculations, and found to perform well. It can predict the surface area seen by a sensor of known FOV pointing in any direction when placed at any point in space above a specified urban surface structure. Moreover, SUM can predict the view factors of the roof, wall and ground facets seen and whether they are sunlit or shaded at any location and time of day. SUM can be used to determine the optimal placement and orientation of remote sensors to study urban radiation emissions; if the facet temperatures are known or modelled it can calculate the average temperature of the system, and it can determine the directional variation of temperature (anisotropy) due to any particular surface-sensor-sun geometric combination. The present surface geometry used in SUM is relatively simple, but there is scope to make it increasingly realistic.     

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.     

Stable Isotope Ratios: Hurricane Olivia

The oxygen and hydrogen isotopic compositions of rains from HurricaneOlivia (1994) in the eastern Pacific were measured. The rains werecollected on 24 and 25 September during airplane flights conducted at anelevation of 3 km. Hurricane Olivia peaked in intensity to a category-4storm between the two dates. Isotope ratios of rains from HurricaneOlivia were markedly lower ( ¹⁸O = –13.9to –28.8) than that of rain collected from a thunderstormat an elevation of 2.3 km outside the influence of Olivia (¹⁸O = –3.8). A distinct decrease in isotoperatios from the first day to the next ( ¹⁸O =–18.4 to –21.9) in Hurricane Olivia wasattributed to decreased updraft velocities and outflow aloft. Thisshifted the isotopic water mass balance so that fewer hydrometeors werelifted and more ice descended to flight level. A decrease in the averagedeuterium excess from the first day to the next (d = 15.5 to 7.1)was attributed to an increase in the relative humidity of the watervapor `source' area. We hypothesize that the `source' region for therain was in the boundary layer near the storm center and that becausethe hurricane was at peak intensity prior to the second day the relative humidity was higher.     

Effects of heat and water vapor transport on eddy covariance measurement of CO2 fluxes

Flux densities of carbon dioxide were measured over an arid, vegetation-free surface by eddy covariance techniques and by a heat budget-profile method, in which CO₂ concentration gradients were specified in terms of mixing ratios. This method showed negligible fluxes of CO₂, consistent with the bareness of the experimental site, whereas the eddy covariance measurements indicated large downward fluxes of CO₂. These apparently conflicting observations are in quantitative agreement with the results of a recent theory which predicts that whenever there are vertical fluxes of sensible or latent heat, a mean vertical velocity is developed. This velocity causes a mean vertical convective mass flux (= _cw for CO₂, in standard notation). The eddy covariance technique neglects this mean convective flux and measures only the turbulent flux _c w. Thus, when the net flux of CO₂ is zero, the eddy covariance method indicates an apparent flux which is equal and opposite to the mean convective flux, i.e., _c w = – _c w. Corrections for the mean convective flux are particularly significant for CO₂ because _cw and _c w are often of similar magnitude. The correct measurement of the net CO₂ flux by eddy covariance techniques requires that the fluxes of sensible and latent heat be measured as well.     

AN EVALUATION OF DISPERSION COEFFICIENTS FOR USE IN AIR QUALITY MODELS

Standard deviations of concentration in horizontal andvertical directions i.e. y andz have been estimated by using fivedifferent schemes based on empirical(due to Pasquill and Briggs)schemes and sophisticated methods(due to Irwin, Draxler, Taylor, Hanna et al.). The fiveschemes are discussed atlength. The purpose of this study is to make use ofmeteorological observations whichare routinely available, to test all the above methods andintercompare the resultswith one another and observations so that the sensitivityof each schemeunder various atmospheric stability conditions could beassessed. It has beenfound that the existing schemes are good enough to providereasonable estimates ofdispersion coefficient (y) during highly unstableconditions (Pasquill stability classes A and B). However, thesame is not true for the case when the stability increasesfrom C to F and turbulencedecreases, specifically during stable atmospheric conditions,when the observedvalues were found to be much higher than all the existingschemes. This suggests thatwhile we continue to use the current methods of estimatingthe dispersion parameters,a rigorous search is required for methods which give predictionswhich are close-to-realityduring such conditions which are represented by lowlevels (in terms of magnitude)of atmospheric turbulence leading to higher levelsof pollution.As one of the sophisticated methods requiresthe use of v and w (standard deviationsof wind velocity fluctuation in y and z directions),these have been estimated andvalidated with observed data (field experiments conductedby EPRI at Kincaid).Statistical evaluation of v and wbased on performance measures indicate a goodperformance of the parameterisations adopted in thisstudy. In the case of w duringunstable conditions a comparison of three differentschemes with observations is made.     

Modelling of turbulence over low hills and varying surface roughness

Neutral surface layer flow over low hills and varying surface roughness is considered with emphasis on closure schemes in relation to the prediction of turbulence quantities. The equations are linearised, Fourier transformed in the two horizontal directions and solved by means of a finite difference method in the vertical. Three closure schemes are. employed, namely mixing length, E- and e-- closure where E, and indicate that differential equations are used for turbulent kinetic energy, dissipation rate and shear stress. Model calculations are compared with experimental data for the step in roughness problem and for the Askervein hill. The mean flow results turn out to be relatively insensitive to the closure scheme. The shear stress and the dimensionless shear, however, are much better predicted with the E- equations than with mixing length closure. In the outer layer of the hill problem, advection of shear stress becomes important. An equation for is needed here.     

In Search of Warmer Climates? The Impact of Climate Change on Flows of British Tourists

This paper investigates the impact of climate change on the chosendestinations of Britishtourists. Destinations are characterised in terms of attractors includingclimate variables, traveland accommodation costs. These variables are used to explain the currentobserved pattern ofoverseas travel in terms of a model based upon the idea of utilitymaximisation. The approachpermits the trade-offs between climate and holiday expenditure to be analysedand effectivelyidentifies the optimal climate for generating tourism. The findings are usedto predict the impactof various climate change scenarios on popular tourist destinations.     

On Trajectory Curvature as a Selection Criterion for Valid Lagrangian Stochastic Dispersion Models

Recently Wilson and Flesch (Boundary-Layer Meteorology, 84, 411-426, 1997) suggested that the average increment d z to the orientation = arctan(w/u) of the Lagrangian velocity-fluctuation vector can be used to distinguish the better Lagrangian stochastic models within the well-mixed class. Here it is demonstrated that the specification of d z constitutes neither a sufficient or universally applicable criterion to distinguish the better Lagrangian stochastic models within the well-mixed class. The hypothesis made by Wilson and Flesch that Lagrangian stochastic models with /P_E irrotational are zero-spin models, having d z=0, is proven     

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     

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.     

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.     

Reactions of alkoxy radicals under atmospheric conditions: The relative importance of decomposition versus reaction with O2

The reactions of alkoxy radicals determine to a large extent the products formed during the atmospheric degradations of emitted organic compounds. Experimental data concerning the decompositions, 1,5-H shift isomerizations and reactions with O₂ of several classes of alkoxy radicals are inconsistent with literature estimations of their absolute or relative rate constants. An alternative, although empirical, method for assessing the relative importance under atmospheric conditions of the reactions of alkoxy radicals with O₂ versus decomposition was derived. This estimation method utilizes the differences in the heats of reaction, (H)=(H_{decomposition}–H_{O 2 reaction}), between these two reactions pathways. For (H)[22–0.5(H_{O 2 reaction})], alkoxy radical decomposition dominates over the reaction with O₂ at room temperature and atmospheric pressure of air, while for (H)[25-0.5(H_{O 2 reaction})], the O₂ reaction dominates over decomposition (where the units of H are in kcal mol^–1). The utility and shortcomings of this approach are discussed. It is concluded that further studies concerning the reactions of alkoxy radicals are needed.     

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.     

Three-dimensional modeling of synthetic cold fronts interacting with northern Alpine foehn

Summary The effect of the Alpine orography on prototype cold fronts approaching from the west is investigated by three-dimensional numerical model simulations. The numerical experiments cover a range of parameter constellations which govern the prefrontal environment of the front. Especially, the appearance and intensity of prefrontal northern Alpine foehn varies from case to case.The behaviour of a cold front north of the Alps depends much on the prefrontal condition it encounters. It is found that prefrontal foehn can either accelerate or retard the approaching front.An important feature is the pressure depression along the northern Alpine rim that results from the southerly foehn flow. In cases where this depression compensates the eastward directed pressure gradient associated with the largescale flow, the front tends to accelerate and the foehn breaks down as soon as the front passes. In contrast, the foehn prevents the front from a rapid eastward propagation if it is connected with a strong southerly wind component.No-foehn experiments are performed for comparison, where either the mountains are removed, or the static stability is set to neutral. Also shown are effects of different crossfrontal temperature contrasts.List of Symbols c _F propagation speed of a front - x, y horizontal grid spacing (cartesian system) - , horizontal grid spacing (geographic system) - t time step - z vertical grid spacing (cartesian system) - cross-frontal potential temperature difference - _i potential temperature step at an inversion - E turbulent kinetic energy - f Coriolis parameter - FGP frontogenesis parameter (see section 2.2) - g gravity acceleration (g=9.81 m s^–2) - vertical gradient of potential temperature - h terrain elevation (above MSL) - h _i height of an inversion (h _i =1000 m MSL) - H height of model lid (H=9000 m MSL) - K _M exchange coefficient of momentum - K _H exchange coefficient of heat and moisture - longitude - N Brunt-Väisäla-frequency - p pressure - Exner function (=T/) - latitude - q _v specific humidity - R _d gas constant of dry air (R _d =287.06 J kg^–1 K^–1) - density of dry air - t time - T temperature - potential temperature - TFP thermal front parameter (see section 2.2) - u, v, w cartesian wind components - u _g ,v _g geostrophic wind components - horizontal wind vector - x, y, z cartesian coordinates Abbreviations GND (above) ground level - MSL (above) mean sea level - UTC universal time coordinated With 20 Figures     

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.     

Analytical solutions for the Ekman layer

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

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