Evaluation of Earthquake Magnitude Estimation and Event Detection Thresholds for Real-Time GNSS Networks: Examples from Recent Events Captured by the Network of the Americas

2020 ◽  
Vol 91 (3) ◽  
pp. 1628-1645 ◽  
Author(s):  
Kathleen M. Hodgkinson ◽  
David J. Mencin ◽  
Karl Feaux ◽  
Charles Sievers ◽  
Glen S. Mattioli

Abstract Several studies have shown that real-time (RT) Global Navigation Satellite Systems (GNSS) measurements can provide an estimate of an earthquake’s moment magnitude (Mw) using scaling laws that relate peak ground displacements (PGDs) and hypocentral distance with Mw. In this study, we use data from GNSS stations operated by UNAVCO as part of the National Science Foundation Geodetic Facility for the Advancement of Geoscience (GAGE) that comprises the Network of the Americas (NOTA) to show that precise point positioning (PPP) solutions distributed in RT during five recent earthquakes could be used to calculate Mw rapidly and reliably. We analyze solutions distributed by UNAVCO during the 8 September 2017 Mw 8.2 Tehuantepec, Mexico, 10 January 2018 Mw 7.5 Great Swan Island, Honduras, 23 January 2018 Mw 7.9 Gulf of Alaska earthquake, and the 4 July 2019 Mw 6.4 and 6 July 2019 Mw 7.1 Ridgecrest, California, earthquakes. We find that RT-GNSS Mw estimates consistent with Advanced National Seismic System Comprehensive Earthquake Catalog values are available tens of seconds to a few minutes after the event origin time. The speed with which an estimate is available is dependent on the proximity to the epicenter of the closest NOTA stations. The results demonstrate that RT-GNSS networks could be used to mitigate the problem of magnitude saturation observed in seismic-based earthquake early warning (EEW) systems. RT-GNSS effectively expands the spectrum of events for which a seismic EEW system can provide accurate warnings of impending ground shaking and provides an independent verification of seismically derived magnitudes. We also analyze the RT-PPP solutions from more than 800 RT-GNSS stations to determine the ambient-noise levels of each NOTA station and combine that with PGD magnitude scaling laws to construct regional network sensitivity maps for the NOTA network. Such maps may be used to determine “blind spots” or regions of lower sensitivity in RT-GNSS networks under consideration for EEW.


2018 ◽  
Vol 52 (1) ◽  
pp. 98 ◽  
Author(s):  
Athanassios Ganas ◽  
Nikoletta Andritsou ◽  
Chrysanthi Kosma ◽  
Panagiotis Argyrakis ◽  
Varvara Tsironi ◽  
...  

We describe and make available a dataset of 64 data points of Global Positioning System (GPS) displacements for significant, shallow earthquakes in Greece during the period 1997-2017. The displacement data can be used by earthquake geologists, engineers and seismologists in an effort to better understand the faulting process, the rupture mechanics, the pattern of ground-motions, and in engineering applications. We include recordings from GNSS networks at near-source to regional distances (2–132 km) for 11 earthquakes between global CMT moment magnitudes (Mw) 5.5 and 6.9. We also model the magnitude scaling properties of peak ground horizontal displacements (PGD and PGD-S) for these events using L1-norm minimisation regression. Our data indicate an almost linear attenuation of seismic strain with distance for this range of seismic magnitudes. We developed a set of relationships based on PGD (in cm) and distance to hypocentre R (in km), which may be used for the rapid estimation of the earthquake magnitude in near real-time.MwPGD = [LOG(PGD) + 8.2849]/(1.6810 – 0.2453LOGR)MwPGD-S = [LOG(PGD-S) + 8.0839]/(1.6793 – 0.2447LOGR)



2021 ◽  
Vol 13 (11) ◽  
pp. 2179
Author(s):  
Pedro Mateus ◽  
Virgílio B. Mendes ◽  
Sandra M. Plecha

