Age of the source of the Jarrafa gravity and magnetic anomalies offshore Libya and its geodynamic implications

The interpretation of the Jarrafa magnetic and gravity highs, NW Libyan offshore, suggests that it may be caused by a body of high-density and high magnetization. Analysis of their power spectra indicates two groups of sources at: (1) 2.7 km depth, probably related to the igneous rocks, some of which were penetrated in the JA-1 borehole, (2) 5 km depth, corresponding to the top of the causative body and (3) 10 km depth, probably referring to the local basement depth. The boundary analysis derived from applied horizontal gradient to both gravity and magnetic data reveals lineaments many of which can be related to geological structures (grabens, horsts and faults).The poor correlation between pseudogravity fields for induced magnetization and observed gravity fields strongly suggests that the causative structure has a remanent magnetization (D = −16°, I = 23°) of Early Cretaceous age, fitting with the opening of the Neo Tethys 3 Ocean.Three-dimensional interpretation techniques indicate that the magnetic source of the Jarrafa magnetic anomaly has a magnetization intensity of 0.46 A/m, which is required to simulate the amplitude of the observed magnetic anomaly. The magnetic model shows that it has a base level at 15 km.The history of the area combined with the analysis and interpretation of the gravity and magnetic data suggests that: (1) the source of the Jarrafa anomaly is a mafic igneous rock and it may have formed during an Early Cretaceous extensional phase and (2) the Jarrafa basin was left-laterally sheared along the WNW Hercynian North Graben Fault Zone, during its reactivation in the Early Cretaceous.     

Interpretation of gravity and magnetic data in the Phlegraean fields geothermal area, Naples, Italy

Available gravity and magnetic data of the Phlegraean Fields geothermal area, Naples, Italy, have been interpreted and the obtained structural models discussed in the light of the other available geological, volcanological and geophysical data.On the basis of the results of a previous seismic reflection survey in the Gulf of Naples and in the Pozzuoli Bay, which delineated a basement characterized by a seismic velocity of 4–6 km/s, it has been possible to evaluate the gravity anomaly connected with the morphology of this horizon ( = 2.7 g/cm³).The residual anomaly map, obtained after subtraction of the regional long-wavelength components relative to mantle and deep crustal structures and the computed components relative to the above-mentioned seismic basement, shows up as a circular low with an amplitude of 10 mgal centred in the Pozzuoli Bay. This gravity low has been interpreted as due to the occurrence, in the centre of Pozzuoli Bay, of light (Δ = −0.2 g/cm³) material with a maximum thickness of about 2 km. However, a contribution to the anomaly due to a narrow magmatic body intruded in the basement, as suggested by volcanological and ground deformation data, cannot be excluded.The aeromagnetic map of the Phlegraean Fields is characterized by three main anomalies which have been fitted by superficial tridimensional parallelepipedic bodies, schematically representing lava flows and domes. Their anomalies have been subsequently subtracted from the observed field, obtaining as a residual a large anomaly centred in the southwestern area of the Pozzuoli Bay. It has been interpreted as being due to a lowmagnetized body which, taking into account the thermal state of the area, should represent that part of the pyroclastic sequence which has lost part of its magnetization by thermo-chemical alteration.     

利用2009年日全食的精细重力观测探寻“引力异常”

日全食期间在是否存在"引力异常"现象是一个关系到引力本质的重大问题,也是一个十分具有现实意义的科学问题.近一个世纪来,许多科学工作者作了大量的观测和实验,研究日全食期间的"引力异常"现象.肯定的、否定的结果并存,更多的是难以确定的结果.在学习、总结和发展了前人工作的基础之上,本文利用7台LaCoste-Romberg重力仪和1台SG-053超导重力仪,对2009年7月22日日全食期间的重力变化进行了精细的观测.研究利用了目前最好的观测环境和观测仪器,无论观测规模还是观测质量在近百年来日全食的重力观测研究中都是首屈一指的.在经过详细的数据处理和图像分析后,可以确定,日全食期间不存在现有观测能力之上的引力异常,以往发现的异常现象受剧烈的气象变化影响可能性极大.另外,在处理重力观测数据时,本文合理的利用了最小二乘多项式拟合法去除理论潮汐值,使计算得以简化,结论更加可靠.     

