Multi‐criterial validation of TOPMODEL in a mountainous catchment

The need for powerful validation methods for hydrological models including the evaluation of internal stages and spatially distributed simulations has often been emphasized. In this study a multi‐criterial validation scheme was used for validation of TOPMODEL, a conceptual semi‐distributed rainfall–runoff model. The objective was to test TOPMODEL's capability of adequately representing dominant hydrological processes by simple conceptual approaches. Validation methods differed in the type of data used, in their target and in mode. The model was applied in the humid and mountainous Brugga catchment (40 km²) in south‐west Germany. It was calibrated by a Monte Carlo method based on hourly runoff data. Additional information for validation was derived from a recession analysis, hydrograph separation with environmental tracers and from field surveys, including the mapping of saturated areas. Although runoff simulations were satisfying, inadequacies of the model structure compared with the real situation with regard to hydrological processes in the study area were found. These belong mainly to the concept of variable contributing areas for saturation excess overland flow and their dynamics, which were overestimated by the model. The simple TOPMODEL approach of two flow components was found to be insufficient. The multi‐criterial validation scheme enables not only to demonstrate limitations with regard to process representation, but also to specify where and why these limitations occur. It may serve as a valuable tool for the development of physically sound model modifications. Copyright © 1999 John Wiley & Sons, Ltd.     

A distributed monthly hydrological model for integrating spatial variations of basin topography and rainfall

Hydrological models at a monthly time‐scale are important tools for hydrological analysis, such as in impact assessment of climate change and regional water resources planning. Traditionally, monthly models adopt a conceptual, lumped‐parameter approach and cannot account for spatial variations of basin characteristics and climatic inputs. A large requirement for data often severely limits the utility of physically based, distributed‐parameter models. Based on the variable‐source‐area concept, we considered basin topography and rainfall to be two major factors whose spatial variations play a dominant role in runoff generation and developed a monthly model that is able to account for their influences in the spatial and temporal dynamics of water balance. As a hybrid of the Xinanjiang model and TOPMODEL, the new model is constructed by innovatively making use of the highly acclaimed simulation techniques in the two existing models. A major contribution of this model development study is to adopt the technique of implicit representation of soil moisture characteristics in the Xinanjiang model and use the TOPMODEL concept to integrate terrain variations into runoff simulation. Specifically, the TOPMODEL topographic index ln(a/tanβ) is converted into an index of relative difficulty in runoff generation (IRDG) and then the cumulative frequency distribution of IRDG is used to substitute the parabolic curve, which represents the spatial variation of soil storage capacity in the Xinanjiang model. Digital elevation model data play a key role in the modelling procedures on a geographical information system platform, including basin segmentation, estimation of rainfall for each sub‐basin and computation of terrain characteristics. Other monthly data for model calibration and validation are rainfall, pan evaporation and runoff. The new model has only three parameters to be estimated, i.e. watershed‐average field capacity WM, pan coefficient η and runoff generation coefficient α. Sensitivity analysis demonstrates that runoff is least sensitive to WM and, therefore, it can be determined by a prior estimation based on the climate and soil properties of the study basin. The other two parameters can be determined using optimization methods. Model testing was carried out in a number of nested sub‐basins of two watersheds (Yuanjiang River and Dongjiang River) in the humid region in central and southern China. Simulation results show that the model is capable of describing spatial and temporal variations of water balance components, including soil moisture content, evapotranspiration and runoff, over the watershed. With a minimal requirement for input data and parameterization, this terrain‐based distributed model is a valuable contribution to the ever‐advancing technology of hydrological modelling. Copyright © 2006 John Wiley & Sons, Ltd.     

End member and Bayesian mixing models consistently indicate near-surface flowpath dominance in a pristine humid tropical rainforest

