Nonlinear multidimensional scaling and visualization of earthquake clusters over space, time and feature space

International audienceWe present a novel technique based on a multi-resolutional clustering and nonlinear multi-dimensional scaling of earthquake patterns to investigate observed and synthetic seismic catalogs. The observed data represent seismic activities around the Japanese islands during 1997-2003. The synthetic data were generated by numerical simulations for various cases of a heterogeneous fault governed by 3-D elastic dislocation and power-law creep. At the highest resolution, we analyze the local cluster structures in the data space of seismic events for the two types of catalogs by using an agglomerative clustering algorithm. We demonstrate that small magnitude events produce local spatio-temporal patches delineating neighboring large events. Seismic events, quantized in space and time, generate the multi-dimensional feature space characterized by the earthquake parameters. Using a non-hierarchical clustering algorithm and nonlinear multi-dimensional scaling, we explore the multitudinous earthquakes by real-time 3-D visualization and inspection of the multivariate clusters. At the spatial resolutions characteristic of the earthquake parameters, all of the ongoing seismicity both before and after the largest events accumulates to a global structure consisting of a few separate clusters in the feature space. We show that by combining the results of clustering in both low and high resolution spaces, we can recognize precursory events more precisely and unravel vital information that cannot be discerned at a single resolution

Spacial distribution and time variation in seismicity around the Antarctic Plate-Indian Ocean region

Spacial distribution and time variations in seismicity around the Antarctic Plate-Indian Ocean (0-160°E , 20-80°S ) is evaluated based on the data compiled at the global seismological centers since 1964. Seismicity in the oceanic area of the Antarctic plate, along with ridges/transform faults between the Indian-Australian Plate, represents characteristic features before and after large earthquakes, such as the Balleny Earthquake on March 25, 1998. Seismicity in the aseismic ridges immediately east of the Australia-Antarctic Discordance (AAD) increased the year before occurrence of the Balleny Earthquake. Long period variations in seismicity during more than 30 years imply several characteristic time periods of increase in cumulative seismicity, associated with dynamic distribution of tectonic stress in space and time between adjacent plates. Time variations in seismicity around the Balleny Earthquake region, in particular, drastically changed before and after the main shock in March 1998. The recent distribution of hypocenters in this area appears to extend toward Wilks Land, followed by the excitation of local events beneath the continental ice sheet of Antarctica. Seismic activities of these areas might reflect far-field tectonic stress in the lithosphere around the region

Global observation of vertical-CLVD earthquakes at active volcanoes

Some of the largest and most anomalous volcanic earthquakes have non-double-couple focal mechanisms. Here, we investigate the link between volcanic unrest and the occurrence of non-double-couple earthquakes with dominant vertical tension or pressure axes, known as vertical compensated-linear-vector-dipole (vertical-CLVD) earthquakes. We determine focal mechanisms for 313 target earthquakes from the standard and surface wave catalogs of the Global Centroid Moment Tensor Project and identify 86 shallow 4.3 ≤ MW ≤ 5.8 vertical-CLVD earthquakes located near volcanoes that have erupted in the last ~100 years. The majority of vertical-CLVD earthquakes occur in subduction zones in association with basaltic-to-andesitic stratovolcanoes or submarine volcanoes, although vertical-CLVD earthquakes are also located in continental rifts and in regions of hot spot volcanism. Vertical-CLVD earthquakes are associated with many types of confirmed or suspected eruptive activity at nearby volcanoes, including volcanic earthquake swarms as well as effusive and explosive eruptions and caldera collapse. Approximately 70% of all vertical-CLVD earthquakes studied occur during episodes of documented volcanic unrest at a nearby volcano. Given that volcanic unrest is underreported, most shallow vertical-CLVD earthquakes near active volcanoes are likely related to magma migration or eruption processes. Vertical-CLVD earthquakes with dominant vertical pressure axes generally occur after volcanic eruptions, whereas vertical-CLVD earthquakes with dominant vertical tension axes generally occur before the start of volcanic unrest. The occurrence of these events may be useful for identifying volcanoes that have recently erupted and those that are likely to erupt in the future

