Model and parameter uncertainty in IDF relationships under climate change

Quantifying distributional behavior of extreme events is crucial in hydrologic designs. Intensity Duration Frequency (IDF) relationships are used extensively in engineering especially in urban hydrology, to obtain return level of extreme rainfall event for a specified return period and duration. Major sources of uncertainty in the IDF relationships are due to insufficient quantity and quality of data leading to parameter uncertainty due to the distribution fitted to the data and uncertainty as a result of using multiple GCMs. It is important to study these uncertainties and propagate them to future for accurate assessment of return levels for future. The objective of this study is to quantify the uncertainties arising from parameters of the distribution fitted to data and the multiple GCM models using Bayesian approach. Posterior distribution of parameters is obtained from Bayes rule and the parameters are transformed to obtain return levels for a specified return period. Markov Chain Monte Carlo (MCMC) method using Metropolis Hastings algorithm is used to obtain the posterior distribution of parameters. Twenty six CMIP5 GCMs along with four RCP scenarios are considered for studying the effects of climate change and to obtain projected IDF relationships for the case study of Bangalore city in India. GCM uncertainty due to the use of multiple GCMs is treated using Reliability Ensemble Averaging (REA) technique along with the parameter uncertainty. Scale invariance theory is employed for obtaining short duration return levels from daily data. It is observed that the uncertainty in short duration rainfall return levels is high when compared to the longer durations. Further it is observed that parameter uncertainty is large compared to the model uncertainty. (C) 2015 Elsevier Ltd. All rights reserved.

Impact of river runoff into the ocean on Indian summer monsoon

Rivers of the world discharge about 36000 km 3 of freshwater into the ocean every year. To investigate the impact of river discharge on climate, we have carried out two 100 year simulations using the Community Climate System Model (CCSM3), one including the river runoff into the ocean and the other excluding it. When the river discharge is shut off, global average sea surface temperature (SST) rises by about 0.5 degrees C and the Indian Summer Monsoon Rainfall (ISMR) increases by about 10% of the seasonal total with large increase in the eastern Bay of Bengal and along the west coast of India. In addition, the frequency of occurrence of La Nina-like cooling events in the equatorial Pacific increases and the correlation between ISMR and Pacific SST anomalies become stronger. The teleconnection between the SST anomalies in the Pacific and monsoon is effected via upper tropospheric meridional temperature gradient and the North African-Asian Jet axis.

Modeling on Runoff and Sediment Yield in Small Watersheds

The prediction and estimate of water and soil loss is fundamental important for understanding the effect of the spatial heterogeneity of underlying surfaces and preventing ecological environment deterioration. In this paper, a dynamic model of runoff and sediment yield in small watersheds is established. The proposed model includes three components: runoff generation caused by rainfall, soil erosion on hillslopes by overland flow, and runoff concentration and sediment transport on watersheds. Applying the proposed model, the runoff and sediment yield processes in a typical catchment on the loess plateau was estimated, which exhibited a good agreement between predicted results and observation.

Model structure and uncertainty for stochastic non-point source modelling applications

The uncertainty associated with a rainfall-runoff and non-point source loading (NPS) model can be attributed to both the parameterization and model structure. An interesting implication of the areal nature of NPS models is the direct relationship between model structure (i.e. sub-watershed size) and sample size for the parameterization of spatial data. The approach of this research is to find structural limitations in scale for the use of the conceptual NPS model, then examine the scales at which suitable stochastic depictions of key parameter sets can be generated. The overlapping regions are optimal (and possibly the only suitable regions) for conducting meaningful stochastic analysis with a given NPS model. Previous work has sought to find optimal scales for deterministic analysis (where, in fact, calibration can be adjusted to compensate for sub-optimal scale selection); however, analysis of stochastic suitability and uncertainty associated with both the conceptual model and the parameter set, as presented here, is novel; as is the strategy of delineating a watershed based on the uncertainty distribution. The results of this paper demonstrate a narrow range of acceptable model structure for stochastic analysis in the chosen NPS model. In the case examined, the uncertainties associated with parameterization and parameter sensitivity are shown to be outweighed in significance by those resulting from structural and conceptual decisions. © 2011 Copyright IAHS Press.

分布式水文模型的误差分析

