2008年汶川地震近断裂区域变形场的空间分布

应用Yoshimitsu Okada及Steketee静态断裂位错模型及目前汶川地震研究成果,理论上反演龙门山近断裂区域在汶川地震中从震源深处至地表的变形(应变)场分布.通过矢量合成及坐标变换方法将Yoshimitsu Okada单一静态断裂位错模型应用于龙门山断裂带错裂研究中,重点研究汶川地震中因龙门山中央主断裂中、北段和山前断裂中段共同错动形成的应变场.模拟揭示了应变场平面和垂直分量的空间分布特征.研究发现:应变场垂直分量的幅值在震源深处大于地表,但在地表分布范围较大;应变场沿断裂走向分量的幅值和影响范围由震源深度至地表均逐渐增加;应变场垂直于断裂走向分量的变化情况类似于垂直向应变场;震源深度附近的大变形可能为龙门山断裂在地震中强烈挤压形成的塑性变形.研究表明从震源至地表,所有应变场分量从断裂两侧附近位置至远离断裂,幅值均迅速衰减;其中地表部分的最大变形与现有计算成果近似,变形场的变化特征与地表科考近似,断裂两盘变形方向与震后InSAR资料揭示的情况一致.     

Strain variation with progressive deformation in basement-involved trishear structures

Basement-involved structures associated with reverse, vertical and normal faults commonly involve non-parallel shear within a triangular deformation (trishear) zone located on the front limbs of the structures. Deformation within the trishear zone is characterized by shear gradients and an associated decrease in the dips of the beds in stratigraphically higher units. Geometric models suggest that the layer-parallel strain within the trishear zone depends on the type of fault (normal, reverse, or vertical), the dip and throw of the fault, the dip of the anticlinal or synclinal axial surfaces, and the distance of any unit above the initial tip of the trishear zone, located at the basement-sediment contact. At any given location, reverse faults typically show increasing layer parallel shortening, followed by decreasing layer parallel shortening and a transition to extension, with increasing throw. The transition from contraction to extension occurs at lower values of throw for stratigraphically lower units and also for faults with smaller dips. Vertical and normal faults exhibit increasing layer-parallel extension of all units with increasing throw, with larger extension for stratigraphically lower units. Experimental models suggest that the trishear zone can expand with increasing fault throw. The strain within the trishear zones is accommodated largely by secondary faults, which are rotated with progressive deformation. The strain variations in the experiments closely mimic those predicted by the geometric models for reverse, vertical, and normal faults.     

海原形变观测资料异常分析

对海原台CZB-1型竖直摆钻孔倾斜仪记录到的地表倾斜方向的异常变化进行了分析,同时对该台YRY-4型分量应变仪依据自检方法表现出的压缩及伸张性异常进行了提取,得出以上两测项观测资料在宁夏南部发生的3次Ms3．0级以上有感地震的2次地震前出现了一定的异常信息,与地震赢救蚺蕴。     

Characterization of landslide ground surface kinematics from terrestrial laser scanning and strain field computation

Assessment and mitigation of the risk induced by landslide activation need an appropriate phenomenon investigation, to obtain useful information about the failure processes. The first step is the complete kinematics characterization of the landslide ground surface, by evaluating the involved displacement and deformation patterns. A dense displacement field can be obtained from comparison of a series of multi-temporal observations performed by means of terrestrial laser scanning. Subsequently, the strain field can be computed from displacement vectors. In this paper, a modified least square technique is employed to compute the strain on the nodes of a regular grid (2D approach) or on the points of a digital terrain model (3D approach). Such a computation takes into account the displacements, their spatial distribution, as well as the measurement and modelling errors. A scale factor is introduced in order to emphasize the contributions of the experimental points on the basis of their distance from each computation point, and to recognize possible scale-depending behaviours. This method has been implemented in Matlab and applied on two landslides located in the northeastern Italian Alps (Lamosano and Perarolo di Cadore). The experiments show that different kinematics can be recognized, and the presence and influence of eventual discontinuities can be revealed.     

