Linking urban water balance and energy balance models to analyse urban design options

Using a water balance modelling framework, this paper analyses the effects of urban design on the water balance, with a focus on evapotranspiration and storm water. First, two quite different urban water balance models are compared: Aquacycle which has been calibrated for a suburban catchment in Canberra, Australia, and the single-source urban evapotranspiration-interception scheme (SUES), an energy-based approach with a biophysically advanced representation of interception and evapotranspiration. A fair agreement between the two modelled estimates of evapotranspiration was significantly improved by allowing the vegetation cover (leaf area index, LAI) to vary seasonally, demonstrating the potential of SUES to quantify the links between water sensitive urban design and microclimates and the advantage of comparing the two modelling approaches. The comparison also revealed where improvements to SUES are needed, chiefly through improved estimates of vegetation cover dynamics as input to SUES, and more rigorous parameterization of the surface resistance equations using local-scale suburban flux measurements. Second, Aquacycle is used to identify the impact of an array of water sensitive urban design features on the water balance terms. This analysis confirms the potential to passively control urban microclimate by suburban design features that maximize evapotranspiration, such as vegetated roofs. The subsequent effects on daily maximum air temperatures are estimated using an atmospheric boundary layer budget. Potential energy savings of about 2% in summer cooling are estimated from this analysis. This is a clear ‘return on investment’ of using water to maintain urban greenspace, whether as parks distributed throughout an urban area or individual gardens or vegetated roofs.

Modelling the activated sludge flocculation process combining laser light diffraction particle sizing and population balance modelling (PBM)

A technique based on laser light diffraction is shown to be successful in collecting on-line experimental data. Time series of floc size distributions (FSD) under different shear rates (G) and calcium additions were collected. The steady state mass mean diameter decreased with increasing shear rate G and increased when calcium additions exceeded 8 mg/l. A so-called population balance model (PBM) was used to describe the experimental data, This kind of model describes both aggregation and breakage through birth and death terms. A discretised PBM was used since analytical solutions of the integro-partial differential equations are non-existing. Despite the complexity of the model, only 2 parameters need to be estimated: the aggregation rate and the breakage rate. The model seems, however, to lack flexibility. Also, the description of the floc size distribution (FSD) in time is not accurate.

Evaluation of modelled regional water balance using lake status data: a comparison of 6ka simulations with the NCAR CCM

Variations in lake area and depth reflect climatically induced changes in the water balance of overflowing as well as closed lakes. A new global data base of lake status has been assembled, and is used to compare two simulations for 6 ka (6000 yr ago) made with successive R15 versions of the NCAR Community Climate Model (CCM). Simulated water balance was expressed as anomalies of annual precipitation minus evaporation (P-E); observed water balance as anomalies of lake status. Comparisons were made visually, by comparing regional averages, and by a statistic that compares the signs of simulated P-E anomalies (smoothly interpolated to the lake sites) with the status anomalies. Both CCM0 and CCM1 showed enhanced Northern-Hemisphere monsoons at 6 ka. Both underestimated the effect, but CCM1 fitted the spatial patterns better. In the northern mid- and high-latitudes the two versions differed more, and fitted the data less satisfactorily. CCM1 performed better than CCM0 in North America and central Eurasia, but not in Europe. Both models (especially CCM0) simulated excessive aridity in interior Eurasia. The models were systematically wrong in the southern mid-latitudes. Problems may have been caused by inadequate treatment of changes in sea-surface conditions in both models. Palaeolake status data will continue to provide a benchmark for the evaluation of modelling improvements.

Population balance modelling of granulation with a physically based coalescence kernel

It was previously published by the authors that granules can either coalesce through Type I (when granules coalesce by viscous dissipation in the surface liquid layer before their surfaces touch) or Type II (when granules are slowed to a halt during rebound, after their surfaces have made contact) (AIChE J. 46 (3) (2000) 529). Based on this coalescence mechanism, a new coalescence kernel for population balance modelling of granule growth is presented. The kernel is constant such that only collisions satisfying the conditions for one of the two coalescence types are successful. One constant rate is assigned to each type of coalescence and zero is for the case of rebound. As the conditions for Types I and II coalescence are dependent on granule and binder properties, the coalescence kernel is thus physically based. Simulation results of a variety of binder and granule materials show good agreement with experimental data. (C) 2002 Elsevier Science Ltd. All rights reserved.

Nitrate leaching through climatologic water balance in a fertigated coffee plantation

