Conceptual modelling to assess how the interplay of hydrological connectivity, catchment storage and tracer dynamics controls nonstationary water age estimates

Acknowledgements We would like to gratefully acknowledge the data provided by SEPA, Iain Malcolm. Mark Speed, Susan Waldron and many MSS staff helped with sample collection and lab analysis. We thank the European Research Council (project GA 335910 VEWA) for funding and are grateful for the constructive comments provided by three anonymous reviewers.

Avaliação de incertezas na calibração automática do modelo SWMM

Este trabalho é referente ao desenvolvimento de um calibrador multiobjetivo automático do modelo SWMM (Storm Water Management Model), e avaliação de algumas fontes de incertezas presentes no processo de calibração, visando à representação satisfatória da transformação chuva-vazão. O código foi escrito em linguagem C, e aplica os conceitos do método de otimização multiobjetivo NSGAII (Non Dominated Sorting Genetic Algorithm) com elitismo controlado, além de utilizar o código fonte do modelo SWMM para a determinação das vazões simuladas. Paralelamente, também foi criada uma interface visual, para melhorar a facilidade de utilização do calibrador. Os testes do calibrador foram aplicados a três sistemas diferentes: um sistema hipotético disponibilizado no pacote de instalação do SWMM; um sistema real de pequenas dimensões, denominado La Terraza, localizado no município de Sierra Vista, Arizona (EUA); e um sistema de maiores dimensões, a bacia hidrográfica do Córrego do Gregório, localizada no município de São Carlos (SP). Os resultados indicam que o calibrador construído apresenta, em geral, eficiência satisfatória, porém é bastante dependente da qualidade dos dados observados em campo e dos parâmetros de entrada escolhidos pelo usuário. Foi demonstrada a importância da escolha dos eventos utilizados na calibração, do estabelecimento de limites adequados nos valores das variáveis de decisão, da escolha das funções objetivo e, principalmente, da qualidade e representatividade dos dados de monitoramento pluvio e fluviométrico. Conclui-se que estes testes desenvolvidos contribuem para o entendimento mais aprofundado dos processos envolvidos na modelagem e calibração, possibilitando avanços na confiabilidade dos resultados da modelagem.

Efficient accommodation of local minima in watershed model calibration

The Gauss-Marquardt-Levenberg (GML) method of computer-based parameter estimation, in common with other gradient-based approaches, suffers from the drawback that it may become trapped in local objective function minima, and thus report optimized parameter values that are not, in fact, optimized at all. This can seriously degrade its utility in the calibration of watershed models where local optima abound. Nevertheless, the method also has advantages, chief among these being its model-run efficiency, and its ability to report useful information on parameter sensitivities and covariances as a by-product of its use. It is also easily adapted to maintain this efficiency in the face of potential numerical problems (that adversely affect all parameter estimation methodologies) caused by parameter insensitivity and/or parameter correlation. The present paper presents two algorithmic enhancements to the GML method that retain its strengths, but which overcome its weaknesses in the face of local optima. Using the first of these methods an intelligent search for better parameter sets is conducted in parameter subspaces of decreasing dimensionality when progress of the parameter estimation process is slowed either by numerical instability incurred through problem ill-posedness, or when a local objective function minimum is encountered. The second methodology minimizes the chance of successive GML parameter estimation runs finding the same objective function minimum by starting successive runs at points that are maximally removed from previous parameter trajectories. As well as enhancing the ability of a GML-based method to find the global objective function minimum, the latter technique can also be used to find the locations of many non-global optima (should they exist) in parameter space. This can provide a useful means of inquiring into the well-posedness of a parameter estimation problem, and for detecting the presence of bimodal parameter and predictive probability distributions. The new methodologies are demonstrated by calibrating a Hydrological Simulation Program-FORTRAN (HSPF) model against a time series of daily flows. Comparison with the SCE-UA method in this calibration context demonstrates a high level of comparative model run efficiency for the new method. (c) 2006 Elsevier B.V. All rights reserved.

Quantifying simulator discrepancy in discrete-time dynamical simulators

