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

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

基于集合经验模态分解法的重力极潮提取与研究

重力极潮,即极移引起的重力变化,是地球由于惯性离心力变化导致地球形变的综合反映,其观测和研究有助于了解地球在长周期频段的响应,约束地球形变、地球内部构造和物理参数.区别于以往研究,本文采用了集合经验模态分解方法（Ensemble Empirical Mode Decomposition,EEMD）从全球5个超导重力仪台站连续重力观测资料中提取重力极潮、消除仪器漂移及极潮频段以外的大部分噪声信号,该方法所用的基函数是基于自身信号获得的;在此基础上,利用最小二乘匹配法分离极潮中的周年项和Chandler项,估算Chandler周期的极潮潮汐因子.结果表明,基于EEMD得到的极潮潮汐因子与前人基于其他消除仪器漂移方法（数学模型或小波分析）得到的结果相符合,精度相当,但由于这种方法是自适应的,因此本文结果更能反映实际物理过程.     

Temporal variability of dimethylsulfide and dimethylsulfoniopropionate in the Sargasso Sea

Vertical profiles of dimethylsulfide (DMS) and β-dimethylsulfoniopropionate, particulate (pDMSP) and dissolved (dDMSP), were measured biweekly in the upper 140 m of the Sargasso Sea (32°10′N, 64°30′W) during 1992 and 1993. DMS and pDMSP showed strong, but different, seasonal patterns; no distinct intra-annual pattern was observed for dDMSP. During winter, concentrations of DMS were generally less than 1 nmol l⁻¹ at all depths, dDMSP was less than 3 nmol l⁻¹ and pDMSP was less than 8 nmol l⁻¹. In spring, concentrations of both dDMSP and pDMSP rose, on a few occasions up to 20 nmol l⁻¹ in the dissolved pool and up to 27 nmol l⁻¹ in the particulate pool. These increases, due to blooms of DMSP-containing phytoplankton, resulted in only minor increases in DMS concentrations (up to 4 nmol l⁻¹). Throughout the summer, the concentrations of DMS continued to increase, reaching a maximum in August of 12 nmol l⁻¹ (at 30 m depth). There was no concomitant summer increase in dDMSP or pDMSP. The differences among the seasonal patterns of DMS, dDMSP, and pDMSP suggest that the physical and biological processes involved in the cycling of DMS change with the seasons. There is a correlation between the concentration of DMS and temperature in this data set, as required by some of the climate feedback models that have been suggested for DMS. A full understanding of the underlying processes controlling DMS is required to determine if the temperature-DMS pattern is of significance in the context of global climate change.     

地球液核的动力学效应研究进展

简要介绍了液核动力学研究的方法及液核动力学效应检测的进展和研究结果,着重介绍了超导重力仪在液核动力学研究中的作用。基于初始参考地球模型（PREM）,采用球对称、非自转、完全弹性和各向同性地球的弹性引力形变理论研究了液核动力学扰动导致的地球固体部分的形变和重力位扰动。根据重力潮汐观测中的近周日共振特征,利用国际超导重力仪观测资料研究了地球的自由核章动（FCN）,精密确定了有关的共振参数,其中FCN的本征周期为429.0（424.3,4 433.7）恒星日,品质因子为 9543（6405,18714）,复共振强度为（－6.10±0.20,－0.01±0.20）×10^－4°/h。最近,我们还估计了全球地球动力学观测网中全球分布的14台超导重力仪21个长期、连续重力观测序列的"积谱密度"以检测固体内核的平动振荡运动。     

Analysis of the tidal gravity observations in Three Gorges reservoir area

The comparative observations of the tidal gravity in Three Gorges reservoir area, recorded with the LCR-ET20 gravimeter from Institute of Geodesy and Geophysics, Chinese Academy of Sciences and the DZW-9 gravimeter from Institute of Seismology, China Earthquake Administration, are introduced. High-accuracy tidal gravity pa-rameters in this area are obtained, and the scale factor of the DZW-9 gravimeter is also estimated at a value of ?(756.06±0.05)×10?8 m?s?2V?1. The observed residuals of the corresponding instruments are investigated. Nu-merical results indicate that the long-term drift of the DZW-9 gravimeter appears a linear characteristic, and the observed accuracy is of the same order comparing with that of the LCR-ET20 gravimeter. The results given in the paper can provide with an effective reference model of the tidal gravity correction to the ground based and space geodesy.     

特大型精密工程SSC测量投影方式的选择

目前世界上最大的工程超导超级对撞机(SSC)正在建造之中,作者参与这项大型精密工程测量的方案设计和部分咨询工作。根据工程的特点,设计测量投影方式是其中一项工作。本文讨论所设计的双重正形投影(stereographic double conformal projection)以及各种观测量投影归化中的问题。     

超导重力技术在探讨核幔边界黏性特征中的初步应用

