Water Cycle Model Research Articles

Evapotranspiration and water yield of a pine-broadleaf forest are not altered by long-term atmospheric [CO2 ] enrichment under native or enhanced soil fertility.

Changes in evapotranspiration (ET) from terrestrial ecosystems affect their water yield (WY), with considerable ecological and economic consequences. Increases in surface runoff observed over the past century have been attributed to increasing atmospheric CO2 concentrations resulting in reduced ET by terrestrial ecosystems. Here, we evaluate the water balance of a Pinus taeda (L.) forest with a broadleaf component that was exposed to atmospheric [CO2 ] enrichment (ECO2 ; +200ppm) for over 17years and fertilization for 6years, monitored with hundreds of environmental and sap flux sensors on a half-hourly basis. These measurements were synthesized using a one-dimensional Richard's equation model to evaluate treatment differences in transpiration (T), evaporation (E), ET, and WY. We found that ECO2 did not create significant differences in stand T, ET, or WY under either native or enhanced soil fertility, despite a 20% and 13% increase in leaf area index, respectively. While T, ET, and WY responded to fertilization, this response was weak (<3% of mean annual precipitation). Likewise, while E responded to ECO2 in the first 7years of the study, this effect was of negligible magnitude (<1% mean annual precipitation). Given the global range of conifers similar to P. taeda, our results imply that recent observations of increased global streamflow cannot be attributed to decreases in ET across all ecosystems, demonstrating a great need for model-data synthesis activities to incorporate our current understanding of terrestrial vegetation in global water cycle models.

Full Coupling Between the Atmosphere, Surface, and Subsurface for Integrated Hydrologic Simulation

AbstractAn ever increasing community of earth system modelers is incorporating new physical processes into numerical models. This trend is facilitated by advancements in computational resources, improvements in simulation skill, and the desire to build numerical simulators that represent the water cycle with greater fidelity. In this quest to develop a state‐of‐the‐art water cycle model, we coupled HydroGeoSphere (HGS), a 3‐D control‐volume finite element surface and variably saturated subsurface flow model that includes evapotranspiration processes, to the Weather Research and Forecasting (WRF) Model, a 3‐D finite difference nonhydrostatic mesoscale atmospheric model. The two‐way coupled model, referred to as HGS‐WRF, exchanges the actual evapotranspiration fluxes and soil saturations calculated by HGS to WRF; conversely, the potential evapotranspiration and precipitation fluxes from WRF are passed to HGS. The flexible HGS‐WRF coupling method allows for unique meshes used by each model, while maintaining mass and energy conservation between the domains. Furthermore, the HGS‐WRF coupling implements a subtime stepping algorithm to minimize computational expense. As a demonstration of HGS‐WRF's capabilities, we applied it to the California Basin and found a strong connection between the depth to the groundwater table and the latent heat fluxes across the land surface.

Optimal Allocation Method of Irrigation Water from River and Lake by Considering the Field Water Cycle Process

At present, the shortage of water resources has become a serious constraint to the further development of social economy. The optimal allocation of multi-water resources is valuable for agricultural water management in arid and semi-arid regions. However, traditional deterministic programming does not solve the complex water resources allocation in irrigation systems. Furthermore, previous allocation methods of irrigation water seldom considered the water cycle process, especially for multi-sources of irrigation. In this study, we develop an inexact fuzzy stochastic simulation-optimization programming (IFSSOP) model for the irrigation water optimal allocation of two water sources. The model combines the crop water model and the field water cycle model with an uncertainty optimization model, which considers the contribution of groundwater to crop water consumption. As a case study, the developed model is used in an arid area with two irrigation water sources: a river and a lake. Accordingly, the total optimal allocation irrigation water amounts of river and lake water under different violation probabilities in various hydrological years were obtained. By comparing the IFSSOP model with the IFSSOP model without considering the contribution of shallow groundwater (IFSSOP-NG model), it can be shown that the system benefits of the developed model are higher. With the lake water source from flood water, the region can save 30–34% of the river water, maintaining the original crop water deficit irrigation ratio. Consequently, application of the IFSSOP model in irrigation scheduling will provide effective water allocation patterns to save more water in an arid region with shallow groundwater.

Validation of terrestrial water storage variations as simulated by different global numerical models with GRACE satellite observations

