基于EEMD分解的气压改正方法研究

利用超导重力仪观测数据精确测定低于1 mHz的地球自由振荡简正模式的分裂频率,是在不与任何弹性系数发生联系的情况下改善一维密度模型的有效方法.但在该频段台站局部气压变化对重力观测数据的影响成为主要干扰来源,且具有频率依赖特性,因此精细地开展气压改正成为利用超导重力数据检测低频自由振荡信号的必要手段.本文基于EEMD方法,提出了一种具有频率依赖特性的气压改正方法.该方法将重力观测和气压变化分解成处于不同频段的本征模态函数,并在相应频段上分别进行重力-气压变化的回归分析,计算得到具有频率依赖特性的气压导纳值,精细地消除气压变化对重力观测的影响,并以此对微弱低频地球自由振荡信号开展高分辨率分析.基于本文提出的气压改正方法,利用大地震后的超导重力数据检测了频率小于1.5 mHz的低频地球自由振荡及其频谱分裂现象.研究结果表明：利用该方法进行气压改正后检测得到的各简正模具有更高的信噪比,估计的本征频率误差水平明显降低,获得的基频球型振荡₀S₂和₀S₃以及一阶球型振荡₁S₂的分裂谱峰的估计精度更高,同时还检测到了部分环型振荡在重力观测中的耦合现象.对低频地球振荡的高分辨率检测结果验证了基于EEMD分解提出的气压改正方法的有效性,同时再次证明了超导重力仪观测数据在低频地球自由振荡检测中的优势.

     

琼中台重力非潮汐变化气压与海潮负荷改正

本文对琼中台连续重力观测数据进行收集整理并处理,基于处理后的数据,进行了潮汐分析和非潮汐分析。潮汐分析采用VAV调和分析方法;非潮汐分析则分别进行了零漂改正、固体潮改正、气压改正和海潮改正。其中,零漂改正采用一般多项式拟合零漂的方法;气压改正采用VAV软件;海潮改正运用SPOTL程序,以NAO.99b潮汐模型计算了琼中台海潮负荷值。最终获得了改正后的琼中台重力非潮汐变化,结果表明琼中台的重力气压导纳值为-0.34×10^-8m/s²/mbar,气压改正幅度约为10×10^-8m/s²,海潮改正幅度约为5×10^-8m/s²。改正后,琼中台重力非潮汐变化数据,比仅进行零漂固体潮改正后的重力非潮汐变化数据中的潮汐信号更加微弱,说明进行海潮改正后的效果是明显的,该方法可进一步去除其中的潮汐信号。     

武汉超导重力仪观测最新结果和海潮模型研究

利用武汉台站GWR_C032超导重力仪观测资料,在对原始数据进行有效预处理的基础上作调和分析,获得反映地球内部介质特征的重力潮汐参数.基于卫星测高技术和有限元方法同时考虑验潮站数据作约束条件获得的多个全球海潮模型,利用负荷理论和数值褶积积分技术计算了重力负荷,对周日和半日频段内的重力潮汐参数实施负荷改正,提出了“负荷改正有效性”概念,研究了全球海潮模型适应性.数值结果说明,海潮改正的有效性高达91%（O1,NAO99）和92%（M2,ORI96）.基于11个海潮模型对主波（O1,K1,M2和S2）的负荷改正说明平均有效性为（86%,70%,73%和84%）,振幅因子与理论模型间的差异分别从（212%,155%,116%和080%）降到（031%,039%,034%和008%）,同时还说明利用NAO99和ORI96全球海潮模型能获得比其他模型更佳的负荷改正效果.文章还利用国际地球动力学计划网络其他7个台站的超导重力仪观测研究了全球海潮模型的适定性问题,结果说明不同模型中不同潮波具有明显的区域特点,早期构制的SCW80全球海潮模型仍可作为大地测量研究中的重要参考模型.     

GS型重力仪面板常数测定的研究

系统分析了GS - 1 5 - 2 1 6重力仪多次在重力基线场标定的结果 ,并与厂家给出的面板常数进行对比 ,认为GS型重力仪的面板常数存在长期变化 ,应定期重新标定 ;检修或工作状态的变化也会引起面板常数的改变 ;面板常数是仪器测程的函数 ,必须分段标定以便拟合一条常数曲线 .     

FG5绝对重力仪

介绍了ＦＧ５型绝对重力仪的发展历史和一些技术特点，并报道了ＦＧ５／１１２型绝对重力仪在我国的观测结果。     

绝对重力仪测量结果的振动划线处理方法初步研究

