{"title":"Optimal mathematical and statistical models to estimate vertical crustal movements using satellite altimetry and tide gauge data","authors":"H. Iz, T. Y. Yang, C. Shum, C. Kuo","doi":"10.1515/jogs-2019-0014","DOIUrl":"https://doi.org/10.1515/jogs-2019-0014","url":null,"abstract":"Abstract Knowledge of vertical crustal movement is fundamental to quantify absolute sea level changes at tide gauge locations as well as for satellite altimetry calibration validations. While GPS measurements at collocated tide gauge stations fulfill this need, currently only few hundred tide gauge stations are equipped with GPS, and their measurements do not span a long period of time. In the past, several studies addressed this problem by calculating relative and geocentric trends from the tide gauge and satellite altimetry measurements respectively, and then difference the two trends to calculate the rate of changes at the tide gauge stations. However, this approach is suboptimal. This study offers an optimal statistical protocol based on the method of condition equations with unknown parameters. An example solution demonstrates the proposed mathematical and statistical models’ optimality in estimating vertical crustal movement and its standard error by comparing them with the results of current methods. The proposed model accounts for the effect of autocorrelations in observed tide gauge and satellite altimetry sea level time series, adjusts observed corrections such as inverted barometer effects, and constraints tide gauge and satellite altimeter measurement to close. The new model can accommodate estimating other systematic effects such as pole tides that are not eliminated by differencing.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2019-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"84611530","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}
{"title":"Evaluation of Gravity Data Derived from Global Gravity Field Models Using Terrestrial Gravity Data in Enugu State, Nigeria","authors":"O. I. Apeh, E. C. Moka, V. Uzodinma","doi":"10.1515/jogs-2018-0015","DOIUrl":"https://doi.org/10.1515/jogs-2018-0015","url":null,"abstract":"Abstract Spherical harmonic expansion is a commonly applied mathematical representation of the earth’s gravity field. This representation is implied by the potential coeffcients determined by using elements/parameters of the field observed on the surface of the earth and/or in space outside the earth in the spherical harmonic expansion of the field. International Centre for Gravity Earth Models (ICGEM) publishes, from time to time, Global Gravity Field Models (GGMs) that have been developed. These GGMs need evaluation with terrestrial data of different locations to ascertain their accuracy for application in those locations. In this study, Bouguer gravity anomalies derived from a total of eleven (11) recent GGMs, using sixty sample points, were evaluated by means of Root-Mean-Square difference and correlation coeficient. The Root-Mean-Square differences of the computed Bouguer anomalies from ICGEMwebsite compared to their positionally corresponding terrestrial Bouguer anomalies range from 9.530mgal to 37.113mgal. Additionally, the correlation coe_cients of the structure of the signal of the terrestrial and GGM-derived Bouguer anomalies range from 0.480 to 0.879. It was observed that GECO derived Bouguer gravity anomalies have the best signal structure relationship with the terrestrial data than the other ten GGMs. We also discovered that EIGEN-6C4 and GECO derived Bouguer anomalies have enormous potential to be used as supplements to the terrestrial Bouguer anomalies for Enugu State, Nigeria.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2018-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"85148800","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}
{"title":"Optimal biased Kriging: Homeogram tapering and applications to geoid undulations in Korea","authors":"B. Schaffrin, T. Bae, Y. Felus","doi":"10.1515/jogs-2018-0016","DOIUrl":"https://doi.org/10.1515/jogs-2018-0016","url":null,"abstract":"Abstract This article studies the Optimal Biased Kriging (OBK) approach which is an alternative geostatistical method that gives up the unbiasedness condition of Ordinary Kriging (OK) to gain an improved Mean Squared Prediction Error (MSPE). The system of equations for the optimal linear biased Kriging predictor is derived and itsMSPE is compared with that of Ordinary Kriging. A major impediment in implementing this system of equations and performing Kriging interpolation with massive datasets is the inversion of the spatial coherency matrix. This problem is investigated and a novel method, called “homeogram tapering”, which exploits spatial sorting techniques to create sparse matrices for efficient matrix inversion, is described. Finally, as an application, results from experiments performed on a geoid undulation dataset from Korea are presented. A precise geoid is usually the indispensable basis for meaningful hydrological studies over wide areas. These experiments use the theory presented here along with a relatively new spatial coherency measure, called the homeogram, also known as the non-centered covariance function.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2018-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"85697120","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}
{"title":"Why and how to predict sea level changes at a tide gauge station with prediction intervals","authors":"H. Iz","doi":"10.1515/jogs-2018-0012","DOIUrl":"https://doi.org/10.1515/jogs-2018-0012","url":null,"abstract":"Abstract Predicting sea level rise is essential for current climate discussions. Empirical models put in use to monitor and analyze sea level variations observed at globally distributed tide gauge stations during the last decade can provide reliable predictions with high resolution. Meanwhile, prediction intervals, an alternative to confidence intervals, are to be recognized and deployed in sea level studies. Predictions together with their prediction intervals, as demonstrated in this study, can quantify the uncertainty of a single future observation from a population, instead of the uncertainty of a conceivable average sea level namely a confidence interval, and it is thereby, better suited for coastal risk assessment to guide policy development for mitigation and adaptation responses.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2018-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"82241679","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}
{"title":"On the geoid and orthometric height vs. quasigeoid and normal height","authors":"L. Sjöberg","doi":"10.1515/jogs-2018-0011","DOIUrl":"https://doi.org/10.1515/jogs-2018-0011","url":null,"abstract":"Abstract The geoid, but not the quasigeoid, is an equipotential surface in the Earth’s gravity field that can serve both as a geodetic datum and a reference surface in geophysics. It is also a natural zero-level surface, as it agrees with the undisturbed mean sea level. Orthometric heights are physical heights above the geoid,while normal heights are geometric heights (of the telluroid) above the reference ellipsoid. Normal heights and the quasigeoid can be determined without any information on the Earth’s topographic density distribution, which is not the case for orthometric heights and geoid. We show from various derivations that the difference between the geoid and the quasigeoid heights, being of the order of 5 m, can be expressed by the simple Bouguer gravity anomaly as the only term that includes the topographic density distribution. This implies that recent formulas, including the refined Bouguer anomaly and a difference between topographic gravity potentials, do not necessarily improve the result. Intuitively one may assume that the quasigeoid, closely related with the Earth’s surface, is rougher than the geoid. For numerical studies the topography is usually divided into blocks of mean elevations, excluding the problem with a non-star shaped Earth. In this case the smoothness of both types of geoid models are affected by the slope of the terrain,which shows that even at high resolutions with ultra-small blocks the geoid model is likely as rough as the quasigeoid model. In case of the real Earth there are areas where the quasigeoid, but not the geoid, is ambiguous, and this problem increases with the numerical resolution of the requested solution. These ambiguities affect also normal and orthometric heights. However, this problem can be solved by using the mean quasigeoid model defined by using average topographic heights at any requested resolution. An exact solution of the ambiguity for the normal height/quasigeoid can be provided by GNSS-levelling.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2018-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90977748","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}
H. Koivula, J. Kuokkanen, Simo Marila, S. Lahtinen, Tuukka Mattila
{"title":"Assessment of sparse GNSS network for network RTK","authors":"H. Koivula, J. Kuokkanen, Simo Marila, S. Lahtinen, Tuukka Mattila","doi":"10.1515/jogs-2018-0014","DOIUrl":"https://doi.org/10.1515/jogs-2018-0014","url":null,"abstract":"Abstract We tested the accuracy and usability of a sparse GNSS reference station network for network RTK (NRTK) using the Finnish permanent GNSS network FinnRef. We modified the configuration of the FinnRef network stations used in NRTK computation. This allowed us to perform the test both inside and outside of the network area using different NRTK methods and two different RTK receivers. In the test area the average distance between the FinnRef stationswas 160 km. As a comparison,we tested also with the commercial Trimnet and HxGN SmartNet positioning services operated by Geotrim Oy and Leica Geosystems Finland, respectively. Tests showed that the horizontal and vertical rms of Trimnet servicewas 16mmand 40 mm, and of HxGN SmartNet service 23mmand 48 mm. The best rms for the sparse NLS (National Land Survey of Finland) Service was 22 mm and 56 mm. These results indicate that a good NRTK solution can be achieved with a sparser network than typically used. This study also indicates, that the methods for NRTK processing can also affect the quality of the solution.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2018-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"77262213","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}
{"title":"Sea level accelerations at globally distributed tide gauge stations during the satellite altimetry era","authors":"H. Iz, C. Shum, C. Kuo","doi":"10.1515/jogs-2018-0013","DOIUrl":"https://doi.org/10.1515/jogs-2018-0013","url":null,"abstract":"Abstract This observational study reports that several globally distributed tide gauge stations exhibit a propensity of statistically significant sea level accelerations during the satellite altimetry era. However, the magnitudes of the estimated tide gauge accelerations during this period are systematically and noticeably smaller than the global mean sea level acceleration reported by recent analyses of satellite altimetry. The differences are likely to be caused by the interannual, decadal and interdecadal sea level variations, which are modeled using a broken trend model with overlapping harmonics in the analyses of tide gauge data but omitted in the analysis of satellite altimetry.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2018-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78383815","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}
{"title":"Fast computation of sine/cosine series coefficients of associated Legendre function of arbitrary high degree and order","authors":"T. Fukushima","doi":"10.1515/jogs-2018-0017","DOIUrl":"https://doi.org/10.1515/jogs-2018-0017","url":null,"abstract":"Abstract In order to accelerate the spherical/spheroidal harmonic synthesis of any function, we developed a new recursive method to compute the sine/cosine series coefficient of the 4π fully- and Schmidt quasi-normalized associated Legendre functions. The key of the method is a set of increasing-degree/order mixed-wavenumber two to four-term recurrence formulas to compute the diagonal terms. They are used in preparing the seed values of the decreasing-order fixed-degree, and fixed-wavenumber two- and three-term recurrence formulas, which are obtained by modifying the classic relations. The new method is accurate and capable to deal with an arbitrary high degree/ order/wavenumber. Also, it runs significantly faster than the previous method of ours utilizing the Wigner d function, say around 20 times more when the maximum degree exceeds 1,000.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2018-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"73550615","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}
{"title":"Geodesic equations and their numerical solution in Cartesian coordinates on a triaxial ellipsoid","authors":"G. Panou, R. Korakitis","doi":"10.1515/jogs-2019-0001","DOIUrl":"https://doi.org/10.1515/jogs-2019-0001","url":null,"abstract":"Abstract In this work, the geodesic equations and their numerical solution in Cartesian coordinates on an oblate spheroid, presented by Panou and Korakitis (2017), are generalized on a triaxial ellipsoid. A new exact analytical method and a new numerical method of converting Cartesian to ellipsoidal coordinates of a point on a triaxial ellipsoid are presented. An extensive test set for the coordinate conversion is used, in order to evaluate the performance of the two methods. The direct geodesic problem on a triaxial ellipsoid is described as an initial value problem and is solved numerically in Cartesian coordinates. The solution provides the Cartesian coordinates and the angle between the line of constant λ and the geodesic, at any point along the geodesic. Also, the Liouville constant is computed at any point along the geodesic, allowing to check the precision of the method. An extensive data set of geodesics is used, in order to demonstrate the validity of the numerical method for the geodesic problem. We conclude that a complete, stable and precise solution of the problem is accomplished.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2018-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"72741609","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}
{"title":"The changing mass of glaciers on the Tibetan Plateau, 2002–2016, using time-variable gravity from the GRACE satellite mission","authors":"Alyson K. Beveridge, C. Harig, F. Simons","doi":"10.1515/jogs-2018-0010","DOIUrl":"https://doi.org/10.1515/jogs-2018-0010","url":null,"abstract":"Abstract The Tibetan Plateau is the largest region of high elevation in the world. The source of water for a number of important rivers, the Himalayan region is vital to the billions of inhabitants of the Asian continent. Over the last fifty years, the climate in the region has warmed more rapidly than anywhere else at the same latitude. Causes and effects, and the geographical details of these alarming warming trends are as yet not fully known. One way of assessing the effects of climate change in this area is to measure the change in glacier volume in the region, but estimates made on the basis of different techniques have not been conclusive to date, and remain difficult to reconcile. We examine the temporal behavior of the mass flux integrated over four distinct groupings of Tibetan glaciers using satellite gravimetry from the Gravity Recovery and Climate Experiment (GRACE). We use a technique known as spatio-spectral localization using spherical Slepian functions to convert global spherical harmonic expansions of the time-dependent geopotential into monthly estimates of mass changes over the Tibetan Plateau. Subsequent reductions are aimed at interpreting this mass change as due to gains or losses in ice mass. We find that (ice) mass has been decreasing on the Tibetan Plateau between 2002 and 2016 but with significant spatial variability throughout the region. Specifically, in the regions of Himalaya, Pamir, Qilian, and Tien Shan, glaciers have been losing ice mass at a rate of −11±3, −1±2, +8±2, and −6±1 Gt/yr, respectively, over the last decade.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2018-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"88096122","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}