Effect of rainfall uncertainty on the performance of physically based rainfall–runoff models

This paper analyses the effect of rain data uncertainty on the performance of two hydrological models with different spatial structures: a semidistributed and a fully distributed model. The study is performed on a small catchment of 19.6 km² located in the north‐west of Spain, where the arrival of low pressure fronts from the Atlantic Ocean causes highly variable rainfall events. The rainfall fields in this catchment during a series of storm events are estimated using rainfall point measurements. The uncertainty of the estimated fields is quantified using a conditional simulation technique. Discharge and rain data, including the uncertainty of the estimated rainfall fields, are then used to calibrate and validate both hydrological models following the generalized likelihood uncertainty estimation (GLUE) methodology. In the storm events analysed, the two models show similar performance. In all cases, results show that the calibrated distribution of the input parameters narrows when the rain uncertainty is included in the analysis. Otherwise, when rain uncertainty is not considered, the calibration of the input parameters must account for all uncertainty in the rainfall–runoff transformation process. Also, in both models, the uncertainty of the predicted discharges increase in similar magnitude when the uncertainty of rainfall input increase.     

Trends in rainfall and runoff in the Blue Nile Basin: 1964–2003

Most of the water from the Nile originates in Ethiopia but there is no agreement on how land degradation or climate change affects the future flow in downstream countries. The objective of this paper is to improve the understanding of future conditions by analysing historical trends. During the period 1964–2003, the average monthly basin‐wide precipitation and monthly discharge data were collected and analysed statistically for two stations in the upper 30% of the Blue Nile Basin and monthly and 10‐day discharge data of one station at the Sudan–Ethiopia border. A rainfall–runoff model examined the causes for observed trends. The results show that, while there was no significant trend in the seasonal and annual basin‐wide average rainfall, significant increases in discharge during the long rainy season (June to September) were observed at all three stations. In the upper Blue Nile, the short rainy season flow (March to May) increased, while the dry season flow (October to February) stayed the same. At the Sudan border, the dry season flow decreased significantly with no change in the short rainy season flow. The difference in response was likely due to the construction of weir in the 1990s at the Lake Tana outlet that affected the upper Blue Nile discharge significantly but affected less than 10% of the discharge at the Sudan border. The rainfall–runoff model reproduced the observed trends, assuming that an additional 10% of the hillsides were eroded in the 40‐year time span and generated overland flow instead of interflow and base flow. Models concerning future trends in the Nile cannot assume that the landscape runoff processes will remain static. Copyright © 2010 John Wiley & Sons, Ltd.     

Understanding coastal wetland hydrology with a new regional‐scale,process‐based hydrological model

Coastal wetlands represent an ecotone between ocean and terrestrial ecosystems, providing important services, including flood mitigation, fresh water supply, erosion control, carbon sequestration, and wildlife habitat. The environmental setting of a wetland and the hydrological connectivity between a wetland and adjacent terrestrial and aquatic systems together determine wetland hydrology. Yet little is known about regional‐scale hydrological interactions among uplands, coastal wetlands, and coastal processes, such as tides, sea level rise, and saltwater intrusion, which together control the dynamics of wetland hydrology. This study presents a new regional‐scale, physically based, distributed wetland hydrological model, PIHM‐Wetland, which integrates the surface and subsurface hydrology with coastal processes and accounts for the influence of wetland inundation on energy budgets and evapotranspiration (ET). The model was validated using in situ hydro‐meteorological measurements and Moderate Resolution Imaging Spectroradiometer (MODIS) ET data for a forested and herbaceous wetland in North Carolina, USA, which confirmed that the model accurately represents the major wetland hydrological behaviours. Modelling results indicate that topographic gradient is a primary control of groundwater flow direction in adjacent uplands. However, seasonal climate patterns become the dominant control of groundwater flow at lower coastal plain and land–ocean interface. We found that coastal processes largely influence groundwater table (GWT) dynamics in the coastal zone, 300 to 800 m from the coastline in our study area. Among all the coastal processes, tides are the dominant control on GWT variation. Because of inundation, forested and herbaceous wetlands absorb an additional 6% and 10%, respectively, of shortwave radiation annually, resulting in a significant increase in ET. Inundation alters ET partitioning through canopy evaporation, transpiration, and soil evaporation, the effect of which is stronger in cool seasons than in warm seasons. The PIHM‐Wetland model provides a new tool that improves the understanding of wetland hydrological processes on a regional scale. Insights from this modelling study provide benchmarks for future research on the effects of sea level rise and climate change on coastal wetland functions and services.     

