大气重力格林函数的理论计算

重力测量中需要扣除大气的影响.大气负荷对重力测量的影响可以分为大气质量变化引起的直接效应和大气负荷引起的地球变形带来的间接效应,大气负荷对重力观测值的直接影响,相对于间接效应量级较大.本文从理论上研究了大气负荷对重力观测的直接影响,仿照Farrell定义的负荷格林函数,引入大气重力格林函数,用来表示大气压变化对于重力观测的直接引力影响.在前人的基础上,本文采用了更为精细的大气模型,考虑大气温度随高程的变化,用离散褶积的方法求得了大气重力格林函数的理论值.实际计算时还要考虑地表温度、台站高程、周围地形等因素的影响,本文讨论了这些因素对大气重力格林函数的影响.考虑地表温度、台站高程、地形改正等各种影响因素以及地球变形引起的间接效应后,对台站周围区域积分即可求得大气变化引起的理论重力信号.     

大气重力信号的理论计算及其检测

基于标准大气定律和大气圆柱体分布模型,本文引进了大气重力格林函数,用离散格积方法求得了大气对重力场观测的影响,对台站高程、周围地形和地表温度变化等因素的影响问题进行了讨论．结果说明台站近区气压变化是大气重力信号的主要贡献者,考虑大气质量负荷引起的弹性地球形变效应后,对距台站0．5°的区域积分获得的大气重力导纳值为-0．3603μGal／hPa,占全球大气变化引起的总信号的90％以上,这一理论模型结果与超导重力仪实测结果相吻合,并能较有效地用于消除重力观测中的气压干扰成分．     

地球物理场观测中的大气效应问题研究

作用在地表的大气质量不仅产生引力位的扰动，而且会引起各种地球物理场的附加效应，本文综述了近年来国内外同行在研究地表位移、重力、倾斜和应变固体潮观测中实施的大气效应改正问题，包括利用台站气压资料建立的实测模型，建立大气重力格林函数，利用全球和区域气象数据计算大气效应，建立大气改正的误差模型等，最后文章讨论了今后研究展望。     

旋转弹性椭球地球模型的固体潮理论值计算

根据Wahr 1981年提出的理论,导出了计算旋转弹性地球模型的重力固体潮、地倾斜固体潮和地面应变固体潮的公式,并在此基础上编写出相应的计算程序.为了显示旋转和扁度对地球模型的重力固体潮、地倾斜固体潮和地面应变固体潮的影响,计算了东经120°不同纬度处的旋转弹性椭球地球模型(1066A模型)和G-B地球模型的重力固体潮、地倾斜固体潮和地面应变固体潮.计算结果表明,旋转和扁度对重力固体潮、地倾斜固体潮和地面应变固体潮的最大幅度分别为1.4×10-8 m/s2、0.2ms和0.5×10-9.      

利用背景噪声互相关观测近对跖点台站间体波格林函数

本文利用秦岭周边和南美地区的部分台站2012—2014年记录的地震事件垂直分量数据,获取了166个台站对之间的互相关函数,从中观测到了近对跖点台站之间的体波格林函数.在此基础上,我们计算了理论地震图的互相关,其结果与实际数据获得的结果很相似,推测格林函数的能量来源应主要是大地震引起的地球自由振荡.更进一步的研究表明,不同的震源机制对于互相关结果中不同震相的振幅有较大影响,走滑断层引起的地震并不能产生较好的互相关结果.通过对比不同位置台站对的互相关结果,本文认为互相关得到的体波格林函数包含了传播路径信息.     

云南地区的背景弹性波场分析

利用云南省地震台网44个台站记录的2008年1-9月连续波形数据进行互相关计算,得到了台站间的格林函数,并获取对应的频散曲线,据此分析了该地区的区域背景弹性波场的来源及分布.研究发现,该地区的区域背景弹性波场有着明显的方向性,15s信号的总体传播方向是从东南向西北,也就是说主要来自云南省的东南方向,据此推测区域背景弹性...     

中国地壳运动观测网络基准站倾斜固体潮观测中的海潮负荷信号改正问题

利用最新的全球海潮模型（Csr3.0,Fes95.2,Tpxo2,Csr4.0和Scw80)和中国近海海潮资料，基于标准地球模型负荷格林函数，采用Agnew的积分格林函数方法研究了倾斜固体潮观测中的海潮负荷效应。文章计算了中国地壳运动观测网络25个基准站8个主要潮流的倾斜负荷，数值结果说明对某一潮波而言，倾斜负荷振幅达10^-9弧度，沿海地区达10^-8弧度或更多。文章构制了北京和上海等8个主要台站的倾斜负荷随时间变化。     