The neutral atmospheric delay is one of the major error sources in Space Geodesy techniques such as Global Navigation Satellite Systems (GNSS), and its modeling for high accuracy applications can be challenging. Improving the modeling of the atmospheric delays (hydrostatic and non-hydrostatic) also leads to a more accurate and precise precipitable water vapor estimation (PWV), mostly in real-time applications, where models play an important role, since numerical weather prediction models cannot be used for real-time processing or forecasting. This study developed an improved version of the Hourly Global Pressure and Temperature (HGPT) model, the HGPT2. It is based on 20 years of ERA5 reanalysis data at full spatial (0.25° × 0.25°) and temporal resolution (1-h). Apart from surface air temperature, surface pressure, zenith hydrostatic delay, and weighted mean temperature, the updated model also provides information regarding the relative humidity, zenith non-hydrostatic delay, and precipitable water vapor. The HGPT2 is based on the time-segmentation concept and uses the annual, semi-annual, and quarterly periodicities to calculate the relative humidity anywhere on the Earth’s surface. Data from 282 moisture sensors located close to GNSS stations during 1 year (2020) were used to assess the model coefficients. The HGPT2 meteorological parameters were used to process 35 GNSS sites belonging to the International GNSS Service (IGS) using the GAMIT/GLOBK software package. Results show a decreased root-mean-square error (RMSE) and bias values relative to the most used zenith delay models, with a significant impact on the height component. The HGPT2 was developed to be applied in the most diverse areas that can significantly benefit from an ERA5 full-resolution model.



Author(s):  
Dino Bindi ◽  
Hoby N. T. Razafindrakoto ◽  
Matteo Picozzi ◽  
Adrien Oth

ABSTRACT We investigate the impact of considering a depth-dependent attenuation model on source parameters assessed through a spectral decomposition. In particular, we evaluate the effect of considering the hypocentral depth as an additional variable for the attenuation model, using as the target the tendency of the average stress drop to increase with depth, as observed in recent studies. We analyze the Fourier spectra of S-wave windows for about 1900 earthquakes with a magnitude above 2.5 recorded in the Ridgecrest region, southern California. Two different parameterizations of the attenuation term are implemented in the spectral decomposition, either as a function of the hypocentral distance alone or as a function of both epicentral distance and depth. The comparison of the spectral attenuation curves shows that, although the hypocentral model describes, on average, the range of values spanned by the attenuation curve for different depths, systematic differences with distance, depth, and frequency are observed. These differences are transferred to the source spectra and, in turn, to the source parameters extracted from the best-fitting ω−2 models. In particular, stress drops for events deeper than 7 km are, on average, almost double even when depth is introduced explicitly in the attenuation model. The increase of stress drop with depth is confirmed also after accounting for the increase of the shear velocity with depth, which absorbs about 30%–40% of the total increase. Moreover, a qualitative comparison with a model for the gradient of the effective normal stress confirms the reliability of the observed trend. Finally, the coherent spatial patterns shown by a simplified 2D tomographic representation of the spectral residuals highlights the impact on ground-shaking variability of the lateral variability of the crustal attenuation properties in the region.



2010 ◽  
Vol 10 (7) ◽  
pp. 1617-1627 ◽  
Author(s):  
A. Y. Babeyko ◽  
A. Hoechner ◽  
S. V. Sobolev

Abstract. We present the GITEWS approach to source modeling for the tsunami early warning in Indonesia. Near-field tsunami implies special requirements to both warning time and details of source characterization. To meet these requirements, we employ geophysical and geological information to predefine a maximum number of rupture parameters. We discretize the tsunamigenic Sunda plate interface into an ordered grid of patches (150×25) and employ the concept of Green's functions for forward and inverse rupture modeling. Rupture Generator, a forward modeling tool, additionally employs different scaling laws and slip shape functions to construct physically reasonable source models using basic seismic information only (magnitude and epicenter location). GITEWS runs a library of semi- and fully-synthetic scenarios to be extensively employed by system testing as well as by warning center personnel teaching and training. Near real-time GPS observations are a very valuable complement to the local tsunami warning system. Their inversion provides quick (within a few minutes on an event) estimation of the earthquake magnitude, rupture position and, in case of sufficient station coverage, details of slip distribution.



Author(s):  
Maryna Sapachova

