Effects of initial soil water content and saturated hydraulic conductivity variability on small watershed runoff simulation using LISEM

Abstract

Soil water content (θ) and saturated hydraulic conductivity (K_s) vary in space. The objective of this study was to examine the effects of initial soil water content (θ_i) and K_s variability on runoff simulations using the LImburg Soil Erosion Model (LISEM) in a small watershed in the Chinese Loess Plateau, based on model parameters derived from intensive measurements. The results showed that the total discharge (TD) and peak discharge (PD) were underestimated when the variability of θ_i and K_s was partially considered or completely ignored compared with those when the variability was fully considered. Time to peak (TP) was less affected by the spatial variability compared to TD and PD. Except for TP in some cases, significant differences were found in all hydrological variables (TD, PD and TP) between the cases in which spatial variability of θ_i or K_s was fully considered and those in which spatial variability was partially considered or completely ignored. Furthermore, runoff simulations were affected more strongly by K_s variability than by θ_i variability. The degree of spatial variability influences on runoff simulations was related to the rainfall pattern and θ_i. Greater rainfall depth and instantaneous rainfall intensity corresponded to a smaller influence of the spatial variability. Stronger effects of the θ_i variability on runoff simulation were found in wetter soils, while stronger effects of the K_s variability were found in drier soils. For accurate runoff simulation, the θ_i variability can be completely ignored in cases of a 1-h duration storm with a return period greater than 10 years, while K_s variability should be fully considered even in the case of a 1-h duration storm with a return period of 20 years.

Editor D. Koutsoyiannis; Associate editor A. Fiori     

Radially convergent groundwater flow in sloping terrain

Abstract

The groundwater flow equation governing the elevation (h) of the steady-state phreatic surface in a sloping aquifer fed by constant recharge over a bi-circular sector is rhh′ ? r ² Bh′ + Pr ² ? PR ² = 0, where r is the radial coordinate, P is a constant involving recharge and aquifer properties, and B is the slope of the aquifer—bedrock boundary. The derived flow equation describes radially convergent flow through a sloping aquifer that discharges to a water body of fixed head. One important simplification is that in which the width of the bi-circular sector is constant, and the draining land becomes a rectangular aquifer. The bi-circular sector and rectangular-strip groundwater flow problems are solved in terms of implicit equations. The solutions for the steady-state phreatic surfaces depend on the ratio of recharge to hydraulic conductivity, the slope of the aquifer-bedrock, and the downstream constant-head boundary. Computational examples illustrate the application of the solutions.     

Hydraulic modelling of Athabasca Vallis,Mars

Abstract

Palaeohydraulic modelling is presented for Athabasca Vallis, the youngest known catastrophic flood channel on Mars. This modelling incorporates three significant advantages over previous modelling of Martian channels: a step-backwater hydraulic model; more accurate topography; and improved flood height indicators. The maximum modelled palaeodischarge is between 1 × 10⁶ and 2 × 10⁶m³s^?1 depending on the friction coefficient selected. An anomalously high palaeostage indicator suggests a region of ponded backwater in the channel in which streamlined forms were created through deposition, with the additional possibility of post-flood subsidence/lowering of the channel slope due to magma extrusion.     

Statistical analysis of the effects on overland flow of spatial variability in soil hydraulic conductivity / Analyse statistique des effets de la variabilité spatiale de la conductivité hydraulique du sol sur l'écoulement de surface

Abstract

The effects of spatial variation of the saturated hydraulic conductivity (K _s ) of the soil on the variation of overland flow were tested by analysing 2000 synthetic rainfall—runoff events, all generated from real, observed rainfall events but with runoff modelled by a two-dimensional distributed model using different spatially variable K _s fields in a small (12 ha) agricultural catchment. The purpose is to determine the influence of spatial variation in K _s on runoff generation. The statistical measures used to describe the variation in the generated K _s were its coefficient of variation and correlation length. Both of these had two levels of typical values obtained from field measurements in other studies. The storms were analysed at a general event level, first using simple graphical and statistical methods and then using analysis of variance (ANOVA). The observed scale of the spatial variation of K _s does cause statistically significant variation in overland flow. The graphical analysis showed that the first flow peak in a multi-event storm had the largest variation and that differences were greater in the rising part of the hydrograph than in its recession. The greatest variation in overland flow was produced by the combination of the greater coefficient of variation and the longer correlation lengths. The smallest variation in overland flow was produced by the combination of the smaller coefficient of variation and the shorter correlation lengths. ANOVA showed that the coefficient of variation and correlation length alone did not explain all the variation of the total flow. ANOVA was not very useful due to the many restrictive assumptions that were not satisfied by the nature of the data and therefore analysis methods with less restrictive assumptions need to be tested.     

Nontectonic signals in reprocessed continuous GPS position time series of CMONOC

