球形电容传感器研究进展

1 研究背景 利用传统平板式电容传感器对曲面结构进行检测时,难以设置合适的位置,而且传感器探头与被测物体的平行度对精度有较大影响.相比之下,球形电容传感器以其非平面检测方式,在各个方向具有相同的传感特性,对特殊曲面结构进行检测时具有较大优势.球形电容传感器发展至今大致为3种:球探头形、差动型以及柔性球形电容传感器.根据...     

宽频带超导重力仪的研制进展

本文报道一种构建宽频带高分辨率超导重力仪（SG）的方法.与已商业化的GWR超导重力仪类似,新型超导重力仪同样使用磁悬浮超导检验质量来构建弹簧振子结构,不同的是弹簧振子的固有频率更高,可以有效覆盖地脉动频带.同时,利用基于超导量子干涉器件（SQUID）的位移传感技术,大幅提高了位移传感灵敏度,以保证仪器在测量带宽扩展之后仍具有与GWR仪器相当的测量分辨率,从而实现对时变重力信号与地球背景噪声信号的同时高分辨观测.宽频带超导重力仪核心部件垂向超导加速度计（VSA）的设计、组装和初步测试已经完成.结果表明,加速度计的噪声本底为8 μGal/Hz^1/2@0.1 Hz,实现了对地脉动信号的高分辨率测量.而为了将仪器噪声本底降低到GWR仪器水平,加速度计的固有频率需要进一步减小,该部分参数优化工作正在进行.本文将对垂向超导加速度计的工作原理、结构和测试结果进行详细介绍.

     

武汉台重力潮汐长期观测结果

采用武汉台超导重力仪(SG C032)14年多的长期连续观测资料,研究了固体地球对二阶和三阶引潮力的响应特征,精密测定了重力潮汐参数,系统研究了最新的固体潮模型和海潮模型在中国大陆的有效性.采用最新的8个全球海潮模型计算了海潮负荷效应,从武汉台SG C032的观测中成功分离出63个2阶潮汐波群和15个3阶潮汐波群信号,3阶潮波涵盖了周日、半日和1/3日三个频段.重力潮汐观测的精度非常高,标准偏差达到1.116 nm·s-2,系统反映了非流体静力平衡、非弹性地球对2阶和3阶引潮力的响应特征.结果表明,现有的武汉国际重力潮汐基准在半日频段非常精确,但在周日频段存在比较明显的偏差,需要进一步精化.对于中国大陆的大地测量观测,固体潮可以采用Dehant等考虑地球内部介质非弹性和非流体静力平衡建立的固体潮理论模型或Xu 等基于全球SG观测建立的重力潮汐全球实验模型作为参考和改正模型,海潮负荷效应应该采用Nao99作为改正模型.     

武汉九峰地震台超导重力仪观测分析研究

连续重力观测和GPS的技术结合能够监测到物质迁移和地壳垂直形变之间的量化关系.和相对重力测量以及绝对重力测量技术相比,其避免了时间分辨率和观测精度低,无法精细描述观测周期内的物质迁移过程问题.本文利用武汉九峰地震台超导重力仪SGC053超过13000 h连续重力观测数据;同址观测的绝对重力仪观测结果;气压数据;周边GPS观测结果;GRACE卫星的时变重力场;全球水储量模型等资料,采用同址观测技术、调和分析法、相关分析方法在扣除九峰地震台潮汐、气压、极移和仪器漂移的基础上,利用重力残差时间序列和GPS垂直位移研究物质迁移和地壳垂直形变之间的量化关系.结果表明:在改正连续重力观测数据的潮汐、气压、极移的影响后,不仅准确观测到2009年的夏秋两季由于水负荷引起的约(6~8)×10-8m·s-2短期的重力变化.而且在扣除2.18×10-8(m·s-2)/a仪器漂移和水负荷的影响后,验证了本地区长短趋势垂直形变和重力变化之间具有一致的负相关性规律.同时长趋势表明该地区地壳处于下沉,重力处于增大过程,增加速率约为1.79×10-8(m·s-2)/a.武汉地区重力梯度关系约为-354×10-8(m·s-2)/m.     

Noise limitations in the core-mode band of superconducting gravimeter data