ObjectiveThe performance of comparative analysis of sensitivity and resultsof detection of avian influenza virus by real time polymerase chainreaction (PCR-RT) and loop-mediated isothermal amplification of thenucleic acids (LAMP) was the main goal of the study.IntroductionAs part of this surveillance study for Avian Influenza both activeand passive surveillance samples were tested using PCR and alsoutilized to validate the LAMP method. Active surveillance samplesinclude pathological material and tracheal and cloacal swabs fromill poultry, which were subsequently assessed for avian influenzaduring diagnosis, and birds collected by hunters. Passive surveillanceincluded environmental samples such as sand and bird faeces.Active surveillance samples were taken mostly from poultry farmsacross Ukraine, where infected birds are required to be diagnosedby State Scientific Research Institute of Laboratory Diagnosticsand Veterinary Sanitary Expertise (SSRILDVSE) by Ukraine Law.Passive surveillance samples were taken primarily during the annualbird migration season. Development of simple, sensitive, and cheapmethods for diagnostics of avian influenza is a very important taskfor practical veterinary medicine. LAMP is one of such methods.The technique is based on isothermal amplification of nucleic acids.It does not require special conditions and equipment (PCR cyclers),therefore it is cheaper in comparison with PCR. Accurate diagnosisis necessary for determining the risk associated with avian influenzain Ukraine and along the Dnipro River during the migratory season.MethodsFor the research, we used PCR-RT commercial kit Bird-Flu-PCR(Ukrzoovetprompostach, Ukraine), LAMP (the protocol has beenoptimized and patented by SSRILDVSE), QIAamp®Viral RNA MiniKit. For the study, we used pathological and biological materials frombirds, which were sent to the SSRILDVSE from all regions of Ukraineaccording to the 2013–2014 State monitoring plan.Set up of the real time PCR reactions and parameters ofamplifications are indicated in the instruction to the kit.The following protocol was used to set up the RT- LAMP: 2.5μL10 X Thermopol buffer, 1 mmol/L betaine, 5 mmol/L MgSO4,1.4 mmol/L - BNTP, 12.5μmol/L SYBR GREEN, 0.5 mmol/LMnCL2, up to 25μL Nuclease-free water, 8 U Bsm DNA polymerase,0.1μM/1 of F3, 0.1μM/1 of B3, 0.8μM/1 of FIP, 0.8μM/1 of BIP,0.4μM/1 of LF, 0.4 of LB, 2μL cDNA.During our work, we used the following optimal temperature andtime for the amplification – 59°C and 60 minutes.The sensitivity of diagnostic kit Bird-Flu-PCR and RT- LAMP wasdetermined by testing cDNA of the reference strain of AIV H5N1,which was provided to us by NSC Institute for Experimental andClinical Veterinary Medicine (Kharkiv, Ukraine). For the standard,we employed concentration in the range of 10.0-0.01 ng/sample.ResultsTable 1.This table shows the reproducibility results obtained by bothmethods. However, taken into account absence of highly pathogenicavian influenza virus circulating in Ukraine during the studied period,it was not possible to confirm these results with protocols of positivesamples.Table 2.It has been established that the sensitivity of PCR-RT kit Bird-Flu-PCR is 0.01 ng/sample for gene M and 0.1 ng/sample for subtypeH5N1.Fig. 1. Visual detection of LAMP products with differentconcentrations of cDNA of avian influenza virus (ng per sample):1 – 10; 2 – 5; 3 – 1.0; 4 – 0.1; 5–7 – 0.01; 8–9 – 0.1; 10 – negative.We have examined the LAMP results using electrophoresis forthe confirmation of visual detection and correct interpretation of theresults (Fig. 2).Fig.2. Electrophoresis results for LAMP products. M –molecular weight marker; 1 – 10.0; 2 – 5.0; 3 – 1.0; 4 – 0.1; 5–7– 0.01; 8 - negative control.It has been established that the sensitivity of LAMP is0.1 ng/sample. Slightly lower sensitivity of LAMP in comparisonto PCR-RT can be explained by visual detection of the products ofthe LAMP reaction.Conclusions1. Sensitivity of both methods is high.2. LAMP is a perspective screening method for the diagnosis ofviral infectious diseases supported by confirmation of positive resultsby PCR-RT.