Comparative analysis of a geomorphology‐based instantaneous unit hydrograph in small mountainous watersheds

A conceptual insytnataneous unit hydrograph (IUH) based on geomorphologival association of linear reservoirs (GR) previously developed by the authors has been compared with other IUH models: a distributed GR variation (GR(v)), the Nash IUH, the Chutha and Dooge IUH, and the Troutman and Karlinger IUH for the analysis of direct runoff hydrographs recorded in three experimental watershed of the north of Spain. The comparison was made through a calibration‐validation process in which a leave‐one‐out cross‐validation method was applied. The results indicate the satisfactory performance of all the models, with the advantage for the GR model of the dependence on only one parameter, which can be identified from the watershed and event characteristics. This property makes its use easier than that of other models. Copyright © 2011 John Wiley & Sons, Ltd.     

Design hydrograph estimation in small and ungauged watersheds: continuous simulation method versus event‐based approach

The proper assessment of design hydrographs and their main properties (peak, volume and duration) in small and ungauged basins is a key point of many hydrological applications. In general, two types of methods can be used to evaluate the design hydrograph: one approach is based on the statistics of storm events, while the other relies on continuously simulating rainfall‐runoff time series. In the first class of methods, the design hydrograph is obtained by applying a rainfall‐runoff model to a design hyetograph that synthesises the storm event. In the second approach, the design hydrograph is quantified by analysing long synthetic runoff time series that are obtained by transforming synthetic rainfall sequences through a rainfall‐runoff model. These simulation‐based procedures overcome some of the unrealistic hypotheses which characterize the event‐based approaches. In this paper, a simulation experiment is carried out to examine the differences between the two types of methods in terms of the design hydrograph's peak, volume and duration. The results conclude that the continuous simulation methods are preferable because the event‐based approaches tend to underestimate the hydrograph's volume and duration. Copyright © 2011 John Wiley & Sons, Ltd.     

Spatiotemporal analysis of climate variability (1971–2010) in spring and summer on the Loess Plateau,China

By using the Variable Infiltration Capacity model with Palmer Drought Severity Index (VIC‐PDSI) model and Standardized Precipitation Index (SPI), spatiotemporal trends of climate variation during the main growing seasons for plants of Loess Plateau between 1971 and 2010 were detected and characterized. The VIC‐PDSI model is established by combining the VIC model with PDSI. The simulation results and the grids system of VIC were applied to substitute for the two‐layer bucket‐type model to do the hydrological accounting, which could improve the physical mechanism of PDSI and expand its application range. Our results suggest that the climate of the study area has experienced a drying and warming trend during the past four decades. Apart from some individual years and regions, there was a perpetuation of water deficit over the Plateau both in spring and summer. The drought frequency increased from southeast to northwest in spring, while the drought frequency decreased from southeast to northwest in summer. The climate in the southern part of the Loess Plateau, accounting for 23.3% of the study region, showed a significant drying and warming trend in spring over the past four decades. The climate variability detected by VIC‐PDSI model shows good agreement with that monitored by SPI. Since a large part of the study region frequently suffered from water shortage during the main growing seasons for plants, people living in such drought‐prone areas should take measures to prevent the negative effects on agricultural production, reforestation, and regional food security caused by drought. Copyright © 2013 John Wiley & Sons, Ltd.     

Extreme precipitation time trends in Ontario, 1960–2010

