The role of space-based observation in understanding and responding to active tectonics and earthquakes

On the Patterns and Scaling Properties of the 2021–2022 Arkalochori Earthquake Sequence (Central Crete, Greece) Based on Seismological, Geophysical and Satellite Observations

The 27 September 2021 damaging mainshock (Mw6.0) close to Arkalochori village is the strongest earthquake that was recorded during the instrumental period of seismicity in Central Crete (Greece). The mainshock was preceded by a significant number of foreshocks that lasted nearly four months. Maximum ground subsidence of about 18 cm was estimated from InSAR processing. The aftershock sequence is located in an almost NE-SW direction and divided into two main clusters, the southern and the northern ones. The foreshock activity, the deformation area, and the strongest aftershocks are located within the southern cluster. Based on body-wave travel times, a 3-D velocity model was developed, while using combined space and ground-based geodetic techniques, the co-seismic ground deformation is presented. Moreover, we examined the co-seismic static stress changes with respect to the aftershocks’ spatial distribution during the major events of the foreshocks, the Mw = 6.0 main event as well as the largest aftershock. Both the foreshock and the aftershock sequences obey the scaling law for the frequency-magnitude distribution as derived from the framework of non-extensive statistical physics (NESP). The aftershock production rate decays according to the modified Omori scaling law, exhibiting various Omori regimes due to the generation of secondary aftershock sequences. The analysis of the inter-event time distribution, based on NESP, further indicates asymptotic power-law scaling and long-range correlations among the events. The spatiotemporal evolution of the aftershock sequence indicates triggering by co-seismic stress transfer, while its slow migration towards the outer edges of the area of the aftershocks, related to the logarithm of time, further indicates a possible afterslip.

Satellite-Observed Thermal Anomalies and Deformation Patterns Associated to the 2021, Central Crete Seismic Sequence

Nowadays, there has been a growing interest in understanding earthquake forerunners, i.e., anomalous variations that are possibly associated with the complex process of earthquake evolution. In this context, the Robust Satellite Technique was coupled with 10 years (2012–2021) of daily night-time MODIS-Land Surface Temperature remote sensing data to detect thermal anomalies likely related to the 27 September 2021, strong onshore earthquake of magnitude Mw6.0 occurring near the Arkalochori village in Central Crete, Greece. Eight intense (signal-to-noise ratio > 3) and infrequent, quite extensive, and temporally persistent thermal signal transients were detected and characterized as pre-seismic anomalies, while one thermal signal transient was identified as a co-seismic effect on the day of the main tectonic event. The thermal anomalies dataset was combined with tectonic parameters of Central Crete, such as active faults and fault density, seismogenic zones and ground displacement maps produced using Sentinel-1 satellite imagery and the Interferometric Synthetic Aperture Radar technique. Regarding the thermal anomaly of 27 September, its greatest portion was observed over the footwall part of the fault where a significant subsidence up to 20 cm exists. We suggest that the thermal anomalies are possibly connected with gas release which happens due to stress changes and is controlled by the existence of tectonic lines and the density of the faults, even if alternative explanations could not be excluded.

Ground Displacements Estimation through GNSS and Geometric Leveling: A Geological Interpretation of the 2016–2017 Seismic Sequence in Central Italy

Between August 2016 and January 2017, a very energetic seismic sequence induced substantial horizontal and vertical ground displacements in the Central Italian Apennines. After this event, the Italian Military Geographical Institute (IGM), owner and manager of the Italian geodetic networks, executed several topographic surveys in the earthquake area in order to update the coordinates of vertices belonging to the IGM95 geodetic network. The measurements began in the areas where the most significant deformation occurred: the localities of Amatrice and Accumoli, in the Rieti Province, and the area covering Norcia and Castelluccio, in the Province of Perugia, all the way to Visso (Province of Macerata). The activities described in this paper focused on the updated measurement of the IGM95 network points through GNSS and the restatement of extensive parts of the high precision geometric lines that were levelled until reaching stable zones. This unprecedented amount of data was used for a new geological interpretation of the seismic sequence, which confirms some of the previous hypotheses of the scientific community. In the analyzed territory, the latest estimate of the geodetic position points has allowed for an accurate determination of the east and the north and of the altitude components of the displacement induced by the earthquake through a comparison with the previous coordinates. The results confirm that the seismicity was induced by normal faults system activity. Still, they also indicate the possible influence of a significant regional thrust that conditioned the propagation of the seismicity in the area. The obtained maps of the displacement are coherent with other geodetic works and with a rupture propagation driven by the documented geotectonic structure.

