Dense Global Navigation Satellite System Research Articles

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.

Constructing a Regional Ionospheric TEC Model in China with Empirical Orthogonal Function and Dense GNSS Observation

Using Global Navigation Satellite Systems (GNSS) observation data for developing a high-precision ionospheric Total Electron Content (TEC) model is one of the essential subjects in ionospheric physics research and the application of satellite navigation correction. In this study, we integrate the Empirical Orthogonal Function (EOF) method with the TEC data provided by the Center for Orbit Determination in Europe (CODE), and observed by the dense GNSS receivers operated by the Crustal Movement Observation Network of China (CMONOC) to construct a regional ionospheric TEC model over China. The EOF analysis of CODE TEC in China from 1998 to 2010 shows that the first-order EOF component accounts for 90.3813% of the total variation of the ionospheric TEC in China. Meanwhile, the average value of CODE TEC is consistent with the spatial and temporal distribution characteristics of the first-order EOF base function, which mainly reflects the latitude and diurnal variations of TEC in China. The first-order coefficient after EOF decomposition shows an obvious 11-year period and semi-annual variations. The maximum amplitude of semi-annual variation mainly appears in March and October, which is closely associated with the variation in geographical longitude, the semi-annual change of the low-latitude electric field, and the ionospheric fountain effect. The second-order coefficient has an evident annual variation, the minimum amplitude mainly occurs in March, August, and September, and the amplitude values in the high solar activity years are more significant than those in the low solar activity years. The third-order coefficient mainly shows the characteristics of annual variation, and the fourth-order coefficient shows the noticeable semi-annual and annual variations. The third and fourth-order coefficients are both modulated by the solar activity index F10.7. The ionospheric TEC model in China, driven by CMONOC real-time GNSS observation data, can better reflect the latitude, local time and seasonal variation characteristics of ionospheric TEC over China. In particular, it can clearly show the spring and autumn asymmetry of ionospheric TEC in the low latitudes. The root mean square error of the absolute error between the model and the actual observation is mainly distributed around 2.45 TECU (1 TECU = 1016 electrons/m2). The values of the TEC model constructed in this study are closer to the actual observed values than those of the CODE TEC in China.

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.

Investigating the ERA5-Based PWV Products and Identifying the Monsoon Active and Break Spells with Dense GNSS Sites in Guangxi, China

Precipitable water vapor (PWV) with high precision and high temporal resolution estimated by Global Navigation Satellite System (GNSS) is widely used in atmospheric research and weather forecasting. However, most previous works are not consensual concerning the characteristics of the PWV at different time scales and the identification of active and break spells during summ er monsoon climate in Guangxi, China. Taking radiosonde (RS) observations as reference, a strong correlation (R > 0.97) exists between GNSS PWV and RS PWV with a mean root mean square error (RMSE) of 2.68 mm. The annual, seasonal, monthly, and diurnal PWV variations of three years (2017, 2018 and 2020) over Guangxi in were comprehensively investigated using 104 GNSS stations and the fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) Atmospheric Reanalysis (ERA5). The mean annual bias and RMSE between GNSS PWV and ERA5 PWV are −1.04 mm and 2.63 mm, respectively. The monthly bias and RMSE range are −0.77 to 3.87 mm, 1.32 to 4.45 mm, and the daily range is −1.41 to 1.07 mm and 1.11 to 5.02 mm, respectively. Additionally, the adopted average standardized rainfall anomaly criteria also identified 7/7/3 active spells and 5/3/7 break spells during the summer monsoon (June–September) from 2017 to 2020, respectively. During the three-year period, the daily amplitude ranges for active spells varied from 1.41 to 2.49 mm, 0.69 to 5.4 mm, and 0.88 to 1.41 mm, while the ranges for break spells were 2.45 to 6.76 mm, 1.66 to 8.17 mm, and 1.48 to 2.99 mm, respectively. The results show a superior performance of GNSS PWV compared to ERA5 PWV in Guangxi, and the maximum, minimum and occurrence time of PWV anomaly vary slightly with the season and the topography of stations. Despite temperature primarily exhibiting a negative correlation with rainfall, acting as a dampener, a positive correlation remains evident between PWV and rainfall. Therefore, densely distributed GNSS stations exhibit excellent capabilities in quantifying atmospheric water vapor and facilitating real-time monitoring of small and medium-scale weather phenomena.