分布式水文模型以其具有明确物理意义的参数结构和对空间分异性的全面反映,能够准确详尽地描述和模拟流域内真实的降水径流过程而被广泛需求和关注。在模拟土地利用、土地覆盖、水土流失等各种变化过程的水文响应,面源污染、陆面过程、气候变化影响评价等诸多领域都有广泛的应用。模型的预报精度和误差至关重要,决定了模型的应用和推广。在分析分布式水文模型建立和验证过程的基础上,提出了模型的4类误差来源:被排除在外的因素引起的误差,实测历史记录资料的随机或系统误差,参数误差和模型结构误差,讨论了各类误差的分析与计算方法,为模型的发展和成长提供了依据。

Effects of soil type on leaching and runoff transport of rare earth elements and phosphorous in laboratory experiments

Through leaching experiments and simulated rainfall experiments, characteristics of vertical leaching of exogenous rare earth elements (REEs) and phosphorus (P) and their losses with surface runoff during simulated rainfall in different types of soils (terra nera soil, cinnamon soil, red soil, loess soil, and purple soil) were investigated. Results of the leaching experiments showed that vertical transports of REEs and P were relatively low, with transport depths less than 6 cm. The vertical leaching rates of REEs and P in the different soils followed the order of purple soil > terra nera soil > red soil > cinnamon soil > loess soil. Results of the simulated rainfall experiments (83 mm h(-1)) revealed that more than 92% of REEs and P transported with soil particles in runoff. The loss rates of REEs and P in surface runoff in the different soil types were in the order of loess soil > terra nera soil > cinnamon soil > red soil > purple soil. The total amounts of losses of REEs and P in runoff were significantly correlated.

Scale-wise evolution of rainfall probability density functions fingerprints the rainfall generation mechanism

© 2010 by the American Geophysical Union.The cross-scale probabilistic structure of rainfall intensity records collected over time scales ranging from hours to decades at sites dominated by both convective and frontal systems is investigated. Across these sites, intermittency build-up from slow to fast time-scales is analyzed in terms of heavy tailed and asymmetric signatures in the scale-wise evolution of rainfall probability density functions (pdfs). The analysis demonstrates that rainfall records dominated by convective storms develop heavier-Tailed power law pdfs toward finer scales when compared with their frontal systems counterpart. Also, a concomitant marked asymmetry build-up emerges at such finer time scales. A scale-dependent probabilistic description of such fat tails and asymmetry appearance is proposed based on a modified q-Gaussian model, able to describe the cross-scale rainfall pdfs in terms of the nonextensivity parameter q, a lacunarity (intermittency) correction and a tail asymmetry coefficient, linked to the rainfall generation mechanism.

Assessing the combined effect of land management and wildfire on runoff and soil erosion in north central Portugal

Strong and sometimes extreme responses in runoff and soil erosion following wildfires have been reported worldwide. However, in the case of North-Central Portugal, little research had been carried out regarding the hydrologic and erosive impacts of several land management activities in recently burnt areas (such as ground preparation, post-fire logging or post-fire mitigation treatments). This study aims to assess post-fire runoff and soil erosion response on Eucalypt and Maritime pine plantations during the first, second and third years following wildfires. The effect of several pre-fire ground preparation operations (ploughed down-slope, contour ploughed and inclined terraces), post-fire logging activities (on both the eucalypt and pine plantations), as well as the application of hydromulch (a post-fire emergency treatment) on overland flow and soil erosion were compared to burnt but undisturbed and untreated areas. The intensive monitoring of runoff, soil erosion and selected soil properties served to determine the main factors involved in post-fire runoff and soil erosion and their spatial and temporal variation. Soil water repellency deserved special attention, due to its supposed important role for overland flow generation. Repeated rainfall simulation experiments (RSE’s), micro-scale runoff plots and bounded sediment fences were carried out and/or installed immediately after the wildfire on seven burnt slopes. Micro-scale runoff plots results under natural rainfall conditions were also compared to the RSE’s results, which was useful for assessing the representativeness of the data obtained with artificial rainfall. The results showed comparable runoff coefficient (20-60%) but lower sediment losses (125-1000 g m-2) than prior studies in Portugal, but especially outside Portugal. Lower sediment losses were related with the historic intensive land use in the area. In evaluating these losses, however, the shallowness and stoniness of the soils, as well as the high organic matter fraction of the eroded sediments (50%) must not be overlooked. Sediment limited erosion was measured in all the ploughed sites, probably due to the time since ploughing (several years). The disturbance of the soil surface cover due to post-fire logging and wood extraction substantially increased sediment losses at both the pine and eucalypt sites. Hydromulch effectiveness in reducing the runoff (70%) and sediment losses (83%) was attributed to the protective high coverage provided by hydromulch. The hydromulch significantly affected the soil cover and other soil properties and these changes also reduced the soil erosion risk. The rainfall amount was the main factor explaining the variance in runoff. However, a shift from rainfall amount to rainfall intensity was detected when either the surface cover or the infiltration capacity (hydrophilic conditions) increased. Sediment losses were controlled by rainfall intensity and surface cover. The role of soil water repellency on runoff generation was not consistent; the overall repellency levels alone were not enough to assess its hydrological impact. Soil water repellency explained runoff generation in the specific-sites model better than in the overall model. Additionally, soil moisture content was a better predictor for soil water repellency than antecedent rainfall. The natural rainfall results confirmed that RSE’s were able to capture the specific sediment losses and its organic matter content as well as the differences between the ploughed and unploughed sites. Repeated RSE’s also captured the seasonal variations in runoff and sediment losses attributed to soil water repellency. These results have implications for post-fire soil erosion modelling and soil conservation practices in the region, or areas with the same land use, climate and soil characteristics. The measured sediment loss, as well as the increasing frequency of ploughing in recently burnt and unburnt eucalypt stands, suggests ploughing is not an effective as a soil conservation measure. Logging activities with less impact are recommended in order to maintain the forest litter protecting the soil surface. Due to its high effectiveness in reducing runoff and soil erosion, hydromulch is recommended for highly sensitive and vulnerable areas.