针对自主研发的高精度绝对重力仪样机测量结果中引入的振动误差, 提出了一种全新的振动抑制算法．本算法利用一簇斜率相同的平行直线与记录到的参考棱镜振动曲线相交, 通过提取这些交点的时间信息, 确定落体下落轨迹上进行重力加速度计算的时间和位移坐标． 然后, 利用最小二乘算法获取测点的绝对重力加速度值．对不同重力加速度值的适用性和不同类型振动干扰的抑制效果模拟计算发现, 对于优势频率大于0.2 Hz, 峰峰值小于500 counts的干扰, 算法可以有效提高测量结果的系统偏差和测量精度．最后, 实际实验样机的测量结果表明, 本算法使最终测量精度提高了一倍, 为绝对重力仪的小型化和流动测量提供了软件技术支持。      

应变固体潮资料与气压的相关分析及其改正的一种方法

通过对应变资料和气压资料的频谱分析，将固体潮频段的信息分别提取出来，再进行相关分析，二者显示出了很强的线性相关性。据此，对应变观测资料进行了气压改正，使应变固体潮调和分析的各项精度指标有了明显的提高，指出该方法的实用效果，同时，对气压改正中存在的问题也进行了讨论     

关于绝对重力仪的发展和我们的思考

介绍了近年来出现的新型绝对重力仪，认为该仪器是控制地球重力场信息的重要工具，特别对FG5型的如何改进及我国绝对重力仪的如何发展，提出了建议。     

Signal path-curving caused by atmospheric refraction

When a signal propagates in the atmosphere, a small variation will occur in its propagation direction because of refraction. The signal path will be a curve other than a straight-line. The characteristics of signal direction variation in propagation in the atmosphere are detailed, new and more precise path-curving correction formulas are derived on the basis of a new atmospheric layer-built model.     

Some effects of intra-red radiation on the dynamics of the atmospheric boundary layer

An approximate infra-red radiation scheme, employing essentially the cooling to space approximation, is included in a one-dimensional model of the atmospheric boundary layer. The approximate scheme is found not to produce significant errors in the behaviour of the dynamical model when integrated over a few hours. Radiative cooling is shown to be important in the development of a clear well-mixed layer which is capped by an essentially dry region; in particular, convetive instability is induced and this enhances the rate of entrainment of dry air. The development of fog is found to require sound models of both radiative transfer and turbulent diffusion.     

基于GPS和大气传输模型的InSAR大气改正方法研究

利用GPS数据改正InSAR大气影响中,GPS站网的低空间密度是限制改正精度的一个主要因素.文中引入大气传输模型（ATM）来考虑大气状态在时间上的演化,并兼顾了风向的估计,把GPS-ZWD（GPS湿延迟）时间序列观测值转换为空间上分布比较稠密的GPS-ZWD网络,然后利用该稠密网络在空间上内插产生水汽延迟图,来更好地模拟InSAR影像获取时刻的水汽场,提高大气改正精度.实验结果表明,在研究地区有限的数据条件下（即只有6个GPS点的情况下）,GPS+ATM算法在重现大气信号能力方面比单纯的使用影像获取时刻的GPS-ZWD数据要强,并且其对长波的大气误差去除更明显;用来做实验的三幅差分干涉图中,GPS+ATM算法对两幅以长波信号为主的干涉图中的大气影响分别降低了21.7%和22.6%,比仅使用SAR过境时刻GPS-ZWD数据时的结果分别改进了4.5%和8.7%.而对以短波信号为主的干涉图没有明显提高.     

几何校正对暗像元算法及离水辐亮度反演精度的影响:以太湖为例

几何校正将引起几何校正前后像元的像素值发生变化,进而导致暗像元的辐亮度在几何校正前后存在偏差.本文以太湖为研究水区,以2003年10月28日的Landsat/TM影像、2008年7月24日的CBERS影像和2009年4月25 日的HJ-1B影像作为数据基础,研究与探讨了几何校正对暗像元大气校正算法以及离水辐亮度反演精度...     

The effects of small water bodies on the atmospheric heat and water budgets over the MacKenzie River Basin