Characterization of the co-seismic pattern and slip distribution of the February 06, 2023, Kahramanmaraş (Turkey) earthquakes (Mw 7.7 and Mw 7.6) with a dense GNSS network

Two consecutive earthquakes with the magnitudes of Mw 7.7 and 7.6 (February 06, 2023) occurred on the East Anatolian Fault Zone (EAFZ) segments and unfortunately resulted in significant devastation to human life and cities in Turkey and Syria. In this study, we aimed to analyse the co-seismic displacements and fault slip distributions of these seismic events. Our unique high-spatial-resolution Global Navigation Satellite System (GNSS) network (comprising 73 permanent GNSS stations and 40 campaign observation sites), providing the recent geodetic dataset for the region, allows better constraint of the co-seismic surface displacements and slip distributions of both earthquakes. The three largest total displacements were identified as 466 cm, 362 cm, and 360 cm. The Fault interactions along the EAFZ were obvious during the consecutive earthquakes. The ruptures mainly occurred in the left-lateral components of the fault segments, with the maximum slips of 7.25 m and 9.43 m for the first event along the EAFZ and the second event on the Çardak Fault, respectively.

Revisiting the Variation of the Ionospheric Irregularities in the Low Latitude Region of China Based on Small Regional Geodetic GNSS Station Network

AbstractThe rate of total electron content index (ROTI) and loss of lock (LoL) from global navigation satellite system (GNSS) observations are an important data source for ionospheric scintillation study. However, there are certain limitations of these data in ionospheric scintillation study, and clearing these limitations is important to understand the results from these data. In this paper, the variation of the ionospheric scintillation is revisited based on LoL and ROTI from a spatial dense GNSS network in the low latitude region of China. Via a priori knowledge of scintillation morphology, the method to determine the baseline ROTI for discriminating quiet and disturbed conditions is provided, and the abnormal data of ROTI and LoL unrelated to ionosphere are found and eliminated. Results show that the data examination is necessary. Using the qualified ROTI and LoL, the morphological variations of ionospheric irregularities are revisited. The temporal variation of ROTI and LoL are generally consistent, but there are some discrepancies. The maximum LoL occurrence is at about 21:00 LT, it is between 20:00 and 23:00 LT for ROTI. For spatial distribution, both parameters reach maximum occurrence in the southwest direction, but LoL is more concentrated. The azimuth range of maximum occurrence is 180–190° for LoL and 210–220° for ROTI. Statistically, the correlation between LoL occurrence and ROTI value from observation of single GNSS station is relatively good when ROTI <5 TECU/min, but it becomes worse when ROTI >5 TECU/min. However, their correlation can be greatly improved when data are averaged among all the stations.

Road extraction from low-cost GNSS-device dense trajectories

ABSTRACT This paper proposes a method for road centerline extraction from dense Global Navigation Satellite System (GNSS) trajectories, collected by using low-cost GNSS-devices, i.e. smartphones. The proposed method basically consists in generating a frequency image by tracking the GNSS trajectories and then by applying the Steger line detector to the generated image to extract the road centerlines. The main motivation of using the Steger algorithm is its capability to detect lines with sub-pixel accuracy. To evaluate the obtained results, reference road centerlines are manually extracted from a georeferenced orthomosaic. The experiments performed demonstrate the high potential of applying the Steger line detector to frequency images, generated by using dense GPS (Global Positioning System) trajectories. The completeness and correctness values for the accomplished experiments were 98% and 99%, respectively. Additionally, the RMSE (Root Mean Square Error) ranged from 0.63 m to 2.39 m, or approximately 1/16 to 1/4 of the expected accuracy (about 10 m) of a point determined by the Single-Point Positioning (SPP) method, which is the GNSS positioning method usually employed by smartphones.

Multiwave Structure of Traveling Ionospheric Disturbances Excited by the Tonga Volcanic Eruptions Observed by a Dense GNSS Network in China

AbstractWe used dense global navigation satellite system data from China to track the propagation of traveling ionospheric disturbances (TIDs) triggered by the 2022 January 15 Tonga volcanic eruption. We identified two TIDs originating from the eruption. One, which has been reported widely by a number of recent investigations, had a velocity of ∼361 m/s. However, another long‐distance propagating TID with a velocity of ∼264 m/s has not been widely discussed. The velocities of these TIDs coincide with previous simulation results of gravity‐wave L0 and L1 modes. We propose that these TIDs were caused by the L0 and L1 ducted modes of gravity waves excited by the volcanic eruption. However, the L1 mode is usually negligible due to its weak amplitude in comparison with that of the L0 mode. The enormous energy release resulted in a stronger amplitude of the L1 mode, which induced a detectable TID in the ionosphere.

