Seasonal low-degree changes in terrestrial water mass load from global GNSS measurements

Large-scale mass redistribution in the terrestrial water storage (TWS) leads to changes in the low-degree spherical harmonic coefficients of the Earth’s surface mass density field. Studying these low-degree fluctuations is an important task that contributes to our understanding of continental hydrology. In this study, we use global GNSS measurements of vertical and horizontal crustal displacements that we correct for atmospheric and oceanic effects, and use a set of modified basis functions similar to Clarke et al. (Geophys J Int 171:1–10, 2007) to perform an inversion of the corrected measurements in order to recover changes in the coefficients of degree-0 (hydrological mass change), degree-1 (centre of mass shift) and degree-2 (flattening of the Earth) caused by variations in the TWS over the period January 2003–January 2015. We infer from the GNSS-derived degree-0 estimate an annual variation in total continental water mass with an amplitude of \((3.49 \pm 0.19) \times 10^{3}\) Gt and a phase of \(70^{\circ } \pm 3^{\circ }\) (implying a peak in early March), in excellent agreement with corresponding values derived from the Global Land Data Assimilation System (GLDAS) water storage model that amount to \((3.39 \pm 0.10) \times 10^{3}\) Gt and \(71^{\circ } \pm 2^{\circ }\), respectively. The degree-1 coefficients we recover from GNSS predict annual geocentre motion (i.e. the offset change between the centre of common mass and the centre of figure) caused by changes in TWS with amplitudes of \(0.69 \pm 0.07\) mm for GX, \(1.31 \pm 0.08\) mm for GY and \(2.60 \pm 0.13\) mm for GZ. These values agree with GLDAS and estimates obtained from the combination of GRACE and the output of an ocean model using the approach of Swenson et al. (J Geophys Res 113(B8), 2008) at the level of about 0.5, 0.3 and 0.9 mm for GX, GY and GZ, respectively. Corresponding degree-1 coefficients from SLR, however, generally show higher variability and predict larger amplitudes for GX and GZ. The results we obtain for the degree-2 coefficients from GNSS are slightly mixed, and the level of agreement with the other sources heavily depends on the individual coefficient being investigated. The best agreement is observed for \(T_{20}^C\) and \(T_{22}^S\), which contain the most prominent annual signals among the degree-2 coefficients, with amplitudes amounting to \((5.47 \pm 0.44) \times 10^{-3}\) and \((4.52 \pm 0.31) \times 10^{-3}\) m of equivalent water height (EWH), respectively, as inferred from GNSS. Corresponding agreement with values from SLR and GRACE is at the level of or better than \(0.4 \times 10^{-3}\) and \(0.9 \times 10^{-3}\) m of EWH for \(T_{20}^C\) and \(T_{22}^S\), respectively, while for both coefficients, GLDAS predicts smaller amplitudes. Somewhat lower agreement is obtained for the order-1 coefficients, \(T_{21}^C\) and \(T_{21}^S\), while our GNSS inversion seems unable to reliably recover \(T_{22}^C\). For all the coefficients we consider, the GNSS-derived estimates from the modified inversion approach are more consistent with the solutions from the other sources than corresponding estimates obtained from an unconstrained standard inversion.     

The unexpected signal in GRACE estimates of $$C_{20}$$