A New Method for InSAR Stratified Tropospheric Delay Correction Facilitating Refinement of Coseismic Displacement Fields of Small-to-Moderate Earthquakes

Focusing on stratified tropospheric delay correction in the small-amplitude coseismic displacement field of small-to-moderate earthquakes (<Mw 6.5), we develop a Simple-Stratification-Correction (SSC) approach based on the empirical phase-elevation relationship and spatial properties of the troposphere, via an equal-size window segmentation. We validate our SSC method using 23 real earthquakes that occurred from January 2016 to May 2021 with a moment magnitude (Mw) ranging from 4.5 to 6.5. We conclude that SSC performs well according to the amount of reduction in semi-variance and the root-mean-square value. This method primarily focuses on stratification delay correction; thus, it is especially useful in regions with complex terrain, while it can mitigate partial large-scale turbulence signals. We investigate three parameters that are empirically setup in the correction working flow and inspect their optimal settings, when implementing SSC for quick response after earthquake. Our method is ready to be integrated into an operational InSAR processing chain to produce a reliable atmospheric phase screen map, which can also serve as an auxiliary product to quickly and timely quantify stratification delays in coseismic interferograms. Through improved accuracy of the coseismic displacement field, the focal mechanism could be better constrained to facilitate the building and expansion of the geodesy-based earthquake catalogue.

The Campo de Dalias GNSS Network Unveils the Interaction between Roll-Back and Indentation Tectonics in the Gibraltar Arc

The Gibraltar Arc includes the Betic and Rif Cordilleras surrounding the Alboran Sea; it is formed at the northwest–southeast Eurasia–Nubia convergent plate boundary in the westernmost Mediterranean. Since 2006, the Campo de Dalias GNSS network has monitored active tectonic deformation of the most seismically active area on the north coast of the Alboran Sea. Our results show that the residual deformation rates with respect to Eurasia range from 1.7 to 3.0 mm/year; roughly homogenous west-southwestward displacements of the northern sites occur, while the southern sites evidence irregular displacements towards the west and northwest. This deformation pattern supports simultaneous east-northeast–west-southwest extension, accommodated by normal and oblique faults, and north-northwest–south-southeast shortening that develops east-northeast–west-southwest folds. Moreover, the GNSS results point to dextral creep of the main northwest–southeast Balanegra Fault. These GNNS results thus reveal, for the first time, present-day interaction of the roll-back tectonics of the Rif–Gibraltar–Betic slab in the western part of the Gibraltar Arc with the indentation tectonics affecting the eastern and southern areas, providing new insights for improving tectonic models of arcuate orogens.

Assessing the Use of Optical Satellite Images to Detect Volcanic Impacts on Glacier Surface Morphology

Globally, about 250 Holocene volcanoes are either glacier-clad or have glaciers in close proximity. Interactions between volcanoes and glaciers are therefore common, and some of the most deadly (e.g., Nevado del Ruiz, 1985) and most costly (e.g., Eyjafjallajökull, 2010) eruptions of recent years were associated with glaciovolcanism. An improved understanding of volcano-glacier interactions is therefore of both global scientific and societal importance. This study investigates the potential of using optical satellite images to detect volcanic impacts on glaciers, with a view to utilise detected changes in glacier surface morphology to improve glacier-clad volcano monitoring and eruption forecasting. Roughly 1400 optical satellite images are investigated from key, well-documented eruptions around the globe during the satellite remote sensing era (i.e., 1972 to present). The most common observable volcanic impact on glacier morphology (for both thick and thin ice-masses) is the formation of ice cauldrons and openings, often associated with concentric crevassing. Other observable volcanic impacts include ice bulging and fracturing due to subglacial dome growth; localized crevassing adjacent to supraglacial lava flows; widespread glacier crevassing, presumably, due to meltwater-triggered glacier acceleration and advance. The main limitation of using optical satellite images to investigate changes in glacier morphology is the availability of cloud- and eruption-plume-free scenes of sufficient spatial- and temporal resolution. Therefore, for optimal monitoring and eruption prediction at glacier-clad volcanoes, optical satellite images are best used in combination with other sources, including SAR satellite data, aerial images, ground-based observations and satellite-derived products (e.g., DEMs).

