SH Wave Propagation in Viscoelastic Media

When comparing solutions for the propagation of SH waves in plane parallel layered elastic and viscoelastic (anelastic) media, one of the first things that becomes apparent is that in the elastic case the location of the saddle points required to obtain a high frequency approximation are located on the real p axis. This is true of the branch points also. In a viscoelastic medium this is not typical. The saddle point corresponding to an arrival lies in the first quadrant of the complex p-plane as do the branch points. Additionally, in the elastic case the saddle point and branch points lie on a straight line drawn through the origin (the positive real axis in the complex p-plane), while in the viscoelastic case this is generally not the case and the saddle point and branch points lie in such a manner as to indicate the degree of their complex values.In this paper simple SH reflected and transmitted particle displacement arrivals due to a point torque source at the surface in a viscoelastic medium composed of a layer over a half space will be considered. The path of steepest descent defining the saddle point in the first quadrant will be parameterized in terms of a real variable and the high frequency solutions and intermediate analytic results obtained will be used to formulate more specific constraints and observations regarding saddle point location relative to branch point locations in the complex p-plane.As saddle point determination for an arrival is, in general, the solution of a non-linear equation in two unknowns (the real and imaginary parts of the complex saddle point p ₀), which must be solved numerically, the use of analytical methods for investigating this problem type is somewhat limited.Numerical experimentation using well documented solution methods, such as Newton's method, was undertaken and some observations were made. Although fairly basic, they did provide for the design of algorithms for the computation of synthetic traces that displayed more efficient convergence and accuracy than those previously employed. This was the primary motivation for this work and the results from the SH problem may be used with minimal modifications to address the more complicated subject of coupled P-SV wave propagation in viscoelastic media.Another reason for revisiting a problem that has received some attention in the literature was to approach it in a fairly comprehensive manner so that a number of specific observations may be made regarding the location of the saddle point in the complex p-plane and to incorporate these into computer software. These have been found to result in more efficient algorithms for the SH wave propagation and a significant enhancement of the comparable software in the P-SV problem.     

关于大地电磁信号非最小相位性的讨论

由于信号的功率谱仅包含振幅信息而不包含相位信息,因此传统功率谱方法仅能处理最小相位信号,难于适用于非最小相位信号．本文对不同地区的实际大地电磁(MT)信号作了分析,结果表明,不同测点MT信号在Z平面上的零点分布并不全部位于单位圆之内,可以判定MT信号一般具有非最小相位的特点,因此有必要改进传统处理方法,发展新的资料处理方式．实际信号分析,理论和数值模型实验证明,信号的高阶统计量包含其相位信息,因此,高阶统计量方法能够有效处理具有非最小相位特征的MT信号．     

混合二维泊松过程的分解算法及其在提取地震丛集模式中的应用.

将一定范围内的地震数据假设为背景地震和丛集地震的叠加，并同时认为背景地震和丛集地震分别满足不同参数的二维泊松过程．通过引入N阶距离概念，将叠加的二维泊松过程转化为一维的混合密度函数，在对距离阶数进行选择的基础上，最终采用遗传算法进行混合密度分解，以达到提取地震丛集模式的效果．文中将该算法应用于我国西南地区松潘及龙陵主震前丛集地震的提取，并与C值的时间扫描结合，深化了这两次大地震前地震活动图象的认识．      

The algorithm of decomposing superim-posed 2-D Poisson processes and its applica-tion to the extracting earthquake clustering pattern

IntroductionClusteringearthquakesareusuallyconsideredasomensofstrongearthquakesorasignaloftectonicmovement.Thus,theyarenotonlyoneoftheprimaryevidencestopredictearthquakesbutalsoasignificantindicatortorecognizetectonicmovement(MEI,etal,1993;EarthquakePre-dictionandPreventionDepartmentofChinaSeismologicalBureau,1998).Ongeneralconditions,webelievethatclusteringearthquakesexistrelativelytobackgroundearthquakes,howtoeffectivelyseparateonefromtheotherbecomesthekeypointofextractingtheclusteringea…     

Spatial controls on the distribution and dynamics of a marginal snowpack in the Australian Alps

