Hydrological system analysis and modelling of the Kara River basin (West Africa) using a lumped metric conceptual model

This paper discusses the analysis and modelling of the hydrological system of the basin of the Kara River, a transboundary river in Togo and Benin, as a necessary step towards sustainable water resources management. The methodological approach integrates the use of discharge parameters, flow duration curves and the lumped conceptual model IHACRES. A Sobol sensitivity analysis is performed and the model is calibrated by applying the shuffled complex evolution algorithm. Results show that discharge generation in three nested catchments of the basin is affected by landscape physical characteristics. The IHACRES model adequately simulates the rainfall–runoff dynamics in the basin with a mean modified Nash-Sutcliffe efficiency measure of 0.6. Modelling results indicate that parameters controlling rainfall transformation to effective rainfall are more sensitive than those routing the streamflow. This study provides insights into understanding the catchment’s hydrological system. Nevertheless, further investigations are required to better understand detailed runoff generation processes.

EDITOR M.C. Acreman; ASSOCIATE EDITOR N Verhoest     

Geophysical flows under location uncertainty,Part I Random transport and general models

A stochastic flow representation is considered with the Eulerian velocity decomposed between a smooth large scale component and a rough small-scale turbulent component. The latter is specified as a random field uncorrelated in time. Subsequently, the material derivative is modified and leads to a stochastic version of the material derivative to include a drift correction, an inhomogeneous and anisotropic diffusion, and a multiplicative noise. As derived, this stochastic transport exhibits a remarkable energy conservation property for any realizations. As demonstrated, this pivotal operator further provides elegant means to derive stochastic formulations of classical representations of geophysical flow dynamics.     

Geophysical flows under location uncertainty,Part II Quasi-geostrophy and efficient ensemble spreading

Models under location uncertainty are derived assuming that a component of the velocity is uncorrelated in time. The material derivative is accordingly modified to include an advection correction, inhomogeneous and anisotropic diffusion terms and a multiplicative noise contribution. In this paper, simplified geophysical dynamics are derived from a Boussinesq model under location uncertainty. Invoking usual scaling approximations and a moderate influence of the subgrid terms, stochastic formulations are obtained for the stratified Quasi-Geostrophy and the Surface Quasi-Geostrophy models. Based on numerical simulations, benefits of the proposed stochastic formalism are demonstrated. A single realization of models under location uncertainty can restore small-scale structures. An ensemble of realizations further helps to assess model error prediction and outperforms perturbed deterministic models by one order of magnitude. Such a high uncertainty quantification skill is of primary interests for assimilation ensemble methods. MATLAB^® code examples are available online.     

Enhanced identification of a hydrologic model using streamflow and satellite water storage data: A multicriteria sensitivity analysis and optimization approach

Hydrologic model development and calibration have continued in most cases to focus only on accurately reproducing streamflows. However, complex models, for example, the so‐called physically based models, possess large degrees of freedom that, if not constrained properly, may lead to poor model performance when used for prediction. We argue that constraining a model to represent streamflow, which is an integrated resultant of many factors across the watershed, is necessary but by no means sufficient to develop a high‐fidelity model. To address this problem, we develop a framework to utilize the Gravity Recovery and Climate Experiment's (GRACE) total water storage anomaly data as a supplement to streamflows for model calibration, in a multiobjective setting. The VARS method (Variogram Analysis of Response Surfaces) for global sensitivity analysis is used to understand the model behaviour with respect to streamflow and GRACE data, and the BORG multiobjective optimization method is applied for model calibration. Two subbasins of the Saskatchewan River Basin in Western Canada are used as a case study. Results show that the developed framework is superior to the conventional approach of calibration only to streamflows, even when multiple streamflow‐based error functions are simultaneously minimized. It is shown that a range of (possibly false) system trajectories in state variable space can lead to similar (acceptable) model responses. This observation has significant implications for land‐surface and hydrologic model development and, if not addressed properly, may undermine the credibility of the model in prediction. The framework effectively constrains the model behaviour (by constraining posterior parameter space) and results in more credible representation of hydrology across the watershed.     

