长江上游区域蒸发皿蒸发量变化及其对水分循环的影响

利用长江上游最近30年(66个测站)蒸发皿蒸发量和最近50年(90个测站)的7种气象要素,分析了蒸发皿蒸发量的区域变化趋势和影响蒸发皿蒸发量变化的因素;针对7个水文站的年径流量变化,探讨了蒸发皿蒸发量变化后对水分循环的影响.结果表明,长江上游蒸发皿蒸发量的变化可以划分为三个分区,研究区域东西两侧(青藏高原和大巴山一带)为显著减少区,分别命名为RⅠ和RⅡ,中间(云贵高原北部到黄土高原南缘以及由二者包围的四川盆地一带)为显著增大区,命名为RⅢ区.影响区域蒸发皿蒸发量变化的原因各有不同,青藏高原一带(RⅠ区)蒸发皿蒸发量减少的原因可归结于太阳辐射强度和风动力扰动减弱所致.大巴山一带(RⅡ区)减少原因是太阳辐射强度、风动力扰动强度、湿度条件都在显著下降所引起的.云贵高原到四川盆地一带(RⅢ区)蒸发皿蒸发量增加是环境气温强烈升高,导致其上空大气水汽含量显著减少,大气很干燥,引发蒸发过程加强所致.蒸发皿蒸发量发生变化的直接后果就是导致水分循环强弱发生变化,对于RⅠ区,尽管蒸发皿蒸发量减少,由于降水量和径流量增加的作用,这一区域的水分循环有所加强.在RⅡ区,降水量、径流量和蒸发量都在减少,因此RⅡ区水分循环显著减弱.在RⅢ区,降水量和径流量同时减少,而蒸发量增大,水量消耗增大,因此RⅢ区水分循环有减弱趋势.     

On the computation of the electrical potential inside a horizontally‐layered half‐space

A new formulation is proposed for the electrical potential developed inside a horizontally‐layered half‐space for a direct current point‐source at the surface. The recursion formula for the kernel coefficient in the potential integral is simpler than the generally used two‐coefficient recursion. The numerical difficulties that may occur during the computation of the integrals and near the source axis are examined and solutions are proposed. The set of equations permits a stable and accurate computation of the tabular potential everywhere in the medium.     

Modelling the bulk flow of a bedrock‐constrained,multi‐channel reach of the Mekong River,Siphandone, southern Laos