Cluster Analysis of Dense GNSS Velocity Field Reveals Characteristics Associated With Regional Tectonics in New Zealand

AbstractNew Zealand is located in the converging zone between the Australian and Pacific Plates. The Pacific Plate is descending beneath the Australian Plate in the north of New Zealand, while the Australian Plate is descending beneath the Pacific Plate in the south. The contrastive tectonics is connected by the Alpine Fault, the transform fault with NE‐SW strike along most of the length of New Zealand. Recently, previous studies showed that a cluster analysis of dense Global Navigation Satellite System (GNSS) velocity field objectively yields major velocity discontinuities, which may reflect an important tectonic boundary. We applied two different types of clustering algorithms: hierarchical agglomerative clustering (HAC) and Euler pole clustering (EPC), to the dense GNSS horizontal velocity data in New Zealand. Both methods agree that the primary kinematic boundary lies from the northern Wairarapa region to the southwestern end of the Alpine Fault. The HAC introduces a cluster boundary along the North Island Dextral Fault Belt whereas EPC emphasizes the importance of the Taupo Volcanic Zone (TVZ). Since EPC fits the data reasonably well, it indicates the TVZ is an important boundary. This result displays the limitation of the flat earth assumption. The Chatham Islands is identified as an independent cluster by the HAC, but not by the EPC. Further cluster splits illustrate that the relative motion between the two plates is accommodated by distributed deformation in the Marlborough Fault Zone in the South Island whereas the North Island is characterized by local block rotations.

Rapid Source Determination of the 2021 Mw 7.4 Maduo Earthquake by a Dense High-Rate GNSS Network in the Tibetan Plateau

Abstract High-rate Global Navigation Satellite Systems (GNSS) observations from an integrated network in the Tibetan plateau were used for rapid source determination of the 2021 Mw 7.4 Maduo earthquake. Constrained by high-rate displacements from nine stations near the source, the magnitude was essentially determined as early as 20 s after the origin time (OT), whereas the 39-s-long rupture propagated only approximately halfway, and together with the epicenter location and rupture length, they were completely determined 29 s after the OT or 10 s before the rupture completion. In addition, to test the capability and feasibility of current GNSS networks for rapid source determination in the Tibetan plateau, two experiments were further implemented to invert source parameters by real observations at farther sites and synthetic observations induced by an Mw 7.0 modeled earthquake. The results demonstrate that current GNSS networks are capable and feasible to provide rapid source determination and earthquake early warning in approximately 30 s after the OT for large earthquakes (Mw 7.0+) based on current station density, spatial coverage and communication conditions in the Tibetan plateau.

Updated absolute gravity rate of change associated with glacial isostatic adjustment in Southeast Alaska and its utilization for rheological parameter estimation

In Southeast Alaska (SE-AK), rapid ground uplift of up to 3 cm/yr has been observed associated with post-Little Ice Age glacial isostatic adjustment (GIA). Geodetic techniques such as global navigation satellite system (GNSS) and absolute gravimetry have been applied to monitor GIA since the last 1990s. Rheological parameters for SE-AK were determined from dense GNSS array data in earlier studies. However, the absolute gravity rate of change observed in SE-AK was inconsistent with the ground uplift rate, mainly because few gravity measurements from 2006 to 2008 resulted in imprecise gravity variation rates. Therefore, we collected absolute gravity data at six gravity points in SE-AK every June in 2012, 2013, and 2015, and updated the gravity variation rate by reprocessing the absolute gravity data collected from 2006 to 2015. We found that the updated gravity variation rate at the six gravity points ranged from −2.05 to −4.40 upmuGal/yr, and its standard deviation was smaller than that reported in the earlier study by up to 88 %. We also estimated the rheological parameters under the assumption of the incompressible Earth to explain the updated gravity variation rate, and their optimal values were determined to be 55 km and 1.2 times 10^{19} Pa s for lithospheric thickness and upper mantle viscosity, respectively. These optimal values are consistent with those independently obtained from GNSS observations, and this fact indicates that absolute gravimetry can be one of the most effective methods in determining sub-surface structural parameters associated with GIA accurately. Moreover, we utilized the gravity variation rates for estimating the ratio of gravity variation to vertical ground deformation at the six gravity points in SE-AK. The viscous ratio values were obtained as −0.168 and −0.171 upmuGal/mm from the observed data and the calculated result, respectively. These ratios are greater (in absolute) than those for other GIA regions (−0.15 to −0.16 upmuGal/mm in Antarctica and Fennoscandia), because glaciers in SE-AK have melted more recently than in other regions.Graphical

