IHMS—Integrated Hydrological Modelling System. Part 1. Hydrological processes and general structure

A newly Integrated Hydrological Modelling System (IHMS) has been developed to study the impact of changes in climate, land use and water management on groundwater and seawater intrusion (SWI) into coastal areas. The system represents the combination of three models, which can, if required, be run separately. It has been designed to assess the combined impact of climate, land use and groundwater abstraction changes on river, drainage and groundwater flows, groundwater levels and, where appropriate, SWI. The approach is interdisciplinary and reflects an integrated water management approach. The system comprises three packages: the Distributed Catchment Scale Model (DiCaSM), MODFLOW (96 and 2000) and SWI models. In addition to estimating all water balance components, DiCaSM, produces the recharge data that are used as input to the groundwater flow model of the US Geological Survey, MODFLOW. The latter subsequently generates the head distribution and groundwater flows that are used as input to the SWI model, SWI. Thus, any changes in land use, rainfall, water management, abstraction, etc. at the surface are first handled by DiCaSM, then by MODFLOW and finally by the SWI. The three models operate at different spatial and temporal scales and a facility (interface utilities between models) to aggregate/disaggregate input/output data to meet a desired spatial and temporal scale was developed allowing smooth and easy communication between the three models. As MODFLOW and SWI are published and in the public domain, this article focuses on DiCaSM, the newly developed unsaturated zone DiCaSM and equally important the interfacing utilities between the three models. DiCaSM simulates a number of hydrological processes: rainfall interception, evapotranspiration, surface runoff, infiltration, soil water movement in the root zone, plant water uptake, crop growth, stream flow and groundwater recharge. Input requirements include distributed data sets of rainfall, land use, soil types and digital terrain; climate data input can be either distributed or non‐distributed. The model produces distributed and time series output of all water balance components including potential evapotranspiration, actual evapotranspiration, rainfall interception, infiltration, plant water uptake, transpiration, soil water content, soil moisture (SM) deficit, groundwater recharge rate, stream flow and surface runoff. This article focuses on details of the hydrological processes and the various equations used in DiCaSM, as well as the nature of the interface to the MODFLOW and SWI models. Furthermore, the results of preliminary tests of DiCaSM are reported; these include tests related to the ability of the model to predict the SM content of surface and subsurface soil layers, as well as groundwater levels. The latter demonstrates how the groundwater recharge calculated from DiCaSM can be used as input into the groundwater model MODFLOW using aggregation and disaggregation algorithms (built into the interface utility). SWI has also been run successfully with hypothetical examples and was able to reproduce the results of some of the original examples of Bakker and Schaars ( 2005 ). In the subsequent articles, the results of applications to different catchments will be reported. Copyright © 2010 John Wiley & Sons, Ltd.     

近10年来中国农村居民点用地时空特征

利用20世纪90年代中国1:10万土地利用动态变化数据,利用单元自动机和人工神经网络模型对中国农村居民点用地进行了区划。在此基础上,研究了90年代中国农村居民点用地动态变化时空格局。研究表明,农村居民点用地受区域位置、经济发展和国家政策的影响。90年代前5年,受沿海开放政策和经济高速发展的影响,长江三角洲、珠江三角洲和华北平原农村居民点用地扩展较快,中西部地区农村居民点用地扩展较慢;90年代后5年,在沿海地区经济发展速度放慢和国家实施耕地资源保护条例的情况下,沿海地区农村居民点用地扩张受到抑制,西部地区扩展加快。90年代中国经济体制处于由计划经济向市场经济转轨的过程中,农村居民点用地的动态变化格局受到经济发展和土地利用政策的深刻影响。     

Sub‐grid scale parameterization of hillslope runoff and erosion processes for catchment‐scale models of semi‐arid landscapes

The processes of hillslope runoff and erosion are typically represented at coarse spatial resolution in catchment‐scale models due to computational limitations. Such representation typically fails to incorporate the important effects of topographic heterogeneity on runoff generation, overland flow, and soil erosion. These limitations currently undermine the application of distributed catchment models to understand the importance of thresholds and connectivity on hillslope and catchment‐scale runoff and erosion, particularly in semi‐arid environments. This paper presents a method for incorporating high‐resolution topographic data to improve sub‐grid scale parameterization of hillslope overland flow and erosion models. Results derived from simulations conducted using a kinematic wave overland flow model at 0.5 m spatial resolution are used to parameterize the depth–discharge relationship in the overland flow model when applied at 16 m resolution. The high‐resolution simulations are also used to derive a more realistic parameterization of excess flow shear stress for use in the 16 m resolution erosion model. Incorporating the sub‐grid scale parameterization in the coarse‐resolution model (16 m) leads to improved predictions of overland flow and erosion when evaluated using results derived from high‐resolution (0.5 m) model simulations. The improvement in performance is observed for a range of event magnitudes and is most notable for erosion estimates due to the non‐linear dependency between the rates of erosion and overland flow. Copyright © 2013 John Wiley & Sons, Ltd.     

A model for distributed GIUH‐based flow routing on natural and anthropogenic hillslopes