Applications of Satellite Radar Imagery for Hazard Monitoring: Insights from Australia

Earth observation (EO) satellites facilitate hazard monitoring and mapping over large-scale and remote areas. Despite Synthetic Aperture Radar (SAR) satellites being well-documented as a hazard monitoring tool, the uptake of these data is geographically variable, with the Australian continent being one example where the use of SAR data is limited. Consequently, less is known about how these data apply in the Australian context, how they could aid national hazard monitoring and assessment, and what new insights could be gleaned for the benefit of the international disaster risk reduction community. The European Space Agency Sentinel-1 satellite mission now provides the first spatially and temporally complete global SAR dataset and the first opportunity to use these data to systematically assess hazards in new locations. Using the example of Australia, where floods and uncontrolled bushfires, earthquakes, resource extraction (groundwater, mining, hydrocarbons) and geomorphological changes each pose potential risks to communities, we review past usage of EO for hazard monitoring and present a suite of new case studies that demonstrate the potential added benefits of SAR. The outcomes provide a baseline understanding of the potential role of SAR in national hazard monitoring and assessment in an Australian context. Future opportunities to improve national hazard identification will arise from: new SAR sensing capabilities, which for Australia includes a first-ever civilian EO capability, NovaSAR-1; the integration of Sentinel-1 SAR with other EO datasets; and the provision of standardised SAR products via Analysis Ready Data and Open Data Cubes to support operational applications.

Complete Atomic Oxygen and UV Protection for Polymer and Composite Materials in a Low Earth Orbit

With the realization of larger and more complex space installations, an increase in the surface area exposed to atomic oxygen (AO) and ultraviolet (UV) effects is expected, making structural integrity of space structures essential for future development. In a low Earth orbit (LEO), the effects of AO and UV degradation can have devastating consequences for polymer and composite structures in satellites and space installations. Composite materials such as carbon fiber-reinforced polymer (CFRP) or polymer materials such as polyetherimide and polystyrene are widely used in satellite construction for various applications including structural components, thermal insulation, and importantly radio frequency (RF) assemblies. In this paper, we present a multilayered material protection solution, a multilayered protection barrier, that mitigates the effects of AO and UV without disrupting the functional performance of tested assemblies. This multilayered protection barrier deposited via a custom-built plasma-enhanced chemical vapor deposition (PECVD) system is designed so as to deposit all necessary layers without breaking vacuum to maximize the adhesion to the surface of the substrate and to ensure no pinhole erosion is present. In the multilayer solution, a moisture and outgassing barrier (MOB) is coupled with an AO and UV capping layer to provide complete protection.

PSI Clustering for the Assessment of Underground Infrastructure Deterioration

Remote sensing images find application in several different domains, such as land cover or land usage observation, environmental monitoring, and urbanization. This latter field has recently witnessed an interesting development with the use of remote sensing for infrastructural monitoring. In this work, we present an analysis of Sentinel-1 images, which were used to monitor the Italian provinces of Bologna and Modena located at the Emilia Region Apennines foothill. The goal of this study was the development of a machine learning-based detection system to monitor the deterioration of public aqueduct infrastructures based on Persistent Scatterer Interferometry (PSI). We evaluated the land deformation over a temporal range of five years; these series feed a k-means clustering algorithm to separate the pixels of the region according to different deformation patterns. Furthermore, we defined the critical areas as those areas where different patterns collided or overlapped. The proposed approach provides an informative tool for the structural health monitoring of underground infrastructures.

Evaluating Potential Ground Subsidence Geo-Hazard of Xiamen Xiang’an New Airport on Reclaimed Land by SAR Interferometry