For science applications of the gravity recovery and climate experiment (GRACE) monthly solutions, the GRACE estimates of \(C_{20}\) (or \(J_{2}\)) are typically replaced by the value determined from satellite laser ranging (SLR) due to an unexpectedly strong, clearly non-geophysical, variation at a period of \(\sim \)160 days. This signal has sometimes been referred to as a tide-like variation since the period is close to the perturbation period on the GRACE orbits due to the spherical harmonic coefficient pair \(C_{22}/S_{22}\) of S2 ocean tide. Errors in the S2 tide model used in GRACE data processing could produce a significant perturbation to the GRACE orbits, but it cannot contribute to the \(\sim \)160-day signal appearing in \(C_{20}\). Since the dominant contribution to the GRACE estimate of \(C_{20}\) is from the global positioning system tracking data, a time series of 138 monthly solutions up to degree and order 10 (\(10\times 10\)) were derived along with estimates of ocean tide parameters up to degree 6 for eight major tides. The results show that the \(\sim \)160-day signal remains in the \(C_{20}\) time series. Consequently, the anomalous signal in GRACE \(C_{20}\) cannot be attributed to aliasing from the errors in the S2 tide. A preliminary analysis of the cross-track forces acting on GRACE and the cross-track component of the accelerometer data suggests that a temperature-dependent systematic error in the accelerometer data could be a cause. Because a wide variety of science applications relies on the replacement values for \(C_{20}\), it is essential that the SLR estimates are as reliable as possible. An ongoing concern has been the influence of higher degree even zonal terms on the SLR estimates of \(C_{20}\), since only \(C_{20}\) and \(C_{40}\) are currently estimated. To investigate whether a better separation between \(C_{20}\) and the higher-degree terms could be achieved, several combinations of additional SLR satellites were investigated. In addition, a series of monthly gravity field solutions (\(60\times 60\)) were estimated from a combination of GRACE and SLR data. The results indicate that the combination of GRACE and SLR data might benefit the resonant orders in the GRACE-derived gravity fields, but it appears to degrade the recovery of the \(C_{20}\) variations. In fact, the results suggest that the poorer recovery of \(C_{40}\) by GRACE, where the annual variation is significantly underestimated, may be affecting the estimates of \(C_{20}\). Consequently, it appears appropriate to continue using the SLR-based estimates of \(C_{20}\), and possibly also \(C_{40}\), to augment the existing GRACE mission.     

Impact of ambient temperature on spring-based relative gravimeter measurements

In this paper, we investigate the impact of ambient temperature changes on the gravity reading of spring-based relative gravimeters. Controlled heating experiments using two Scintrex CG5 gravimeters allowed us to determine a linear correlation (R \(^{2}>\) 0.9) between ambient temperature and gravity variations. The relation is stable and constant for the two CG5 we used: ?5 nm/s\(^{2}/^\circ \)C. A linear relation is also seen between gravity and residual sensor temperature variations (R \(^{2}>\) 0.75), but contrary to ambient temperature, this relation is neither constant over time nor similar between the two instruments. The linear correction of ambient temperature on the controlled heating time series reduced the standard deviation at least by a factor of 2, to less than 10 nm/s\(^{2}\). The laboratory results allowed for reprocessing the data gathered on a field survey that originally aimed to characterize local hydrological heterogeneities on a karstic area. The correction of two years of monthly CG5 measurements from ambient temperature variations halved the standard deviation (from 62 to 32 nm/s\(^{2}\)) and led us to a better hydrological interpretation. Although the origin of this effect is uncertain, we suggest that an imperfect control of the sensor temperature may be involved, as well as a change of the properties of an electronic component.     

A high-precision digital astrogeodetic traverse in an area of steep geoid gradients close to the coast of Perth,Western Australia

We present results from a new vertical deflection (VD) traverse observed in Perth, Western Australia, which is the first of its kind in the Southern Hemisphere. A digital astrogeodetic QDaedalus instrument was deployed to measure VDs with \({\sim }\)0.2\(''\) precision at 39 benchmarks with a \({{\sim }}1~\hbox {km}\) spacing. For the conversion of VDs to quasigeoid height differences, the method of astronomical–topographical levelling was applied, based on topographical information from the Shuttle Radar Topography Mission. The astronomical quasigeoid heights are in 20–30 mm (RMS) agreement with three independent gravimetric quasigeoid models, and the astrogeodetic VDs agree to 0.2–0.3\(''\) (north–south) and 0.6–0.9\(''\) (east–west) RMS. Tilt-like biases of \({\sim }1\,\,\hbox {mm}\) over \({\sim }1\,\,\hbox {km}\) are present for all quasigeoid models within \({\sim }20\,\,\hbox {km}\) of the coastline, suggesting inconsistencies in the coastal zone gravity data. The VD campaign in Perth was designed as a low-cost effort, possibly allowing replication in other Southern Hemisphere countries (e.g., Asia, Africa, South America and Antarctica), where VD data are particularly scarce.     

