Estimation of possible climate change impacts on water availability, extreme flow events and soil moisture in the Goulburn and Ovens Basins, Victoria

A rainfall-runoff model (IHACRES) is applied on a daily timestep to a large area of the state of Victoria, Australia. Successful calibrations of this dynamic lumped parameter model were performed for 5 rivers contributing streamflow to the Ovens Basin, and for 9 rivers of the Goulburn Basin. This is the first application of the model on such a scale, involving two basins where the total drainage area of the catchments modelled is about 6,500 km². The models were tested by simulation over the entire common period of observation for the 14 catchments under consideration. The results show that the models closely simulate the observed streamflow.The effect of historical climate variability on streamflow was investigated. The models were used for estimation of the potential impact of climatic change on water availability for irrigation for different climate scenarios developed in the Division of Atmospheric Research, CSIRO. This allows conditional estimates to be made of water supply in these basins for the periods 2030 and 2070 under current vegetation conditions. Projecting the future hydrologic regime in this region is extremely important, in particular for supporting irrigation management of the Basin.The problem of estimating the impact of climate change on the probability of extreme events of the hydrological regime was analysed. Flood frequency was found to increase for the scenarios providing the maximum amount of water; to 50% at 2030 and 100% at 2070. The probability of flood events for the dry scenarios rapidly decreases for these dates. Drought frequency, as defined by a soil wetness index, increased 35% for the dry scenario at 2030 and 80% for this scenario at 2070.     

REGIONAL DEPENDENCY OF PRECIPITATION-ALTITUDE RELATIONSHIP IN THE SWISS ALPS

The spatial distribution of the regional long-term rate of precipitation change, P, with altitude in Switzerland is analysed. In order to eliminate the bias of point precipitation measurement, the analysis is based on precipitation sums corrected for the systematic error of precipitation measurement, particularly the wind-induced error and wetting losses. The P values indicate a large spatial heterogeneity. They show different local and regional patterns. In this respect, the effect of the main alpine ridge dividing the Swiss territory into northern and southern parts is crucial. Water budget computations indicate that the regional P values have a tendency to overestimate areal precipitation. The possible reasons are the redistribution of precipitation by wind and its accumulation in the valleys and, at the same time, the biased precipitation networks preferring location of gauges also in the valleys.     

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.     

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.     

Transport-Dissipation Analytical Solutions to the E-∈Turbulence Model and their Role in Predictions of the Neutral ABL

E- turbulence model predictions of the neutralatmospheric boundary layer (NABL) are reinvestigated to determine thecause for turbulence overpredictions found in previous applications. Analytical solutions to the coupled E and equations for the case of steady balance between transport and dissipation terms, the dominant balance just below the NABL top, are derived. It is found that analytical turbulence profiles laminarizeat a finite height only for values of closure parameter ratio c ₂ /_e equal toor slightly greater than one, with laminarization as z for greater . The point = 2 is additionally foundthat where analytical turbulent length scale (l) profilesmade a transition from ones ofdecreasing ( < 2) to increasing ( > 2)values with height. Numerically predicted profiles near the NABL topare consistent with analytical findings. The height-increasingvalues of l predicted throughout the NABL with standard values ofclosure parameters thus appear a consequence of 2.5(> 2), implied by these values (c ₂ = 1.92, _{= 1.3}, _{e = 1}). Comparison of numericalpredictions with DNS data shows that turbulence overpredictions obtained with standard-valued parameters are rectifiedby resetting and _e to 1.1 and 1.6, respectively, giving, with c ₂ = 1.92, 1.3, and laminarization of the NABL's cappingtransport-dissipation region at a finite height.     

The impact of climate change on the river rhine: A scenario study

This paper concerns the impact of human-induced global climate change on the River Rhine discharge. For this purpose a model for climate assessment, named ESCAPE, is coupled to a water balance model, named RHINEFLOW. From climate scenarios, changes in regional annual water availability and seasonal discharge in the River Rhine Basin are estimated. The climate scenarios are based on greenhouse gases emissions scenarios. An assessment is made for best guess seasonal discharge changes and for changes in frequencies of low and high discharges in the downstream reaches of the river. In addition, a quantitative estimation of the uncertainties associated with this guess is arrived at.The results show that the extent and range of uncertainty is large with respect to the best guess changes. The uncertainty range is 2–3 times larger for the Business-as-Usual than for the Accelerated Policies scenarios. This large range stems from the doubtful precipitation simulations from the present General Circulation Models. This scenario study showed the precipitation scenarios to be the key-elements within the present range of reliable climate change scenarios.For the River Rhine best guess changes for annual water availability are small according to both scenarios. The river changes from a present combined snow-melt-rain fed river to an almost entirely rain fed river. The difference between present-day large average discharge in winter and the small average discharge in autumn should increase for all scenarios. This trend is largest in the Alpine part of the basin. Here, winter discharges should increase even for scenarios forecasting annual precipitation decreases. Summer discharge should decrease. Best guess scenarios should lead to increased frequencies of both low and high flow events in the downstream (Dutch) part of the river. The results indicate changes could be larger than presently assumed in worst case scenarios used by the Dutch water management authorities.     

