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

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

大气变化对位移、重力和倾斜观测影响的理论计算

利用弹性地球负荷理论和流体静力学平衡假设条件下的干空气及其垂直温度梯度分布,理论计算了大气压力作用下地表位移、重力和倾斜格林函数．详细讨论了不同地球模型、台站高程、周边地形、地表大气温度和水汽含量等因素对计算结果的影响．数值结果说明大气效应弹性项影响主要来自于远离台站的区域;引力项对重力的影响主要贡献来自于台站近区,对倾斜的影响来自于远离台站的区域．各种因素对大气重力和倾斜格林函数的影响分别在01°和1°角距范围内,但其综合影响较为复杂．     

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

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

水质量分布变化对重力观测的影响

利用全球7000多个陆地气象台站观测资料、TOPEX/Poseidon卫星测高资料、多层海水温盐度数据和负荷格林函数,采用数值积分法,计算了陆地水储量及海水质量分布变化引起我国20个测站的重力变化.结果表明,我国一些测站观测的重力受水质量迁移影响超过3μGal,可见在应用重力仪观测资料研究地球动力学时,有必要作水负荷改正.另外,文中还对武汉测站的超导重力仪观测的残差和水负荷的计算结果作了比较,并讨论了不同的地球模型及气象数据对计算测站重力的影响.     

不同地球模型对中国陆区大气负荷位移建模的影响

目前进行质量负荷位移建模时通常使用的是全球或区域平均的一维地球模型，这类地球模型无法反映地壳结构的局部差异，如地壳厚度、地震波波速等.本文利用地表气压再分析数据产品NCEP/NCAR Reanalysis 1(NCEP R-1)和不同地球模型(Preliminary reference Earth model(PREM)、AK135、STW105和CRUST1.0)计算的负荷格林函数计算中国陆区的大气负荷垂向形变，结合GPS实测位移数据，评估基于不同地球模型计算得到的地表形变的差异.通过比较PREM与不同地球模型的建模结果发现，一维地球模型之间的建模差异较小，主要集中在气压变化较大的东部地区和西北地区，气压变化较小的青藏高原地区的建模差异则更小.相反，PREM与CRUST1.0模型在青藏高原地区的建模差异最大，在东部地区的建模差异最小.PREM与STW105、AK135和CRUST1.0在中国陆区的最大位移建模差异分别为0.05、0.34和0.44 mm.不同地球模型的建模差异是由负荷格林函数不同产生的，而负荷格林函数在近场受上地壳厚度的影响较大.对于局部地壳结构差异明显的青藏高原地区...     

大气潮对重力固体潮影响的研究

通过气压变化对重力观测的影响和大气潮对重力固体潮的负荷效应二方面的研究,得出将ＧＳ型重力仪置于密封装置,能使气压的影响降低１个数量级。     

浅谈银川台洞温变化对重力观测资料的影响

通过对银川台重力观测资料进行差分变化对比,与理论固体潮及室温变化进行对比分析,阐述了2007年台站重力观测山洞改造期间引起重力变化的主要因素。研究结果表明,山洞改造对重力观测的直接影响不显著,引起资料内在精度下降的主要原因是改造时保温措施不当,引起洞温迅速变化,继而影响重力观测。     

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

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

温度变化对我国GPS台站垂直位移的影响

地表温度变化可以引起GPS台站上安装GPS天线的地表水泥墩内部温度变化,还可以通过热传导的方式引起GPS台站基岩温度变化,从而引起GPS台站垂直位移变化.在中国区域,由GPS台站基岩温度变化引起GPS台站垂直位移变化的周年振幅最大可以达到1 mm;在长江以北地区,此周年振幅一般大于0.5 mm.在我国地壳运动观测网络中的23个GPS基准站中,温度变化对GPS台站垂直位移总影响的周年振幅最大值为2.8 mm,其中13个GPS基准站垂直位移的周年振幅变化大于1 mm.因此,温度变化是引起GPS台站垂直位移周年变化的一个不可忽视因素.     

不同产状断层错动的地表重力变化和形变

根据点源位错引起的地表重力变化和形变的理论公式,应用数值算法求得半无限空间任意产状断层位错在地表的重力变化和形变．计算结果表明,产状和位错方式是决定断层的地表重力和形变效应的基本因素．铲式断层活动的地表重力和形变效应与矩形断层相比有显著差异．断层位错引起的重力变化、垂直形变、视垂直形变的图象特征相似．根据重力变化和形变的特殊图象特征,有可能从观测资料中提取断层活动引起的重力变化和形变信息,判识和估计发震断层的几何学和运动学特征．      

