A Possible Detection of the 26 December 2004 Great Sumatra-Andaman Islands Earthquake with Solution Products of the International GNSS Service

The main goal of this work is to critically review the IGS solution products and Precise Point Positioning (PPP) in order to demonstrate their potential to contribute to studies of large earthquakes such as the one that devastated Southeast Asia on December 26th, 2004. In view of a possible detection of the Mw 9.0 Sumatra-Andaman Islands Earthquake of December 26, 2004, position solutions, ranging from intervals of years to one second, of four International GNSS Service (IGS) stations within 3000 km of the epicenter were examined. The IGS combined, cumulative solution product (IGS04P51), consisting of epoch and station velocity solutions and based on data spans of several years prior to the earthquake, was used as a reference. Four IGS combined weekly position solutions (igs04P1301-4), two weeks before and after the earthquake, were utilized for the weekly solution resolution. PPP static and kinematic solutions with IGS Final combined orbits and clocks were used for the mean daily and instantaneous 5-min and 1-sec epoch solutions, respectively. The most significant changes, detected by both weekly and daily solutions occurred in longitude. The nearest IGS station ntus, about 1000 km east of the epicenter, moved westward about 15 mm, while the more distant Indian station iisc (∼ 2300 km NW from the epicenter), shifted about 15 mm eastward. In spite of position errors caused by interpolation of the 5-min IGS clocks, the 1-sec solutions, based on separate data sets, available only for two stations (iisc, dgar), still showed seismic surface waves, in particular at the Indian station iisc. Precise daily IGS combined polar motion and length-of-day products, after correcting for the atmospheric effects, also likely detected, statistically significant, anomalistic excitations on December 26, 2004 that could be caused by this great earthquake.     

Stress‐dependent elastic wave velocity of microfractured sandstone

A model for the stress‐dependent elastic wave velocity response of fractured rock mass is proposed based on experimental evidence of stress‐dependent fracture normal and shear stiffness. Previously proposed models and previous experimental studies on stress‐dependent fracture stiffness have been reviewed to provide a basis for the new model. Most of the existing stress‐dependent elastic wave velocity models are empirical, with model parameters that do not have clear physical meanings. To propose the new model, the rock mass is assumed to have randomly oriented microscopic fractures. In addition, the characteristic length of microfractures is assumed to be sufficiently short compared to the rock mass dimensions. The macroscopic stress‐dependent elastic wave velocity response is assumed to be attributed to the stress dependency of fracture stiffness. The stress‐dependent fracture normal stiffness is defined as a generalized power law function of effective normal stress, which is a modification of the Goodman's model. On the other hand, the stress dependency of fracture shear stiffness is modeled as a linear function of normal stress based on experimental data. Ultrasonic wave velocity responses of a dry core sample of Berea sandstone were tested at effective stresses ranging from 2 to 55 MPa. Visual observation of thin sections obtained from the Berea sandstone confirms that the assumptions made for microstructure of rock mass model are appropriate. It is shown that the model can describe the stress‐dependent ultrasonic wave velocity responses of dry Berea sandstone with a set of reasonable material parameter values. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

Modeling discontinuous rock mass based on smoothed finite element method

This paper aims at developing a method for modeling rock mass with preexisting multiple discontinuities within the framework of the smoothed finite element method (SFEM). The discontinuity is simulated by an interface element with zero thickness, the stiffness matrix of which are derived explicitly based on the SFEM. An elastic damage constitutive relation with residual strength is introduced in order to describe the nonlinear mechanical behavior of the discontinuities. The computation codes of the present method were developed. The present method has been verified to be a sound approach for modeling discontinuous rock mass, inheriting the advantages of the SFEM.     

基于重力全张量数据的Parker-Oldenburg算法研究

相比重力数据,重力张量数据通常包含更多的异常信息。本文根据重力数据与重力张量数据的关系,利用位场转化技术,将重力张量数据应用于传统的Parker-Oldenburg密度界面反演算法中。通过模型试验,证明了在网格间距较大或者数据存在一定噪音时,使用本文算法进行反演能得到更好的效果。实验结果说明利用重力张量数据可以有效地提高密度界面反演的分辨率。     

大同地区3次地震序列构造应力张量变化特征的研究

利用1989~2002年间3次大同地震序列中共计700多次中小地震的震源机制解资料,应用Gephart(1990)的应力张量反演方法研究了这3次地震序列的构造应力张量的总体变化特征和时序变化特征,研究发现3次主震发生前震源及附近地区的构造应力作用较强,主震发生前后,应力方向存在细节差异,但是最大主压应力方向与华北地区的构造应力场方向基本一致,只有1999年震前阶段的最大主压应力方向为226°(SW向),分析认为这可能是华北地区构造应力场与大同地区局部构造应力场相互作用的结果。     

Large-time response of ice cover to a load moving along a frozen channel

Unsteady response of an ice cover in a channel with vertical walls is studied for large times. The ice deflection is caused by a load moving along the frozen channel at a constant speed. The ice cover is modelled as a thin elastic plate clamped to the walls of the channel. The time-dependent problem is solved by using the Fourier transform along the channel and the method of separating variables. In the system moving along the channel together with the load, the large-time deflection of the ice cover consists of steady deflection and standing waves in front and behind the load. The number of waves, their frequencies and wavenumbers depend on the speed of the load and the values of the critical speeds for the channel. The number of the waves and their amplitudes are calculated for a given load and its speed. The maximum stress in the ice as a function of the load speed is estimated.     

Effective medium modeling to assess gas hydrate and free-gas evident from the velocity structure in the Makran accretionary prism, offshore Pakistan

Seismic properties of sediments are strongly influenced by pore fluids. Stiffness of unconsolidated marine sediment increases with the presence of gas hydrate and decreases with the presence of gas. A strong bottom-simulating reflector (BSR) observed on a seismic profile in the Makran accretionary prism, offshore Pakistan, indicates the presence of gas hydrate and free-gas across the BSR. Elastic properties of gas depend largely on pressure and temperature. We, therefore, first determine the elastic modulus of gas at pressure and temperature calculated at the BSR depth in the study region. The interval velocities derived from the seismic data are interpreted by the effective medium theory, which is a combination of self-consistent approximation and differential effective medium theories, together with a smoothing approximation, for assessment of gas hydrate and free-gas. The results show the saturations of gas hydrate and free-gas as 22 and 2.4% of pore space, respectively, across the BSR.     

旋转椭球面上的应变与转动张量表达

以旋转椭球体面上某点为原点建立一个大地坐标单位活动坐标架. 通过平移, 使活动坐标架的原点与以椭球中心为原点的笛卡尔单位标架的原点相重合. 然后再通过两次标架旋转, 使活动坐标架与笛卡尔单位标架完全重合. 本文给出了使两个单位标架相重合的转换关系式, 以及该点位移在两个单位标架中的坐标转换式; 在此基础上, 考虑该点的位移及活动坐标架皆为该点大地坐标的函数, 经复杂推导, 分别给出了该点位移向量的微分在大地坐标系中的分量以及该点分别沿坐标曲线的弧微分表达式, 继而导出了该点的位移梯度矩阵; 最后推导出了椭球坐标系的应变张量与转动张量表达式, 并对转动张量的几何含义进行了较详细的解释, 且采用曲面理论对球面与椭球面的应变张量间的内在关系进行了讨论.