The land reclaimed from the seaside may have a long-term subsidence trend, which poses a potential geohazard in the future land use. Xiamen Xiang’an New Airport (XXNA) is built on reclaimed land since 2016. Based on the spaceborne Sentinel-1 data between January 2018 to April 2019 and the time series interferometric synthetic aperture radar (InSAR) technique, this paper analyzed the reclaimed land subsidence evolution at XXNA in this period. InSAR measurements show that XXNA is suffering from severe subsidence, mainly in three regions because of the earth and sand compacting. By analyzing the spatial subsidence characterizations of the main subsiding areas combined with historical land reclamation and future land use planning, we find the potential threat of subsidence to future land use. Correlation between subsidence and the period of reclamation was found, indicating that the consolidation and compression in dredger fill is the main cause of subsidence. By combining subsidence monitoring results with different land use types and adopting the Expectation (Ex) and Entropy (En) methods, we analyzed the key area with potential subsidence geo-hazard. This work shows that with SAR interferometry, it is possible to find the large area ground subsidence in the airport reclaimed area. The areas with potential subsidence geo-hazards are consistent with the deep reclaimed earth, which means high subsidence risk in the future.

LiCSAR: An Automatic InSAR Tool for Measuring and Monitoring Tectonic and Volcanic Activity

Space-borne Synthetic Aperture Radar (SAR) Interferometry (InSAR) is now a key geophysical tool for surface deformation studies. The European Commission’s Sentinel-1 Constellation began acquiring data systematically in late 2014. The data, which are free and open access, have global coverage at moderate resolution with a 6 or 12-day revisit, enabling researchers to investigate large-scale surface deformation systematically through time. However, full exploitation of the potential of Sentinel-1 requires specific processing approaches as well as the efficient use of modern computing and data storage facilities. Here we present Looking Into Continents from Space with Synthetic Aperture Radar (LiCSAR), an operational system built for large-scale interferometric processing of Sentinel-1 data. LiCSAR is designed to automatically produce geocoded wrapped and unwrapped interferograms and coherence estimates, for large regions, at 0.001° resolution (WGS-84 coordinate system). The products are continuously updated at a frequency depending on prioritised regions (monthly, weekly or live update strategy). The products are open and freely accessible and downloadable through an online portal. We describe the algorithms, processing, and storage solutions implemented in LiCSAR, and show several case studies that use LiCSAR products to measure tectonic and volcanic deformation. We aim to accelerate the uptake of InSAR data by researchers as well as non-expert users by mass producing interferograms and derived products.

Spatial forecasting of seismicity provided from Earth observation by space satellite technology

Understanding the controls on the distribution and magnitude of earthquakes is required for effective earthquake forecasting. We present a study that demonstrates that the distribution and size of earthquakes in Italy correlates with the steady state rate at which the Earth’s crust moves. We use a new high-resolution horizontal strain rate (S) field determined from a very dense velocity field derived from the combination of Global Navigation Satellite System (GNSS) and satellite radar interferometry from two decades of observations. Through a statistical approach we study the correlation between the S and the magnitude of M ≥ 2.5 earthquakes that occurred in the same period of satellite observations. We found that the probability of earthquakes occurring is linked to S by a linear correlation, and more specifically the probability that a strong seismic event occurs doubles with the doubling of S. It also means that lower horizontal strain rate zone can have as large earthquakes as high horizontal strain rate zones, just with a reduced probability. The work demonstrates an independent and quantitative tool to spatially forecast seismicity.

LiCSBAS: An Open-Source InSAR Time Series Analysis Package Integrated with the LiCSAR Automated Sentinel-1 InSAR Processor

For the past five years, the 2-satellite Sentinel-1 constellation has provided abundant and useful Synthetic Aperture Radar (SAR) data, which have the potential to reveal global ground surface deformation at high spatial and temporal resolutions. However, for most users, fully exploiting the large amount of associated data is challenging, especially over wide areas. To help address this challenge, we have developed LiCSBAS, an open-source SAR interferometry (InSAR) time series analysis package that integrates with the automated Sentinel-1 InSAR processor (LiCSAR). LiCSBAS utilizes freely available LiCSAR products, and users can save processing time and disk space while obtaining the results of InSAR time series analysis. In the LiCSBAS processing scheme, interferograms with many unwrapping errors are automatically identified by loop closure and removed. Reliable time series and velocities are derived with the aid of masking using several noise indices. The easy implementation of atmospheric corrections to reduce noise is achieved with the Generic Atmospheric Correction Online Service for InSAR (GACOS). Using case studies in southern Tohoku and the Echigo Plain, Japan, we demonstrate that LiCSBAS applied to LiCSAR products can detect both large-scale (>100 km) and localized (~km) relative displacements with an accuracy of <1 cm/epoch and ~2 mm/yr. We detect displacements with different temporal characteristics, including linear, periodic, and episodic, in Niigata, Ojiya, and Sanjo City, respectively. LiCSBAS and LiCSAR products facilitate greater exploitation of globally available and abundant SAR datasets and enhance their applications for scientific research and societal benefit.

