Exploring the effects of hillslope-channel link dynamics and excess rainfall properties on the scaling structure of peak-discharge

Several studies revealed that peak discharges (Q) observed in a nested drainage network following a runoff-generating rainfall event exhibit power law scaling with respect to drainage area (A) as Q(A) = αA^θ. However, multiple aspects of how rainfall-runoff process controls the value of the intercept (α) and the scaling exponent (θ) are not fully understood. We use the rainfall-runoff model CUENCAS and apply it to three different river basins in Iowa to investigate how the interplay among rainfall intensity, duration, hillslope overland flow velocity, channel flow velocity, and the drainage network structure affects these parameters. We show that, for a given catchment: (1) rainfall duration and hillslope overland flow velocity play a dominant role in controlling θ, followed by channel flow velocity and rainfall intensity; (2) α is systematically controlled by the interplay among rainfall intensity, duration, hillslope overland flow velocity, and channel flow velocity, which highlights that it is the combined effect of these factors that controls the exact values of α and θ; and (3) a scale break occurs when runoff generated on hillslopes runs off into the drainage network very rapidly and the scale at which the break happens is determined by the interplay among rainfall duration, hillslope overland flow velocity, and channel flow velocity.     

Effects of hillslope geometry on spatial infiltration using the TOPMODEL and SCS-CN models

ABSTRACT

The estimation of infiltration is a main issue in runoff simulation. The geometry of hillslopes (plan shape and profile curvature) may affect the responses, as well as infiltration over the hillslopes. In this study, the equations of TOPMODEL (a topography-based model) were applied to complex hillslopes to develop the complex TOPMODEL. This model was coupled with the SCS-CN (Soil Conservation Service Curve Number) model to examine the effects of geometry on infiltration and derive a saturation excess-based curve number (CN). The effects of plan shape and profile curvature upon the spatial distribution of CN and infiltration were studied. The results show that convergent hillslopes have 15.4% less infiltration and divergent hillslopes have 7.8% more infiltration than parallel ones. The infiltration over concave hillslopes is 13.5% lower and infiltration over convex hillslopes 5.8% higher than for straight ones. The degree of convergence/divergence has a greater effect on the CN compared to that of profile curvature.     

An integral-balance nonlinear model to simulate changes in soil moisture,groundwater and surface runoff dynamics at the hillslope scale

We present a system of ordinary differential equations (ODEs) capable of reproducing simultaneously the aggregated behavior of changes in water storage in the hillslope surface, the unsaturated and the saturated soil layers and the channel that drains the hillslope. The system of equations can be viewed as a two-state integral-balance model for soil moisture and groundwater dynamics. Development of the model was motivated by the need for landscape representation through hillslopes and channels organized following stream drainage network topology. Such a representation, with the basic discretization unit of a hillslope, allows ODEs-based simulation of the water transport in a basin. This, in turn, admits the use of highly efficient numerical solvers that enable space–time scaling studies. The goal of this paper is to investigate whether a nonlinear ODE system can effectively replicate observations of water storage in the unsaturated and saturated layers of the soil. Our first finding is that a previously proposed ODE hillslope model, based on readily available data, is capable of reproducing streamflow fluctuations but fails to reproduce the interactions between the surface and subsurface components at the hillslope scale. However, the more complex ODE model that we present in this paper achieves this goal. In our model, fluxes in the soil are described using a Taylor expansion of the underlying storage flux relationship. We tested the model using data collected in the Shale Hills watershed, a 7.9-ha forested site in central Pennsylvania, during an artificial drainage experiment in August 1974 where soil moisture in the unsaturated zone, groundwater dynamics and surface runoff were monitored. The ODE model can be used as an alternative to spatially explicit hillslope models, based on systems of partial differential equations, which require more computational power to resolve fluxes at the hillslope scale. Therefore, it is appropriate to be coupled to runoff routing models to investigate the effect of runoff and its uncertainty propagation across scales. However, this improved performance comes at the expense of introducing two additional parameters that have no obvious physical interpretation. We discuss the implications of this for hydrologic studies across scales.     