The nearly nine-year continuous GPS data collected since 1 March 1999 from the Crustal Motion Observation Network of China (CMONOC) were consistently analyzed. Most of the nonlinear movements in the cumulative position time series produced by CMONOC data center disappeared; and more accurate vertical terms and tectonic signals were extracted. Displacements caused by atmospheric pressure loading, nontidal ocean loading, soil moisture mass loading, and snow cover mass loading using the National Centers for Environmental Prediction (NCEP) Reanalysis I/II models and Estimation of the Circulation and Climate of the Ocean (ECCO) data can explain most of the vertical annual terms at many stations, while only parts can be explained at Lhasa and southern coastal sites, indicating that there are some deformation mechanisms that are still unknown or not modeled accurately. The remarkable differences in vertical position time series for short-baseline sites reveal that GPS stations can be greatly affected by local factors; and attention should be paid when explaining observed GPS velocity vectors.     

Motion of animated streamlets appears to surpass their graphical alterations in human visual detection of vector field maxima

ABSTRACT

Animations have become a frequently utilized illustration technique on maps but changes in their graphical loading remain understudied in empirical geovisualization and cartographic research. Animated streamlets have gained attention as an illustrative animation technique and have become popular on widely viewed maps. We conducted an experiment to investigate how altering four major animation parameters of animated streamlets affects people’s reading performance of field maxima on vector fields. The study involved 73 participants who performed reaction-time tasks on pointing maxima on vector field stimuli. Reaction times and correctness of answers changed surprisingly little between visually different animations, with only a few occasional statistical significances. The results suggest that motion of animated streamlets is such a strong visual cue that altering graphical parameters makes only little difference when searching for the maxima. This leads to the conclusion that, for this kind of a task, animated streamlets on maps can be designed relatively freely in graphical terms and their style fitted to other contents of the map. In the broader visual and geovisual analytics context, the results can lead to more generally hypothesizing that graphical loading of animations with continuous motion flux could be altered without severely affecting their communicative power.     

Plant drought tolerance trait is the key parameter in improving the modeling of terrestrial transpiration in arid and semi-arid regions

陆面过程蒸腾作用的模拟制约着天气,气候降水预测的精确度.近几十年来,为了更好地描述植被蒸腾的水力约束,陆面过程模式发展了基于植物性状的植物水力胁迫方案.然而,我们对于植物性状在蒸腾模拟中的地位仍然缺乏了解,植物性状对蒸腾的重要性仍需进一步量化.本研究利用Morris方法评估植物性状参数在通用陆面模式植物水力胁迫方案(CoLM-P₅₀HS)中的重要性,针对17种植物性状,筛选出最为重要的:耐旱性状(P₅₀),气孔性状,和光合作用性状.在12个FLUXNET站点中,参数的重要性由归一化敏感度来衡量.P₅₀的重要性随着降水的减少而增加,而气孔性状和光合作用性状的重要性则随着降水的减少而减少.在干旱或半干旱地区,P₅₀比气孔性状和光合作用性状更重要,这意味着当植物经常经历干旱时,水力安全策略比植物生长策略更关键.而耐旱性状的巨大变异性进一步暗示了多种植物水力安全策略的共存.忽视P₅₀的变异性可能会对陆面过程模式蒸腾作用的模拟造成严重误差.因此,为了更好地表示植物水力功能的变异性,需要增加对耐旱性状的观测并耦合到陆面模式中.     

土壤饱和导水率空间预测的不确定性分析

当土壤转换函数应用于土壤水力性质估计时,对于预测值的不确定性往往容易被忽视。为了有针对性地提出减少这种不确定性的方法和措施,提高土壤转换函数的实际应用能力,以两种现有的土壤转换函数(Vereecken和HYPRES模型)为例,将其应用于山东省平度市土壤饱和导水率的空间预测,并利用拉丁超立方抽样(LHS)方法对预测结果的不确定性进行了分析。结果表明,饱和导水率空间预测的不确定性主要来源于土壤基本性质的空间插值误差和土壤转换函数自身的预测误差。当Vereecken模型应用于饱和导水率空间预测时,预测结果的不确定性主要由土壤基本性质空间插值误差所决定,土壤转换函数预测误差的影响较小,而HYPRES模型则是受二者的双重影响。     

土壤水分对土壤参数的敏感性及其参数优化方法研究

利用2008年1月1日至2009年9月31日黑河流域阿柔冻融观测站的气象和土壤水分数据,采用基于方差的多参数敏感性分析方法研究通用陆面模型(Common Land Model,CoLM)模拟的土壤水分对土壤质地(砂土和黏土)的敏感性,进而采用SCE-UA参数优化算法分别优化土壤质地和土壤水力参数,分析不同优化策略对土壤水分模拟结果的影响。研究结果表明,浅层土壤水分对土壤质地较为敏感,敏感性系数达到了0.45以上,并且砂土含量对土壤水分的影响更为显著;利用SCE-UA算法优化土壤质地或土壤水力参数都可以有效地提高土壤水分的模拟精度,优化土壤水力参数易产生"异参同效"现象,而优化土壤质地能够使土壤水力参数的取值范围更加合理。     

水力类泥石流冲出规模模拟实验

在物源条件一定的情况下,通过实验模拟,研究水力类泥石流在不同流速、流量和纵比降组合条件下的沟道启动、揭底机理和运动特征。用Herschel-ulkley模型分析泥石流运动时的受力情况,整理实验数据获得侵蚀高度与纵比降的关系和流量与揭底距离的关系。