Green's functions and response functions in geochemical modelling

The use of Green's functions and response functions in modelling atmospheric and oceanic tracers is discussed. It is shown that low-resolution Green's function representations can describe tracer distributions more accurately than low-resolution models. The Green's function formalism is also used to discuss the problem of determining surface CO₂ sources from observations of surface CO₂ concentrations.     

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.     

Gravity reduction with three-dimensional atmospheric pressure data for precise ground gravity measurements

The redistribution of air masses induces gravity variations (atmospheric pressure effect) up to about 20 μgal. These variations are disturbing signals in gravity records and they must be removed very carefully for detecting weak gravity signals. In the past, different methods have been developed for modelling of the atmospheric pressure effect. These methods use local or two-dimensional (2D) surface atmospheric pressure data and a standard height-dependent air density distribution. The atmospheric pressure effect is consisting of the elastic deformation and attraction term. The deformation term can be well modelled with 2D surface atmospheric pressure data, for instance with the Green's function method. For modelling of the attraction term, three-dimensional (3D) data are required. Results with 2D data are insufficient.From European Centre for Medium-Range Weather Forecasts (ECMWF) 3D atmospheric pressure data are now available. The ECMWF data used here are characterised by a spacing of Δ and Δλ = 0.5°, 60 pressure levels up to a height of 60 km and an interval of 6 h. These data are used for modelling of the atmospheric attraction term. Two attraction models have been developed based on the point mass attraction of air segments and the gravity potential of the air masses. The modelling shows a surface pressure-independent part of gravity variations induced by mass redistribution of the atmosphere in the order of some μgal. This part can only be determined by using 3D atmospheric pressure data. It has been calculated for the Vienna Superconducting Gravimeter site.From this follows that the gravity reduction can be improved by applying the 3D atmospheric attraction model for analysing long-periodic tidal waves including the polar tide. The same improvement is expected for reduction of long-term absolute gravity measurements or comparison of gravity measurements at different seasonal times. By using 3D atmospheric pressure data, the gravity correction can be improved up to some μgal.     

A spectrum of the acoustic-gravity waves associated with Auroral electro-jet

This model concerns the analysis of Aurora initiated travelling pressure waves in an isothermal atmosphere. Electro-dynamic Lorentz Force associated with auroral electric current density during the periods of geomagnetic activity is invoked as a possible exciting source.The dispersion phenomena in auroral induced acoustic gravity modes in the earth's atmosphere are examined and various cut-off frequencies analysed.Finally, an attempt is made to obtain a far field representation of the forced oscillations by means of Green's Function technique. Therefrom, the spectral amplitudes of the ground-level oscillations are computed. Incorporating various scaling factors, it is deduced that these amplitude components are in reasonable agreement with the results of recent measurements.     

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.     

基于OH全天空气辉成像仪观测的中国低纬地区的重力波传播统计特征

利用位于海南富克(19.5°N,109.1°E)和广西桂平(23.4°N,110.1°E)两个台站两年多的OH全天空气辉成像仪观测数据,对中国低纬地区的重力波传播统计特征进行了研究.从富克和桂平的气辉成像观测中, 分别提取了65和86个重力波事件.研究结果表明,观测水平波长,观测周期和水平相速度分别集中分布在10~35 km, 4~14 min和20~90 m·s^-1范围.重力波传播方向,在夏季表现出很强的东北方向传播.然而,在冬季主要沿东南和西南方向传播. 同时,结合流星雷达风场观测和TIMED/SABER卫星的温度数据,也发现在中层-低热层中传播的大多数重力波表现为耗散传播.且低层-中层大气中背景风场的滤波作用和多普勒频移可能对纬向方向传播的重力波产生的各向异性起到重要的调制作用.然而,经向方向传播的重力波产生的各向异性可能同时被低层大气中波源的非均匀分布以及潮汐变化所影响.     

Annual effect of cosmic rays on ionization in the high-latitude troposphere

The seasonal effect of the daily variations in the cosmic ray intensity on the conductivity of the Earth-high-conductivity layer column has been analyzed based on the observations of the cosmic ray intensity, atmospheric current, and electric field vertical component, performed from summer 2006 to spring 2007 at Apatity station. The method for correcting the measurements of the atmospheric current and electric field vertical component under complex tropospheric conditions by numerically simulating the spatial structure of the current and field lines in the observation region has been proposed. It has been indicated that cosmic rays are the main source of ions in the winter polar lower atmosphere and are responsible for the type of daily variations in the conductivity, whereas the daily variations in the atmospheric current more depends on the conductivity rather than on the vertical electric field.