In recent years, several claims of detection of weak harmonic signals in the sub-tidal band of frequencies of high-quality gravimeter data have been made. Here, we review our attempts at confirming Smylie's claim of detection of the Slichter mode of inner-core oscillation using the same four data sets he and his colleagues used. We further examine the homogeneous 2 year data sets obtained from the superconducting gravimeters sited in Strasbourg and at Cantley, Quebec. We show that the power spectra of tidal-reduced, pressure and slew-corrected gravity residuals from these data are indistinguishable from a random walk process except in narrow bands dominated by residual earth tides and by harmonics of the diurnal thermal-atmospheric tide. Such a random ‘brown’ residual could result from mechanical instabilities and electronic noise in the instruments themselves, from site-specific tectonic noise, from local hydrological variations and pier instabilities or it could derive from unaccounted for variability in the atmospheric pressure in the vicinity of the instruments. We argue that the homogeneous 2 year data sets show no evidence whatever of Slichter harmonics even though our preprocessing methods reduce ‘apparent noise levels’ in the core-mode band by almost an order magnitude in comparison with the four data sets originally employed by Smylie et al. and, in their first attempted confirmation of his claimed discovery, by Jensen et al.     

Parameters of the fluid core resonance inferred from superconducting gravimeter data

We have estimated the parameters of fluid core resonance (FCR) due to the nearly diurnal free wobble of the Earth's core based on the superconducting gravimeter (SG) data obtained at the following four observation sites; Esashi and Matsushiro in Japan, Canberra in Australia and Membach in Belgium. By fitting the tidal admittances normalized with the O1 wave at each site to a model of the damped harmonic oscillator, we obtained values of 429.66 ± 1.43 sidereal days, 9350–10,835, −4.828E−4 ± 3.4E−6, −3.0E−5 ± 4.5E−6 for the eigenperiod, the Q-value and the real and imaginary parts of the resonance strength, respectively. Our values obtained from only using the gravity data are very consistent with those inferred from the VLBI nutation data. Our study strongly indicates that the systematic difference between two estimations from the gravity and the nutation in particular for the Q-value, which has been shown in previous works, is mainly caused by the inaccurate correction for the ocean tide effects. The error in the ocean tide correction is discussed based on the comparison among four global ocean tide models; Schwiderski model (1980), NAO.99b (Matsumoto et al., 2000), CSR4.0 (Eanes and Bettadpur, 1994) and GOT99.2b (Ray, 1999).     

Detection of inertial gravity oscillations in the Earth's core with a superconducting gravimeter at Brussels

As stated by Crossley and Rochester the discovery and identification of even a few undertone periods would provide valuable gross constraints on the stability profile and hence, the thermal regime of the core. Three years of continuous measurements of gravity variations with a superconducting gravimeter at Brussels offer a possibility to detect such oscillations in that part of the spectrum where tidal contributions can be eliminated. Spectral peaks of 10–15 ngal amplitude have been detected around a period of 13.9 h which were enhanced immediately after the two deep earthquakes of magnitude 7.2 which occurred during the observation interval.     

Least squares self-coherence for sub-nGal signal detection in the superconducting gravimeter records

In this paper, we use superconducting gravimeter (SG) data recorded at three stations of the global geodynamics project (GGP) network, with good geographical distribution, to search for possible significant peaks in the gravity spectrum that are in the assumed period range of the Slichter triplet. Seven-year long series from Cantley (Canada), and four-year long series from both Canberra (Australia) and Moxa (Germany) stations are used. First, a solid Earth and ocean tide model is subtracted from the data, followed by a local atmospheric pressure correction based on a frequency-, and location-dependent admittance estimated by the least squares response method. Subsequently, the residual series are filtered with a Parzen-based bandpass filter with a passband (12 h–78 s). A sub-nGal detection level is confirmed by injecting an artificial sine wave of different amplitudes. The Least Squares Self-Coherency spectrum shows the existence of many apparently statistically significant peaks at the 95% confidence level in the band (3–8 h). Although a few peaks are close to the claimed Slichter periods in previous research, the large number of candidate peaks may be related to other mechanisms such as global pressure variations, or hydrology.     

Results from 44 months of observations with a superconducting gravimeter at Moxa/Germany

