Global Navigation Satellite System Network Research Articles

The Postseismic Deformation of the 6 July 2019 Mw 7.1 Ridgecrest Earthquake from Burst Overlap Interferometry, InSAR, and GNSS

Abstract The July 2019 Ridgecrest earthquake sequence consists of an Mw 6.4 foreshock and an Mw 7.1 mainshock, which ruptured a complex set of orthogonal faults in the eastern California shear zone. We measure the co- and postseismic deformation associated with this sequence using the Burst Overlap Interferometry (BOI) technique in addition to the commonly used Interferometric Synthetic Aperture Radar (InSAR). The BOI technique, which provides displacement in the satellite’s along-track direction, offers important information on the postseismic deformation that cannot be measured by traditional InSAR and is only sparsely measured by the Global Navigation Satellite System networks. The BOI data reveal up to 4 cm displacement in the along-track direction, 10 km north of the northern tip of the seismic rupture, and up to 3 cm displacement along the coseismically active faults. These results rule out the possibility of significant shallow afterslip near the mainshock hypocenter, suggesting that the previously suggested poroelastic rebound is likely to be the cause for the significant postseismic line of sight deformation near the mainshock rupture. Based on the aftershocks’ moment tensor distribution, surface rupture, and simple forward modeling, we propose that the postseismic deformation north of the Ridgecrest rupture is caused by an aseismic slip along a north-trending normal fault of the Ridgecrest rupture that was induced by the Ridgecrest earthquake. Furthermore, we observed that both deformation and seismic activity decay slower over time as the distance from the Coso geothermal area increases. This decay is influenced by the mechanical properties of the crust, which are affected by the increased heat flow at Coso and thus suppress deformation and seismicity, ultimately controlling their temporal evolution.

Probing primordial black holes and dark matter clumps in the Solar System with gravimeter and Global Navigation Satellite Systems networks

Probing primordial black holes and dark matter clumps in the Solar System with gravimeter and Global Navigation Satellite Systems networks

Preliminary Results on Deformations of the Central North Anatolian Fault Zone

The North Anatolian Fault Zone (NAFZ) in the north of Türkiye is considered one of the most active faults on earth and have caused many major earthquakes in the last century. Surface ruptures by those earthquakes indicate that this tectonic feature consists of many segments, which reveal different characteristics in terms of geometry and mechanism. NAFZ lies along between Bingöl in the eastern Turkey and Saros Gulf in the west. This study focuses on the central part of the NAFZ between Niksar and Ilgaz provinces. In addition to the main branch, the central section of the NAFZ also contains two sub-branches extending southward: Merzifon-Esençay and Sungurlu Faults. In the literature, there are several studies for the region, however a comprehensive study mainly focusing on the main branch of the NAFZ and its sub-branches along the central NAFZ has not been conducted yet, which would fully reveal the fault mechanism and seismic potential in the region. The aim of this study is to derive an up-to-date high-spa-tial-resolution Global Navigation Satellite Systems (GNSS) velocity field for the central NAFZ and to correlate with the strain accumulations along the faults. Based on this, a high-spatial-resolution GNSS network consisting of 60 sites was established enclosing the central NAFZ. The previously archived GNSS data for the same observation sites were collected from the different studies in the region and a new GNSS campaign measurement was conducted between October 2023 and February 2024. All GNSS data were processed using GAMIT/GLOBK software, and the results indicate the GNSS velocities ranging from 3-30 mm/year and their associated uncertainties of maximum ± 1.5 mm/year with respect to the Eurasian tectonic plate.

Empirical orthogonal function based modelling of ionosphere using Turkish GNSS network

