scholarly journals Adjoint-based direct data assimilation of GNSS time series for optimizing frictional parameters and predicting postseismic deformation following the 2003 Tokachi-oki earthquake

2020 ◽  
Vol 72 (1) ◽  
Author(s):  
Masayuki Kano ◽  
Shin’ichi Miyazaki ◽  
Yoichi Ishikawa ◽  
Kazuro Hirahara

Abstract Postseismic Global Navigation Satellite System (GNSS) time series followed by megathrust earthquakes can be interpreted as a result of afterslip on the plate interface, especially in its early phase. Afterslip is a stress release process accumulated by adjacent coseismic slip and can be considered a recovery process for future events during earthquake cycles. Spatio-temporal evolution of afterslip often triggers subsequent earthquakes through stress perturbation. Therefore, it is important to quantitatively capture the spatio-temporal evolution of afterslip and related postseismic crustal deformation and to predict their future evolution with a physics-based simulation. We developed an adjoint data assimilation method, which directly assimilates GNSS time series into a physics-based model to optimize the frictional parameters that control the slip behavior on the fault. The developed method was validated with synthetic data. Through the optimization of frictional parameters, the spatial distributions of afterslip could roughly (but not in detail) be reproduced if the observation noise was included. The optimization of frictional parameters reproduced not only the postseismic displacements used for the assimilation, but also improved the prediction skill of the following time series. Then, we applied the developed method to the observed GNSS time series for the first 15 days following the 2003 Tokachi-oki earthquake. The frictional parameters in the afterslip regions were optimized to A–B ~ O(10 kPa), A ~ O(100 kPa), and L ~ O(10 mm). A large afterslip is inferred on the shallower side of the coseismic slip area. The optimized frictional parameters quantitatively predicted the postseismic GNSS time series for the following 15 days. These characteristics can also be detected if the simulation variables can be simultaneously optimized. The developed data assimilation method, which can be directly applied to GNSS time series following megathrust earthquakes, is an effective quantitative evaluation method for assessing risks of subsequent earthquakes and for monitoring the recovery process of megathrust earthquakes.

2020 ◽  
Author(s):  
Masayuki Kano ◽  
Shin'ichi Miyazaki ◽  
Yoichi Ishikawa ◽  
Kazuro Hirahara

Abstract PostseismicGlobal Navigation Satellite System (GNSS) time seriesfollowed by megathrust earthquakes can be interpreted as a result of afterslip on the plate interface especially in its early phase. Afterslip is a stress release process accumulated by adjacent coseismic slip andcan be considered a recovery process for future events during earthquake cycles. Spatio-temporal evolution of afterslip often triggers subsequent earthquakes through stress perturbation. Therefore, it is important toquantitativelycapture the spatio-temporal evolution of afterslip and related postseismic crustal deformation and to predict their future evolution with a physics-based simulation. We developedanadjoint data assimilation method, which directly assimilates GNSS time series into a physics-based model to optimize the frictional parameters that control the slip behavior on the fault.The developed method was validated with synthetic data. Through the optimization of frictional parameters, the spatial distributions of afterslip can be roughly reproduced but not in detail if the observation noise is included. The optimization of frictional parameters provides not only the reproduction ofpostseismic displacements used for the assimilation but also the improvement in the prediction skill of the following time series. Then, we appliedthe developed method to the observed GNSS time series for the first 15 d following the 2003 Tokachi-oki earthquake. The frictional parameters in the afterslip regions were optimized to A-B ~ O(10 kPa), A ~ O(100 kPa), and L ~ O(10 mm). The large afterslip is inferred on the shallower side of the coseismic slip area. The optimized frictional parameters quantitatively predicted the postseismicGNSS time series for the following 15 d. These characteristics can be also detected if the simulation variables can be simultaneously optimized. The developed data assimilation method, which can be directly applied to GNSS time series following megathrust earthquakes, isan effective quantitative evaluation method for assessing risks of subsequent earthquakes and for monitoring the recovery process of megathrust earthquakes.


2020 ◽  
Author(s):  
Masayuki Kano ◽  
Shin'ichi Miyazaki ◽  
Yoichi Ishikawa ◽  
Kazuro Hirahara