噪声地震学方法及其应用

基于背景噪声的地震方法发展迅速,已广泛应用于全球和区域地球内部结构研究、浅地表地质调查及油气田勘探开发.本文简要介绍了背景噪声的来源,回顾了噪声地震学的发展历程.给出了基于背景噪声的全波场和面波格林函数恢复的公式,较为详细综述了噪声源的分布和记录台站间距对格林函数恢复的影响.讨论了两台站互相关法和空间自相关法获取面波频散特性的区别与理论连接.对基于噪声的面波层析成像法、程函方程层析成像法、空间自相关法的原理进行了总结.介绍了噪声地震学方法在各领域特别是浅地表方面的应用现状.最后简要展望了噪声地震学的发展前景.     

大气负载引起地表位移:对重力和基线测量的影响

大气质量对地壳加载并使之形变.通过将逐日的全球气压数据和质量负载的格林函数进行摺积求和,我们怡算大气负载(包括那些与短期大气扰动有关的大气负载)对地表刚点定位测量和地表重力观测的时间依赖关系.我们茸出了无海洋地球和具有一负压海洋(inverted barometer ocean)的地球的响应.负压海洋显著地影响近海岸台站的负载响应.峰峰垂直位移常常是15一20mm,伴随有3一6卜Gal的扰动.基线变化可能大到20mm或更大.这种扰动在高纬度区和冬季达到最大.这些振幅的变化与Rabbel和Zs。hau(1985)的结果一致.他们用距某一点半径有关的高斯函数模拟了气压扰动.使用台距约I00okm以内测点的现有气压数据就足以计算地壳的形变.单是局部气压数据是不够的.Rabbel和Zsc五-au为使用局部气压和区域平均气压算得的这些位移位所假设的改正值,对正好在内地的测点证明是精确的,但一般地说在距海岸几百公里之内,改正值是不适当的.     

球体位错理论计算程序的总体设计与具体实现

本文详细介绍了球体位错理论计算程序的总体设计思想、 各类配套文件的具体内涵以及各类输出文件的物理含义, 同时介绍了程序的使用方法和注意事项, 为读者独立使用该程序提供参考。 球体位错理论计算程序主要由三部分组成: ① 位错格林函数计算程序, 基于具体的球对称地球模型提供离散的二维格林函数数值框架; ② 积分计算程序, 对离散的格林函数数值框架进行双二次样条插值运算, 并对四类独立位错源对应的格林函数进行适当组合, 从而计算出任意位置任意类型震源在地表产生的同震变形(含位移、 应变、 重力变化和大地水准面变化); ③ 辅助文件, 用于提供发震断层模型和计算点位信息。 一般情况下, 读者不需要理解位错格林函数计算程序和积分计算程序, 只需要对辅助文件提供的信息进行针对性改动, 就可以计算目标地震在目标观测站引起的同震变形。     

Loading effects in Metsähovi from the atmosphere and the Baltic Sea

Loading by atmosphere and by the Baltic Sea cause gravity change at Metsähovi, located 15 km from the open sea. Gravity is changed by both the Newtonian attraction of the loading mass and by the crustal deformation. We have performed loading calculations using appropriate Green's function for both gravity and deformation, for both atmospheric and Baltic loading. The loading by atmosphere has been computed using a detailed surface pressure field from high resolution limited area model (HIRLAM) for north Europe up to 10° distances. Baltic Sea level is modelled using tide gauge records. Calculations show that 1 m of uniform layer of water corresponds to 31 nm s⁻² in gravity and −11 mm in height. Modelled loading is compared with observations of the superconducting gravimeter T020 for years 1994–2002. The combination of HIRLAM and a tide gauge record decreases RMS of gravity residuals by 14% compared to single admittance in air pressure corrections without sea level data. Regression of gravity residuals on the tide gauge record at Helsinki (at 30 km distance) gives a gravity effect of 26 nm s⁻² m⁻¹ for Baltic loading.The gravity station is co-located with a permanent GPS station. We have also associated the loading effects of the atmosphere and of the Baltic Sea with temporal height variations. The range of modelled vertical motion due to air pressure was 46 mm and that due to sea level 18 mm. The total range was 38 mm. The effects of the Baltic Sea and of the atmosphere partly cancel each other, since at longer periods the inverse barometer assumption is valid. Regression of the modelled height on local air pressure gives −0.37 mm hPa⁻¹, corresponding approximately to width 6° for pressure system.We have tested the models using one year of daily GPS data. Multilinear regression on local air pressure and sea level in Helsinki gives the coefficient −0.34 mm hPa⁻¹ for pressure, and −11 mm m⁻¹ for sea level. These match model values. Loading by air pressure and Baltic Sea explains nearly 40% of the variance of daily GPS height solutions.     

