青海钻孔应变映震效能研究

利用青海省6套YRY-4分量钻孔应变仪2007—2015年的观测资料,将该地区符合扫描范围的地震从震前异常和同震特征两方面进行研究。结果显示:湟源、德令哈分量钻孔应变差应变曲线地震前大多有短期超差异常;分量钻孔应变仪的同震响应幅度与震级、震中距都有关。同一地震,多台仪器多数情况同震响应幅度与震中距成反比,少数成正比;一般单台仪器同震响应幅度与地震震级成正比关系,震级越大同震响应幅度越大;分量钻孔应变仪的同震响应延迟时间大多与地震震中距成正比,震中距越远同震响应延迟时间越长;分量钻孔应变仪同震响应持续时间与震级基本成正比,震级越大同震持续时间越长。     

Travel time computations using a compact eikonal equation for vertical transverse isotropic media

Eikonal solvers often have stability problems if the velocity model is mildly heterogeneous. We derive a stable and compact form of the eikonal equation for P‐wave propagation in vertical transverse isotropic media. The obtained formulation is more compact than other formulations and therefore computationally attractive. We implemented ray shooting for this new equation through a Hamiltonian formalism. Ray tracing based on this new equation is tested on both simple as well as more realistic mildly heterogeneous velocity models. We show through examples that the new equation gives travel times that coincide with the travel time picks from wave equation modelling for anisotropic wave propagation.     

Using the applied element method for modelling calcium silicate brick masonry subjected to in‐plane cyclic loading

The response of calcium silicate unreinforced masonry construction to horizontal cyclic loading has recently become the focus of experimental and numerical research, given its extensive use in some areas of the world that are now exposed to induced earthquakes (eg, north of the Netherlands). To assess the seismic behaviour of such construction, a relatively wide range of modelling methodologies are available, amongst which the discrete elements approach, which takes into account the intrinsic heterogeneity of a brick‐mortar assembly, can probably be deemed as the most appropriate computational procedure. On the other hand, however, since discrete elements numerical methods are based on a discontinuum domain, often they are not able to model every stage of the structural response adequately, and because of the high computational burden required, the analysis scale should be chosen carefully. The applied element method is a relatively recent addition to the discrete elements family, with a high potential for overcoming the aforementioned limitations or difficulties. Initially conceived to model blast events and concrete structures, its use in the earthquake engineering field is, of late, increasing noticeably. In this paper, the use of the applied element method to model the in‐plane cyclic response of calcium silicate masonry walls is discussed and scrutinised, also through the comparison with experimental results of in‐plane cyclic shear‐compression tests on unreinforced masonry walls.     

Morphodynamic assessment of side channel systems using a simple one‐dimensional bifurcation model and a comparison with aerial images

Side channel construction is a common intervention applied to increase a river's conveyance capacity and to increase its ecological value. Past modelling efforts suggest two mechanisms affecting the morphodynamic change of a side channel: (1) a difference in channel slope between the side channel and the main channel and (2) bend flow just upstream of the bifurcation. The objective of this paper was to assess the conditions under which side channels generally aggrade or degrade and to assess the characteristic timescales of the associated morphological change. We use a one‐dimensional bifurcation model to predict the development of side channel systems and the characteristic timescale for a wide range of conditions. We then compare these results to multitemporal aerial images of four side channel systems. We consider the following mechanisms at the bifurcation to be important for side channel development: sediment diversion due to the bifurcation angle, sediment diversion due to the transverse bed slope, partitioning of suspended load, mixed sediment processes such as sorting at the bifurcation, bank erosion, deposition due to vegetation, and floodplain sedimentation. There are limitations to using a one‐dimensional numerical model as it can only account for these mechanisms in a parametrized manner, but the model reproduces general behaviour of the natural side channels until floodplain‐forming processes become important. The main result is a set of stability diagrams with key model parameters that can be used to assess the development of a side channel system and the associated timescale, which will aid in the future design and maintenance of side channel systems. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Modelling soil moisture in a high‐latitude landscape using LiDAR and soil data