Seasonal snowpacks in marginal snow environments are typically warm and nearly isothermal, exhibiting high inter‐ and intra‐annual variability. Measurements of snow depth and snow water equivalent were made across a small subalpine catchment in the Australian Alps over two snow seasons in order to investigate the extent and implications of snowpack spatial variability in this marginal setting. The distribution and dynamics of the snowpack were found to be influenced by upwind terrain, vegetation, solar radiation, and slope. The role of upwind vegetation was quantified using a novel parameter based on gridded vegetation height. The elevation range of the catchment was relatively modest (185 m), and elevation impacted distribution but not dynamics. Two characteristic features of marginal snowpack behaviour are presented. Firstly, the evolution of the snowpack is described in terms of a relatively unstable accumulation state and a highly stable ablation state, as revealed by temporal variations in the mean and standard deviation of snow water equivalent. Secondly, the validity of partitioning the snow season into distinct accumulation and ablation phases is shown to be compromised in such a setting. Snow at the most marginal locations may undergo complete melt several times during a season and, even where snow cover is more persistent, ablation processes begin to have an effect on the distribution of the snowpack early in the season. Our results are consistent with previous research showing that individual point measurements are unable to fully represent the variability in the snowpack across a catchment, and we show that recognising and addressing this variability are particularly important for studies in marginal snow environments.     

Spatiotemporal variability in hydrochemistry of shallow groundwater in a small pre‐alpine catchment: The importance of landscape elements

Topography and landscape characteristics affect the storage and release of water and, thus, groundwater dynamics and chemistry. Quantification of catchment scale variability in groundwater chemistry and groundwater dynamics may therefore help to delineate different groundwater types and improve our understanding of which parts of the catchment contribute to streamflow. We sampled shallow groundwater from 34 to 47 wells and streamflow at seven locations in a 20‐ha steep mountainous catchment in the Swiss pre‐Alps, during nine baseflow snapshot campaigns. The spatial variability in electrical conductivity, stable water isotopic composition, and major and trace ion concentrations was large and for almost all parameters larger than the temporal variability. Concentrations of copper, zinc, and lead were highest at sites that were relatively dry, whereas concentrations of manganese and iron were highest at sites that had persistent shallow groundwater levels. The major cation and anion concentrations were only weakly correlated to individual topographic or hydrodynamic characteristics. However, we could distinguish four shallow groundwater types based on differences from the catchment average concentrations: riparian zone‐like groundwater, hillslopes and areas with small upslope contributing areas, deeper groundwater, and sites characterized by high magnesium and sulfate concentrations that likely reflect different bedrock material. Baseflow was not an equal mixture of the different groundwater types. For the majority of the campaigns, baseflow chemistry most strongly resembled riparian‐like groundwater for all but one subcatchment. However, the similarity to the hillslope‐type groundwater was larger shortly after snowmelt, reflecting differences in hydrologic connectivity. We expect that similar groundwater types can be found in other catchments with steep hillslopes and wet areas with shallow groundwater levels and recommend sampling of groundwater from all landscape elements to understand groundwater chemistry and groundwater contributions to streamflow.     

整点气温缺测的插补方法研究及其初步应用

长期连续完整的历史气温资料是震前气温异常判别研究的重要数据基础。本文考虑了参考站与缺测站之间的距离,建立改进的线性回归模型。利用该模型插补缺测和错误的气温整点值数据,在一定程度上解决了长期连续观测数据缺测的情况。通过对收集的唐山观测站气温整点值数据进行插补,并应用插补完整的数据分析研究了2012年5月28日唐山4.8级地震前兆异常。结果表明:①插补值与其前后观测值衔接吻合,插补后完整连续数据符合夏高冬低的年变规律;②插补误差在±0.5℃范围内的比例为60.2%,在±0.8℃范围内的比例为80.3%,其误差绝对值大于1.0℃的比例为9.6%,平均绝对误差为0.84℃,插补值与观测值的相关系数大部分在0.9以上;③从3月27日起出现增温异常,特别是震前2天增温幅度约8℃。     

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.     

多分量地震资料处理解释技术研究

随着油田勘探开发的不断深入,多波多分量地震资料的应用也越来越广泛.本文结合实际资料,对多波多分量地震资料的处理和解释方法进行了探讨和研究,主要针对多波资料的静校正、叠加成像、转换波资料的共转换点的求取方法等技术进行了分析研究,开发了综合利用多波资料进行地震属性和各向异性研究的方法,形成了一套较完善的多分量地震资料的处理、解释技术系列.经过在胜利油田罗家地区的初步应用,表明了预测结果与钻井资料相吻合.     

在SeisDps实时系统中用零点、极点计算ML的过程

在SeisDps地震处理软件的实时系统中,不是用系统的幅频特性计算震级ML,而是利用台站的零点、极点和增益,将记录的数字地震波形仿真成伍德-安德森地震仪记录的波形,计算具有3 大优点的ML震级,本文主要介绍用零点、极点和增益计算ML的过程。