2018 ◽  
Vol 52 (3) ◽  
pp. 100-108 ◽  
Author(s):  
Takeshi Nakamura ◽  
Narumi Takahashi ◽  
Kensuke Suzuki

AbstractThe deployment of real-time permanent ocean-bottom seismic and tsunami observatories is significant for disaster mitigation and prevention during the occurrence of large subduction earthquakes near trough areas. On April 1, 2016, a moderate-sized suboceanic earthquake occurred beneath Dense Oceanfloor Network System for Earthquakes and Tsunamis (DONET) stations that were recently deployed in deep ocean-bottom areas near the Nankai Trough in southwest Japan. P-waves arrived at the ocean-bottom station within 4 s of the origin time, which was 6 and 13 s earlier than the arrival of P- and S-waves at a land station in the coastal area, respectively; this implies earlier detection of strong motion than at land stations. However, the waveforms are amplified by sediment layers and even contaminated with acceleration offsets at some stations, which would lead to overestimations during source investigations. Such amplification and offset did not occur at a borehole station connected to DONET. The amplifications caused by the sediment layers and the offset were found to have a considerable spatial variation, not only between the DONET stations and land and borehole stations but also among the DONET stations, implying that the amplitude evaluation could be unstable. Therefore, procedures for correcting or suppressing the amplification and offset problem are required for conducting waveform analyses, such as magnitude estimations and source modeling, during large subduction earthquakes.



2019 ◽  
Vol 37 (3) ◽  
pp. 429-446 ◽  
Author(s):  
Michal Kačmařík ◽  
Jan Douša ◽  
Florian Zus ◽  
Pavel Václavovic ◽  
Kyriakos Balidakis ◽  
...  

Abstract. An analysis of processing settings impacts on estimated tropospheric gradients is presented. The study is based on the benchmark data set collected within the COST GNSS4SWEC action with observations from 430 Global Navigation Satellite Systems (GNSS) reference stations in central Europe for May and June 2013. Tropospheric gradients were estimated in eight different variants of GNSS data processing using precise point positioning (PPP) with the G-Nut/Tefnut software. The impacts of the gradient mapping function, elevation cut-off angle, GNSS constellation, observation elevation-dependent weighting and real-time versus post-processing mode were assessed by comparing the variants by each to other and by evaluating them with respect to tropospheric gradients derived from two numerical weather models (NWMs). Tropospheric gradients estimated in post-processing GNSS solutions using final products were in good agreement with NWM outputs. The quality of high-resolution gradients estimated in (near-)real-time PPP analysis still remains a challenging task due to the quality of the real-time orbit and clock corrections. Comparisons of GNSS and NWM gradients suggest the 3∘ elevation angle cut-off and GPS+GLONASS constellation for obtaining optimal gradient estimates provided precise models for antenna-phase centre offsets and variations, and tropospheric mapping functions are applied for low-elevation observations. Finally, systematic errors can affect the gradient components solely due to the use of different gradient mapping functions, and still depending on observation elevation-dependent weighting. A latitudinal tilting of the troposphere in a global scale causes a systematic difference of up to 0.3 mm in the north-gradient component, while large local gradients, usually pointing in a direction of increasing humidity, can cause differences of up to 1.0 mm (or even more in extreme cases) in any component depending on the actual direction of the gradient. Although the Bar-Sever gradient mapping function provided slightly better results in some aspects, it is not possible to give any strong recommendation on the gradient mapping function selection.



Sensors ◽  
2019 ◽  
Vol 19 (15) ◽  
pp. 3376 ◽  
Author(s):  
Yuan Du ◽  
Guanwen Huang ◽  
Qin Zhang ◽  
Yang Gao ◽  
Yuting Gao

