Distributed jobs on EGEE Grid infrastructure for an Earth science application: moment tensor computation at the centroid of an earthquake

The Centroid Moment Tensor is a good summary of the fault rupture happening during an earthquake. We here describe a method for its determination in case of large earthquakes (with magnitude greater than 6.5) recorded world-widely by permanent digital seismic stations networks, as well as how deploying it on a Grid computing facility finally enabled its fast completion. Two different approaches of this process used over the last years are herein detailed and discussed. Some important impacts of this work are then highlighted and future trends are foreseen.     

频率域激电信号采集平台研究

为提高传统频率域激电信号接收机的通用性,设计了一种可采集和处理变频信号、双频信号和伪随机多频信号等多种频率域激电信号的采集平台。该平台以STM32微控制器为核心设计硬件电路,实现低通滤波、程控放大、模数转换以及数据通信等功能;在LabVIEW中设计上位机,进行参数设置、数字滤波、频谱分析、幅值测量、求值计算等处理,得到幅频率等目标参数,并具有波形显示与回放、数据存储与查看等功能;设计了阻容网络模拟实验,3种激电法幅频率实测值与理论值绝对误差均小于0.5%。实验表明,平台能够灵活选取不同方法进行采集和处理,测量精度较高,在资源勘查以及煤矿水害探测等方面具有良好的应用前景。      

Source parameters and focal mechanisms of local earthquakes: Single broadband observatory at ISM Dhanbad

A three-component broadband seismograph is in operation since January 2007 at the Indian School of Mines (ISM) campus. We have used the broadband seismograms of two local earthquakes M <3 recorded by this single station to illustrate its efficacy in understanding the source processes and tectonics in Dhanbad area. Source parameters and fault plane solutions are obtained through waveform inversion. It is observed that these two earthquakes occurred in the lower crust at a depth of 26 km by strike slip faulting. North-south compressional and east-west tensional stresses are dominant in the area, and the lower crust is the source area for the local earthquakes.     

Statistics for Modeling Heavy Tailed Distributions in Geology: Part I. Methodology

Many natural phenomena exhibit size distributions that are power laws or power law type distributions. Power laws are specific in the sense that they can exhibit extremely long or heavy tails. The largest event in a sample from such distribution usually dominates the underlying physical or generating process (floods, earthquakes, diamond sizes and values, incomes, insurance). Often, the practitioner is faced with the difficult problem of predicting values far beyond the highest sample value and designing his system either to profit from them, or to protect against extreme quantiles. In this paper, we present a novel approach to estimating such heavy tails. The estimation of tail characteristics such as the extreme value index, extreme quantiles, and percentiles (rare events) is shown to depend primarily on the number of extreme data that are used to model the tail. Because only the most extreme data are useful for studying tails, thresholds must be selected above which the data are modeled as power laws. The mean square error (MSE) is used to select such thresholds. A semiparametric bootstrap method is developed to study estimation bias and variance and to derive confidence limits. A simulation study is performed to assess the accuracy of these confidence limits. The overall methodology is applied to the Harvard Central Moment Tensor catalog of global earthquakes.     

TR method (TRM): A separation and stress inversion method for heterogeneous fault-slip data driven by Andersonian extensional and compressional stress regimes

The recently proposed TR method (TRM) which uses the slip preference of the faults to separate heterogeneous fault-slip data in extensional and compressional Andersonian stress regimes, is enhanced so as to determine stress tensors with the use of the Wallace–Bott slip criterion. Published natural fault-slip data from the extensional region of Tympaki, Crete, Greece and artificial fault-slip data modeled from the Chelungpu thrust, activated during the 1999 Chi–Chi earthquake in Taiwan, have been used as case studies. In the first case, the fault-slip data previously considered as homogeneous might actually be of heterogeneous origin as they determine two distinct stress tensors that both fit well with the neotectonic faulting deformation of the region. In the second case, where the fault-slip data belonging to three different subsets are of low diversity, the TRM succeeds in defining the driving stress tensors. The Misfit Angle minimization criterion can adequately separate the fault-slip data between two subsets when the percentage of the “Stress Tensor Discriminator Faults” is higher than approximately 70%.     