The impacts of forest conversion on runoff generation in the tropics have received much interest, but scientific progress is still hampered by challenging fieldwork conditions and limited knowledge about runoff mechanisms. Here, we assessed the runoff generation, flow paths and water source dynamics of a pristine rainforest catchment in Costa Rica using end member mixing analysis (EMMA) and a Bayesian mixing model (MixSIAR). Geochemical tracer data collected over a 4-week field campaign were combined with tritium data used to assess potential deeper groundwater flow pathways to the perennial stream. The streamflow composition was best captured using three end-members, namely throughfall, shallow (5–15 cm) and deeper (15–50 cm) soil water. We estimated the end-member contributions to the main stream and two tributaries using the two mixing approaches and found good agreement between results obtained from EMMA and MixSIAR. The system was overwhelmingly dominated by near-surface sources, with little evidence for deeper and older groundwater as tritium-derived baseflow mean transit time was between 2.0 and 4.4 years. The shallow soil flow pathway dominated streamflow contributions in the main stream (median 39% and 49% based on EMMA and MixSIAR, respectively), followed by the deeper soil (32% and 31%) and throughfall (25% and 19%). The two tributaries had even greater shallow soil water contributions relative to the main stream (83% and 74% for tributary A and 42% and 63% for tributary B). Tributary B had no detectable deep soil water contribution, reflecting the morphology of the hillslope (steeper slopes, shallower soils and lower vegetation density compared to hillslope A). Despite the short sampling campaign and associated uncertainties, this study allowed to thoroughly assess runoff generation mechanisms in a humid tropical catchment. Our results also provide a first comparison of two increasingly used mixing models and suggest that EMMA and MixSIAR yield comparable estimates of water source partitioning in this tropical, volcanic rainforest environment.     

Spatially distributed tracer‐aided modelling to explore water and isotope transport,storage and mixing in a pristine,humid tropical catchment

Rapidly transforming headwater catchments in the humid tropics provide important resources for drinking water, irrigation, hydropower, and ecosystem connectivity. However, such resources for downstream use remain unstudied. To improve understanding of the behaviour and influence of pristine rainforests on water and tracer fluxes, we adapted the relatively parsimonious, spatially distributed tracer‐aided rainfall–runoff (STARR) model using event‐based stable isotope data for the 3.2‐km² San Lorencito catchment in Costa Rica. STARR was used to simulate rainforest interception of water and stable isotopes, which showed a significant isotopic enrichment in throughfall compared with gross rainfall. Acceptable concurrent simulations of discharge (Kling–Gupta efficiency [KGE] ~0.8) and stable isotopes in stream water (KGE ~0.6) at high spatial (10 m) and temporal (hourly) resolution indicated a rapidly responding system. Around 90% of average annual streamflow (2,099 mm) was composed of quick, near‐surface runoff components, whereas only ~10% originated from groundwater in deeper layers. Simulated actual evapotranspiration (ET) from interception and soil storage were low (~420 mm/year) due to high relative humidity (average 96%) and cloud cover limiting radiation inputs. Modelling suggested a highly variable groundwater storage (~10 to 500 mm) in this steep, fractured volcanic catchment that sustains dry season baseflows. This groundwater is concentrated in riparian areas as an alluvial–colluvial aquifer connected to the stream. This was supported by rainfall–runoff isotope simulations, showing a “flashy” stream response to rainfall with only a moderate damping effect and a constant isotope signature from deeper groundwater (~400‐mm additional mixing volume) during baseflow. The work serves as a first attempt to apply a spatially distributed tracer‐aided model to a tropical rainforest environment exploring the hydrological functioning of a steep, fractured‐volcanic catchment. We also highlight limitations and propose a roadmap for future data collection and spatially distributed tracer‐aided model development in tropical headwater catchments.     

Coupling the Xinanjiang model to a kinematic flow model based on digital drainage networks for flood forecasting

The Xinanjiang model, which is a conceptual rainfall‐runoff model and has been successfully and widely applied in humid and semi‐humid regions in China, is coupled by the physically based kinematic wave method based on a digital drainage network. The kinematic wave Xinanjiang model (KWXAJ) uses topography and land use data to simulate runoff and overland flow routing. For the modelling, the catchment is subdivided into numerous hillslopes and consists of a raster grid of flow vectors that define the water flow directions. The Xinanjiang model simulates the runoff yield in each grid cell, and the kinematic wave approach is then applied to a ranked raster network. The grid‐based rainfall‐runoff model was applied to simulate basin‐scale water discharge from an 805‐km² catchment of the Huaihe River, China. Rainfall and discharge records were available for the years 1984, 1985, 1987, 1998 and 1999. Eight flood events were used to calibrate the model's parameters and three other flood events were used to validate the grid‐based rainfall‐runoff model. A Manning's roughness via a linear flood depth relationship was suggested in this paper for improving flood forecasting. The calibration and validation results show that this model works well. A sensitivity analysis was further performed to evaluate the variation of topography (hillslopes) and land use parameters on catchment discharge. Copyright © 2009 John Wiley & Sons, Ltd.     