Abstract. Estimates of terrestrial water storage (TWS) variations from the Gravity Recovery and Climate Experiment (GRACE) satellite mission are used to assess the accuracy of four global numerical model realizations that simulate the continental branch of the global water cycle. Based on four different validation metrics, we demonstrate that for the 31 largest discharge basins worldwide all model runs agree with the observations to a very limited degree only, together with large spreads among the models themselves. Since we apply a common atmospheric forcing data set to all hydrological models considered, we conclude that those discrepancies are not entirely related to uncertainties in meteorologic input, but instead to the model structure and parametrization, and in particular to the representation of individual storage components with different spatial characteristics in each of the models. TWS as monitored by the GRACE mission is therefore a valuable validation data set for global numerical simulations of the terrestrial water storage since it is sensitive to very different model physics in individual basins, which offers helpful insight to modellers for the future improvement of large-scale numerical models of the global terrestrial water cycle.

Incorporation of a dynamic root distribution into CLM4.5: Evaluation of carbon and water fluxes over the Amazon

Roots are responsible for the uptake of water and nutrients by plants and have the plasticity to dynamically respond to different environmental conditions. However, most land surface models currently prescribe rooting profiles as a function only of vegetation type, with no consideration of the surroundings. In this study, a dynamic rooting scheme, which describes root growth as a compromise between water and nitrogen availability, was incorporated into CLM4.5 with carbon–nitrogen (CN) interactions (CLM4.5-CN) to investigate the effects of a dynamic root distribution on eco-hydrological modeling. Two paired numerical simulations were conducted for the Tapajos National Forest km83 (BRSa3) site and the Amazon, one using CLM4.5-CN without the dynamic rooting scheme and the other including the proposed scheme. Simulations for the BRSa3 site showed that inclusion of the dynamic rooting scheme increased the amplitudes and peak values of diurnal gross primary production (GPP) and latent heat flux (LE) for the dry season, and improved the carbon (C) and water cycle modeling by reducing the RMSE of GPP by 0.4 g C m-2 d-1, net ecosystem exchange by 1.96 g C m-2 d-1, LE by 5.0 W m-2, and soil moisture by 0.03 m3 m-3, at the seasonal scale, compared with eddy flux measurements, while having little impact during the wet season. For the Amazon, regional analysis also revealed that vegetation responses (including GPP and LE) to seasonal drought and the severe drought of 2005 were better captured with the dynamic rooting scheme incorporated.

Effect of Catchment-Scale Green Roof Deployment on Stormwater Generation and Reuse in a Tropical City

Low-impact development (LID) comprises a broad spectrum of stormwater management technologies for mitigating the impacts of urbanization on hydrological processes. Among these technologies, green roofs are one of the most adopted solutions, especially in densely populated metropolitan areas, where roofs take up a significant portion of the impervious surfaces and land areas are scarce. While the in situ hydrological performance of green roofs—i.e., reduction of runoff volume and peak discharge—is well addressed in literature, less is known about their impact on stormwater management and reuse activities at a catchment or city scale. This study developed an integrated urban water cycle model (IUWCM) to quantitatively assess the effect of uniform green roof deployment (i.e., 25, 50, and 100% conversion of traditional roofs) over the period 2009–2011 in the Marina Reservoir catchment, a 100-km2, highly urbanized area located in the heart of Singapore. The IUWCM consists of two components: (1) a physically based model for extensive green roofs integrated within a one-dimensional numerical hydrological-hydraulic catchment model linked with (2) an optimization-based model describing the operation of the downstream, stormwater-fed reservoir. The event-based hydrological performance of green roofs varied significantly throughout the simulation period with a median of about 5% and 12% for the catchment scale reduction of runoff volume and peak discharge (100% conversion of traditional roofs). The high variability and lower performance (with respect to temperate climates) are strongly related to the tropical weather and climatic conditions—e.g., antecedent dry weather period and maximum rainfall intensity. Average annual volume reductions were 0.6, 1.2, and 2.4% for the 25, 50, and 100% green roof scenarios, respectively. The reduction of the stormwater generated at the catchment level through green roof implementation had a positive impact on flood protection along Marina Reservoir shores and the energy costs encountered when operating the reservoir. Vice versa, the drinking water supply, which depends on the amount of available stormwater, decreased due to the evapotranspiration losses from green roofs. Better performance in terms of stormwater reuse could only be obtained by increasing the time of concentration of the catchment. This may be achieved through the combination of green roofs with other LID structures.