Refined models of gravitational potential energy compared with stress and strain rate patterns in Iberia

This study examines the role of gravitational potential energy (GPE) in generating second-order (spatial scale ∼10² km) variations in the Iberia stress and strain-rate patterns. We present a new map of present-day strain rate field derived from the secular velocity field computed using all available continuously operating Global Navigation Satellite Systems (GNSS) stations in Iberia. The estimated strain rate field is generally consistent with the tectonic framework of the Iberian region, even though sporadic sharp local variations downgrade its correlation with the regional stress patterns. Many of the sharp spatial variations in the strain rate map are consistent with local changes of deformation style determined by prevailing faults. To obtain a more accurate estimate of GPE we use new data on the structure of the crust and apply a thin sheet approach using a 3-D definition of deviatoric stress. The GPE is derived from two isostatically compensated models (GPEd and GPEe compensated by density and elevation adjustment, respectively) and from the truncated geoid (GPEg). The GPE stresses are then summed with the first-order stress field due to the Eurasia–Nubia (EU–NU) convergence and the results compared with both the stress and strain rate data. In agreement with previous studies, we find that the GPE does not significantly change the NW–SE average direction of the most compressive stress (SHmax) imposed by the EU–NU collision, its main effect being to cause spatially changing stress regimes. From the analysis of the different GPE models we find: (1) in the Pyrenees, the tectonic forces have a secondary role when compared to the GPE. In this region, the model that best correlates with observations is the one emphasizing the role of surface elevation as a source of GPE (GPEe); (2) in the Iberian Chain and the Betics, the GPE imposes NE–SW extension consistent with a strike-slip regime and is equally (GPEe) or more (GPEg) important than the tectonic forces. In these regions, both deep heterogeneities associated with mantle convection and elevation are important sources of GPE; (3) in western Iberia, the GPE differences work against dominant tectonic forces by reducing the SHmax magnitude. The GPEg model is the one that best predicts the average strike-slip regime in Galicia; and finally (4) in the Gulf of Cadiz the gravitational potential stresses have a minor role and the style of deformation is clearly controlled by the tectonic forces.     

基于Chebyshev自褶积组合窗的有限差分算子优化方法

有限差分法广泛应用于地震波数值模拟、成像和波形反演中,差分数值解的精度直接影响着地震成像和反演的效果.因为有限差分算子可以通过截断伪谱法的空间褶积序列得到,而截断窗函数的属性影响有限差分算子逼近微分算子的精度.具体地讲,窗函数的幅值响应的主瓣和旁瓣决定了有限差分算子逼近的精度,主瓣越窄,旁瓣衰减越大,则有限差分算子逼近微分算子的精度越高,更好地压制数值频散.基于此认识,本文提出了一种基于Chebyshev自褶积组合窗截断逼近的有限差分算子优化方法.Chebyshev自褶积组合窗的主瓣较窄,且旁瓣衰减大,其可通过只调节三个参数,更直观和可视化地控制主瓣和旁瓣的形状,改变有限差分算子逼近微分算子的精度;该窗函数截断逼近的有限差分算子不仅有较大的谱覆盖范围,而且精度误差波动较小,这表明低阶的差分算子可以达到高阶算子的精度,且逼近误差更稳定;从经济上来讲,将有效地减少模拟计算花费,提高计算效率.     

2000年中国地区地磁场长期变化的区域特征

本文用中国地磁场复测点和地磁台站资料,建立了2000年代中国参考地磁场长期变化模型CGRF-SV2000.模型显示,中国地区地磁场变化比较平缓,X、Y、Z、H、D、I、F七个地磁要素的“无符号平均年变率”分别为122 nT/a、82 nT/a、438 nT/a、118 nT/a、096(′)/a、299(′)/a、224 nT/a,比国际参考地磁场IGRF给出的全球年变率约小1/3到1/2.各地磁要素的变化显示,地磁北极正在向中国移近,或者说,中国正在向高地磁纬度方向移动,平均移动速度约为3(′)/a.磁偏角变化还显示,中国地区东西部偏角差异继续扩大.作为检验和对比研究,本文利用第8代国际参考地磁场（IGRF）模型,分析了全球地磁场长期变化的时空特征,讨论了全球长期变模型IGRF_SV与中国长期变模型CGRF_SV的异同点.对比分析表明,中国地区地磁场的长期变化与全球长期变化总趋势基本符合,但是,CGRF_SV也表现出一些特有的局部异常特征.     