TOPMODEL: A critique

TOPMODEL (a TOPography based hydrological MODEL) is now 20 years old and has been the subject of numerous applications to a wide variety of catchments. This paper represents a critical review of some of the issues involved in application of the TOPMODEL concepts, including the basic assumptions involved; the derivation of topographic index distributions from digital terrain data; additional model components; meaning and calibration of the model parameters; and issues involved in model validation and predictive uncertainty. The aim is to provoke a thoughtful approach to hydrological modelling and the interaction of modelling and field work. Some recommendations are made for future modelling practice. © 1997 John Wiley & Sons, Ltd.     

Review of field observations of runoff generation in the tropics

Over the last thirty years, French hydrologists have observed flow characteristics at the outlets of small tropical watersheds, characteristics which vary between arid zones, humid zones and equatorial rainforest zones. The main physical phenomena related to the various types of runoff generation have been observed and analysed. The generation of runoff and subsurface runoff following tropical rainstorms depends on soil surface formations, texture and structure as well as vegetation cover, geological substrata, weathered material and climate.

Analysis is also made of experimental rainfall-runoff relationships, the shape of flood hydrographs and the conditions necessary for the occurrence of steady base flow. Scale effects between experimental plot and watershed have been studied to identify any effect this may have upon runoff generation and flow. Tables detail watershed characteristics and relationships between these and rainfall and runoff.     

A MATLAB implementation of TOPMODEL

The MATLAB SIMULINK programming language is applied to the TOPMODEL rainfall–runoff model. SIMULINK requires a good recognition of model dynamics, which has been achieved here in a version based on the first TOPMODEL (Beven and Kirkby, 1979). Introducing the topographic index distribution in a vector form allows the generalization and simplification of the SIMULINK structure. The SIMULINK version of TOPMODEL has a very easy to understand graphical representation, which shows, in a straightforward way, all the physical interactions that take place in the model. Moreover, owing to its modular structure it is easy to add new and/or develop old submodels, depending on the available data and the goal of the modelling. In the example given here TOPMODEL was extended by two submodels representing the soil moisture and evaporation distribution in the catchment. Preparation of the data and presentation of the results is done in MATLAB. Discharge predictions and spatial patterns of hydrological response are demonstrated for a separate validation period. © 1997 John Wiley & Sons, Ltd.     

Modelling non-stationary water ages in a tropical rainforest: A preliminary spatially distributed assessment

Pristine tropical forests play a critical role in regional and global climate systems. For a better understanding of the eco-hydrology of tropical “evergreen” vegetation, it is essential to know the partitioning of water into transpiration and evaporation, runoff and associated water ages. For this purpose, we evaluated how topography and vegetation influence water flux and age dynamics at high temporal (hourly) and spatial (10 m) resolution using the Spatially Distributed Tracer-Aided Rainfall-Runoff model for the tropics (STARRtropics). The model was applied in a tropical rainforest catchment (3.2 km²) where data were collected biweekly to monthly and during intensive monitoring campaigns from January 2013 to July 2018. The STARRtropics model was further developed, incorporating an isotope mass balance for evapotranspiration partitioning into transpiration and evaporation. Results exhibited a rapid streamflow response to rainfall inputs (water and isotopes) with limited mixing and a largely time-invariant baseflow isotope composition. Simulated soil water storage showed a transient response to rainfall inputs with a seasonal component directly resembling the streamflow dynamics which was independently evaluated using soil water content measurements. High transpiration fluxes (max 7 mm/day) were linked to lower slope gradients, deeper soils and greater leaf area index. Overall water partitioning resulted in 65% of the actual evapotranspiration being driven by vegetation with high transpiration rates over the drier months compared to the wet season. Time scales of water age were highly variable, ranging from hours to a few years. Stream water ages were conceptualized as a mixture of younger soil water and slightly older, deeper soil water and shallow groundwater with a maximum age of roughly 2 years during drought conditions (722 days). The simulated soil water ages ranged from hours to 162 days and for shallow groundwater up to 1,200 days. Despite the model assumptions, experimental challenges and data limitation, this preliminary spatially distributed model study enhances knowledge about the water ages and overall young water dominance in a tropical rainforest with little influence of deeper and older groundwater.     