攀西地区重力场特征及地壳密度结构

攀西地区位于峨眉山大火成岩省中西部,构造和岩浆特征显著,地震活动强烈.通过对野外重力测量得到的云县—会东和普洱—七甸两条剖面的高精度重力观测数据进行处理和分析,构建了沿剖面的二维地壳密度结构,其中普洱—七甸剖面与孟连—马龙宽角地震剖面部分位置重合.同时结合区域重力异常特征及下地壳视密度填图结果,得到如下初步认识:红河断裂带是南北地震带南段地区重要的构造分界线,断裂带南北向密度结构和莫霍面分布形态存在较大差异,沿走向构造变化.云县—会东剖面上大姚—会东段下地壳底部存在密度较高的壳幔过渡层,结合研究区下地壳底部壳幔过渡层的密度分布特征,认为该过渡层不是攀西裂谷下的"裂谷垫",而是由岩浆底侵作用造成的.     

基于重力异常反演的三维地震波速度结构分析方法研究

利用重力异常反演测试三维地震波速度结构,存在解不唯一、可靠性不高的问题。将面波反演充分融合到重力异常反演方程中,降低传统反演方法的非唯一性,并提升可靠性。以川滇地区为例,采用融合后的重力异常反演方法分析三维地震波速度结构。通过速度和密度的关系转换,得到对应的重力异常数据。由于面波频射数据主要对地震波横波速度敏感,因此将重力异常数据和初始横波速度相连,依据地震波速度和岩石密度之间的关系,获取重力异常反演方程,用于分析速度结构。选取21.6°~34.2°N、97.1°~105.9°E范围内的川滇地区活动块体作为实验数据,经过实验分析发现：使用该方法迭代反演川滇地区地壳上地幔顶部横波速度,重力异常数据和面波频射数据的残差值分别是6.24 mGal和0.027 km/s,实际拟合效果较好;分析该地区不同深度切面横波速度发现,在24 km深度处,上地壳中含有相对低速层,在44 km深度处,中下地壳中存在低速层;且该方法分析川滇地区三维地震波速度结构解的分辨率较高。     

gPhone重力仪的面波频段响应实测研究

以中国大陆构造环境监测网络昆明台和恩施台gPhone相对重力仪连续重力潮汐观测数据为基础,研究了gPhone重力仪在1 mHz以上频段的高频响应。从瑞利面波角度获得gPhone重力仪的高频响应,并且通过与同址观测的STS-1地震仪LHZ分量数据进行对比,验证了gPhone重力仪高频观测结果的可靠性。对比从gPhone重力仪和STS-1地震仪观测数据中提取到的面波波形和群速度频散曲线,发现昆明台两类仪器观测到面波信号的振幅和相位都较为一致,而恩施台仅振幅较为一致,相位上存在较明显差异,gPhone重力仪记录的面波信号在各频段存在不同的时间延迟。用两类仪器观测数据获得了大地震激发的自由振荡,结果表明两类仪器观测到的基频球型模态自由振荡的频率和振幅都吻合较好,进一步验证了gPhone重力仪对高频频段信号振幅响应的可靠性。以上研究结果表明:利用gPhone重力仪能够准确地观测到大地震激发的面波和自由振荡等高频信号的振幅,但在记录信号的相位信息时,有些仪器会有相位偏移产生,如果研究中需要考虑信号的相位,则必须获得仪器相位偏移量,再进行仪器相位校正。     

蓟县地震台重力仪观测数据质量与水文响应分析

