Towards a cm-geoid for Hungary: Recent efforts and results

Some steps were taken recently for Hungary aiming at the determination of geoid heights with a cm-accuracy. The present HGTUB98 gravimetric solution was based on terrestrial gravity data, height data and the EGM96 geopotential model, and was computed with the 1D Spherical FFT method. The gravity data were used in the area 45.5 ° ≤ϑ ≤ 49 °, 16 ° ≤ λ ≤ 23 °, the resolution of the grid was 30″ × 50″. The DTM used had a resolution of 1 km × 1 km.Our solution was evaluated using GPS/levelling data at 340 and 308 points respectively and at 138 vertical deflection points. We have compared our solution to the European EGG97 geoid solution, the gravimetric solution HGR97B developed by A. Kenyeres and the litospheric geoid solution by G. Papp. We have correlated our recent HGTUB98 solution to the Moho model of Central Europe. The comparison with GPS/levelling yielded respectively an accuracy of ±8.7 cm and ±4.4 cm (in terms of standard deviation) when a linear trend was removed. The comparison of the 1D planar FFT solution for the deflections of the vertical with 138 astrogeodetic deflections yielded an accuracy (in terms of standard deviation) of ±0.62″ and ±0.52″ for ξ and η, respectively.     

The high‐resolution gravimetric geoid of North Iberia: NIBGEO

The gravimetric geoid computed in the northern part of Iberia, is presented in this paper. This computation has been performed considering two study windows fitted to the areas with higher density of gravity data, to reduce the computation errors associated to the scarcity of gravity data, as much as possible. The bad influence of a bathymetry with poorer resolution than the topography is also reduced considering the smallest marine area possible. Moreover, the computation of this gravimetric model is based on the most recent geopotential model: EIGEN‐GL04C (obtained in 2006). The method used in the computation of the new gravimetric geoid has been the Stokes integral in convolution form. The terrain correction has been applied to the gridded gravity anomalies, to obtain the corresponding reduced anomalies. Also the indirect effect has been taken into account. Thus, a new geoid model has been calculated and it is provided as a data grid in the Geodetic Reference System of 1980, distributed for the northern part of Iberia from 40 to 44 degrees of latitude and ?10 to 4 degrees of longitude, on a 161 × 561 regular grid with a mesh size of 1.5′ × 1.5′. This new geoid and the previous geoid Iberian Gravimetric Geoid 2005, are compared with the geoid undulations measured for eight points of the European Vertical Reference Network (EUVN) on Iberia. The new geoid shows an improvement in precision and reliability, fitting the geoidal heights of these EUVN points with more accuracy than the previous geoid. Moreover, this new geoid has a smaller standard deviation (12.6 cm) than that obtained by any previous geoid developed for the Iberian area up to date. This geoid obtained for the northern part of Iberia will complement the previously obtained geoid for South Spain and the Gibraltar Strait area; both geoids jointly will give a complete picture of the geoid for Spain and the Gibraltar Strait area. This new model will be useful for orthometric height determination by GPS over this study area, because it will allow orthometric height determination in the mountains and remote areas, in which levelling has many logistic problems. This new model contributes to our knowledge of the geoid, but the surrounding areas must be better known to constrain the lithospheric and mantle models.     

全球碳酸锂贸易格局与供应危机传播研究

随着新能源产业的快速发展, 全球碳酸锂贸易量逐渐增长, 但其出口集中于少数国家, 这意味着全球碳酸锂供应存在一定的风险。为了深入了解全球碳酸锂贸易格局演变、供应危机传播路径及影响, 本文构建了2000—2021年全球碳酸锂贸易网络和级联失效模型, 概述了主要出口国贸易流向和贸易格局演变过程, 并模拟了不同供应危机来源的雪崩规模、传播路径及其影响。我们得到: (1)20多年来全球碳酸锂贸易集中度逐渐增加, 出口前三位国家保持稳定, 而主要进口大国变化明显, 贸易格局已由南美供应全球转变为南美供应东亚; (2)碳酸锂贸易供应危机的影响范围逐渐扩大, 间接传播逐渐占据主导地位; (3)智利供应危机对中日韩的影响程度逐渐扩大, 但对欧洲的影响时间最长; 阿根廷和中国的雪崩规模在上升、传播次数增加, 而美国则相反。本文将贸易格局演变与危机传播结合, 明确了全球碳酸锂主要贸易国智利、阿根廷、中国、德国等作为供应危机来源的影响程度及其传播路径。以上分析结果将有利于维护全球碳酸锂贸易稳定和主要贸易国采取应对策略缓解供应危机。     

