Statistical identification of orographic effects in the regional analysis of extreme rainfall

Regional models of extreme rainfall must address the spatial variability induced by orographic obstacles. However, the proper detection of orographic effects often depends on the availability of a well‐designed rain gauge network. The aim of this study is to investigate a new method for identifying and characterizing the effects of orography on the spatial structure of extreme rainfall at the regional scale, including where rainfall data are lacking or fail to describe rainfall features thoroughly. We analyse the annual maxima of daily rainfall data in the Campania region, an orographically complex region in Southern Italy, and introduce a statistical procedure to identify spatial outliers in a low order statistic (namely the mean). The locations of these outliers are then compared with a pattern of orographic objects that has been a priori identified through the application of an automatic geomorphological procedure. The results show a direct and clear link between a particular set of orographic objects and a local increase in the spatial variability of extreme rainfall. This analysis allowed us to objectively identify areas where orography produces enhanced variability in extreme rainfall. It has direct implications for rain gauge network design criteria and has led to promising developments in the regional analysis of extreme rainfall. Copyright © 2015 John Wiley & Sons, Ltd.     

Critical rainfall statistics for predicting watershed flood responses: rethinking the design storm concept

Recent advances have been made to modernize estimates of probable precipitation scenarios; however, researchers and engineers often continue to assume that rainfall events can be described by a small set of event statistics, typically average intensity and event duration. Given the easy availability of precipitation data and advances in desk‐top computational tools, we suggest that it is time to rethink the ‘design storm’ concept. Design storms should include more holistic characteristics of flood‐inducing rain events, which, in addition to describing specific hydrologic responses, may also be watershed or regionally specific. We present a sensitivity analysis of nine precipitation event statistics from observed precipitation events within a 60‐year record for Tompkins County, NY, USA. We perform a two‐sample Kolmogorov–Smirnov (KS) test to objectively identify precipitation event statistics of importance for two related hydrologic responses: (1) peak outflow from the Six Mile Creek watershed and (2) peak depth within the reservoir behind the Six Mile Creek Dam. We identify the total precipitation depth, peak hourly intensity, average intensity, event duration, interevent duration, and several statistics defining the temporal distribution of precipitation events to be important rainfall statistics to consider for predicting the watershed flood responses. We found that the two hydrologic responses had different sets of statistically significant parameters. We demonstrate through a stochastic precipitation generation analysis the effects of starting from a constrained parameter set (intensity and duration) when predicting hydrologic responses as opposed to utilizing an expanded suite of rainfall statistics. In particular, we note that the reduced precipitation parameter set may underestimate the probability of high stream flows and therefore underestimate flood hazard. Copyright © 2016 John Wiley & Sons, Ltd.     

The use of the exponent K(q) function to delimit homogeneous regions in regional frequency analysis of extreme annual daily rainfall

The regional frequency analysis of extreme annual rainfall data is a useful methodology in hydrology to obtain certain quantile values when no long data series are available. The most crucial step in the analysis is the grouping of sites into homogeneous regions. This work presents a new grouping criterion based on some multifractal properties of rainfall data. For this purpose, a regional frequency analysis of extreme annual rainfall data from the Maule Region (Chile) has been performed. Daily rainfall data series of 53 available stations have been studied, and their empirical moments scaling exponent functions K(q) have been obtained. Two characteristics parameters of the K(q) functions (γ_max and K(0)) have been used to group the stations into three homogeneous regions. Only five sites have not been possible to include into any homogenous regions, being the local frequency analysis of extreme daily rainfall the most appropriate method to be used at these locations. Copyright © 2014 John Wiley & Sons, Ltd.     

Three dimensional numerical simulation of residential building on shrink–swell soils in response to climatic conditions

Shrink–swell soils can cause distresses in buildings, and every year, the economic loss associated with this problem is huge. This paper presents a comprehensive system for simulating the soil–foundation–building system and its response to daily weather conditions. Weather data include rainfall, solar radiation, air temperature, relative humidity, and wind speed, all of which are readily available from a local weather station or the Internet. These data are used to determine simulation flux boundary conditions. Different methods are proposed to simulate different boundary conditions: bare soil, trees, and vegetation. A coupled hydro‐mechanical stress analysis is used to simulate the volume change of shrink–swell soils due to both mechanical stress and water content variations. Coupled hydro‐mechanical stress‐jointed elements are used to simulate the interaction between the soil and the slab, and general shell elements are used to simulate structural behavior. All the models are combined into one finite element program to predict the entire system's behavior. This paper first described the theory for the simulations. A site in Arlington, Texas, is then selected to demonstrate the application of the proposed system. Simulation results are shown, and a comparison between measured and predicted movements for four footings in Arlington, Texas, over a 2‐year period is presented. Finally, a three‐dimensional simulation is made for a virtual residential building on shrink–swell soils to identify the influence of various factors. Copyright © 2015 John Wiley & Sons, Ltd.     