对天津蓟县地震台2套相对重力仪2018—2019年观测数据进行潮汐分析和提取非潮汐重力残差处理,研究水文变化所引起的重力响应。结果表明,蓟县地震台2套重力仪精度较高;GS-15重力仪重力残差与水文变化间存在较好的相关性,具体表现为降水时重力减少,降水结束后重力增大,水位下降时重力减少,水位上升时重力增大,但其间都存在延迟对应关系。     

Gravity control network at airports in nigeria

A network of gravity of control bases located at most of the airports and airfields in Nigeria has been established as a means of controlling future gravimetric studies in Nigeriaand of integrating existing gravity surveys into a unified whole for the purpose of producing gravity anomaly maps for the country. The network will also provide a means of checking and adjusting gravimeter calibration differences to a common standard.The LaCoste and Romberg gravimeters, models G446, G464 and G468, used for the observations were transported between the base stations by means of a six-seater fixed-wing aircraft. The scale adopted of the network was the same as that of gravimeter G464 to which the other two meters were referenced. The absolute gravity values observed ranged from 978184.882 ± 0.014 mGal at Nguru to 977844.379 ± 0.020 mGal at Jos, yielding a gravity range of 340.503 mGal. The accuracy of the individual absolute gravity value is better than 0.030 mGal using a datum value of 978120.933 mGal for site 103401 estlished at the Kano International Airport. The 34 observation sites are illustrated with sufficient detail to permit reoccupation within 0.50 m to 1.0 m.     

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.     

气压变化对石英弹簧重力仪读数的影响

根据我们在低压仓中对国产ZSM—Ⅲ型(简称Z)、美国Wodnre(简称W)及加拿大Scintrex CG—2型(简称C)8台重力仪进行气压实验的结果,气压变化和重力仪读数呈极好的线性关系,但每台仪器有不同的回归系数,因此在以地震预报为目的的高精度重力测量中,应该对野外观测的结果进行气压改正。用高差所建立的重力基线场的数值是否准确有进一步考虑的必要。     

新疆精河6.6级地震周边地区密度构造、均衡异常以及岩石圈挠曲机理

