R. Pail, J. Bamber, R. Biancale, R. Bingham, C. Braitenberg, A. Eicker, F. Flechtner, T. Gruber, A. Güntner, G. Heinzel, M. Horwath, L. Longuevergne, J. Müller, I. Panet, H. Savenije, S. Seneviratne, N. Sneeuw, T. V. van Dam, B. Wouters
Abstract As changes in gravity are directly related to mass variability, satellite missions observing the Earth’s time varying gravity field are a unique tool for observing mass transport processes in the Earth system, such as the water cycle, rapid changes in the cryosphere, oceans, and solid Earth processes, on a global scale. The observation of Earth’s gravity field was successfully performed by the GRACE and GOCE satellite missions, and will be continued by the GRACE Follow-On mission. A comprehensive team of European scientists proposed the next-generation gravity field mission MOBILE in response to the European Space Agency (ESA) call for a Core Mission in the frame of Earth Explorer 10 (EE10). MOBILE is based on the innovative observational concept of a high-low tracking formation with micrometer ranging accuracy, complemented by new instrument concepts. Since a high-low tracking mission primarily observes the radial component of gravity-induced orbit perturbations, the error structure is close to isotropic. This geometry significantly reduces artefacts of previous along-track ranging low-low formations (GRACE, GRACE-Follow-On) such as the typical striping patterns. The minimum configuration consists of at least two medium-Earth orbiters (MEOs) at 10000 km altitude or higher, and one low-Earth orbiter (LEO) at 350-400 km. The main instrument is a laser-based distance or distance change measurement system, which is placed at the LEO. The MEOs are equipped either with passive reflectors or transponders. In a numerical closed-loop simulation, it was demonstrated that this minimum configuration is in agreement with the threshold science requirements of 5 mm equivalent water height (EWH) accuracy at 400 km wavelength, and 10 cm EWH at 200 km. MOBILE provides promising potential future perspectives by linking the concept to existing space infrastructure such as Galileo next-generation, as future element of the Copernicus/Sentinel programme, and holds the potential of miniaturization even up to swarm configurations. As such MOBILE can be considered as a precursor and role model for a sustained mass transport observing system from space.
{"title":"Mass variation observing system by high low inter-satellite links (MOBILE) – a new concept for sustained observation of mass transport from space","authors":"R. Pail, J. Bamber, R. Biancale, R. Bingham, C. Braitenberg, A. Eicker, F. Flechtner, T. Gruber, A. Güntner, G. Heinzel, M. Horwath, L. Longuevergne, J. Müller, I. Panet, H. Savenije, S. Seneviratne, N. Sneeuw, T. V. van Dam, B. Wouters","doi":"10.1515/jogs-2019-0006","DOIUrl":"https://doi.org/10.1515/jogs-2019-0006","url":null,"abstract":"Abstract As changes in gravity are directly related to mass variability, satellite missions observing the Earth’s time varying gravity field are a unique tool for observing mass transport processes in the Earth system, such as the water cycle, rapid changes in the cryosphere, oceans, and solid Earth processes, on a global scale. The observation of Earth’s gravity field was successfully performed by the GRACE and GOCE satellite missions, and will be continued by the GRACE Follow-On mission. A comprehensive team of European scientists proposed the next-generation gravity field mission MOBILE in response to the European Space Agency (ESA) call for a Core Mission in the frame of Earth Explorer 10 (EE10). MOBILE is based on the innovative observational concept of a high-low tracking formation with micrometer ranging accuracy, complemented by new instrument concepts. Since a high-low tracking mission primarily observes the radial component of gravity-induced orbit perturbations, the error structure is close to isotropic. This geometry significantly reduces artefacts of previous along-track ranging low-low formations (GRACE, GRACE-Follow-On) such as the typical striping patterns. The minimum configuration consists of at least two medium-Earth orbiters (MEOs) at 10000 km altitude or higher, and one low-Earth orbiter (LEO) at 350-400 km. The main instrument is a laser-based distance or distance change measurement system, which is placed at the LEO. The MEOs are equipped either with passive reflectors or transponders. In a numerical closed-loop simulation, it was demonstrated that this minimum configuration is in agreement with the threshold science requirements of 5 mm equivalent water height (EWH) accuracy at 400 km wavelength, and 10 cm EWH at 200 km. MOBILE provides promising potential future perspectives by linking the concept to existing space infrastructure such as Galileo next-generation, as future element of the Copernicus/Sentinel programme, and holds the potential of miniaturization even up to swarm configurations. As such MOBILE can be considered as a precursor and role model for a sustained mass transport observing system from space.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":"4 1","pages":"48 - 58"},"PeriodicalIF":1.3,"publicationDate":"2019-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"75830225","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract The Central Andes experienced major earthquake (Mw =8.2) in April 2014 in a region where the giant 1877 earthquake (Mw=8.8) occurred. The 2014 Iquique earthquake did not break the entire seismic gap zones as previously predicted. Geodetic and seismological observations indicate a highly coupled plate interface. To assess the locking mechanism of plate interfaces beneath Central Andes, a 2.5-D gravity model of the crust and upper mantle structure of the central segment of the subduction zone was developed based on terrestrial and satellite gravity data from the LAGEOS, GRACE and GOCE satellite missions. The densities and major structures of the gravity model are constrained by velocity models from receiver function and seismic tomography. The gravity model defined details of crustal and slab structure necessary to understand the cause of megathrust asperity generation. The densities of the upper and lower crust in the fore-arc (2970 – 3000 kg m−3) are much higher than the average density of continental crust. The high density bodies are interpreted as plutonic or ophiolitic structures emplaced onto continental crust. The plutonic or ophiolitic structures may be exerting pressure on the Nazca slab and lock the plate interfaces beneath the Central Andes subduction zone. Thus, normal pressure exerted by high density fore-arc structures and buoyancy force may control plate coupling in the Central Andes. However, this interpretation does not exclude other possible factors controlling plate coupling in the Central Andes. Seafloor roughness and variations in pore-fluid pressure in sediments along subduction channel can affect plate coupling and asperity generation.
2014年4月,在1877年发生8.8级大地震的地区,安第斯山脉中部发生了8.2级大地震。2014年的伊基克地震并没有像之前预测的那样破坏整个地震间隙带。大地测量和地震观测显示了一个高度耦合的板块界面。为了评估安第斯中部板块界面的锁锁机制,基于LAGEOS、GRACE和GOCE卫星任务的地面和卫星重力数据,建立了俯冲带中段地壳和上地幔结构的2.5维重力模型。重力模型的密度和主要结构受接收函数和地震层析成像的速度模型的约束。重力模型定义了地壳和板块结构的细节,这对于理解巨型逆冲隆起的成因是必要的。弧前上下地壳密度(2970 ~ 3000 kg m−3)远高于大陆地壳的平均密度。高密度体被解释为侵位在大陆地壳上的深成或蛇绿构造。深成或蛇绿构造可能对纳斯卡板块施加压力,并锁住了中安第斯俯冲带下的板块界面。因此,高密度弧前构造施加的正压力和浮力可能控制安第斯山脉中部的板块耦合。然而,这种解释并不排除控制安第斯山脉中部板块耦合的其他可能因素。海底粗糙度和俯冲通道沉积物孔隙流体压力的变化可以影响板块耦合和凹凸不平的产生。
{"title":"Plate Coupling Mechanism of the Central Andes Subduction: Insight from Gravity Model","authors":"R. Mahatsente","doi":"10.1515/jogs-2019-0002","DOIUrl":"https://doi.org/10.1515/jogs-2019-0002","url":null,"abstract":"Abstract The Central Andes experienced major earthquake (Mw =8.2) in April 2014 in a region where the giant 1877 earthquake (Mw=8.8) occurred. The 2014 Iquique earthquake did not break the entire seismic gap zones as previously predicted. Geodetic and seismological observations indicate a highly coupled plate interface. To assess the locking mechanism of plate interfaces beneath Central Andes, a 2.5-D gravity model of the crust and upper mantle structure of the central segment of the subduction zone was developed based on terrestrial and satellite gravity data from the LAGEOS, GRACE and GOCE satellite missions. The densities and major structures of the gravity model are constrained by velocity models from receiver function and seismic tomography. The gravity model defined details of crustal and slab structure necessary to understand the cause of megathrust asperity generation. The densities of the upper and lower crust in the fore-arc (2970 – 3000 kg m−3) are much higher than the average density of continental crust. The high density bodies are interpreted as plutonic or ophiolitic structures emplaced onto continental crust. The plutonic or ophiolitic structures may be exerting pressure on the Nazca slab and lock the plate interfaces beneath the Central Andes subduction zone. Thus, normal pressure exerted by high density fore-arc structures and buoyancy force may control plate coupling in the Central Andes. However, this interpretation does not exclude other possible factors controlling plate coupling in the Central Andes. Seafloor roughness and variations in pore-fluid pressure in sediments along subduction channel can affect plate coupling and asperity generation.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":"7 1","pages":"13 - 21"},"PeriodicalIF":1.3,"publicationDate":"2019-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78638272","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract The paper presents a new approach to the airborne vector gravimetry problem. The idea of the approach is to take into account spatial correlation of the gravity field to improve observability of horizontal components of the gravity disturbance vector (GDV). We consider the GDV determination problem given airborne data at a set of parallel survey lines assuming that lines are flown in the same direction at a constant height above the reference ellipsoid. We use a 2-D random field model for the gravity field at the flight height. The random field is governed by two autoregressive equations (one in the direction along the lines, the other across the lines). Then we pose the estimation problem simultaneously for the GDV horizontal components and systematic errors of an inertial navigation system at all the lines simultaneously. The developed estimation algorithm is based on 2D Kalman filtering and smoothing techniques. Numerical results obtained from simulated data processing showed improved accuracy of the gravity horizontal component determination.
{"title":"Two-dimensional Kalman filter approach to airborne vector gravimetry","authors":"V. Vyazmin, Y. Bolotin","doi":"10.1515/jogs-2019-0009","DOIUrl":"https://doi.org/10.1515/jogs-2019-0009","url":null,"abstract":"Abstract The paper presents a new approach to the airborne vector gravimetry problem. The idea of the approach is to take into account spatial correlation of the gravity field to improve observability of horizontal components of the gravity disturbance vector (GDV). We consider the GDV determination problem given airborne data at a set of parallel survey lines assuming that lines are flown in the same direction at a constant height above the reference ellipsoid. We use a 2-D random field model for the gravity field at the flight height. The random field is governed by two autoregressive equations (one in the direction along the lines, the other across the lines). Then we pose the estimation problem simultaneously for the GDV horizontal components and systematic errors of an inertial navigation system at all the lines simultaneously. The developed estimation algorithm is based on 2D Kalman filtering and smoothing techniques. Numerical results obtained from simulated data processing showed improved accuracy of the gravity horizontal component determination.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":"21 1","pages":"87 - 96"},"PeriodicalIF":1.3,"publicationDate":"2019-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90983053","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract Atmospheric water vapour, a major component in weather systems serves as the main source for precipitation, provides latent heat which helps maintain the earth’s energy balance and a major parameter in Numerical Weather Prediction (NWP) models. An observational technique based on the Global Navigation Satellite System (GNSS) has made it possible to easily retrieve Precipitable Water (PW) at station’s antenna position with very high spatial and temporal variabilities. GNSS techniques are superior to ground-based and balloons sensors in terms of accuracy, ease of use, wider coverage and easier assimilation into NWP models. This study sought to use prediction models using daily observational data from Four (4) International GNSS Service stations in West Africa. The best prediction model can be used in cases of station outages and to predict PW over data poor regions using computed Zenith Tropospheric Delays (ZTD). gLAB software was used to process the stations’ data in Precise Point Positioning mode and PW were retrieved using station’s temperature and pressure values. Computed PW were compared against Total Column Water Vapour from ERA-Interim Reanalysis data in 2016. Correlation coefficient (R2) values ranging from 0.947 — 0.995 were obtained for the four stations. With computed PW’s, three regression models were tested to find the best-fit with PW as the dependent variable and ZTD being the independent variable. The quadratic model gave the highest R2 and lowest RMSE values as against the linear and exponential models. Time series forecasts models such as moving average, autoregressive, exponential smoothing and autoregressive integrated moving average were also employed. The forecasts results were compared against ZTD with autoregressive model reporting the highest R2 and lowest RMSE amongst the forecast models developed.
{"title":"Application of GNSS derived precipitable water vapour prediction in West Africa","authors":"Akwasi Acheampong, K. Obeng","doi":"10.1515/jogs-2019-0005","DOIUrl":"https://doi.org/10.1515/jogs-2019-0005","url":null,"abstract":"Abstract Atmospheric water vapour, a major component in weather systems serves as the main source for precipitation, provides latent heat which helps maintain the earth’s energy balance and a major parameter in Numerical Weather Prediction (NWP) models. An observational technique based on the Global Navigation Satellite System (GNSS) has made it possible to easily retrieve Precipitable Water (PW) at station’s antenna position with very high spatial and temporal variabilities. GNSS techniques are superior to ground-based and balloons sensors in terms of accuracy, ease of use, wider coverage and easier assimilation into NWP models. This study sought to use prediction models using daily observational data from Four (4) International GNSS Service stations in West Africa. The best prediction model can be used in cases of station outages and to predict PW over data poor regions using computed Zenith Tropospheric Delays (ZTD). gLAB software was used to process the stations’ data in Precise Point Positioning mode and PW were retrieved using station’s temperature and pressure values. Computed PW were compared against Total Column Water Vapour from ERA-Interim Reanalysis data in 2016. Correlation coefficient (R2) values ranging from 0.947 — 0.995 were obtained for the four stations. With computed PW’s, three regression models were tested to find the best-fit with PW as the dependent variable and ZTD being the independent variable. The quadratic model gave the highest R2 and lowest RMSE values as against the linear and exponential models. Time series forecasts models such as moving average, autoregressive, exponential smoothing and autoregressive integrated moving average were also employed. The forecasts results were compared against ZTD with autoregressive model reporting the highest R2 and lowest RMSE amongst the forecast models developed.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":"19 1","pages":"41 - 47"},"PeriodicalIF":1.3,"publicationDate":"2019-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91296971","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract The topographic potential bias at geoid level is the error of the analytically continued geopotential from or above the Earth’s surface to the geoid. We show that the topographic potential can be expressed as the sum of two Bouguer shell components, where the density distribution of one is spherical symmetric and the other is harmonic at any point along the normal to a sphere through the computation point. As a harmonic potential does not affect the bias, the resulting topographic bias is that of the first component, i.e. the spherical symmetric Bouguer shell. This implies that the so-called terrain potential is not likely to contribute significantly to the bias. We present three examples of the geoid bias for different topographic density distributions.
{"title":"The topographic bias in gravimetric geoid determination revisited","authors":"L. Sjöberg","doi":"10.1515/jogs-2019-0007","DOIUrl":"https://doi.org/10.1515/jogs-2019-0007","url":null,"abstract":"Abstract The topographic potential bias at geoid level is the error of the analytically continued geopotential from or above the Earth’s surface to the geoid. We show that the topographic potential can be expressed as the sum of two Bouguer shell components, where the density distribution of one is spherical symmetric and the other is harmonic at any point along the normal to a sphere through the computation point. As a harmonic potential does not affect the bias, the resulting topographic bias is that of the first component, i.e. the spherical symmetric Bouguer shell. This implies that the so-called terrain potential is not likely to contribute significantly to the bias. We present three examples of the geoid bias for different topographic density distributions.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":"9 1","pages":"59 - 64"},"PeriodicalIF":1.3,"publicationDate":"2019-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"77734559","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract The signal content and error level of recent GOCE-based high resolution gravity field models is assessed by means of signal degree variances and comparisons to independent GNSS-levelling geoid heights. The signal of the spherical harmonic series of these models is compared to the pre-GOCE EGM2008 model in order to identify the impact of GOCE data, of improved surface and altimetric gravity data and of modelling approaches. Results of the signal analysis show that in a global average roughly 80% of the differences are due to the inclusion of GOCE satellite information, while the remaining 20% are contributed by improved surface data. Comparisons of the global models to GNSS-levelling derived geoid heights demonstrate that a 1 cm geoid from the global model is feasible, if there is a high quality terrestrial gravity data set available. For areas with less good coverage an accuracy of several centimetres to a decimetre is feasible taking into account that GOCE provides now the geoid with a centimetre accuracy at spatial scales of 80 to 100 km. Comparisons with GNSS-levelling geoid heights also are a good tool to investigate possible systematic errors in the global models, in the spirit levelling and in the GNSS height observations. By means of geoid height differences and geoid slope differences one can draw conclusions for each regional data set separately. These conclusions need to be considered for a refined analysis e.g. to eliminate suspicious GNSS-levelling data, to improve the global modelling by using full variance-covariance matrices and by consistently weighting the various data sources used for high resolution gravity field models. The paper describes the applied procedures, shows results for these geoid height and geoid slope differences for some regional data sets and draws conclusions about possible error sources and future work to be done in this context.
{"title":"Signal and error assessment of GOCE-based high resolution gravity field models","authors":"T. Gruber, M. Willberg","doi":"10.1515/jogs-2019-0008","DOIUrl":"https://doi.org/10.1515/jogs-2019-0008","url":null,"abstract":"Abstract The signal content and error level of recent GOCE-based high resolution gravity field models is assessed by means of signal degree variances and comparisons to independent GNSS-levelling geoid heights. The signal of the spherical harmonic series of these models is compared to the pre-GOCE EGM2008 model in order to identify the impact of GOCE data, of improved surface and altimetric gravity data and of modelling approaches. Results of the signal analysis show that in a global average roughly 80% of the differences are due to the inclusion of GOCE satellite information, while the remaining 20% are contributed by improved surface data. Comparisons of the global models to GNSS-levelling derived geoid heights demonstrate that a 1 cm geoid from the global model is feasible, if there is a high quality terrestrial gravity data set available. For areas with less good coverage an accuracy of several centimetres to a decimetre is feasible taking into account that GOCE provides now the geoid with a centimetre accuracy at spatial scales of 80 to 100 km. Comparisons with GNSS-levelling geoid heights also are a good tool to investigate possible systematic errors in the global models, in the spirit levelling and in the GNSS height observations. By means of geoid height differences and geoid slope differences one can draw conclusions for each regional data set separately. These conclusions need to be considered for a refined analysis e.g. to eliminate suspicious GNSS-levelling data, to improve the global modelling by using full variance-covariance matrices and by consistently weighting the various data sources used for high resolution gravity field models. The paper describes the applied procedures, shows results for these geoid height and geoid slope differences for some regional data sets and draws conclusions about possible error sources and future work to be done in this context.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":"58 2 1","pages":"71 - 86"},"PeriodicalIF":1.3,"publicationDate":"2019-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"77754139","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract The present study computes B-W extension model (extended Bursa-Wolf model) coordinate transformation parameters from World Geodetic System 1984 (WGS-84) to the Everest datum namely Everest (1830) and Everest (1956) using records of coordinate measurements from Global Positioning System (GPS) observable across Nepal region. Synthetic or modeled coordinates were determined by using the Artificial Neural Network (ANN) and Singular Value Decomposition (SVD) methods. We studied 9-transformation parameters with the help of the ANN technique and validated the outcomes with the SVD method. The comparative analysis of the ANN, as well as SVD methods, was done with the observed output following one way ANOVA test. The analysis showed that the null hypothesis for both datums were acceptable and suggesting all models statistically significantly equivalent to each other. The outcomes from this study would complement a relatively better understanding of the techniques for coordinate transformation and precise coordinate assignment while assimilating data sets from different resources.
{"title":"Coordinate transformation parameters in Nepal by using neural network and SVD methods","authors":"K. Ansari, P. Gyawali, P. Pradhan, Kwan-Dong Park","doi":"10.1515/jogs-2019-0003","DOIUrl":"https://doi.org/10.1515/jogs-2019-0003","url":null,"abstract":"Abstract The present study computes B-W extension model (extended Bursa-Wolf model) coordinate transformation parameters from World Geodetic System 1984 (WGS-84) to the Everest datum namely Everest (1830) and Everest (1956) using records of coordinate measurements from Global Positioning System (GPS) observable across Nepal region. Synthetic or modeled coordinates were determined by using the Artificial Neural Network (ANN) and Singular Value Decomposition (SVD) methods. We studied 9-transformation parameters with the help of the ANN technique and validated the outcomes with the SVD method. The comparative analysis of the ANN, as well as SVD methods, was done with the observed output following one way ANOVA test. The analysis showed that the null hypothesis for both datums were acceptable and suggesting all models statistically significantly equivalent to each other. The outcomes from this study would complement a relatively better understanding of the techniques for coordinate transformation and precise coordinate assignment while assimilating data sets from different resources.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":"28 4 1","pages":"22 - 28"},"PeriodicalIF":1.3,"publicationDate":"2019-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87048367","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract For the calculation of gravity disturbance in the Earth’s external gravity field, the Stokes-Pizzetti integral is a commonly used method. However, when the target point approaches the Earth’s surface, such problems as singularity and discontinuity arise due to the Stokes kernel structure itself. To settle the problems, firstly the reason for singularity and discontinuity was discussed, and then modification was made to the integral formula, by which the singularity at the surface point is eliminated. Finally the non-singular integral formulas for the calculation of disturbing gravity were derived. In numerical experiments, an area in China was selected to test the modified formula. Numerical results show that the modified formula performs much better than classical Stokes-Pizzetti integral formula when dealing with the calculation of the radial component of gravity disturbance near the Earth’s surface.
{"title":"The modified integral method for the determination of gravity disturbance near the Earth’s surface","authors":"D. Zhao, S. Li, Q. Wang, Z. Gong","doi":"10.1515/jogs-2019-0016","DOIUrl":"https://doi.org/10.1515/jogs-2019-0016","url":null,"abstract":"Abstract For the calculation of gravity disturbance in the Earth’s external gravity field, the Stokes-Pizzetti integral is a commonly used method. However, when the target point approaches the Earth’s surface, such problems as singularity and discontinuity arise due to the Stokes kernel structure itself. To settle the problems, firstly the reason for singularity and discontinuity was discussed, and then modification was made to the integral formula, by which the singularity at the surface point is eliminated. Finally the non-singular integral formulas for the calculation of disturbing gravity were derived. In numerical experiments, an area in China was selected to test the modified formula. Numerical results show that the modified formula performs much better than classical Stokes-Pizzetti integral formula when dealing with the calculation of the radial component of gravity disturbance near the Earth’s surface.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":"9 1","pages":"65 - 70"},"PeriodicalIF":1.3,"publicationDate":"2019-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90313111","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract For the modelling and determination of the Earth’s external gravity potential as well as its second-order radial derivatives in the space near sea surface, the surface layer integral method was discussed in the paper. The reasons for the applicability of the method over sea surface were discussed. From the original integral formula of disturbing potential based on the surface layer method, the expression of the radial component of the gravity gradient tensor was derived. Furthermore, an identity relation was introduced to modify the formula in order to reduce the singularity problem. Numerical experiments carried out over the marine area of China show that, the modi-fied surface layer integral method effectively improves the accuracy and reliability of the calculation of the second-order radial gradient component of the disturbing potential near sea surface.
{"title":"The surface layer integral method for the modelling of the radial gravity gradient component over sea surface","authors":"D. Zhao, Z. Gong, J. Feng","doi":"10.1515/jogs-2019-0012","DOIUrl":"https://doi.org/10.1515/jogs-2019-0012","url":null,"abstract":"Abstract For the modelling and determination of the Earth’s external gravity potential as well as its second-order radial derivatives in the space near sea surface, the surface layer integral method was discussed in the paper. The reasons for the applicability of the method over sea surface were discussed. From the original integral formula of disturbing potential based on the surface layer method, the expression of the radial component of the gravity gradient tensor was derived. Furthermore, an identity relation was introduced to modify the formula in order to reduce the singularity problem. Numerical experiments carried out over the marine area of China show that, the modi-fied surface layer integral method effectively improves the accuracy and reliability of the calculation of the second-order radial gradient component of the disturbing potential near sea surface.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":"52 1","pages":"127 - 132"},"PeriodicalIF":1.3,"publicationDate":"2019-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"80254860","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract The AntGrav project, funded by the German Research Foundation (DFG) has the main objective to homogenize and optimize Antarctic gravity field information. Within this project an evaluation procedure is needed to inspect all different kind of gravity field surveys available in Antarctica. In this paper a suitable methodology is proposed. We present an approach for fast 3D gravity point data reduction in different spectral bands. This is achieved through pre-calculating a fine 3D mesh of synthesized gravity functionals over the entirety of the Antarctic continent, for which two different global models are used: the combined satellite model GOCO05s for the long-wavelength part, and the topographic model Earth2014 for the shorter wavelengths. To maximize the applicability separate meshes are calculated for different spectral bands in order to specifically reduce a certain band or a selected combination. All meshes are calculated for gravity anomalies as well as gravity disturbances. Utilizing these meshes, synthesized gravity data at arbitrary positions is computed by conventional 3D interpolation methods (e.g. linear, cubic or spline). It is shown that the applied approach can reach a worst-case interpolation error of less than 1 mGal. Evaluation results are presented for the AntGG grid and exemplary for the in-situ measurements of the AGAP and BAS-LAND campaigns. While general properties, large-scale errors and systematic effects can usually be detected, small-scale errors (e.g. of single points) are mostly untraceable due to the uncertainties within the topographic model.
{"title":"Evaluation of terrestrial and airborne gravity data over Antarctica – a generic approach","authors":"P. Zingerle, R. Pail, M. Scheinert, T. Schaller","doi":"10.1515/jogs-2019-0004","DOIUrl":"https://doi.org/10.1515/jogs-2019-0004","url":null,"abstract":"Abstract The AntGrav project, funded by the German Research Foundation (DFG) has the main objective to homogenize and optimize Antarctic gravity field information. Within this project an evaluation procedure is needed to inspect all different kind of gravity field surveys available in Antarctica. In this paper a suitable methodology is proposed. We present an approach for fast 3D gravity point data reduction in different spectral bands. This is achieved through pre-calculating a fine 3D mesh of synthesized gravity functionals over the entirety of the Antarctic continent, for which two different global models are used: the combined satellite model GOCO05s for the long-wavelength part, and the topographic model Earth2014 for the shorter wavelengths. To maximize the applicability separate meshes are calculated for different spectral bands in order to specifically reduce a certain band or a selected combination. All meshes are calculated for gravity anomalies as well as gravity disturbances. Utilizing these meshes, synthesized gravity data at arbitrary positions is computed by conventional 3D interpolation methods (e.g. linear, cubic or spline). It is shown that the applied approach can reach a worst-case interpolation error of less than 1 mGal. Evaluation results are presented for the AntGG grid and exemplary for the in-situ measurements of the AGAP and BAS-LAND campaigns. While general properties, large-scale errors and systematic effects can usually be detected, small-scale errors (e.g. of single points) are mostly untraceable due to the uncertainties within the topographic model.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":"2 1","pages":"29 - 40"},"PeriodicalIF":1.3,"publicationDate":"2019-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"79201612","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}