When making predictions with complex simulators it can be important to quantify the various sources of uncertainty. Errors in the structural specification of the simulator, for example due to missing processes or incorrect mathematical specification, can be a major source of uncertainty, but are often ignored. We introduce a methodology for inferring the discrepancy between the simulator and the system in discrete-time dynamical simulators. We assume a structural form for the discrepancy function, and show how to infer the maximum-likelihood parameter estimates using a particle filter embedded within a Monte Carlo expectation maximization (MCEM) algorithm. We illustrate the method on a conceptual rainfall-runoff simulator (logSPM) used to model the Abercrombie catchment in Australia. We assess the simulator and discrepancy model on the basis of their predictive performance using proper scoring rules. This article has supplementary material online. © 2011 International Biometric Society.

A three dimensional numerical investigation of hillslope flow processes

Physically based distributed models of catchment hydrology are likely to be made available as engineering tools in the near future. Although these models are based on theoretically acceptable equations of continuity, there are still limitations in the present modelling strategy. Of interest to this thesis are the current modelling assumptions made concerning the effects of soil spatial variability, including formations producing distinct zones of preferential flow. The thesis contains a review of current physically based modelling strategies and a field based assessment of soil spatial variability. In order to investigate the effects of soil nonuniformity a fully three dimensional model of variability saturated flow in porous media is developed. The model is based on a Galerkin finite element approximation to Richards equation. Accessibility to a vector processor permits numerical solutions on grids containing several thousand node points. The model is applied to a single hillslope segment under various degrees of soil spatial variability. Such variability is introduced by generating random fields of saturated hydraulic conductivity using the turning bands method. Similar experiments are performed under conditions of preferred soil moisture movement. The results show that the influence of soil variability on subsurface flow may be less significant than suggested in the literature, due to the integrating effects of three dimensional flow. Under conditions of widespread infiltration excess runoff, the results indicate a greater significance of soil nonuniformity. The recognition of zones of preferential flow is also shown to be an important factor in accurate rainfall-runoff modelling. Using the results of various fields of soil variability, experiments are carried out to assess the validity of the commonly used concept of `effective parameters'. The results of these experiments suggest that such a concept may be valid in modelling subsurface flow. However, the effective parameter is observed to be event dependent when the dominating mechanism is infiltration excess runoff.

(Figure 2) Mean discharge per day (total and baseflow part) of the San Francisco River and daily precipitation at station ECPL (Planta), in 2007

We investigated controls on the water chemistry of a South Ecuadorian cloud forest catchment which is partly pristine, and partly converted to extensive pasture. From April 2007 to May 2008 water samples were taken weekly to biweekly at nine different subcatchments, and were screened for differences in electric conductivity, pH, anion, as well as element composition. A principal component analysis was conducted to reduce dimensionality of the data set and define major factors explaining variation in the data. Three main factors were isolated by a subset of 10 elements (Ca2+, Ce, Gd, K+, Mg2+, Na+, Nd, Rb, Sr, Y), explaining around 90% of the data variation. Land-use was the major factor controlling and changing water chemistry of the subcatchments. A second factor was associated with the concentration of rare earth elements in water, presumably highlighting other anthropogenic influences such as gravel excavation or road construction. Around 12% of the variation was explained by the third component, which was defined by the occurrence of Rb and K and represents the influence of vegetation dynamics on element accumulation and wash-out. Comparison of base- and fast flow concentrations led to the assumption that a significant portion of soil water from around 30 cm depth contributes to storm flow, as revealed by increased rare earth element concentrations in fast flow samples. Our findings demonstrate the utility of multi-tracer principal component analysis to study tropical headwater streams, and emphasize the need for effective land management in cloud forest catchments.

Desempenho de uma lagoa de infiltração na absorção de cheias e na recarga de aquífero

A type of macro drainage solution widely used in urban areas with predomi-nance of closed catchments (basins without outlet) is the implementation of detention and infiltration reservoirs (DIR). This type of solution has the main function of storing surface runoff and to promote soil infiltration and, consequently, aquifer recharge. The practice is to avoid floods in the drainage basin low-lying areas. The catchment waterproofing reduces the distributed groundwater recharge in urban areas, as is the case of Natal city, RN. However, the advantage of DIR is to concentrate the runoff and to promote aquifer recharge to an amount that can surpass the distributed natu-ral recharge. In this paper, we proposed studying a small urban drainage catchment, named Experimental Mirassol Watershed (EMW) in Natal, RN, whose outlet is a DIR. The rainfall-runoff transformation processes, water accumulation in DIR and the pro-cess of infiltration and percolation in the soil profile until the free aquifer were mod-eled and, from rainfall event observations, water levels in DIR and free aquifer water level measurements, and also, parameter values determination, it is was enabled to calibrate and modeling these combined processes. The mathematical modeling was carried out from two numerical models. We used the rainfall-runoff model developed by RIGHETTO (2014), and besides, we developed a one-dimensional model to simu-late the soil infiltration, percolation, redistribution soil water and groundwater in a combined system to the reservoir water balance. Continuous simulation was run over a period of eighteen months in time intervals of one minute. The drainage basin was discretized in blocks units as well as street reaches and the soil profile in vertical cells of 2 cm deep to a total depth of 30 m. The generated hydrographs were transformed into inlet volumes to the DIR and then, it was carried out water balance in these time intervals, considering infiltration and percolation of water in the soil profile. As a re-sult, we get to evaluate the storage water process in DIR as well as the infiltration of water, redistribution into the soil and the groundwater aquifer recharge, in continuous temporal simulation. We found that the DIR has good performance to storage excess water drainage and to contribute to the local aquifer recharge process (Aquifer Dunas / Barreiras).

Modelling storage-driven connectivity between landscapes and riverscapes : towards a simple framework for long-term ecohydrological assessment

Acknowledgements: We thank Iain Malcolm of Marine Scotland Science for access to data from the Girnock and the Scottish Environment Protection Agency for historical stage-discharge relationships. CS contributions on this paper were in part supported by the NERC/JPI SIWA project (NE/M019896/1).

Measuring eco-roof mitigation on flash floods via GIS simulation

Purpose: The rapid and ongoing expansion of urbanised impervious areas could lead to more frequent flood inundation in urban flood-prone regions. Nowadays, urban flood inundation induced by rainstorm is an expensive natural disaster in many countries. In order to reduce the flooding risk, eco-roof systems (or green roof systems) could be considered as an effective mechanism of mitigating flooding disasters through their rainwater retention capability. However, there is still a lack of examining the stormwater management tool. The purpose of this paper is to evaluate the effects on flooding disaster from extensive green roofs. Design/methodology/approach: Based on geographical information system (GIS) simulation, this research presents a frame of assessing eco-roof impacts on urban flash floods. The approach addresses both urban rainfall-runoff and underground hydrologic models for traditional impervious and green roofs. Deakin University’s Geelong Waurn Ponds campus is chosen as a study case. GIS technologies are then utilised to visualise and analyse the effects on flood inundation from surface properties of building roofs. Findings: The results reveal that the eco-roof systems generate varying degrees of mitigation of urban flood inundation with different return period storms. Originality/value: Although the eco-roof technology is considered as an effective stormwater management tool, it is not commonly adopted and examined in urban floods. This study will bring benefits to urban planners for raising awareness of hazard impacts and to construction technicians for considering disaster mitigation via roof technologies. The approach proposed here could be used for the disaster mitigation in future urban planning.

Water Demand and Allocation in the Mara River Basin, Kenya/Tanzania in the Face of Land Use Dynamics and Climate Variability

The Mara River Basin (MRB) is endowed with pristine biodiversity, socio-cultural heritage and natural resources. The purpose of my study is to develop and apply an integrated water resource allocation framework for the MRB based on the hydrological processes, water demand and economic factors. The basin was partitioned into twelve sub-basins and the rainfall runoff processes was modeled using the Soil and Water Assessment Tool (SWAT) after satisfactory Nash-Sutcliff efficiency of 0.68 for calibration and 0.43 for validation at Mara Mines station. The impact and uncertainty of climate change on the hydrology of the MRB was assessed using SWAT and three scenarios of statistically downscaled outputs from twenty Global Circulation Models. Results predicted the wet season getting more wet and the dry season getting drier, with a general increasing trend of annual rainfall through 2050. Three blocks of water demand (environmental, normal and flood) were estimated from consumptive water use by human, wildlife, livestock, tourism, irrigation and industry. Water demand projections suggest human consumption is expected to surpass irrigation as the highest water demand sector by 2030. Monthly volume of water was estimated in three blocks of current minimum reliability, reserve (>95%), normal (80–95%) and flood (40%) for more than 5 months in a year. The assessment of water price and marginal productivity showed that current water use hardly responds to a change in price or productivity of water. Finally, a water allocation model was developed and applied to investigate the optimum monthly allocation among sectors and sub-basins by maximizing the use value and hydrological reliability of water. Model results demonstrated that the status on reserve and normal volumes can be improved to ‘low’ or ‘moderate’ by updating the existing reliability to meet prevailing demand. Flow volumes and rates for four scenarios of reliability were presented. Results showed that the water allocation framework can be used as comprehensive tool in the management of MRB, and possibly be extended similar watersheds.

Modeling on Runoff Concentration Caused by Rainfall on Hillslopes and Application in Maoping Slope

Based on the fact that the concentration flowlines of overland flow depend on the surface landform of hillslope, a kinematic wave model was developed for simulating runoff generation and flow concentration caused by rainfall on hillslopes. The model-simulated results agree well with experimental observations. Applying the model to the practical case of Maoping slope, we obtained the characteristics of runoff generation and infiltration on the slope. Especially, the simulated results adequately reflected the confluent pattern of surface runoff, which offers a scientific foundation for designing the drainage engineering on the Maoping slope.

Modeling of Hillslope Runoff and Soil Erosion at Rainfall Events Using Cellular Automata Approach

National Science Fund for Distinguished Young Scholars of China [40225004]; National Natural Science Foundation of China [40471048]

Rainfall and runoff water quality of the Pang and Lambourn, tributaries of the River Thames, south-eastern England

The water quality of rainfall and runoff is described for two catchments of two tributaries of the River Thames, the Pang and Lambourn. Rainfall chemistry is variable and concentrations of most determinands decrease with increasing volume of catch probably due to 'wash out' processes. Two rainfall sites have been monitored, one for each catchment. The rainfall site on the Lambourn shows higher chemical concentrations than the one for the Pang which probably reflects higher amounts of local inputs from agricultural activity, Rainfall quality data at a long-term rainfall site on the Pang (UK National Air Quality Archive) shows chemistries similar to that for the Lambourn site. but with some clear differences. Rainfall chemistries show considerable variation on an event-to-event basis. Average water quality concentrations and flow-weighted concentrations as well as fluxes vary across the sites, typically by about 30%. Stream chemistry is much less variable due to the main Source of water coming from aquifer sources of high storage. The relationship between rainfall and runoff chemistry at the catchment outlet is described in terms of the relative proportions of atmospheric and within-catchment sources. Remarkably, in view of the quantity of agricultural and sewage inputs to the streams, the catchments appear to be retaining both P and N.

Characterization of Soil Surface Roughness Effects on Runoff and Soil Erosion Rates under Simulated Rainfall

This study aims at identifying the influence of soil surface roughness from small to large aggregates (random roughness) on runoff and soil loss and to investigate the interaction with soil surface seal formation. Bulk samples of a silty clay loam soil were sieved to four aggregate-size classes of 3 to 12, 12 to 20, 20 to 45, 45 to 100 mm, and packed in soil trays set at a 5% slope. Rainfall simulations using an oscillating nozzle simulator were conducted for 90 min at an average rainfall intensity of 50.2 mm h(-1). Soil surface roughness was measured using an instantaneous profile laser scanner and surface sealing was studied by macroscopic analysis of epoxy impregnated soil samples. The rainfall simulations revealed longer times to initiate runoff with increasing soil surface roughness. For random roughness levels up to 6 mm, a decrease in final runoff rate with increasing roughness was observed. This can be attributed to a decreased breakdown of the larger roughness elements on rougher surfaces, thus keeping infiltration rate high. For a random roughness larger than 6 mm, a greater final runoff rate was observed. This was caused by the creation of a thick depositional seal in the concentrated flow areas, thus lowering the infiltration rates. Analysis of impregnated soil sample blocks confirmed the formation of a structural surface seal on smooth surfaces, whereas thick depositional seals were visible in the depressional areas of rougher surfaces. Therefore, from our observations it can be learned that soil surface roughness as formed by the presence of different aggregate sizes reduces runoff but that its effect diminishes due to aggregate breakdown and the formation of thick depositional seals in the case of rough soil surfaces. Sediment concentration increased with increasing soil surface roughness, due to runoff concentration in flow paths. Nevertheless, final soil loss rates were comparable for all soil roughness categories, indicating that random roughness is only important in influencing runoff rates and the time to initiate runoff, but not in influencing sediment export through soil loss rates.