Attempts to reduce the number of parameters in distributed rainfall–runoff models have not yet resulted in a model that is accurate for both natural and anthropogenic hillslopes. We take on the challenge by proposing a distributed model for overland flow and channel flow based on a combination of a linear response time distribution and the hillslope geomorphologic instantaneous unit hydrograph (GIUH), which can be calculated with only a digital elevation model and a map with field boundaries and channel network as input. The spatial domain is subdivided into representative elementary hillslopes (REHs) for each of which we define geometric and flow velocity parameters and compute the GIUH. The catchment GIUH is given by the sum of all REH responses. While most distributed models only perform well on natural hillslopes, the advantage of our approach is that it can also be applied to modified hillslopes with for example a rectangular drainage network and terrace cultivation. Tests show that the REH‐GIUH approach performs better than classical routing functions (exponential and gamma). Simulations of four virtual hillslopes suggest that peak flow at the catchment outlet is directly related to drainage density. By combining the distributed flow routing model with a lumped‐parameter infiltration model, we were also able to demonstrate that terrace cultivation delays the response time and reduces peak flow in comparison to the same hillslope, but with a natural stream network. The REH‐GIUH approach is a first step in the process of coupling distributed hydrological models to erosion and water quality models at the REH (associated with agricultural management) and at the catchment scale (associated with the evaluation of the environmental impact of human activities). It furthermore provides a basis for the development of models for large catchments and urban or peri‐urban catchments. Copyright © 2013 John Wiley & Sons, Ltd.     

Field investigation and modelling of coupled stream discharge and shallow water‐table dynamics in a small Mediterranean catchment (Sardinia)

In the semi‐arid Mediterranean environment, the rainfall–runoff relationships are complex because of the markedly irregular patterns in rainfall, the seasonal mismatch between evaporation and rainfall, and the spatial heterogeneity in landscape properties. Watersheds often display considerable non‐linear threshold behavior, which still make runoff generation an open research question. Our objectives in this context were: to identify the primary processes of runoff generation in a small natural catchment; to test whether a physically based model, which takes into consideration only the primary processes, is able to predict spatially distributed water‐table and stream discharge dynamics; and to use the hydrological model to increase our understanding of runoff generation mechanisms. The observed seasonal dynamics of soil moisture, water‐table depth, and stream discharge indicated that Hortonian overland‐flow was negligible and the main mechanism of runoff generation was saturated subsurface‐flow. This gives rise to base‐flow, controls the formation of the saturated areas, and contributes to storm‐flow together with saturation overland‐flow. The distributed model, with a 1D scheme for the kinematic surface‐flow, a 2D sub‐horizontal scheme for the saturated subsurface‐flow, and ignoring the unsaturated flow, performed efficiently in years when runoff volume was high and medium, although there was a smoothing effect on the observed water‐table. In dry years, small errors greatly reduced the efficiency of the model. The hydrological model has allowed to relate the runoff generation mechanisms with the land‐use. The forested hillslopes, where the calibrated soil conductivity was high, were never saturated, except at the foot of the slopes, where exfiltration of saturated subsurface‐flow contributed to storm‐flow. Saturation overland‐flow was only found near the streams, except when there were storm‐flow peaks, when it also occurred on hillslopes used for pasture, where soil conductivity was low. The bedrock–soil percolation, simulated by a threshold mechanism, further increased the non‐linearity of the rainfall–runoff processes. Copyright © 2013 John Wiley & Sons, Ltd.     

Simulation of stream flow components in a mountainous catchment in northern Thailand with SWAT,using the ANSELM calibration approach

Highland agriculture is intensifying rapidly in South‐East Asia, leading to alarmingly high applications of agrochemicals. Understanding the fate of these contaminants requires carefully planned monitoring programmes and, in most cases, accurate simulation of hydrological pathways into and through water bodies. We simulate run‐off in a steep mountainous catchment in tropical South‐East Asia. To overcome calibration difficulties related to the mountainous topography, we introduce a new calibration method, named A Nash–Sutcliffe Efficiency Likelihood Match (ANSELM), that allows the assignment of optimal parameters to different hydrological response units in simulations of stream discharge with the Soil and Water Assessment Tool (SWAT) hydrological model. ANSELM performed better than the Parasol calibration tool built into SWAT in terms of model efficiency and computation time. In our simulation, the most sensitive model parameters were those related to base flow generation, surface run‐off generation, flow routing and soil moisture change. The coupling of SWAT with ANSELM yielded reasonable simulations of both wet‐season and dry‐season storm hydrographs. Nash–Sutcliffe model efficiencies for daily stream flow during two validation years were 0.77 and 0.87. These values are in the upper range or even higher than those reported for other SWAT model applications in temperate or tropical regions. The different flow components were realistically simulated by SWAT, and showed a similar behaviour in all the study years, despite inter‐annual climatic differences. The realistic partitioning of total stream flow into its contributing components will be an important factor for using this hydrological model to simulate solute transport in the future. Copyright © 2014 John Wiley & Sons, Ltd.     

Inter‐comparison of hydro‐climatic regimes across northern catchments: synchronicity,resistance and resilience