探地雷达波的相位参数及其应用

从探地雷达的实测波形出发,利用雷达反射波的特性并通过Ｈｉｌｂｅｒｔ变换导出计算相位的公式。在此基础上,借助球面波的表达式分析雷达波相位参数在均匀介质与非均匀介质听变化特征,并利用地基土层地质界面识别与地下土洞判识的实测了来进一步说明打位参数法的优点。     

Non-double-couple mechanisms of microearthquakes induced during the 2000 injection experiment at the KTB site, Germany: A result of tensile faulting or anisotropy of a rock?

Moment tensors of microearthquakes induced during the 2000 injection experiment at the KTB deep drilling borehole at a depth level of 5.4 km are studied. A family of 37 most reliable moment tensors contains significant non-double-couple (non-DC) components. The DC is on average 60% and the non-DC is 40%. Fault plane solutions computed from the DC part show preferred strike-slip mechanisms with small normal or reverse components. A predominant azimuth of P and T axes is in the range of N320°–340°E and of N230°–250°E, respectively. The non-DC components contain both the isotropic (ISO) and compensated linear vector dipole (CLVD) components. The mean value of ISO is 1.5%, the mean value of CLVD is − 5.7%. The predominantly negative CLVD components are inconsistent with the concept of the non-DC mechanisms as a result of tensile faulting due to fluid injection into the rock. The main origin of the non-DC components is probably anisotropy in the focal area. The other origins are errors produced by mismodelling of the medium when calculating the Green functions, and numerical errors produced by noise and limitations of input data. Adopting four alternative models of anisotropy obtained by other seismic measurements at the KTB, we have employed the non-DC components for estimating an optimum orientation of anisotropy in the focal area. The optimum orientation of the symmetry plane of anisotropy is nearly vertical with a strike of N335°–340°E. This strike coincides well with the strike of 330° typical for many major lithological units and faults and with the orientation of the transversely isotropic model inferred by other authors. After removing the anisotropy effects from the moment tensors by calculating the source tensors, the distribution of ISO is significantly narrowed. This indicates predominantly shear, but not tensile faulting.     

Seismological study of the central Ecuadorian margin: Evidence of upper plate deformation

A seismic study of a segment of the convergent margin of Ecuador is presented. During the SISTEUR campaign a network of 24 Ocean Bottom Seismometers (OBS) was deployed on the Carnegie Ridge, one line along the main axes of the ridge and two lines across the strike of the edge of the ridge, during one month. This marine network was complemented with a land network of 20 stations distributed in two lines: one parallel to the margin and the other perpendicular to it.The seismic event recorded by these networks, were located using different crustal models defined from the wide-angle seismic data modeling. Relative location techniques were used to improve earthquake locations. Seismogram waveform modeling allowed us to constrain hypocentral location for events farther than ～50 km from the network. This modeling also provided additional information to constrain the focal mechanisms of these events. The upper limit of the Interplate Seismogenic Zone (ISZ) is estimated to be at a 10 km depth in the region. The background seismic activity of the upper plate provided new insights:1) A seismic cluster that reaches the base of the overriding plate is linked to the Jipijapa-Portoviejo fault. The reactivation of this Quaternary fault is confirmed by focal mechanisms that provide rupture planes parallel to its superficial projection (N10°–N25°).2) The focal mechanisms presented in this study are compatible with a homogeneous regional stress field corresponding to an E–W to ESE–WNW compression and an NNE–SSW extension. The presence of strike-slip deformation, with a reverse component, corresponds to the NNE escape of the North Andean Block. Normal faulting accommodating this movement suggests that this part of the North Andean Block cannot be considered as a rigid block.     