Multiresolution tectonic features over the Earth inferred from a wavelet transformed geoid

Geoid signals give information about the underlying density structure and can be used to locate the source depth of the mass anomalies. Wavelet analysis allows a multiresolution analysis of the signal and permits one to zoom into a specific area bounded by a particular length scale. The ability of wavelets to resolve the geoid signal into individual wavelength components without losing the spatial information makes this method superior to the more common spherical harmonic method. The wavelet analysis allows one to zoom into a specific area and look at the regional geology. We have used a wavelet transform of the geoid to study the regional geology of Japan and the Philippine Plate, South America, Europe, North America, East Africa and the Middle East, India and the Himalayas, China and Southeast Asia, and Australia. By filtering the Earth’s geoid anomalies with 2-D Gaussian wavelets at various horizontal length scales, one can detect the subduction zones along South America, the Aleutians, and the western Pacific; the Himalayas; the Zagros Mountains; the Mid-Atlantic ridge; and the island chains of the mid-Pacific. We have processed geoid data with a horizontal resolution down to approximately 200 km. Using an adjustable wavelet, one can detect structures that can only be picked up visually with much higher resolution spherical harmonic gravity data. We have also looked at the wavelength at which the maximum signal occurs over a range of scales. This method, known as E-max and k-max, is especially effective for detecting plate tectonic boundaries and ancient suture zones along with areas of strong non-isostatic gravitational potential due to high differential stress. These areas are likely to be at high risk of earthquakes. These methods will be especially useful to future studies of the geoid potentials of other planets, such as Mars and Venus, since they will allow careful studies of the regional geology variations with geoid data of the resolution available from satellites.     

地面重力数据结合GPS/水准数据精化兴城物探实验区大地水准面

区域大地水准面的确定是GPS测量常需解决的问题。目前确定大地水准面的方法主要包括重力法、GPS水准几何法及组合法,其中组合法因其精度和可靠性都较高,常用于计算高精度区域大地水准面。高精度的大地水准面模型是组合法确定区域大地水准面的关键。在我国,EGM2008全球重力场模型精度和分辨率均高于此前的所有模型,研究基于该模型的组合法大地水准面精化具有重要的实践意义。笔者以吉林大学兴城教学实习基地物探实验区为例,基于实测重力数据、EGM2008重力场模型和GPS水准数据,采用组合法精化了区域大地水准面,比较了组合法大地水准面模型和无重力实测数据的几何法大地水准面模型的精度差异,分析了该方法在物探测量中的适用性。结果表明,实验区组合法大地水准面模型精度最高达到1.2 cm,并且误差分布区间较小,总体上精度和可靠性高于对比的几何方法,并且组合法和几何法获取的两种大地水准面模型均能满足大比例尺物探测量要求。EGM2008模型精度较高,故平坦地区使用组合法时,高密度的实测重力数据可能带来高频扰动,有可能降低EGM2008重力场模型本身的精度,所以重力数据采集过程中要顾及重力点的密度和空间分布。本文方法更适用于地形复杂的地区。     

Climate downscaling over southern South America for present-day climate (1970-1989) using the MM5 model. Mean,interannual variability and internal variability