Results for more than 42 months of observations with the superconducting gravimeter CD-034 at the Geodynamic Observatory Moxa are discussed. Moxa observatory is one of the newer stations within the ‘Global Geodynamics Project’ (GGP). A special feature of the gravimeter at Moxa is its dual sensor system; differences in the results obtained from the two sensor recordings are generally well within the standard deviations of the tidal analysis. One significant difference concerns the slightly different drift rates of 31 and 49.5 nm/s² per year for upper and lower sensor; both sensor drifts can be fitted by a linear function. We find that the noise levels are close to the ‘New Low Noise Model’ for the seismic-modes and are also low in the tidal bands. Due to this low noise, Moxa is a station well suited to search for small geodynamic signals. The long-period variation in the gravity residuals correlates well with the polar motion.The difference signal between the two sensor recordings has a peak-to-peak amplitude of about 6 nm/s² and shows systematic variations. Its spectrum is characterised by instrumental noise between 0.2 and 0.4 cph. The noise level of the difference and of the sum of the two residual datasets are clearly lower, respectively, higher than the noise contents of the gravity residuals themselves. This is a strong indication for the existence of broadband signals common to the two residual datasets, leading to the conjecture that the reduction of environmental effects is still not sufficient.Our results once more emphasize the necessity to correct the data for barometric pressure effects when analyzing the data for seismic modes. The reduction visibly increases the signal-to-noise ratio in the low frequencies of the mode band and helps to avoid misinterpreations of peaks. Besides the well known barometric pressure influence we can establish hydrological effects in the data which are probably caused by soil moisture and groundwater table variations as well as by batch-wise water movement within the weathering layer. As the major part of the observatory surroundings is above gravimeter level, an anticorrelation between hydrological and gravity changes is observed. In addition, it can be shown that global hydrological effects reach an order of magnitude that makes it necessary to consider these effects when analyzing long-period signals like polar motion. Vice versa these effects are large enough to be detectable in the gravity data. A first joint analysis of five datasets from the GGP network shows no indications for signals related to the Slichter triplet or core modes.     

Least-squares self-coherency analysis of superconducting gravimeter records in search for the Slichter triplet

We develop a new approach for the spectral analysis of the superconducting gravimeter data to search for the spheroidal oscillation ₁S₁ of the Earth solid inner core. The new method, which we call least-squares (LS) self-coherency analysis, is based on the product of the least-squares spectra of segments of the time series under consideration. The statistical foundation of this method is presented in the new least-squares product spectrum theorem that establishes rigorously confidence levels for detecting significant peaks. We apply this approach along with a number of other innovative ideas to a 6-year long gravity series collected at the Canadian Superconducting Gravimeter Installation (CSGI) in Cantley, Canada, by splitting it into 72 statistically independent monthly records. Each monthly record is analysed spectrally and all monthly LS spectra are multiplied to construct the self-coherency spectrum of the 6-year gravity series. The self-coherency spectrum is then used to detect significant peaks in the band 3-7 h at various significant levels with the aim to identify a triplet of periods associated with the rotational/ellipsoidal splitting of ₁S₁ (Slichter triplet). From all the Slichter periods predicted by various researchers so far, Smylie's triplet appears to be the most supported one, albeit very weakly, both, before and after the atmospheric pressure effect is removed from the series. Using the viscous splitting law [Smylie, D.E., 1992. The inner core translational triplet and the density near Earth's center. Science 255, 1678-1682] as guide, we can also see one interesting and statistically significant triplet with periods A = {4.261 h, 4.516 h, 4.872 h}, which changes slightly to A′ = {4.269 h, 4.516 h, 4.889 h} after the atmospheric pressure correction is applied to the gravity series.     

Incessant excitation of the Earth's free oscillations: global comparison of superconducting gravimeter records

Records of superconducting gravimeters (SGs) at Canberra (Australia), Esashi (Japan), Metsähovi (Finland) and Syowa Stations (Antarctica) were analyzed to search for further evidence of background free oscillations of the Earth. Spectrograms for 1-year period and averaged power spectra for seismically quiet periods were obtained for each of the stations. Anomalous features of the oscillations observed at Syowa Station, such as an apparent seasonal variation and a high intensity at frequencies between 3 and 4 mHz, were absent at the other SG stations. Among the SG stations used in this study, the background free oscillations were detected most consistently and distinctly at Canberra, where the noise level was comparable to that at the IDA quietest station, while that at Syowa Station was close to the critical limit for detecting the oscillations. The background free oscillations provide a good reference to evaluate the noise level in the milliHertz band.     

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.     

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

利用武汉台站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全球海潮模型仍可作为大地测量研究中的重要参考模型.     

Complexity as a streamflow metric of hydrologic alteration