Nitrate losses from soil profiles by leaching should preferentially be monitored during high rainfall events and during irrigation when fertilizer nitrogen applications are elevated. Using a climatologic water balance, based on the models of Thornthwaite and Penman Monteith for potential evapotranspiration, drainage soil water fluxes below the root zone were estimated in a fertigated coffee crop. Soil solution extraction at the depth of 1 m allowed the calculation of nitrate leaching. The average nitrate concentration in soil solution for plots that received nitrogen by fertigation at a rate of 400 kg ha-1, was 5.42 mg L-1, surpassing the limit of the Brazilian legislation of 10.0 mg L-1, only during one month. For plots receiving 800 kg ha-1 of nitrogen, the average was 25.01 mg L-1, 2.5 times higher than the above-mentioned limit. This information indicates that nitrogen rates higher than 400 kg ha-1 are potentially polluting the ground water. Yearly nitrate amounts of leaching were 24.2 and 153.0 kg ha-1 for the nitrogen rates of 400 and 800 kg ha-1, respectively. The six times higher loss indicates a cost/benefit problem for coffee fertigations above 400 kg ha-1.

Impacts of deforestation on water balance components of a watershed on the Brazilian East Coast

The Brazilian East coast was intensely affected by deforestation, which drastically cut back the original biome. The possible impacts of this process on water resources are still unknown. The purpose of this study was an evaluation of the impacts of deforestation on the main water balance components of the Galo creek watershed, in the State of Espírito Santo, on the East coast of Brazil. Considering the real conditions of the watershed, the SWAT model was calibrated with data from 1997 to 2000 and validated for the period between 2001 and 2003. The calibration and validation processes were evaluated by the Nash-Sutcliffe efficiency coefficient and by the statistical parameters (determination coefficient, slope coefficient and F test) of the regression model adjusted for estimated and measured flow data. After calibration and validation of the model, new simulations were carried out for three different land use scenarios: a scenario in compliance with the law (C1), assuming the preservation of PPAs (permanent preservation areas); an optimistic scenario (C2), which considers the watershed to be almost entirely covered by native vegetation; and a pessimistic scenario (C3), in which the watershed would be almost entirely covered by pasture. The scenarios C1, C2 and C3 represent a soil cover of native forest of 76, 97 and 0 %, respectively. The results were compared with the simulation, considering the real scenario (C0) with 54 % forest cover. The Nash-Sutcliffe coefficients were 0.65 and 0.70 for calibration and validation, respectively, indicating satisfactory results in the flow simulation. A mean reduction of 10 % of the native forest cover would cause a mean annual increase of approximately 11.5 mm in total runoff at the watershed outlet. Reforestation would ensure minimum flows in the dry period and regulate the maximum flow of the main watercourse of the watershed.

Water Balance Components in Covered and Uncovered Soil Growing Irrigated Muskmelon

ABSTRACT Knowledge of the terms (or processes) of the soil water balance equation or simply the components of the soil water balance over the cycle of an agricultural crop is essential for soil and water management. Thus, the aim of this study was to analyze these components in a Cambissolo Háplico (Haplocambids) growing muskmelon (Cucumis melo L.) under drip irrigation, with covered and uncovered soil, in the municipality of Baraúna, State of Rio Grande do Norte, Brazil (05º 04’ 48” S, 37º 37’ 00” W). Muskmelon, variety AF-646, was cultivated in a flat experimental area (20 × 50 m). The crop was spaced at 2.00 m between rows and 0.35 m between plants, in a total of ten 50-m-long plant rows. At points corresponding to ⅓ and ⅔ of each plant row, four tensiometers (at a distance of 0.1 m from each other) were set up at the depths of 0.1, 0.2, 0.3, and 0.4 m, adjacent to the irrigation line (0.1 m from the plant row), between two selected plants. Five random plant rows were mulched using dry leaves of banana (Musa sp.) along the drip line, forming a 0.5-m-wide strip, which covered an area of 25 m2 per of plant row with covered soil. In the other five rows, there was no covering. Thus, the experiment consisted of two treatments, with 10 replicates, in four phenological stages: initial (7-22 DAS - days after sowing), growing (22-40 DAS), fruiting (40-58 DAS) and maturation (58-70 DAS). Rainfall was measured with a rain gauge and water storage was estimated by the trapezoidal method, based on tensiometer readings and soil water retention curves. For soil water flux densities at 0.3 m, the tensiometers at the depths of 0.2, 0.3, and 0.4 m were considered; the tensiometer at 0.3 m was used to estimate soil water content from the soil water retention curve at this depth, and the other two to calculate the total potential gradient. Flux densities were calculated through use of the Darcy-Buckingham equation, with hydraulic conductivity determined by the instantaneous profile method. Crop actual evapotranspiration was calculated as the unknown of the soil water balance equation. The soil water balance method is effective in estimating the actual evapotranspiration of irrigated muskmelon; there was no significant effect of soil coverage on capillary rise, internal drainage, crop actual evapotranspiration, and muskmelon yield compared with the uncovered soil; the transport of water caused by evaporation in the uncovered soil was controlled by the break in capillarity at the soil-atmosphere interface, which caused similar water dynamics for both management practices applied.

