Mapping tectonic strain in the central Alpine-Himalayan Belt with Sentinel-1 InSAR and GNSS observations

Author(s):  
Chris Rollins ◽  
Tim Wright ◽  
Jonathan Weiss ◽  
Andrew Hooper ◽  
Richard Walters ◽  
...  

<p>Geodetic measurements of crustal deformation provide crucial constraints on a region’s tectonics, geodynamics and seismic hazard. However, such geodetic constraints have traditionally been hampered by poor spatial and/or temporal sampling, which can result in ambiguities about how the lithosphere accommodates strain in space and time, and therefore where and how often earthquakes might occur. High-resolution surface deformation maps address this limitation by imaging (rather than presuming or modelling) where and how deformation takes place. These maps are now within reach for the Alpine-Himalayan Belt thanks to the COMET-LiCSAR InSAR processing system, which performs large-scale automated processing and time-series analysis of Sentinel-1 InSAR data. Expanding from our work focused on Anatolia, we are combining LiCSAR products with GNSS data to generate high-resolution maps of tectonic strain rates across the central Alpine-Himalayan Belt. Then, assuming that the buildup rate of seismic moment (deficit) from this geodetically-derived strain is balanced over the long term by the rate of moment release in earthquakes, we pair these strain rate maps with seismic catalogs to estimate the recurrence intervals of large, moderate and small earthquakes throughout the region. We also use arguments from dislocation modeling to identify two key signatures of a locked fault in a strain rate field, allowing us to convert the strain maps to “effective fault maps” and assess the contribution of individual fault systems to crustal deformation and seismic hazard. Finally, we address how to expand these approaches to the Alpine-Himalaya Belt as a whole.</p>

2020 ◽  
Author(s):  
Chris Rollins ◽  
Tim Wright ◽  
Jonathan Weiss ◽  
Andrew Hooper ◽  
Richard Walters

<p>Geodetic measurements of crustal deformation rates can provide important constraints on a region’s earthquake hazard that purely seismicity-based hazard models may miss. For example, geodesy might show that strain (or a deficit of seismic moment) is accumulating faster than the total rate at which known earthquakes have released it, implying that the long-term hazard may include larger earthquakes with long recurrence intervals (and/or temporal increases in seismicity rates). Conversely, the moment release rate in recent earthquakes might surpass the geodetic moment buildup rate, suggesting that the long-term-average earthquake activity and hazard may in fact may be more quiescent than might be estimated using the earthquake history alone. Such geodetic constraints, however, have traditionally been limited by poor spatial and/or temporal sampling, resulting in ambiguities about how the lithosphere accommodates strain in space and time that can bias estimates of the resulting hazard. High-resolution deformation maps address this limitation by imaging (rather than presuming and/or modelling) where and how deformation takes place. These maps are now within reach for the Alpine-Himalayan Belt – one of the most populous and seismically hazardous regions on Earth – thanks to the COMET-LiCSAR InSAR processing system, which performs large-scale automated processing and timeseries analysis of Sentinel-1 data provided by the EU’s Copernicus programme. We are pairing LiCSAR products with GNSS data to generate high-resolution maps of interseismic surface motion (velocity) and strain rate for the Anatolia region. Here we quantitively investigate what these strain rate distributions imply for seismic hazard in this region, using two approaches in parallel.</p><p>First, building on previous work, we develop a fully probability-based method to pair geodesy and seismic catalogs to estimate the recurrence times of large, moderate and small earthquakes in a given region. We assume that earthquakes 1) obey a power-law magnitude-frequency distribution up to a maximum magnitude and 2) collectively release seismic moment at the same rate that we estimate it is accumulating from the strain rate maps. Iterating over various magnitude-frequency distributions and their governing parameters, and formally incorporating uncertainties in moment buildup rate and the magnitudes of recorded earthquakes, we build a probabilistic long-term-average earthquake model for Anatolia as a whole, including the most likely maximum earthquake magnitude. Second, we estimate how seismic hazard may vary from place to place within Anatolia. Using insights from dislocation models, we identify two key signatures of a locked fault in a strain rate field, allowing us to convert the newly developed strain maps to “effective fault maps.” Additionally, we explore how characteristics of earthquake magnitude-frequency distributions may scale with the rate of strain (or moment) buildup, and what these scaling relations imply for the distribution of hazard in Anatolia, using the seismic catalog to evaluate these hypotheses. We also explore the implications of our findings for seismic hazard and address how to expand these approaches to the Alpine-Himalaya Belt as a whole.</p>


2016 ◽  
Vol 11 (2) ◽  
pp. 265-271 ◽  
Author(s):  
Takahiro Maeda ◽  
◽  
Hiroyuki Fujiwara

We have developed a data mining system of parallel distributed processing system which is applicable to the large-scale and high-resolution numerical simulation of ground motion by transforming into ground motion indices and their statistical values, and then visualize their values for the seismic hazard information. In this system, seismic waveforms at many locations calculated for many possible earthquake scenarios can be used as input data. The system utilizes Hadoop and it calculates the ground motion indices, such as PGV, and statistical values, such as maximum, minimum, average, and standard deviation of PGV, by parallel distributed processing with MapReduce. The computation results are being an output as GIS (Geographic Information System) data file for visualization. And this GIS data is made available via the Web Map Service (WMS). In this study, we perform two benchmark tests by applying three-component synthetic waveforms at about 80,000 locations for 10 possible scenarios of a great earthquake in Nankai Trough to our system. One is the test for PGV calculation processing. Another one is the test for PGV data mining processing. A maximum of 10 parallel processing are tested for both cases. We find that our system can hold the performance even when the total tasks is larger than 10. This system can enable us to effectively study and widely distribute to the communities for disaster mitigation since it is built with data mining and visualization for hazard information by handling a large number of data from a large-scale numerical simulation.


