Using hysteretic behaviour and hydrograph classification to identify hydrological function across the “hillslope–depression–stream” continuum in a karst catchment

In cockpit karst landscapes, fluxes from upland areas contribute large volumes of water to low-lying depressions and stream flow. Hydrograph hysteresis and similarity between monitoring sites is important for understanding the space–time variability of hydrologic responses across the “hillslope–depression–stream” continuum. In this study, the hysteretic feature of hydrographs was assessed by characterizing the loop-like relationships between responses at upstream sites relative to subsurface discharge at the outlet of a small karst catchment. A classification of hydrograph responses based on the multi-scale smoothing Kernel -derived distance classifies the hydrograph responses on the basis of similarities between hillslope and depression sites, and those at the catchment outlet. Results demonstrate that the temporal and spatial variability of hydrograph hysteresis and similarity between hillslope flow and outlet stream flow can be explained by the local heterogeneity of depression aquifer. Large depression storage deficits emerging in the highly heterogeneous aquifer produce strong hysteresis and multiple relationships of upstream hydrographs relative to the outlet subsurface discharge. In contrast, when depression storage deficits are filled during consecutive rainfall events, depression hydrographs at the high permeability sites are almost synchronous or exhibit a monotonous function with the hydrographs at the outlet. This reduced hydrograph hysteresis enhances preferential flow paths in fractured rocks and conduits that can accelerate the hillslope flow to the outlet. Therefore, classification of hydrograph similarities between any upstream sites and the catchment outlet can help to identify the dominant hydrological functions in the heterogeneous karst catchment.     

A new numerical approach to the hyperstatic reaction method for segmental tunnel linings

This paper proposes a numerical approach to the hyperstatic reaction method (HRM) for the analysis of segmental tunnel linings. The influence of segmental joints has been considered directly using a fixity ratio that is determined on the basis of the rotational stiffness. The parameters necessary for the calculation are presented. A specific implementation has been developed using a FEM framework. This code is able to consider the three‐dimensional (3D) effect of segment joints in successive rings on the tunnel lining behaviour. The present HRM allows one to take an arbitrary distribution of segment joints along the tunnel boundary into consideration. In addition, the rotational stiffness of segment joints has been simulated using nonlinear behaviour, as it is closer to the true behaviour of a joint than linear or bilinear behaviour. The numerical results of three hypotheses on ring interaction, which allow the 3D effect of a segmental tunnel lining to be taken into account, have been compared with data obtained from the shield‐driven tunnel of the Bologna–Florence high‐speed railway line project. The numerical results presented in the paper show that the proposed HRM can be used to effectively estimate the behaviour of a segmental tunnel lining. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical and theoretical research on stress distribution in the loosening zone of the trapdoor problem

The traditional theory of soil arching effect was developed on the assumption that stress distribution in the loosening zone is uniform. However, because of the deflection of principal stress' direction, the stress distribution in the loosening zone is actually ununiform. For the evaluation of principal stress axis deflection and stress redistribution, a discrete element method numerical model of trapdoor problem is established for the simulation of soil arching effect. Based on the numerical results, an arc shape of major principal stress trajectory and uniform horizontal stress distribution at the same depth of the loosening zone are adopted. An analytical model is raised to estimate the average loosening earth pressure acting on the trapdoor and stress distribution in the loosening zone at a limit state. In addition, comparison studies are carried out between the predictions of the proposed solutions and discrete element method numerical results as well as available model test results, thereby validating the accuracy of the proposed theoretical model. Both numerical and theoretical results indicate that the vertical stress distribution in the loosening zone is obviously ununiform. The load acting in the middle of loosening zone is transferred toward two sides so that the vertical stress distribution in loosening zone is concave.     

Borehole failure analysis in a sandstone under anisotropic stresses

Borehole failure under anisotropic stresses in a sandstone is analyze numerically for various borehole sizes using a nonlinear elastic–plastic constitutive model for a Cosserat continuum. Borehole failure is identified as macroscopic failure of the borehole through the development of shear bands and breakouts. The results compare well both qualitatively and quantitatively with experimental results from polyaxial tests on Red Wildmoor sandstone. They show that the hole size effect of the borehole failure strength is independent of the far‐field stress anisotropy and follows a ? power law of the hole size. A similar scale effect equation with a ? power law is proposed for the scale effect of the maximum plastic shear strain at failure. This equation can be useful for better predicting hole‐size‐dependent failure with standard codes based on classical continua. The effect of stress anisotropy on the borehole failure stress is found to be independent of the hole size. The failure stress decreases linearly to 40% as the stress anisotropy increases. However, the maximum plastic shear strain at failure is stress anisotropy independent and therefore the critical plastic shear strain for failure is only hole‐size dependent. Copyright © 2009 John Wiley & Sons, Ltd.     

Influence of heterogeneity on unsaturated hydraulic properties: (1) local heterogeneity and scale effect