The impact of climate change on the behaviour of intensity–duration–frequency curves is critical to the estimation of design storms, and thus to the safe design of drainage infrastructure. The present study develops a regional time trend methodology that detects the impact of climate change on extreme precipitation from 1960 to 2010. The regional time trend linear regression method is fitted to different durations of annual maximum precipitation intensities derived from multiple sites in Ontario, Canada. The results show the relationship between climate change and increased extreme precipitation in this province. The regional trend analysis demonstrates, under nonstationary conditions arising from climate change, that the intensity of extreme precipitation increased decennially between 1.25% for the 30‐min storm and 1.82% for the 24‐h storm. A comparison of the results with a regional Mann–Kendall test validates the found regional time‐trend results. The results are employed to extrapolate the intensity–duration–frequency curves temporally and spatially for future decades across the province. The results of the regional time trend assessment help with the establishment of new safety margins for infrastructure design in Ontario. Copyright © 2016 John Wiley & Sons, Ltd.     

On the role of spatial resolution on snow estimates using a process‐based snow model across a range of climatology and elevation

Hydrological processes in mountainous settings depend on snow distribution, whose prediction accuracy is a function of model spatial scale. Although model accuracy is expected to improve with finer spatial resolution, an increase in resolution comes with modelling costs related to increased computational time and greater input data and parameter information. This computational and data collection expense is still a limiting factor for many large watersheds. Thus, this work's main objective is to question which physical processes lead to loss in model accuracy with regard to input spatial resolution under different climatic conditions and elevation ranges. To address this objective, a spatially distributed snow model, iSnobal, was run with inputs distributed at 50‐m—our benchmark for comparison—and 100‐m resolutions and with aggregated (averaged from the fine to the large resolution) inputs from the 50‐m model to 100‐, 250‐, 500‐, and 750‐m resolution for wet, average, and dry years over the Upper Boise River Basin (6,963 km²), which spans four elevation bands: rain dominated, rain–snow transition, and snow dominated below treeline and above treeline. Residuals, defined as differences between values quantified with high resolution (>50 m) models minus the benchmark model (50 m), of simulated snow‐covered area (SCA) and snow water equivalent (SWE) were generally slight in the aggregated scenarios. This was due to transferring the effects of topography on meteorological variables from the 50‐m model to the coarser scales through aggregation. Residuals in SCA and SWE in the distributed 100‐m simulation were greater than those of the aggregated 750 m. Topographic features such as slope and aspect were simplified, and their gradient was reduced due to coarsening the topography from the 50‐ to 100‐m resolution. Therefore, solar radiation was overestimated, and snow drifting was modified and caused substantial SCA and SWE underestimation in the distributed 100‐m model relative to the 50‐m model. Large residuals were observed in the wet year and at the highest elevation band when and where snow mass was large. These results support that model accuracy is substantially reduced with model scales coarser than 50 m.     

Measurements of the initiation of post‐wildfire runoff during rainstorms using in situ overland flow detectors

Overland flow detectors (OFDs) were deployed in 2012 on a hillslope burned by the 2010 Fourmile Canyon fire near Boulder, Colorado, USA. These detectors were simple, electrical resistor‐type instruments that output a voltage (0–2·5 V) and were designed to measure and record the time of runoff initiation, a signal proportional to water depth, and the runoff hydrograph during natural convective rainstorms. Initiation of runoff was found to be spatially complex and began at different times in different locations on the hillslope. Runoff started first at upstream detectors 56% of the time, at the mid‐stream detectors 6%, and at the downstream detectors 38% of the time. Initiation of post‐wildfire runoff depended on the time‐to‐ponding, travel time between points, and the time to fill surface depression storage. These times ranged from 0·5–54, 0·4–1·1, and 0·2–14 minutes, respectively, indicating the importance of the ponding process in controlling the initiation of runoff at this site. Time‐to‐ponding was modeled as a function of the rainfall acceleration (i.e. the rate of change of rainfall intensity) and either the cumulative rainfall at the start of runoff or the soil–water deficit. Measurements made by the OFDs provided physical insight into the spatial and temporal initiation of post‐wildfire runoff during unsteady flow in response to time varying natural rainfall. They also provided data that can be telemetered and used to determine critical input parameters for hydrologic rainfall–runoff models. Copyright © 2015 John Wiley & Sons, Ltd.     