The general nature of bulk flow within bedrock single‐channel reaches has been considered by several studies recently. However, the flow structure of a bedrock‐constrained, large river with a multiple channel network has not been investigated previously. The multiple channel network of the Siphandone wetlands in Laos, a section of the Mekong River, was modelled using a steady one‐dimensional hydraulic model. The river network is characterized by a spatially‐varying channel‐form leading to significant changes in the bulk flow properties between and along the channels. The challenge to model the bulk flow in such a remote region was the lack of ideal boundary conditions. The flow models considered both low flow, high inbank and overbank flows and were calibrated using SPOT satellite sensor imagery and limited field data concerning water levels. The application of the model highlighted flow characteristics of a large multi‐channel network and also further indicated the field data that would be required to properly characterize the flow field empirically. Important results included the observation that adjacent channels within the network had different water surface slopes for the same moments in time; thus calibration data for modelling similar systems needs to account for these significant local differences. Further, the in‐channel hydraulic roughness coefficient strongly varied from one cross‐section to the next (Manning's ‘n’ range: 0·01 to 0·10). These differences were amplified during low flow but persisted in muted form during high discharges. Copyright © 2012 John Wiley & Sons, Ltd.     

Runout analysis of a large rockfall in the Dolomites/Italian Alps using LIDAR derived particle sizes and shapes

LiDAR data were used to quantify and analyse a rockfall event which occurred in 2003 in the Western Dolomites (Italian Alps). In addition to previously existing airborne laserscanning (ALS) data, high resolution terrestrial laserscanning (TLS) data were collected. By using the original point clouds, the volume, axial ratio and runout length of single boulders as well as the surface roughness in the runout zone of the rockfall were derived. The total volume of the rockfall event of approximately 10 000 m³ was estimated by a reconstruction of the pre‐event surface at the detachment zone. The analysis of the laser scanning data of the accumulation zone revealed a power law scaling for boulder volumes larger than 8 m³. The dependence of runout length on boulder volume is complex; it is moderated by particle sphericity. In addition, we used ALS and TLS data to derive the spatial distribution of surface roughness on the talus cone. TLS allow for more accurate roughness mapping than ALS data, but for most applications the point density of ALS data seems to be sufficiently high to derive measures of roughness. Different sampling approaches for plane fitting on the scale of 5 m did not show significant effects besides the computational time. The results of our analyses provide important perspectives for rockfall modelling and process understanding with potential applications in both ‘applied’ (natural hazards) and ‘pure’ geomorphological research. Copyright © 2012 John Wiley & Sons, Ltd.     

Testing bedload transport equations with consideration of time scales

Bedload transport is known to be a highly fluctuating temporal phenomenon, even under constant (mean) flow conditions, as a consequence of stochasticity, bedform migration, grain sorting, hysteresis, or sediment supply limitation. Because bedload transport formulas usually refer to a single mean transport value for a given flow condition, one can expect that prediction accuracy (when compared to measurements) will depend on the amplitude and duration of fluctuations, which in turn depend on the time scale used for observations. This paper aims to identify how the time scale considered can affect bedload prediction. This was done by testing 16 common bedload transport formulas with four data sets corresponding to different measurement period durations: (i) highly fluctuating (quasi‐)instantaneous field measurements; (ii) volumes accumulated at the event scale on two small alpine gravel‐bed rivers, potentially affected by seasonal fluctuations; (iii) volumes accumulated at the interannual scale in a meandering gravel bed river, thought to be weakly subject to fluctuations; (iv) time‐integrated flume measurements with nearly uniform sediments. The tests confirmed that the longer the measurement period, the better the precision of the formula's prediction interval. They also demonstrate several consequential limitations. Most threshold formulas are no longer valid when the flow condition is below two times the threshold condition for the largest elements' motion on the bed surface (considering D₈₄). In such conditions, equations either predict zero transport, or largely overestimate the real transport, especially when D₈₄ is high. There is a need for new sediment data collected with highly reliable techniques such as recording slot bedload samplers to further investigate this topic. Copyright © 2012 John Wiley & Sons, Ltd.     

Improved estimators of correlation and R2 for skewed hydrologic data

ABSTRACT

The coefficient of determination R² and Pearson correlation coefficient ρ = R are standard metrics in hydrology for the evaluation of the goodness of fit between model simulations and observations, and as measures of the degree of dependence of one variable upon another. We show that the standard product moment estimator of ρ, termed r, while well-behaved for bivariate normal data, is upward biased and highly variable for bivariate non-normal data. We introduce three alternative estimators of ρ which are nearly unbiased and exhibit much less variability than r for non-normal data. We also document remarkable upward bias and tremendous increases in variability associated with r using both synthetic data and daily streamflow simulations from 905 calibrated rainfall–runoff models. We show that estimators of ρ = R accounting for skewness are needed for daily streamflow series because they exhibit high variability and skewness compared to, for example, monthly/annual series, where r should perform well.     

Spatial scale effect on the upper soil effective parameters of a distributed hydrological model

Nonlinear dynamics and spatial variability in hydrological systems make the formulation of scaling theories difficult. Therefore, the development of knowledge related to scale effects, scaling techniques, parameterization and linkages of parameters across scales is highly relevant. The main purpose of this work is to analyse the spatial effect of the static storage capacity parameter H_u and the saturated hydraulic conductivity parameter k_s from microscale (sub‐grid level) to mesoscale (grid level) and its implication to the definition of an optimum cell size. These two parameters describe the upper soil water characteristics in the infiltration process conceptualization of the TETIS hydrological model. At microscale, the spatial heterogeneity of H_u and k_s was obtained generating random parameter fields through probability distribution functions and a spatial dependence model with pre‐established correlation lengths. The effective parameters at mesoscale were calculated by solving the inverse problem for each parameter field. Results indicate that the adopted inverse formulation allows transferring the nonlinearity of the system from microscale to the mesoscale via non‐stationary effective parameters. Their values at each cell and time step are in the range of zero to the mean value of the parameter at microscale. The stochastic simulations showed that the variance of the estimated effective parameters decreases when the ratio between mesoscale cell size and correlation length at microscale increases. For a ratio greater than 1, we found cell sizes having the characteristics of a representative elementary area (REA); in such case, the microscale variability pattern did not affect the system response at mesoscale. Copyright © 2011 John Wiley & Sons, Ltd.     

基于多模型的西宁市人口规模预测

人口是反映国情、国力基本情况的重要指标,是区域研究所必须考虑的重要因素之一。合理、准确地预测城市人口规模,是城市与区域规划中首先要考虑的基本问题,也是保证规划科学性与可实施性的关键性前提。以西宁市2000-2011年历年总人口为样本数据,分别构建了一元线性回归模型、马尔萨斯模型、logistic模型及GM（1,1）模型,并进行模型检验。结果表明：（1）模型均通过模型精度检验且精度较高,GM（1,1）模型拟合度最高,均误差达到0.004%,马尔萨斯模型拟合度最低,为-1.440 8%;（2）分析模型预测精度差异产生原因及适用性,表明深入、准确地分析样本数据特征,恰当选择分析方法对于控制人口预测精度尤为重要。由于西宁市2000-2011年人口样本数据在2005及2009年数据存在波动性,破坏了其与一元线性回归模型及马尔萨斯模型的拟合度,导致在4种模型中,Logistic及GM（1,1）模型预测精度较高,而GM（1,1）模预测精度最高,所以采用GM（1,1）模型进行西宁市人口预测,得到西宁市人口预测的最终结果：2012年西宁市总人口将达到225.89×10⁴人,2015年将达到233.39×10⁴人,2020年将达到246.37×10⁴人。从结果看,未来9 a西宁市人口将呈现持续平稳增长的态势,但随着时间推进人口增长速度将逐渐下降。     

中国天山山区235 a气候变化及降水趋势预测

  利用4个降水重建序列与天山山区10个气象站降水变化的响应关系,重建了天山山区近235 a来气候变化的年降水序列,重建序列的方差解释量达到47.7%。对过去235 a天山山区气候变化中降水变化的特征分析发现：降水大致经历了7干7湿的变化阶段,其中偏湿年份为124 a,多于偏干年份;天山山区的降水量以2.1 a、3.0 a、5.8 a、6.0 a的高频变化和24~25 a的低频变化周期最为显著。研究表明,近百年天山山区气候变化中的降水变化分为两个阶段,20世纪初到80年代,降水量逐渐减少,导致冰川退缩、河流径流量减少、湖泊萎缩、沙漠面积增大、植被覆盖率减少;80年代后期至今,降水量迅速增加,导致河川径流量增加、湖泊面积增大和沙尘暴日数减少,生态环境趋于好转,但同时造成洪水和地质灾害频繁。采用基于方差分析的周期叠加方法对天山山区未来23 a气候变化中降水的演变趋势进行预测,结果表明：未来23 a气候变化中天山山区以湿润为主,期间有部分年份可能较常年干旱。     

云南近50 年来的气候变化

利用1961—2008 年云南省124 站及其相邻省市区的37 个地面站的逐日平均气温、降水量资料,客观分析了考虑气温垂直递减率的气温细网格数据。在此基础上,分析研究了云南近50 年来的气候变化。(1)云南近50 年来的气温变化,趋势上与北半球和全国一致,但气候变暖的速率相对缓慢;(2) 云南气候带面积有明显变化,北热带和南亚热带面积增加,而中亚热带、北亚热带和温带的面积减少;(3) 云南冬季是气温上升幅度最大的季节,然后依次是秋季、夏季和春季,虽然近50 年来云南大部分区域气温呈上升趋势,但少数区域气温却呈下降趋势,这些降温区主要集中在低海拔河谷地区;(4) 云南近50 年来降水量年平均变化不大,但季节和空间分布的变化却比较明显,雨季和主汛期降水量呈下降趋势,干季呈上升趋势;在空间分布上,滇东地区呈一致性下降趋势,滇中呈一致性上升趋势;滇西和滇南地区降水量的增减趋势呈交错分布。