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.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

An original semi‐discrete approach to assess gas conductivity of concrete structures

For civil engineering structures with a tightness role, structural permeability is a key issue. In this context, this paper presents a new proposition of a numerical modelling of leakage rate through a cracked concrete structure undergoing mode I cracking. The mechanical state of the material, considered in the framework of continuum mechanics based on finite element modelling, is described by means of the stress‐based nonlocal damage model which takes into account the stress state and provides realistic local mechanical fields. A semi‐discrete method based on the strong discontinuity approach to estimate crack opening is then considered in the post‐treatment phase. Using a Poiseuille's like relation, the coupling between the mechanical state of the material and its dry gas conductivity is performed. For validation purposes, an original experimental campaign is conducted on a dry concrete disc loaded in a splitting setup. During the loading, gas conductivity and digital image correlation analysis are performed. The comparison with the 3D experimental mechanical global response highlights the performance of the mechanical model. The comparison between crack openings measured by digital image correlation and estimated by the strong discontinuity method shows a good agreement. Finally, the results of the semi‐discrete approach coupled with the gas conductivity compared with experimental data show a good estimation of the structural conductivity. Consequently, if the mechanical problem is well modelled at the global scale, then the proposed approach provides good estimation of gas conductivity. Copyright © 2016 John Wiley & Sons, Ltd.     

Can 234U–230Th dating be used to date large semi‐arid tufas? Challenges from a study in the Naukluft Mountains,Namibia

The large, extensive tufa deposits of the semi‐arid Naukluft Mountains, Namibia are potentially important palaeoenvironmental indicators in an area with few proxy records. Tufas are reliable indicators of increased moisture availability, and have been shown to be amenable to ²³⁴U–²³⁰Th dating, although two challenges are detrital contamination and open‐system behaviour. Densely cemented tufa facies are good candidates for dating, minimising these problems. We report attempts to date five densely‐cemented units, which are only found rarely within the Naukluft deposits. We applied a detailed methodology using multiple subsample analysis, measurement of insoluble residues, application of ‘isochron’ mixing lines, and attempted open‐systems modelling, alongside observations of micromorphology and cathodoluminescence in order to assess the validity of any obtained dates. Surprisingly, densely cemented tufas were found not always to be suitable for dating. Two units contained detrital contamination, which could not be corrected for using a single leachate correction or ‘isochron’ methods. Two units contained ‘excess ²³⁰Th’. This could result under a closed‐system if initial (²³⁴U/²³⁸U) was sufficiently high. Alternatively this may be the result of open‐system behaviour, and loss of uranium, or incorporation of initial unsupported ²³⁰Th, which render samples unsuitable for ²³⁴U–²³⁰Th dating. Micromorphological appearance and cathodoluminescence behaviour are used to explore these possibilities. This study exemplifies the need for careful sample selection, and highlights the importance of analysing multiple subsamples from any tufa sample. The detailed methodology applied proves to be a powerful tool for identifying the range of problems that can be encountered when selecting suitable candidate samples for successful dating. It also shows that semi‐arid tufa sequences may contain very little material suitable for dating. A reliable age of c 80 ka was obtained for a banded unit within a large fluvial barrage, with less reliable dates suggesting tufa deposition during times since >350 ka through to the late Holocene. Copyright © 2010 John Wiley & Sons, Ltd.     

Hydrological response of a Brazilian semi‐arid catchment to different land use and climate change scenarios: a modelling study