西江水文干旱历时与强度的遭遇概率分析

利用西江下游马口水文站1959 2009年月径流量数据计算径流干旱指数,经游程理论提取了水文干旱特征值.应用Copula函数分析水文干旱强度和历时之间的联合概率分布.对构建的干旱历时和强度联合分布模式进行分析,结果表明:(1)径流干旱历时和强度之间具有高关联性,秩相关系数达0.617;(2)三参数Weibull分布较好地描述了干旱历时和强度的边缘分布特征;(3)经拟合优度检验结果优选的干旱历时和强度之间的较优连接函数为Archimedean类的Gumbel-Hougaard Copula函数;(4)5～10年重现期和20年重现期的水文干旱分别达到了重旱级别和特旱级别;(5)干旱历时和强度之间的遭遇概率可为特定干旱历时与水文干旱级别或特定干旱强度与干旱历时之间的对应关系提供概率意义上的干旱特征诊断与预测.     

Impulse response of the kinematic diffusion model as a probability distribution of hydrologic samples with zero values

It is hypothesized that the unit impulse response of a linearized kinematic diffusion (KD) model is a probability distribution suitable for frequency analysis of hydrologic samples with zero values. Such samples may include data on monthly precipitation in dry seasons, annual low flow, and annual maximum peak discharge observed in arid and semiarid regions. The hypothesized probability distribution has two parameters, which are estimated using the methods of moments (MOM) and maximum likelihood (MLM). Also estimated are errors in quantiles for MOM and MLM. The distribution shows an equivalency of MOM and MLM with respect to the mean value—an important property for ML-estimation in the case of the unknown true distribution function. The hypothesized KD distribution is tested on 44 discharge data series and compared with the Muskingum-like (M-like) probability distribution function. A comparison of empirical distribution with KD and M-like distributions shows that MOM better reproduces the upper tail of the distribution, while MLM is more robust for higher sample values and more conditioned on the value of the probability of the zero value event. The KD-model is suitable for frequency analysis of short samples with zero values and it is more universal than the M-like model as its modal value cannot be only equaled to zero value but also to any positive value.     

Effects of climatic variability on the hydrologic response of a freshwater watershed

A generalized watershed model was used to evaluate the effects of global climate changes on the hydrologic responses of freshwater ecosystems. The Enhanced Trickle Down (ETD) model was applied to W-3 watershed located near Danville, Vermont. Eight years of field data was used to perform model calibration and verification and the results were presented in Nikolaidis et al., (1993). Results from the Goddard Institute for Space Studies (GISS) and the Geophysical Fluid Dynamics Laboratory (GFDL) general circulation models which simulated the doubling of present day atmospheric CO₂ scenarios were used to perform the hydrologic simulations for the W-3 watershed. The results indicate that the W-3 watershed will experience increases in annual evapotranspiration and decreases in annual outflow and soil moisture. Stochastic models that simulate collective statistical properties of meteorological time series were developed to generate data to drive the ETD model in a Monte-Carlo fashion for quantification of the uncertainty in the model predictions due to input time series. This coupled deterministic and stochastic model was used to generate probable scenarios of future hydrology of the W-3 watershed. The predicted evapotranspiration and soil moisture under doubling present day atmospheric CO₂ scenarios exceed the present day uncertainty due to input time series by a factor greater than 2. The results indicate that the hydrologic response of the W-3 watershed will be significantly different than its present day response. The Enhanced Trickle Down model can be used to evaluate land surface feedbacks and assessing water quantity management in the event of climate change.     

Selection of ground motion time series and limits on scaling