The study investigates ionospheric total electron content (TEC) variation over Turkey from the five selected global navigation satellite system (GNSS) stations situated in diverse parts of Turkey. The geomagnetic indices are used and observed TEC are modeled with the technique known as Empirical orthogonal function (EOF). It is valuable to note that the correlation coefficient between observed GNSS TEC values and EOF TEC values varies from 0.8020 to 0.9394. The root means square error (RMSE) values between observed GNSS TEC values and EOF TEC lie between 3.1665 TECU to 4.4220 TECU. These results show that the EOF model performs quite well in the Turkish region and can present the model TEC variations perfectly. Finally, these GNSS observed and EOF-predicted TEC values along with geomagnetic indices are studied with the tropospheric wind speed. The results showed that both observed and modeled TEC have very low correlations with tropospheric wind speed and do not provide any significant value. Hence, we concluded that the ionospheric region is not affected by tropospheric wind speed. It happens because the tropospheric wind speed is a matter of the lower troposphere and its atmospheric pressure while the ionosphere is far from the earth and depends upon the number of free electrons.

Enhancing Atmospheric Monitoring Capabilities: A Comparison of Low- and High-Cost GNSS Networks for Tropospheric Estimations

Global Navigation Satellite System (GNSS) signals experience delays when passing through the atmosphere due to the presence of free electrons in the ionosphere and air density in the non-ionized part of the atmosphere, known as the troposphere. The Precise Point Positioning (PPP) technique demonstrates highly accurate positioning along with Zenith Tropospheric Delay (ZTD) estimation. ZTD estimation is valuable for various applications including climate modelling and determining atmospheric water vapor. Current GNSS network resolutions are not completely sufficient for the scale of a few kilometres that regional climate and weather models are increasingly adopting. The Centipede-RTK network is a low-cost option for increasing the spatial resolution of tropospheric monitoring. This study is motivated by the question of whether low-cost GNSS networks can provide a viable alternative without compromising data quality or precision. This study compares the performance of the low-cost Centipede-RTK network in calculating the Zenith Tropospheric Delay (ZTD) to that of the existing EUREF Permanent Network (EPN), using two alternative software packages, RTKLIB demo5 version and CSRS-PPP version 3, to ensure robustness and software independence in the findings. This investigation indicated that the ZTD estimations from both networks are almost identical when processed by the CSRS-PPP software, with the highest mean difference being less than 3.5 cm, confirming that networks such as Centipede-RTK could be a reliable option for dense precise atmospheric monitoring. Furthermore, this study revealed that the Centipede-RTK network, when processed using CSRS-PPP, provides ZTD estimations that are very similar and consistent with the EUREF ZTD product values. These findings suggest that low-cost GNSS networks like Centipede-RTK are viable for enhancing network density, thus improving the spatial resolution of tropospheric monitoring and potentially enriching climate modelling and weather prediction capabilities, paving the way for broader application and research in GNSS meteorology.

GSeisRT: A Continental BDS/GNSS Point Positioning Engine for Wide-Area Seismic Monitoring in Real Time

Precise coseismic displacements in earthquake/tsunamic early warning are necessary to characterize earthquakes in real time in order to enable decision-makers to issue alerts for public safety. Real-time global navigation satellite systems (GNSSs) have been a valuable tool in monitoring seismic motions, allowing permanent displacement computation to be unambiguously achieved. As a valuable tool presented to the seismic community, the GSeisRT software developed by Wuhan University (China) can realize multi-GNSS precise point positioning with ambiguity resolution (PPP-AR) and achieve centimeter-level to sub-centimeter-level precision in real time. While the stable maintenance of a global precise point positioning (PPP) service is challenging, this software is capable of estimating satellite clocks and phase biases in real time using a regional GNSS network. This capability makes GSeisRT especially suitable for proprietary GNSS networks and, more importantly, the highest possible positioning precision and reliability can be obtained. According to real-time results from the Network of the Americas, the mean root mean square (RMS) errors of kinematic PPP-AR over a 24 h span are as low as 1.2, 1.3, and 3.0 cm in the east, north, and up components, respectively. Within the few minutes that span a typical seismic event, a horizontal displacement precision of 4 mm can be achieved. The positioning precision of the GSeisRT regional PPP/PPP-AR is 30%–40% higher than that of the global PPP/PPP-AR. Since 2019, GSeisRT has successfully recorded the static, dynamic, and peak ground displacements for the 2020 Oaxaca, Mexico moment magnitude (Mw) 7.4 event; the 2020 Lone Pine, California Mw 5.8 event; and the 2021 Qinghai, China Mw 7.3 event in real time. The resulting immediate magnitude estimates have an error of around 0.1 only. The GSeisRT software is open to the scientific community and has been applied by the China Earthquake Networks Center, the EarthScope Consortium of the United States, the National Seismological Center of Chile, Institute of Geological and Nuclear Sciences Limited (GNS Science Te Pū Ao) of New Zealand, and the Geospatial Information Agency of Indonesia.