Fusion of Spatially Heterogeneous GNSS and InSAR Deformation Data Using a Multiresolution Segmentation Algorithm and Its Application in the Inversion of Slip Distribution

The fusion of global navigation satellite system (GNSS) and interferometric synthetic aperture radar (InSAR) deformation data can leverage the advantages of GNSS high temporal resolution and InSAR high spatial resolution, and obtain more abundant deformation data for constraints on geophysical structural and mechanical parameters. Existing studies seldom consider the spatial heterogeneity of largescale deformation data, which easily leads to obvious spatial aggregation of errors in the results of fusion. Here, we propose a novel multiresolution segmentation fusion (MRSF) method that uses a multiresolution segmentation algorithm to automatically classify the spatial heterogeneity of InSAR deformation data with similar deformation characteristics. We applied the MRSF method to the fusion of GNSS and InSAR deformation data covering the central valley aquifer system (CVAS) in southern California to verify its precision and robustness. Results show that the MRSF method can accurately reflect spatiotemporal evolution characteristics of displacement data and reliably estimate deformation for the times and locations of missing data. We then tested this method for geophysical parameter estimation by constructing three different sets of data, including dense GNSS sites, sparse GNSS sites, and sparse GNSS sites fused with InSAR data using MRSF, to invert the slip distribution of the Cascadia subduction zone. Results show that the inverted slip of the fused InSAR and GNSS data is comparable to that of the dense GNSS sites. Therefore, the MRSF method can obtain deformation results with high precision and high spatiotemporal resolution and effectively compensate for the lack of data caused by sparse GNSS sites during the geophysical inversion process.

On the Potential of Rapid Moment Magnitude Estimation for Strong Earthquakes in Sichuan–Yunnan Region, China, Using Real-Time CMONOC GNSS Observations

Abstract High-rate (≥1 Hz) real-time Global Navigation Satellite Systems (GNSS) data have been now incorporated into earthquake early warning systems where dense GNSS network exists, like the United States west coast, Japan, and so on. In this manuscript, we exploit the potential of rapid moment magnitude (Mw) estimation in Sichuan–Yunnan region using GNSS records from Crustal Movement Observation Network of China (CMONOC) for seven magnitude (M) ≥ 6 events from 2011 to 2020. We find that currently it is still quite challenging to get reliable Mw using the peak ground displacement or peak ground velocity approach, which is caused by the sparse stations and probably the inaccurate scaling law itself. By contrast, the Mw estimations from focal mechanism grid search are all consistent with Global Centroid Moment Tensor catalog, even by one station. Considering the CMONOC network is going to be densified in the coming decade, it is expected that high-rate GNSS waveforms will play a more important role in fast postevent response.

High-precision GNSS PWV retrieval using dense GNSS sites and in-situ meteorological observations for the evaluation of MERRA-2 and ERA5 reanalysis products over China