On the representation of frequency spectra in meteorology

The common representation of frequency spectra in meteorology and climatology is discussed. It is pointed out that this representation is misleading since spectral peaks and spectral gaps are obtained even when the spectrum density is monotonously decreasing in the whole frequency range. A plea is made for using the spectrum distribution function, F() which gives an unambiguous picture of the distribution of variance with frequency.     

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.     

On the inner-layer scale height of boundary-layer flow over low hills

Jackson and Hunt's (1975) equation for the depth of the inner layer of flows over low hills does not depend on any closure assumption as contrarily supposed in literature. This equation contains a constant which can arbitrarily be specified. It is suggested that this inner-layer constant should be determined from experimental data. A preliminary check with some data from the Askervein experiment suggests that Jackson and Hunt's equation fits these data almost as well as Jensen's equation provided that fitted inner-layer constants are used.     

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.     

The effect of synoptic-scale advection on the performance of the Priestley-Taylor evaporation formula

Measurements of the appropriate parameters for the calculation of the latent heat flux over a black spruce forest in northern Quebec were carried out in August, 1980. Values of the Priestley-Taylor parameter, , were derived by exploiting the Bowen-ratio-energy-balance (BREB) technique. Derived values of are then related to synoptic-scale warm air advection, derived from surface synoptic charts and tephigrams of the planetary boundary layer. It is found that when warm air advection is present, > 1.26, especially when the surface is wet. When advection enhancement is removed, however, values of approach unity. A new approach to calculating the latent heat flux, when warm air advection is present, is therefore proposed.     

Averaging Intervals For Spectral Analysis Of Nonstationary Turbulence

We formulate a method for determining the smallest time interval Tover which a turbulence time series can be averaged to decompose it intoinstantaneous mean and random components. From the random part the method defines the optimal interval (or averaging window) AW over which this part should be averaged to obtain the instantaneous spectrum. Both T and AW vary randomly with time and depend on physical properties of the turbulence. T also depends on the accuracy of the measurements and is thus independent of AW. Interesting features of the method are its real-time capability and the non-equality between AW and T.     

Vegetation responses and feedbacks to climate: a review of models and processes

Vegetation changes both in stationary and changing climates. Such changes can significantly affect hydrological and climate dynamics. Probabilistic, inferential, empirical, statistical, threshold, ecophysiological, and mechanistic vegetation models provide tools and ideas to construct coupled climate and vegetation schemes to study climate/vegetation feedbacks. Their logic is discussed, typical applications are presented, and their usefulness is assessed. Developing coupled climate and vegetation schemes also implies tackling scaling issues explicitly. Just as the Courant-Friedrichs-Lewy (CFL) criterion guarantees that information is not transferred faster through space than time in climate models, information should be transmitted fast enough in vegetation models for the landscape to register vegetation responses. To guarantee that this is the case, a migration criterion, or m criterion, is proposed. The CFL criterion and the m criterion set formal constraints on the design of coupled atmosphere and vegetation schemes. In particular, the ratio of climate and vegetation space scales should be approximately five orders of magnitude less than the ratio of climate and vegetation time scales.     

Simulating the causal elements of urban heat islands

A comparison was made between the resultant surface temperatures and sensible heat fluxes of building interfaces calculated by steady-state and transient (implicit) methods. Both procedures used identical environmental (summer and winter) input. For exterior conditions, the results indicated that the correlation between the two methods is sufficiently large, enabling them to be used interchangeably for the spatial analysis of urban canopy layers of entire cities. Using a steady-state approach as a surrogate for unsteady conditions, computer resources can be saved up to a factor of ten. An urban energy budget model (URBAN 3) has been used to demonstrate that the distribution of sensible heat flux and net longwave radiation — the prime causes of urban heat-island generation — was far from the homogeneity assumed in many macroscopic models or even some street-level studies. The individual emanations of reradiation and sensible heat flux showed different diurnal and spatial patterns. Under the input scenarios used, daytime heat islands assumed a doughnut shape in the inner city. It is believed that many of the lower boundary conditions used in macroscopic numerical models are inadequate in light of this study.     

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.     

Near-field dispersion from instantaneous sources in the surface layer