Soil moisture has a pronounced effect on earth surface processes. Global soil moisture is strongly driven by climate, whereas at finer scales, the role of non‐climatic drivers becomes more important. We provide insights into the significance of soil and land surface properties in landscape‐scale soil moisture variation by utilizing high‐resolution light detection and ranging (LiDAR) data and extensive field investigations. The data consist of 1200 study plots located in a high‐latitude landscape of mountain tundra in north‐western Finland. We measured the plots three times during growing season 2016 with a hand‐held time‐domain reflectometry sensor. To model soil moisture and its temporal variation, we used four statistical modelling methods: generalized linear models, generalized additive models, boosted regression trees, and random forests. The model fit of the soil moisture models were R² = 0.60 and root mean square error (RMSE) 8.04 VWC% on average, while the temporal variation models showed a lower fit of R² = 0.25 and RMSE 13.11 CV%. The predictive performances for the former were R² = 0.47 and RMSE 9.34 VWC%, and for the latter R² = 0.01 and RMSE 15.29 CV%. Results were similar across the modelling methods, demonstrating a consistent pattern. Soil moisture and its temporal variation showed strong heterogeneity over short distances; therefore, soil moisture modelling benefits from high‐resolution predictors, such as LiDAR based variables. In the soil moisture models, the strongest predictor was SAGA (System for Automated Geoscientific Analyses) wetness index (SWI), based on a 1 m² digital terrain model derived from LiDAR data, which outperformed soil predictors. Thus, our study supports the use of LiDAR based SWI in explaining fine‐scale soil moisture variation. In the temporal variation models, the strongest predictor was the field‐quantified organic layer depth variable. Our results show that spatial soil moisture predictions can be based on soil and land surface properties, yet the temporal models require further investigation. Copyright © 2017 John Wiley & Sons, Ltd.     

On the importance of soil moisture in calibration of rainfall–runoff models: two case studies

Streamflow modelling results from the GR4H and PDM hydrological models were evaluated in two Australian sub-catchments, using (1) calibration to streamflow and (2) joint-calibration to streamflow and soil moisture. Soil moisture storage in the models was evaluated against soil moisture observations from field measurements. The PDM had the best performance in terms of both streamflow and soil moisture estimations during the calibration period, but was outperformed by GR4H during validation. It was also shown that the soil moisture estimation was improved significantly by joint-calibration for the case where streamflow and soil moisture estimations were poor. In other cases, addition of the soil moisture constraint did not degrade the results. Consequently, it is recommended that GR4H be used, in preference to the PDM, in the foothills of the Murrumbidgee catchment or other Australian catchments with semi-arid to sub-humid climate, and that soil moisture data be used in the calibration process.     

川滇块体东边界主要断裂带现今运动特征分析

基于2009年以来的GPS观测数据,利用块体模型和GPS剖面方法分别计算川滇块体东边界主要断裂带的滑动速度,并结合跨断裂带的区域应变时间序列分析断裂带现今的运动特征。结果表明:从速度场变化来看,2013—2015期的速度场在川滇块体东北部有东向增加的微弱变化;从滑动速率结果来看,鲜水河北段的左旋走滑运动有所增强,拉张运动有所增加;小江断裂带的左旋走滑运动普遍有微弱的增强;从去掉线性的区域应变时间序列结果来看,小江断裂带南段主张应变在2014年底出现了趋势性转折,值得进一步关注。     

2017年四川九寨沟7.0级地震前地震应变场分析

以地震应变场作为地震活动的变量,通过自然正交函数展开方法,计算2017年8月8日四川九寨沟7.0级地震前的地震应变场,提取出震前时间因子的异常变化。计算结果发现应变场前4个时间因子在震前1~3年的中短期异常并不显著,只有第2个和第4个时间因子震前有小幅度的异常变化,分析其原因,时间因子可能受到2008年汶川8.0级地震和2013年岷县漳县6.6级地震前大幅度异常的影响。与时间因子对应的空间等值线形成局部应变高值异常危险区,可能是2013年岷县漳县6.6级地震和2017年九寨沟7.0级地震空间异常的中短期特征。再对比九寨沟7.0级地震前后与松潘-平武7.2级双震的空间异常随时间演变,分析异常发展变化模式的差异,最后应用地震资料的累积频次从物理角度解释异常形成的机制。     

多方向不规则波传播变形数值模拟

在推广的缓坡方程数学模型基础上建立了多方向不规则波数学模型,综合考虑了波浪折射、绕射、反射、底摩擦和风能输入等因素。基于线性波浪理论,将波浪方向谱在频率和方向上按等能量分割法离散后,分别计算各组成波的传播变形,再计算合成波要素。缓坡方程数学模型采用改进的ADI法求解,计算效率高,稳定性好。采用椭圆形浅滩不规则波模型试验结果和单突堤不规则波绕射理论解对数学模型进行了验证,数值模拟结果和试验值及理论解符合良好。利用该模型进行了某港港内波浪折射、绕射和反射的联合数值模拟,给出了合理的港内波高分布。