Precise and efficient evaluation of gravimetric quantities at arbitrarily scattered points in space

Gravimetric quantities are commonly represented in terms of high degree surface or solid spherical harmonics. After EGM2008, such expansions routinely extend to spherical harmonic degree 2190, which makes the computation of gravimetric quantities at a large number of arbitrarily scattered points in space using harmonic synthesis, a very computationally demanding process. We present here the development of an algorithm and its associated software for the efficient and precise evaluation of gravimetric quantities, represented in high degree solid spherical harmonics, at arbitrarily scattered points in the space exterior to the surface of the Earth. The new algorithm is based on representation of the quantities of interest in solid ellipsoidal harmonics and application of the tensor product trigonometric needlets. A FORTRAN implementation of this algorithm has been developed and extensively tested. The capabilities of the code are demonstrated using as examples the disturbing potential T, height anomaly \(\zeta \), gravity anomaly \(\Delta g\), gravity disturbance \(\delta g\), north–south deflection of the vertical \(\xi \), east–west deflection of the vertical \(\eta \), and the second radial derivative \(T_{rr}\) of the disturbing potential. After a pre-computational step that takes between 1 and 2 h per quantity, the current version of the software is capable of computing on a standard PC each of these quantities in the range from the surface of the Earth up to 544 km above that surface at speeds between 20,000 and 40,000 point evaluations per second, depending on the gravimetric quantity being evaluated, while the relative error does not exceed \(10^{-6}\) and the memory (RAM) use is 9.3 GB.     

A neural network model for predicting weighted mean temperature

Water vapor is an important element of the Earth’s atmosphere, and most of it concentrates at the bottom of the troposphere. Knowledge of the water vapor measured by Global Navigation Satellite Systems (GNSS) is an important direction of GNSS research. In particular, when the zenith wet delay is converted to precipitable water vapor, the weighted mean temperature \(T_\mathrm{m}\) is a variable parameter to be determined in this conversion. The purpose of the study is getting a more accurate \(T_\mathrm{m}\) model for global users by a combination of two different characteristics of \(T_\mathrm{m}\) (i.e., the \(T_\mathrm{m}\) seasonal variations and the relationships between \(T_\mathrm{m}\) and surface meteorological elements). The modeling process was carried out by using the neural network technology. A multilayer feedforward neural network model (the NN) was established. The NN model is used with measurements of only surface temperature \(T_\mathrm{S}\). The NN was validated and compared with four other published global \(T_\mathrm{m}\) models. The results show that the NN performed better than any of the four compared models on the global scale.     

Geodetic methods to determine the relativistic redshift at the level of 10$$^{}$$ in the context of international timescales: a review and practical results

The frequency stability and uncertainty of the latest generation of optical atomic clocks is now approaching the one part in \(10^{18}\) level. Comparisons between earthbound clocks at rest must account for the relativistic redshift of the clock frequencies, which is proportional to the corresponding gravity (gravitational plus centrifugal) potential difference. For contributions to international timescales, the relativistic redshift correction must be computed with respect to a conventional zero potential value in order to be consistent with the definition of Terrestrial Time. To benefit fully from the uncertainty of the optical clocks, the gravity potential must be determined with an accuracy of about \(0.1\,\hbox {m}^{2}\,\hbox {s}^{-2}\), equivalent to about 0.01 m in height. This contribution focuses on the static part of the gravity field, assuming that temporal variations are accounted for separately by appropriate reductions. Two geodetic approaches are investigated for the derivation of gravity potential values: geometric levelling and the Global Navigation Satellite Systems (GNSS)/geoid approach. Geometric levelling gives potential differences with millimetre uncertainty over shorter distances (several kilometres), but is susceptible to systematic errors at the decimetre level over large distances. The GNSS/geoid approach gives absolute gravity potential values, but with an uncertainty corresponding to about 2 cm in height. For large distances, the GNSS/geoid approach should therefore be better than geometric levelling. This is demonstrated by the results from practical investigations related to three clock sites in Germany and one in France. The estimated uncertainty for the relativistic redshift correction at each site is about \(2 \times 10^{-18}\).     