Estimate of Groundwater Recharge Based on Water Balance in The Unsaturated Soil Zone

ABSTRACT Groundwater management depends on the knowledge on recharge rates and water fluxes within aquifers. The recharge is one of the water cycle components most difficult to estimate. As a result, despite the chosen method, the estimates are subject to uncertainties that can be identified by means of comparison with other approaches. In this study, groundwater recharge estimates based on the water balance in the unsaturated zone is assessed. Firstly, the approach is evaluated by comparing the results with those of another method. Then, the estimates are used as inputs in a transient groundwater flow model in order to assess how the water table would respond to the obtained recharges rates compared to measured levels. The results suggest a good performance of the adopted approach and, despite some inherent limitations, it has advantages over other methods since the data required are easier to obtain.

Differential growth responses to water balance of coexisting deciduous tree species are linked to wood density in a Bolivian tropical dry forest

A seasonal period of water deficit characterizes tropical dry forests (TDFs). There, sympatric tree species exhibit a diversity of growth rates, functional traits, and responses to drought, suggesting that each species may possess different strategies to grow under different conditions of water availability. The evaluation of the long-term growth responses to changes in the soil water balance should provide an understanding of how and when coexisting tree species respond to water deficit in TDFs. Furthermore, such differential growth responses may be linked to functional traits related to water storage and conductance. We used dendrochronology and climate data to retrospectively assess how the radial growth of seven coexisting deciduous tree species responded to the seasonal soil water balance in a Bolivian TDF. Linear mixed-effects models were used to quantify the relationships between basal area increment and seasonal water balance. We related these relationships with wood density and sapwood production to assess if they affect the growth responses to climate. The growth of all species responded positively to water balance during the wet season, but such responses differed among species as a function of their wood density. For instance, species with a strong growth response to water availability averaged a low wood density which may facilitate the storage of water in the stem. By contrast, species with very dense wood were those whose growth was less sensitive to water availability. Coexisting tree species thus show differential growth responses to changes in soil water balance during the wet season. Our findings also provide a link between wood density, a trait related to the ability of trees to store water in the stem, and wood formation in response to water availability.

Weather variables, water balance, growth, and agro industrial yield of sugarcane

The objective of this study was to evaluate the meteorological variables, water deficiency, growth, and agro-industrial yield of sugarcane varieties: RB72454, RB863129, RB867515, RB92579, RB93509, RB931003, RB951541, and RB971755, in rainfed crop in two harvests in the Rio Largo-AL region. The meteorological variables were obtained in an automatic station and water balance was done by Thornthwaite & Mather method. During the study period, the air temperature ranged from 16.6 to 35.9 ºC. In the first production cycle rained 1,806 mm and the crop evapotranspiration was 1,775 mm. In the second cycle, the rainfall totaled 1,632 mm and the crop evapotranspiration was 1,290 mm. The average water excess of two production cycles was 689 mm and the water deficit totaled 665 mm. The average agricultural productivity in the plant was 86.8 t ha-1, in the first ratoon was 75.2 t ha-1 and the agro-industrial yield average was 12.9 and 10.9 tons of sugar per hectare in the plant and first ratoon, respectively. The air temperature was not limiting to the growth of sugarcane and the rainfall was higher than the crop evapotranspiration, but due to poor distribution of the rains there was water deficit. The most productive varieties were RB93509, RB92579, and RB863129.

Application of artificial neural networks as an alternative to volumetric water balance in drip irrigation management in watermelon crop