Brazilian semi‐arid regions are characterized by water scarcity, vulnerability to desertification, and climate variability. The investigation of hydrological processes in this region is of major interest not only for water planning strategies but also to address the possible impact of future climate and land‐use changes on water resources. A hydrological distributed catchment‐scale model (DiCaSM) has been applied to simulate hydrological processes in a small representative catchment of the Brazilian northeast semi‐arid region, and also to investigate the impact of climate and land‐use changes, as well as changes associated with biofuel/energy crops production. The catchment is part of the Brazilian network for semi‐arid hydrology, established by the Brazilian Federal Government. Estimating and modelling streamflow (STF) and recharge in semi‐arid areas is a challenging task, mainly because of limitation in in situ measurements, and also due to the local nature of some processes. Direct recharge measurements are very difficult in semi‐arid catchments and contain a high level of uncertainty. The latter is usually addressed by short‐ and long‐time‐scale calibration and validation at catchment scale, as well as by examining the model sensitivity to the physical parameters responsible for the recharge. The DiCaSM model was run from 2000 to 2008, and streamflow was successfully simulated, with a Nash–Sutcliffe (NS) efficiency coefficient of 0·73, and R² of 0·79. On the basis of a range of climate change scenarios for the region, the DiCaSM model forecasted a reduction by 35%, 68%, and 77%, in groundwater recharge (GWR), and by 34%, 65%, and 72%, in streamflow, for the time spans 2010–2039, 2040–2069, and 2070–2099, respectively, could take place for a dry future climate scenario. These reductions would produce severe impact on water availability in the region. Introducing castor beans to the catchment would increase the GWR and streamflow, mainly if the caatinga areas would be converted into castor beans production. Changing an area of 1000 ha from caatinga to castor beans would increase the GWR by 46% and streamflow by 3%. If the same area of pasture is converted into castor beans, there would be an increase in GWR and streamflow by 24% and 5%, respectively. Such results are expected to contribute towards environmental policies for north‐east Brazil (NEB), and to biofuel production perspectives in the region. Copyright © 2010 John Wiley & Sons, Ltd.     

Gully‐head erosion processes on a semi‐arid valley floor in Kenya: a case study into temporal variation and sediment budgeting

A three year monitoring programme of gully‐head retreat was established to assess the significance of sediment production in a drainage network that expanded rapidly by gully‐head erosion on the low‐angled alluvio‐lacustrine Njemps Flats in semi‐arid Baringo District, Kenya. This paper discusses the factors controlling the large observed spatial and temporal variation in gully‐head retreat rates, ranging from 0 to 15 m a^?1. The selected gullies differed in planform and in runoff‐contributing catchment area but soil material and land use were similar. The data were analysed at event and annual timescales. The results show that at annual timescale rainfall amount appears to be a good indicator of gully‐head retreat, while at storm‐event timescale rainfall distribution has to be taken into account. A model is proposed, including only rainfall (P) and the number of dry days (DD) between storms: which explains 56 per cent of the variation in retreat rate of the single‐headed gully of Lam1. A detailed sediment budget has been established for Lam1 and its runoff‐contributing area (RCA). By measuring sediment input from the RCA, the sediment output by channelized flow and linear retreat of the gully head for nine storms, it can be seen that erosion shifts between different components of the budget depending on the duration of the dry period (DD) between storms. Sediment input from the RCA was usually the largest component for the smaller storms. The erosion of the gully head occurred as a direct effect of runoff falling over the edge (GH_waterfall) and of the indirect destabilization of the adjacent walls by the waterfall erosion and by saturation (GH_mass/storage). The latter component (GH_mass/storage) was usually much larger that the former (GH_waterfall). The sediment output from the gully was strongly related to the runoff volume while the linear retreat, because of its complex behaviour, was not. Overall, the results show that the annual retreat is the optimal timescale to predict retreat patterns. More detailed knowledge about relevant processes and interactions is necessary if gully‐head erosion is to be included in event‐based soil erosion models. Copyright © 2001 John Wiley & Sons, Ltd.     

Evaporation of perennial semi‐arid woodland in southeastern Australia is adapted for irregular but common dry periods