Leaf transpiration of drought tolerant plant can be captured by hyperspectral reflectance using PLSR analysis

Abstract: A clear understanding of plant transpiration is a crucial step for water cycle and climate modeling, especially for arid ecosystems in which water is one of the major constraints. Traditional field measurements of leaf scale transpiration are always time-consuming and often unfeasible in the context of large spatial and temporal scales. This study focused on a dominant native plant in the arid land of central Asia, Haloxylon ammondendron, with the aim of deriving the leaf-scale transpiration through hyperspectral reflectance using Partial Least Squares Regression (PLSR) analysis. The results revealed that the PLSR model based on the first-order derivative spectra at wavelengths selected through stepwise regression analysis can closely trace leaf transpiration with a high accuracy (R2 = 0.78, RMSE = 1.62 µmol g-1 s-1). The accuracy is also relatively stable even at a spectral resolution of 10 nm, which is very close to the bandwidths of several running satellite-borne hyperspectral sensors such as Hyperion. The results also proved that the first-order derivative spectra within the shortwave infrared (SWIR) domain, especially at 2435, 2440, 2445, and 2470 nm, were critical for PLSR models to predict leaf transpiration. These findings highlight a promising strategy for developing remote sensing methods to potentially characterize transpiration at broad scales.

Coupling carbon allocation with leaf and root phenology predicts tree–grass partitioning along a savanna rainfall gradient

Abstract. The relative complexity of the mechanisms underlying savanna ecosystem dynamics, in comparison to other biomes such as temperate and tropical forests, challenges the representation of such dynamics in ecosystem and Earth system models. A realistic representation of processes governing carbon allocation and phenology for the two defining elements of savanna vegetation (namely trees and grasses) may be a key to understanding variations in tree–grass partitioning in time and space across the savanna biome worldwide. Here we present a new approach for modelling coupled phenology and carbon allocation, applied to competing tree and grass plant functional types. The approach accounts for a temporal shift between assimilation and growth, mediated by a labile carbohydrate store. This is combined with a method to maximize long-term net primary production (NPP) by optimally partitioning plant growth between fine roots and (leaves + stem). The computational efficiency of the analytic method used here allows it to be uniquely and readily applied at regional scale, as required, for example, within the framework of a global biogeochemical model.We demonstrate the approach by encoding it in a new simple carbon–water cycle model that we call HAVANA (Hydrology and Vegetation-dynamics Algorithm for Northern Australia), coupled to the existing POP (Population Orders Physiology) model for tree demography and disturbance-mediated heterogeneity. HAVANA-POP is calibrated using monthly remotely sensed fraction of absorbed photosynthetically active radiation (fPAR) and eddy-covariance-based estimates of carbon and water fluxes at five tower sites along the North Australian Tropical Transect (NATT), which is characterized by large gradients in rainfall and wildfire disturbance. The calibrated model replicates observed gradients of fPAR, tree leaf area index, basal area, and foliage projective cover along the NATT. The model behaviour emerges from complex feedbacks between the plant physiology and vegetation dynamics, mediated by shifting above- versus below-ground resources, and not from imposed hypotheses about the controls on tree–grass co-existence. Results support the hypothesis that resource limitation is a stronger determinant of tree cover than disturbance in Australian savannas.

Water and the city: exploring links between urban growth and water demand management

Urban water management is currently understood as a socio-technical problem, including both technologies and engineering interventions as well as socioeconomic dimensions and contexts vis-à-vis both end users and institutions. In this framework, perhaps the most important driver of urban water demand, at the intersection between engineering, social and economic domains, is urban growth. This paper examines aspects of the interplay between the dynamics of urban growth and the urban water cycle. Specifically, a cellular automata urban growth model is re-engineered to provide growth patterns at the level of detail needed by an urban water cycle model. The resulting toolkit is able to simulate spatial changes in urban areas while simultaneously estimating their water demand impact under different water demand management scenarios, with an emphasis on distributed technologies whose applicability depends on urban form. The method and tools are tested in the case study of Mesogeia, Greece, and conclusions are drawn, regarding both the performance of the urban growth model and the effectiveness of different urban water management practices.

The application of a Web-geographic information system for improving urban water cycle modelling.