Effects of the catchment runoff coefficient on the performance of TOPMODEL in rainfall–runoff modelling

Effects of the catchment runoff coefficient on the performance of TOPMODEL in simulating catchment rainfall–runoff relationships are investigated in this paper, with an aim to improve TOPMODEL's simulation efficiency in catchments with a low runoff coefficient. Application of TOPMODEL in the semi‐arid Yihe catchment, with an area of 2623 km² in the Yellow River basin of China, produced a Nash–Sutcliffe model efficiency of about 80%. To investigate how the catchment runoff coefficient affects the performance of TOPMODEL, the whole observed discharge series of the Yihe catchment is multiplied with a larger‐than‐unity scale factor to obtain an amplified discharge series. Then TOPMODEL is used to simulate the amplified discharge series given the original rainfall and evaporation data. For a set of different scale factors, TOPMODEL efficiency is plotted against the corresponding catchment runoff coefficient and it is found that the efficiency of TOPMODEL increases with the increasing catchment runoff coefficient before reaching a peak (e.g. about 90%); after the peak, however, the efficiency of TOPMODEL decreases with the increasing catchment runoff coefficient. Based on this finding, an approach called the discharge amplification method is proposed to enhance the simulation efficiency of TOPMODEL in rainfall–runoff modelling in catchments with a low runoff coefficient. Copyright © 2004 John Wiley & Sons, Ltd.     

Power function decay of hydraulic conductivity for a TOPMODEL‐based infiltration routine

TOPMODEL rainfall‐runoff hydrologic concepts are based on soil saturation processes, where soil controls on hydrograph recession have been represented by linear, exponential, and power function decay with soil depth. Although these decay formulations have been incorporated into baseflow decay and topographic index computations, only the linear and exponential forms have been incorporated into infiltration subroutines. This study develops a power function formulation of the Green and Ampt infiltration equation for the case where the power n = 1 and 2. This new function was created to represent field measurements in the New York City, USA, Ward Pound Ridge drinking water supply area, and provide support for similar sites reported by other researchers. Derivation of the power‐function‐based Green and Ampt model begins with the Green and Ampt formulation used by Beven in deriving an exponential decay model. Differences between the linear, exponential, and power function infiltration scenarios are sensitive to the relative difference between rainfall rates and hydraulic conductivity. Using a low‐frequency 30 min design storm with 4·8 cm h^?1 rain, the n = 2 power function formulation allows for a faster decay of infiltration and more rapid generation of runoff. Infiltration excess runoff is rare in most forested watersheds, and advantages of the power function infiltration routine may primarily include replication of field‐observed processes in urbanized areas and numerical consistency with power function decay of baseflow and topographic index distributions. Equation development is presented within a TOPMODEL‐based Ward Pound Ridge rainfall‐runoff simulation. Copyright © 2006 John Wiley & Sons, Ltd.     

Two‐dimensional physically based finite element runoff model for small agricultural watersheds: II. Model testing and field application

Mathematical models are useful analysis tools to understand problems in watersheds associated with runoff, and to find solutions through land use changes and best management practices. However, before a model is applied in the field, it must be tested and checked to ensure that the model represents the real world adequately. In this paper, a two‐dimensional physically based finite element runoff model ROMO2D has been verified and validated by comparing the model output with analytic solution under simplified conditions, published data, and field measurements. Calibration of the model was done manually through a multi‐objective calibration procedure, using observed field data. Before going for field validation/application of ROMO2D, analysis was carried out to determine the optimum number of finite elements into which the watershed should be discretized and the size of the time step. A sensitivity analysis of the model was performed using the observed values of watershed parameters. The model was applied to a 1·45 ha agricultural watershed located in the Shiwalik foothills (India) to simulate runoff. The results demonstrated the potential of the model to simulate runoff from small agricultural watersheds for individual storm events with reasonable accuracy. Copyright © 2008 John Wiley & Sons, Ltd.     

