A Loading Correction Model for GPS Measurements Derived from Multiple-Data Combined Monthly Gravity

Time-dependent loading deformations of the Earth’s surface, due to nontidal changes in the atmosphere, ocean, land water/ice, etc., contribute significantly to the seasonal and secular Global Positioning System (GPS) site displacements, especially for the up component. While loading deformations derived from general circulation model (GCM) outputs are usually used to correct loading signals in the GPS site displacements, this study aims to provide a loading correction model based on the multiple-data combined monthly gravity products LDCmgm90. We have adopted GPS measurements from 249 IGS reference frame stations and 3 different GCM-based loading models to test the reliability of the LDCmgm90 model. Compared to the GCM-based models, the LDCmgm90 loading correction is more effective in attenuating seasonal (especially annual) loading signals and can bring more significant improvements to most stations for both the data-trend-removed and the data-trend-retained cases. Thus, we have validated the LDCmgm90 model from the loading aspect and proved it to be a reliable loading-correction model for GPS displacements. The relatively better secular loading signals provided by the LDCmgm90 loading model may provide us a chance to study the long-term, nonloading signals in GPS data.

Estimation of crustal deformation parameters and strain build-up in Northwest Himalaya using GNSS data measurements

GPS measurements have proved extremely useful in quantifying strain accumulation rate and assessing seismic hazard in a region. Continuous GPS measurements provide estimates of secular motion used to understand the earthquake and other geodynamic processes. GNSS stations extending from the South of India to the Higher Himalayan region have been used to quantify the strain build-up rate in Central India and the Himalayan region to assess the seismic hazard potential in this realm. Velocity solution has been determined after the application of Markov noise estimated from GPS time series data. The recorded GPS data are processed along with the closest International GNSS stations data for estimation of daily basis precise positioning. The baseline method has been used for the estimation of the linear strain rate between the two stations. Whereas the principal strain axes, maximum shear strain, rotation rate, and crustal shortening rate has been calculated through the site velocity using an independent approach; least-square inversion approach-based triangulation method. The strain rate analysis estimated by the triangulation approach exhibits a mean value of extension rate of 26.08 nano-strain/yr towards N131°, the compression rate of –25.38 nano-strain/yr towards N41°, maximum shear strain rate of 51.47 nano-strain/yr, dilation of –37.57 nano-strain/yr and rotation rate of 0.7°/Ma towards anti-clockwise. The computed strain rate from the Baseline method and the Triangulation method reports an extensive compression rate that gradually increases from the Indo-Gangetic Plain in South to Higher Himalaya in North. The slip deficit rate between India and Eurasia Plate in Kumaun Garhwal Himalaya has been computed as 18±1.5 mm/yr based on elastic dislocation theory. Thus, in this study, present-day surface deformation rate and interseismic strain accumulation rate in the Himalayan region and the Central Indian region have been estimated for seismic hazard analysis using continuous GPS measurements.

Greenland ice velocity maps from the PROMICE project

We present the Programme for Monitoring of the Greenland Ice Sheet (PROMICE) Ice Velocity product (https://doi.org/10.22008/promice/data/sentinel1icevelocity/greenlandicesheet, Solgaard and Kusk, 2021), which is a time series of Greenland Ice Sheet ice velocity mosaics spanning September 2016 through to the present. The product is based on Sentinel-1 synthetic aperture radar data and has a 500 m grid spacing. A new mosaic is available every 12 d and spans two consecutive Sentinel-1 cycles (24 d). The product is made available within ∼ 10 d of the last acquisition and includes all possible 6 and 12 d pairs within the two Sentinel-1A cycles. We describe our operational processing chain from data selection, mosaicking, and error estimation to final outlier removal. The product is validated against in situ GPS measurements. We find that the standard deviation of the difference between satellite- and GPS-derived velocities (and bias) is 20 m yr−1 (−3 m yr−1) and 27 m yr−1 (−2 m yr−1) for the components in an eastern and northern direction, respectively. Over stable ground the values are 8 m yr−1 (0.1 m yr−1) and 12 m yr−1 (−0.6 m yr−1) in an eastern and northern direction, respectively. This is within the expected values; however, we expect that the GPS measurements carry a considerable part of this uncertainty. We investigate variations in coverage from both a temporal and spatial perspective. The best spatial coverage is achieved in winter due to the comprehensive data coverage by Sentinel-1 and high coherence, while summer mosaics have the lowest coverage due to widespread melt. The southeast Greenland Ice Sheet margin, along with other areas of high accumulation and melt, often has gaps in the ice velocity mosaics. The spatial comprehensiveness and temporal consistency make the product ideal both for monitoring and for studying ice-sheet-wide and glacier-specific ice discharge and dynamics of glaciers on seasonal scales.

Slip bursts during coalescence of slow slip events in Cascadia

<p>Both laboratory experiments and dynamic simulations suggest that earthquakes can be preceded by a precursory phase of slow slip. Observing processes leading to an acceleration or spreading of slow slip along faults is therefore key to understand the dynamics potentially leading to seismic ruptures. Here, we use continuous GPS measurements of the ground displacement to image the daily slip along the fault beneath Vancouver Island during a slow slip event in 2013. We image the coalescence of three originally distinct slow slip fronts merging together. We show that during coalescence phases lasting for 2 to 5 days, the rate of energy (moment) release significantly increases. This observation supports the view proposed by theoretical and experimental studies that the coalescence of slow slip fronts is a possible mechanism for initiating earthquakes.</p>