Real-time kinematic (RTK) positioning is a satellite navigation technique that is widely used to enhance the precision of position data obtained from global navigation satellite systems (GNSS). This technique can reduce or eliminate significant correlation errors via the enhancement of the base station observation data. However, observations received by the base station are often interrupted, delayed, and/or discontinuous, and in the absence of base station observation data the corresponding positioning accuracy of a rover declines rapidly. With the strategies proposed till date, the positioning accuracy can only be maintained at the centimeter-level for a short span of time, no more than three min. To address this, a novel asynchronous RTK method (that addresses asynchronous errors) that can bridge significant gaps in the observations at the base station is proposed. First, satellite clock and orbital errors are eliminated using the products of the final precise ephemeris during post-processing or the ultra-rapid precise ephemeris during real-time processing. Then the tropospheric error is corrected using the Saastamoinen model and the asynchronous ionospheric delay is corrected using the carrier phase measurements from the rover receiver. Finally, a straightforward first-degree polynomial function is used to predict the residual asynchronous error. Experimental results demonstrate that the proposed approach can achieve centimeter-level accuracy for as long as 15 min during interruptions in both real-time and post-processing scenarios, and that the accuracy of the real-time scheme can be maintained for 15 min even when a large systematic error is projected in the U direction.



Sensors ◽  
2020 ◽  
Vol 20 (13) ◽  
pp. 3676 ◽  
Author(s):  
Tao Ni ◽  
Wenhang Li ◽  
Dingxuan Zhao ◽  
Zhifei Kong

Autonomous vehicles can achieve accurate localization and real-time road information perception using sensors such as global navigation satellite systems (GNSSs), light detection and ranging (LiDAR), and inertial measurement units (IMUs). With road information, vehicles can navigate autonomously to a given position without traffic accidents. However, most of the research on autonomous vehicles has paid little attention to road profile information, which is a significant reference for vehicles driving on uneven terrain. Most vehicles experience violent vibrations when driving on uneven terrain, which reduce the accuracy and stability of data obtained by LiDAR and IMUs. Vehicles with an active suspension system, on the other hand, can maintain stability on uneven roads, which further guarantees sensor accuracy. In this paper, we propose a novel method for road profile estimation using LiDAR and vehicles with an active suspension system. In the former, 3D laser scanners, IMU, and GPS were used to obtain accurate pose information and real-time cloud data points, which were added to an elevation map. In the latter, the elevation map was further processed by a Kalman filter algorithm to fuse multiple cloud data points at the same cell of the map. The model predictive control (MPC) method is proposed to control the active suspension system to maintain vehicle stability, thus further reducing drifts of LiDAR and IMU data. The proposed method was carried out in outdoor environments, and the experiment results demonstrated its accuracy and effectiveness.



Sensors ◽  
2021 ◽  
Vol 21 (6) ◽  
pp. 1956
Author(s):  
Natalia Wielgocka ◽  
Tomasz Hadas ◽  
Adrian Kaczmarek ◽  
Grzegorz Marut

Global Navigation Satellite Systems (GNSS) have revolutionized land surveying, by determining position coordinates with centimeter-level accuracy in real-time or up to sub-millimeter accuracy in post-processing solutions. Although low-cost single-frequency receivers do not meet the accuracy requirements of many surveying applications, multi-frequency hardware is expected to overcome the major issues. Therefore, this paper is aimed at investigating the performance of a u-blox ZED-F9P receiver, connected to a u-blox ANN-MB-00-00 antenna, during multiple field experiments. Satisfactory signal acquisition was noticed but it resulted as >7 dB Hz weaker than with a geodetic-grade receiver, especially for low-elevation mask signals. In the static mode, the ambiguity fixing rate reaches 80%, and a horizontal accuracy of few centimeters was achieved during an hour-long session. Similar accuracy was achieved with the Precise Point Positioning (PPP) if a session is extended to at least 2.5 h. Real-Time Kinematic (RTK) and Network RTK measurements achieved a horizontal accuracy better than 5 cm and a sub-decimeter vertical accuracy. If a base station constituted by a low-cost receiver is used, the horizontal accuracy degrades by a factor of two and such a setup may lead to an inaccurate height determination under dynamic surveying conditions, e.g., rotating antenna of the mobile receiver.



Sign in / Sign up

Export Citation Format

Share Document