TOPMODEL在珠江流域布柳河流域的应用及其与新安江模型的比较

TOPMODEL是一种以地形为基础的半分布式流域水文模型.对珠江流域布柳河流域的DEM信息进行处理,提取流域的水系、子流域边界、地形指数及水流路径距离的分布,将TOPMODEL应用于该流域中.另外将新安江模型也应用于该流域进行比较.此外,分析了两种模型结构差异所带来的模拟功能差异.两种模型模拟结果精度差异不大,而TOPMODEL实现了空间产流面积分布的可视化.     

Analytical solution to a bias in the TOPMODEL framework balance

The increasing need for distributed hydrological modelling leads to an intense use of spatially distributed predictions of physically based models, such as TOPMODEL as addressed here. The ability of these models to reproduce the internal behaviour of catchments physically is increasingly tested through field experiments (geochemical investigation, distributed measurements network, etc.). This paper will show that, in the case of TOPMODEL, an implicit approximation remains in the classic derivation of the equations that consists in neglecting the surface of saturated areas with respect to the total surface of the catchment. This simplifying, though unnecessary, approximation leads to a systematic underestimation of the catchment water storage deficit and to divergence in the water budget accounting. This may also significantly change the predicted ratio between subsurface and surface water fluxes in the total discharge. An analytical solution is suggested that leads to water balance accounting which is better defined, and more consistent in comparison with field water storage recording. It is expected that this work will ensure more accurate TOPMODEL predictions, consistent with the assumptions of the model. This will then improve the interpretation of comparisons between results of simulation and field experiments. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling the hydrological response of mediterranean catchments,Prades, Catalonia. The use of distributed models as aids to hypothesis formulation

This paper attempts to extend the physical arguments underlying the distributed TOPMODEL concepts in an application to the strongly seasonal contributing area responses in two adjacent small mediterranean catchments in the Prades region of Catalonia, Spain. A perceptual model of hydrological response in these catchments is used to suggest possible modifications of the model in a hypothesis testing framework, including an attempt to modify the topographic index approach to reflect the expansion of the effective area of subsurface flow during the wetting-up sequence. It is found that slight improvements in modelling efficiency are possible but that different model parameter distributions are appropriate for different parts of the record. The model was much more successful for the catchment producing the higher runoff volumes. © 1997 John Wiley & Sons, Ltd.     

Dune formation on the humid tropical sector of the North Queensland Coast,Australia

Several previous attempts have been made to explain the apparent poor development of coastal dunes in the humid tropics in terms of lack of wind energy, failure of sand supply to the shoreline, excessive climatic wetness, salt crust formation on beaches, and the character of tropical back-beach vegetation. However, recent published reports indicate that coastal dune occurrences are more common in the humid tropics than was formerly thought, throwing suspicion on the idea that environmental conditions militate against dune formation in these areas as a whole. Evidence from the humid tropical sector of the North Queensland coast suggests that the poor development of dunes in this area primarily reflects poor sediment sorting in the beach and nearshore zone and low wind energy at the shoreline due to the nature of the coastal orientation and physiography in relation to the prevailing southeasterly winds. These limiting factors are not unique to humid tropical climates.     

Sensitivity analysis of extended TOPMODEL for agricultural watersheds equipped with tile drains

An extension of TOPMODEL was developed for rainfall–runoff simulation in agricultural watersheds equipped with tile drains. Tile drain functions are incorporated into the framework of TOPMODEL. Nine possible flow generation scenarios are suggested for tile-drained watersheds and applied in the modelling procedure. In the model development, two methods of simulation of the flow in the unsaturated zone were compared: the traditional, physically based storage approach and a new approach using a transfer function. A regionalized sensitivity analysis was used to determine the sensitivity of parameters and to compare the behaviour of the transfer function with that of the simple storage-related formulation. The number of accepted combinations of parameter values, on average, was higher for the transfer function approach than when using a Monte Carlo method of parameter estimation. Since the rainfall–runoff response pattern tends to vary seasonally, seven events distributed throughout a year were used in the sensitivity analysis to investigate the seasonal variation of the hydrological characteristics. © 1997 John Wiley & Sons, Ltd.     