Assessing soil erosion due to land use change at the Alqueva reservoir surrounding area

Tese de doutoramento, Ciências do Mar, da Terra e do Ambiente, Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2015

Developing a desertification indicator system for a small Mediterranean catchment: a case study from the Serra de Rodes, Alt Empordà, Catalunya, NE Spain

La desertificació és un problema de degradació de sòls de gran importància en regions àrides, semi-àrides i sub-humides, amb serioses conseqüències ambientals, socials i econòmiques com a resultat de l'impacte d'activitats humanes en combinació amb condicions físiques i medi ambientals desfavorables (UNEP, 1994). L'objectiu principal d'aquesta tesi va ser el desenvolupament d'una metodologia simple per tal de poder avaluar de forma precisa l'estat i l'evolució de la desertificació a escala local, a través de la creació d'un model anomenat sistema d'indicators de desertificació (DIS). En aquest mateix context, un dels dos objectius específics d'aquesta recerca es va centrar en l'estudi dels factors més importants de degradació de sòls a escala de parcel.la, comportant un extens treball de camp, analisi de laboratori i la corresponent interpretació i discussió dels resultats obtinguts. El segon objectiu específic es va basar en el desenvolupament i aplicació del DIS. L'àrea d'estudi seleccionada va ser la conca de la Serra de Rodes, un ambient típic Mediterràni inclòs en el Parc Natural del Cap de Creus, NE Espanya, el qual ha estat progressivament abandonat pels agricultors durant el segle passat. Actualment, els incendis forestals així com el canvi d'ús del sòl i especialment l'abandonament de terres són considerats els problemes ambientals més importants a l'àrea d'estudi (Dunjó et al., 2003). En primer lloc, es va realitzar l'estudi dels processos i causes de la degradació dels sòls a l'àrea d'interés. En base a aquest coneixement, es va dur a terme la identificació i selecció dels indicadors de desertificació més rellevants. Finalment, els indicadors de desertificació seleccionats a escala de conca, incloent l'erosió del sòl i l'escolament superficial, es van integrar en un model espaial de procés. Ja que el sòl és considerat el principal indicador dels processos d'erosió, segons la FAO/UNEP/UNESCO (1979), tant el paisatge original així com els dos escenaris d'ús del sòl desenvolupats, un centrat en el cas hipotétic del pas d'un incendi forestal, i l'altre un paisatge completament cultivat, poden ser ambients classificats sota baixa o moderada degradació. En comparació amb l'escenari original, els dos escenaris creats van revelar uns valors més elevats d'erosió i escolament superficial, i en particular l'escenari cultivat. Per tant, aquests dos hipotètic escenaris no semblen ser una alternativa sostenible vàlida als processos de degradació que es donen a l'àrea d'estudi. No obstant, un ampli ventall d'escenaris alternatius poden ser desenvolupats amb el DIS, tinguent en compte les polítiques d'especial interés per la regió de manera que puguin contribuir a determinar les conseqüències potencials de desertificació derivades d'aquestes polítiques aplicades en aquest escenari tan complexe espaialment. En conclusió, el model desenvolupat sembla ser un sistema força acurat per la identificació de riscs presents i futurs, així com per programar efectivament mesures per combatre la desertificació a escala de conca. No obstant, aquesta primera versió del model presenta varies limitacions i la necessitat de realitzar més recerca en cas de voler desenvolupar una versió futura i millor del DIS.