This paper considers the near-field dispersion of an ensemble of tracer particles released instantaneously from an elevated source into an adiabatic surface layer. By modelling the Lagrangian vertical velocity as a Markov process which obeys the Langevin equation, we show analytically that the mean vertical drift velocity w(t) is w()=bu _*(1–e ^–(1+)), where is time since release (nondimensionalized with the Lagrangian time scale at the source), b Batchelor's constant, and u _*, the friction velocity. Hence, the mean height and mean depth of the ensemble are calculated. Although the derivation is formally valid only when 1, the predictions for w, mean height and mean depth are consistent in the downstream limit ( 1) with surface-layer Lagrangian similarity theory and with the diffusion equation. By comparing the analytical predictions with numerical, randomflight solutions of the Langevin equation, the analytical predictions are shown to be good approximations at all times, both near-field and far-field.     

Small-Scale Variability in the Coastal Atmospheric Boundary Layer

The influence of the main large-scale wind directions on thermally driven mesoscale circulations at the Baltic southwest coast, southeast of Sweden, is examined. The aim of the study is to highlight small-scale alterations in the coastal atmospheric boundary layer. A numerical three-dimensional mesoscale model is used in this study, which is focused on an overall behaviour of the coastal jets, drainage flows, sea breezes, and a low-level eddy-type flow in particular. It is shown that synoptic conditions, together with the moderate terrain of the southeast of Sweden (max. height h₀ 206 m), governs the coastal mesoscale dynamics triggered by the land-sea temperature difference T. The subtle nature of coastal low-level jets and sea breezes is revealed; their patterns are dictated by the interplay between synoptic airflow, coastline orientation, and T.The simulations show that coastal jets typically occur during nighttime and vary in height, intensity and position with respect to the coast; they interact with downslope flows and the background wind. For the assigned land surface temperature (varying ±8 K from the sea temperature) and the opposing constant geostrophic wind 8 m s^-1, the drainage flow is more robust to the opposing ambient flow than the sea breeze later on. Depending on the part of the coast under consideration, and the prevailing ambient wind, the sea breeze can be suppressed or enhanced, stationary at the coast or rapidly penetrating inland, locked up in phase with another dynamic system or almost independently self-evolving. A low-level eddy structure is analyzed. It is governed by tilting, divergence and horizontal advection terms. The horizontal extent of the coastal effects agrees roughly with the Rossby radius of deformation.     

Triggering of atmospheric circulations by moisture inhomogeneities of the earth's surface

A two-dimensional mesoscale soil-atmosphere model is used to simulate the triggering of atmospheric convection by horizontally varying soil water content. The variation is periodic with a wavelength between 4 and 40 km, which is considered a realistic scale for the variation of land surface characteristics. Three stages of convection can be clearly discerned: a short initial stage when convection sets in and where the size of the conective cells is determined by , a mature stage with well developed cells whose size is still determined by , and a decay/transformation stage, characterized by the formation of narrow regions of strong updrafts and wide regions of moderate downdrafts, independent of . Parameters relevant for the transition are given, and the importance of the feedback between soil and atmosphere is demonstrated. The dependence of convective parameters, e.g., height of the convective layer, vertical velocity and fluxes of heat and moisture on is investigated. The calculations of the mature stage are compared with the predictions of a linear model.     

Stable boundary layers: Observations,models and variability part I: Modelling and measurements

Selected field measurements of evening stable boundary layers are presented in detail comparable with published Large Eddy Simulation results. Such models appear to match idealized theories more closely than do some boundary-layer observations. Any attempt to compare detailed observations with idealized models therefore highlights the variability of the real boundary layer.Here direct turbulence measurements across the stable boundary layer from a heterogeneous and an ideal site are contrasted. Recommendations are made for the information needed to distinguish heterogeneous and ideal cases.The companion paper (Part II) discusses further the issues of data, analysis in the presence of variability, and the effects of averaging over heterogeneous terrain.Part of UK Meteorological Office Atmospheric Process Research Division.     

Increasing cloud in a warming world

Cloud amount records for the U.S.A. have been analyzed in the context of the warming world analogue model described by Lough et al. (1983). Cloud amount increases over practically the entire U.S.A. in all seasons. This result considerably strengthens the more tentative conclusion of Henderson-Sellers (1986) that cloud amount increases over Europe in the same warming world scenario. These results are in contrast to the few numerical model predictions of cloud changes in warming world experiments. A possible, rather tantalizing, conclusion is that current GCM cloud prediction schemes tend to enhance temperature increases through cloud-climate feedback whereas the historical data could suggest a negative feedback. Part, possibly all, of this difference may be the result of the fundamental distinction between the two experimental scenarios: the equilibrium change modelled by GCMs as compared to the smaller transient change represented by the historical analogue. On the other hand the current real-world experiment is, itself, a transient change in boundary and atmospheric conditions. At the least, surface-observed cloudiness seems to offer a useful and complementary data source with which to examine one aspect of the performance of numerical climate models.