We explore the potential of using a complexity measure from statistical physics as a streamflow metric of basin-scale hydrologic alteration. The complexity measure that we employ is a non-trivial function of entropy. To determine entropy, we use the so-called permutation entropy (PE) approach. The PE approach is desirable in this case since it accounts for temporal streamflow information and it only requires a weak form of stationarity to be satisfied. To compute the complexity measure and assess hydrologic alteration, we employ daily streamflow records from 22 urban basins, located in the metropolitan areas of the cities of Baltimore, Philadelphia, and Washington DC, in the United States. We use urbanization to represent hydrologic alteration since urban basins are characterized by varied and often pronounced human impacts. Based on our application of the complexity measure to urban basins, we find that complexity tends to decline with increasing hydrologic alteration while entropy rises. According to this evidence, heavily urbanized basins tend to be temporally less complex (less ordered or structured) and more random than basins with low urbanization. This complexity loss may have important implications for stream ecosystems whose ability to provide ecosystem services depend on the flow regime. We also find that the complexity measure performs better in detecting alteration to the streamflow than more conventional metrics (e.g., variance and median of streamflow). We conclude that complexity is a useful streamflow metric for assessing basin-scale hydrologic alteration.     

Loading effect of a self-consistent equilibrium ocean pole tide on the gravimetric parameters of the gravity pole tides at superconducting gravimeter stations

The gravimetric parameters of the gravity pole tide are the amplitude factor δ, which is the ratio of gravity variations induced by polar motion for a real Earth to variations computed for a rigid one, and the phase difference κ between the observed and the rigid gravity pole tide. They can be estimated from the records of superconducting gravimeters (SGs). However, they are affected by the loading effect of the ocean pole tide. Recent results from TOPEX/Poseidon (TP) altimeter confirm that the ocean pole tide has a self-consistent equilibrium response. Accordingly, we calculate the gravity loading effects as well as their influence on the gravimetric parameters of gravity pole tide at all the 26 SG stations in the world on the assumption of a self-consistent equilibrium ocean pole tide model. The gravity loading effect is evaluated between 1 January 1997 and 31 December 2006. Numerical results show that the amplitude of the gravity loading effect reaches 10⁻⁹ m s⁻², which is larger than the accuracy (10⁻¹⁰ m s⁻²) of a SG. The gravimetric factor δ is 1% larger at all SG stations. Then, the contribution of a self-consistent ocean pole tide to the pole tide gravimetric parameters cannot be ignored as it exceeds the current accuracy of the estimation of the pole tide gravity factors. For the nine stations studied in Ducarme et al. [Ducarme, B., Venedikov, A.P., Arnoso, J., et al., 2006. Global analysis of the GGP superconducting gravimeters network for the estimation of the pole tide gravimetric amplitude factor. J. Geodyn. 41, 334–344.], the mean of the modeled tidal factors δ_m = 1.1813 agrees very well with the result of a global analysis δ_CH = 1.1816 ± 0.0047 in that paper. On the other hand, the modeled phase difference κ_m varies from −0.273° to 0.351°. Comparing to the two main periods of the gravity pole tide, annual period and Chandler period, κ_m is too small to be considered. Therefore, The computed time difference κ_L induced by a self-consistent ocean pole tide produces a negligible effect on κ_m. It confirms the results of Ducarme et al., 2006, where no convincing time difference was found in the SG records.     

Analysis of observations with dual sensor superconducting gravimeters

Among the 21 superconducting gravimeters presently operating worldwide four instruments exist that are equipped with two vertically aligned sensor units. Three of the instruments are installed in Germany (Bad Homburg, Moxa, Wettzell) and one in South Africa (Sutherland). Comparisons of the data sets obtained with the dual sensor systems yield information on instrumental effects and sensitivity as well as on the efficiency of reductions of environmental effects applied to the data. The latter is an important constraint when looking for small geodynamic signals like Slichter and core modes or aperiodic variations.From analyses of the two data sets of each instrument a small but significant difference of 1-3% in the response of the sensor units on barometric pressure variations is found. Likewise, the records of lower and upper sensor vary slightly but not systematically with regard to the noise levels in the different frequency ranges. The tidal analyses yield an agreement of the tidal parameters generally well within the standard deviations determined from the least squares adjustment in the tidal analysis. The deviations are in the range between 0×10⁻⁴ and 3×10⁻⁴ for the amplitude factor and the phases differ between 0.0005° and 0.01° for the four main tidal constituents O1, K1, M2, and S2.The comparison of the gravity residuals of the two sensors with each other as well as with their sum and difference in the time and frequency domain shows the existence of identical signals in the records of the two sensors in the whole range of observation. This probably means that either the environmental reductions applied are not sufficient or there are additional disturbing effects in the data which have not been taken care of yet. From the study it emerges that it is not possible to get entirely rid of the tidal signals in the data. This is probably also due to the fact that despite reductions the data sets contain additional signals and slightly different noise at tidal frequencies which affect the result of the tidal analysis.     