Abstract Postseismic Global Navigation Satellite System (GNSS) time series followed by megathrust earthquakes can be interpreted as a result of afterslip on the plate interface, especially in its early phase. Afterslip is a stress release process accumulated by adjacent coseismic slip and can be considered a recovery process for future events during earthquake cycles. Spatio-temporal evolution of afterslip often triggers subsequent earthquakes through stress perturbation. Therefore, it is important to quantitatively capture the spatio-temporal evolution of afterslip and related postseismic crustal deformation and to predict their future evolution with a physics-based simulation. We developed an adjoint data assimilation method, which directly assimilates GNSS time series into a physics-based model to optimize the frictional parameters that control the slip behavior on the fault. The developed method was validated with synthetic data. Through the optimization of frictional parameters, the spatial distributions of afterslip could roughly (but not in detail) be reproduced if the observation noise was included. The optimization of frictional parameters reproduced not only the postseismic displacements used for the assimilation but also improved the prediction skill of the following time series. Then, we applied the developed method to the observed GNSS time series for the first 15 days following the 2003 Tokachi-oki earthquake. The frictional parameters in the afterslip regions were optimized to A-B ~ O(10 kPa), A ~ O(100 kPa), and L ~ O(10 mm). A large afterslip is inferred on the shallower side of the coseismic slip area. The optimized frictional parameters quantitatively predicted the postseismic GNSS time series for the following 15 days. These characteristics can also be detected if the simulation variables can be simultaneously optimized. The developed data assimilation method, which can be directly applied to GNSS time series following megathrust earthquakes, is an effective quantitative evaluation method for assessing risks of subsequent earthquakes and for monitoring the recovery process of megathrust earthquakes.


2020 ◽  
Author(s):  
Masayuki Kano ◽  
Shin'ichi Miyazaki ◽  
Yoichi Ishikawa ◽  
Kazuro Hirahara

Abstract Postseismic Global Navigation Satellite System (GNSS) time series followed by megathrust earthquakes can be interpreted as a result of afterslip on the plate interface, especially in its early phase. Afterslip is a stress release process accumulated by adjacent coseismic slip and can be considered a recovery process for future events during earthquake cycles. Spatio-temporal evolution of afterslip often triggers subsequent earthquakes through stress perturbation. Therefore, it is important to quantitatively capture the spatio-temporal evolution of afterslip and related postseismic crustal deformation and to predict their future evolution with a physics-based simulation. We developed an adjoint data assimilation method, which directly assimilates GNSS time series into a physics-based model to optimize the frictional parameters that control the slip behavior on the fault. The developed method was validated with synthetic data. Through the optimization of frictional parameters, the spatial distributions of afterslip could roughly (but not in detail) be reproduced if the observation noise was included. The optimization of frictional parameters reproduced not only the postseismic displacements used for the assimilation but also improved the prediction skill of the following time series. Then, we applied the developed method to the observed GNSS time series for the first 15 days following the 2003 Tokachi-oki earthquake. The frictional parameters in the afterslip regions were optimized to A-B ~ O (10 kPa), A ~ O (100 kPa), and L ~ O (10 mm). A large afterslip is inferred on the shallower side of the coseismic slip area. The optimized frictional parameters quantitatively predicted the postseismic GNSS time series for the following 15 days. These characteristics can also be detected if the simulation variables can be simultaneously optimized. The developed data assimilation method, which can be directly applied to GNSS time series following megathrust earthquakes, is an effective quantitative evaluation method for assessing risks of subsequent earthquakes and for monitoring the recovery process of megathrust earthquakes.


Author(s):  
M. Béjar-Pizarro ◽  
P. Ezquerro Martín ◽  
G. Herrera ◽  
R. Tomás ◽  
C. Guardiola-Albert ◽  
...  

Abstract. The Tertiary detritic aquifer of Madrid (TDAM), with an average thickness of 1500 m and a heterogeneous, anisotropic structure, supplies water to Madrid, the most populated city of Spain (3.2 million inhabitants in the metropolitan area). Besides its complex structure, a previous work focused in the north-northwest of Madrid city showed that the aquifer behaves quasi elastically trough extraction/recovery cycles and ground uplifting during recovery periods compensates most of the ground subsidence measured during previous extraction periods (Ezquerro et al., 2014). Therefore, the relationship between ground deformation and groundwater level through time can be simulated using simple elastic models. In this work, we model the temporal evolution of the piezometric level in 19 wells of the TDAM in the period 1997–2010. Using InSAR and piezometric time series spanning the studied period, we first estimate the elastic storage coefficient (Ske) for every well. Both, the Ske of each well and the average Ske of all wells, are used to predict hydraulic heads at the different well locations during the study period and compared against the measured hydraulic heads, leading to very similar errors when using the Ske of each well and the average Ske of all wells: 14 and 16 % on average respectively. This result suggests that an average Ske can be used to estimate piezometric level variations in all the points where ground deformation has been measured by InSAR, thus allowing production of piezometric level maps for the different extraction/recovery cycles in the TDAM.


