Which catchment characteristics control the temporal dependence structure of daily river flows?

Hydrological classification systems seek to provide information about the dominant processes in the catchment to enable information to be transferred between catchments. Currently, there is no widely agreed‐upon system for classifying river catchments. This paper develops a novel approach to classifying catchments based on the temporal dependence structure of daily mean river flow time series, applied to 116 near‐natural ‘benchmark’ catchments in the UK. The classification system is validated using 49 independent catchments. Temporal dependence in river flow data is driven by the flow pathways, connectivity and storage within the catchment and can thus be used to assess the influence catchment characteristics have on moderating the precipitation‐to‐flow relationship. Semi‐variograms were computed for the 116 benchmark catchments to provide a robust and efficient way of characterising temporal dependence. Cluster analysis was performed on the semi‐variograms, resulting in four distinct clusters. The influence of a wide range of catchment characteristics on the semi‐variogram shape was investigated, including: elevation, land cover, physiographic characteristics, soil type and geology. Geology, depth to gleyed layer in soils, slope of the catchment and the percentage of arable land were significantly different between the clusters. These characteristics drive the temporal dependence structure by influencing the rate at which water moves through the catchment and/or the storage in the catchment. Quadratic discriminant analysis was used to show that a model with five catchment characteristics is able to predict the temporal dependence structure for un‐gauged catchments. This method could form the basis for future regionalisation strategies, as a way of transferring information on the precipitation‐to‐flow relationship between gauged and un‐gauged catchments. © 2014 The Authors. Hydrological Processes by published by John Wiley & Sons, Ltd.     

热水条件下花岗质糜棱岩的摩擦滑动实验研究

为了探讨大陆地壳断层深部的力学性质,我们选择了采自红河断裂带的糜棱岩作为实验样品,进行热水条件下的高温高压摩擦滑动实验.实验在一个以气体为介质的高温高压三轴实验系统中进行.实验条件是：有效正应力为200 MPa;孔隙水压为30 MPa（在400 ℃到600 ℃之间为超临界水条件）;温度为100 ℃到600 ℃;轴向加载的速率范围从0.04 μm/s到0.2 μm/s再到1 μm/s.实验结果表明：（1）当温度小于300 ℃时,糜棱岩的摩擦强度随着温度的上升而增大;当温度大于300 ℃时,糜棱岩的摩擦强度随着温度的上升而减小.这种趋势和以往花岗岩的摩擦滑动数据基本一致;（2）糜棱岩在200 ℃和400 ℃时表现为速度弱化,其余温度下为速度强化;（3）糜棱岩与已有花岗岩的摩擦滑动数据并不完全一致;（4）花岗质糜棱岩速度弱化向速度强化转变的温度在430 ℃附近,以此我们可以推测：在变形机制为摩擦滑动的深部条件下,地震成核的深度范围可以比以往的估计更深.     

A search for atmospheric effects on gravity at different time and space scales

We investigate in this paper various approaches to correct gravity changes for the effect of atmospheric pressure changes. Two specific locations are considered: Strasbourg (France) as mid-latitude station, where regular pressure fronts occur and Djougou (Benin) as equatorial station with large thermally driven S₁ and S₂ waves of planetary extent. We first review the classical approaches based on a constant or frequency-dependent admittance using only local pressure and gravity data. We consider then a model of atmospheric loading and show the barometric admittance in terms of elastic, Newtonian and total load, as a function of the distance from the station. We consider both a 2D pressure model (surface loading) and a 2.5D model, where the density decreases with height (standard atmosphere). Assuming horizontal advection in the atmospheric dynamics, we convert this spatially dependent admittance into a frequency-dependent admittance. Using global pressure data from European Centre for Medium-Range Weather Forecasts (ECMWF) at about 12 km spatial resolution and 3 h sampling, we compute the model-predicted pressure admittance for Djougou and Strasbourg and we simulate the frequency dependence inferred from gravity and pressure observations below 4 cycle per day. A long gravity and pressure data set (1996–2013) from Strasbourg is used to investigate the low frequency part of the pressure admittance while a common 2.5 year data set (August 2010–February 2013) for Strasbourg and Djougou is then analyzed to investigate the high frequency part of the admittance. In both cases, our results are in close agreement with the predictions inferred from an atmospheric 2.5D loading model with a distance–time relationship due to horizontal advection. The frequency dependence of the barometric admittance is explained by the competing contributions of Newtonian attraction and elastic surface deformation according to the distance from the gravimeter. In the far field (low frequencies), the magnitude of the admittance decreases with frequency because of the combined elasticity effect and Newtonian attraction (when the atmosphere is below the horizon) while, on the contrary, in the near field (high frequencies), elasticity becomes negligible and the pressure admittance mainly decreases with increasing frequency because of the decreasing attraction effect of the atmospheric masses inside the cylindrical pressure cell centered on the sensor location of decreasing radius. In the last part, we show that there is variability in time in the pressure admittance for both stations.     