Soil property variation mapping through data mining of soil category maps

Spatial information on soil properties is an important input to hydrological models. In current hydrological modelling practices, soil property information is often derived from soil category maps by the linking method in which a representative soil property value is linked to each soil polygon. Limited by the area‐class nature of soil category maps, the derived soil property variation is discontinuous and less detailed than high resolution digital terrain or remote sensing data. This research proposed dmSoil, a data‐mining‐based approach to derive continuous and spatially detailed soil property information from soil category maps. First, the soil–environment relationships are extracted through data mining of a soil map. The similarity of the soil at each location to different soil types in the soil map is then estimated using the mined relationships. Prediction of soil property values at each location is made by combining the similarities of the soil at that location to different soil types and the representative soil property values of these soil types. The new approach was applied in the Raffelson Watershed and Pleasant Valley in the Driftless Area of Wisconsin, United States to map soil A horizon texture (in both areas) and depth to soil C horizon (in Pleasant Valley). The property maps from the dmSoil approach capture the spatial gradation and details of soil properties better than those from the linking method. The new approach also shows consistent accuracy improvement at validation points. In addition to the improved performances, the inputs for the dmSoil approach are easy to prepare, and the approach itself is simple to deploy. It provides an effective way to derive better soil property information from soil category maps for hydrological modelling. Copyright © 2014 John Wiley & Sons, Ltd.     

Including spatial distribution in a data‐driven rainfall‐runoff model to improve reservoir inflow forecasting in Taiwan

Multi‐step ahead inflow forecasting has a critical role to play in reservoir operation and management in Taiwan during typhoons as statutory legislation requires a minimum of 3‐h warning to be issued before any reservoir releases are made. However, the complex spatial and temporal heterogeneity of typhoon rainfall, coupled with a remote and mountainous physiographic context, makes the development of real‐time rainfall‐runoff models that can accurately predict reservoir inflow several hours ahead of time challenging. Consequently, there is an urgent, operational requirement for models that can enhance reservoir inflow prediction at forecast horizons of more than 3 h. In this paper, we develop a novel semi‐distributed, data‐driven, rainfall‐runoff model for the Shihmen catchment, north Taiwan. A suite of Adaptive Network‐based Fuzzy Inference System solutions is created using various combinations of autoregressive, spatially lumped radar and point‐based rain gauge predictors. Different levels of spatially aggregated radar‐derived rainfall data are used to generate 4, 8 and 12 sub‐catchment input drivers. In general, the semi‐distributed radar rainfall models outperform their less complex counterparts in predictions of reservoir inflow at lead times greater than 3 h. Performance is found to be optimal when spatial aggregation is restricted to four sub‐catchments, with up to 30% improvements in the performance over lumped and point‐based models being evident at 5‐h lead times. The potential benefits of applying semi‐distributed, data‐driven models in reservoir inflow modelling specifically, and hydrological modelling more generally, are thus demonstrated. Copyright © 2012 John Wiley & Sons, Ltd.     

A new magnitude category disaggregation approach for temporal high‐resolution rainfall intensities

Simulation of quick runoff components such as surface runoff and associated soil erosion requires temporal high‐resolution rainfall intensities. However, these data are often not available because such measurements are costly and time consuming. Current rainfall disaggregation methods have shortcomings, especially in generating the distribution of storm events. The objectives of this study were to improve point rainfall disaggregation using a new magnitude category rainfall disaggregation approach. The procedure is introduced using a coupled disaggregation approach (Hyetos and cascade) for multisite rainfall disaggregation. The new procedure was tested with ten long‐term precipitation data sets of central Germany using summer and winter precipitation to determine seasonal variability. Results showed that dividing the rainfall amount into four daily rainfall magnitude categories (1–10, 11–25, 26–50, >50 mm) improves the simulation of high rainfall intensity (convective rainfall). The Hyetos model category approach (Hyetos_Cat) with seasonal variation performs representative to observed hourly rainfall compared with without categories on each month. The mean absolute percentage accuracy of standard deviation for hourly rainfall is 89.7% in winter and 95.6% in summer. The proposed magnitude category method applied with the coupled Hyetos_Cat–cascade approach reproduces successfully the statistical behaviour of local 10‐min rainfall intensities in terms of intermittency as well as variability. The root mean square error performance statistics for disaggregated 10‐min rainfall depth ranges from 0.20 to 2.38 mm for summer and from 0.12 to 2.82 mm for the winter season in all categories. The coupled stochastic approach preserves the statistical self‐similarity and intermittency at each magnitude category with a relatively low computational burden. Copyright © 2014 John Wiley & Sons, Ltd.     