Boundary crossing and non‐linear theory in earth‐system sciences – a proof of concept based on tsunami and post‐eruption scenarios on Java Island,Indonesia

Although volcanic eruptions are well‐known to be the trigger of some weather and climatic changes, land‐cover changes by pyroclastic‐flows and lahars do not get this recognition, neither do major hazards such as tsunami. These two earth processes are even lesser considered as being able to modify other earth processes they are not directly connected to, such as landslides or river discharge in non‐connected basins more than a hundred kilometres away. In this contribution the authors argue that these ideas are mainly driven by the process of being ‘educated’ in a single academic discipline and once put to the test interactions and retroactions between earth processes and atmospheric processes are far more reaching than commonly thought. For this study, the site of Java Island (Indonesia) was chosen to conduct (1) an analysis of a major tsunami impact – in the same area as the 2006 Java tsunami and (2) an analysis of the post‐eruption impacts of Merapi Volcano after a major eruption – excluding any ejecta in the atmosphere for the sake of the demonstration. The atmospheric feedback simulations were conducted using the regional climate model (RegCM‐4) with calibration from weather stations in Java Island. As a result, both simulations have proven that large scale deposits of pyroclasts (not introducing the ejectas sent in the atmosphere) and tsunamis can have outstanding impacts on the atmospheric situation and the bio‐geomorphologic evolution of the landscape in the following weeks to months. Interestingly enough these impacts are not limited to the area impacted by the earth process and the effect are not linear in time as they work following thresholds. These rainfalls ‘tele‐impacts’ are important enough to, in turn, modify earth‐surface processes in areas remote from the original phenomenon. This system acts in the same manner as a famous butterfly in Africa that could trigger a hurricane on the other side of the Atlantic Ocean. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling atmospheric and induced non-tidal oceanic loading contributions to surface gravity and tilt measurements

We investigate the contribution of atmospheric and its induced non-tidal oceanic loading effects on surface time-varying gravity and tilt measurements for several stations in Western Europe. The ocean response to pressure forcing can be modelled accordingly to the inverted barometer, i.e. assuming that air pressure variations are fully compensated by static sea height changes, or using ocean general circulation models. We validate two runs of the HUGO-m barotropic ocean model by comparing predicted sea surface height variations with hundred tide-gauge measurements along the European coasts. We then show that global surface pressure field, as well as a barotropic high-resolution ocean model forced by air pressure and winds allow in most cases a significant reduction of the variance of gravity residuals and, to a smaller extends tilt residuals.We finally show that precise gravity measurements with superconducting gravimeters allow the observation of large storm surges, occurring in the North Sea, even for inland stations. However, we also confirm that the continental hydrology contribution cannot be neglected. Thanks to their specific sensitivity feature, only tiltmeters closest to the coast can clearly detect the loading due to these storm surges.     

精河6.6级地震前精河水管仪、伸缩仪异常特征分析

2017年8月9日精河发生6.6级地震。本文首先利用相关、回归分析研究了气温对水管仪、伸缩仪的影响,并进行了干扰的定量剔除;其次,通过回归残差分析研究地震前水管仪、伸缩仪的异常特征,认为精河6.6级地震前存在3项异常:①2016年5月22日起,水管仪NS分量S倾幅度0.53";②2017年3月24日起,EW分量E倾幅度0.28";③2016年9月10日起,伸缩仪NS分量出现压缩现象。     

A search for atmospheric effects on gravity at different time and space scales