人口集聚区吸纳人口迁入的影响因素——以东莞市为例

人口集聚区的地域特征差异对人口迁移行为产生较大影响,从地域特征出发研究人口集聚区对于迁入人口的作用机制可为地域功能区划和配置公共服务设施提供参考。本文以人口迁移的推拉理论为基础,构建人口集聚区对迁入人口拉力作用的指标体系。以东莞市为案例,综合运用主成分分析和多元回归分析方法,识别人口集聚区对迁入人口的作用力以及各类拉力的作用强度。结果显示：① 东莞市各乡镇人口集聚程度差异较大,且空间分布不均衡;在1995-2000年全市人口迅速集聚期,区域人口集聚度主要与地区就业岗位增加、人均收入水平提高、人口迁移的路径依赖相关。② 该市以劳动密集型加工业为主的产业特色给迁入人口提供了更多的就业机会;其相对较高的工资水平,对外来人口形成了强大的吸引力。③ 人口以路径依赖的形式进行链式迁移可以减低迁移过程中的风险和成本。④ 对人口集聚作用由强到弱依次为：路径依赖拉力、就业岗位拉力、收入水平拉力;综合的区域拉力与人口集聚度呈三次曲线关系。     

The inversion of density structure by graphic processing unit (GPU) and identification of igneous rocks in Xisha area

Organic reefs, the targets of deep-water petroleum exploration, developed widely in Xisha area. However, there are concealed igneous rocks undersea, to which organic rocks have nearly equal wave impedance. So the igneous rocks have become interference for future exploration by having similar seismic reflection characteristics. Yet, the density and magnetism of organic reefs are very different from igneous rocks. It has obvious advantages to identify organic reefs and igneous rocks by gravity and magnetic data. At first, frequency decomposition was applied to the free-air gravity anomaly in Xisha area to obtain the 2D subdivision of the gravity anomaly and magnetic anomaly in the vertical direction. Thus, the distribution of igneous rocks in the horizontal direction can be acquired according to high-frequency field, low-frequency field, and its physical properties. Then, 3D forward modeling of gravitational field was carried out to establish the density model of this area by reference to physical properties of rocks based on former researches. Furthermore, 3D inversion of gravity anomaly by genetic algorithm method of the graphic processing unit (GPU) parallel processing in Xisha target area was applied, and 3D density structure of this area was obtained. By this way, we can confine the igneous rocks to the certain depth according to the density of the igneous rocks. The frequency decomposition and 3D inversion of gravity anomaly by genetic algorithm method of the GPU parallel processing proved to be a useful method for recognizing igneous rocks to its 3D geological position. So organic reefs and igneous rocks can be identified, which provide a prescient information for further exploration.     

岩石强度时效性演化模型