A procedure to select time series for use in non-linear analyses that are intended to result in an average response of the non-linear system is proposed that is not based simply on magnitude, distance, and spectral shape. A simple model of a yielding system is used as a proxy for the non-linear behavior of a more complicated yielding system. As an example, Newmark displacements are used as a proxy for more complex slope-stability models. The candidate scaled time series are evaluated to find those that yield a response of the simple non-linear system that is near the expected response for the design event. Those scaled time series with responses near the expected value are selected as the optimum time series for defining average response even if the scale factors are larger than commonly accepted (e.g. scale factors >factor of 2).     

Functional evaluation of a simplified scaling method for assessing the spatial variability of soil hydraulic properties at the hillslope scale

Abstract

Mapping soil hydraulic parameters with traditional scaling methods that use laboratory-determined hydraulic characteristics (the LAB method) is not always feasible as it involves expensive, time-consuming and sophisticated measurements on soil samples collected in several locations of the study area. An alternative scaling method (the AP method) has been recently proposed to indirectly retrieve the soil hydraulic properties following the Arya-Paris physico-empirical pedotransfer function, which makes use of particle-size distribution and bulk density values. In this synthetic study we verify the performance of the AP method from a functional perspective, by evaluating the differences in the simulated soil water budget through a Monte Carlo approach. Notwithstanding that the AP method can provide soil hydraulic property patterns with faster experimental procedures and minor costs, we observe significant bias in the predicted spatially-averaged soil water budget due to a poor parametric calibration of the AP method and an imprecise identification of the spatial correlation structure of the AP-estimated scaling factors.

Citation Nasta, P., Romano, N., and Chirico, G.B., 2013. Functional evaluation of a simplified scaling method for assessing the spatial variability of soil hydraulic properties at the hillslope scale. Hydrological Sciences Journal, 58 (5), 1059–1071.     

Simulated effect of a forest road on near‐surface hydrologic response and slope stability

Forest management practices often result in significant changes to hydrologic and geomorphic responses at or near the earth's surface. A well‐known, but not fully tested, hypothesis in hillslope hydrology[sol ]geomorphology is that a near‐surface permeability contrast, caused by the surface compaction associated with forest roads, can result in diverted subsurface flow paths that produce increased up‐slope pore pressures and slope failure. The forest road focused on in this study is located in a steep forested, zero‐order catchment within the H. J. Andrews Experimental Forest (Oregon). A three‐phase modelling effort was employed to test the aforementioned hypothesis: (i) two‐dimensional (vertical slice), steady‐state, heterogeneous, saturated subsurface flow simulations at the watershed scale for establishing the boundary conditions for the catchment‐scale boundary‐value problem in (ii); (ii) two‐dimensional (vertical slice), transient, heterogeneous, variably saturated subsurface flow simulations at the catchment scale for estimating near‐surface hydrologic response and pore pressure distributions; and (iii) slope stability analyses, using the infinite slope approach, driven by the pore pressure distributions simulated in (ii), for assessing the impact of the forest road. Both observed and hypothetical rainfall events are used to drive the catchment‐scale simulations. The results reported here support the hypothesis that a forest road can have an effect on slope stability. The permeability contrast associated with the forest road in this study led to a simulated altering of slope‐parallel subsurface flow with increased pore pressures up‐slope of the road and, for a large rainfall event, a slope failure prediction. Copyright © 2005 John Wiley & Sons, Ltd.     

自动提取反应谱特征参数的人工神经网络方法

本以反应谱的特征提取为主题,讨论了被估计目标特征的自动提取问题,中给出了反应谱特征参数提取的神经网络模型,该应用模型通过学习积累有关特征的知识,然后利用已学到的知识自动完成新目标相应特征的提取,体现出神经网络设计智能系统的潜力。     

A procedure for drainage network identification from geomorphology and its application to the prediction of the hydrologic response