We investigate in this paper various approaches to correct gravity changes for the effect of atmospheric pressure changes. Two specific locations are considered: Strasbourg (France) as mid-latitude station, where regular pressure fronts occur and Djougou (Benin) as equatorial station with large thermally driven S₁ and S₂ waves of planetary extent. We first review the classical approaches based on a constant or frequency-dependent admittance using only local pressure and gravity data. We consider then a model of atmospheric loading and show the barometric admittance in terms of elastic, Newtonian and total load, as a function of the distance from the station. We consider both a 2D pressure model (surface loading) and a 2.5D model, where the density decreases with height (standard atmosphere). Assuming horizontal advection in the atmospheric dynamics, we convert this spatially dependent admittance into a frequency-dependent admittance. Using global pressure data from European Centre for Medium-Range Weather Forecasts (ECMWF) at about 12 km spatial resolution and 3 h sampling, we compute the model-predicted pressure admittance for Djougou and Strasbourg and we simulate the frequency dependence inferred from gravity and pressure observations below 4 cycle per day. A long gravity and pressure data set (1996–2013) from Strasbourg is used to investigate the low frequency part of the pressure admittance while a common 2.5 year data set (August 2010–February 2013) for Strasbourg and Djougou is then analyzed to investigate the high frequency part of the admittance. In both cases, our results are in close agreement with the predictions inferred from an atmospheric 2.5D loading model with a distance–time relationship due to horizontal advection. The frequency dependence of the barometric admittance is explained by the competing contributions of Newtonian attraction and elastic surface deformation according to the distance from the gravimeter. In the far field (low frequencies), the magnitude of the admittance decreases with frequency because of the combined elasticity effect and Newtonian attraction (when the atmosphere is below the horizon) while, on the contrary, in the near field (high frequencies), elasticity becomes negligible and the pressure admittance mainly decreases with increasing frequency because of the decreasing attraction effect of the atmospheric masses inside the cylindrical pressure cell centered on the sensor location of decreasing radius. In the last part, we show that there is variability in time in the pressure admittance for both stations.     

地球自转的天文地球动力学效应研究(Ⅰ)--引力位和重力的变化

详细讨论了地球自转的天文地球动力学效应，地极移动和日长变化导致地球引力位系数产生时变特性并引起重力的摄动，根据理论力学的基本概念，导出了由于地球自转变化引起的地球引力位系数变化、重力摄动、垂线偏差和地球形变的表达式，并定量地研究了极移和日长变化对测站重力观测值和地球形变的影响，建议在高精度的空间大地测量中要顾及到地球自转变化引起的一系列效应。     

Clear evidence for the sign-reversal of the pressure admittance to gravity near 3 mHz

The influence of the atmosphere on gravity measurements by Newtonian attraction and vertical displacements due to surface loading is well known and studied especially at frequencies below 1 mHz. Less work has been done at the higher seismic frequencies where inertial effects come into play. The sensor mass of a vertical accelerometer responds to several forces caused by the atmosphere, on global, regional and local scales. For the seismic frequencies the atmosphere above the observation site has the largest influence. A simple “Gedankenexperiment” demonstrates that the gravitational effect on one hand and the free air and inertial effects due to deformation on the other have opposite sign. Since the inertial effect is strongly frequency-dependent there should be a crossover-frequency where the pressure admittance to gravity changes sign. Simple analytical models clearly show this property near frequencies of a few mHz and are used here to amplify the variance reduction of the gravity residuals. The crossover-frequency depends on the properties of the models of the atmospheric phenomena and the elasticity of the Earth’s crust. Therefore in reality it must be expected to vary in time and space.We show one extremely clear example of this sign-reversal detected in the gravity data from the superconducting gravimeter GWR-C025 in Vienna among other examples which demonstrate the reality of this effect.     

甘肃祁连山主动源不同枪组激发效果对比

大容量气枪震源激发探测因具绿色环保、传播距离远和高度重复性等特点成为地下介质探测的新方法,利用其进行地下介质物性变化高精度探测,则对不同激发条件下气枪传播信号的一致性提出要求。本文选取了祁连山主动源不同枪组（4支和3支气枪）在激发环境一致情况下分别激发时各远台波形记录,通过信号叠加、互相关及频谱分析等方法,对比分析了各台站记录的不同枪组激发信号,及其与震中距的关系,并讨论了不同枪组激发对地下介质走时变化观测的影响。结果表明：（1）不同枪组的激发信号波形相关度随震中距增大而降低,信号震相随激发能量有微弱的变化,震中距超过50 km后其变化几乎分辨不出;（2）远台接收信号的能量大小随主动源激发能量变化,各台站能量变化主要体现在优势频率范围内;（3）接收台站台基越坚固不同激发能量均方根振幅比越小;（4）叠加一定数量的波形可以有效提高远台记录信号的信噪比,增加气枪激发能量则相对更容易提高信号信噪比。