Research in urban water management has experienced a transition from traditional model applications to modelling water cycles as an integrated part of urban areas. This includes the interlinking of models of many research areas (e.g. urban development, socio-economy, urban water management). The integration and simulation is realized in newly developed frameworks (e.g. DynaMind and OpenMI) and often assumes a high knowledge in programming. This work presents a Web based urban water management modelling platform which simplifies the setup and usage of complex integrated models. The platform is demonstrated with a small application example on a case study within the Alpine region. The used model is a DynaMind model benchmarking the impact of newly connected catchments on the flooding behaviour of an existing combined sewer system. As a result the workflow of the user within a Web browser is demonstrated and benchmark results are shown. The presented platform hides implementation specific aspects behind Web services based technologies such that the user can focus on his main aim, which is urban water management modelling and benchmarking. Moreover, this platform offers a centralized data management, automatic software updates and access to high performance computers accessible with desktop computers and mobile devices.

Application of Imitated Ecosystems in Urban Waterlogging Prevention

In the constantly advancing of urbanization, waterlogging disasters causing a serious impact on urban living, economic development, harmony and stability, etc.In order to solve these issues brought by urban waterlogging , this paper will combine the actual urban construction, accord the characteristics and causes of waterlogging, draw lessons from the water cycle model of natural ecosystems, ,and apply the imitation of ecosystems in urban waterlogging Disaster Prevention, to explore an efficient and orderly waterlogging prevention system, and strive to solve the problem fundamentally.

Exploring the relationship between peatland net carbon balance and apparent carbon accumulation rate at century to millennial time scales

Each year, a peatland has an annual net carbon balance ( NCB), which can be positive (net uptake), zero or negative. Over centuries to millennia, this NCB accumulates as a peat profile. Contemporary peatlands can be sampled (cored), and the past apparent carbon accumulation rate ( aCAR) can be determined as the quantity of peat carbon in any particular dated interval down the core profile. We use a process-based peatland carbon and water cycle model to compare peatland annual NCB during millennia of peat accumulation to the contemporary estimate of aCAR, resulting from this accumulation. Integrating over the entire profile, the accumulated NCB must equal the aCAR, but for shorter time intervals, these two quantities can diverge. A climate variation/perturbation that leads to persistent, slow carbon loss or negligible carbon gain through enhanced decomposition will necessarily reduce the aCAR for time periods before the climate variation/perturbation occurred. This can compromise peatland climate–carbon balance relationships inferred from joint analysis of peat cores and paleoclimate reconstructions.

Research on Driving Mechanism of Socio-Economic Factors for Water Cycle

The paper selects water price as a potentially effective socio-economic factor for the mitigation of conflict between socio-economic development and regional water resources in the Weihe River basin on basis of the driving mechanism research of socio-economic factors for water cycle, and evaluates the variation of socio-economic and hydrological variables by comprehensive simulation of economic system and water resources system by coupling of economic model CGE, water allocation model ROWAS and distributed water cycle model WEP. The results indicate that raising water price would bring down socio-economic development and regional water demand to some extent, but regional water resources tensions just reduce so slightly that corresponding schemes turn out to be not sufficient to deal with water crisis in the Weihe River basin.

Impacts of Water Consumption Structure Adjustment on Social-Economic Development and Water Cycle

Water consumption structure adjustment is considered to be effective to alleviate regional water crisis for sustainable socio-economic development and rational water resources utilization. The paper studies mechanism of interactions between socio-economic system and water resources system by coupling of economic model CGE, water allocation model ROWAS and distributed water cycle model WEP to discover the impacts of water consumption structure adjustment on regional socio-economic development and water cycle in the Weihe River basin of China. The results indicate that water consumption structure adjustment transferring water consumption from agriculture to industry and services would distinctly promote economic growth and resident welfare with surface water flow increasing and groundwater level declining. As a result, water supply structure adjustment and integrated regulation of surface water and groundwater would be essential for regional harmonious development.

Water cycle modelling of peri-urban hydroecological systems: a case study of the South Creek catchment, Australia