Identifying channel initiation points is central to geomorphology and hydrology as they relate morphology, climate, and soil properties at the boundary between different surface runoff paths. Since catchment response is strongly influenced by the dynamics of water movement on the hillslope and in the channel, rainfall-runoff modeling is one of the fields in which improving the identification of channel initiation can lead to benefits. Among the various filtering criteria that can be found in the literature for channel recognition from digital elevation models, the one using contributing area and topographic slope shows interesting features in this context. Nevertheless, the area-slope criterion has been poorly applied. This is mainly due to the difficulties in objectively defining appropriate threshold values. This study proposes a new procedure to assess the area-slope threshold value. The resulting channel network is then used as input to a semi-distributed, event-based rainfall-runoff model able to describe severe rainfall events in small, steep basins. This model accounts for network and hillslope contributions to the total dispersion in the routing process, a key factor in determining the main features of the hydrologic response. In a geomorphologically homogeneous region, the set of model parameters shows interesting invariance properties with respect to storm and basin characteristics.     

A procedure on ground motion selection and scaling for nonlinear response of simple structural systems

This study presents a ground-motion selection and scaling methodology that preserves the basic seismological features of the scaled records with reduced scatter in the nonlinear structural response. The methodology modifies each strong-motion recording with known fundamental seismological parameters using the estimations of ground-motion prediction equations for a given target hazard level. It provides robust estimations on target building response through scaled ground motions and calculates the dispersion about this target. This alternative procedure is not only useful for record scaling and selection but, upon its further refinement, can also be advantageous for the probabilistic methods that assess the engineering demand parameters for a given target hazard level. Case studies that compare the performance of the proposed procedure with some other record selection and scaling methods suggest its usefulness for building performance assessment and loss models. Copyright © 2012 John Wiley & Sons, Ltd.     

1999年11月29日岫岩5.4级地震序列震源参数测定及标度关系分析

根据Andrews谱积分的方法, 采用近震源Brune圆盘模型, 测定了1999年11月29日辽宁岫岩M_S5.4地震序列的震源参数。 结果表明, 得到辐射能量与震级的关系, 与古登堡-里克特给出的关系基本一致; 地震矩与震级的关系中斜率和截距都小于陈培善等给出的全球结果; 视应力与震级呈半对数线性关系; 在双对数坐标下视应力随地震矩的增加而增加; 地震矩随拐角频率的7次方衰减, 在地震序列的不同阶段各参数的拟合关系斜率也不完全相同。     

Effects of hillslope topography on runoff response in a small catchment in the Fudoji Experimental Watershed,central Japan

We investigated the role of different hillslope units with different topographic characteristics on runoff generation processes based on field observations at two types of hillslopes (0·1 ha): a valley‐head (a convergent hillslope) and a side slope (a planar hillslope), as well as at three small catchments having two types of slopes with different drainage areas ranging from 1·9 to 49·7 ha in the Tanakami Mountains, central Japan. We found that the contribution of the hillslope unit type to small catchment runoff varied with the magnitude of rainfall. When the total amount of rainfall for a single storm event was < 35 mm, runoff in the small catchment was predominantly generated from the side slope. As the amount of rainfall increased (>35 mm), the valley‐head also began to contribute to the catchment runoff, adding to runoff from the side slope. Although the direct runoff from the valley‐head was greater than that from the side slope, the contribution from the side slope was quantitatively greater than that from the valley‐head due to the proportionally larger area occupied by the side slope in the small catchment. The storm runoff responses of the small catchments reflected the change in the runoff components of each hillslope unit as the amount of rainfall increased and rainfall patterns changed. However, similar runoff responses were found for the small catchments with different areas. The similarity of the runoff responses is attributable to overlay effects of different hillslope units and the similar composition ratios of the valley‐head and side slope in the catchments. This study suggests that the relative roles of the valley‐head and side slope are important in runoff generation and solute transport as the catchment size increases from a hillslope/headwater to a small catchment. Copyright © 2011 John Wiley & Sons, Ltd.     

Simulated effect of a forest road on near‐surface hydrologic response: redux