Multicriteria sensitivity analysis as a diagnostic tool for understanding model behaviour and characterizing model uncertainty

Complex hydrological models are being increasingly used nowadays for many purposes such as studying the impact of climate and land‐use change on water resources. However, building a high‐fidelity model, particularly at large scales, remains a challenging task, due to complexities in model functioning and behaviour and uncertainties in model structure, parameterization, and data. Global sensitivity analysis (GSA), which characterizes how the variation in the model response is attributed to variations in its input factors (e.g., parameters and forcing data), provides an opportunity to enhance the development and application of these complex models. In this paper, we advocate using GSA as an integral part of the modelling process by discussing its capabilities as a tool for diagnosing model structure and detecting potential defects, identifying influential factors, characterizing uncertainty, and selecting calibration parameters. Accordingly, we conduct a comprehensive GSA of a complex land surface–hydrology model, Modélisation Environmentale–Surface et Hydrologie (MESH), which combines the Canadian land surface scheme with a hydrological routing component, WATROUTE. Various GSA experiments are carried out using a new technique, called Variogram Analysis of Response Surfaces, for alternative hydroclimatic conditions in Canada using multiple criteria, various model configurations, and a full set of model parameters. Results from this study reveal that, in addition to different hydroclimatic conditions and SA criteria, model configurations can also have a major impact on the assessment of sensitivity. GSA can identify aspects of the model internal functioning that are counter‐intuitive and thus help the modeller to diagnose possible model deficiencies and make recommendations for improving development and application of the model. As a specific outcome of this work, a list of the most influential parameters for the MESH model is developed. This list, along with some specific recommendations, is expected to assist the wide community of MESH and Canadian land surface scheme users, to enhance their modelling applications.     

Uncertainty and sensitivity in a bank stability model: implications for estimating phosphorus loading

Eutrophication of aquatic ecosystems is one of the most pressing water quality concerns in the United States and around the world. Bank erosion has been largely overlooked as a source of nutrient loading, despite field studies demonstrating that this source can account for the majority of the total phosphorus load in a watershed. Substantial effort has been made to develop mechanistic models to predict bank erosion and instability in stream systems; however, these models do not account for inherent natural variability in input values. To quantify the impacts of this omission, uncertainty and sensitivity analyses were performed on the Bank Stability and Toe Erosion Model (BSTEM), a mechanistic model developed by the US Department of Agriculture – Agricultural Research Service (USDA‐ARS) that simulates both mass wasting and fluvial erosion of streambanks. Generally, bank height, soil cohesion, and plant species were found to be most influential in determining stability of clay (cohesive) banks. In addition to these three inputs, groundwater elevation, stream stage, and bank angle were also identified as important in sand (non‐cohesive) banks. Slope and bank height are the dominant variables in fluvial erosion modeling, while erodibility and critical shear stress had low sensitivity indices; however, these indices do not reflect the importance of critical shear stress in determining the timing of erosion events. These results identify important variables that should be the focus of data collection efforts while also indicating which less influential variables may be set to assumed values. In addition, a probabilistic Monte‐Carlo modeling approach was applied to data from a watershed‐scale sediment and phosphorus loading study on the Missisquoi River, Vermont to quantify uncertainty associated with these published results. While our estimates aligned well with previous deterministic modeling results, the uncertainty associated with these predictions suggests that they should be considered order of magnitude estimates only. Copyright © 2016 John Wiley & Sons, Ltd.     

土壤水环境中有机污染物生物处理过程模型研究

生物处理是一种经济有效处理土壤水环境中有机污染物的手段,本文在研究土地生物处理过程的基础上,建立了综合描述有机污染物在土壤－水－微生物系统中扩散、吸附／解吸、屏蔽和生物降解过程的数学模型。为确定模型中各参数在模型计算中的作用和相对重要性,进行了参数灵敏度分析,预计数学模型可以定量预测有机污染物进行土地生物处理所需的要时间和程度,为构建土地生物处理工程提供参考。     

全球温升1.5℃和2.0℃情景下中国降雨诱发地质灾害危险性和人口暴露度评估研究