Stability Analysis of GNSS Stations Affected by Samos Earthquake

An earthquake cycle can cause meters of displacement on the surface and at Global Navigation Satellite System (GNSS) stations. This study focuses on the identification of GNSS stations that have significant displacement because of a Mw 7.0 earthquake near Samos Island on 30 October 2020. The S-transformation method is used to examine 3D, 2D and 1D coordinate systems along with threshold and statistical test approaches. The highest coseismic offset among the 21 GNSS stations is displayed by SAMO, and CESM, MNTS, IZMI and IKAR also experience significant displacement. Significantly displaced stations are successfully identified in both 3D and 2D analyses. In the up component, SAMO is the only unstable station. The coordinate S-transformation method can be used in detecting unstable points in a GNSS network and provide valuable information about the effects of an earthquake on GNSS stations.

THE EFFECT OF SEASONAL VARIATION ON GNSS ZENITH TROPOSPHERIC DELAY

Abstract. The Global Navigation Satellite System (GNSS) signal experiences delays caused by the atmosphere, leading to the lengthening of the geometric path of the ray, commonly referred to as tropospheric delay. This delay is a significant source of error in GNSS positioning, contributing to a bias in the height component of several centimeters, even when meteorological data are simultaneously recorded and used in tropospheric models. In this study, considering seasonal variations, we investigated the impact of tropospheric delay on the GNSS height component. GNSS stations, part of the Turkish RTK CORS Network known as TUSAGA-Active (Turkish National Permanent GNSS Network Active), covered different heights over the 2014–2019 period. Daily coordinates of GNSS stations and tropospheric zenith delay were obtained through the GAMIT/GLOBK software solution.In the study, temperature, pressure, and relative humidity data of meteorological stations at different heights were converted to mean sea level. By using these values, interpolation estimates were made for the continuous GNSS stations in the same region with the IDW method. The most significant delay in GNSS signals occurs in July and August. This effect, which causes periodic changes in the zenith delay, varies inversely with the station's height. With the increase in the amount of water vapor in the atmosphere in parallel with the rise in the temperature in the summer months, it is seen that the stations at a low height are more exposed to the tropospheric effect than the stations at higher heights. In addition, GNSS stations' reduced meteorological values (temperature, pressure and relative humidity) show that the zenith delay values changed directly proportional to the temperature and inversely proportional to the pressure and relative humidity.

Drought‐Induced Vertical Displacements and Water Loss in the Po River Basin (Northern Italy) From GNSS Measurements

AbstractWe study vertical ground displacement time series from Global Navigation Satellite System (GNSS) stations to measure deformation associated with hydrological drought in the Po river basin. Focusing on interannual trend changes, rather than seasonal (annual) components, we found a clear spatially correlated deformation signal that is temporally (anti)correlated with changes in the Po river level and the SPEI‐12 drought index, with stations moving upward during periods of river/index level decrease and vice versa. In the 2021–2022 time span, which culminated in the most severe drought of the last two centuries, we estimate the amount and spatial distribution of water loss in the basin and its surroundings. Excluding the seasonal signals, between January 2021 and August 2022, the GNSS stations underwent uplift, up to 7 mm, which corresponds to ∼70 Gtons of water loss. Compared to Global Land Data Assimilation System and Gravity Recovery and Climate Experiment estimates, GNSS results show a similar temporal evolution of water content but a more heterogeneous distribution of values. We show that continuous GNSS networks provide an effective way to monitor multiannual trend changes in water storage even in small water basins and serve as a reliable indicator of drought severity.