2020 ◽  
Vol 12 (17) ◽  
pp. 2850 ◽  
Author(s):  
Simon Daout ◽  
Andreas Steinberg ◽  
Marius Paul Isken ◽  
Sebastian Heimann ◽  
Henriette Sudhaus

Inferring the geometry and evolution of an earthquake sequence is crucial to understand how fault systems are segmented and interact. However, structural geological models are often poorly constrained in remote areas and fault inference is an ill-posed problem with a reliability that depends on many factors. Here, we investigate the geometry of the Mw 6.3 2008 and 2009 Qaidam earthquakes, in northeast Tibet, by combining InSAR time series and teleseismic data. We conduct a multi-array back-projection analysis from broadband teleseismic data and process three overlapping Envisat tracks covering the two earthquakes to extract the spatio-temporal evolution of seismic ruptures. We then integrate both geodetic and seismological data into a self-consistent kinematic model of the earthquake sequence. Our results constrain the depth and along-strike segmentation of the thrust-faulting sequence. The 2008 earthquake ruptured a ∼32° north-dipping fault that roots under the Olongbulak pop-up structure at ∼12 km depth and fault slip evolved post-seismically in a downdip direction. The 2009 earthquake ruptured three south-dipping high-angle thrusts and propagated from ∼9 km depth to the surface and bilaterally along the south-dipping segmented 55–75° high-angle faults of the Olonbulak pop-up structure that displace basin deformed sedimentary sequences above Paleozoic bedrock. Our analysis reveals that the inclusion of the post-seismic afterslip into modelling is beneficial in the determination of fault geometry, while teleseismic back-projection appears to be a robust tool for identifying rupture segmentation for moderate-sized earthquakes. These findings support the hypothesis that the Qilian Shan is expanding southward along a low-angle décollement that partitions the oblique convergence along multiple flower and pop-up structures.


Sensors ◽  
2020 ◽  
Vol 20 (3) ◽  
pp. 946 ◽  
Author(s):  
Jin Zhang ◽  
Cheng Wu ◽  
Yiming Wang

Abnormal falls in public places have significant safety hazards and can easily lead to serious consequences, such as trampling by people. Vision-driven fall event detection has the huge advantage of being non-invasive. However, in actual scenes, the fall behavior is rich in diversity, resulting in strong instability in detection. Based on the study of the stability of human body dynamics, the article proposes a new model of human posture representation of fall behavior, called the “five-point inverted pendulum model”, and uses an improved two-branch multi-stage convolutional neural network (M-CNN) to extract and construct the inverted pendulum structure of human posture in real-world complex scenes. Furthermore, we consider the continuity of the fall event in time series, use multimedia analytics to observe the time series changes of human inverted pendulum structure, and construct a spatio-temporal evolution map of human posture movement. Finally, based on the integrated results of computer vision and multimedia analytics, we reveal the visual characteristics of the spatio-temporal evolution of human posture under the potentially unstable state, and explore two key features of human fall behavior: motion rotational energy and generalized force of motion. The experimental results in actual scenes show that the method has strong robustness, wide universality, and high detection accuracy.


2021 ◽  
Vol 7 (5) ◽  
pp. 938-949
Author(s):  
Miao Yahui ◽  
Jiang Ce ◽  
Gu Zifan ◽  
Xi Zenglei

In recent years, the world tobacco industry has shown a trend of continuous growth in market demand and increasingly concentrated tobacco market, which brings new opportunities to the tobacco agriculture development. However, the modern tobacco agriculture is bound to occupy a large amount of agricultural land, which will have a certain impact on the distribution of urban construction land (UCL). Therefore, reasonable and effective modern tobacco agricultural planning plays a vital role in promoting modern agricultural and urban development.In this paper, we took Baiyangdian basin as the study area, applying the long time series nighttime light (NTL) data to extract UCL and analysing the impact of modern tobacco agriculture planning on the spatio-temporal evolution of UCL. Firstly, we used a power function model to fit the two kinds of international mainstream NTL data to form a long time series NTL dataset from 1992 to 2018 to make the NTL data comparable over a long time. Then, the threshold segmentation method based on land use data calibration was applied to extract the UCL in Baiyangdian basin in 2000, 2010 and 2018, and to analyze its spatio-temporal evolution characteristics and patterns by combining the landscape metrics and gravity model. The results show that the UCL in Baiyangdian basin shows the expansion trend centered on the main urban area of Baoding City, land intensification degree has increased, and the modern tobacco agriculture planning has a profound impact on the spatio-temporal distribution of UCL. Our study will provide technical support and experience for the scientific modern tobacco agriculture planning.


Sign in / Sign up

Export Citation Format

Share Document