Stochastic modelling of rainfall from satellite data

Satellite-based rainfall monitoring is widely used for climatological studies because of its full global coverage but it is also of great importance for operational purposes especially in areas such as Africa where there is a lack of ground-based rainfall data. Satellite rainfall estimates have enormous potential benefits as input to hydrological and agricultural models because of their real time availability, low cost and full spatial coverage. One issue that needs to be addressed is the uncertainty on these estimates. This is particularly important in assessing the likely errors on the output from non-linear models (rainfall-runoff or crop yield) which make use of the rainfall estimates, aggregated over an area, as input. Correct assessment of the uncertainty on the rainfall is non-trivial as it must take account of • the difference in spatial support of the satellite information and independent data used for calibration • uncertainties on the independent calibration data • the non-Gaussian distribution of rainfall amount • the spatial intermittency of rainfall • the spatial correlation of the rainfall field This paper describes a method for estimating the uncertainty on satellite-based rainfall values taking account of these factors. The method involves firstly a stochastic calibration which completely describes the probability of rainfall occurrence and the pdf of rainfall amount for a given satellite value, and secondly the generation of ensemble of rainfall fields based on the stochastic calibration but with the correct spatial correlation structure within each ensemble member. This is achieved by the use of geostatistical sequential simulation. The ensemble generated in this way may be used to estimate uncertainty at larger spatial scales. A case study of daily rainfall monitoring in the Gambia, west Africa for the purpose of crop yield forecasting is presented to illustrate the method.

Evaluation of two hybrid metric-conceptual models, for simulating phosphorus transfer from agricultural land in the river enborne, a lowland UK catchment

Across Europe, elevated phosphorus (P) concentrations in lowland rivers have made them particularly susceptible to eutrophication. This is compounded in southern and central UK by increasing pressures on water resources, which may be further enhanced by the potential effects of climate change. The EU Water Framework Directive requires an integrated approach to water resources management at the catchment scale and highlights the need for modelling tools that can distinguish relative contributions from multiple nutrient sources and are consistent with the information content of the available data. Two such models are introduced and evaluated within a stochastic framework using daily flow and total phosphorus concentrations recorded in a clay catchment typical of many areas of the lowland UK. Both models disaggregate empirical annual load estimates, derived from land use data, as a function of surface/near surface runoff, generated using a simple conceptual rainfall-runoff model. Estimates of the daily load from agricultural land, together with those from baseflow and point sources, feed into an in-stream routing algorithm. The first model assumes constant concentrations in runoff via surface/near surface pathways and incorporates an additional P store in the river-bed sediments, depleted above a critical discharge, to explicitly simulate resuspension. The second model, which is simpler, simulates P concentrations as a function of surface/near surface runoff, thus emphasising the influence of non-point source loads during flow peaks and mixing of baseflow and point sources during low flows. The temporal consistency of parameter estimates and thus the suitability of each approach is assessed dynamically following a new approach based on Monte-Carlo analysis. (c) 2004 Elsevier B.V. All rights reserved.

The effect of soil phosphorus on particulate phosphorus in land runoff

Accumulation of surplus phosphorus (P) in the soil and the resulting increased transport of P in land runoff contribute to freshwater eutrophication. The effects of increasing soil P (19–194 mg Olsen-P (OP) kg−1) on the concentrations of particulate P (PP), and sorption properties (Qmax, k and EPCo) of suspended solids (SS) in overland flow from 15 unreplicated field plots established on a dispersive arable soil were measured over three monitoring periods under natural rainfall. Concentrations of PP in plot runoff increased linearly at a rate of 2.6 μg litre−1 per mg OP kg−1 of soil, but this rate was approximately 50% of the rate of increase in dissolved P (< 0.45 μm). Concentrations of SS in runoff were similar across all plots and contained a greater P sorption capacity (mean + 57%) than the soil because of enrichment with fine silt and clay (0.45–20 μm). As soil P increased, the P enrichment ratio of the SS declined exponentially, and the values of P saturation (Psat; 15–42%) and equilibrium P concentration (EPCo; 0.7–5.5 mg litre−1) in the SS fell within narrower ranges compared with the soils (6–74% and 0.1–10 mg litre−1, respectively). When OP was < 100 mg kg−1, Psat and EPCo values in the SS were smaller than those in the soil and vice-versa, suggesting that eroding particles from soils with both average and high P fertility would release P on entering the local (Rosemaund) stream. Increasing soil OP from average to high P fertility increased the P content of the SS by approximately 10%, but had no significant (P > 0.05) effect on the Psat, or EPCo, of the SS. Management options to reduce soil P status as a means of reducing P losses in land runoff and minimizing eutrophication risk may therefore have more limited effect than is currently assumed in catchment management.