Hydrochemical variations during flood pulses in the south‐west China peak cluster karst: impacts of CaCO3–H2O–CO2 interactions

High‐resolution measurements of rainfall, water level, pH, conductivity, temperature and carbonate chemistry parameters of groundwater at two adjacent locations within the peak cluster karst of the Guilin Karst Experimental Site in Guangxi Province, China, were made with different types of multiparameter sonde. The data were stored using data loggers recording with 2 min or 15 min resolution. Waters from a large, perennial spring represent the exit for the aquifer's conduit flow, and a nearby well measures water in the conduit‐adjacent, fractured media. During flood pulses, the pH of the conduit flow water rises as the conductivity falls. In contrast, and at the same time, the pH of groundwater in the fractures drops, as conductivity rises. As Ca²⁺ and HCO₃^? were the dominant (>90%) ions, we developed linear relationships (both r² > 0·91) between conductivity and those ions, respectively, and in turn calculated variations in the calcite saturation index (SI_C) and CO₂ partial pressure (P) of water during flood pulses. Results indicate that the P of fracture water during flood periods is higher than that at lower flows, and its SI_C is lower. Simultaneously, P of conduit water during the flood period is lower than that at lower flows, and its SI_C also is lower. From these results we conclude that at least two key processes are controlling hydrochemical variations during flood periods: (i) dilution by precipitation and (ii) water–rock–gas interactions. To explain hydrochemical variations in the fracture water, the water–rock–gas interactions may be more important. For example, during flood periods, soil gas with high CO₂ concentrations dissolves in water and enters the fracture system, the water, which in turn has become more highly undersaturated, dissolves more limestone, and the conductivity increases. Dilution of rainfall is more important in controlling hydrochemical variations of conduit water, because rainfall with higher pH (in this area apparently owing to interaction with limestone dust in the lower atmosphere) and low conductivity travels through the conduit system rapidly. These results illustrate that to understand the hydrochemical variations in karst systems, considering only water–rock interactions is not sufficient, and the variable effects of CO₂ on the system should be evaluated. Consideration of water–rock–gas interactions is thus a must in understanding variations in karst hydrochemistry. Copyright © 2004 John Wiley & Sons, Ltd.     

Moisture transport and seasonal variations in the stable isotopic composition of rainfall in Central American and Andean Páramo during El Niño conditions (2015–2016)

High‐elevation tropical grassland systems, called Páramo, provide essential ecosystem services such as water storage and supply for surrounding and lowland areas. Páramo systems are threatened by climate and land use changes. Rainfall generation processes and moisture transport pathways influencing precipitation in the Páramo are poorly understood but needed to estimate the impact of these changes, particularly during El Niño conditions, which largely affect hydrometeorological conditions in tropical regions. To fill this knowledge gap, we present a stable isotope analysis of rainfall samples collected on a daily to weekly basis between January 2015 and May 2016 during the strongest El Niño event recorded in history (2014–2016) in two Páramo regions of Central America (Chirripó, Costa Rica) and the northern Andes (Cajas, south Ecuador). Isotopic compositions were used to identify how rainfall generation processes (convective and orographic) change seasonally at each study site. Hybrid Single Particle Lagrangian Integrated Trajectory model (HYSPLIT) air mass back trajectory analysis was used to identify preferential moisture transport pathways to each Páramo site. Our results show the strong influence of north‐east trade winds to transport moisture from the Caribbean Sea to Chirripó and the South American low‐level jet to transport moisture from the Amazon forest to Cajas. These moisture contributions were also related to the formation of convective rainfall associated with the passage of the Intertropical Convergence Zone over Costa Rica and Ecuador during the wetter seasons and to orographic precipitation during the transition and drier seasons. Our findings provide essential baseline information for further research applications of water stable isotopes as tracers of rainfall generation processes and transport in the Páramo and other montane ecosystems in the tropics.