Real-time modeling of transient crustal deformation through the quantification of uncertainty deduced from GNSS data

Abstract

Abstract We propose a new method for real-time uncertainty monitoring of earthquake and volcano source models using data from the global navigation satellite system (GNSS) and explore its application concerning observation station placement. The Geospatial Information Authority of Japan operates two main types of GNSS earth observation network system (GEONET) coordinates for crustal deformation monitoring on different time scales: post-processing analysis values and real-time GEONET analysis system for rapid deformation monitoring (REGARD). REGARD uses the Markov Chain Monte Carlo (MCMC) method termed real-time automatic uncertainty estimation of a coseismic single rectangular model using GNSS data (RUNE) for single rectangular fault model estimation to handle uncertainty. Thus far, no GNSS monitoring system can automatically detect transient crustal deformation events, such as volcanic activity and earthquake swarms, on timescales of a day or less. We extended RUNE and developed a core program for a new monitoring system for earthquake and volcanic source models and their uncertainties. Our program achieved automatic and stable MCMC utilization for rectangular fault, dike, Mogi, and spheroid models by increasing the computational speed, improving search efficiency, and adjusting hyperparameters. The program automatically determines the standard deviation of the likelihood function assuming a normal distribution with weights for each observation station. The calculation time was within 15 s for 1 × 106 samples on a standard 1U server. We assessed the reliability of the developed method using synthetic and observed GNSS data from the 2015 Sakurajima volcanic event. The results were consistent with the assumed model and previous studies and indicated an advantage in automatically quantifying uncertainty in a short computation time. Based on MCMC samples, we developed a new visualization algorithm to indicate areas on a map in which the number of observation stations should be expanded. We assessed the reliability using data from the 2023 Noto Peninsula earthquake [Mj 6.5]. The results indicate that the algorithm is helpful in studying the placement of stations. The above model extensions and their application are essential to achieve a rapid quantitative understanding of disaster events near urban areas and for utilizing this information in emergency response activities. Graphical Abstract

Keywords

• GNSS
• Real time
• Bayesian inversion
• Uncertainty estimation
• MCMC
• REGARD