Vertical Land Motion Due To Present‐Day Ice Loss From Greenland's and Canada's Peripheral Glaciers

AbstractGreenland's bedrock responds to ongoing ice loss with an elastic vertical land motion (VLM) that is measured by Greenland's Global Navigation Satellite System (GNSS) Network (GNET). The measured VLM also contains other contributions, including the long‐term viscoelastic response of the Earth to the deglaciation of the last glacial period. Greenland's ice sheet (GrIS) produces the most significant contribution to the total VLM. The contribution of peripheral glaciers (PGs) from both Greenland (GrPGs) and Arctic Canada (CanPGs) has not carefully been accounted for in previous GNSS analyses. This is a significant concern, since GNET stations are often closer to PGs than to the ice sheet. We find that, PGs produce significant elastic rebound, especially in North and East Greenland. Across these regions, the PGs produce up to 32% of the elastic rebound. For a few stations in the North, the VLM from PGs is larger than that due to the GrIS.

Earth Observation Capacity of the National Geoinformation Center of Bulgaria as Part of the Tools for Geo- and Anthropogenic-Hazard Management in EPOS

Exploring the nature of the Earth system and its multi-scale architecture requires an integrated symbiotic approach not only to understand our planet but also to reveal the hazards and footprints of their events on human life. Following the conceptual model of the European Plate Observing System (EPOS) we integrate a variety of research to obtain a multilateral use of geoinformation for scientific and practical purposes. In 2018, we started building a Bulgarian distributed scientific infrastructure called the National Geoinformation Center, aiming to provide data, products, and expertise to predict and prevent natural and anthropogenic risks and disasters, as well as environmental changes. The present paper demonstrates the capabilities of seismological, Global Navigation Satellite System (GNSS), geomagnetic and multi-scale laboratories (MSL) in Bulgaria that participate to the respective thematic core services of EPOS. More than 40 seismological stations register and transmit data in real time through the operational center of the Bulgarian Seismological Network. Data from 18 GNSS stations, part of the National GNSS network, are officially registered in the EPOS GNSS Data Gateway. Data and products are available via the EPOS GNSS Product Portal for interested users. Panagjurishte geomagnetic observatory (PAG) provides real time variations and long-period data series for the Earth’s magnetic field. The Palaeomagnetic Laboratory at the National Institute of Geophysics, Geodesy and Geography (NIGGG) participates in MSL with its unique facilities for paleo-, archeo- and environmental magnetism studies. Four case studies are presented to illustrate best practices for utilizing the data collected. The analysis of seismicity in a region of a long-time exploited salt deposit is upgraded with GNSS data to reveal the reason for crustal movements in the region. The results show that the majority of the seismic events took place in the upper part of the Earth’s crust at a depth of up to 5 km, where the salt body is buried. The observed significant local deformations of the Earth’s surface are very likely related to the technological process of deposit exploitation. A recent model of the geomagnetic declination over the Bulgarian territory is also presented. It tracks the variation of the regional magnetic field from 2015 to 2020 and adds short-wavelength signals from the secular network observations. Last but not least, the application of a mineral magnetic laboratory analysis along the depth of a Holocene soil profile is demonstrated. Results reveal the soil evolution in response to environmental changes adding a new magnetic proxy parameter for more precise identification and characterization of the soil formation processes.

Water Storage Variations Recovered from Global Navigation Satellite System Network Using Spatial Constraints: A Case Study of the Contiguous United States