Modeling the strain localization around an underground gallery with a hydro-mechanical double scale model; effect of anisotropy

This paper concerns the application of a finite element squared approach for modelling hydromechanical coupling in the simulation of gallery excavation in the context of radioactive waste repositories. The micromechanics of Callovo-Oxfordian claystone is modelled at the microscale, taking into account the interaction of different mechanical constituents and its interaction with pore fluid. In a framework of computational homogenization, the micromechanical behaviour is coupled to the macroscale boundary value problem of a poromechanical continuum with local second gradient paradigm. The simulations concern several cases of the “Transverse action” benchmark by Andra, in the context of which the model is used.     

The strain seismograms of P- and S-waves of a local event recorded by four-gauge borehole strainmeter

At a sampling rate of 100 samples per second,the YRY-4 four-gauge borehole strainmeters(FGBS) are capable of recording transient strains caused by seismic waves such as P and S waves or strain seismograms. At such a high sampling rate, data from the YRY-4 strainmeters demonstrate fairly satisfactory self-consistency. The strain tensor seismograms demonstrate the senses of motion of P waves, that is, the type of seismic wave travels in the direction of the maximum normal strain change. The observed strain patterns of S waves significantly differ from those of P waves and should contain information about the source mechanism. Spectrum analysis shows that the strain seismograms are consistent with conventional broadband seismograms from the same site.     

Hydrogeochemical study in the Main Ethiopian Rift: new insights to the source and enrichment mechanism of fluoride

The central Main Ethiopian Rift suffers a severe water quality problem, characterized by an anomalously high fluoride (F) content that causes an endemic fluorosis disease. The current study, conducted in the Ziway–Shala lakes basin, indicates that the F content exceeds the permissible limit for drinking prescribed by the World Health Organization (WHO; 1.5 mg/l) in many important wells (up to 20 mg/l), with even more extreme F concentration in hot springs and alkaline lakes (up to 97 and 384 mg/l respectively). The groundwater and surface water from the highlands, typically characterized by low total dissolved solids (TDS) and Ca (Mg)–HCO₃ hydrochemical facies, do not show high F content. The subsequent interaction of these waters with the various rocks of the rift valley induces a general increase of the TDS, and a variation of the chemical signature towards Na–HCO₃ compositions, with a parallel enrichment of F. The interacting matrixes are mainly rhyolites consisting of volcanic glass and only rare F-bearing accessory minerals (such as alkali amphibole). Comparing the abundance and the composition of the glassy groundmass with other mineral phases, it appears that the former stores most of the total F budget. This glassy material is extremely reactive, and its weathering products (i.e. fluvio/volcano-lacustrine sediments) further concentrate the fluoride. The interaction of these “weathered/reworked” volcanic products with water and carbon dioxide at high pH causes the release of fluoride into the interacting water. This mainly occurs by a process of base-exchange softening with the neo-formed clay minerals (i.e. Ca–Mg uptake by the aquifer matrix, with release of Na into the groundwater). This is plausibly the main enrichment mechanism that explains the high F content of the local groundwater, as evidenced by positive correlation between F, pH, and Na, and inverse correlation between F and Ca (Mg). Saturation indices (SI) have been calculated (using PHREEQC-2) for the different water groups, highlighting that the studied waters are undersaturated in fluorite. In these conditions, fluoride cannot precipitate as CaF₂, and so mobilizes freely without forming other complexes. These results have important implications for the development of new exploitation strategies and accurate planning of new drilling sites. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.