A mission planning method for multi-satellite wide area observation

Space-based Earth observation is now developing rapidly because of the advantages in coverage, convenience, and flexibility. However, it is difficult for one single satellite to realize the quick observation of wide areas. In order to catch all the significant information of a wide area in a short time, multi-satellite observation mission would be proposed. In this article, a mission planning method for online multi-satellite wide area observation is established to serve future multi-satellite observation missions. Firstly, a method for area division is proposed, and the whole area is divided into subareas. Then the multi-satellite observation path planning is realized by a strategy of path deduction. After that, a remaining time allocation method to maximize the observation gazing time of each subarea is proposed. Finally, the algorithms for the whole mission planning process are provided. Numerical simulations show that the mission planning method is able to ensure the complete coverage of different wide target areas, with high reliability and low computational complexity.

Analogy model test by experimental simulation of the subsidence zone associated with the 3 September 2017 North Korean nuclear test

The surveillance of nuclear tests relies on the seismic network and the interferometric synthetic aperture radar (InSAR) data. Using the surface displacements such as subsidence craters caused by a nuclear test is an effective tool to monitor and estimate source characteristics of nuclear tests. In this study, based on the scaling laws of the process of subsidence crater formation, we developed an explosive model test in a vacuum chamber using the yield and buried depth estimated by Wang et al. to simulate the surface subsidence zone and collapse the crater associated with the 3 September 2017 North Korean nuclear test, and compared with their research results. The explosive simulation apparatus independently developed by the Army Engineering University of PLA is also introduced. The simulated results indicate that the radii of the subsidence zone and collapse crater are ∼257 m and ∼154 m, respectively, which are close to the empirical formula and InSAR observations. The explosive model test in a vacuum chamber could help to characterize the surface displacements induced by nuclear tests and provide an effective method for nuclear monitoring and characterization around the world.

The State of Remote Sensing Capabilities of Cascading Hazards over High Mountain Asia

Cascading hazard processes refer to a primary trigger such as heavy rainfall, seismic activity, or snow melt, followed by a chain or web of consequences that can cause subsequent hazards influenced by a complex array of preconditions and vulnerabilities. These interact in multiple ways and can have tremendous impacts on populations proximate to or downstream of these initial triggers. High Mountain Asia (HMA) is extremely vulnerable to cascading hazard processes given the tectonic, geomorphologic, and climatic setting of the region, particularly as it relates to glacial lakes. Given the limitations of in situ surveys in steep and often inaccessible terrain, remote sensing data are a valuable resource for better understanding and quantifying these processes. The present work provides a survey of cascading hazard processes impacting HMA and how these can be characterized using remote sensing sources. We discuss how remote sensing products can be used to address these process chains, citing several examples of cascading hazard scenarios across HMA. This work also provides a perspective on the current gaps and challenges, community needs, and view forward towards improved characterization of evolving hazards and risk across HMA.

The Spatial and Temporal Influence of Cloud Cover on Satellite-Based Emergency Mapping of Earthquake Disasters

The ability to rapidly access optical satellite imagery is now an intrinsic component of managing the disaster response that follows a major earthquake. These images provide synoptic data on the impacts, extent, and intensity of damage, which is essential for mitigating further losses by feeding into the response coordination. However, whilst the efficiency of the response can be hampered when cloud cover limits image availability, spatio-temporal variations in cloud cover have never been considered as part of the design of effective disaster mapping. Here we show how annual variations in cloud cover may affect our capacity to respond rapidly throughout the year and consequently contribute to overall earthquake risk. We find that on a global scale when accounting for cloud, the worst time of year for an earthquake disaster is between June and August. During these months, 40% of the global population at risk from earthquakes are obscured from optical satellite view for >3 consecutive days. Southeastern Asia is particularly strongly affected, accounting for the majority of the population at risk from earthquakes that could be obscured by cloud in every month. Our results demonstrate the importance of the timing of earthquakes in terms of our capacity to respond effectively, highlighting the need for more intelligent design of disaster response that is not overly reliant on optical satellite imagery.