In the work reported here the comprehensive physics‐based Integrated Hydrology Model (InHM) was employed to conduct both three‐ and two‐dimensional (3D and 2D) hydrologic‐response simulations for the small upland catchment known as C3 (located within the H. J. Andrews Experimental Forest in Oregon). Results from the 3D simulations for the steep unchannelled C3 (i) identify subsurface stormflow as the dominant hydrologic‐response mechanism and (ii) show the effect of the down‐gradient forest road on both the surface and subsurface flow systems. Comparison of the 3D results with the 2D results clearly illustrates the importance of convergent subsurface flow (e.g. greater pore‐water pressures in the hollow of the catchment for the 3D scenario). A simple infinite‐slope model, driven by subsurface pore‐water pressures generated from the 3D and 2D hydrologic‐response simulations, was employed to estimate slope stability along the long‐profile of the C3 hollow axis. As expected, the likelihood of slope failure is underestimated for the lower pore pressures from the 2D hydrologic‐response simulation compared, in a relative sense, to the higher pore pressures from the 3D hydrologic response simulation. The effort reported herein provides a firm quantitative foundation for generalizing the effects that forest roads can have on near‐surface hydrologic response and slope stability at the catchment scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Effects of experimental uncertainty on the calculation of hillslope flow paths

Measurement uncertainty is a key hindrance to the quantification of water fluxes at all scales of investigation. Predictions of soil‐water flux rely on accurate or representative measurements of hydraulic gradients and field‐state hydraulic conductivity. We quantified the potential magnitude of errors associated with the parameters and variables used directly and indirectly within the Darcy – Buckingham soil‐water‐flux equation. These potential errors were applied to a field hydrometric data set collected from a forested hillslope in central Singapore, and their effect on flow pathway predictions was assessed. Potential errors in the hydraulic gradient calculations were small, approximately one order of magnitude less than the absolute magnitude of the hydraulic gradients. However, errors associated with field‐state hydraulic conductivity derivation were very large. Borehole (Guelph permeameter) and core‐based (Talsma ring permeameter) techniques were used to measure field‐saturated hydraulic conductivity. Measurements using these two approaches differed by up to 3\9 orders of magnitude, with the difference becoming increasingly marked within the B horizon. The sensitivity of the shape of the predicted unsaturated hydraulic conductivity curve to ±5% moisture content error on the moisture release curve was also assessed. Applied moisture release curve error resulted in hydraulic conductivity predictions of less than ±0\2 orders of magnitude deviation from the apparent conductivity. The flow pathways derived from the borehole saturated hydraulic conductivity approach suggested a dominant near‐surface flow pathway, whereas pathways calculated from the core‐based measurements indicated vertical percolation to depth. Direct tracer evidence supported the latter flow pathway, although tracer velocities were approximately two orders of magnitude smaller than the Darcy predictions. We conclude that saturated hydraulic conductivity is the critical hillslope hydrological parameter, and there is an urgent need to address the issues regarding its measurement further. Copyright © 2000 John Wiley & Sons, Ltd.     

Assessing the relative contributions of the flood tide and the ebb tide to tidal channel network dynamics

Flood and ebb currents provide different contributions to the initiation and evolution of tidal channel networks, generating diverse network structures and channel cross-sections. In order to separate the effects of these contributions, a physical model of a sloping tidal-flat basin was set up in the laboratory. Depending on the degree of tidal asymmetry imposed offshore, either flood or ebb currents can be enhanced. The experimental results show that the ebb current has a higher capability to initiate and shape tidal networks than the flood current. Headward erosion is mainly induced by the ebb flow. The slightly inclined flat surface tends to reduce the energy of the flood current and to enhance the ebb current, thus prolonging the duration of morphodynamic activity as well as sediment motion. Overall, flood-dominated tides favour the formation of small-scale channel branches in the upper basin zone, while long lasting ebb-dominated tides result in more complex, wider and deeper tidal networks. © 2019 John Wiley & Sons, Ltd.