Simulating current global river runoff with a global hydrological model: model revisions, validation, and sensitivity analysis

Global hydrological models (GHMs) model the land surface hydrologic dynamics of continental-scale river basins. Here we describe one such GHM, the Macro-scale - Probability-Distributed Moisture model.09 (Mac-PDM.09). The model has undergone a number of revisions since it was last applied in the hydrological literature. This paper serves to provide a detailed description of the latest version of the model. The main revisions include the following: (1) the ability for the model to be run for n repetitions, which provides more robust estimates of extreme hydrological behaviour, (2) the ability of the model to use a gridded field of coefficient of variation (CV) of daily rainfall for the stochastic disaggregation of monthly precipitation to daily precipitation, and (3) the model can now be forced with daily input climate data as well as monthly input climate data. We demonstrate the effects that each of these three revisions has on simulated runoff relative to before the revisions were applied. Importantly, we show that when Mac-PDM.09 is forced with monthly input data, it results in a negative runoff bias relative to when daily forcings are applied, for regions of the globe where the day-to-day variability in relative humidity is high. The runoff bias can be up to - 80% for a small selection of catchments but the absolute magnitude of the bias may be small. As such, we recommend future applications of Mac-PDM.09 that use monthly climate forcings acknowledge the bias as a limitation of the model. The performance of Mac-PDM.09 is evaluated by validating simulated runoff against observed runoff for 50 catchments. We also present a sensitivity analysis that demonstrates that simulated runoff is considerably more sensitive to method of PE calculation than to perturbations in soil moisture and field capacity parameters.

Hillslope scale surface runoff, sediment, and nutrient losses associated with tramline wheelings

Research on arable sandy loam and silty clay loam soils on 4° slopes in England has shown that tramlines (i.e. the unseeded wheeling areas used to facilitate spraying operations in cereal crops) can represent the most important pathway for phosphorus and sediment loss from moderately sloping fields. Detailed monitoring over the October–March period in winters 2005–2006 and 2006–2007 included event-based sampling of surface runoff, suspended and particulate sediment, and dissolved and particulate phosphorus from hillslope segments (each ∼300–800 m2) established in a randomized block design with four replicates of each treatment at each of two sites on lighter and heavier soils. Experimental treatments assessed losses from the cropped area without tramlines, and from the uncropped tramline area, and were compared to losses from tramlines which had been disrupted once in the autumn with a shallow tine. On the lighter soil, the effects of removal or shallow incorporation of straw residues was also determined. Research on both sandy and silty clay loam soils across two winters showed that tramline wheelings represented the dominant pathway for surface runoff and transport of sediment, phosphorus and nitrogen from cereal crops on moderate slopes. Results indicated 5·5–15·8% of rainfall lost as runoff, and losses of 0·8–2·9 kg TP ha−1 and 0·3–4·8 t ha−1 sediment in tramline treatments, compared to only 0·2–1·7% rainfall lost as runoff, and losses of 0·0–0·2 kg TP ha−1 and 0·003–0·3 t ha−1 sediment from treatments without tramlines or those where tramlines had been disrupted. The novel shallow disruption of tramline wheelings using a tine once following the autumn spray operation consistently and dramatically reduced (p < 0·001) surface runoff and loads of sediment, total nitrogen and total phosphorus to levels similar to those measured in cropped areas between tramlines. Results suggest that options for managing tramline wheelings warrant further refinement and evaluation with a view to incorporating them into spatially-targeted farm-level management planning using national or catchment-based agri-environment policy instruments aimed at reducing diffuse pollution from land to surface water systems. Copyright © 2010 John Wiley & Sons, Ltd.