Spatial heterogeneity is ubiquitous in nature, which may significantly affect the soil hydraulic property curves. The models of a closed‐form functional relationship of soil hydraulic property curves (e.g. VG model or exponential model) are valid at point or local scale based on a point‐scale hydrological process, but how do scale effects of heterogeneity have an influence on the parameters of these models when the models are used in a larger scale process? This paper uses a two‐dimensional variably saturated flow and solute transport finite element model (VSAFT2) to simulate variations of pressure and moisture content in the soil flume under a constant head boundary condition. By changing different numerical simulation block sizes, a quantitative evaluation of parameter variations in the VG model, resulting from the scale effects, is presented. Results show that the parameters of soil hydraulic properties are independent of scale in homogeneous media. Parameters of α and n in homogeneous media, which are estimated by using the unsaturated hydraulic conductivity curve (UHC) or the soil water retention curve (WRC), are identical. Variations of local heterogeneities strongly affect the soil hydraulic properties, and the scale affects the results of the parameter estimations when numerical experiments are conducted. Furthermore, the discrepancy of each curve becomes considerable when moisture content becomes closer to a dry situation. Parameters estimated by UHC are totally different from the ones estimated by WRC. Copyright © 2012 John Wiley & Sons, Ltd.     

基于遥感的GEF海河项目实施效果评价

GEF海河项目于2011年完工,该项目在华北地区16个县市实施,惠及2 000多万人。本文从ET、地下水超采量削减、污染物排放量控制、地表水水资源量与利用量以及出入境水质水量5个方面对GEF海河项目实施效果进行了评价。结果表明,全流域ET值减少了68.4 mm,COD和NH3-N入河总量分别减少了36.8%和8.5%、入海总量分别减少了46.1%和36.4%,地表水资源量增加了8.94%,总用水量减少了1.99×108m3,生态环境用水量增加了7.77×108m3,河湖环境明显改善。     

Mesoscale numerical simulation study of warm fog dissipation by salt particles seeding

Based on the dynamic framework of WRF and Morrison 2-moment explicit cloud scheme, a salt-seeding scheme was developed and used to simulate the dissipation of a warm fog event during 6–7 November 2009 in the Beijing and Tianjin area. The seeding effect and its physical mechanism were studied. The results indicate that when seeding fog with salt particles sized 80 μm and at a quantity of 6 gm~(-2) at the fog top, the seeding effect near the ground surface layer is negative in the beginning period, and then a positive seeding effect begins to appear at 18 min, with the best effect appearing at 21 min after seeding operation. The positive effect can last about 35 min. The microphysical mechanism of the warm fog dissipation is because of the evaporation due to the water vapor condensation on the salt particles and coalescence with salt particles.The process of fog water coalescence with salt particles contributed mostly to this warm fog dissipation. Furthermore, two series of sensitivity experiments were performed to study the seeding effect under different seeding amounts and salt particles sizes. The results show that seeding fog with salt particles sized of 80 μm can have the best seeding effect, and the seeding effect is negative when the salt particle size is less than 10 μm. For salt particles sized 80 μm, the best seeding effect, with corresponding visibility of 380 m, can be achieved when the seeding amount is 30 g m~(-2).     

A modeling study of effective radiative forcing and climate response due to tropospheric ozone

This study simulates the effective radiative forcing(ERF) of tropospheric ozone from 1850 to 2013 and its effects on global climate using an aerosol–climate coupled model, BCC AGCM2.0.1 CUACE/Aero, in combination with OMI(Ozone Monitoring Instrument) satellite ozone data. According to the OMI observations, the global annual mean tropospheric column ozone(TCO) was 33.9 DU in 2013, and the largest TCO was distributed in the belts between 30°N and 45°N and at approximately 30°S; the annual mean TCO was higher in the Northern Hemisphere than that in the Southern Hemisphere;and in boreal summer and autumn, the global mean TCO was higher than in winter and spring. The simulated ERF due to the change in tropospheric ozone concentration from 1850 to 2013 was 0.46 W m~(-2), thereby causing an increase in the global annual mean surface temperature by 0.36℃, and precipitation by 0.02 mm d~(-1)(the increase of surface temperature had a significance level above 95%). The surface temperature was increased more obviously over the high latitudes in both hemispheres, with the maximum exceeding 1.4?C in Siberia. There were opposite changes in precipitation near the equator,with an increase of 0.5 mm d~(-1)near the Hawaiian Islands and a decrease of about-0.6 mm d~(-1)near the middle of the Indian Ocean.     

太原市地表热环境时空格局分析

针对鲜有资源型城市的地表热环境研究这一情况,该文基于Landsat8遥感数据,反演了太原市在2013—2015年间共3期的地表温度。提出一个考虑空间格局结构的城市热岛效应指数,在3期数据中,太原市区整体热岛效应指数分别为9.1%、9.8%和8.9%;同时,根据城市功能和路网结构将太原市区划分为6个分区,从分区级层面量化分析了地表热环境的时空分异特征规律,结果表明,太原市区地表温度整体呈现空间分散化特征,不同分区也呈现出不同的时空格局演变过程;最后,选取40个地表温度实测样本点验证反演精度,最大误差绝对值为1.1℃,最小误差绝对值为0.2℃,均方根误差为0.7℃。