2020 ◽  
Author(s):  
Jonathan Weiss ◽  
Richard Walters ◽  
Tim Wright ◽  
Yu Morishita ◽  
Milan Lazecky ◽  
...  

<p>Geodetic measurements of small rates of interseismic crustal motion made at high spatial resolutions and over large areas are crucial for understanding the earthquake cycle, characterizing spatial variations in lithosphere rheology and fault frictional properties, illuminating the mechanics of large-scale continental deformation, and improving earthquake hazard models. The densification of Global Navigation Satellite System (GNSS) networks has provided an unprecedented view of the kinematics of deforming regions. However, large gaps in spatial coverage still hamper our ability to fully characterize patterns of surface deformation. Interferometric synthetic aperture radar (InSAR), and the Sentinel-1 (S-1) satellites in particular, have the potential to overcome this obstacle by providing spatially continuous measurements of surface motions, without instruments on the ground, with precision approaching that obtained from GNSS, and at a resolution of a few tens of meters. In order to manage and process the large data volumes produced by S-1, we have developed open-source, automated workflows to efficiently generate interferograms and line-of-sight (LOS) velocity fields. These outputs are valuable for a range of applications, from earthquake rapid-response to investigating human-induced ground-level changes. In this talk, we demonstrate our ability to measure plate-scale interseismic deformation using data from the first ~5 years of the S-1 mission. We estimate LOS velocities for the Anatolian microplate, an area encompassing ~800,000 km<sup>2</sup> and including the highly seismogenic North Anatolian Fault Zone (NAFZ). By combining S-1 InSAR and GNSS data, we create high-resolution surface velocity and strain rate fields for the region, which illuminate horizontal deformation patterns dominated by the westward motion of Anatolia relative to Eurasia, localized strain accumulation along the NAFZ, and rapid vertical signals associated with groundwater extraction. Relatively low levels of strain characterize other active regions including the East and Central Anatolian Fault Zones and the Western Anatolian Extensional Province. We find that GNSS data alone are insufficient for characterizing key details of the strain rate field that are critical for understanding the relationship between strain accumulation and release in earthquakes. We highlight two important results stemming from our work including probabilistic estimates of the recurrence times of earthquakes of varying magnitudes for the region and a new NAFZ locking distribution that shows close correspondence to the surface rupture extents of large 20<sup>th</sup> century earthquakes.</p>


Author(s):  
A. V. Crewe ◽  
M. Ohtsuki

We have assembled an image processing system for use with our high resolution STEM for the particular purpose of working with low dose images of biological specimens. The system is quite flexible, however, and can be used for a wide variety of images.The original images are stored on magnetic tape at the microscope using the digitized signals from the detectors. For low dose imaging, these are “first scan” exposures using an automatic montage system. One Nova minicomputer and one tape drive are dedicated to this task.The principal component of the image analysis system is a Lexidata 3400 frame store memory. This memory is arranged in a 640 x 512 x 16 bit configuration. Images are displayed simultaneously on two high resolution monitors, one color and one black and white. Interaction with the memory is obtained using a Nova 4 (32K) computer and a trackball and switch unit provided by Lexidata.The language used is BASIC and uses a variety of assembly language Calls, some provided by Lexidata, but the majority written by students (D. Kopf and N. Townes).


2014 ◽  
Vol 31 (2) ◽  
Author(s):  
Mariela Gabioux ◽  
Vladimir Santos da Costa ◽  
Joao Marcos Azevedo Correia de Souza ◽  
Bruna Faria de Oliveira ◽  
Afonso De Moraes Paiva

Results of the basic model configuration of the REMO project, a Brazilian approach towards operational oceanography, are discussed. This configuration consists basically of a high-resolution eddy-resolving, 1/12 degree model for the Metarea V, nested in a medium-resolution eddy-permitting, 1/4 degree model of the Atlantic Ocean. These simulations performed with HYCOM model, aim for: a) creating a basic set-up for implementation of assimilation techniques leading to ocean prediction; b) the development of hydrodynamics bases for environmental studies; c) providing boundary conditions for regional domains with increased resolution. The 1/4 degree simulation was able to simulate realistic equatorial and south Atlantic large scale circulation, both the wind-driven and the thermohaline components. The high resolution simulation was able to generate mesoscale and represent well the variability pattern within the Metarea V domain. The BC mean transport values were well represented in the southwestern region (between Vitória-Trinidade sea mount and 29S), in contrast to higher latitudes (higher than 30S) where it was slightly underestimated. Important issues for the simulation of the South Atlantic with high resolution are discussed, like the ideal place for boundaries, improvements in the bathymetric representation and the control of bias SST, by the introducing of a small surface relaxation. In order to make a preliminary assessment of the model behavior when submitted to data assimilation, the Cooper & Haines (1996) method was used to extrapolate SSH anomalies fields to deeper layers every 7 days, with encouraging results.


Sign in / Sign up

Export Citation Format

Share Document