The gravity field model IGGT_R1 based on the second invariant of the GOCE gravitational gradient tensor

Based on tensor theory, three invariants of the gravitational gradient tensor (IGGT) are independent of the gradiometer reference frame (GRF). Compared to traditional methods for calculation of gravity field models based on the gravity field and steady-state ocean circulation explorer (GOCE) data, which are affected by errors in the attitude indicator, using IGGT and least squares method avoids the problem of inaccurate rotation matrices. The IGGT approach as studied in this paper is a quadratic function of the gravity field model’s spherical harmonic coefficients. The linearized observation equations for the least squares method are obtained using a Taylor expansion, and the weighting equation is derived using the law of error propagation. We also investigate the linearization errors using existing gravity field models and find that this error can be ignored since the used a-priori model EIGEN-5C is sufficiently accurate. One problem when using this approach is that it needs all six independent gravitational gradients (GGs), but the components \(V_{xy}\) and \(V_{yz}\) of GOCE are worse due to the non-sensitive axes of the GOCE gradiometer. Therefore, we use synthetic GGs for both inaccurate gravitational gradient components derived from the a-priori gravity field model EIGEN-5C. Another problem is that the GOCE GGs are measured in a band-limited manner. Therefore, a forward and backward finite impulse response band-pass filter is applied to the data, which can also eliminate filter caused phase change. The spherical cap regularization approach (SCRA) and the Kaula rule are then applied to solve the polar gap problem caused by GOCE’s inclination of \(96.7^{\circ }\). With the techniques described above, a degree/order 240 gravity field model called IGGT_R1 is computed. Since the synthetic components of \(V_{xy}\) and \(V_{yz}\) are not band-pass filtered, the signals outside the measurement bandwidth are replaced by the a-priori model EIGEN-5C. Therefore, this model is practically a combined gravity field model which contains GOCE GGs signals and long wavelength signals from the a-priori model EIGEN-5C. Finally, IGGT_R1’s accuracy is evaluated by comparison with other gravity field models in terms of difference degree amplitudes, the geostrophic velocity in the Agulhas current area, gravity anomaly differences as well as by comparison to GNSS/leveling data.     

A first modeling of dynamic and static crustal strain field from near-field dilatation measurements: example of the 2013 $$M_w$$ 6.2 Ruisui earthquake,Taiwan

We analyze the high-resolution dilatation data for the October 2013 \(M_w\) 6.2 Ruisui, Taiwan, earthquake, which occurred at a distance of 15–20 km away from a Sacks–Evertson dilatometer network. Based on well-constrained source parameters (\(\hbox {strike}=217^\circ \), \(\hbox {dip}=48^\circ \), \(\hbox {rake}=49^\circ \)), we propose a simple rupture model that explains the permanent static deformation and the dynamic vibrations at short period (\(\sim \)3.5–4.5 s) for most of the four sites with less than 20 % of discrepancies. This study represents a first attempt of modeling simultaneously the dynamic and static crustal strain using dilatation data. The results illustrate the potential for strain recordings of high-frequency seismic waves in the near-field of an earthquake to add constraints on the properties of seismic sources.     

On the consistency of the current conventional EOP series and the celestial and terrestrial reference frames

