These technologies provide critical insights into the geodynamics of lithospheric plates and enable the detailed characterization of both interplate and intraplate tectonic domains, illuminating complex tectonic processes that govern seismicity and crustal strain. Furthermore, the synergy between different geodetic platforms supports real-time hazard assessment and early warning systems for earthquakes, tsunamis, and volcanic eruptions. This multidimensional approach enhances both scientific understanding of Earth’s dynamic behavior and societal resilience by informing evidence-based disaster risk reduction strategies.
This session invites contributions employing geodetic, geophysical, geological, and seismotectonic data to investigate active deformation zones, including intraplate volcanic settings. We especially welcome interdisciplinary studies that integrate geodesy, seismology, tectonics, and geophysics to better constrain strain accumulation, plate motions, and lithospheric deformation. The overarching goal is to highlight recent geodetic advances and explore their implications for understanding lithospheric geodynamics and enhancing hazard preparedness.
Orals: Thu, 7 May, 10:45–12:30 | Room 0.96/97
Fault zones release tectonic strain through a combination of seismic and aseismic slip. Creeping fault segments may have typically less elastic strain energy accumulated available to rupture in an earthquake compared to largely locked sections. However, where creeping fault segments transition into locked ones, stress rates are the highest as the slip deficit and stored elastic strain show large spatial along-fault gradient.
The North Anatolian Fault Zone (NAFZ) in Türkiye hosts two prominent creeping segments: 1) the Ismetpasa segment on the central part of the NAFZ, which appears devoid of micro-seismicity down to magnitude M=1.8 in the regional catalogs, but hosted the nucleation of two M>7 earthquakes in 1943 and 1944; and 2) the western portion of the submarine Main Marmara Fault, which poses a high seismic risk due to its proximity to the Istanbul metropolitan region. Some of the largest earthquakes of the instrumental era (2019 M5.8 and 2025 M6.2) close to Istanbul nucleated at the eastern edge of this partially creeping segment.
In this study, we combine near-fault dense seismo-geodetic deployments, with deep-learning seismicity catalogs to investigate the role of aseismic deformation in driving the seismicity and controlling the source properties along those creeping segments of the NAFZ. At the Ismetpasa segment, we present the first evidence of significant microseismicity on and up to ~5km off the main NAFZ fault branch (mostly ML < 2) surrounding the creeping patches. This microseismicity is likely driven by the aseismic slip on the main fault plane. We interpret this seismic activity as the signature of a weak, damaged fault zone surrounding the tip of the ruptures of the M>7 1943 and 1944 events.
In the Marmara region, we show a series of eastward propagating M>5 events and a gradual eastward unlocking of the Main Marmara Fault over the last ~15 years. Seismic activity progresses from creeping toward transitional segments and is currently arriving at the locked Princes Islands segment south of Istanbul, which has the potential to host a M~7 earthquake. These findings highlight the role of aseismic slip in modulating the available shear stress and elastic stored energy, which, in turn, control the nucleation and arrest of large ruptures. Our results also illustrate the importance of monitoring fault systems including multi-disciplinary instrumentation that enables capturing the entire frequency band from slow to fast slip.
How to cite: Martínez-Garzón, P., Becker, D., Jolivet, R., Jara, J., Cakir, Z., Chen, X., Nunez-Jara, S., Kartal, R. F., Türker, E., Dresen, G., Ben-Zion, Y., Cotton, F., Kadirioglu, F. T., Kilic, T., and Bohnhoff, M.: Seismic potential of creeping segments of the North Anatolian Fault: insights from seismo-geodetic deployments and deep learning catalogs , EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-8784, https://doi.org/10.5194/egusphere-egu26-8784, 2026.
Satellite geodesy provides critical insights into tectonic deformation, fault activity, and seismic hazard. However, in regions of widespread continental deformation, observational coverage has until recently relied on sparse GNSS measurements, limiting the resolution of short-wavelength deformation features. By integrating InSAR we can greatly improve the resolution, and we have recently constructed a transnational velocity field for the entire Alpine-Himalayan Belt at 1 km spacing, from over 222,000 Sentinel-1 SAR images (2016–2024) and a new compilation of GNSS velocities [Elliott et al., in review]. This dataset spans more than 11,000 km from southwestern Europe to eastern China, covering over 20 million km², and is referenced consistently to the Eurasian frame.
From these velocities, we derive horizontal strain rates, providing near-continuous deformation mapping across the planet’s largest actively deforming region. Results reveal a bimodal pattern of tectonic strain, which is concentrated along major faults in some regions but distributed across broader zones in others. Vertical motions, in contrast, exhibit shorter-wavelength signals dominated by non-tectonic processes, particularly groundwater depletion.
Satellite geodesy also provides critical insights into volcanic deformation and hazard, and we have processed InSAR data for the ~1300 subaerial volcanoes most likely to erupt. Scale is less of an issue for volcanoes, with volcanic activity usually confined to within 40 km of each volcanic centre, but timeliness is important for hazard monitoring, and we process data in near-real time form a subset of volcanoes. For historical analyses we have integrated our InSAR results with local GNSS networks [Bedon et al, in prep], but it remains a challenge to incorporate GNSS from multiple disparate networks for ongoing monitoring on a global basis.
The spatial resolution of InSAR measurements is better than GNSS by orders of magnitude, but inclusion of GNSS is key for two reasons: firstly, for tying InSAR to a global reference frame and secondly, to provide a third component of the velocity field, which allows the full 3-D field to be constrained. However, the combination leads to very different resolutions in the north-south direction, constrained predominantly by GNSS, and the east-west direction, where InSAR dominates. When estimating the strain rate this leads to non-localisation of strain for north-south trending strike-slip faults and east-west trending dip-slip faults but also leads to short wavelength shear strain (e.g., from near-surface creep) being wrongly attributed to dilatation on faults of any orientation [Fang et al., 2024].
We are addressing this issue in two ways. Firstly, by inclusion of along-track velocity estimates from Sentinel-1 burst overlap regions [Nergizci et al., 2024] and secondly by the addition of InSAR velocity measurements from NISAR. The left-looking nature of NISAR acquisitions will provide two more independent velocity measurement vectors that will enable full 3-D estimation at high resolution. Whilst the accuracy in the north-south direction will be ~4 times worse than in the east-west direction, the improvement in resolution will be by orders of magnitude.
References
Elliott et al. (in review). Preprint: doi:10.31223/X5GX6B.
Fang et al (2024). doi:10.1029/2024GL111199.
Nergizci et al. (2024). doi:10.1016/j.procs.2024.06.401.
How to cite: Hooper, A., Elliott, J., Fang, J., Lazecký, M., Wright, T., Espin Bedon, P., Muhammet Nergizci, M., Maghsoudi, Y., Ou, Q., Payne, J., Novoa Lizama, C., Rollins, C., and Wang, D.: Large-scale high-resolution deformation of tectonic and volcanic regions, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-18724, https://doi.org/10.5194/egusphere-egu26-18724, 2026.
Prior to the Tajogaite eruption (19/09-13/12/2021), several unrest signals were detected on La Palma starting in 2017. In particular, seismic swarms associated with variations in gas emissions and ground deformation were documented between 2017 and the onset of the 2021 eruptive process (Torres-González et al., 2020; Fernández et al., 2021). However, the spatio-temporal relationship between the 2021 eruptive activity on La Palma and the regional tectonic framework of the Canary Archipelago has not yet been thoroughly investigated.
The joint analysis of regional- and local-scale continuous GNSS data from permanent open-access stations across the Canary Islands allowed the identification of anomalous inflation on La Palma beginning in 2014. This anomalous signal is temporally correlated with the cessation of sill intrusion events that drove the long-term uplift of El Hierro Island following its submarine eruption (2012–2014). Our analysis explores whether this inflation can be directly related to the Tajogaite eruption—suggesting that post-eruptive processes at El Hierro may have induced dilatation within La Palma’s magmatic plumbing system—or whether it reflects a broader regional uplift associated with deep intrusions and lateral magma transport across the Canary Archipelago.
References
Torres-González, P.A., Luengo-Oroz, N., Lamolda, H. et al. (2020) Unrest signals after 46 years of quiescence at Cumbre Vieja, La Palma, Canary Islands, J. Volcanol. Geotherm. Res., 392, 106757, https://doi.org/10.1016/j.jvolgeores.2019.106757.
Fernández, J., Escayo, J., Hu, Z. et al. (2021) Detection of volcanic unrest onset in La Palma, Canary Islands, evolution and implications. Sci. Rep., 11, 2540. https://doi.org/10.1038/s41598-021-82292-3.
How to cite: Charco, M., Pallero, J. L. G., Santamaría-Gómez, Á., and González, P. J.: Pre-Eruptive Inflation on La Palma (2014–2021): GNSS Evidence for Regional Magmatic Interactions in the Canary Islands, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-20387, https://doi.org/10.5194/egusphere-egu26-20387, 2026.
The Quaternary Olkaria volcanic complex is a high-temperature geothermal system in the Central Kenya Rift located in a structural transition zone between the Mt. Longonot and Mt. Eburru eruptive centers. This region is subject to protracted crustal deformation and ground subsidence. Vertical ground motion is associated with crustal-scale processes, including magmatic intrusion and eruption, as well as normal and transfer faulting. In addition, the region has been the subject of fluid extraction associated with geothermal energy production. In the past 20 years, increased re-injection strategies were implemented in order to slow down drawdown and thereby mitigate against exponential ground subsidence. Although an extensive benchmark network for monitoring surface deformation was established in 1983, the lack of subsequent precise leveling surveys has necessitated the use of state-of-the-art geodetic techniques to quantify the magnitude and temporal evolution of ground deformation to better understand the roles of tectonic and anthropogenically induced land-surface changes.
Here we present new radar interferometry observations spanning the past decade, combining Sentinel-1 data (2016–2026) and TerraSAR-X data (2024–2026), to constrain vertical surface motion at high spatial and temporal resolution. These InSAR time series are complemented by measurements from a local GNSS network installed on and around the Olkaria dome. Our results show that rapid subsidence observed since 2016 slowed markedly around 2020 and has since largely stagnated. High-resolution X-band and persistent scatterer C-band data reveal that localized subsidence persists near fluid-extraction sites, whereas regional subsidence rates in the area of the volcanic complex have decreased by approximately an order of magnitude. Independent, statistically robust GNSS time-series analyses support these observations. We further assess different InSAR processing strategies and highlight the critical importance of rigorous atmospheric correction due to the influence of high seasonal moisture availability.
Overall, our analysis indicates that vertical land-surface deformation in the Olkaria region at annual timescales is primarily driven by deep-seated magmatic processes, while geothermal energy production has only contributed to localized subsidence through fluid extraction.
How to cite: Bookhagen, B., Kimata, J., Omenda, P., Saitet, D., and Strecker, M.: GNSS and InSAR observations at the Olkaria geothermal site in the central Kenya Rift System, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-8227, https://doi.org/10.5194/egusphere-egu26-8227, 2026.
The 2021 Mw 7.4 Maduo earthquake, rupturing within the Bayan Har block, provides a unique opportunity to investigate long-term postseismic deformation and lateral rheological heterogeneity of the lithosphere in northeastern Tibet. Here, we present four years (2021–2025) of postseismic deformation derived from continuous GNSS and InSAR observations to characterize the long-term deformation pattern and its underlying rheological controls. The GNSS and InSAR time series reveal a sustained growth of cumulative postseismic deformation, with both deformation amplitude and spatial extent progressively increasing with time. Maximum line-of-sight deformation inferred from InSAR reaches ~80 mm four years after the earthquake, consistent with the horizontal displacement magnitudes recorded by GNSS. The deformation exhibits a clear temporal transition, characterized by rapid growth during the first 1–2 years, followed by a substantially reduced but persistent deformation rate during years 3–4. During the later stage, localized regions exhibit additional deformation of up to ~5 mm, indicating that postseismic processes continue to operate over multi-year timescales. Modeling of the early postseismic deformation indicates that first-year displacements are jointly controlled by afterslip and viscoelastic relaxation, whereas the contribution from poroelastic rebound is negligible. Rheological inversion of the early postseismic deformation constrains optimal steady-state viscosities of 2–5 × 10¹⁹ Pa s for the lower crust and 3–10 × 10¹⁹ Pa s for the upper mantle, indicating the presence of a mechanically weak lower crust beneath the Bayan Har block. By incorporating the full four-year deformation time series, we further identify pronounced lateral variations in postseismic deformation behavior across the East Kunlun fault. South of the fault, the long-term deformation decay is broadly consistent with a weak lower crust characterized by viscosities on the order of ~10¹⁹ Pa s. In contrast, north of the fault, systematic spatial and temporal misfits between observations and homogeneous rheological models require a substantially stronger lower crust, with effective viscosities on the order of ~10²¹ Pa s. These results indicate that the East Kunlun fault represents a first-order rheological boundary separating laterally contrasting lithospheric domains in northeastern Tibet, and highlight the critical role of long-term GNSS and InSAR observations in resolving lateral rheological heterogeneity that cannot be captured by short-term postseismic data alone.
How to cite: Li, J., Chen, Y., Zhang, Z., Zhan, W., and Deng, Z.: Constraints on Lateral Rheological Heterogeneity in Northeastern Tibet from Long-Term GNSS and InSAR Observations following the 2021 Maduo Earthquake, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-9031, https://doi.org/10.5194/egusphere-egu26-9031, 2026.
Margins between tectonic plates host most large earthquakes recorded in the lithosphere. Over periods of tens to hundreds of years, relative plate motions along portions of crustal seismogenic faults promote the slow accrual of stress (i.e., the inter-seismic stress gain) that is later suddenly released via earthquakes (i.e., the co-seismic stress drop) – a process generally referred to as seismic cycle. Virtually all models of seismic hazard assessment assume that whole-plate motions (i.e., motions that are adequately described via Euler vectors) remain steady over the seismic cycle, and that the impact of inter- and co-seismic stress variations is solely crustal deformation in the vicinity of seismogenic faults. From the standpoint of plate dynamics, however, plate-margin stress variations during the seismic cycle generate torques that may be comparable in magnitude to those associated with viscous stresses at the lithosphere/asthenosphere interface, which resist plate motions. On this basis, it is plausible to hypothesize that whole-plate motions may be susceptible to temporal variations over the seismic cycle. The availability of progressively longer and denser GNSS position time series measured at sites located inside several tectonic plates indeed favor testing such a hypothesis. Here I will show results from recent studies that analyze publicly available GNSS data and infer temporal variations of the motions of several tectonic plates. These changes appear consistent with the torque variations associated with inter- or co-seismic phases of large earthquakes occurred along their margins. I will speculate on whether the link between whole-plate motions and the seismic cycle is robust enough to draw any additional information in the context of models of seismic hazard assessment.
How to cite: Iaffaldano, G.: Testing whole-plate motion steadiness over the seismic cycle, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-12754, https://doi.org/10.5194/egusphere-egu26-12754, 2026.
We present multi-disciplinary case studies for multi-year deployments of an autonomous subsea sensor logging node with acoustic communication capability.
Errors associated with single-point observations of tectonic plate motion, seafloor subsidence and geostrophic current circulation are typically on the same order as expected annual rates of those phenomena.
Hardware platforms suffer from settling rates that contaminate signal for the first few months of deployment. These also require platforms to remain in the same 3D position without recovery and re-deployment to achieve a continuous time series.
Sonardyne Fetch systems have a 10-year lifetime and are used as multi-disciplinary platforms for Acoustic-Ranging, GNSS-A, PIES and self-calibrating (Ambient-Zero-Ambient/AZA) pressure logging. Acoustic modems allow regular data offload without hardware recovery.
Here, we present case studies from academia and industry across disciplines and incorporating manned and unmanned data recovery platforms. We also present projects under active investigation to form permanent near-real-time data communication to pre-existing cabled infrastructure. This will expand the footprint and potential for relocation of subsea observatories while minimising logistical and environmental impact.
How to cite: Reis, W.: Long-term Platforms for Subsea Acoustic Ranging, Geodetic Monitoring and Ocean Circulation: Past and Future Applications, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-20042, https://doi.org/10.5194/egusphere-egu26-20042, 2026.
High-precision seafloor geodesy such as GNSS-A positioning or seabed transponder networks is critically dependent on acoustic travel time measurements to recover seafloor benchmarks; however, in stratified ocean acoustic ranging, sound-speed-profile variations invalidate the approximately linear travel-time-to-range mapping that commonly holds for the terrestrial case. Simplified sound-speed models therefore yield biases that could be aliased into deformation time series, leading to errors in inference of plate motions and lithospheric deformation. Furthermore, underwater transponders deployed within the geodetic network are unable to maintain strict time synchronization in cold, high-pressure deep-ocean environments, imposing an added challenge on the usage of travel-time ranging.
To overcome these limitations, we develop an equivalent-gradient framework with a closed-form delay–range relationship and represent synchronization imperfections by a lumped time-bias term, enabling joint recovery of seafloor transponder position(s) and the bias. Specifically, let k = 1,...,K index the reception epochs along a moving surface vessel trajectory; sk denotes the GNSS-referenced vessel position and tk the recorded one-way arrival timestamp from a fixed seafloor transponder. We then form inter-epoch TDOA measurements that eliminate the unknown transmit epoch and reduce the problem to estimating a reference one-way delay τ0 together with the transponder location u. Under the equivalent-gradient framework, travel time is efficiently mapped to an slant range dk = Req (τ; ξk), , where ξk collects the equivalent-gradient parameters derived from the layered SSP, yielding the range–geometry constraints dk^2 = u − s^2. A squared-difference with respect to a reference epoch leads to a stable pseudo-linear regression:
This yields a WLS closed-form initializer followed by weighted Gauss–Newton refinement. An SDR-based global initializer is also developed, offering complementary insight into the problem’s geometry. The approach accommodates different acoustic link geometries (e.g., ship-to-seafloor and AUV-to-seafloor) and can exploit identifiable multipath (e.g., surface-reflected arrivals) for additional constraints. Monte-Carlo simulations under realistic stratified SSPs provide a controlled assessment of performance and robustness, showing that the proposed method substantially reduces range bias and improves seafloor position recovery relative to constant-sound-speed and single-gradient baselines, while remaining stable under SSP mismatch.
We further present an underwater acoustic transponder prototype integrating a chip-scale atomic clock (CSAC) and an FPGA-based multi-channel parallel clock disciplining subsystem.Sea trials in the South China Sea validate the end-to-end design and demonstrate representative ranging results, confirming kilometer-scale capability and stable real-time performance under in situ conditions. Overall, the proposed approach improves the fidelity of seafloor positioning time series and strengthens geodetic constraints on ilithospheric deformation and related earthquake hazard assessment.
How to cite: Liu, Z., Zhou, F., Zhao, J., Ji, X., and Cao, C.: Equivalent-Gradient Sound-Speed Correction and Joint Time-Bias Estimation for Stratified-Ocean Acoustic Ranging in Seafloor Geodesy, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-3939, https://doi.org/10.5194/egusphere-egu26-3939, 2026.
About 80% of shallow (depth < 60 km) earthquakes with Mw > 6.5 worldwide occur on offshore faults, highlighting the need for seafloor geodetic measurements to monitor strain accumulation before, during, and after major events. The current state-of-the-art technique for measuring horizontal seafloor deformation is GNSS/Acoustic (GNSS/A), which provides episodic measurements of absolute seafloor motion with centimeter-level precision at any depth. However, widespread application of GNSS/A remains limited by three main constraints: (1) high operational cost; (2) the inability to leave acoustic transponders in shallow waters (< 500 m) because of trawling activity; and (3) the need for long acquisition sessions to average out poorly modeled sound-speed variability in the water column. Here we present a new approach suitable for shallow water (< 300 m) that potentially enables centimeter-level seafloor geodesy at reduced cost and with shorter acquisition times. The method combines high-resolution optical imaging of the seafloor acquired by low-cost autonomous underwater vehicles (AUVs) with GNSS/A surveys. Acoustic beacons are used as ground control points, analogous to aerial photogrammetry, allowing georeferencing of the optical mosaics in a global reference frame. Natural markers such as rocks, reefs, and outcrops can then be re-imaged over time to measure displacement. Compared to classical GNSS/A, this approach uses acoustic beacons only during the survey, enabling multiple seafloor points to be monitored within a single experiment using a limited number of transponders. We will present results from a proof-of-concept experiment conducted in autumn 2025 near Toulon, southern France.
How to cite: Reveneau, H., Nocquet, J.-M., Royer, J.-Y., Furst, S., Varais, N., Verdun, J., Coulombier, T., Ballu, V., Eceiza, A., Sladen, A., and Sakic, P.: Hybrid Optical–GNSS/Acoustic Method for Centimeter-Precision Seafloor Geodesy in Shallow Water, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-10362, https://doi.org/10.5194/egusphere-egu26-10362, 2026.
Posters on site: Thu, 7 May, 16:15–18:00 | Hall X2
Satellite geodesy has become a powerful and reliable approach for investigating of geophysical processes, including tectonic, seismic, and volcanic activity. In seismotectonic research, the Global Navigation Satellite System (GNSS) provides precise estimates of plate motion parameters and offers valuable insights into present-day crustal deformation. The determination of Euler pole parameters and plate velocities represents a fundamental problem in global tectonics and a critical component of many geodynamic studies. This study aims to determine and analyze the parameters that characterize tectonic plate motion, specifically the rotation pole position (Φ, λ) and the angular rotation rate (Ω). The present-day kinematics of major tectonic plates are investigated using hundreds of geodetic observations provided by the Nevada Geodetic Laboratory (NGL) over the period 2000–2025. The analysis is based on GNSS station positions and velocities expressed in the International Terrestrial Reference Frame 2020 (ITRF2020). The estimation of Euler pole parameters for five major tectonic plates resulted in the following preliminary findings: for the African plate (𝜆𝐴𝑓=−81.97°±0.37; 𝜑𝐴𝑓=50.08°±0.14 and 𝛺𝐴𝑓=0.2665°𝑀𝑦𝑟⁄±0.0010);for the Eurasian plate (𝜆𝐸𝑢=−99.51°±0.98; 𝜑𝐸𝑢=54.61°±0.62 and 𝛺𝐸𝑢=0.2581°𝑀𝑦𝑟⁄±0.0016); for the North American plate (𝜆𝑁𝐴=−88.70°±0.32; 𝜑𝑁𝐴=−8.18°±0.42 and 𝛺𝑁𝐴=0.1863°𝑀𝑦𝑟⁄±0.0009 ); for South American plate (𝜆𝑆𝐴=−128.17°±0.92; 𝜑𝑆𝐴=−19.12°±0.35 and 𝛺𝑆𝐴=0.1178°𝑀𝑦𝑟⁄±0.0015) ; and for the Australian plate (𝜆𝐴𝑢=37.87°±0.22; 𝜑𝐴𝑢=32.84°±0.15 and 𝛺𝐴𝑢=0.6331°𝑀𝑦𝑟⁄±0.0007).The results were compared with the NNR-MORVEL56 plate motion model, providing key information for geodynamic modeling and insights into present-day tectonic behavior and seismic hazard.
Keywords
GNSS Stations; Tectonic Plates Velocities; Euler pole Parameters; ITRF2020.
How to cite: Allal, S. H., Hasni, K., and Dekkiche, H.: Recent Kinematics Analysis of Major Tectonic Plates from GNSS Positions and Velocities in the ITRF2020, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-10128, https://doi.org/10.5194/egusphere-egu26-10128, 2026.
North-Eastern Italy is of particular interest in tectonics as it lies on the northernmost edge of the convergent margin between Eurasia and the Adria microplate, influencing regional deformation and seismicity. The Friuli Venezia Giulia Deformation Network (FReDNet) was established in the area in 2002 to monitor crustal deformation and contribute to regional seismic hazard assessment.
Tunini et al. (2024) described the time series spanning two decades from GNSS stations located in north-eastern Italy and surrounding areas, as well as the resulting velocity field in the ITRF14 reference frame. The documented dataset was obtained by processing GPS observations with GAMIT/GLOBK software version 10.71. The time series, estimated using the same procedure, are collected daily and stored as part of a long-term monitoring project, with annual updates of velocity solutions computed.
In this study, we present a multi-constellation solution and evaluate the differences compared to the GPS-only solution. We also update the solution to the new implementation of the International Terrestrial Reference Frame (ITRF2020). We then analyse the impact of the new reference system on the characteristics of the time series and the velocities. We also include additional stations that have recently become available in the solution.
We acknowledge the CINECA award under the ISCRA initiative, for the availability of high performance computing resources and support (IscraC IsCd5_MGNSS20).
Reference:
Tunini, L., Magrin, A., Rossi, G., and Zuliani, D.: Global Navigation Satellite System (GNSS) time series and velocities about a slowly convergent margin processed on high-performance computing (HPC) clusters: products and robustness evaluation, Earth Syst. Sci. Data, 16, 1083–1106, https://doi.org/10.5194/essd-16-1083-2024, 2024.
How to cite: Magrin, A., Tunini, L., Zuliani, D., and Rossi, G.: Adria-Eurasia collision front: Multi constellation GNSS data Processing in ITRF2020 reference frame, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-11657, https://doi.org/10.5194/egusphere-egu26-11657, 2026.
Ocean Bottom Pressure gauges (OBPs) are devices installed on the seafloor to continuously measure ocean bottom pressure. They are expected to detect vertical seafloor crustal deformation caused by transient tectonic events. However, OBP data contain pressure changes originating from various sources other than tectonic events, such as ocean tides, non-tidal fluctuations due to meteorological and oceanic variations, and instrumental drift. We need to remove these non-tectonic components to identify the pressure changes related to vertical crustal deformation. However, the timescale of non-tidal fluctuations (ocean noise) is similar to that of transient tectonic signals, such as slow slip events (SSEs), which makes it difficult to separate these components.
Otsuka et al. (2023) applied principal component analysis (PCA) to seafloor pressure data obtained by an OBP array to detect transient events, assuming that pressure changes due to transient events cause temporal fluctuations in the PCA results. To verify the performance of this method, they applied PCA to synthetic data and confirmed that the method can detect temporal changes in the composition of principal components (PCs). In the present study, we apply this method to OBP data obtained before the 2011 Tohoku earthquake, which includes transient events resulting from aseismic slip, as reported by Ito et al. (2013), to verify whether we can identify the event through temporal changes in the PCs decomposed from the OBP data. We performed PCA on de-tided OBP data covering about four months before the Tohoku earthquake using a short sliding time window, to examine temporal variations in the PCs.
Changes in the PCs were evaluated using the normalized inner product (NIP) of the eigenvector of each PC (Otsuka et al., 2023), which measures the difference in the direction of the vector. We expect the NIPs to be stable if the OBP data do not contain any transient events, whereas evident changes in the NIPs of more than one PC would occur when transient pressure variations are included in the data. In the present study, the NIPs of PC1 (the most significant component) and PC2 (the second most significant component) remained stable over time, whereas the NIPs of PC3 and PC4 began to decrease as the moving window (60-day length) approached late January 2011 and remained low for about 40 days. Based on a comprehensive analysis of seismic and geodetic data, including the same OBP data, Ito et al. (2013) reported that an SSE lasted for about 40 days from late January 2011. The NIP changes detected in the present study may correspond to pressure changes due to this transient event.
How to cite: Yamada, T., Hino, R., Kubota, T., Otsuka, H., and Ohta, Y.: An Attempt to Detect Transient Crustal Deformation from Ocean Bottom Pressure Gauge Data Using Principal Component Analysis, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-3717, https://doi.org/10.5194/egusphere-egu26-3717, 2026.
Precise seafloor geodetic measurements are essential for understanding plate tectonics, earthquake cycles, and volcanic activity. Seafloor pressure gauges provide a key tool for monitoring vertical changes in the seafloor elevation. However, distinguishing millimeter-scale tectonic signals from instrumental drift and environmental noise remains a fundamental challenge in deep-ocean observing.
This presentation evaluates the performance and data integrity of self calibrating bottom pressure recorders (BPR) deployed by Ocean Networks Canada (ONC) offshore Vancouver Island. Instruments like the RBR BPRZero and the Sonardyne FETCH AZA utilize Ambient-Zero-Ambient (AZA) in-situ calibration mechanisms to quantify sensor drift. The AZA method involves switching a pressure gauge from ambient (seafloor) pressure to atmospheric pressure within the instrument's housing. By comparing this internal pressure reading to an accurate barometer also measuring internal pressure, the drift can be precisely determined and a calibration function applied.
A forensic analysis of the dataset reveals that, while the high-resolution pressure measurements capture true genuine environmental signals, they are also significantly contaminated by instrumental artifacts that are partially related to the measurement approach itself. In this study we report findings from preliminary deployments and discuss the methodological challenges encountered, proposing mitigation strategies for future applications.
How to cite: Schlesinger, A., Heesemann, M., Dexter, J., Leconte, J.-M., Davis, E., Sun, T., Kreimer, N., and Aghaei, O.: Assessing Data Integrity in Seafloor Geodesy: An Analysis of Self-Calibrating Pressure Data Collected by Ocean Networks Canada, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-1585, https://doi.org/10.5194/egusphere-egu26-1585, 2026.
Continuous observations of seafloor water pressure, which enable detection of vertical seafloor motion, are an essential technique in seafloor geodetic observations. For such geodetic applications, observation instruments equipped with Digiquartz sensors manufactured by Paroscientific Ltd. are commonly used because of their extremely high accuracy, resolution, and stability. Downsizing of observation instruments are important for making dense network allowing detection of subtle crustal deformation; however, the relatively high power consumption of Digiquartz makes it difficult to downsize self-pop-up mobile observation systems.
Our group has used RBRduo3, temperature–pressure loggers as auxiliary instruments to monitor environmental conditions in other seafloor observations, such as seafloor acoustic ranging. Given that these products are compact, lightweight, and capable of easily acquiring continuous records over periods exceeding one year, we investigated their potential applicability to seafloor geodetic observations. To this end, we conducted parallel observations with ocean-bottom pressure recorders (OBPRs) equipped with high-precision Digiquartz sensors and RBR loggers.
Pressure records from RBR loggers exhibit large transient variations immediately after deployment on the seafloor, with amplitudes up to about an order of magnitude larger than the typical transients observed in Digiquartz sensors. However, except for approximately the first three days after installation, this behavior can be well approximated by a time-dependent function combining an exponential term and a linear term, and no other irregular fluctuations are observed.
Results from parallel observations at sites with water depths exceeding 5,000 m show that, aside from the initial post-deployment transients, pressure time series obtained by RBR loggers agree well with those from Digiquartz sensors. In contrast, at shallower sites with water depths less than about 2,000 m, the pressure time series from the two instruments differ substantially. The time series of the pressure differences closely resembles the temperature time series, suggesting that these discrepancies arise from insufficient temperature correction of the RBR pressure data.
Assuming that apparent pressure variations caused by temperature changes dominate the short-period components of the pressure fluctuations recorded by the RBR logger, we estimated a coefficient for temperature correction by minimizing the power of the fluctuations. Using this coefficient, we removed the temperature-correlated component over the entire frequency band. As a result, we obtained pressure time series that agree with those from the Digiquartz sensors within approximately 0.2 hPa in the parallel observations. This demonstrates that, with appropriate temperature correction, it is possible to obtain pressure variation data from RBR loggers that are comparable in quality to those from Digiquartz sensors.
By taking advantage of their compact size and low power consumption, RBR loggers could be applied as add-on instruments to ocean bottom seismometers and similar seafloor observation instruments, and are expected to contribute to an increased number of observation points in mobile seafloor observation networks.
How to cite: Hino, R., Suzuki, S., Sato, M., Yamada, T., and Ohta, Y.: Quality assessment of seafloor pressure data from RBR loggers for applications to seafloor geodesy, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-6150, https://doi.org/10.5194/egusphere-egu26-6150, 2026.
The Cascadia Subduction Zone (CSZ) is a convergent margin extending from Northern California to Northern Vancouver Island that is capable of creating M9+ earthquakes. Direct seafloor geodetic observations near the deformation front and the locked zone are required to constrain possible rupture scenarios. Until recently, observations of plate motion have been limited to land-based stations, which are insufficient to resolve the variability of potential rupture scenarios. Consequently, different model scenarios yield varied outcomes for hazard assessments and the development of effective mitigation strategies. Offshore Oregon and Washington (Central Cascadia), this critical data gap is being addressed by GNSS-Acoustic (GNSS-A) observations that recently provided seafloor geodetic evidence indicating near-full locking on the shallow megathrust (DeSanto et al., 2025).
To the north, offshore Vancouver Island, the Northern Cascadia Subduction Zone Observatory (NCSZO) project, primarily funded by the Canada Foundation for Innovation (CFI) and operated by Ocean Networks Canada (ONC), provides seafloor geodetic observations that will constrain model scenarios and offer the opportunity to observe along-strike variations. The NCSZO complements ONC’s NEPTUNE cabled seafloor observatory, which provides real-time data from seismometers, bottom pressure recorders, CORK borehole observatories, and other sensors relevant to seafloor geodesy. The NCSZO is composed of two main offshore components: a GNSS-Acoustic (GNSS-A) seafloor geodesy network with seven stations and a Deformation Front Laboratory providing pressure and tilt measurements across the deformation front. This presentation will provide an overview of the NCSZO and will highlight the first results from the GNSS-A stations, which consist of accurately located seafloor benchmarks. These benchmark locations are monitored via tens of thousands of acoustic interrogations during yearly surveys utilizing an autonomous Wave Glider.
Following the completion of a fourth observation campaign in 2025, we start to see meaningful results from several sites that provide constraints on the locking of the Juan de Fuca plate with the overriding North American plate. These first direct GNSS-A measurements in Northern Cascadia are a significant step towards improving the reliability of regional earthquake and tsunami hazard and mitigation models.
How to cite: Heesemann, M., Hutchinson, J., Sun, T., Wang, K., Davis, E., Bailly, N., Schlesinger, A., and Trenaman, F.: Near-Trench Seafloor Geodesy: GNSS-Acoustic Plate Motion Measurements at the Northern Cascadia Subduction Zone, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-13992, https://doi.org/10.5194/egusphere-egu26-13992, 2026.
Ocean-bottom pressure (OBP) gauges are a key tool in seafloor geodesy, providing time-series observations of vertical seafloor motion by continuously recording pressure variations associated with changes in water depth. However, OBP time series contain multiple superimposed signals in addition to crustal deformation, including tides, instrumental drift, and non-tidal oceanic fluctuations driven by atmospheric and oceanic processes. In particular, non-tidal oceanographic fluctuations have time scales comparable to long-term crustal deformation, making it difficult to separate them and accurately estimate long-term seafloor deformation.
Ocean-model-based corrections have been proposed to isolate and remove non-tidal oceanographic components, offering a physically based approach. However, model–observation mismatches can leave large residuals after correction, making it harder to identify deformation-related signals. As a fundamental step toward understanding these model–observation discrepancies, this study quantifies the physical processes controlling OBP variability using an ocean general circulation model named OFES2 (Sasaki et al., 2020), which does not assimilate oceanic observations.
As a first step, we assess the extent to which OFES2 can account for the observed OBP variability by directly comparing modeled and observed time series. The modeled OBP is compared with observed OBP records from pressure gauges deployed off northeastern Japan. Our results show that OFES2 reproduces the seasonal cycle of the observed OBP time series to some extent, whereas agreement at periods of several days to about a month remains limited in both amplitude and phase.
To investigate the physical processes controlling the modeled OBP variability, we decompose the modeled OBP into four components: atmospheric pressure loading, sea surface height variability, horizontal density advection, and vertical advection, and quantify the relative contribution of each component. The decomposition indicates that atmospheric pressure loading and sea surface height variability dominate the modeled OBP fluctuations, while horizontal density advection contributes to part of the seasonal variability. Moreover, correlation analyses using the time-derivative of the observed OBP and the decomposed pressure components reveal episodic enhancements in correlation with the horizontal density advection term, suggesting that its contribution can temporarily increase during specific periods.
We will extend the analysis by comparing OFES2 with longer-term OBP observations and discussing in more detail the contributions of each component to seafloor pressure variability.
How to cite: Hagihara, T., Ohta, Y., Wang, S., Sasaki, Y., Hino, R., and Otsuka, H.: Process Decomposition of Ocean-Bottom Pressure Variability: What OFES2 Reproduces and Misses, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-8768, https://doi.org/10.5194/egusphere-egu26-8768, 2026.
The Cascadia subduction zone represents a seismic hazard to the Pacific Northwest region of North America, yet the state of fault locking near the deformation front, which could cause a devastating tsunami upon rupturing, remains poorly understood due to limited offshore observations along the subduction zone. In this study, we present the first seafloor geodetic measurements of the horizontal deformation rates on the accretionary prism from an array of four Global Navigation Satellite System-Acoustic (GNSS-Acoustic) sites surveyed from 2016-2022. These GNSS-Acoustic sites, despite resting on the North American plate, show velocities that are a significant fraction of the subducting Juan de Fuca plate velocity. In contrast, the continuous GNSS stations along the Oregon coast are moving at velocities <1 cm/yr relative to the North American Plate. Locking models constrained by these offshore velocities show that the subduction zone interface near the deformation front must be nearly locked offshore Oregon. To satisfy both the onshore and offshore geodetic observations, the locked zone must be relatively narrow and only minimal aseismic creep is permissible at the deformation front. These results suggest that appreciable elastic strain has accumulated near the deformation front, which elevates the potential for tsunamigenesis along this portion of the subduction zone.
How to cite: DeSanto, J., Schmidt, D., Zumberge, M., Sasagawa, G., and Chadwell, C. D.: Near full locking along the shallow megathrust of the Cascadia subdduction zone identified from seven years of GNSS-Acoustic observations, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-11756, https://doi.org/10.5194/egusphere-egu26-11756, 2026.
The Main Marmara Fault (MMF) keeps unruptured since the last Marmara earthquake in 1766, while the most part of its extension along the North Anatolian Fault has released its strain through M7 class earthquakes progressively occurred during the last century. Revealing the coupling condition MMF is quite important because it controls magnitude of the next earthquake expected to be occurred.
Since the MMF is covered with sea water, acoustic ranging technique is required to monitor the fault behavior in geodetic measure. So far, two experiments were carried out to measure the creep rate of the MMF using acoustic extensometers; no significant creep at the west of Central High (Sakic et al., 2016), while nearly a half of the strain is released by creep at the Western High (Yamamoto et al., 2019). We further conducted the experiment at the west edge of Kumburgaz Basin to fill the spatial gap of the two experiments.
The experiment started from May 2017 installing five extensometers across a valley formed by MMF activity. Acoustic ranging was successfully operated between nine pairs out of ten combinations. The longest one is up to 3 km and the shortest one is just 0.5 km. The observed data were extracted remotely from a ship via acoustic communication while the measurement continues. We already extracted the data for nearly two years from the beginning to April 2019, just before the interruption due to COVID-19. Change rates in baseline length are evaluated using roundtrip times, which are converted into distance using sound speed corrected for in situ temperature of sea water. Combining change rates and crossing angle of their baselines with MMF, nearly 10-15 mm/yr of right-lateral creep is expected at the site.
Obtained result taken together the past experiments, the west half of MMF is partially creeping and the east half is rocked, which indicates that MMF still has a potential for over M7 earthquake. This distribution creep is consistent with electromagnetic imaging beneath the Sea of Marmara (Kaya-Eken et al., 2025). Our extensometers were still working at least the time of the Mw6.2 Marmara earthquake in April 2025, which occurred just beneath our site as a series of eastward progressive intermediate ruptures (Martinez-Garzon et al., 2025). We are expecting further data analysis after 2019.
How to cite: Kido, M., Takahashi, N., Yamamoto, Y., Özener, H., and Kaneda, Y.: Creep of the Main Marmara Fault at the west of the Kumburgaz Basin observed by acoustic extensometers, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-21272, https://doi.org/10.5194/egusphere-egu26-21272, 2026.
Deformation measurements from space-borne Interferometric Synthetic Aperture Radar (InSAR) resolve geodynamic signals at multiple spatial scales and are commonly integrated with ground-based measurements from the Global Navigation Satellite System (GNSS). The integration process requires accounting for reference frame discrepancies and measurement uncertainties. In this study, we propose a new approach, CoInSAR, which uses least squares collocation to transform InSAR displacement time series into the GNSS reference frame while correcting residual tropospheric errors from turbulent atmospheric effects that typically persist in InSAR data. In our approach, we construct the observation vector at each epoch from the displacement differences between GNSS and InSAR, which comprises the reference frame difference, measurement noise, and residual tropospheric errors. We represent the measurement noise via data variances and derive the stochastic model for residual tropospheric errors using an empirical covariance function estimated from InSAR displacements. By accounting for these stochastic components, we use least squares collocation to estimate the transformation parameter between the two reference frames, interpolate corrections for the residual tropospheric errors, and generate the integrated displacements. To assess the performance of our method, we applied CoInSAR to measure land subsidence in the San Joaquin Valley, California, and seismic deformation in Chile. Our results show high agreement between GNSS and CoInSAR time series and variance reduction in regions outside the GNSS network. Moreover, CoInSAR-based deformation estimates are not only consistent with physics-based models but also capture small-scale deformation features, highlighting CoInSAR’s potential to improve the modeling of geodynamic signals in regions with sparse GNSS coverage.
How to cite: Wang, C.-Y., Gómez, D., and Figueroa, M.: A least squares collocation approach to integrate InSAR and GNSS observations: CoInSAR, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-12697, https://doi.org/10.5194/egusphere-egu26-12697, 2026.
While technologies such as satellite imaging are used in addition to seismic data for large magnitude earthquakes to confirm fault mechanisms, they are significantly less considered for the analysis of smaller sized events. However, constraining the nature of small to moderate magnitude seismic events with a sparse seismic network presents a significant challenge for monitoring agencies. The limitation in seismic coverage as well as low signal-to-noise levels measured in seismic data imply significant uncertainties in accurately estimating the source parameters (i.e., epicentral position, depth, magnitude and mechanism: earthquake, explosion, collapse). Quantifying and reducing these uncertainties becomes paramount, especially for the monitoring of shallow underground nuclear tests.
Here, we explore the advantages of combining seismic data and satellite borne Synthetic-Aperture-Radar-Interferometry (InSAR) techniques to recover the source parameters of superficial geophysical events using moment tensor techniques. We present the work we have undertaken for different types of shallow seismic events including earthquakes and collapses. From the acquisition of satellite images and seismic data to the comparison of the source solutions provided by both datasets, we explore the strengths and weaknesses of each approach that still need to be understood, and we propose joint approaches where possible. In the future, joint seismic-InSAR full moment tensor inversions may lead the way in the monitoring of regions with low seismic coverage.
How to cite: Guilhem Trilla, A., De Boever, H., Burgos, G., and Pinel-Puyssegur, B.: Improving Seismic Source Characterization of Moderate Magnitude Events at Regional Distance With Seismic and InSAR Data, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-18873, https://doi.org/10.5194/egusphere-egu26-18873, 2026.
Spectacular natural hazards such as large earthquakes or volcanic eruptions are accompanied by smaller-amplitude processes that produce millimeter-scale ground deformation. Although subtle, these signals provide critical insight into the physical state of the system and their associated hazards. Interferometric Synthetic Aperture Radar (InSAR) offers the spatial resolution required to observe such deformation, but its exploitation over long time spans remains challenging due to centimeter-scale noise and systematic biases. Here, we demonstrate how state-of-the-art interferometric processing combined with a Kalman Filter–based Time Series analysis (KFTS) enables the extraction of millimeter-scale deformation from 10 years of Sentinel-1A/B data.
We present two case studies: the Chaman fault system (Pakistan–Afghanistan) and the Natron rift (northern Tanzania) in the East African Rift. Careful step-by-step corrections of the interferograms include tropospheric and ionospheric corrections, azimuth shift compensation, and rigorous assessment of closure phase biases. Measurement uncertainties derived from coherence are propagated within the KFTS time-series inversion, allowing iterative estimation of phase evolution with associated uncertainties.
In the Chaman fault zone, we detect aseismic deformation characterized by fault creep rates of about 5 mm/yr, as well as a slow-slip event with ~1 cm of cumulative displacement resolved using combined ascending and descending geometries. In northern Tanzania, we resolve long-term rift opening of only a few millimeters per year between 2015 and 2025, consistent with GNSS campaign measurements. We further assess the potential of Independent Component Analysis (ICA) for InSAR signal separation and discuss current limitations imposed by residual noise and vegetation-related biases.
How to cite: Dalaison, M., Drique, L., Jolivet, R., Grandin, R., Pinel-Puysségur, B., Navarrete, I., Calais, E., and Olive, J.-A.: Resolving millimeter-scale tectonic deformation in decade-long InSAR time series : from long-term rifting to slow-slip events, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-17353, https://doi.org/10.5194/egusphere-egu26-17353, 2026.
Monitoring Volcanic Deformation Using InSAR: An Optimised InSAR Time-Series Approach for Seasonally Snow-Covered Volcanoes
1Tianyuan Zhu*, 1Juliet Biggs, 1Alison Rust, 2Milan Lazecký, 3Loreto Cordova
1School of Earth Sciences, University of Bristol, Bristol, United Kingdom
2COMET, School of Earth and Environment, University of Leeds, Leeds, United Kingdom
3Servicio Nacional de Geología y Minería (SERNAGEOMIN), Santiago, Chile
*tianyuan.zhu@bristol.ac.uk
Satellite-based Interferometric Synthetic Aperture Radar (InSAR) has been widely used for monitoring volcanic deformation, especially since Sentinel-1 launched in 2014, providing an unprecedented volume of routinely acquired, open-access data. Automated systems now continuously process interferograms and regularly update deformation time series, providing a valuable dataset for monitoring volcanoes globally. However, seasonal snow leads to coherence loss and subsequent unwrapping errors in interferograms, causing gaps in the network of the automated time-series analysis and reducing deformation accuracy. As ~41% of subaerial Holocene volcanoes exhibit seasonal snow cover (with snow persistence of 7-90%), optimising InSAR processing for seasonally snow-covered volcanoes would substantially improve monitoring active volcanoes, especially in high-latitude and high-altitude areas.
In this study, we developed an optimised InSAR time-series processing workflow using MODIS 8-Day Snow Product, which has been successfully applied to Laguna del Maule (LdM), a caldera with strong seasonal snow cover in Chile (Snow Persistence=51%). At LdM, the default product from the LiCSBAS auto-processing system underestimates the average line-of-sight deformation by 28% at GNSS station MAU2 between 10/2014 and 06/2023. To improve the accuracy of time series, we adapt the LiCSBAS time-series processing strategy using the quantified relationship between MODIS 8-Day Snow Products and Sentinel-1 InSAR coherences. The optimised workflow, including an algorithm based on Graph Theory for network selection, reduced data requirements by ~90% and LiCSBAS processing time by ~80%, while improving the accuracy of the LiCSBAS-processed deformation to match GNSS observations.
Vegetation is another crucial factor in coherence loss, and cloud cover affects optical satellite data. Using MODIS products, we also show that over 50 of 484 seasonally snow-covered volcanoes have lower Normalized Difference Vegetation Index (NDVI) and cloud-obscured duration than LdM (NDVI=0.16; Cloud Duration=144 days), confirming that seasonal snow is their dominant source of coherence loss and MODIS products are applicable.
We applied our validated workflow to seasonally snow-covered volcanoes across a range of environments with different cloud and vegetation cover to produce a long-term deformation (2014–present) using Sentinel-1 data. The optimised workflow has implications for the accuracy and efficiency of global volcano monitoring, improving the quality of modelling and forecasting.
How to cite: Zhu, T.: Monitoring Volcanic Deformation Using InSAR: An Optimised InSAR Time-Series Approach for Seasonally Snow-Covered Volcanoes, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-14378, https://doi.org/10.5194/egusphere-egu26-14378, 2026.
Global Navigation Satellite System (GNSS) observations are widely used to monitor volcanic deformation through continuous measurements of surface displacements. Detecting subtle volcano-tectonic signals relies on stable and repeatable position estimates over long time periods, which can be difficult to achieve in complex atmospheric environments.
In geodetic GNSS processing, station coordinates are estimated jointly with tropospheric parameters, including zenith delays and, in many processing strategies, horizontal tropospheric gradients. These parameters are known to be coupled with the vertical component of the position and to reflect both large-scale and locally anisotropic atmospheric conditions. At volcanic summits, tropospheric variability is often enhanced by strong orographic effects, frequent cloud formation, intense precipitation, and, in tropical regions, persistently high atmospheric humidity.
In this study, we investigate the behaviour and temporal variability of tropospheric parameters estimated for GNSS stations installed near the summits of La Soufrière de Guadeloupe and Montagne Pelée in Martinique, two active volcanoes located in the tropical Lesser Antilles. The analysis focuses on zenith tropospheric delays, horizontal gradients, and their consistency over time, as well as on diagnostic indicators derived from GNSS phase residuals.
GNSS data are processed using two independent geodetic software packages, GINS and GipsyX, enabling a comparative assessment of tropospheric estimates and residual patterns obtained under different processing strategies. The GNSS-derived tropospheric parameters are examined in conjunction with observations from nearby meteorological stations and with the ECMWF ERA5 reanalysis, providing an external reference for the observed atmospheric variability.
This work presents an initial investigation of tropospheric modelling at volcanic summits in tropical environments. It discusses possible implications for the stability and interpretation of GNSS position time series used in volcanic deformation monitoring.
How to cite: Sakic, P. and Nahmani, S.: Investigation of the effect of tropospheric delay on the quality of GNSS time series in a volcanological context, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-20324, https://doi.org/10.5194/egusphere-egu26-20324, 2026.
Steep relief and nearby ocean influenced weather patterns at Teide–Pico Viejo volcano create strong, elevation‑correlated tropospheric delays that can obscure millimetric ground motion in C‑band InSAR. We applied a robust workflow to isolate true deformation at Teide (2015–2025) by combining multi‑geometry Sentinel‑1 time series. Our pipeline integrates (i) ascending/descending TOPS stacks processed with SBAS method; (ii) numerical weather model corrections (e.g., ERA5‑based slant delays using PyAPS) to remove water‑vapor structure; and (iii) Common-mode filtering to minimize atmospheric residuals due a single deformation reference area (Mohammadnia et al., 2025). Finally, we recovered 3‑D deformation fields by combining information from 3 line‑of‑sight geometries. The atmospheric-corrected solutions substantially suppress topography‑correlated variance and reveal coherent, low‑amplitude deformation that otherwise would have been misinterpreted as larger magnitude volcanic deformation. The emergent pattern is dominated by slow, millimeters per year, deformation since 2022. Signals are superimposed with centimetric seasonal vertical signals. Our results demonstrate that rigorous atmospheric corrections are essential to recover sub-centimeter deformation spanning multiple years at high‑relief volcanoes to isolate magmatic-hydrothermal pressurization signals.
References:
Mohammadnia, M., Yip, M.W., Webb, A.A.G., González, P.J. (2025) Spontaneous transient summit uplift at Taftan volcano (Makran subduction arc) imaged using an InSAR common-mode filtering method, Geophysical Research Letters, doi:10.1029/2025GL114853
Acknowledgements: We thank Spanish Agencia Estatal de Investigación project PID2022-139159NB-I00 (Volca-Motion) funded by MCIN/AEI/10.13039/501100011033 and “FEDER Una manera de hacer Europa”.
How to cite: Gonzalez, P. J., Mohammadnia, M., Boulesteix, T., Webb, A. A. G., and Charco, M.: Atmospherically corrected deformation at a prominent‑topography volcano: Teide-Pico Viejo volcano, Tenerife, Canary Islands (2015–2025), EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-10208, https://doi.org/10.5194/egusphere-egu26-10208, 2026.
