The first absolute gravity measurements in Indonesia

For the purposes of the calibration of the superconducting gravimeter (SG) in Bandung and the establishment of the absolute gravity (AG) points, we carried out AG measurements for the first time in Indonesia in November 2002. The measurements in Bandung were conducted between November 15th and 20th by means of a FG5 (#210), and 14,520 effective drops were obtained. The gravity value newly determined at the AG point in Bandung is 977976701.2 μgal (1 μgal = 10⁻⁸ ms⁻²) and the scale factor for the SG is −52.22 μgal/V. We also established another AG point in Yogyakarta near Merapi volcano and carried out AG measurements in Yogyakarta between November 22nd and 26th. The gravity value determined for this station is 978203093.5 μgal.     

AGrav—The New International Absolute Gravity Database of BGI and BKG and its benefit for the Global Geodynamics Project (GGP)

A new database for absolute gravity (AG) measurements has been implemented at BGI and BKG and is operational now for storing absolute gravity data either in the form of metadata or as detailed measurement results. The database development was proposed by the IGFS (International Gravity Field Service) and is expected to have a great importance for the GGOS (Global Geodetic Observing System) initiative. This database will provide an overview about AG stations and observations and by this improve the cooperation between gravity groups and foster the combination with other geodetic observation techniques. The international community of absolute gravimeter users is asked to contribute to this database.In addition to its primary purposes, demonstration of the global site distribution and information about available observations, the database could also provide an important contribution to the Global Geodynamics Project (GGP). Precise repeated absolute gravity measurements at the superconducting gravimeter (SG) sites are necessary for the determination of SG drift parameters and can be used for checking SG instrument calibration factors. The AGrav database is capable of storing the necessary AG observations at the SG location in detail up to the “single drop level” and provides this information for the combination with SG time series. An example for a selected station is presented. It is proposed to establish an interface between the AGrav and GGP databases.     

武汉九峰站地下水变化对重力场观测的影响

本文从地下水渗透过程的物理机制出发,采用一维水动力学模拟,利用井水位和降雨数据模拟计算了武汉九峰站附近的土壤含水率变化,在此基础上估计了地下水变化导致的重力效应.其峰对峰变化幅度达到15.94μGal,说明当利用精密重力观测研究长周期效应时实施台站地下水改正的必要性;频域分析表明,地下水重力效应在周年频段上的振幅最大,说明地下水变化对重力的最大影响来自季节性变化.对比模拟计算的地下水重力效应和经过潮汐、大气、极移等改正后的高精度超导重力残差(峰对峰变化幅度为12.73μGal),发现两者在时域和频域均具有良好的一致性,说明超导重力残差信号主要来源于局部地区地下水的变化,同时也验证了本文使用的水动力学模拟方法的正确性.     

Hydrological signals in height and gravity in northeastern Italy inferred from principal components analysis

This work describes a study of GPS heights, gravity and hydrological time series collected by stations located in northeastern Italy. During the last 12 years, changes in the long-term behaviors of the GPS heights and gravity time series are observed. In particular, starting in 2004–2005, a height increase is observed over the whole area. The temporal and spatial variability of these parameters has been studied as well as those of key hydrological variables, namely precipitation, hydrological balance and water table by using the Empirical Orthogonal Functions (EOF) analysis. The coupled variability between the GPS heights and the hydrological balance and precipitation data has been investigated by means of the Singular Value Decomposition (SVD) approach. Significant common patterns in the spatial and temporal variability of these parameters have been recognized. In particular, hydrology-induced variations are clearly observable starting in 2002–2003 in the southern part of the Po Plain for the longest time series, and from 2004–2005 over the whole area. These findings, obtained by means of purely mathematical approaches, are supported by sound physical interpretation suggesting that the climate-related fluctuations in the regional/local hydrological regime are one of the main contributors to the observed variations. A regional scale signal has been identified in the GPS station heights; it is characterized by the opposite behavior of the southern and northern stations in response to the hydrological forcing. At Medicina, in the southern Po Plain, the EOF analysis has shown a marked common signal between the GPS heights and the Superconducting Gravimeter (SG) data both over the long and the short period.     

Gravity field variations from superconducting gravimeters for GRACE validation

The network of superconducting gravimeters (SG) of the ‘Global Geodynamics Project’ (GGP) offers the unique opportunity to supplement and validate the gravity field variations derived from the GRACE satellite mission. Because of the different spatial and temporal resolutions of the gravity data a combination of all datasets can be used to retrieve a maximum of information regarding mass transfers especially related to hydrology which is deployable as constraint for hydrological modelling.For a consistent combination of the datasets the gap between terrestrial data of superconducting and absolute gravimeters (AG) and from satellite data has to be bridged. A successful combination of SG and AG data could be realized for several stations which resulted in time series of the highest accuracy and long-term stability.In principle, the same reductions applied to GRACE data have to be taken into account for the terrestrial data. The separation of local hydrological effects in SG observations is crucial for the comparison with satellite-derived gravity data. It is shown that even for stations with a hydrological challenging situation such as Moxa/Germany local hydrology-induced effects can be successfully modelled.Currently, the study focuses on Europe with its dense and long-term observation network. Regarding the consistency of the SG gravity variations they are representative for a larger region. From a comparison with GRACE-derived gravity field changes, and the variations due to hydrological models a principle good agreement emerges.     

Accurate determination of calibration factor for tidal gravity observation of a GWR-superconducting gravimeter

Introduction The superconducting gravimeters (SGs) are widely used to observe the change in temporal and spatial gravity field by many countries along the world, considering their characteristics as of the high-precision, good continuity and stability. By using the temporal variation of the tidal grav-ity obtained from the global network of SGs, the Global Geophysics Project (GGP) aims at some hot problems in geophysics and geodynamics (SUN, HSU, 1997). The original observations of the c…     

GWR超导重力仪潮汐观测标定因子的精密测定

超导重力仪被普遍认为是迄今为止观测精度高、连续性和稳定性最好的相对重力仪，但是在获得真实重力场变化之前，必须利用格值（标定因子）对这类仪器的直接输出量（电压变化）进行标定．研究表明，格值的准确度将影响到观测数据的后续分析和资料的解释，本文利用两期（每期３天）FG５绝对重力仪在武汉国际重力潮汐基准站上的同址测量资料和该台站已知潮汐参数，详细研究了超导重力仪标定因子及其精度．      

Tackling mass redistribution phenomena by time-dependent GRACE- and terrestrial gravity observations

Time variable gravity field models derived from the satellite mission GRACE have been demonstrated to be consistent with water mass variations in the global hydrological cycle. Independent observations are provided by terrestrial measurements. In order to achieve a maximum of reliability and information gain, ground-based gravity observations may be deployed for comparison with the gravity field variations derived from the GRACE satellite mission. In this context, the data of the network of superconducting gravimeters (SG) of the ‘Global Geodynamics Project’ (GGP) are of particular interest. This study is focused on the dense SG network in Central Europe with its long-term gravity observations. It is shown that after the separation and reduction of local hydrological effects in the SG observations especially for subsurface stations, the time-variable gravity signals from GRACE agree well with the terrestrial observations from the SG station cluster.Station stability of the SG sites with respect to vertical deformations was checked by GNSS based observations. Most of the variability can be explained by loading effects due to changes in continental water storage, and, in general, the stability of all stations has been confirmed.From comparisons based on correlation and coherence analyses in combination with the root mean square (RMS) variability of the time series emerges, that the maximum correspondence between the SG and GRACE time series is achieved when filtering the GRACE data with Gaussian filters of about 1000 km filter length, which is in accordance with previous publications.Empirical Orthogonal Functions (EOF) analysis was applied to the gravity time series in order to identify common characteristic spatial and temporal patterns. The high correspondence of the first modes for GRACE and SG data implies that the first EOF mode represents a large-scale (Central European) time-variable gravity signal seen by both the GRACE satellites and the SG cluster.     

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.     

Integration of micro-gravity and geodetic data to constrain shallow system mass changes at Krafla Volcano, N Iceland

New and previously published micro-gravity data are combined with InSAR data, precise levelling and GPS measurements to produce a model for the processes operating at Krafla volcano, 20 years after its most recent eruption. The data have been divided into two periods: from 1990 to 1995 and from 1996 to 2003 and show that the rate of deflation at Krafla is decaying exponentially. The net micro-gravity change at the centre of the caldera is shown, using the measured free air gradient, to be −85 μGal for the first and −100 μGal for the second period. After consideration of the effects of water extraction by the geothermal power station within the caldera, the net gravity decreases are −73±17 μGal for the first and −65±17 μGal for the second period. These decreases are interpreted in terms of magma drainage. Following a Mogi point source model, we calculate the mass decrease to be ∼2×10¹⁰ kg/year reflecting a drainage rate of ∼0.23 m³/s, similar to the ∼0.13 m³/s drainage rate previously found at Askja volcano, N. Iceland. Based on the evidence for deeper magma reservoirs and the similarity between the two volcanic systems, we suggest a pressure-link between Askja and Krafla at deeper levels (at the lower crust or the crust-mantle boundary). After the Krafla fires, co-rifting pressure decrease of a deep source at Krafla stimulated the subsequent inflow of magma, eventually affecting conditions along the plate boundary in N. Iceland, as far away as Askja. We anticipate that the pressure of the deeper reservoir at Krafla will reach a critical value and eventually magma will rise from there to the shallow magma chamber, possibly initiating a new rifting episode. We have demonstrated that by examining micro-gravity and geodetic data, our knowledge of active volcanic systems can be significantly improved.Editorial responsibility: A. Harris     

Influence of soil consolidation and thermal expansion effects on height and gravity variations

The daily GPS height series of the Medicina station were analyzed for the period July 1996–September 2001. The station is located in the middle Po Plain on fine-grained alluvial deposits. A seasonal oscillation in the order of 18 mm (peak-to-peak amplitude) is present in the data. This crustal deformation has been modeled by including variations in the atmospheric, oceanic and hydrologic mass. The vertical positions can also be affected significantly by soil consolidation. Geotechnical parameters derived by in situ tests and laboratory analyses of the clayey soil collected at Medicina allowed the estimate of the soil settlement relevant to the seasonal oscillation of the surficial water table. Thermal expansion of the geodetic monument has to be taken into account in the case of high-precision vertical positioning. In this work models both for the soil consolidation and the thermal expansion effects are provided. The continuous gravity observations collected at Medicina by means of a superconducting gravimeter also exhibit a marked seasonal oscillation, which has been interpreted as the sum of loading and Newtonian attraction effects, as well as of the contribution due to soil consolidation. Especially the study concerning the soil consolidation effect has allowed a better insight on the seasonal vertical movements occurring at the Medicina station by providing quantitative information on soil behavior due to change of effective pressures. The results can be applied to those stations characterized by similar fine-grained soils and surficial hydrogeology.     

Recent geophysical investigation at Somma-Vesuvio volcanic complex

Activity at Somma-Vesuvio volcanic area in southern Italy is monitored by seismic stations and periodic geodetic and gravity surveys. The seismic network, which consists at present of four vertical stations and one three-component station, recorded an increase in earthquake activity in 1978 and between November 1988 and March 1989. During the later activity, earthquakes were located in a cluster about 3 km beneath the summit of the volcano. Two tide gauges, two tiltmeters and a recording gravimeter are also operating at Somma-Vesuvio. Yearly levelling surveys are conducted along several closed routes that extend from as much as 6 km from the base of the volcano to the summit area. Survey results reveal no significant ground movement since 1959, except for a slight subsidence around the rim of the summit crater. Gravity changes have been larger than the expected 10 μGal uncertainty of the measurements. The lack of contemporary elevation changes implies that the observed gravity changes are the result of a slight change in density structure. The cone of Somma-Vesuvio has been very stable for the last few decades, showing no indications of a buildup to activity. The lack of surface movement should rule out a magma-supply rate to this volcano at the historic eruptive rate of 0.002 km³/yr.     

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

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

2011年日本东北地区太平洋近海地震对亚洲板块及韩国大地控制网的影响分析

本研究分析了2011年3月11日发生的M_w9.0日本东北地区太平洋近海地震对亚洲地区和韩国国内GPS卫星常年跟踪站的位移影响.为此,利用了日本东北地区太平洋近海地震发生前后两周(2011年3月4日到3月18日)的GPS站点数据,包括震中附近地区(韩国,中国,中国台湾地区,日本和俄罗斯)55个GPS卫星常年跟踪站和284个IGS 全球跟踪站,并采用GAMIT/GLOBK软件进行处理和平差,估算出所有GPS站点的同震形变.结果显示,日本东北地区太平洋近海地震引起的同震形变影响在亚洲地区比较明显,包括日本和附近国家,距离震中2702 km的中国武汉(WUHN)站也观测到同震形变.为精确分析日本东北地区太平洋近海地震对韩国国家大地控制网的影响,通过GAMIT/GLOBK软件计算出韩国GPS卫星常年跟踪站之间的基线长度变形,并分析出弹性变形量.结果表明:大部分GPS站点均向震中方向膨胀,且向震中的垂直方向收缩.由日本东北地区太平洋近海地震导致的最大剪应变达到韩国国家大地控制网年均变形率的约7倍,对韩国的地壳产生14.5~57.7 mm的水平位移,并导致韩国国家大地控制网产生弹性变形.因此,在不及时更新维护韩国国家大地控制网的情况下,GPS测量成果将会发生最大20 mm的位置误差.     

Glacial isostatic adjustment: New developments from advanced observing systems and modeling

The measurement of glacial isostatic adjustment (GIA) is one of the key ways in which geophysicists probe the long-term mantle rheology and Pleistocene ice history. GIA models are also tied to global and regional relative sea-level (RSL) histories, to 20th century tide-gauge (TG) data and to space and terrestrial geodetic measurements. Two new types of observation are related to the high-resolution space–gravity data recovered from the Gravity and Climate Experiment (GRACE) satellite pair and the soon-to-be launched Gravity and Ocean Circulation Experiment (GOCE) with on-board three-component gradiometer. Gravity mapping has the unique capability of isolating those regions that lack isostatic equilibrium. When coupled with other space and terrestrial geodetic measurements, such as those of the Global Positioning System (GPS) networks and with multi-decade terrestrial gravity data, new constraints on GIA are in the offing and should soon illuminate new interpretations of ice-sheet history and mantle response. GIA studies also incorporate space-based altimetry data, which now provide multi-decadal coverage over continents, oceans and lakes. As we are approaching 72 monthly solutions of GRACE gravity coefficients for determining the Earth's secular component of gravity change over the continents, a new issue has surfaced: the problem of relying on interannual hydrological modeling to determine the hydrological contribution to the linear trend in the gravity field. Correctly extracting this contribution is germane to using the GIA-driven component for modeling solid-Earth and paleo-climatic parameters.Seismic and heat-flux-based models of the Earth's interior are emerging with ever higher levels of sophistication regarding material strength (or viscosity). A basic question raised is: how good are traditional Newtonian and non-Newtonian viscosity models that only allow radial variations of Earth parameters? In other words: under what circumstances must this assumption be abandoned for joint interpretations of new and traditional data sets. In this short review we summarize the issues raised in the papers forming this special issue (SI) dedicated to GIA.     

Temporal Gravity Variations near Shrinking Vatnajökull Ice Cap, Iceland

Repeated gravity measurements were carried out from 1991 until 1999 at sites SE of Vatnajökull, Iceland, to estimate the mass flow and deformation accompanying the shrinking of the ice cap. Published GPS data show an uplift of about 13 ± 5 mm/a near the ice margin. A gravity decrease of –2 ± 1 μGal/a relative to the Höfn base station, was observed for the same sites. Control measurements at the Höfn station showed a gravity decrease of –2 ± 0.5 µGal/a relative to the station RVIK 5473 at Reykjavík (about 250 km from Höfn). This is compatible, as a Bouguer effect, with a 10 ± 3 mm/a uplift rate of the IGS point at Höfn and an uplift rate of ~20 mm/a near the ice margin. Although the derived gravity change rates at individual sites have large uncertainties, the ensemble of the rates varies systematically and significantly with distance from the ice. The relationship between gravity and elevation changes and the shrinking ice mass is modelled as response to the loading history. The GPS data can be explained by 1-D modelling (i.e., an earth model with a 15-km thick elastic lithosphere and a 7·10¹⁷ Pa·s asthenosphere viscosity), but not the gravity data. Based on 2-D modelling, the gravity data favour a low-viscosity plume in the form of a cylinder of 80 km radius and 10¹⁷ to 10¹⁸ Pa·s viscosity below a 6 km-thick elastic lid, embedded in a layered PREM-type earth, although the elevation data are less well explained by this model. Strain-porosity-hydrology effects are likely to enhance the magnitude of the gravity changes, but need verification by drilling. More accurate data may resolve the discrepancies or suggest improved models.     

Analysis and comparison of the tidal gravity observations obtained with LCR-ET20 spring gravimeter

Introduction The tidal variation signals of the solid Earth recorded with the modern measuring techniques can be used to provide effectively the means to cognize the Earths motion, deformation and structure character, especially to monitor the material motion in the Earths interior. In the recent 20 years, the distribution character of the global gravity field has made great progress since the continuous and deep theoretical studies in geophysics and geodesy domains. The successful manufactur…     

Combined discrete and continuous gravity observations at Mount Etna

Systematic investigation of discrete gravity measurements has continued at Mount Etna since 1986. The network now covers an area of 400 km² with about 70 stations 0.5–3 km apart. Mass redistributions occurring at depths ranging between about 8 km below sea level and a few hundred metres below the surface (magma level changes within the shallower parts of the feeding conduits) have been identified from these data. Conventional (discrete) microgravity monitoring on a network of stations furnishes only instantaneous states of the mass distribution at continuously active systems. In order to obtain information on the rate at which the volcanic processes (and thus mass transfers) occur, three stations for continuously recording gravity where installed on Mount Etna in 1998. A 16-month long sequence from one of the continuously running stations (PDN, located 2 km from the active northeast crater at the summit of Etna volcano) is presented. After removing the effects of Earth Tide and tilt, the correlation of the residual gravity sequence with simultaneous recordings of meteorological parameters acquired at the same station was analysed. Once the meteorological effects have also been removed, continuous gravity changes are within 10 μGal of gravity changes measured using conventional microgravity observations at sites very close to the continuous station. This example shows how discrete and continuous gravity observations can be used together at active volcanoes to get a fuller and more accurate picture of the spatial and temporal characteristics of volcanic processes.     

Temporal variations of gravity in the Aswan region of Egypt

The Kalabsha region at the northwest side of Lake Aswan has been a subject of seismological and geodetic investigations since 1981, after the 5.5 M earthquake on 14 November. The subsequent seismic activity was related to the effect of water-level oscillations in the Lake. Geodetic measurements provided indications of regional stress and of the local character of vertical movements. As crustal dynamics represents a potential hazard with respect to the size of the lake, other techniques have been applied to the monitoring system. Repeated gravity test measurements were already made in the late 1980s and in 1994, but only in 1997 was a new epoch of long-term investigation initiated. The measured data showed sufficient accuracy of approx. 5 μGal. Already the first difference of the two campaigns of the years 1997 and 2000 showed a few important features. First of all, the temporal non-tidal changes of gravity exceed the level of confidence significantly. The changes were compared with older available geodetic data, and it was concluded, that the gravity changes were not related directly to the vertical movements of the surface. On the contrary, a water loading effect may be expected on the eastern side of the network, especially in the Kalabsha local net. The development of the stress field is considered to be the main source of observed gravity changes. In both the Kalabsha and Seiyal nets the changes differ according to the relative positions of the particular measuring points with respect to the faults.     

鲜水河断裂带北段GPS测量及其运动特征

鲜水河GPS监测网由18个沿鲜水河断裂带分布的点组成．1991年对该网进行了首次观测,1996和2000年进行了两期整网复测．本利用该网199l和2000年的二期GPS观测资料结合川西地区20世纪70年代以来的常规大地测量资料,计算给出了这一强震活动带现今运动图像：由测区的北西虾拉沱测点到南东的乾宁测点水平位移速率逐渐增大,而断层两盘的相对运动不明显．表明利用短边GPS测量可有效监视活动断裂带两侧一定区域的地壳运动．对比有形变测量资料以来川西地区6级以上地震分布的特点,认为可根据鲜水河断裂带断层两盘的相对运动和跨越断层的地壳整体运动的强弱这两种运动图像来研究区域强震的活动特征．