Investigating the role of the Itoigawa-Shizuoka tectonic line towards the evolution of the Northern Fossa Magna rift basin

AbstractThe Itoigawa-Shizuoka tectonic line (ISTL) fault system is considered to have one of the highest probabilities for a major inland earthquake occurrence in the whole of Japan. It is a complex fault system with the dip directions of the local fault segments changing from north to south between an east-dipping low-angle thrust fault, a strike slip fault and a west-dipping thrust fault. The tectonic relations between the different parts of the fault system and the surrounding geological units are yet to be fully explained. This study aims to reveal the juncture of the northern and central parts of the ISTL and investigate its contribution towards the shaping of the Northern Fossa Magna rift basin. We conducted 3 deployments of 1 or 2 linear arrays of seismic stations across the central and northern ISTL regions and observed local micro-earthquakes for a period of 3 years. Each deployment recorded continuous waveform data for approximately 3 months. Using arrival times of 1193 local earthquakes, we jointly determined earthquake locations and a 3D velocity model, applying the tomography method. We were able to image the regional crustal structures from the surface to a depth of 20km with a spatial resolution of 5km. Subsequently, we used the obtained 3D velocity model to relocate the background local seismicity from 2003 to 2009. The juncture of the northern and central parts of the ISTL was well constrained by our results. The depth extension of the northern parts of the ISTL fault segments follows the bottom of the Miocene Northern Fossa Magna rift basin (NFM) and forms an east-dipping low-angle fault. In contrast, the central parts of the ISTL fault segments are estimated to lie along the eastern boundary of the Matsumoto basin forming an oblique strike slip fault (Fig. 1)

Ambient seismic noise tomography of the southern East Sea (Japan Sea) and the Korea Strait

Group velocity maps were derived for the southern East Sea (Japan Sea) and the Korea Strait (Tsushima Strait) for the 5–36 s period range, which is sensitive to shear wave velocities of the crust and the uppermost mantle. Images produced in our study enhance our understanding of the tectonic evolution of a continental margin affected by subducting oceanic slabs and a colliding continental plate. The seismic structure of the study area has not been described well because seismic data for the region are scarce. In this study, we applied the ambient noise tomography technique that does not rely on earthquake data. We calculated ambient noise cross-correlations recorded at station pairs of dense seismic networks located in the regions surrounding the study area, such as the southern Korean Peninsula and southwestern part of the Japanese Islands. We then measured the group velocity dispersion curves of the fundamental mode Rayleigh waves from cross-correlograms and constructed 2-D group velocity maps reflecting group velocity structure from the upper crust to uppermost mantle. The results show that three distinct anomalies with different characteristics exist. Anomalies are located under the Ulleung Basin (UB), the boundary of the Basin, and the area between Tsushima Island and the UB. 1-D velocity models were obtained by inversion of dispersion curves that represent vertical variations of shear wave velocity at locations of three different anomalies. The 1-D velocity models and 2-D group velocity maps of lateral variations in shear wave group velocities show that the high velocity anomaly beneath the UB originates from crustal thinning and mantle uplift. Confirming the exact causes of two low velocity anomalies observed under the UB boundary and between Tsushima Island and the UB is difficult because additional information is unavailable. However, complex fault systems, small basins formed by faulting, and deep mantle flow can be possible causes of the existence of low velocity anomalies in the region.This work was funded by the Korean Meteorological Administration and Development Program under Grant CATER 2012-5051

Summary: Global Seismology and the Polar Region

“Polar Seismology” has been developed since the International Geophysical Year (IGY 1957–1958) and contributed significantly to global seismology in particular through the big project of the International Polar Year (IPY 2007–2008). At present, in the first stage of the twenty-first century, “polar regions” play an important role to monitor and understand the drastic variations in the Earth’s system as well as to advance the inter-disciplinary studies of the interactions among multispheres within the system