Measurements of water vapour flux from semi‐arid perennial woodland (mallee) were made for 3 years using eddy covariance instrumentation. There have been no previous long‐term, detailed measures of water use in this ecosystem. Latent energy flux (LE) on a half hourly basis was the measure of the combined soil and plant evaporation, ‘evapotranspiration’ (E_LE) of the site. Aggregation over 3 years of the site measured rain (1136 mm) and the estimated evaporation (794 mm) suggests that 342 mm or 30% of rain had moved into or past the root zone of the vegetation. Above average rainfall during 2011 and the first quarter of 2012 (633 mm, 15 months) would likely have been the period during which significant groundwater recharge occurred. At times immediately after rainfall, E_LE rates were the same or exceeded estimates of potential E calculated from a suitably parameterized Penman–Monteith (E_PMo) equation. Apparent free water E from plant interception and soil evaporation was about 2.3 mm and lasted for 1.3 days following rainfall in summer, while in autumn, E was 5.1 mm that lasted over 5.4 days. The leaf area index (LAI) needed to adjust a wind function calibrated Penman equation (E_PMe) to match the E_LE values could be back calculated to generate seasonal change in LAI from 0.12 to 0.46 and compared well with normalized difference vegetation index; r = 0.38 and p = 0.0213* and LAI calculated from digital cover photography. The apparently conservative response of perennial vegetation evaporation to available water in these semi‐arid environments reinforces the conclusion that these ecosystems use this mechanism to survive the reasonably common dry periods. Plant response to soil water availability is primarily through gradual changes in leaf area. Copyright © 2015 John Wiley & Sons, Ltd.     

Streamflow estimation using optimal regional dependency function

The determination of spatial dependency of regionalized variable (ReV) is important in engineering studies. Regional dependency function that leads to calculation of weighting coefficients is required in order to make regional or point‐wise estimations. After obtaining this dependency function, it is possible to complete missing records in the time series and locate new measurement station. Also determination of regional dependency function is also useful to understand the regional variation of ReV. Point Cumulative Semi‐Variogram (PCSV) is another methodology to understand the regional dependency of ReV related to the magnitude and the location. However, this methodology is not useful to determine the weighting coefficient, which is required to make regional and point‐wise estimations. However, in Point Semi‐Variogram (PSV) proposed here, weighting coefficient depends on both magnitude and location. Although the regional dependency function has a fluctuating structure in PSV approach, this function gradually increases with distance in PCSV. The study area is selected in Mississippi river basin with 38 streamflow stations used for PCSV application before. It is aimed to compare two different geostatistical models for the same data set. PSV method has an ability to determine the value of variable along with optimum number of neighbour stations and influence radius. PSV and slope PSV approaches are compared with the PCSV. It was shown that slope slope point semi‐variogram (SPSV) approaches had relative error below 5%, and PSV and PCSV methods revealed relative errors below 10%. Copyright © 2009 John Wiley & Sons, Ltd.     

弧形防波堤波浪力试验研究

在圆筒形防波堤和V形防波堤的基础上,结合离岸堤后形成V形和半圆形连续出现的韵律海岸的地貌平面形态特点,提出一种前墙为连续的半圆筒形防波堤形式。通过在波浪水槽内进行规则波物理模型试验,探究这一新型防波堤的波压力分布规律及波高、周期、水深等因素对波压力的影响。将试验结果与海港水文公式和合田良实公式计算的理论值对比分析,给出了以合田良实公式的折减系数来拟合新型弧形防波堤波浪总水平力的计算公式。结果表明:新型弧形防波堤上的波压力随波高、周期的增大而增大,其水平波浪总力比同等尺度直墙少10%左右。     

半潜式起重平台系泊系统设计与优化研究

近年来,国内兴起了大力发展半潜式起重平台的势头,系泊系统的设计与优化是该类型平台开发中的关键问题之一。采用时域方法对半潜式起重平台系泊系统进行耦合动力分析,研究系泊系统主要参数对其动力响应特性的影响,并给出主要参数的设计优化流程,在此基础上,以能够满足安全校核且性能较优的系泊方案为目标,对该半潜式起重平台的系泊系统进行优化设计,并研究单根系泊缆发生破断对整个系泊系统的影响。结果表明:基于主要参数对系泊性能的影响规律,经优化设计得到的系泊方案,预留的安全欲度适中,能够保障半潜式起重平台在恶劣海况下的生存能力,同时具备较好的经济性。