Search

Seismic fault surfaces have complex geometries over a broad-scale range. The Hurst exponent (H) is an important parameter characterizing the complexity of fault geometries. This exponent of seismic faults has only been estimated at outcrops and from surface traces of large earthquakes. This is because geometry data of subsurface faults usually have large uncertainties and many gaps. This study examined the applicability of the lifting scheme, which is one of the wavelet transform methods, to estimate H of irregularly sampled geometry of subsurface faults. We analysed the surface geometry of the subducting oceanic plate at the Nankai trough, Japan, which is part of the fault plane of interplate earthquakes. The geometries of the subducting plate were estimated along six survey lines by integrating seismic refraction and reflection surveys. Two-way traveltimes of reflected waves from the plate surface, which were measured from the reflection survey data, were converted to depth using the velocity structure estimated by refraction survey. The intervals between sampled points were irregular because the reflected waves were obscured or invisible in some places. The Hurst exponents were estimated from the scale dependence of the wavelet coefficients that were derived by the lifting scheme without interpolation. Analyses of the synthetic data simulating the irregularly sampled plate geometries indicated that the lifting scheme yields stable but largely biased estimates of H , especially for small H (<0.5). We introduced the empirical bias correction to achieve an unbiased estimation of the exponent. The analysis of plate geometry at the Nankai trough was conducted at narrow-scale ranges with consideration for the accuracy of velocity structures. We may conclude that H of five survey lines across the trough axis is > 0.8, and that of a line along trough is > 0.7. These estimated exponents had large uncertainties due to analyses at narrow-scale ranges, but were close to the estimates in studies of the surface traces of large earthquakes. [ABSTRACT FROM AUTHOR]

The process by which continental crust has formed is not well understood, though such crust mostly forms at convergent plate margins today. It is thus imperative to study modern intra-oceanic arcs, such as those common in the western Pacific Ocean. New seismic studies along the representative Izu-Bonin intra-oceanic arc provide unique along-strike images of arc crust and uppermost mantle to complement earlier, cross-arc lithospheric profiles. These reveal two scales (1000-10 km scale) of variations, one at the scale of the Izu versus Bonin (thick versus thin) arc crust, the other at the intervolcano (~50 km) scale. These images show that: (1) the bulk composition of the Izu-Bonin arc crust is more mafic than typical continental crust, (2) the middle crust with seismic velocities similar to continental crust is predominantly beneath basaltic arc volcanoes, (3) the bulk composition beneath basaltic volcanoes changes little at thick and thin arc segments, and (4) a process to return lower crustal components to the mantle, such as delamination, is required for an arc crust to evolve into continental crust. Continued thickening of the Izu-Bonin crust, accompanied by delamination of lowermost crust, can yield velocity structure of typical continental crust. Keywords: intra-oceanic arc, continental crust, crustal evolution, seismic imaging, Izu-Bonin arc. doi: 10.1130/G23901A.1

A new high-resolution velocity model of the Mariana arc-backarc system obtained from active-source seismic profiling demonstrates velocity variations within the arc middle and lower crusts of intermediate to felsic and mafic compositions. The characteristics of the oceanic-island-arc crust are a middle crust with velocity of ~6 km/s, laterally heterogeneous lower crust with velocities of ~7 km/s, and unusually low mantle velocities. Petrologic modeling suggests that the volume of the lower crust, composed of restites and olivine cumulates after the extraction of the middle crust, should be significantly larger than is observed, suggesting that a part of the lower crust, especially the cumulates, is seismically a part of the mantle. Keywords: seismic structure, intra-oceanic island arc, crustal growth, Izu-Bonin-Mariana arc.

The offshore part of the North Anatolian Fault (NAF) beneath the Marmara Sea is a well-known seismic gap for future M > 7 earthquakes in the sense that more than 250 yr have passed since the last major earthquake in the Central Marmara region. Although many studies discussed the seismic potential for the future large earthquake in this region on the basis of historical record, geodetic and geological observations, it is difficult to evaluate the actual situation on the seismic activity and structure along the NAF beneath the Marmara Sea due to the lack of ocean bottom seismic observations. Using ocean bottom seismometer observations, an assessment of the location of possible asperities that could host an expected large earthquake is undertaken based on heterogeneities in the microseismicity distribution and seismic velocity structure. Specifically, seismic tomography and precise hypocentre estimations are conducted using offshore seismic data whose recording period is 11 months. About five times more microearthquakes are detected with respect to events recorded in a land-based catalogue. A comparison with previously published results from offshore observation data suggests that the seismicity pattern had not changed from 2014 September to 2017 May. The location accuracy of microearthquakes is greatly improved from only the land-based earthquake catalogue, particularly for depth direction. There are several aseismic and inactive zones of microearthquake, and the largest one is detected using land-based seismic observation, whereas other zones are newly detected via offshore observations. The obtained velocity model shows a strong lateral contrast, with two changing points. The western changing point corresponds to a segmentation boundary, where the dip angle of the NAF segments changed. High-velocity zones from tomographic images are characterized by low seismicity eastward of the segment boundary. To the east of 28.50°E, the high-velocity zone becomes thicker in the depth direction and is characterized by low seismicity. Although the low seismic activity alone could be interpreted as both strong coupling and fully creeping, the high-velocity features at the same can be concluded that these zones are consist of brittle material and strong coupling. From comparison with other geodetic and seismic studies, we interpret these zones as locked zones that had been ruptured by the past large earthquakes and could be ruptured by future ones. These zones might accumulate strain since the main shock rupture associated with the 1766 May M s 7.3 earthquake, the latest major earthquake in this region. [ABSTRACT FROM AUTHOR]

To decrease human and economic damage owing to earthquakes, it is necessary to discover signals preceding earthquakes. We focus on the concept of "early warning signals" developed in bifurcation analysis, in which an increase in the variances of variables precedes its transition. If we can treat earthquakes as one of the transition phenomena that moves from one state to the other state, this concept is useful for detecting earthquakes before they start. We develop a covariance matrix from multi-channel time series data observed by an observatory on the seafloor and calculate the first eigenvalue and corresponding eigenstate of the matrix. By comparing the time dependence of the eigenstate to some past earthquakes, it is shown that the contribution from specific observational channels to the eigenstate increases before earthquakes, and there is a case in which the eigenvalue increases as predicted in early warning signals. This result suggests the first eigenvalue and eigenstate of multi-channel data are useful to identify signals preceding earthquakes. [ABSTRACT FROM AUTHOR]

Three-dimensional (3D) modeling of tsunami events is intended to promote tsunami safety. However, the developed 3D modeling methods based on Computational Fluid Dynamics and photorealistic particle visualization have some weaknesses, such as not being similar to the original environment, not measuring the wave's end point, and low image accuracy. The method for 3D modeling of tsunamis that results from this research can fulfil those weaknesses because it has advantages, such as being able to predict the end point of waves, similar to the original environment, and the height and area of inundation. In addition, the method produces more detailed and sharper spatial data. Modeling in this research is conducted using Agisoft Metashape Professional software to a produce 3D orthomosaic from pictures taken with Unmanned Aerial Vehicle (UAV) technique or drone (photogrammetry), and 3ds max software is used for wave simulation. We take a sample of an area in Cilacap, Indonesia that was impacted by the 2006 southwest coast tsunamis and may be vulnerable to future big megathrust earthquakes and tsunamis. The results could be used to provide several benefits, such as the creation of evacuation routes and the determination of appropriate locations for building shelters. [ABSTRACT FROM AUTHOR]