Precipitable water vapor (PWV) product with high accuracy and high spatiotemporal resolution is important for climate change research. Hourly PWV products with a high spatiotemporal resolution can be provided by global navigation satellite system (GNSS) and global reanalysis data. The National Aeronautics and Space Administration (NASA) and the European Center for Medium-Range Weather Forecasts (ECMWF) have recently released their global reanalysis of the Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2) and the fifth-generation ECMWF Reanalysis (ERA5), respectively. A comprehensive evaluation of these two PWV reanalysis products over China has yet to be fully investigated. In this study, the PWV products of MERRA-2 and ERA5 reanalysis over China have been systematically evaluated based on 339 GNSS sites and in-situ ground meteorological observations from year 2016 to 2018. To mitigate the uncertainties of this evaluation, three enhanced vertical correction models for temperature, pressure and PWV that we developed previously are adopted. Four weighted mean temperature (Tm) models are established for the four regions of China (North China, South China, Tibet Plateau, and Northwest China), and they exhibit better accuracy than the other classical Tm models in China. The evaluation results of the MERRA-2- and ERA5-derived PWV show that the correlation coefficient, biases, and root mean square error (RMSE) in China are 0.98/0.99, −0.51/0.38 mm, and 2.50/1.99 mm, respectively. They both exhibit relatively large difference when applied to South China. For the same grid spatial resolution, ERA5 PWV product exhibits similar performance with MERRA-2. The monthly-mean-reanalysis-derived PWV exhibits obvious seasonality and has the largest difference in summer. ERA5 PWV can provide more details of the spatial pattern compared to MERRA-2 PWV. The diurnal anomaly variation for ERA5 is more consistent with GNSS data, and some differences at certain moments during autumn and winter have been observed. These results show that ERA5 PWV and MERRA-2 PWV exhibits excellent applicability in China and can help us to efficiently use PWV products for climate research in China.

Joint Inversion of GNSS and GRACE for Terrestrial Water Storage Change in California.

Global Navigation Satellite System (GNSS) vertical displacements measuring the elastic response of Earth's crust to changes in hydrologic mass have been used to produce terrestrial water storage change (∆TWS) estimates for studying both annual ∆TWS as well as multi‐year trends. However, these estimates require a high observation station density and minimal contamination by nonhydrologic deformation sources. The Gravity Recovery and Climate Experiment (GRACE) is another satellite‐based measurement system that can be used to measure regional TWS fluctuations. The satellites provide highly accurate ∆TWS estimates with global coverage but have a low spatial resolution of ∼400 km. Here, we put forward the mathematical framework for a joint inversion of GNSS vertical displacement time series with GRACE ∆TWS to produce more accurate spatiotemporal maps of ∆TWS, accounting for the observation errors, data gaps, and nonhydrologic signals. We aim to utilize the regional sensitivity to ∆TWS provided by GRACE mascon solutions with higher spatial resolution provided by GNSS observations. Our approach utilizes a continuous wavelet transform to decompose signals into their building blocks and separately invert for long‐term and short‐term mass variations. This allows us to preserve trends, annual, interannual, and multi‐year changes in TWS that were previously challenging to capture by satellite‐based measurement systems or hydrological models, alone. We focus our study in California, USA, which has a dense GNSS network and where recurrent, intense droughts put pressure on freshwater supplies. We highlight the advantages of our joint inversion results for a tectonically active study region by comparing them against inversion results that use only GNSS vertical deformation as well as with maps of ∆TWS from hydrological models and other GRACE solutions. We find that our joint inversion framework results in a solution that is regionally consistent with the GRACE ∆TWS solutions at different temporal scales but has an increased spatial resolution that allows us to differentiate between regions of high and low mass change better than using GRACE alone.

Hydrological drought characterization based on GNSS imaging of vertical crustal deformation across the contiguous United States

Continuous Global Navigation Satellite System (GNSS) measurements allow us to track subtle elastic crustal deformation in the response to hydrological mass variations and provide an additional tool to independently characterize hydrological extremes (e.g., droughts and floods). In this study, we develop a time-varying GNSS imaging strategy that depends on the principal component analysis of GNSS-sensed vertical crustal displacement (VCD) in 2006–2020 and the monthly images of hydrology-induced deformation are generated for drought characterization across the contiguous United States. The first 12 principal components are selected in our time-varying imaging system, which account for 85% of the data variance. Considering that surface water loads are inversely correlated with the induced elastic vertical motions, we reverse the signs of the GNSS-imaged time series in all grids in subsequent studies (referred to as negative VCD (NVCD)). The GNSS-NVCD data generally correlate well with the water estimates from the Gravity Recovery and Climate Experiment (GRACE) and North American Land Data Assimilation System (NLDAS). Using the GNSS-imaged gridded NVCD products, we produce a GNSS-based drought severity index (GNSS-DSI) based on the climatological methodology, which is implemented by standardizing the GNSS NVCD anomalies that deviate from climatological normal. In most regions, strong linear correlations are accessible for GNSS-DSI relative to GRACE-DSI and the self-calibrating Palmer Drought Severity Index (scPDSI). The new drought monitoring tool, which is based solely on GNSS-measured vertical positions, is used for hydrological drought characterization (onset, end, duration, magnitude, intensity, and recovery); it succeeds in identifying well-documented historical droughts from the US drought monitor (USDM). Our study presents a new drought characterization framework using solely GNSS-measured hydrological loading displacements from a dense GNSS network, which has great potential to strengthen operational drought monitoring and assessment.