依据EIGEN-6C4重力模型和ETOPO1高程模型数据,围绕新疆精河6.6级地震展开岩石圈均衡与挠曲机理研究,得到如下结论：（1）震中附近的布格与自由空气重力异常分别为-221和-92mGal(10~(-5 )m·s~(-2)）,震中位于重力异常高梯度带上;（2）震中周边地区地壳厚度约为50km,密度结构总体变化平缓,东西方向地壳厚度变化较小,但自南向北地壳厚度逐渐变薄,精河6.6级地震初始破裂发生在上中地壳分界面附近;（3）震中附近岩石圈承载的垂向构造应力为20MPa左右,震中位于岩石圈垂向构造应力极大值附近的高梯度带上;（4）地震周边地区岩石圈有效弹性厚度最优解为26km,加载比最优解为F_1=1,F_2=F_3=0,表明该区域岩石圈相对坚硬,且导致岩石圈变形的初始加载全部来自地表．     

Some critical factors for engineering and environmental microgravity investigations

The gravity method is one of the geophysical tools used for engineering and environmental investigations where the detection of cavities, karst phenomena, subsoil irregularities, or landfills is essential. In many cases, deep or small-scale heterogeneities generating low-amplitude anomalies have to be detected and the reliability of further interpretation requires highly accurate measurements, carefully corrected for any quantifiable disturbing effects. The purpose of this study is to investigate the factors likely to limit measurement quality and how to make improvements.Calibrations of a Scintrex gravimeter were made between French relative and absolute base stations, and the relative uncertainties on the calibration factors were estimated for these links. Ranging from 10⁻³, for calibration on an old gravity net, to 10⁻⁴, for a high amplitude absolute base line, this accuracy will be generally sufficient for microgravity surveys.Continuous gravity recordings of Scintrex gravimeters, installed at the same stable site, enabled the estimation of the stability and accuracy of the instruments and revealed that some of the time variations of g measurements, such as instrumental drift, tidal effects and seismic noise, are not entirely removed by standard processing procedures. The accuracy of corrected gravity measurements is mainly limited by inadequate corrections of tidal effects and by a poor estimation of ocean loading effects. In comparison with residual defaults in tidal corrections, instrumental and seismic noises are taken more properly into account by statistical data processing.In field operation, residual tidal effects are generally integrated into an experimental terrain drift estimated on the basis of frequent repeated measurements. A differential gravity approach, based on a fixed gravimeter reference whose recordings are used to correct measurements made with a mobile gravimeter, has also been investigated at a test site. Compared to standard processing, this method can help improve repeatability of gravity measurements.Microgravity surveys in the urban environment require effective and accurate consideration of the effects of infrastructures, nearby buildings and basements, as well as those of topography, in the vicinity of a gravity station. Correction procedures, applied at the same experimental site, where gravity points are located close to buildings, walls and basement slope, appear to have almost totally eliminated these disturbances.     

Effective barometric admittance and gravity residuals

In the analysis of surface gravity signals that may originate from the Earth's core, the step of correcting for the atmospheric pressure fluctuations is one that must be done carefully. We apply two techniques for determining the local, or effective, barometric admittance function between simultaneous observations of surface gravity and pressure. The first is a frequency domain fit that computes the admittance on a band-by-band basis. Using data from both the Canadian and French superconducting gravimeters we determined that the magnitude of the local, or background, admittance increases smoothly and monotonically from about 0.2 μgal mbar⁻¹ at long periods (> 10 days) to about 0.35 μgal mbar⁻¹ at frequencies greater than 3 cycles per day (c.p.d.); the phase lag is within a few degrees of 180°. By comparison, the effective admittances of the large-scale harmonics of the solar heating tide (S₁---S₇) are much smaller, between 0.1 and 0.3 μgal mbar⁻¹, for most of the harmonics of a day. In the second approach we fit a symmetrical time domain admittance function having lengths between 1 and 19 h using both a standard least-squares fit to a white noise residual and a new, and clearly superior, fit assuming a brown noise residual. Both time and frequency domain approaches give comparable results and contribute to a significant lowering of the residual level in non-tidal bands.     

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.     

Some space gravity formulas and the dimensions and the mass of the Earth

Summary Adopting thePizzetti-Somigliana method and using elliptic integrals we have obtained closed formulas for the space gravity field in which one of the equipotential surfaces is a triaxial ellipsoid. The same formulas are also obtained in first approximation of the equatorial flattening avoiding the use of the elliptic integrals. Using data from satellites and Earth gravity data the gravitational and geometric bulge of the Earth's equator are computed. On the basis of these results and on the basis of recent gravity data taken around the equator between the longitudes 50° to 100° E, 155° to 180° E, and 145° to 180° W, we question the advantage of using a triaxial gravity formula and a triaxial ellipsoid in geodesy. Closed formulas for the space field in which a biaxial ellipsoid is an equipotential surface are also derived in polar coordinates and its parameters are specialized to give the international gravity formula values on the international ellipsoid. The possibility to compute the Earth's dimensions from the present Earth gravity data is the discussed and the value ofMG=(3.98603×10²⁰ cm³ sec^–2) (M mass of the Earth,G gravitational constant) is computed. The agreement of this value with others computed from the mean distance Earth-Moon is discussed. The Legendre polinomials series expansion of the gravitational potential is also added. In this series the coefficients of the polinomials are closed formulas in terms of the flattening andMG.Publication Number 327, and Istituto di Geodesia e Geofisica of Università di Trieste.     

ANALYSIS ON STRUCTURAL FEATURES:INSIGHTS FROM SATELLITE-DERIVED GRAVITY MODEL AND SEISMIC PATTERN

In this paper, based on a large number of cumulative observational data from the seismic monitoring network in China, we grid the research area to calculate the density values at each grid node and convert the qualitative earthquake epicenter distribution to quantitative seismic pattern. Minimum magnitude of completeness(M_C)is determined by magnitude-rank analysis, which provides lower limit earthquake and original time. New satellite-derived gravity model v23.1, which is based on satellites CryoSat-2 and Jason-1 data, is used to determine the Bouguer gravity anomaly derived from free-air gravity anomaly and elevation database sets SRTM30, and ultimately, the complete Bouguer correction is obtained. In this paper, the Xingtai earthquake zone and Tanlu fault zone (Anhui segment) are selected for case study. Bouguer gravity anomaly presents a NE-trending U-shaped narrow strip in the Xingtai earthquake zone, and its location is consistent with Shulu Fault Basin. Grid density value contours are restricted by the U-shaped strip, and the extreme value of seismic activity density lies in the bottom of the U-shaped strip as shown in the cross section. The results of Bouguer gravity anomaly and upward continuations to the different heights show good linearity and gradient in the Tanlu fault zone (Anhui segment); and both long-axis direction of seismic pattern and nodal plane strike of seismogenic fault from focal mechanism solutions trend NNE. In short, the Tanlu fault zone(Anhui segment)is a large deep-seated fault that still has the ability to control seismic activity along it. Based on the measured gravity and magmatic data, using the edge detection TDX method to interpret the concealed boundary of the Anqing M4.8 earthquake near the Tanlu fault, and combining with the results from deep seismic reflection profiles of the study area, we discussed the causative fault of the Anqing earthquake.     

The application of combined gravity and seismic data formation separation for revealing deep structure

In some oilfields with 3D seismic data, the deeper structure cannot be observed due to poor quality deep seismic data. Layer stripping using both seismic and gravity data is a solution for this problem but it cannot get satisfactory results because the horizontal variations in formation density are ignored. We present a variable-density formation separation technique to address this problem. Based on 3D seismic depth data and laterallyvariable density derived from 3D seismic velocity data, the upper formation gravity effect is calculated by forward modeling and removed from the Bouguer gravity. The formation-separated gravity anomaly with variable density is obtained, which mainly reflects the deeper geological structure. In block XX of North Africa, the shallow formations seismic data is excellent but the data at the top of basement is poor. The formation-separated gravity anomaly processed under the control of 3D seismic data fits well with the known seismic interpretation and wells. It makes the geological interpretation more reliable.     

重力仪检定基线的建立及其误差讨论

1978年5月使用了18台精度较高的重力仪在南京和北京分别建立了地震系统的重力检定场.本文结合以上成果和有关重力比较测定的资料, 从分析格值变化的各种因素及其产生的误差入手, 说明仪器格值标定的重要性和具体要求, 提出了选建地震系统重力检定场的方法技术和减小格值误差、提高测量精度的可能途径.      

重力剖面和激电测深联合勘测在涡河断裂研究中的应用

为查明涡河断裂的位置、性质等参数,选择重力剖面和激电测深联合勘测方法,对指定区域进行物探勘测工作。结果表明:由密度差异引起的重力异常位置和岩石电阻率、极化率异常位置基本吻合,说明以重力剖面和激电测深为手段的联合勘测是查找断裂的有效方法。     

重力-地震联合变密度剥层技术研究深层构造

在某些已经完成三维地震勘探工作的油气田区,因深层地震资料不好,深层构造难以得到解决。虽然结合地震资料与重力资料的重力剥皮技术得到开发和使用,但由于它没有考虑地层密度的横向变化,因此未能获得满意的处理效果。为了改善该类地区的重力勘探效果,本文提出了变密度剥层方法,它利用各地层的三维地震深度资料和三维地震层速度转换的具横向变化的密度资料,通过正演其重力效应,并从布格重力异常中去除,从而得到主要反映深部地质构造信息的变密度剥层重力异常。在北非某区块,除基底顶面外,各层地震资料良好,通过重力-地震联合变密度剥层处理得到的剥层重力异常,与已知地震和钻井资料吻合良好,为深部基底研究奠定了良好的基础。