Precise transformation between the celestial reference frames (CRF) and terrestrial reference frames (TRF) is needed for many purposes in Earth and space sciences. According to the Global Geodetic Observing System (GGOS) recommendations, the accuracy of positions and stability of reference frames should reach 1 mm and 0.1 mm year\(^{-1}\), and thus, the Earth Orientation Parameters (EOP) should be estimated with similar accuracy. Different realizations of TRFs, based on the combination of solutions from four different space geodetic techniques, and CRFs, based on a single technique only (VLBI, Very Long Baseline Interferometry), might cause a slow degradation of the consistency among EOP, CRFs, and TRFs (e.g., because of differences in geometry, orientation and scale) and a misalignment of the current conventional EOP series, IERS 08 C04. We empirically assess the consistency among the conventional reference frames and EOP by analyzing the record of VLBI sessions since 1990 with varied settings to reflect the impact of changing frames or other processing strategies on the EOP estimates. Our tests show that the EOP estimates are insensitive to CRF changes, but sensitive to TRF variations and unmodeled geophysical signals at the GGOS level. The differences between the conventional IERS 08 C04 and other EOP series computed with distinct TRF settings exhibit biases and even non-negligible trends in the cases where no differential rotations should appear, e.g., a drift of about 20 \(\upmu \)as year\(^{-1 }\)in \(y_{\mathrm{pol }}\) when the VLBI-only frame VTRF2008 is used. Likewise, different strategies on station position modeling originate scatters larger than 150 \(\upmu \)as in the terrestrial pole coordinates.     

Rectangular rotation of spherical harmonic expansion of arbitrary high degree and order

In order to move the polar singularity of arbitrary spherical harmonic expansion to a point on the equator, we rotate the expansion around the y-axis by \(90^{\circ }\) such that the x-axis becomes a new pole. The expansion coefficients are transformed by multiplying a special value of Wigner D-matrix and a normalization factor. The transformation matrix is unchanged whether the coefficients are \(4 \pi \) fully normalized or Schmidt quasi-normalized. The matrix is recursively computed by the so-called X-number formulation (Fukushima in J Geodesy 86: 271–285, 2012a). As an example, we obtained \(2190\times 2190\) coefficients of the rectangular rotated spherical harmonic expansion of EGM2008. A proper combination of the original and the rotated expansions will be useful in (i) integrating the polar orbits of artificial satellites precisely and (ii) synthesizing/analyzing the gravitational/geomagnetic potentials and their derivatives accurately in the high latitude regions including the arctic and antarctic area.     

The first Australian gravimetric quasigeoid model with location-specific uncertainty estimates

We describe the computation of the first Australian quasigeoid model to include error estimates as a function of location that have been propagated from uncertainties in the EGM2008 global model, land and altimeter-derived gravity anomalies and terrain corrections. The model has been extended to include Australia’s offshore territories and maritime boundaries using newer datasets comprising an additional \({\sim }\)280,000 land gravity observations, a newer altimeter-derived marine gravity anomaly grid, and terrain corrections at \(1^{\prime \prime }\times 1^{\prime \prime }\) resolution. The error propagation uses a remove–restore approach, where the EGM2008 quasigeoid and gravity anomaly error grids are augmented by errors propagated through a modified Stokes integral from the errors in the altimeter gravity anomalies, land gravity observations and terrain corrections. The gravimetric quasigeoid errors (one sigma) are 50–60 mm across most of the Australian landmass, increasing to \({\sim }100\) mm in regions of steep horizontal gravity gradients or the mountains, and are commensurate with external estimates.     

Contribution of solar radiation and geomagnetic activity to global structure of 27-day variation of ionosphere