Global Navigation Satellite System (GNSS) vertical displacements are widely used to infer terrestrial water storage (TWS) variations. The traditional Laplacian inversion requires dedicated efforts to determine the optimal parameters, which has an important effect on the spatial patterns. In this study, we develop a new GNSS inversion method with flexible spatial constraints. One major merit is that the new method only requires loose boundary conditions rather than optimal parameters. A closed-loop simulation shows that the inversion using spatial constraints is improved by 7–21% compared with the Laplacian constraints. We apply this method to 18 watersheds across the Contiguous United States (CONUS) to infer daily TWS variations from January 2018 to August 2022. The results show that the amplitudes of monthly TWS time series from the spatial and Laplacian constraints are comparable to the Gravity Recovery and Climate Experiment (GRACE) Follow-On (GFO) in 16 watersheds. Furthermore, the standard deviation between the spatial constraints and GFO is at the same level as that between the Laplacian constraints and GFO. We also extract the daily TWS variations caused by heavy precipitation events in California. Our results demonstrate that spatial constraint inversion supplements the existing constraint strategies of GNSS inversion in hydrogeodesy; therefore, spatial constraint inversion can be an alternative tool for GNSS inversion.

ADVANCED ANALYSIS TOOLS FOR THE EUROPEAN GROUND MOTION SERVICE DATA

Abstract. This paper is focused on Advanced Differential Interferometric SAR (A-DInSAR). In the first part, the paper describes the European Ground Motion Service (EGMS), a new service of the Copernicus Land Monitoring Service. The EGMS provides consistent, regular, standardized, harmonized and reliable information regarding natural and anthropogenic ground motion over Europe. It is derived from A-DInSAR analysis of Copernicus Sentinel-1 data. The EGMS includes three main products: Basic, i.e. basic deformation maps delivered for individual and consistent frames of image stacks; Calibrated, i.e. deformation map integrated into a reference Global Navigation Satellite System network; and Ortho, which includes the East-West and Up-Down deformation, obtained by combining the ascending and descending results. The second part of the paper describes examples of data analysis. The first one is a clustering of PSI results, which highlights the so-called Active Deformation Areas (ADAs). The second one is an automatic classification of the ADAs. The third analysis concerns the differential deformation maps, which are useful to highlight areas where high deformation gradients occur. Such gradients are often associated to damages in buildings and infrastructures.

IONOSPHERIC TEMPORAL-SPATIAL CORRELATION ANALYSIS USING GNSS NETWORKS OVER CHINA

Abstract. The Global Navigation Satellite System (GNSS) is well recognized as a valuable tool for studying ionospheric variation, providing a means to measure variations in the Total Electron Content (TEC) of the ionosphere. In this manuscript, we focus specifically on the use of GNSS networks to study ionospheric diurnal variation over China, highlighting some of the key findings and challenges in this area. We begin by describing the basic principles of GNSS-based TEC measurements. We then analyze ionospheric diurnal variation using about 320 GNSS stations in China, including the Southern, Middle, and Northern areas. We also discuss some of the challenges associated with GNSS-based TEC measurement and analysis, such as the effects of receiver and satellite biases and understanding of spatial–temporal TEC variation. Finally, we highlight the sunset enhancement and scintillation of GNSS-based ionospheric diurnal variation research. Overall, this abstract provides a valuable overview of the use of GNSS networks for studying ionospheric diurnal variation and its potential applications for improving space weather forecasting and mitigating the effects of ionospheric disturbances on GNSS-based navigation systems.

Optimizing Boundary Conditions in GNSS Tomography: A Continuous 7‐Month Case Study in the Amazon

AbstractDifferent estimates of the regional water vapor scale height, taken from ERA5 reanalysis, in situ observations and the direct optimization of retrieved water vapor profiles in GNSS tomography, are found to have major impact in the performance of tomographic inversions, with the better results displaying mean errors comparable to radiosondes. The analysis uses 7 months of GNSS (Global Navigation Satellite Systems) observations in the Amazon Dense GNSS Network near Manaus, Brazil, in 2011–2012, to compute a time series of water vapor profiles, with a tomographic technique capable of producing quasi‐instantaneous inversions with minimal external data or constraints. Results compare very well with 12‐hourly in situ radiosondes, especially in the lower troposphere above 2 km, and its daily‐to‐seasonal variability compares well with WRF (Weather Research and Forecasting) convective‐permitting simulations driven by ERA5 boundary conditions, suggesting that GNSS tomography may be an important source of atmospheric water vapor data for different applications.

