利用SV分量接收函数反演地壳横波速度结构

详细讨论了远震体波SV分量接收函数的特点及其在反演地壳S波速度结构中的优势．与径向接收函数类似，SV分量接收函数可通过对远震体波的SV分量直接反褶积P分量获得．研究分析表明:与径向接收函数相比，SV分量接收函数的振幅随震中距的变化更加稳定，波形简单且突出了对结构最敏感的PS转换波信息．理论数值实验显示:在反演地壳S波速度结构时，SV分量接收函数比径向接收函数具有更好的收敛性．作为实例，利用SV分量接收函数反演方法反演了海拉尔台下的S波速度结构．      

接收函数方法的研究综述

从远震体波波形数据中提取的P波接收函数和S波接收函数已经成为研究台站下方地壳上地幔速度间断面的最有效的方法之一.P波接收函数已经被广泛应用于获取地壳内部S波速度结构、地壳厚度及物质成分组成、地幔过渡带的厚度变化以及岩石圈地幔的间断面等,而S波接收函数是P波接收函数的一个很好的补充,因为在Moho和地幔过渡带之间的深度范围内,Sp转换波比来自浅部间断面的多次波到达早,而在P波接收函数中,相同深度范围的间断面的Ps转换波往往被来自浅部间断面的多次波干扰或者淹没,因此S波接收函数是目前获取岩石圈地幔深度范围内速度间断面结构(如Moho和LAB)的比较有效的方法,比面波观测具有更高的分辨率.本文详细阐述了P波接收函数和S波接收函数的方法原理以及在地壳上地幔速度间断面的研究中所采用的研究思路.     

接收函数复谱比的最大或然性估计及非线性反演

利用Ｓｈａｍｗａｙ和Ｄｅｒ的多道最大或然性反褶积原理，讨论了从单台三分量远震Ｐ波波形中分离接收函数径向与垂向分量复谱比的方法．根据Ｔａｒａｎｔｏｌａ的波形反演理论发展了接收函数复谱比的非线性反演方法．本文的反演方法除了拟合接收函数的复谱比之外，还需拟合时间域中接收函数垂向与径向分量的初至振幅比．合成波形数据反演的结果表明本文方法的反演结果不依赖初始模型．利用该方法和ＣＤＳＮ台网兰州台记录的宽频带远震Ｐ波波形数据，研究了该台站的接收函数及其随方位角的变化．接收函数复谱比非线性反演给出了该台站下方１４０ｋｍ深度岩石层的Ｓ波速度结构，得到兰州台下方地壳厚度为５４ｋｓ，岩石层厚度为９４ｋｍ．     

接收函数方法及其新的进展

远震P波波形数据中包含了大量地震台站下方地壳和上地幔速度间断面所产生的PS转换波及其多次反射波的信息，由此提取的接收函数是了解地壳上地幔速度细结构的重要步骤手段之一．最近几年，接收函数和面波联合反演方法获得了较大的成功，两种分别对两种波形进行拟合，对反演的速度结构提供了有效的约束．地震勘探中的一些成熟技术被引进接收函数的数据处理，使其可用于地壳和上地幔主要速度界面的侧向变化研究．为增强接收函数的信噪比，将不同事件提取的接收函数进行分类，按方位角进行叠加，可以反应一些速度界面的横向变化．接收函数方法可用于PS转换震相的剪切波分裂研究．泊松比是推测地球内部物质构成的有效参数之一．接收函数方法分离出的转换波为获取泊松比提供了一条行之有效的便捷途径，单一的地震台远震记录中转换波可以估计点位下面地壳Vp/Vs值．另一方面，一些最新的反演技术被引入接收函数的反演．例如，格子收索方法力图解决接收函数反演中的不稳定性和非唯一性问题，格子搜索设计比较容易地结合先验约束并保证搜索是全空间的，避免了使用任何初始模型．该方法能保证在格子的间隔和参数的限度内获得全局最小，该方法被用于中东、南美及临近地区的深部结构．相邻算法是基于计算几何的概念而构建的一种非线性反演方法．该方法避免了此前一些方法(遗传算法、模拟退火)的一些缺陷，如大量样本的舍弃、过多参数的引进等，该算法对我国的五大连池和腾冲火山区的深部结构的反演，以及滇西地区的深部动力学研究获得了良好的效果．     

P波和S波接收函数的贝叶斯联合反演

S波接收函数对于研究岩石圈速度结构具有重要价值. 本文利用合成地震图技术研究了S波接收函数的动力学特征. 在接收函数非线性复谱比反演方法的基础上,发展了基于贝叶斯理论的P波和S波接收函数的非线性联合反演方法. 结果表明:(1)适用于S波接收函数反演的震中距范围约为55°~80°,S波接收函数反演要求所用远震事件的震级大于5级; (2)与陡变的岩石圈底部界面(LAB)相比,梯度带类型LAB上生成的S_LP转换波相对较弱,台站下方的沉积盖层有助于相对增强S_LP震相; (3)由于S波接收函数径向分量不符合δ脉冲,不依赖于等效震源假定的三分量接收函数多道最大或然性反褶积方法更适合S波接收函数的估计;(4)数值检验的结果表明,在初始模型速度参数偏离真实模型20%的情况下,本文的方法能够预测300 km深度范围内的P波和S波速度结构;(5)观测数据的反演结果表明,由于P波接收函数低频分量相对不足,本文的联合反演方法对于大于100 km深度上地幔的S波速度结构约束相对较弱.     