基于地层倾斜规律的防斜打直钻井技术研究

地层自然倾斜、钻具组合、钻井参数是影响钻井井斜方位变化的主要因素。根据地层层面空间网格数据,建立了地层倾角、倾向的计算方法,并针对典型防斜打直钻具组合和钻井参数之间的关系进行了力学分析。应用结果表明,井眼实钻轨迹与设计轨迹基本相同,形成了利用地层自然倾斜规律,在钻井自然中靶的同时,释放钻压提高钻速的防斜打直优快钻井技术,对于高陡构造油气藏的高效开发具有重要的理论意义和工程实用价值。     

A unified earthquake catalogue for South Asia covering the period 1900–2014

Seismicity analysis is very much pertinent for Indian subcontinent and its adjoining region which is seismically active including many great earthquakes associated with collision and subduction tectonics in the northern, north-eastern part of the subcontinent and in the Andaman and Nicobar Island. An earthquake catalogue has been generated for South Asia covering the period 1900–2014 by compiling the records of earthquake occurrences from International Seismological Center, Global Centroid Moment Tensor (GCMT), US Geological Survey, India Meteorological Department and published literature. The uniform magnitude scaling in moment magnitude M _W,GCMT is achieved through connecting relationships between different magnitude types. These relationships are derived by orthogonal standard regression analysis on available data pairs. The derived relationships have been compared with the existing equations already reported by others. The catalogue is subsequently subjected to a seismicity declustering algorithm to identify the foreshocks, main-shocks and aftershocks. The catalogue thus compiled is envisaged to be a useful resource for seismotectonic and seismic hazard studies in the region.

     

Toward real-time regional earthquake simulation II: Real-time Online earthquake Simulation (ROS) of Taiwan earthquakes