气候变化情景下极端降水事件的频次和强度预估呈增加趋势,这会导致全球部分地区极端降雨诱发地质灾害风险的增加。本文基于中国降雨诱发地质灾害易发性模型和不同地貌分区的累积事件降雨量-降雨历时阈值曲线,采用最新的CMIP6全球气候模式多模式集合结果,基于全球温升目标情景的视角,从地质灾害空间易发性和发生频次两方面,探讨温升情景下中国地质灾害危险性的可能变化及其对暴露人口的潜在影响。结果表明,CMIP6多模式集合预估的多年平均降水在温升1.5℃和2.0℃情景下相比基准时期可能增加5.4%～9.5%,导致中等至极高地质灾害易发区范围预估增加0.33%～0.74%,由于预估的极端降水事件增加,地质灾害发生频次预估增加7.0%～11.2%,进一步综合未来人口空间分布,潜在地质灾害暴露人口可能增加6.20亿人次（18.90%）和4.26亿人次（12.97%）。各地貌分区未来情景下地质灾害危险性预估增加且存在显著的空间异质性,温升2.0℃情景下中等至极高易发性范围相比基准时期增加0.71%～1.28%,地质灾害发生频次预估增加1.2%～15.6%,其中,青藏高原区地质灾害危险性增加最明显。综合考虑未来人口...     

Geomorphic effectiveness: a linear concept in a non‐linear world

Geomorphic effectiveness has been an influential concept in geomorphology since its introduction by Reds Wolman and John Miller in 1960. It provided a much needed framework to assess the significance of an event by comparing event magnitude to the resultant geomorphic effects. Initially, this concept was applied primarily in river channels, under the linear assumption that geomorphic responses to similarly sized flood events will be consistent. Numerous authors have since attempted to quantify a direct, proportional relationship between event magnitude and different forms of geomorphic response in a variety of geomorphic settings. In doing so, these investigations applied an array of metrics that were difficult to compare across different spatiotemporal scales, and physiographic and geomorphic environments. Critically, the emergence of other geomorphic concepts such as sensitivity, connectivity, thresholds, and recovery has shown that relationships between causes (events) and geomorphic effects (responses) are often complex and non‐linear. This paper disentangles the complex historical development of the geomorphic effectiveness concept and reviews the utility of various metrics for quantifying effectiveness. We propose that total energy (joules) is the most appropriate metric to use for quantifying the magnitude of disturbance events (cause) and volumetric sediment flux associated with landform modification is the most appropriate metric for quantifying geomorphic effects. While both metrics are difficult to quantify, they are the only ones which facilitate comparison across a range of spatiotemporal scales (comparability) in a variety of geomorphic environments (flexibility). The geomorphic effectiveness concept can continue to be useful provided that geomorphologists use flexible and comparable metrics. Today, geomorphologists are better prepared to consider the influence of non‐linear processes on determinations of geomorphic effectiveness, allowing investigators to not only determine if a disturbance event was effective but also to explain why or why not. Copyright © 2016 John Wiley & Sons, Ltd.     

Sensitivity analysis of EUROSEM using Monte Carlo simulation I: hydrological,soil and vegetation parameters

Knowledge about model uncertainty is essential for erosion modelling and provides important information when it comes to parameterizing models. In this paper a sensitivity analysis of the European soil erosion model (EUROSEM) is carried out using Monte Carlo simulation, suitable for complex non‐linear models, using time‐dependent driving variables. The analysis revealed some important characteristics of the model. The variability of the static output parameters was generally high, with the hydrologic parameters being the most important ones, especially saturated hydraulic conductivity and net capillary drive followed by the percentage basal area for the hydrological and vegetation parameters and detachability and cohesion for the soil erosion parameters. Overall, sensitivity to vegetation parameters was insignificant. The coefficient of variation for the sedigraph was higher than for the hydrograph, especially from the beginning of the rainstorm and up to the peak, and may explain difficulties encountered when trying to match simulated hydrographs and sedigraphs with observed ones. The findings from this Monte Carlo simulation calls for improved within‐storm modelling of erosion processes in EUROSEM. Information about model uncertainty will be incorporated in a new EUROSEM user interface. Copyright © 2000 John Wiley & Sons, Ltd.