The higher mid‐latitudes of the Northern Hemisphere are particularly sensitive to climate change as small differences in temperature determine frozen ground status, precipitation phase, and the magnitude and timing of snow accumulation and melt. An international inter‐catchment comparison program, North‐Watch, seeks to improve our understanding of the sensitivity of northern catchments to climate change by examining their hydrological and biogeochemical responses. The catchments are located in Sweden (Krycklan), Scotland (Mharcaidh, Girnock and Strontian), the United States (Sleepers River, Hubbard Brook and HJ Andrews) and Canada (Catamaran, Dorset and Wolf Creek). This briefing presents the initial stage of the North‐Watch program, which focuses on how these catchments collect, store and release water and identify ‘types’ of hydro‐climatic catchment response. At most sites, a 10‐year data of daily precipitation, discharge and temperature were compiled and evaporation and storage were calculated. Inter‐annual and seasonal patterns of hydrological processes were assessed via normalized fluxes and standard flow metrics. At the annual‐scale, relations between temperature, precipitation and discharge were compared, highlighting the role of seasonality, wetness and snow/frozen ground. The seasonal pattern and synchronicity of fluxes at the monthly scale provided insight into system memory and the role of storage. We identified types of catchments that rapidly translate precipitation into runoff and others that more readily store water for delayed release. Synchronicity and variance of rainfall–runoff patterns were characterized by the coefficient of variation (cv) of monthly fluxes and correlation coefficients. Principal component analysis (PCA) revealed clustering among like catchments in terms of functioning, largely controlled by two components that (i) reflect temperature and precipitation gradients and the correlation of monthly precipitation and discharge and (ii) the seasonality of precipitation and storage. By advancing the ecological concepts of resistance and resilience for catchment functioning, results provided a conceptual framework for understanding susceptibility to hydrological change across northern catchments. Copyright © 2010 John Wiley & Sons, Ltd.     

Modelling rating curves using remotely sensed LiDAR data

Accurate stream discharge measurements are important for many hydrological studies. In remote locations, however, it is often difficult to obtain stream flow information because of the difficulty in making the discharge measurements necessary to define stage‐discharge relationships (rating curves). This study investigates the feasibility of defining rating curves by using a fluid mechanics‐based model constrained with topographic data from an airborne LiDAR scanning. The study was carried out for an 8m‐wide channel in the boreal landscape of northern Sweden. LiDAR data were used to define channel geometry above a low flow water surface along the 90‐m surveyed reach. The channel topography below the water surface was estimated using the simple assumption of a flat streambed. The roughness for the modelled reach was back calculated from a single measurment of discharge. The topographic and roughness information was then used to model a rating curve. To isolate the potential influence of the flat bed assumption, a ‘hybrid model’ rating curve was developed on the basis of data combined from the LiDAR scan and a detailed ground survey. Whereas this hybrid model rating curve was in agreement with the direct measurements of discharge, the LiDAR model rating curve was equally in agreement with the medium and high flow measurements based on confidence intervals calculated from the direct measurements. The discrepancy between the LiDAR model rating curve and the low flow measurements was likely due to reduced roughness associated with unresolved submerged bed topography. Scanning during periods of low flow can help minimize this deficiency. These results suggest that combined ground surveys and LiDAR scans or multifrequency LiDAR scans that see ‘below’ the water surface (bathymetric LiDAR) could be useful in generating data needed to run such a fluid mechanics‐based model. This opens a realm of possibility to remotely sense and monitor stream flows in channels in remote locations. Copyright © 2012 John Wiley & Sons, Ltd.     

Potential of time‐lapse photography of snow for hydrological purposes at the small catchment scale

Time‐lapse photography provides an attractive source of information about snow cover characteristics, especially at the small catchment scale. The objective of this study was to design and test a monitoring system, which allows multi‐resolution observations of snow cover characteristics. The main aim was to simultaneously investigate the spatio‐temporal patterns of snow cover, snow depth and snowfall interception in the area very close to the camera, and the spatio‐temporal patterns of snow cover in the far range. The multi‐resolution design was tested at three sites in the eastern part of the Austrian Alps (Hochschwab‐Rax region). Digital photographs were taken at hourly time steps between 6:00 and 18:00 in the period November, 2004 to December, 2006. The results showed that the time‐lapse photography allows effective mapping of the snow depths at high temporal resolution in the region close to the digital camera at many snow stake locations. It is possible to process a large number of photos by using an automatic procedure for accurate snow depth readings. The digital photographs can also be used to infer the settling characteristics of the snow pack and snow interception during the day. Although it is not possible to directly estimate the snow interception mass, the photos may indeed give very useful information on the snow processes on and beneath the forest canopy. The main advantage of using time‐lapse photography in the far range of the digital camera is to observe the spatio‐temporal patterns of snow cover over different landscape configurations. The results illustrate that digital photographs can be very useful for parameterising processes such as sloughing on steep slopes, snow deposition in gullies and snow erosion on mountain ridges in a distributed snow model. Copyright © 2011 John Wiley & Sons, Ltd.