Severe soil erosion during a 3‐day exceptional rainfall event: combining modelling and field data for a fallow cereal field

Exceptional rainfall events cause significant losses of soil, although few studies have addressed the validation of model predictions at field scale during severe erosive episodes. In this study, we evaluate the predictive ability of the enhanced Soil Erosion and Redistribution Tool (SERT‐2014) model for mapping and quantifying soil erosion during the exceptional rainfall event (~235 mm) that affected the Central Spanish Pyrenees in October 2012. The capacity of the simulation model is evaluated in a fallow cereal field (1.9 ha) at a high spatial scale (1 × 1 m). Validation was performed with field‐quantified rates of soil loss in the rills and ephemeral gullies and also with a detailed map of soil redistribution. The SERT‐2014 model was run for the six rainfall sub‐events that made up the exceptional event, simulating the different hydrological responses of soils with maximum runoff depths ranging between 40 and 1017 mm. Predicted average and maximum soil erosion was 11 and 117 Mg ha^?1 event^?1, respectively. Total soil loss and sediment yield to the La Reina gully amounted to 16.3 and 9.0 Mg event^?1. These rates are in agreement with field estimations of soil loss of 20.0 Mg event^?1. Most soil loss (86%) occurred during the first sub‐event. Although soil accumulation was overestimated in the first sub‐event because of the large amount of detached soil, the enhanced SERT‐2014 model successfully predicted the different spatial patterns and values of soil redistribution for each sub‐event. Further research should focus on stream transport capacity. Copyright © 2014 John Wiley & Sons, Ltd.     

A long‐term study of stable isotopes as tracers of processes governing water flow and quality in a lowland river basin: the upper Thames,UK

A long‐term study of O, H and C stable isotopes has been undertaken on river waters across the 7000‐km² upper Thames lowland river basin in the southern UK. During the period, flow conditions ranged from drought to flood. A 10‐year monthly record (2003–2012) of the main River Thames showed a maximum variation of 3‰ (δ¹⁸O) and 20‰ (δ²H), although interannual average values varied little around a mean of –6.5‰ (δ¹⁸O) and –44‰ (δ²H). A δ²H/δ¹⁸O slope of 5.3 suggested a degree of evaporative enrichment, consistent with derivation from local rainfall with a weighted mean of –7.2‰ (δ¹⁸O) and –48‰ (δ²H) for the period. A tendency towards isotopic depletion of the river with increasing flow rate was noted, but at very high flows (>100 m³/s), a reversion to the mean was interpreted as the displacement of bank storage by rising groundwater levels (corroborated by measurements of specific electrical conductivity). A shorter quarterly study (October 2011–April 2013) of isotope variations in 15 tributaries with varying geology revealed different responses to evaporation, with a well‐correlated inverse relationship between Δ¹⁸O and baseflow index for most of the rivers. A comparison with aquifer waters in the basin showed that even at low flow, rivers rarely consist solely of isotopically unmodified groundwater. Long‐term monitoring (2003–2007) of carbon stable isotopes in dissolved inorganic carbon (DIC) in the Thames revealed a complex interplay between respiration, photosynthesis and evasion, but with a mean interannual δ¹³C‐DIC value of –14.8 ± 0.5‰, exchange with atmospheric carbon could be ruled out. Quarterly monitoring of the tributaries (October 2011–April 2013) indicated that in addition to the aforementioned factors, river flow variations and catchment characteristics were likely to affect δ¹³C‐DIC. Comparison with basin groundwaters of different alkalinity and δ¹³C‐DIC values showed that the origin of river baseflow is usually obscured. The findings show that long‐term monitoring of environmental tracers can help to improve the understanding of how lowland river catchments function. Copyright © NERC 2015. Hydrological Processes © 2015 John Wiley & Sons, Ltd.