在一定的外界荷载作用下,岩石材料的强度会随时间不断降低,在工程上常表现为由于强度的弱化而导致岩体的突发性破坏行为以及岩体的长期强度等不同现象和问题。因此,岩石强度的时间效应研究对岩体工程施工安全和工程长期稳定性都具有重要的指导意义。在对大量试验结果分析的基础上,总结了岩石强度演化规律与内在机制,发现岩石强度的损失与岩石本身的固有强度、所受应力状态（屈服接近度）以及应力作用的持续时间有关,在此基础上,提出了以岩石凝聚力和内摩擦角随上述因素变化为核心的岩石强度时效性演化模型。通过对多组试验数据的模拟分析表明,提出的岩石强度时效性演化模型能很好地描述岩石强度的演化规律,可以较好地定量反映不同应力水平下岩石强度弱化过程,且物理意义明确,为岩石强度的时间效应研究开拓了一个新的思路和方法。     

Spectral analysis of a railway track in contact with a multilayered poroviscoelastic soil subjected to a harmonic moving load

The paper proposes an alternative approach to dispersion curves to analyse resonant phenomena in the context of wave propagation induced by a harmonic load moving over a railway track coupled to a multilayered poroviscoelastic medium. It consists of a shading representation of the spectral response in the frequency–wavenumber domain. Unlike the usual dispersion curves, it can take into account both the damping and the coupled system. Specific attention is paid to the first mode contribution. Critical ranges are determined from this shading representation: (i) for the load frequency when the load moves at a fixed constant speed, and (ii) for the load speed when the load frequency is constant. Both sub-Rayleigh and super-Rayleigh regimes are studied.     

The inversion of density structure by graphic processing unit （GPU） and identification of igneous rocks in Xisha area

Organic reefs, the targets of deep-water petroleum exploration, developed widely in Xisha area. However, there are concealed igneous rocks undersea, to which organic rocks have nearly equal wave impedance. So the igneous rocks have become interference for future exploration by having similar seismic reflection characteristics. Yet, the density and magnetism of organic reefs are very different from igneous rocks. It has obvious advantages to identify organic reefs and igneous rocks by gravity and magnetic data. At first, frequency decomposition was applied to the free-air gravity anomaly in Xisha area to obtain the 2D subdivision of the gravity anomaly and magnetic anomaly in the vertical direction. Thus, the distribution of igneous rocks in the horizontal direction can be acquired according to high-frequency field, low-frequency field, and its physical properties. Then, 3D forward modeling of gravitational field was carried out to establish the density model of this area by reference to physical properties of rocks based on former researches. Furthermore, 3D inversion of gravity anomaly by genetic algorithm method of the graphic processing unit(GPU) parallel processing in Xisha target area was applied, and 3D density structure of this area was obtained. By this way, we can confine the igneous rocks to the certain depth according to the density of the igneous rocks. The frequency decomposition and 3D inversion of gravity anomaly by genetic algorithm method of the GPU parallel processing proved to be a useful method for recognizing igneous rocks to its 3D geological position. So organic reefs and igneous rocks can be identified, which provide a prescient information for further exploration.     

基于非线性时频分析的地震和爆破识别

利用非线性赵-阿特拉斯-马克斯(Zhao-Atlas-Marks)时频分析方法对地震和爆破进行分析,结果显示:爆破的时频谱的时频聚集性比地震的时频聚集性较好;地震P波的平均频率大于爆破P波的平均频率;地震S波的平均时频中心相对于爆破分布在时间—频率面左上角。     

基于频域衰减的时域全波形反演

时域全波形反演由于采用了全频段信息,因此在迭代过程中不同波长的信息不能由低到高的逐步重建,极易陷入局部极小值.本文通过分频段的方式,对地震数据做正反傅里叶变换,利用频域指数衰减的方法逐级分离出地震数据中的高频成分,在时域上实现由低频向高频的波形反演,从而降低了反演的非线性,使不同波长的信息得到稳步恢复.同时,在高频成分衰减的过程中,后至波的能量也被削弱,由此也降低了深层反射在初始反演过程中的干扰.整个反演仅增加对数据做正反傅里叶变换过程,相较于混合域反演,无需提取全部波场的相应频率成分.在计算效率方面,利用GPU进行加速,并采用CUDA自带函数库中cufft来提高计算效率.通过对Marmousi模型测试,验证了所述方法的有效性.