Modelling offers a time efficient and cost-effective approach to integrate and analyse different components of complex peri-urban hydroecological systems. This paper evaluates and applies a catchment water balance model (Peri-Urban SimHyd) developed to simulate the complete water cycle of the South Creek catchment, a complex peri-urban catchment in Western Sydney (Australia). Peri-Urban SimHyd combined the modelling of daily rainfall-runoff from both impervious and pervious surfaces with the modelling of monthly reticulated potable water supply, wastewater flows, and surface and groundwater extractions, to simulate the monthly water balance of the catchment. The model was successful (the Nash-Sutcliffe efficiency >0.68) in simulating the monthly total wastewater discharges and streamflows in the South Creek and its main tributary, the Eastern Creek, over a period of 15 years (1992–2006). The model performance was found to be suitable for simulating macro water balance and catchment yield but not for low streamflow hydrographs. Further, the model demonstrated its applicability in assessing the impacts of water cycle management strategies, i.e. ‘replacement flows’ and ‘water smart farms’ scenarios, on monthly and annual streamflows in the South Creek catchment. In general, Peri-Urban SimHyd has a generic structure and can be employed to assess impacts of various ‘what if’ scenarios on water cycle and streamflows in other peri-urban hydroecological systems.

Measurement of δ18O, δ17O, and 17O-excess in Water by Off-Axis Integrated Cavity Output Spectroscopy and Isotope Ratio Mass Spectrometry

Stable isotopes of water have long been used to improve understanding of the hydrological cycle, catchment hydrology, and polar climate. Recently, there has been increasing interest in measurement and use of the less-abundant (17)O isotope in addition to (2)H and (18)O. Off-axis integrated cavity output spectroscopy (OA-ICOS) is demonstrated for accurate and precise measurements δ(18)O, δ(17)O, and (17)O-excess in liquid water. OA-ICOS involves no sample conversion and has a small footprint, allowing measurements to be made by researchers collecting the samples. Repeated (514) high-throughput measurements of the international isotopic reference water standard Greenland Ice Sheet Precipitation (GISP) demonstrate the precision and accuracy of OA-ICOS: δ(18)OVSMOW-SLAP = -24.74 ± 0.07‰ (1σ) and δ(17)OVSMOW-SLAP = -13.12 ± 0.05‰ (1σ). For comparison, the International Atomic Energy Agency (IAEA) value for δ(18)OVSMOW-SLAP is -24.76 ± 0.09‰ (1σ) and an average of previously reported values for δ(17)OVSMOW-SLAP is -13.12 ± 0.06‰ (1σ). Multiple (26) high-precision measurements of GISP provide a (17)O-excessVSMOW-SLAP of 23 ± 10 per meg (1σ); an average of previously reported values for (17)O-excessVSMOW-SLAP is 22 ± 11 per meg (1σ). For all these OA-ICOS measurements, precision can be further enhanced by additional averaging. OA-ICOS measurements were compared with two independent isotope ratio mass spectrometry (IRMS) laboratories and shown to have comparable accuracy and precision as the current fluorination-IRMS techniques in δ(18)O, δ(17)O, and (17)O-excess. The ability to measure accurately δ(18)O, δ(17)O, and (17)O-excess in liquid water inexpensively and without sample conversion is expected to increase vastly the application of δ(17)O and (17)O-excess measurements for scientific understanding of the water cycle, atmospheric convection, and climate modeling among others.

Dynamic mapping of snow-free vegetation and bare soil albedos at global 1 km scale from 10-year analysis of MODIS satellite products

Land surface albedo is considered to be a fundamental quantity for accurate assessment of the surface energy budget. Whereas satellite products usually provide total surface albedo, land surface models often require a separation of snow-free surface albedos of the vegetation canopy and the underlying bare soil. In this paper, visible (0.3–0.7μm) and near infrared (0.7–5.0μm) vegetation and bare soil albedos are retrieved from a decade of MODIS data at global scale by using a Kalman Filter-based method. The retrieved quantities show good agreement with those derived in previous studies. The main advantage of the proposed method over other approaches is that it provides albedo components efficiently throughout the year, thus allowing the seasonal cycles to be captured. This study also exploits the relationship linking the retrieved bare soil albedo to soil moisture in order to reinforce indirectly the quality of the albedo splitting. The proposed method makes it possible to update the bare soil and vegetation albedo fields in the ECOCLIMAP land cover map, the database of physiographic parameters used by Météo-France for climate and NWP modelling. The mean relative difference between the new and the old ECOCLIMAP shortwave albedos (0.3–5.0μm) is −8.5% at the global scale if only albedos lower than 0.5 are considered, and 41.2% otherwise. This update of the ECOCLIMAP database will have numerous NWP and climate applications for carbon, energy, and water cycle modelling.

Spatial and temporal variability of the interception in the natural wetland valley, the lower Biebrza basin case study