This work focuses on evaluating the ability of the MM5 regional model to represent the basic features of present climate over South America. The spatial distribution of seasonal means and the interannual variability, as well as annual cycles for precipitation and near-surface temperature have been evaluated. The internal variability has also been investigated. The analysis has two objectives: one of them is to quantify the dynamic downscaling ability to represent the current climate and the other is to identify critical aspects of the regional climate model in South America in order to interpret the reliability of future projections for the end of the twenty-first century in the A2 scenario of the IPCC Special Report on Emissions Scenarios. In general, the MM5 model is able to reproduce adequately the main general features, seasonal cycle and year-to-year variability of near surface variables over South America. The spatial distribution of temperature is well represented, but some systematic errors were identified, such as an overestimation in central and northern Argentina and an underestimation in the mountainous regions throughout the year. The general structure of precipitation is also well captured by the regional model, although it overestimates the precipitation in the Andean region (specifically in central and southern Chile) in all seasons and underestimates the rainfall over tropical latitudes. The annual cycle of precipitation is adequately represented in the subregions analyzed, but its representation is better over La Plata basin (LPB), Cuyo (CU) and southeastern Pampas (SEP). The annual cycle of mean temperature is well represented, too. The model systematically overestimates the interannual variability of temperature and underestimates the interannual variability of precipitation. From the analyses of interannual and internal variability, as well as the biases, it can be concluded that regardless the season, the simulated precipitation is reliable at subtropical latitudes, Uruguay, southern Brazil and east-central of Argentina, but is less reliable over areas of complex topography. For temperature, the regional model is reliable over subtropical latitudes, Uruguay and the south of Brazil only during winter, but it is less reliable or it is even in the limit of reliability over central and southern Chile all along the year. Therefore, it is concluded that the MM5 model is a useful tool for the generation of regional climate change scenarios and for the evaluation of regional climate change scenarios over southern South America.     

Precision assessment of the orthometric heights determination in northern part of Algeria by combining the GPS data and the local geoid model

The main purpose of this article is to discuss the use of GPS positioning together with a gravimetrically determined geoid, for deriving orthometric heights in the North of Algeria, for which a limited number of GPS stations with known orthometric heights are available, and to check, by the same opportunity, the possibility of substituting the classical spirit levelling. For this work, 247 GPS stations which are homogeneously distributed and collected from the international TYRGEONET project, as well as the local GPS/Levelling surveys, have been used. The GPS/Levelling geoidal heights are obtained by connecting the points to the levelling network while gravimetric geoidal heights were interpolated from the geoid model computed by the Geodetic Laboratory of the National Centre of Spatial Techniques from gravity data supplied by BGI. However, and in order to minimise the discordances, systematic errors and datum inconsistencies between the available height data sets, we have tested two parametric models of corrector surface: a four parameter transformation and a third polynomial model are used to find the adequate functional representation of the correction that should be applied to the gravimetric geoid. The comparisons based on these GPS campaigns prove that a good fit between the geoid model and GPS/levelling data has been reached when the third order polynomial was used as corrector surface and that the orthometric heights can be deducted from GPS observations with an accuracy acceptable for the low order levelling network densification. In addition, the adopted methodology has been also applied for the altimetric auscultation of a storage reservoir situated at 40 km from the town of Oran. The comparison between the computed orthometric heights and observed ones allowed us to affirm that the alternative of levelling by GPS is attractive for this auscultation.     

Recent geoid computations for the Hellenic area

New geoid computations for the Hellenic area are carried out using (a) gravity anomalies for the land area available from old and new data bases, and gravity data for the sea area derived from altimetry and a recent digitization of sea gravity maps, and (b) a 1km × 1km digital terrain model. The EGM96 geopotential model is used as the reference field. In order to assess the quality of the computed geoid heights in the continental area comparisons were carried out with GPS/leveling heights and the recently available European Gravimetric Geoid EGG97. In the sea area the geoid heights were compared with sea surface heights of the recent and more accurate TOPEX/POSEIDON (T/P) altimetry mission. At the end of this article the improvement of the data bases is discussed and some plans for further development in the methodological schedule are pointed out.     

智利北部—阿根廷西北部两类斑岩铜矿成矿方式

智利北部和阿根廷西北部的中新生代斑岩铜矿形成于古生代地体拼贴造山带背景。随着大西洋的张开,南美大陆向西漂移,中新生代期间,南美克拉通块体俯冲到古生代造山带之下形成加厚或双倍地壳。智利北部作为南美活动大陆边缘的组成部分,不断"吞食"向东俯冲的太平洋（纳斯卡）板块,斑岩铜矿成矿作用发生在俯冲板块断离后导致的大规模岩浆活动,并沿再活化岩石圈不连续（先存的古生代拼接带、区域断裂）反复就位,形成安第斯型斑岩铜矿。阿根廷西北部大规模铜（金、钼）成矿与加厚的造山带垮塌有关,大规模成矿受控于造山岩石圈去根、软流圈物质和热上涌引发的大规模岩浆活动。总体而言,智利北部、阿根廷西北部安第斯型和造山带垮塌型斑岩铜矿,乃至南美安第斯山铜（金）矿成矿带形成,与中新生代以来南美大陆向西漂移、大西洋张开事件关系密切。     