Twenty-seven-day variation caused by solar rotation is one of the main periodic effects of solar radiation influence on the ionosphere, and there have been many studies on this periodicity using peak electron density \(\mathrm{N_{m}F_{2}}\) and solar radio flux index F10.7. In this paper, the global electron content (GEC) and observation of Solar EUV Monitor (SEM) represent the whole ionosphere and solar EUV flux, respectively, to investigate the 27-day variation. The 27-day period components of indices \((\hbox {GEC}_{27}\), \(\hbox {SEM}_{27}\), \(\hbox {F10.7}_{27}\), \(\hbox {Ap}_{27})\) are obtained using Chebyshev band-pass filter. The comparison of regression results indicates that the index SEM has higher coherence than F10.7 with 27-day variation of the ionosphere. The regression coefficients of \(\hbox {SEM}_{27 }\) varied from 0.6 to 1.4 and the coefficients of \(\hbox {Ap}_{27}\) varied from \({-}\)0.6 to 0.3, which suggests that EUV radiation seasonal variations are the primary driver for the 27-day variations of the ionosphere for most periods. TEC map grid points on three meridians where IGS stations are dense are selected for regression, and the results show that the contribution of solar EUV radiation is positive at all geomagnetic latitudes and larger than geomagnetic activity in most latitudes. The contribution of geomagnetic activity is negative at high geomagnetic latitude, increasing with decreasing geomagnetic latitudes, and positive at low geomagnetic latitudes. The global structure of 27-day variation of ionosphere is presented and demonstrates that there are two zonal anomaly regions along with the geomagnetic latitudes lines and two peaks in the north of Southeast Asia and the Middle Pacific where \(\hbox {TEC}_{27}\) magnitude values are notably larger than elsewhere along zonal anomaly regions.     

Adding source positions to the IVS combination—First results

The consistent estimation of terrestrial reference frames (TRF), celestial reference frames (CRF) and Earth orientation parameters (EOP) is still an open subject and offers a large field of investigations. Until now, source positions resulting from Very Long Baseline Interferometry (VLBI) observations are not routinely combined on the level of normal equations in the same way as it is a common process for station coordinates and EOPs. The combination of source positions based on VLBI observations is now integrated in the IVS combination process. We present the studies carried out to evaluate the benefit of the combination compared to individual solutions. On the level of source time series, improved statistics regarding weighted root mean square have been found for the combination in comparison with the individual contributions. In total, 67 stations and 907 sources (including 291 ICRF2 defining sources) are included in the consistently generated CRF and TRF covering 30 years of VLBI contributions. The rotation angles \(A_1\), \(A_2\) and \(A_3\) relative to ICRF2 are ?12.7, 51.7 and 1.8 \({\upmu }\) as, the drifts \(D_\alpha \) and \(D_\delta \) are ?67.2 and 19.1 \(\upmu \) as/rad and the bias \(B_\delta \) is 26.1 \(\upmu \) as. The comparison of the TRF solution with the IVS routinely combined quarterly TRF solution shows no significant impact on the TRF, when the CRF is estimated consistently with the TRF. The root mean square value of the post-fit station coordinate residuals is 0.9 cm.     

Performance evaluation of GNSS-TEC estimation techniques at the grid point in middle and low latitudes during different geomagnetic conditions

Global Navigation Satellite Systems (GNSS) have become a powerful tool use in surveying and mapping, air and maritime navigation, ionospheric/space weather research and other applications. However, in some cases, its maximum efficiency could not be attained due to some uncorrelated errors associated with the system measurements, which is caused mainly by the dispersive nature of the ionosphere. Ionosphere has been represented using the total number of electrons along the signal path at a particular height known as Total Electron Content (TEC). However, there are many methods to estimate TEC but the outputs are not uniform, which could be due to the peculiarity in characterizing the biases inside the observables (measurements), and sometimes could be associated to the influence of mapping function. The errors in TEC estimation could lead to wrong conclusion and this could be more critical in case of safety-of-life application. This work investigated the performance of Ciraolo’s and Gopi’s GNSS-TEC calibration techniques, during 5 geomagnetic quiet and disturbed conditions in the month of October 2013, at the grid points located in low and middle latitudes. The data used are obtained from the GNSS ground-based receivers located at Borriana in Spain (40\(^{\circ }\)N, 0\(^{\circ }\)E; mid latitude) and Accra in Ghana (5.50\(^{\circ }\)N, ?0.20\(^{\circ }\)E; low latitude). The results of the calibrated TEC are compared with the TEC obtained from European Geostationary Navigation Overlay System Processing Set (EGNOS PS) TEC algorithm, which is considered as a reference data. The TEC derived from Global Ionospheric Maps (GIM) through International GNSS service (IGS) was also examined at the same grid points. The results obtained in this work showed that Ciraolo’s calibration technique (a calibration technique based on carrier-phase measurements only) estimates TEC better at middle latitude in comparison to Gopi’s technique (a calibration technique based on code and carrier-phase measurements). At the same time, Gopi’s calibration was also found more reliable in low latitude than Ciraolo’s technique. In addition, the TEC derived from IGS GIM seems to be much reliable in middle-latitude than in low-latitude region.     