基于深度学习的接收函数横波速度预测

远震接收函数中包含了大量台站下方速度间断面所产生的Ps转换波及其多次反射波的信息,被广泛应用于反演精细的地壳及上地幔速度结构.但当地壳结构复杂如存在沉积层或高速层时,转换波和多次波震相的到时和振幅发生变化,从而导致接收函数反演存在强非唯一性.深度学习作为一种高效的特征提取方法,能够建立接收函数与横波速度在空间上的映射关系.因此,利用全球模型数据以及高质量观测接收函数建立样本库,设计了利用接收函数预测横波速度的卷积神经网络.测试集结果表明,合成数据预测的横波速度与模型拟合程度较高,在实际数据检测中预测横波速度结果与全局反演的结果基本一致,预测的横波速度间断面与传统H-κ叠加结果也基本一致.利用该方法对琉球海沟地区海底地震数据的接收函数进行反演,获得了台站下方横波速度精细结构,进一步获得了该地区更为精细的地壳结构特征.测试实验及应用表明利用卷积神经网络预测台站下方横波速度的深度学习方法不仅计算效率高,且具有较高的可靠性.     

中国数字地震台网的接收函数及其非线性反演

利用中国数字台网（CDSN）记录的85个远震事件的宽频带P波波形数据和分离接收函数的最大或然性反褶方法,获得了CDSN台网10个台站不同方位的岩石层接收函数.利用这些台站不同方位的平均接收函数和非线性接收函数反演方法,获得了上述各台站下方100km深度范围内的岩石层S波速度结构.结果表明CDSN台网各台站下方的地壳厚度和岩石层速度结构存在明显的差别.其中,拉萨台下方的地壳厚度为66kin,壳幔界面较模糊,而余山台下方的地壳厚度为34km,壳慢界面两侧速度反差明显.刮用本文方法估计的地壳厚度与已有的结果基本一致,由于CDSN台网覆盖了中国大陆的各主要构造单元,本文的结果为研究中国大陆的岩石层S波速度结构及其横向变化提供了新的证据.     

利用宽频带远震体波波形研究青藏高原地壳上地幔速度结构的初步结果

利用根据中美合作研究青藏高原深部结构计划布设在青藏高原上的11个宽频带数字地震仪记录到的远震体波数据,采用接收函数(receiver function)反演的方法,对各台站下面地壳上地幔地震波速度结构进行了研究。台站接收函数是通过将三分向地震记录的两个水平分量旋转合成得到径向分量,然后在频率域除以垂直分量并变换回到时间域得到的,它仅与台站下面介质结构有关,而基本上与震源函数和传播路径无关。为压制噪声干扰,对来自同一方向上一定震中距范围内的远震记录得到的接收函数进行了叠加。采用分层弹性介质中弹性波传播矩阵理论,我们可以计算得到分层介质的理论接收函数以及它对各层弹性参数的偏导数,从而利用迭代线性反演可从观测接收函数得到台站下面的一维速度结构。本文给出了其中3个台,即温泉台、格尔木台和日喀则台的初步结果,它们分别位于高原的中部、北部和南部。从各台的接收函数中都可看到清晰的 Moho 面上的 P-S 转换波震相,其相对直达 P 波的走时延迟分别为:温泉台7.9s(东北方向结果),8.3s(东南方向结果);格尔木台8.2s;日喀则台9.0s,如此大的延迟表明高原地壳的巨厚.      

接收函数方法及研究进展

远震P波波形数据中包含了大量在台站下方地壳上地幔速度间断面所产生的P-S转换波及其多次反射波的信息，是研究台站下方局部区域S速度分布理想的震相，由此产生的接收函数方法是反演台站下方S波速度结构的有效手段。接收函数方法可以通过波形反演拟合接收函数的径向分量，对观测台站下方地球介质的S波速度结构进行估计，也可以通过偏移叠加获得的接收函数道集（地震剖面图）追踪速度间断面。这种方法避免了对天然地震震源及其附近结构混响效应等复杂因素的影响，对S波速度的垂向分布敏感，垂向分辨率高。由于宽频带流动地震台阵的发展，用此方法还可获得研究区域速度结构的横向变化，横向分辨能力主要取决于台站的间距。本文回顾20年来接收函数研究的进展，探讨了方法研究的发展趋势，介绍了对地壳-上地幔结构的部分研究结果。     

地壳分层各向异性介质接收函数及其粒子群反演