旋转椭球型地球的固体地幔与液态地核间相互作用而产生的逆向本征模通常称之为地球自由核章动,自由核章动的品质因子（Q值）能有效反映核幔边界层能量耗散特征,与核幔边界的黏滞度密切相关.本文首次利用全球地球动力学计划网络23个台站27组高密度采样的高精度超导重力仪器观测数据,采用迭积技术,确定了自由核章动参数Q值,进而计算了核幔边界的黏滞系数.数值结果说明获得的核幔边界动力学黏滞系数达到10³ Pa·s量级,与加拿大科学家Smylie等利用VLBI观测资料获得的最新结果一致,这说明重力技术是有效应用于研究地球深内部结构的重要手段之一.     

Hydro-gravimetry in West-Africa: First results from the Djougou (Benin) superconducting gravimeter

The increasing number of hydro-gravimetry studies proves the rising interest of the hydrology community toward this monitoring method. The accuracy of superconducting gravimeters (SG) potentially allows the retrieval of small water storage changes (WSC) down to a few millimeters of equivalent water thickness. However, the importance of corrections applied to SG data to achieve such a precision in gravity residuals should be recalled. The Djougou permanent gravity station presented in this paper and located in northern Benin, West-Africa, provides a good opportunity to review these considerations. This station is equipped since July 2010 with the superconducting gravimeter SG-060 aimed at deriving WSC at different time-scales, daily to inter-annual. In this area, WSC are (1) part of the control system for evapotranspiration (ET) process, a key variable of the West-African monsoon cycle and (2) the state variable for resource management, a critical issue in storage-poor hard rock basement contexts such as in northern Benin. The potential for deriving WSC from time-lapse gravity data partly depends on environmental features such as topography and the instrument shelter. Therefore, this issue is addressed first, with the background idea that such sensitivity analysis should be undertaken before setting up any new instrument. In Djougou, local topography is quite flat leading to a theoretical straightforward relationship between gravity changes and WSC, close to the standard Bouguer value. However, the shelter plays a significant masking role, which is the principal limitation to the retrieval of fast hydrological processes such as ET following a rain event. Several issues concerning classical gravity corrections are also addressed in the paper. These include gap-filling procedures during rain-events and drift estimates for short time series. Special attention is provided to atmospheric corrections, and different approaches are tested: a simple scalar admittance, a filtered scalar admittance, a frequency-dependent admittance and direct atmospheric loading calculations. It is shown that the physically based approach of direct loading calculations performs better in both residual minimization and ET retrieval. Moreover, non-local hydrological effects are investigated and account for about 20% of the gravity residuals. Finally, gravity residuals are briefly analyzed at two distinct time scales: rapid (up to a few days) and seasonal. At the rapid time-scale, it is shown that ET retrieval is hardly achievable given shelter size and state-of-the-art atmospheric corrections. Still, mean values retrieved from this study are in accordance with known values of potential ET and lateral flow. Direct comparison of gravity changes with hydrological data (neutron probe monitoring and water table levels) show some discrepancies, particularly for the hydrological year of 2011, for which all hydrological data show a deficit, but SG and FG5 data do not. This preliminary analysis both provides a basis and call for further hydro-gravity modeling, to comprehensively investigate the water-cycle at the Djougou station.     

On the comparison of tidal gravity parameters with tidal models in central Europe

Two accurately calibrated superconducting gravimeters (SGs) provide high quality tidal gravity records in three central European stations: C025 in Vienna and at Conrad observatory (A) and OSG050 in Pecný (CZ). To correct the tidal gravity factors from ocean loading effects we compared the load vectors from different ocean tides models (OTMs) computed with different software: OLFG/OLMP by the Free Ocean Tides Loading Provider (FLP), ICET and NLOADF. Even with the recent OTMs the mass conservation is critical but the methods used to correct the mass imbalance agree within 0.1 nm/s². Although the different software agrees, FLP probably provides more accurate computations as this software has been optimised. For our final computation we used the mean load vector computed by FLP for 8 OTMs (CSR4, NAO99, GOT00, TPX07, FES04, DTU10, EOT11a and HAMTIDE). The corrected tidal factors of the 3 stations agree better than 0.04% in amplitude and 0.02° in phase. Considering the weighted mean of the three stations we get for O1 δ_c = 1.1535 ± 0.0001, for K1 δ_c = 1.1352 ± 0.0003 and for M2 δ_c = 1.1621 ± 0.0003. These values confirm previous ones obtained with 16 European stations. The theoretical body tides model DDW99/NH provides the best agreement for M2 (1.1620) and MATH01/NH for O1 (1.1540) and K1 (1.1350). The largest discrepancy is for O1 (0.05%). The corrected phase α_c does not differ significantly from zero except for K1 and S2. The calibrations of the two SG's are consistent within 0.025% and agree with Strasbourg results within 0.05%.