A machine learning approach for slow slip event detection using GNSS time-series

Extracting tectonic transient displacements on the Earth’s surface from Global Navigation Satellite System (GNSS) time series remains a challenge, because GNSS station displacements depend on multiple processes occurring simultaneously, along with noise that obscures low-magnitude transient signals. We present a novel method for automatic detection of slow slip events (SSEs) in time series of a GNSS network by training a supervised machine learning (ML) model for classification. The proposed methodology detects both temporally and spatially the signatures of SSEs or regional transients within a GNSS network. The time series of a GNSS network were transformed into grayscale images, from which descriptors, including Bag of Visual Words (BoW) and Extended Local Binary Patterns (ELBP), were extracted. These descriptors served as input features for two distinct ML models: Support Vector Machines (SVM) and Artificial Neural Networks (NN). To train and test the ML classification model, two 3-year synthetic datasets were generated, one with GNSS networks featuring slow slip events (SSEs) of varying location, duration, onset time, and magnitude, and the other without SSEs, resulting in positive and negative sets, respectively. For each GNSS network, an image was created by combining the east and north components of the time series, which have been previously detrended and common mode error filtered. Each image is further divided into sub-images corresponding to 60 days time windows, in order to temporarily detect the existence of a transient. For training and testing, the datasets were separated into 75% for training and 25% for testing, each with 50% positive and 50% negative cases. In the final step, we analyze the positively classified images, representing the time windows in which the classifier detected transients. Within each of these windows, we identify the network’s time series with the highest velocity, indicating the stations and geographic area where the detected transients occurred. The test results demonstrate that both ML models achieved high performance using both ELBP and BoW descriptors as features. Finally, our ML models were validated on a real dataset with a transient signal recorded before the 2014 Iquique earthquake in Chile, and they effectively detected this anomalous signal. The proposed method can effectively detect transient signals related to SSEs with high accuracy, sensitivity, and specificity in both the test and instrumentally recorded datasets.

Active subduction and strain partitioning in western Myanmar revealed by a dense survey GNSS network

The active tectonics of the Myanmar (formerly Burma) microplate in central Myanmar are caused by the oblique collision between the Sunda and Eurasia plates to the east and the Indian plate to the west. The kinematics of this plate boundary zone remain a subject of controversy, due in part to limited geodetic and seismic observations in the region. Here, we report 41 new and 67 updated survey-mode and continuous Global Navigation Satellite System (GNSS) velocities spanning central and western Myanmar, between 15 – 23°N and 93 – 98°E. For the first time, our network allows us to directly observe the occurrence of active convergence along the central Rakhine coast in a manner consistent with active subduction and kinematic coupling of the megathrust, at a rate of 15 – 20 mm/yr. The data also establish a nearly constant strike-slip rate of 20 – 23 mm/yr along the entire length of the Sagaing fault. Using an elastic block-modeling approach, we examine several proposed models for the partitioning of oblique convergence into strike-slip and thrust motion across the Indo-Myanmar Ranges and the Myanmar Central Basin. We find evidence for shortening perpendicular to local topography throughout the Indo-Myanmar Ranges, indicating active convergence and nearly complete strain partitioning between thrust and transform faults north of 19°N, and estimate a rate of moment-deficit accumulation along the megathrust of 0.8−1.5×1019 N m/yr, equivalent to one MW 9.0 earthquake every 250 to 450 years, assuming all elastic strain is released seismically. We also find that a relatively small set of major active crustal faults can explain much of the data, though the data cannot resolve smaller faults that may be active at rates lower than a few mm/yr. Our results have important implications for earthquake and tsunami preparedness across the region, and indicate areas for future study where the geodetic observations remain ambiguous about ongoing fault activity and seismic hazard, particularly within the densely populated Ayeyarwady Delta region and along the southern Rakhine coast.

Imaging the Spatiotemporal Evolution of Plate Coupling With Interferometric Radar (InSAR) in the Hikurangi Subduction Zone