The effects of small water bodies or lakes on the surface sensible and latent heat fluxes and the transport of heat and water vapour in the atmospheric boundary layer (ABL) over the Mackenzie River Basin (MRB) are studied from two cases, which occurred on 2 and 8 June 1999 during the warm season. The synoptic condition for the cases is representative of about 33% of the synoptic situation over the MRB. The two events are simulated using the Canadian mesoscale compressible community (MC2) model. A one‐way nesting grid approach is employed with the highest resolution of 100 m over a domain of 100 km². Experiments were conducted with (LAKE) and without (NOLAKE) the presence of small water bodies, whose size distribution is obtained through an inversion algorithm using information of their linear dimension determined from aircraft measurement of surface temperature during MAGS (the Mackenzie GEWEX (Global Energy and Water Cycle Experiment) Study) in 1999. The water bodies are assumed to be distributed randomly in space with a fractional area coverage of 10% over the MRB. The results show that, in the presence of lakes, the domain‐averaged surface sensible heat flux on 2 June 1999 (8 June 1999) decreases by 9·3% (6·6%). The surface latent heat flux is enhanced by 18·2% (81·5%). Low‐level temperature advection and the lake surface temperature affect the air–land/lake temperature contrast, which in turn controls the sensible heat flux. In the absence of lakes the surface wind speed impacts the latent heat flux, but in the presence of lakes the moisture availability and the atmospheric surface layer stability control the latent heat flux. The enhancement is smaller on 2 June 1999 as a result of a stable surface layer caused by the presence of colder lake temperatures. The domain‐averaged apparent heat source and moisture sink due to turbulent transports were also computed. The results show that, when lakes are present, heating and drying occur in the lowest 100 m from the surface. Above 100 m and within the ABL, there was apparent cooling. However, the apparent moistening profiles reveal that lakes tend to moisten the ABL through transfer of moisture from the lowest 50–100 m layer. Copyright © 2004 John Wiley & Sons, Ltd.     

考虑雷暴活动的全球近地稳态大气电场的精确数学解

全球大气电模式是指大气和地表构成的闭合回路,其中考虑雷暴为主要充电电源,并将雷暴作为一个垂直电偶极子来考虑.本文从电流连续方程出发,采用球谐函数展开法,边界条件采用下边界为地表大气电势,上边界为电离层底部电势,非齐次电流源函数项采用垂直偶极子假设.通过所做的非线性函数变换,得出了考虑雷暴活动的全球电模式的近地稳态解析解,并证明其能展开为有限多项式级数之和.所得结论与过去的理论分析结果和实验观测是一致的.该文的工作是对以前相关工作仅能得到近似解的一种改进,同时本文所得公式可用于今后对大气电参量进行相关的数值计算及模拟之中.     

Adaptability of the ocean and earth tidal models based on global observations of the superconducting gravimeters

As the gravity field is the most primary and direct physical quantity reflecting the density variation of the Earth’s interior and its geodynamic properties under various environmental changes and as the knowledge of the fine structure of the Earth’s interior and its geo-dynamics has a strong impact on space research, grav-ity observations become more and more important in Earth sciences. Therefore based on a global network of superconducting gravimeters (SGs), the Global Geo-dynamic Proje…     

Wavelet filter analysis of atmospheric pressure effects in the long-period seismic mode band

The importance of the reduction of atmospheric pressure effects becomes very clear when investigating seismic normal-mode spectra below 1.5 mHz. The usual simple correction method consists in subtracting a term converted from local atmospheric pressure (pressure multiplied by a frequency-independent admittance) from the gravity record in time domain. Thus, estimating an efficient admittance is the key for an improved correction. Band-pass filters derived from dyadic orthogonal wavelet transform, having narrow pass-bands with good frequency response but without Gibbs phenomenon and causing no phase lag, are very helpful to estimate an efficient admittance, which is both time and frequency-dependent. Processing of high quality superconducting gravimeter (SG) records for the great Sumatra earthquake (Mw = 9.3, Dec 26, 2004) with wavelet filters reveal the three very well resolved splitting singlets of overtone ₂S₁ with a single gravity record after correction with time-dependent and frequency-dependent admittances. We also observe all coupled toroidal modes below 1.5 mHz, except ₀T₅, ₀T₇ and ₁T₁, with good signal-to-noise ratio (SNR); moreover, toroidal modes ₁T₂ and ₁T₃ are for the first time unambiguously revealed in vertical components.     

一个大气边界层非线性三维数值模式及其在太湖流域的应用

本模式用以K半经验理论为基础的闭合方程组求解,引进了地形作用、动量、热量和水汽的侧向扩散和垂直湍流交换以及凝结加热等物理过程,其中湍流交换系数K是地表粗糙度、几何高度、大气稳定度及风速切变等因子的函数。本模式还考虑了山地、平原和水面在地形高度、地表粗糙度、辐射、蒸发及热交换等方面的差异,在下垫面上建立热平衡方程与边界层内控制方程耦合。将模式应用于太湖流域,计算得到的该流域边界层内温、压、湿、风的分布特征与实际情况相似,模式具有一定的实用性。