On the Segmentation of the Cephalonia–Lefkada Transform Fault Zone (Greece) from an InSAR Multi-Mode Dataset of the Lefkada 2015 Sequence

We use Interferometric Synthetic Aperture Radar (InSAR) to study the Cephalonia–Lefkada Transform Fault Zone (CTF) in the Ionian Sea. The CTF separates continental subduction to the north from oceanic subduction to the south, along the Hellenic Subduction Zone. We exploit a rich multi-modal radar dataset of the most recent major earthquake in the region, the 17 November 2015 Mw 6.4 event, and present new surface displacement results that offer additional constraints on the fault segmentation of the area. Based on this dataset, and by exploiting available information of earthquake relocation, we propose a new rupture process for the 2015 sequence, complementary to those published already. Our modelling includes an additional southern fault segment, oblique to the segment related with the mainshock, which indicates that the CTF structure is more complex than previously believed.

Earthquake Environmental Effects of the 1992 MS7.3 Suusamyr Earthquake, Kyrgyzstan, and Their Implications for Paleo-Earthquake Studies

Large pre-historical earthquakes leave traces in the geological and geomorphological record, such as primary and secondary surface ruptures and mass movements, which are the only means to estimate their magnitudes. These environmental earthquake effects (EEEs) can be calibrated using recent seismic events and the Environmental Seismic Intensity Scale (ESI2007). We apply the ESI2007 scale to the 1992 MS7.3 Suusamyr Earthquake in the Kyrgyz Tien Shan, because similar studies are sparse in that area and geological setting, and because this earthquake was very peculiar in its primary surface rupture pattern. We analyze literature data on primary and secondary earthquake effects and add our own observations from fieldwork. We show that the ESI2007 distribution differs somewhat from traditional intensity assessments (MSK (Medvedev-Sponheuer-Karnik) and MM (Modified Mercalli)), because of the sparse population in the epicentral area and the spatial distribution of primary and secondary EEEs. However, the ESI2007 scale captures a similar overall pattern of the intensity distribution. We then explore how uncertainties in the identification of primary surface ruptures influence the results of the ESI2007 assignment. Our results highlight the applicability of the ESI2007 scale, even in earthquakes with complex and unusual primary surface rupture patterns.

Global Earthquake Response with Imaging Geodesy: Recent Examples from the USGS NEIC

The U.S. Geological Survey National Earthquake Information Center leads real-time efforts to provide rapid and accurate assessments of the impacts of global earthquakes, including estimates of ground shaking, ground failure, and the resulting human impacts. These efforts primarily rely on analysis of the seismic wavefield to characterize the source of the earthquake, which in turn informs a suite of disaster response products such as ShakeMap and PAGER. In recent years, the proliferation of rapidly acquired and openly available in-situ and remotely sensed geodetic observations has opened new avenues for responding to earthquakes around the world in the days following significant events. Geodetic observations, particularly from interferometric synthetic aperture radar (InSAR) and satellite optical imagery, provide a means to robustly constrain the dimensions and spatial complexity of earthquakes beyond what is typically possible with seismic observations alone. Here, we document recent cases where geodetic observations contributed important information to earthquake response efforts—from informing and validating seismically-derived source models to independently constraining earthquake impact products—and the conditions under which geodetic observations improve earthquake response products. We use examples from the 2013 Mw7.7 Baluchistan, Pakistan, 2014 Mw6.0 Napa, California, 2015 Mw7.8 Gorkha, Nepal, and 2018 Mw7.5 Palu, Indonesia earthquakes to highlight the varying ways geodetic observations have contributed to earthquake response efforts at the NEIC. We additionally provide a synopsis of the workflows implemented for geodetic earthquake response. As remote sensing geodetic observations become increasingly available and the frequency of satellite acquisitions continues to increase, operational earthquake geodetic imaging stands to make critical contributions to natural disaster response efforts around the world.