We developed a Real-time Online earthquake Simulation system (ROS) to simulate regional earthquakes in Taiwan. The ROS uses a centroid moment tensor solution of seismic events from a Real-time Moment Tensor monitoring system (RMT), which provides all the point source parameters including the event origin time, hypocentral location, moment magnitude and focal mechanism within 2 min after the occurrence of an earthquake. Then, all of the source parameters are automatically forwarded to the ROS to perform an earthquake simulation, which is based on a spectral-element method (SEM). A new island-wide, high resolution SEM mesh model is developed for the whole Taiwan in this study. We have improved SEM mesh quality by introducing a thin high-resolution mesh layer near the surface to accommodate steep and rapidly varying topography. The mesh for the shallow sedimentary basin is adjusted to reflect its complex geometry and sharp lateral velocity contrasts. The grid resolution at the surface is about 545 m, which is sufficient to resolve topography and tomography data for simulations accurate up to 1.0 Hz. The ROS is also an infrastructural service, making online earthquake simulation feasible. Users can conduct their own earthquake simulation by providing a set of source parameters through the ROS webpage. For visualization, a ShakeMovie and ShakeMap are produced during the simulation. The time needed for one event is roughly 3 min for a 70 s ground motion simulation. The ROS is operated online at the Institute of Earth Sciences, Academia Sinica (http://ros.earth.sinica.edu.tw/). Our long-term goal for the ROS system is to contribute to public earth science outreach and to realize seismic ground motion prediction in real-time.     

Reviews in Garnet—Biotite Geothermometer：Its Versions,Inherent Problems of Accuracy and Precision,and Perspective for Further Research

This paper briefly introduces twenty-one versions of garnet-biotite Fe-Mg exchange geothermometry,points out the sources of systematic errors inherent in geothermometry.The sources lie in that:(1) ideal Fe-Mg mixing in garnet and biotite is assumed;(2)non-ideal properties of bot garnet and biotite are not considered completely;(3) minor elements-Mn, Ca in garnet,and Al^Ⅵ,Ti,Mn in biotite are almost neglected;(4)effects of pressure on equilibrium,although little,are almost not taken into consideration,and experiments were conducted at one fixed pressure;(5)the data used for egression analysis is too scarce in quantity;(6)internal consistency of the data used for empirical calibration is not fully guaranteed;(7) the authors take energy parameters W‘s(Margule‘s parametrs)independent of pressure and temperature,though the former are actually the functions of the latter items;and (8)the most important composition-activity relationship is not consistent with the mineral compositions used for calibration.Furthermore,it is believed that empirical calibrations cannot avoid the following problems;(1) effect of retrograde metamorphism;(2)inconsistent data set;(3)effect of pressure on temperature estimation;and (4)effect of Fe^3 calculation.Exact experimental work should be done in order to accurately calibrate the geotermometer in which the Margule‘s parameters should be treated as the functions of pressure,temperature,and mineral compositions.     

目标函数对马斯京根法河道流量演算精度的影响

在用马斯京根法进行河道流量演算时,由于传统的试算法在精度和客观性上的欠缺,目前广泛使用最小二乘法来进行优化计算.在应用最小二乘法时,发现选择不同的目标函数会对最终的流量计算结果的精度产生影响.因此,本文应用了两种目标函数:河槽蓄量误差最小和出流量误差最小,推导了它们在最小二乘意义上的流量演进参数解析式,进而研究了对流量计算精度的影响.对3场洪水过程的模拟结果表明,以出流量误差最小为目标函数所获得的流量计算精度更高:与河槽蓄量误差最小相比,相对平均绝对误差分别降低了4%,25%和25%,说明使用出流量误差最小作为优化的目标函数更为有效.     

A logging method for estimating permeability, velocity of second compressional wave, and electroacoustic constant in electrolyte-saturated porous formations

A new algorithm is suggested for joint calculation of three key parameters of reservoir rocks: permeability, electroacoustic constant, and velocity of second bulk compressional wave (slow wave). The algorithm is advantageous in the use of Stoneley waves which are most stable in acoustoelectric log waveform data.     

Tsunamis in and near Greece and their relation to the earthquake focal mechanisms

The major earthquake-induced tsunamis reliable known to have occurred in and near Greece since antiquity are considered in the light of the recently obtained reliable data on the mechanisms and focal depths of the earthquakes occurring here. (The earthquake data concern the major shocks of the period 1962–1986.) First, concise information is given on the most devastating tsunamis. Then the relation between the (estimated) maximum tsunami intensity and the earthquake parameters (mechanism and focal depth) is examined. It is revealed that the most devastating tsunamis took place in areas (such as the western part of the Corinthiakos Gulf, the Maliakos Gulf, and the southern Aegean Sea) where earthquakes are due to shallow normal faulting. Other major tsunamis were nucleated along the convex side of the Hellenic arc, characterized by shallow thrust earthquakes. It is probably somewhere there (most likely south of Crete) that the region's largest known tsunami occurred in AD 365, claiming many lives and causing extensive devastation in the entire eastern Mediterranean. Such big tsunamis seem to have a return period of well over 1000 years and can be generated by large shallow earthquakes associated with thrust faulting beneath the Hellenic trench, where the African plate subduces under the Euroasian plate. Lesser tsunamis are known in the northernmost part of the Aegean Sea and in the Sea of Marmara, where strike-slip faulting is observed. Finally, an attempt is made to combine the tsunami and earthquake data into a map of the region's main tsunamigenic zones (areas of the sea bed believed responsible for past tsunamis and expected to nucleate tsunamis in the future).     

Source parameters of March 10 and September 13, 2007, United Arab Emirates earthquakes

On March 10 and September 13, 2007 two earthquakes with moment magnitudes 3.66 and 3.94, respectively, occurred in the eastern part of the United Arab Emirates (UAE). The two events were widely felt in the northern Emirates and Oman and were accompanied by a few aftershocks. Ground motions from these events were well recorded by the broadband stations of Dubai (UAE) and Oman seismological networks and provide an excellent opportunity to study the tectonic process and present day stress field acting in this area. In this study, we report the focal mechanisms of the two main shocks by two methods: first motion polarities and regional waveform moment tensor inversion. Our results indicate nearly pure normal faulting mechanisms with a slight strike slip component. We associated the fault plane trending NNE–SSW with a suggested fault along the extension of the faults bounded Bani Hamid area. The seismicity distribution between two earthquake sequences reveals a noticeable gap that may be a site of a future event. The source parameters (seismic moment, moment magnitude, fault radius, stress drop and displacement across the fault) were also estimated from displacement spectra. The moment magnitudes were very consistent with waveform inversion. The recent deployment of seismic networks in Dubai and Oman reveals tectonic activity in the northern Oman Mountains that was previously unknown. Continued observation and analysis will allow for characterization of seismicity and assessment of seismic hazard in the region.     

最大有效力矩准则的理论拓展

在分析"导致变形带内先存面理或层理发生转动的最大有效力矩与先存面理或层理方向有关"的基础上, 对最大有效力矩准则(M_eff=0.5(σ₁-σ₃) Lsin2αsinα)进行理论上的拓展, 提出了可以判定任意方向先存面理最大有效力矩的准则——泛最大有效力矩准则(M_G-eff=0.5(σ₁-σ₃)Lsin2αsin(α-θ)), 其中当先存面理与最大主压应力(σ₁)平行时, 则成为最大有效力矩准则。该准则的理论分析表明:①当先存面理与σ₁平行时, 在σ₁左右两侧±54.7°方向出现2个有效力矩的最大值, 形成共轭的变形带, 钝角(109.4°)对着σ₁方向; ②当先存面理与σ₁斜交时, 在σ₁的另一侧出现1个有效力矩的最大值, 从而只出现一个方向的变形带, 并随着先存面理偏离σ₁方向, 变形带与σ₁的夹角逐渐减小(从θ=0°时的54.7°, 减小到θ=90°时的35.3°), 而与先存面理之间的夹角逐渐增大(从θ=0°时的54.7°, 增加到θ=90°时的125.3°); ③当先存面理与σ₁垂直时, 在σ₁左右两侧± 35.3°方向出现2个有效力矩的最大值, 也形成共轭的变形带, 但锐角(70.6°)对着σ₁方向。在主应变平面上变形带与先存面理方向及变形带剪切方向(左旋或右旋)已知的情况下, 可以确定最大主压应力方向。泛最大有效力矩准则克服了最大有效力矩准则与滑移线理论不相容的问题, 可以解释大多膝褶带非共轭发育等多种现象, 预期在韧性变形域中具有广阔的应用前景。      

时间域瞬变电磁系统的应用软件开发与利用

针对时间域脉冲瞬变电磁仪的PROTEM系统,结合其工作特点编制了一套系列应用软件,功能包括空间各点一次激励磁场强度计算、野外质量检查各种误差值统计计算、野外观测数据预处理及各种转换参量计算。应用该软件对某金属矿区的瞬变电磁资料进行了处理,结果表明程序运行有效、可靠,取得了较好的地质效果。     

Monte-Carlo Method for Coalbed Methane Resource Assessment in Key Coal Mining Areas of China

Monte-Carlo method is used for estimating coalbed methane (CBM) resources in key coal mining areas of China. Monte-Carlo method is shown to be superior to the traditional volumetric method with constant parameters in the calculation of CBM resources. The focus of the article is to introduce the main algorithm and the realization of functions estimated by Monte-Carlo method, including selection of parameters, determination of distribution function, generation of pseudo-random numbers, and evaluation of the parameters corresponding to pseudo-random numbers. A specified software on the basis of Monte-Carlo method is developed using Visual C++ for the assessment of the CBM resources. A case study shows that calculation results using Monte-Carlo method have smaller error range in comparison with those using volumetric method.