EGY-HGM2016: an improved hybrid local geoid model for Egypt based on the combination of GOCE-based geopotential model with gravimetric and GNSS/levelling measurements

An improved hybrid gravimetric geoid model for Egypt, EGY-HGM2016, has been recently computed implementing the least-squares collocation (LSC) method through the remove-compute-restore (RCR) procedure. The computation of EGY-HGM2016 involves different datasets in terms of gravity anomalies determined from the GOCE (gravity field and steady-state ocean circulation explorer)-based global geopotential model (SPW-R4) up to d/o 200 and EGM2008 from d/o 201 to 720 combined with terrestrial gravity datasets in terms of 2140 gravity field anomalies and about 121,480 marine surface gravity anomalies. In addition, orthometric heights from 17 GPS/levelling measurements have been considered during the modelling process to improve the determination of the hybrid gravimetric geoid over the Egyptian region. The EGY-HGM2016 model estimated over Egypt provides geoid heights that are ranging from 7.677 to 21.095 m with a standard deviation (st. dev.) of about 2.534 m in the northwest of the country excluding the involvement of the orthometric heights from GPS/levelling measurements. When the later dataset is considered during the implementation of LSC process, hybrid residual height anomalies ranging from ?1.5 to +0.9 m, with a mean of 0.22 m and a st. dev. of 0.17 m, are obtained. Comparison of the predicted hybrid gravimetric geoid with the corresponding ones obtained from EGM2008, GOCE-based SPW R4 model, and GPS/levelling reveals considerable improvements of our EGY-HGM2016 model over Egypt.     

应用GPS数据与重力数据确定大地水准面

作为GPS/重力边值问题理论及方法的应用,在对GPS/重力方法确定(似)大地水准面的原理进行简要介绍与分析的基础上,利用收集到的N区的600个GPS/重力数据和48个高精度GPS水准数据,计算出该区域的(似)大地水准面。通过拟合法和系统差直接改正法进行的精度分析表明,应用GPS/重力数据结合水准方法确定的该地区(似)大地水准面的精度达到厘米级精度。     

EGMlab, a scientific software for determining the gravity and gradient components from global geopotential models

Nowadays, Global Geopotential Models (GGMs) are used as a routine stage in the procedures to compute a gravimetric geoid. The GGMs based geoidal height also can be used for GPS/levelling and navigation purposes in developing countries which do not have accurate gravimetric geoid models. Also, the GGM based gravity anomaly including the digital elevation model can be used in evaluation and outlier detections schemes of the ground gravity anomaly data. Further, the deflection of vertical and gravity gradients components from the GGMs can be used for different geodetic and geophysical interpretation purposes. However, still a complete and user-friendly software package is not available for universities and geosciences communities. In this article, first we review the procedure for determination of the basic gravity field and gradient components from the GGMs, then general MATLAB based software is presented and applied as a sample case study for determination of gravity components based on the most recent EIGEN-GL04C GRACE model in Sweden. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

On the determination of a new Australian geoid

A new, high precision, high accuracy and high resolution gravimetric geoid of Australia has been produced using most updated data, theory and methodology. This paper presents a concise report of the new Australian geoid determination. Some aspects of the new geoid computation, such as data validation, geoid determination strategies and computational procedures, are described. The relative precision of the new geoid is better than 5 cm for average baseline length of 4km~40km and 18 cm for average baseline length of 120km when compared with three local GPS/levelling networks.     

大地基准的现代化和卫星大地测量新成果--参加国际大地测量协会(IAG)2003年日本札晃大会札记

国际大地测量协会（IAG）2003年 7月在日本札晃举行了近 8天的学术讨论会。现将讨论会中关于我国近邻国家大地基准现代化,低轨卫星（LEO）解算地球重力场的新方法---《一步法》,及GPS数据处理方法的新进展等进行扼要介绍。《国际地面参考框架（ITRF）》一直作为全球的大地测量坐标框架,它目前存在的主要问题是它的框架点是一律被假定作常速运动的,但实际情况与ITRF的这一假定是不一致的。在这次大会上我国的近邻国家如日本、蒙古、新西兰、韩国、马来西亚等报告了本国大地坐标系统、坐标框架的更新和现代化,以及地区大地水准面的精化工作。通常利用二步法从低轨卫星的轨道信息求定地球重力场及其时变,大会上报告的一步法是对这一经典方法的重大改进和突破。GPS数据处理方法的新进展中比较突出的有：顾及观测数据短暂性相关影响的GPS连续运行网坐标及其移动速度的解算方案;含有系统误差的GPS观测数据的解算等。     