Determination of runoff and sediment yield from a small watershed in sub‐humid subtropics using the HSPF model

The Hydrologic Simulation Programme‐Fortran (HSPF), a hydrologic and water quality computer model, was employed for simulating runoff and sediment yield during the monsoon months (June–October) from a small watershed situated in a sub‐humid subtropical region of India. The model was calibrated using measured runoff and sediment yield data for the monsoon months of 1996 and was validated for the monsoon months of 2000 and 2001. During the calibration period, daily‐calibrated runoff had a Nash‐Sutcliffe efficiency (E_NS) value of 0·68 and during the validation period it ranged from 0·44 to 0·67. For daily sediment yield E_NS was 0·71 for the calibration period and it ranged from 0·68 to 0·90 for the validation period. Sensitivity analysis was performed to assess the impact of important watershed characteristics. The model parameters obtained in this study could serve as reference values for model application in similar climatic regions, with practical implications in watershed planning and management and designing best management practices. Copyright © 2007 John Wiley & Sons, Ltd.     

Using a topographic index to distribute variable source area runoff predicted with the SCS curve‐number equation

Because the traditional Soil Conservation Service curve‐number (SCS‐CN) approach continues to be used ubiquitously in water quality models, new application methods are needed that are consistent with variable source area (VSA) hydrological processes in the landscape. We developed and tested a distributed approach for applying the traditional SCS‐CN equation to watersheds where VSA hydrology is a dominant process. Predicting the location of source areas is important for watershed planning because restricting potentially polluting activities from runoff source areas is fundamental to controlling non‐point‐source pollution. The method presented here used the traditional SCS‐CN approach to predict runoff volume and spatial extent of saturated areas and a topographic index, like that used in TOPMODEL, to distribute runoff source areas through watersheds. The resulting distributed CN–VSA method was applied to two subwatersheds of the Delaware basin in the Catskill Mountains region of New York State and one watershed in south‐eastern Australia to produce runoff‐probability maps. Observed saturated area locations in the watersheds agreed with the distributed CN–VSA method. Results showed good agreement with those obtained from the previously validated soil moisture routing (SMR) model. When compared with the traditional SCS‐CN method, the distributed CN–VSA method predicted a similar total volume of runoff, but vastly different locations of runoff generation. Thus, the distributed CN–VSA approach provides a physically based method that is simple enough to be incorporated into water quality models, and other tools that currently use the traditional SCS–CN method, while still adhering to the principles of VSA hydrology. Copyright © 2004 John Wiley & Sons, Ltd.     

Simulation and comparative study of two types of Topographic Index model for a homogeneous mountain catchment

In order to expand the application range of the classic Topographic Index model (TOPMODEL) and develop a more appropriate submodel of hydrological processes for use in the land surface model, two types of TOPMODEL are investigated, one with saturated hydraulic conductivity change with depth obeying exponential law (classical e-TOPMODEL or e-TOPMODEL for short) and the other obeying general power law (general p-TOPMODEL or p-TOPMODEL for short). Using observation date in the Suomo River catchment located in the upper reaches of the Yangtze River, the sensitivity study of the p-TOPMODEL was conducted and the simulated results from the model were examined and evaluated first, and then the results were compared with the results from the e-TOPMODEL to find the similarities and differences between the two types of models. The main conclusions obtained from the above studies are (1) topographic index and its distribution derived from the p-TOPPMODEL for the Suomo Basin are sensitive to changes of parameter n and m; (2) changes of n and m have impacts on the simulation results of various hydrological components (such as daily runoff, monthly averaged runoff, monthly averaged surface runoff and subsurface runoff), but have the weaker impacts on forty-year averaged total runoff; and (3) for the same value of m, the simulated results of e-TOPMODEL display higher surface runoff and lower subsurface runoff than the general p-TOPMODEL does but multi-year averaged total runoffs produced from the two types of TOPMODEL show insignificant difference. The differences between the two types of models indicate that it is necessary to pay close attention to correct selection from different hydrological models for use in land surface model development. The result mentioned above is useful to provide some referential information for the model selection.