Precision irrigation seeks to establish strategies which achieve an efficient ratio between the volume of water used (reduction in input) and the productivity obtained (increase in production). There are several studies in the literature on strategies for achieving this efficiency, such as those dealing with the method of volumetric water balance (VWB). However, it is also of great practical and economic interest to set up versatile implementations of irrigation strategies that: (i) maintain the performance obtained with other implementations, (ii) rely on few computational resources, (iii) adapt well to field conditions, and (iv) allow easy modification of the irrigation strategy. In this study, such characteristics are achieved when using an Artificial Neural Network (ANN) to determine the period of irrigation for a watermelon crop in the Irrigation Perimeter of the Lower Acaraú, in the state of Ceará, Brazil. The Volumetric Water Balance was taken as the standard for comparing the management carried out with the proposed implementation of ANN. The statistical analysis demonstrates the effectiveness of the proposed management, which is able to replace VWB as a strategy in automation.

Multimodel estimate of the global terrestrial water balance: setup and first results

Six land surface models and five global hydrological models participate in a model intercomparison project (WaterMIP), which for the first time compares simulation results of these different classes of models in a consistent way. In this paper the simulation setup is described and aspects of the multi-model global terrestrial water balance are presented. All models were run at 0.5 degree spatial resolution for the global land areas for a 15-year period (1985-1999) using a newly-developed global meteorological dataset. Simulated global terrestrial evapotranspiration, excluding Greenland and Antarctica, ranges from 415 to 586 mm year-1 (60,000 to 85,000 km3 year-1) and simulated runoff ranges from 290 to 457 mm year-1 (42,000 to 66,000 km3 year-1). Both the mean and median runoff fractions for the land surface models are lower than those of the global hydrological models, although the range is wider. Significant simulation differences between land surface and global hydrological models are found to be caused by the snow scheme employed. The physically-based energy balance approach used by land surface models generally results in lower snow water equivalent values than the conceptual degree-day approach used by global hydrological models. Some differences in simulated runoff and evapotranspiration are explained by model parameterizations, although the processes included and parameterizations used are not distinct to either land surface models or global hydrological models. The results show that differences between model are major sources of uncertainty. Climate change impact studies thus need to use not only multiple climate models, but also some other measure of uncertainty, (e.g. multiple impact models).

The surface urban energy and water balance scheme (SUEWS): evaluation in Los Angeles and Vancouver

An urban energy and water balance model is presented which uses a small number of commonly measured meteorological variables and information about the surface cover. Rates of evaporation-interception for a single layer with multiple surface types (paved, buildings, coniferous trees and/or shrubs, deciduous trees and/or shrubs, irrigated grass, non-irrigated grass and water) are calculated. Below each surface type, except water, there is a single soil layer. At each time step the moisture state of each surface is calculated. Horizontal water movements at the surface and in the soil are incorporated. Particular attention is given to the surface conductance used to model evaporation and its parameters. The model is tested against direct flux measurements carried out over a number of years in Vancouver, Canada and Los Angeles, USA. At all measurement sites the model is able to simulate the net all-wave radiation and turbulent sensible and latent heat well (RMSE = 25–47 W m−2, 30–64 and 20–56 W m−2, respectively). The model reproduces the diurnal cycle of the turbulent fluxes but typically underestimates latent heat flux and overestimates sensible heat flux in the day time. The model tracks measured surface wetness and simulates the variations in soil moisture content. It is able to respond correctly to short-term events as well as annual changes. The largest uncertainty relates to the determination of surface conductance. The model has the potential be used for multiple applications; for example, to predict effects of regulation on urban water use, landscaping and planning scenarios, or to assess climate mitigation strategies.

The water balance of northern Africa during the mid-Holocene: an evaluation of the 6ka BP PMIP simulations

Runoff fields over northern Africa (10–25°N, 20°W–30°E) derived from 17 atmospheric general circulation models driven by identical 6 ka BP orbital forcing, sea surface temperatures, and CO2 concentration have been analyzed using a hydrological routing scheme (HYDRA) to simulate changes in lake area. The AGCM-simulated runoff produced six-fold differences in simulated lake area between models, although even the largest simulated changes considerably underestimate the observed changes in lake area during the mid-Holocene. The inter-model differences in simulated lake area are largely due to differences in simulated runoff (the squared correlation coefficient, R2, is 0.84). Most of these differences can be attributed to differences in the simulated precipitation (R2=0.83). The higher correlation between runoff and simulated lake area (R2=0.92) implies that simulated differences in evaporation have a contributory effect. When runoff is calculated using an offline land-surface scheme (BIOME3), the correlation between runoff and simulated lake area is (R2=0.94). Finally, the spatial distribution of simulated precipitation can exert an important control on the overall response.