Detection of spheriodal free oscillation excited by Peru 8.2 M_s earthquake with five international superconducting gravimeter data

In this paper, authors obtain the spectral peaks of the earth free oscillation and check all normal modes from 0S0 to 0S48 accurately, with the Fourier analysis and the maximum entropy spectrum method dealing jointly with six groups of the observational residual data from five international superconducting gravimeter stations. By comparing the observational results in this paper with three former groups of observations or models, authors notice an extra discrepancy between two observational 0S2 modes excited separately by Peru earthquake and Alaska earthquake, which probably mirrors the anisotropy of the Earth's inner core. The analysis on the splitting 1S2 mode shows that the asymmetric factor of rotationwise spectral splitting is possible to be different from that of anti-rotationwise spectral splitting.     

Burris型重力仪性能分析

对Burris型相对重力仪长时间断电后再次供电,分析Burris型重力仪升温过程,并使用2台Burris型重力仪（B95和B101） 2014年7月至2016年12月的实测数据,分析仪器性能。结果表明：① Burris型相对重力仪长时间断电并再次供电后,加热125 min达到恒温点,在升温过程中,读数先大幅减小后小幅增大,绝对变化率逐渐减小,通电110 h后零漂小于3×10^-8 m/（s²·h）,达地震重力观测仪器进网要求;② Burris型重力仪静态零漂存在由负变正的过程,不同仪器所需时间不同;③随着仪器使用年限加长,Burris型重力仪稳定性变差。忽略重力仪突跳读数,可获取测量精度在10×10^-8 m/s²左右的流动重力数据;④ Burris型重力仪混合零漂和静态零漂变化基本一致,呈逐渐增大趋势,但存在差别。     

Wavelet filter analysis of splitting and coupling of seismic normal modes below 1.5 mHz with superconducting gravimeter records after the December 26, 2004 Sumatra earthquake

New generation superconducting gravimeters (SGs), which have been demonstrated to be better than the best seismometers STS-1 at frequencies below 1 mHz, can be accepted as the quietest vertical seismometers for observation of long-period earth free oscillations. Wavelet filtering with narrow band-pass frequency response as shown in this paper is very helpful in removing at- mospheric pressure effects from on gravity records in long-period seismic mode frequency bands. The processing of high quality SG records after the great Sumatra earthquake (Dec. 26, 2004) with wavelet filtering leads to clear observations of all coupled toroidal modes below 1.5 mHz except these for 0T5, 0T7 and 1T1; moreover 1T2 and 1T3 are, for the first time, unambiguously revealed in the vertical components of the free oscillations. The three well-resolved splitting singlets of overtones 2S1 are observed from a single SG record for the first time.     

Impulse response of the kinematic diffusion model as a probability distribution of hydrologic samples with zero values

It is hypothesized that the unit impulse response of a linearized kinematic diffusion (KD) model is a probability distribution suitable for frequency analysis of hydrologic samples with zero values. Such samples may include data on monthly precipitation in dry seasons, annual low flow, and annual maximum peak discharge observed in arid and semiarid regions. The hypothesized probability distribution has two parameters, which are estimated using the methods of moments (MOM) and maximum likelihood (MLM). Also estimated are errors in quantiles for MOM and MLM. The distribution shows an equivalency of MOM and MLM with respect to the mean value—an important property for ML-estimation in the case of the unknown true distribution function. The hypothesized KD distribution is tested on 44 discharge data series and compared with the Muskingum-like (M-like) probability distribution function. A comparison of empirical distribution with KD and M-like distributions shows that MOM better reproduces the upper tail of the distribution, while MLM is more robust for higher sample values and more conditioned on the value of the probability of the zero value event. The KD-model is suitable for frequency analysis of short samples with zero values and it is more universal than the M-like model as its modal value cannot be only equaled to zero value but also to any positive value.