地壳不同深度介质的地震各向异性是研究地壳不同深度范围变形方式的重要依据.鉴于地壳介质的复杂性,如何从远震体波接收函数中提取不同深度的各向异性参数仍是一个有待深入研究的课题.在已有研究的基础上,本文利用广义反射-透射系数矩阵方法计算的合成地震图,研究了复杂地壳分层各向异性介质的接收函数随反方位角(back azimuth)变化及不同层位各向异性参数对接收函数波场的影响,为各向异性介质接收函数的解释提供了新的理论依据.通过引入粒子群优化理论,发展了分层各向异性介质接收函数全局反演算法.数值及观测数据的验证结果表明,在各向同性速度模型确定的前提下,我们的方法能够可靠地提取地壳分层各向异性参数;在反演中引入曲波变换去噪技术,对于正确解析不同层位的各向异性参数具有重要价值.     

Temporal variations and scaling of streamflow and baseflow and their nitrate-nitrogen concentrations and loads

The patterns of temporal variations of precipitation (P), streamflow (SF) and baseflow (BF) as well as their nitrate-nitrogen (nitrate) concentrations (C) and loads (L) from a long-term record (28 years) in the Raccoon River, Iowa, were analyzed using variogram and spectral analyses. The daily P is random but scaling may exist in the daily SF and BF with a possible break point in the scaling at about 18 days and 45 days, respectively. The nitrate concentrations and loads are shown to have a half-year cycle while daily P, SF, and BF have a one-year cycle. Furthermore, there may be a low-frequency cycle of 6–8 years in C. The power spectra of C and L in both SF and BF exhibit fractal 1/f scaling with two characteristic frequencies of half-year and one-year, and are fitted well with the spectrum of the gamma distribution. The nitrate input to SF and BF at the Raccoon watershed seems likely to be a white noise process superimposed on another process with a half-year and one-year cycle.     

燕郊等测点迁移优化与地磁观测研究

为优化地磁观测条件，开展了燕郊、夏垫等测点迁移工作；按照测眯迁移原则与实施技术方案，完整地收集并整理了地质构造与地球物理等方面的基本资料；进行了野外实地勘察，磁场梯度的测量，确定了新测点，在新老测点上进行了较长时间的地磁场对比观测；应用多种方法分析研究了地磁对比观测资料，结果表明，新老测点与有关测点的地磁变化具有良好的一致性，并得到了新老测点之间的地磁数据的按点差。     

东、南洞庭湖的径流、泥沙特征及冲淤规律

通过实地调查并对1957年以来水文、泥沙观测资料做系统分析和计算,探讨东、南洞庭湖出、入湖水量、沙量的年际和年内变化特征,以及长江下荆江段裁弯对湖区径流和泥沙的影响。提出了湖区泥沙汛淤枯冲的变化规律及水位升降与湖区泥沙冲淤的关系;论证了丰、平、枯年湖区淤积严重,面积日益缩小对径流的调节作用正在减弱。     

GIS与防灾减灾

地理信息系统(GIS)在我国已得到广泛的应用。在简要介绍G1S的发展历程、特点和发展趋势的基础上，重点从防灾减灾的角度介绍了GIS在气象灾害、地震灾害、地质灾害等领域的应用成果。     

Assessing Groundwater Depletion and Dynamics Using GRACE and InSAR: Potential and Limitations

In the last decade, remote sensing of the temporal variation of ground level and gravity has improved our understanding of groundwater dynamics and storage. Mass changes are measured by GRACE (Gravity Recovery and Climate Experiment) satellites, whereas ground deformation is measured by processing synthetic aperture radar satellites data using the InSAR (Interferometry of Synthetic Aperture Radar) techniques. Both methods are complementary and offer different sensitivities to aquifer system processes. GRACE is sensitive to mass changes over large spatial scales (more than 100,000 km²). As such, it fails in providing groundwater storage change estimates at local or regional scales relevant to most aquifer systems, and at which most groundwater management schemes are applied. However, InSAR measures ground displacement due to aquifer response to fluid‐pressure changes. InSAR applications to groundwater depletion assessments are limited to aquifer systems susceptible to measurable deformation. Furthermore, the inversion of InSAR‐derived displacement maps into volume of depleted groundwater storage (both reversible and largely irreversible) is confounded by vertical and horizontal variability of sediment compressibility. During the last decade, both techniques have shown increasing interest in the scientific community to complement available in situ observations where they are insufficient. In this review, we present the theoretical and conceptual bases of each method, and present idealized scenarios to highlight the potential benefits and challenges of combining these techniques to remotely assess groundwater storage changes and other aspects of the dynamics of aquifer systems.     

Fault and system stiffnesses and stick-slip phenomena

Summary This paper discusses the influence of system stiffness on the dynamic instability of fault surfaces under laboratory conditions for a number of test modes. In conjunction with shear load stiffness, the normal load stiffness, often neglected, is shown to have a considerable effect on the stick-slip process —its presence or absence and its characteristics. Also appropriate stiffnesses are suggested for an earthquake sequence modeled as a growing dislocation.