A New Method for Large-Scale Landslide Classification from Satellite Radar

Following a large continental earthquake, information on the spatial distribution of triggered landslides is required as quickly as possible for use in emergency response coordination. Synthetic Aperture Radar (SAR) methods have the potential to overcome variability in weather conditions, which often causes delays of days or weeks when mapping landslides using optical satellite imagery. Here we test landslide classifiers based on SAR coherence, which is estimated from the similarity in phase change in time between small ensembles of pixels. We test two existing SAR-coherence-based landslide classifiers against an independent inventory of landslides triggered following the Mw 7.8 Gorkha, Nepal earthquake, and present and test a new method, which uses a classifier based on coherence calculated from ensembles of neighbouring pixels and coherence calculated from a more dispersed ensemble of ‘sibling’ pixels. Using Receiver Operating Characteristic analysis, we show that none of these three SAR-coherence-based landslide classification methods are suitable for mapping individual landslides on a pixel-by-pixel basis. However, they show potential in generating lower-resolution density maps, which are used by emergency responders following an earthquake to coordinate large-scale operations and identify priority areas. The new method we present outperforms existing methods when tested at these lower resolutions, suggesting that it may be able to provide useful and rapid information on landslide distributions following major continental earthquakes.

The rise, collapse, and compaction of Mt. Mantap from the 3 September 2017 North Korean nuclear test

Surveillance of clandestine nuclear tests relies on a global seismic network, but the potential of spaceborne monitoring has been underexploited. We used satellite radar imagery to determine the complete surface displacement field of up to 3.5 meters of divergent horizontal motion with 0.5 meters of subsidence associated with North Korea's largest underground nuclear test. Combining insight from geodetic and seismological remote sensing, we found that the aftermath of the initial explosive deformation involved subsidence associated with subsurface collapse and aseismic compaction of the damaged rocks of the test site. The explosive yield from the nuclear detonation with best-fitting source parameters for 450-meter depth was 191 kilotonnes of TNT equivalent. Our results demonstrate the capability of spaceborne remote sensing to help characterize large underground nuclear tests.

Dynamic strain determination using fibre-optic cables allows imaging of seismological and structural features

Natural hazard prediction and efficient crust exploration require dense seismic observations both in time and space. Seismological techniques provide ground-motion data, whose accuracy depends on sensor characteristics and spatial distribution. Here we demonstrate that dynamic strain determination is possible with conventional fibre-optic cables deployed for telecommunication. Extending recently distributed acoustic sensing (DAS) studies, we present high resolution spatially un-aliased broadband strain data. We recorded seismic signals from natural and man-made sources with 4-m spacing along a 15-km-long fibre-optic cable layout on Reykjanes Peninsula, SW-Iceland. We identify with unprecedented resolution structural features such as normal faults and volcanic dykes in the Reykjanes Oblique Rift, allowing us to infer new dynamic fault processes. Conventional seismometer recordings, acquired simultaneously, validate the spectral amplitude DAS response between 0.1 and 100 Hz bandwidth. We suggest that the networks of fibre-optic telecommunication lines worldwide could be used as seismometers opening a new window for Earth hazard assessment and exploration.

Constant strain accumulation rate between major earthquakes on the North Anatolian Fault

Earthquakes are caused by the release of tectonic strain accumulated between events. Recent advances in satellite geodesy mean we can now measure this interseismic strain accumulation with a high degree of accuracy. But it remains unclear how to interpret short-term geodetic observations, measured over decades, when estimating the seismic hazard of faults accumulating strain over centuries. Here, we show that strain accumulation rates calculated from geodetic measurements around a major transform fault are constant for its entire 250-year interseismic period, except in the ~10 years following an earthquake. The shear strain rate history requires a weak fault zone embedded within a strong lower crust with viscosity greater than ~10²⁰ Pa s. The results support the notion that short-term geodetic observations can directly contribute to long-term seismic hazard assessment and suggest that lower-crustal viscosities derived from postseismic studies are not representative of the lower crust at all spatial and temporal scales.