Study of structures of the sporadic E layer by using dense GNSS network observations

<h3>Abstract</h3> The sporadic E (Es) layer has been known to introduce long-range propagation of aeronautical very high frequency (VHF) navigation beyond the radio horizon and cause potential interference on the navigation system. This study utilizes a rate of total electron content (TEC) index (ROTI) map with dense Global Navigation Satellite System (GNSS) observations for effective Es layer detection. The daytime Es layer shows a well-defined frontal structure when ROTI values are mapped at the typical Es layer height (100 km). A methodology of detecting and characterizing the Es layer frontal structure without manual operation is developed by utilizing the Hough transform. The front direction and drift velocity are successfully derived. Sub-structures in the Es layer front are revealed by analysis using the characteristics of the frontal structure and TEC variation. The developed method is suitable for an automated real-time Es-layer monitoring system in a wide area.

Ionospheric Responses at Low Latitudes to the Annular Solar Eclipse on 21 June 2020

AbstractIn this study, we utilized both ground‐based and space‐borne observations including total electron content (TEC) from Beidou geostationary satellites, two‐dimensional TEC maps from the worldwide dense Global Navigation Satellite System receivers, ionosondes, and in situ electron density (Ne) and electron temperature (Te) from both Swarm and China Seismo‐Electromagnetic Satellite satellites, to investigate the low‐latitude ionospheric responses to the annular solar eclipse on 21 June 2020. The decrease in TEC during the eclipse at low latitudes showed a local time dependence with the largest depletions in the noon and afternoon sectors. It was also found that the TEC depletions at different latitudes in the equatorial ionization anomaly (EIA) region over the East Asian sector cannot solely be explained by the solar flux changes associated with the obscuration rate. The differences in TEC reduction between stations can be more than a factor of 2 at latitudes with the same obscuration rate of over 90%. Compared with TEC variations in the Northern Hemisphere, the TEC also underwent a considerable decrease in the EIA region in the conjugate hemisphere without eclipse shadow. Meanwhile, the hmF2 near the magnetic equator increased around the onset of the eclipse, indicating an enhancement of the eastward equatorial electric field. Furthermore, the TEC decrease during the eclipse in the EIA region in both hemispheres lasted for a long period of more than 7 hr after the eclipse, with a TEC depletion of 2–6 TEC units. The Ne from Swarm and China Seismo‐Electromagnetic Satellite satellites showed a complicated variation after the eclipse, whereas no visible change was observed in Te. The enhanced equatorial electric field, neutral wind changes, and the associated plasma transport act together to generate the observed ionospheric effects at low latitudes during the eclipse. Our results also suggest that the eclipse‐induced perturbations of dynamic processes can continue to impact the ionosphere after the eclipse.

Toward Near‐Field Tsunami Forecasting Along the Cascadia Subduction Zone Using Rapid GNSS Source Models

AbstractOver the past 15 years and through multiple large and devastating earthquakes, tsunami warning systems have grown considerably in their efficacy in providing timely and accurate forecasts to affected communities. However, one part of tsunami warning that still needs improvement is forecasts catered to local, near‐field communities in the time after an earthquake rupture but before coastal inundation. In this study, we test a rapid, Global Navigation Satellite Systems (GNSS)‐driven earthquake characterization model using a large data set of synthetic megathrust ruptures for its near‐field tsunami forecasting potential. We also provide a framework for tsunami forecasting that focuses on the likelihood of exceedance of user‐defined coastal amplitudes that may be of use in the first 15 min following an earthquake. Specifically, we can estimate the earthquake magnitude, without saturation, for 82% of tested ruptures. We can also identify test ruptures as dominantly thrust events, without analyst guidance for 92% of tested thrust ruptures. Finally, modeling the tsunami component of our rapidly estimated fault rupture leads to greater than 80% accuracy in identifying tsunami impact at key coastal amplitude thresholds. This is promising for near‐field warning when the time prior to inundation is limited to tens of minutes. We focus this study on large megathrust ruptures along the quiescent Cascadia subduction zone where there is already a dense GNSS network.