AbstractThe coupling at the interface between tectonic plates is a key geophysical parameter to capture the frictional locking across plate boundaries and provides a means to estimate where tectonic strain is accumulating through time. Here, we use both interferometric radar (InSAR) and Global Navigation Satellite System (GNSS) data to investigate the plate coupling of the Hikurangi subduction zone beneath the North Island of New Zealand, where multiple slow slip cycles are superimposed on the long‐term loading. We estimate the plate coupling across the subduction zone over three multi‐year observational periods targeting different stages of the slow slip cycle. Our results highlight the importance of the observational time period when interpreting coupling maps, emphasizing the temporal variability of plate coupling. Leveraging multiple geodetic data sets, we demonstrate how InSAR provides powerful constraints on the spatial resolution of both plate coupling and slow fault slip, even in a region where a dense GNSS network exists.

High‐Density Integrated GNSS and Hydrologic Monitoring Network for Short‐Scale Hydrogeodesy in High Mountain Watersheds

AbstractWe installed a purpose‐built network of co‐located Global Navigation Satellite System (GNSS) stations and meteorological instrumentation to investigate water storage in a high‐mountain watershed along the Idaho‐Montana border. Twelve GNSS stations are distributed across the Selway‐Lochsa watersheds at approximately 30–40 km spacing, filling a critical observational gap between localized point measurements and regional geodetic and satellite data sets. The unique coupling of geodetic and hydrologic observations in this network enables direct comparison between co‐located GNSS measurements of the elastic response of the solid Earth and local changes in measured water storage. This network is specifically designed to address questions of hydrologic storage and movement at the mountain watershed scale. Here, we describe technical details of the network and its deployment; introduce new hydrologic, meteorologic, and geodetic data sets recorded by the network; process and analyze the source data (e.g., time series of daily three‐dimensional GNSS site positions, removal of non‐hydrologic signals); and characterize basic empirical relationships between water storage, water movement, and GNSS‐inferred surface displacement. The network shows preliminary evidence for spatial differences in displacement resulting from a range of snow loads across elevations, but longer and more complete data records are needed to support these initial findings. We also provide examples of additional scientific applications of this network, including estimations of snow depth and snow water equivalent from GNSS multipath reflectometry. Finally, we consider the challenges, limitations, and opportunities of deploying GNSS and weather stations at high elevations with heavy snowpack and offer ideas for technical improvements.

GNSS visibility and performance implications for the GENESIS mission

The GENESIS mission prepared for launch in 2027 integrates the four space-geodetic techniques on a single spaceborne platform in medium Earth orbit. With its unique observations and alternative tie concepts, the mission aims to contribute to an improved accuracy and homogeneity of future terrestrial reference system realizations. To assess the expected contribution of Global Navigation Satellite System (GNSS) tracking, a comprehensive GNSS coverage analysis is performed based on detailed link-budget simulations, taking into account the best available gain patterns and signal-specific transmit power estimates derived for this work from measurements of a high-gain dish antenna. The benefit of different receiver antenna concepts for the GENESIS spacecraft is assessed, and it is demonstrated that a single-antenna system with either a nadir-looking or side-looking boresight is a viable alternative to the dual-antenna configuration considered in initial mission studies. Compared to terrestrial users and missions in low Earth orbit, GENESIS will collect GNSS signals transmitted at up to two times larger off-boresight angles. Only limited information on the actual transmit antenna phase patterns is presently available in this region, which hampers a quantitative assessment of the expected measurement and orbit determination accuracy. As such, a comprehensive release of manufacturer calibrations is encouraged for all blocks of GPS and Galileo satellites. In parallel, a need for in-flight characterization and calibration of the GNSS transmit antennas for off-boresight angles of up to 30^circ using observations of the GENESIS mission itself is expected. The impact of such calibrations on the overall quality of terrestrial reference frame parameters will need to be assessed in comprehensive simulations of global GNSS network solutions with joint processing of terrestrial and GENESIS GNSS observations.