Factors affecting FFT gravimetric geoid determination precision

The Fast Fourier Transformation (FFT) has become a routine mathematical tool for the refinement of the Earth's gravity field, such as the computation of precise gravimetric geoid and terrain corrections, particularly over a large area. This paper presents ideas and methodologies to evaluate the accuracy of geoid undulation computations using FFT. A global geopotential model is used as a ‘ground truth’ gravity field model to assess the geoid determination precision by using FFT technique. It is demonstrated that special considerations must be given for a high precision FFT gravimetric geoid determination. A maximum of a few decimetres error could be introduced by the FFT algorithm if the gravity anomalies are not long wavelength filtered and/or no zero padding is applied.     

First record of stereospondyls (Tetrapoda,Temnospondyli) in the Upper Triassic of Southern Brazil

Stereospondyls survived the Permo-Triassic extinctions in a refuge probably located in the landmass that nowadays comprises Australia. Subsequently, they radiated to other parts of Pangaea, reaching their highest distribution and diversification during the Early Triassic. An incomplete interclavicle from the Caturrita Formation represents their first record in the Upper Triassic of Brazil. Previously, Upper Triassic South American stereospondyls were restricted to Argentina. This new record reinforces a former hypothesis that suggests the presence of a more diverse stereospondyl fauna in South America during the Late Triassic than previously assumed. Additionally, the presence of a stereospondyl and a phytosaur in the Caturrita Formation reinforces the hypothesis of a change to more humid climatic conditions in the Paraná Basin during the Upper Triassic. The record of Early Triassic stereospondyls in South America suggests that they first colonized Brazil and/or Uruguay, spreading from South Africa during the Early Triassic, subsequently reaching Argentina. Up till now, there is no record of Middle Triassic stereospondyls in either Argentina and Brazil, probably due to either taphonomic bias or insufficient prospecting. Despite the lack of direct evidence, one should not dismiss an earlier stereospondyl colonization of Argentina still during the Early or Middle Triassic.     

关于重力勘查的高度改正应采用何种高程系统的讨论

重力勘查的观测值必须进行正常场与高度改正,为此要测量重力测点的平面位置与高程,尤其对高程测量精度要求甚严。以前测地工作使用经纬仪、水准仪和测距仪等,高程系统按“规范”(1)采用大地水准面的“正高系统”.近来全球卫星定位系统(GPS)开始应用于重力勘查中,GPS测量的高程为WGS-84的“大地高系统”,这样需要把“大地高”换算到“正高”.这样做即要花费许多工作量,又会增加换算中产生的误差。本文通过对地球正常重力场理论公式的分析,认为可以采用GPS测量的大地高进行重力高度改正。     

南美洲铁矿重要成矿区带与成矿规律研究

南美洲铁矿资源丰富、分布广泛、矿床规模大、矿石品位高、类型较齐全。从南美洲大地构造格架分析,分布在克拉通地盾区或其边缘的沉积-变质型铁矿和分布在安第斯海岸造山带的火山成因型铁矿最为重要,沉积-变质型铁矿主要形成于太古宙—元古宙,火山成因型铁矿主要形成于中新生代。综合研究分析,圈出10个铁矿重要成矿区带,主要分布在委内瑞拉的玻利瓦尔州,巴西的阿马帕州、帕拉州、米纳斯吉拉斯州、北里奥格兰德州、巴拉伊巴州、巴伊亚州、南马托格罗索州,玻利维亚的圣克鲁斯州,秘鲁的伊卡—阿雷基帕海岸带和智利的卡拉马—圣地亚哥海岸带等。通过对重要成矿区带内典型矿床的剖析,并采用统计分析方法总结南美洲铁矿成矿时空分布规律,为未来在该区寻找更多的铁矿资源提供参考。     

大地测量在研究青藏地壳运动及其机制中的作用

基于对GPS和重力测量数据的分析,扼要地介绍了重力与内部构造关系,大地水准面的场源与其动力效应等方面的研究内容和方法。总结了利用重力与形变资料研究青藏地区的现代地壳运动及其动力机制的研究结果。