Abstract Spatial and temporal variability of the interception in the natural wetland valley, the lower Biebrza basin case study. The paper presents the research carried out in the lower basin of Biebrza River valley in order to identify interception for natural wetland plant communities. Maximum interception, i.e. the largest amount of water, expressed in millimeters, which can be captured and retained by plant canopy from rainfall is one of the key parameters of the water cycle modeling. Maximum interception was determined based on the difference of the masses of wet and dry fresh plant samples. Collection of plant material samples took place during the five measurement sessions, which began immediately after the flood recedes, and then lasted until the end of the growing season. Interception spatial variability was analyzed on the basis of the results of maximum interception measured for selected plant aggregations in the different sampling points. The obtained values were extrapolated to the area of the lower basin of Biebrza River using vegetation map of the Biebrza National Park. By conducting a test sessions in the five coming months, the maps of the spatial variability also show changes over time. Methodology used in the described tests allowed for obtaining of satisfactory results. They present, in a correct way, variation occurring between the plant aggregations due to their morphology. In most cases the results are consistent with data from the literature. As results of the analysis of spatial variability of the maximum interception, the highest values were found for the plant communities located in the immediate vicinity of the river channel. With the increase of the distance from river towards the valley edges the maximum interception values decrease. These changes can be seen in the form of strips parallel to the river channel, which corresponds to the plant zones. Obtained map of spatial variability of the maximum interception, which is the results of extrapolation of the values assigned to plant communities, has a high correlation with the map resulting from the analysis of satellite images

Integrated urban water modelling using Aquacycle model

The traditional approach to urban water management is primarily based on a supply-oriented approach, where water follows a one-way path from supply to a single use, treatment and then discharged to the environment. The wide acknowledgement of the need to shift towards more sustainable practices has led to the recognition that demand is multi-faceted; delivered supply should correspond to the varying characteristics of users, taking into account their actual needs in terms of quantity, quality and level of reliability. Such an approach can entail the application of reuse and reclamation methods, which allow for the multiple use of water to meet higher to lower quality needs. Integrated urban water cycle modelling can be a powerful tool in estimating the potential for the application of such options, as it considers issues related to water supply, wastewater production and stormwater runoff in a holistic framework. This paper presents preliminary results from the application of an integrated urban water balance model, Aquacycle, in the Greater Athens Area. The Aquacycle model accounts for water pathways by simulating two subsystems of the urban water cycle (the rainfall-run-off network and the water supplywastewater network), and the interactions between them. The software package can also be used for evaluating alternative stormwater and wastewater reuse schemes. In this paper, Aquacycle was used for simulating water use, wastewater production and stormwater drainage in the Greater Athens Area. For applying the model, the region was divided into smaller – cluster – scales, according to specific structural characteristics. Then, the model was calibrated and validated against measured data, so that model outputs can account for the overall complexity of the Athens urban environment. Finally, three alternative scenarios, comprising the application of on-site recycling and reuse technologies were simulated and evaluated on the basis of three sustainability indicators (i.e. reduction of imported freshwater, wastewater and stormwater discharge).

Source to tap urban water cycle modelling

The continuous expansion of urban areas is associated with increased water demand, both for domestic and non-domestic uses. To cover this additional demand, centralised infrastructure, such as water supply and distribution networks tend to become more and more complicated and are eventually over-extended with adverse effects on their reliability. To address this, there exist two main strategies: (a) Tools and algorithms are employed to optimise the operation of the external water supply system, in an effort to minimise risk of failure to cover the demand (either due to the limited availability of water resources or due to the limited capacity of the transmission system and treatment plants) and (b) demand management is employed to reduce the water demand per capita. Dedicated tools do exist to support the implementation of these two strategies separately. However, there is currently no tool capable of handling the complete urban water system, from source to tap, allowing for an investigation of these two strategies at the same time and thus exploring synergies between the two. This paper presents a new version of the UWOT model (Makropoulos et al., 2008), which adopts a metabolism modelling approach and is now capable of simulating the complete urban water cycle from source to tap and back again: the tool simulates the whole water supply network from the generation of demand at the household level to the water reservoirs and tracks wastewater generation from the household through the wastewater system and the treatment plants to the water bodies. UWOT functionality is demonstrated in the case of the water system of Athens and outputs are compared against the current operational tool used by the Water Company of Athens. Results are presented and discussed: The discussion highlights the conditions under which a single source-to-tap model is more advantageous than dedicated subsystem models.