VMF3/GPT3: refined discrete and empirical troposphere mapping functions

Incorrect modeling of troposphere delays is one of the major error sources for space geodetic techniques such as Global Navigation Satellite Systems (GNSS) or Very Long Baseline Interferometry (VLBI). Over the years, many approaches have been devised which aim at mapping the delay of radio waves from zenith direction down to the observed elevation angle, so-called mapping functions. This paper contains a new approach intended to refine the currently most important discrete mapping function, the Vienna Mapping Functions 1 (VMF1), which is successively referred to as Vienna Mapping Functions 3 (VMF3). It is designed in such a way as to eliminate shortcomings in the empirical coefficients b and c and in the tuning for the specific elevation angle of \(3^{\circ }\). Ray-traced delays of the ray-tracer RADIATE serve as the basis for the calculation of new mapping function coefficients. Comparisons of modeled slant delays demonstrate the ability of VMF3 to approximate the underlying ray-traced delays more accurately than VMF1 does, in particular at low elevation angles. In other words, when requiring highest precision, VMF3 is to be preferable to VMF1. Aside from revising the discrete form of mapping functions, we also present a new empirical model named Global Pressure and Temperature 3 (GPT3) on a \(5^{\circ }\times 5^{\circ }\) as well as a \(1^{\circ }\times 1^{\circ }\) global grid, which is generally based on the same data. Its main components are hydrostatic and wet empirical mapping function coefficients derived from special averaging techniques of the respective (discrete) VMF3 data. In addition, GPT3 also contains a set of meteorological quantities which are adopted as they stand from their predecessor, Global Pressure and Temperature 2 wet. Thus, GPT3 represents a very comprehensive troposphere model which can be used for a series of geodetic as well as meteorological and climatological purposes and is fully consistent with VMF3.     

Calibration for the shear strain of 3-component borehole strainmeters in eastern Taiwan through Earth and ocean tidal waveform modeling

