The sensitivity of ocean tide loading displacements to the structure of the upper mantle and crust of Taiwan Island

Ocean tide loading (OTL) displacements are sensitive to the shallow structure of the solid Earth; hence, the high-resolution spatial pattern of OTL displacement can provide knowledge to constrain the shallow Earth structure, especially in coastal areas. In this study, we investigate the sensitivity of the modeled M2 OTL displacement over Taiwan Island to perturbations of three physical quantities, namely, the density, bulk modulus, and shear modulus in the upper mantle and crust. Then, we compare the sensitivity of the modeled M2 OTL displacement to Earth models with the sensitivity to ocean tide models using root mean square (RMS) differences. We compute the displacement Green’s function and OTL displacement relative to the center of mass of the solid Earth (CE) reference frame, analyze the sensitivity to the three physical quantities in the CRUST1.0 model and the Preliminary Reference Earth Model (PREM), and present their spatial patterns. We find that displacement Green’s functions and OTL displacements are more sensitive to the two elastic moduli than the density in the upper mantle and crust. Moreover, their distinctive sensitivity patterns suggest that the three physical quantities might be constrained independently. The specific relationships between the perturbed structural depths and the distance ranges of peak sensitivities from the observation points to the coastline revealed by the shear modulus can mitigate the nonuniqueness problem in inversion. In particular, the horizontal tidal components observed by the Global Positioning System (GPS) can yield better results in inversions than the vertical component owing to the smaller OTL model errors and the higher structural sensitivity (except for the shear modulus in the asthenosphere).

Estimating ocean tide loading displacements with GPS and GLONASS

Ground displacements due to ocean tide loading have previously been successfully observed using Global Positioning System (GPS) data, and such estimates for the principal lunar M2 constituent have been used to infer the rheology and structure of the asthenosphere. The GPS orbital repeat period is close to that of several other major tidal constituents (K1, K2, S2); thus, GPS estimates of ground displacement at these frequencies are subject to GPS systematic errors. We assess the addition of GLONASS (GLObal NAvigation Satellite System) to increase the accuracy and reliability of eight major ocean tide loading constituents: four semi-diurnal (M2, S2, N2, K2) and four diurnal constituents (K1, O1, P1, Q1). We revisit a previous GPS study, focusing on 21 sites in the UK and western Europe, expanding it with an assessment of GLONASS and GPS+GLONASS estimates. In the region, both GPS and GLONASS data have been abundant since 2010.0. We therefore focus on the period 2010.0–2014.0, a span considered long enough to reliably estimate the major constituents. Data were processed with a kinematic precise point positioning (PPP) strategy to produce site coordinate time series for each of three different modes: GPS, GLONASS and GPS+GLONASS. The GPS solution with ambiguities resolved was used as a baseline for performance assessment of the additional modes. GPS+GLONASS shows very close agreement with ambiguity resolved GPS for lunar constituents (M2, N2, O1, Q1) but with substantial differences for solar-related constituents (S2, K2, K1, P1), with solutions including GLONASS being generally closer to model estimates. While no single constellation mode performs best for all constituents and components, we propose to use a combination of constellation modes to recover tidal parameters: GPS+GLONASS for most constituents, except for K2 and K1 where GLONASS (north and up) and GPS with ambiguities resolved (east) perform best.