We propose an approach for calibrating the horizontal tidal shear components [(differential extension (\(\gamma _1\)) and engineering shear (\(\gamma _2\))] of two Sacks–Evertson (in Pap Meteorol Geophys 22:195–208, 1971) SES-3 borehole strainmeters installed in the Longitudinal Valley in eastern Taiwan. The method is based on the waveform reconstruction of the Earth and ocean tidal shear signals through linear regressions on strain gauge signals, with variable sensor azimuth. This method allows us to derive the orientation of the sensor without any initial constraints and to calibrate the shear strain components \(\gamma _1\) and \(\gamma _2\) against \(M_2\) tidal constituent. The results illustrate the potential of tensor strainmeters for recording horizontal tidal shear strain.     

Enriching the semantics of the directed polyline–polygon topological relationships: the $$\mathcal {}$$-intersection matrix

In this paper, we define an intersection matrix for enriching the semantics of the topological relationships between a directed polyline and a polygon. In particular, we propose the \(\mathcal {DLP}\)-intersection matrix which enables us to model the origin and destination points, as well as the right- and left-hand sides of the directed polyline. This matrix overcomes the limitation of the well-known DE-9IM, because it allows the representation of the different dimensions of the intersection results at the same time. Accordingly, the geo-operators have been revised and extended in order to address the notions of right- and left-hand sides of a directed polyline, as well as additional notions related to the orientation of the polyline. The \(\mathcal {DLP}\)-intersection matrix has been implemented by extending the Java Topology Suite methods in order to address the new geo-operators based on the notion of orientation.     

A low-complexity algorithm for fast acquisition of weak DSSS signal in high dynamic environment

The initial acquisition of direct-sequence spread-spectrum (DSSS) signal transmitted in bursts by ground terminals at satellite-borne receiver poses an engineering challenge. We propose a low-complexity acquisition algorithm that is capable of capturing extremely weak DSSS signal in the presence of large Doppler dynamics. The algorithm uses fast Fourier transform (FFT)-based frequency-domain techniques to implement circular correlations between the received signal and the local pseudo-random noise (PRN) code, and it coherently accumulates the correlation results across multiple PRN code periods, to achieve a sufficient signal–noise ratio for reliable acquisition. We invoke another FFT procedure to perform the coherent accumulation and the fine compensation for the residual Doppler frequency offset. To highlight the advantage of the proposed algorithm, we make a complexity comparison among the proposed algorithm and two other benchmark strategies, namely the modified double block zero padding (MDBZP) and two-dimensional exhaustive search (2D-ES). It is shown that the proposed algorithm has the lowest complexity, which is particularly desirable for satellite-borne receivers where the computational resource is limited. The acquisition performance of the proposed algorithm is verified by theoretical analysis and Monte Carlo simulations and compared with that of MDBZP and 2D-ES. Moreover, we have carried out extensive tests on a hardware verification system, and we show the claimed tradeoff between performance and cost is indeed attainable with the suggested algorithm. Numerically, it is found the proposed algorithm can achieve a detection rate of 0.9 and a false alarm rate of \(10^{ - 5}\) at C/N ₀ = 29.5 dBHz over a Doppler frequency offset range of \(\left[ { - 7.5\,{\text{kHz}},7.5\,{\text{kHz}}} \right]\) in floating-point simulation, which coincides with the analytical results. The same performance is achieved at C/N ₀ = 31 dBHz in fixed-point simulation and at C/N ₀ = 31.5 dBHz on a hardware system.     

Model improvements and validation of TerraSAR-X precise orbit determination

The radar imaging satellite mission TerraSAR-X requires precisely determined satellite orbits for validating geodetic remote sensing techniques. Since the achieved quality of the operationally derived, reduced-dynamic (RD) orbit solutions limits the capabilities of the synthetic aperture radar (SAR) validation, an effort is made to improve the estimated orbit solutions. This paper discusses the benefits of refined dynamical models on orbit accuracy as well as estimated empirical accelerations and compares different dynamic models in a RD orbit determination. Modeling aspects discussed in the paper include the use of a macro-model for drag and radiation pressure computation, the use of high-quality atmospheric density and wind models as well as the benefit of high-fidelity gravity and ocean tide models. The Sun-synchronous dusk–dawn orbit geometry of TerraSAR-X results in a particular high correlation of solar radiation pressure modeling and estimated normal-direction positions. Furthermore, this mission offers a unique suite of independent sensors for orbit validation. Several parameters serve as quality indicators for the estimated satellite orbit solutions. These include the magnitude of the estimated empirical accelerations, satellite laser ranging (SLR) residuals, and SLR-based orbit corrections. Moreover, the radargrammetric distance measurements of the SAR instrument are selected for assessing the quality of the orbit solutions and compared to the SLR analysis. The use of high-fidelity satellite dynamics models in the RD approach is shown to clearly improve the orbit quality compared to simplified models and loosely constrained empirical accelerations. The estimated empirical accelerations are substantially reduced by 30% in tangential direction when working with the refined dynamical models. Likewise the SLR residuals are reduced from \(-3\,\pm \,17\) to \(2\,\pm \,13\) mm, and the SLR-derived normal-direction position corrections are reduced from 15 to 6 mm, obtained from the 2012–2014 period. The radar range bias is reduced from \(-10.3\) to \(-6.1\) mm with the updated orbit solutions, which coincides with the reduced standard deviation of the SLR residuals. The improvements are mainly driven by the satellite macro-model for the purpose of solar radiation pressure modeling, improved atmospheric density models, and the use of state-of-the-art gravity field models.