Search

The preface for the article collection "10 years after the 2011 Tohoku earthquake: a milestone of solid earth science" highlights the scientific discoveries made since the earthquake, thanks to near-field observations and existing knowledge about the northeastern Japan arc. The earthquake revealed dynamic frictional behaviors of the subduction interface and provided insights into subduction zone dynamics. The collection of research on the Tohoku earthquake is expected to benefit the understanding of subduction earthquakes worldwide. The preface also mentions various research papers included in the collection, covering topics such as coseismic stress drop, afterslip, seafloor observations, slow earthquake activity, tsunami deposits, paleoseismology, and numerical modeling of the earthquake cycle. Ongoing research on postseismic deformation and historical research on earthquake occurrence are also discussed. The authors express their appreciation to the reviewers and the editorial board for their contributions to the special issue. [Extracted from the article]

Ocean bottom pressure-gauge (OBP) records play an essential role in seafloor geodesy. Oceanographic fluctuations in OBP data, however, pose as a significant noise source in seafloor transient crustal deformation observations, including slow slip events (SSEs), making it crucial to evaluate them quantitatively. To extract the significant fluctuation phenomena common to multiple observation networks, including oceanographic fluctuations and tectonic signals, we applied principal component analysis (PCA) to the 3-year Dense Oceanfloor Network System for Earthquakes and Tsunamis (DONET) OBP time series for 40 stations during 2016–2019. PCA could separate several oceanographic signals based on the characteristics of their spatial distributions, although evident transient tectonic signals could not be confirmed from the observed pressure records during this observed period. The spatial distribution of the first four principal components (PCs) reflected the common component, inclined component along sea depth, longitude component, and parabola-like pattern, respectively. By subtracting each PC (in particular, PC-2 and PC-4) from the time series, we could significantly reduce the sea depth dependence of OBP records, which has been highlighted in several previous studies and is also evident in this region. We interpreted PCs 2–4 as the reflection of the strength and meandering of ocean geostrophic currents based on a comparison with the PC spatial distribution of the numerical oceanographic models. In addition, to evaluate the ability of PCA to separate transient tectonic signal from OBP time series, including oceanographic fluctuations, we conducted a synthetic ramp assuming an SSE by rectangular fault and then applied PCA. The assumed synthetic tectonic signal could be separated from the oceanographic signals and included in the principal component independently depending on its amplitude, suggesting that the spatial distribution of each PC would change if the amplitude of the synthetic signal were sufficiently large. We propose a transient event-detection method based on the spatial distribution difference of a specific PC with or without a tectonic signal. We used the normalized inner product (NIP) between these PCs as the indicator of their similarities. This method can detect transient tectonic signals more significantly than the moment-magnitude scale of 5.9 from OBP records. [ABSTRACT FROM AUTHOR]

The 2011 Tohoku-Oki great earthquake increased the difficulty of evaluating the long-term probability of seismic activity along the Japan Trench because of the unknown impact of the unprecedentedly large slip. In this study, the long-term activity of "Miyagi-ken-Oki earthquakes", an M > 7 earthquake sequence off Miyagi Prefecture, located at the edge of the source area of the Tohoku-Oki earthquake was simulated. We conducted numerical simulations of earthquake generation cycles based on the rate- and state-dependent friction law representing the stress accumulation and release process on the plate interface. We also validated the results based on analyses of observational data, including time intervals and afterslip distributions for the previous Miyagi-ken-Oki earthquakes. The simulation results were then compared with repeating-earthquake-derived interplate slip observations over 30 years. The results showed that the spatial and temporal characteristics of aseismic slips on the plate interface near the M > 7 Miyagi-ken-Oki source were qualitatively reproduced. The time interval between the M ~ 9 earthquake and the first M > 7 earthquake is shorter than the average recurrence interval of M > 7 earthquakes during the latter stage of each M ~ 9 earthquake cycle. In contrast, the interval between the first and the second M > 7 earthquakes is the longest in each M ~ 9 earthquake cycle. The time intervals of the M > 7 earthquakes fluctuated in the early stage compared to those in the latter stage of the M ~ 9 earthquake cycle. These characteristics were associated with the extent of the locked and afterslip areas in and around the source. Hence, monitoring the spatio-temporal distribution of the slip rate in and around the source area during the preparation process of earthquakes occurring in the downdip area off Miyagi Prefecture is very important to assess whether the next M > 7 earthquake is approaching. Furthermore, earthquake cycle simulations combined with earthquake and slow slip monitoring can provide meaningful insights for long-term assessments of great interplate earthquakes. [ABSTRACT FROM AUTHOR]

We investigated temperature records associated with seafloor pressure observations at eight stations that experienced the 2011 Mw 9 Tohoku earthquake near its epicenter. The temperature data were based on the temperature measured inside the pressure transducer. We proposed a method to estimate ambient water temperature from the internal temperature using an equation of heat conduction. The estimated seafloor water temperature showed remarkable anomalies, especially increases several hours after the Mw 9 earthquake. A station of P03 (sea depth of 1.1 km) showed an abrupt temperature increase of + 0.19 °C that occurred ~ 3 h after the earthquake, which lasted for several hours. At stations of GJT3 (sea depth of 3.3 km) and TJT1 (sea depth of 5.8 km), there were abrupt temperature anomalies of + 0.20 °C and + 0.10 °C that began to occur 3–4 h after the earthquake. These anomalies both decayed to their original levels over a few tens of days. During the decay processes, only TJT1 showed several intermittent temperature rises. A water temperature anomaly within + 0.03 °C was found up to ~ 500 m above TJT1 2 weeks after the earthquake. There was no significant anomaly at the other five stations. Processes to cause these seafloor temperature anomalies were discussed. The temperature anomaly of P03 was reasonably caused by a tsunami-generated turbidity current, as also suggested by a previous study. Meanwhile, we proposed a scenario that the abrupt temperature anomalies of GJT3/TJT1 and the intermittent anomalies of TJT1 were caused by warm water discharges from the subseafloor. The pathways of the warm water were probably composed of the branch normal fault between GJT3 and TJT1, the reverse fault near TJT1, the backstop interface, and perhaps reverse faults at the frontal prism. The proposed scenario was almost compatible with other studies based on epicentral observations. We estimated the heat properties of the initial temperature anomalies of GJT3/TJT1. The estimated heat source might be explained by that most of the geothermal fluids trapped in those fault pathways were discharged to the seafloor immediately after the earthquake. The onsets of the subsequent intermittent anomalies of TJT1 were possibly activated by low or falling ocean tidal loading. [ABSTRACT FROM AUTHOR]

Long-term continuous observation of seafloor pressure is effective for detecting seafloor vertical deformations that are associated with transient tectonic phenomena such as slow slip events. Since the aseismic slip event prior to the 2011 Tohoku earthquake, several discoveries have been made on spontaneous slow slip events and various other types of slow earthquake along the Japan and Kuril Trenches. Seafloor observation network for earthquakes and tsunamis along the Japan Trench (S-net) is expected to provide invaluable information on slow slip activities via geodetic signals that are detected by pressure observation. This study inspects the quality of the S-net pressure data in view of seafloor geodesy by comparison with records obtained by more than 100 autonomous ocean bottom pressure recorders (OBPRs) deployed along the Japan Trench. OBPRs have long been standard tools in seafloor geodesy, and the data collected are considered a benchmark in terms of quality. Most of the S-net stations showed noise levels that are considerably higher than those of the OBPRs over periods of more than 2 d. We speculate that a strong correlation between pressure and temperature accounts for much of the long-term noise. In this study, the temperature-dependent fluctuation component was estimated by prediction filtering and removed from the original data, leading to a significant reduction in the noise level at 51 stations, which reached levels almost equivalent to those of OBPRs. Although no significant pressure changes have been identified as associated with the 2018 Boso SSE or repeated tremor bursts in the northern Japan Trench thus far, our findings indicate that these stations are sufficiently sensitive to detect slow slip events occurring nearby. [ABSTRACT FROM AUTHOR]

Temporal seismic observations from pop-up type ocean-bottom seismometers were used to detect tectonic tremors immediately following the 2011 Tohoku-Oki earthquake in the northern periphery of the aftershock area. Near-field observations clearly distinguished tremors from regular earthquakes based on their spectral shape in the frequency band of 1–4 Hz. In addition to tremors accompanied by very low-frequency earthquakes (VLFEs), we detected 130 tremors without known VLFE activity during April–October 2011. The newly detected tremors were in the vicinity of a sequence of small repeating earthquakes, indicating a mixed distribution of tremors and regular interplate earthquakes in the region. Tremor activity was high immediately after the deployment of seismometers and gradually decreased. In addition, the tremor activity fluctuated with two activations with an interval of approximately 90 days, similar to the intervals between tremor bursts after 2016. The results of the study suggest that the observed tremors occurred under the influence of aseismic slip caused by the decaying afterslip of the preceding Tohoku-Oki and Mw 7.4 interplate earthquakes and episodic accelerations with a quasi-periodicity unique to the area. [ABSTRACT FROM AUTHOR]

The 2011 Tohoku-Oki earthquake generated a surprisingly large near-trench slip, and earth scientists have devoted significant attention to understanding why. Some studies proposed special rupture mechanisms, such as extensive dynamic frictional weakening; others simulated this near-trench slip behavior without supposing the extensive dynamic weakening. However, we have not reached a decisive conclusion for this question due to limited spatial near-trench slip resolution. Hence, in this study we use new tsunami data recorded just above the large slip area in addition to offshore and onshore geodetic data to improve the spatial resolution of stress release in the Tohoku-Oki earthquake and quantitatively examine the mechanical state of the plate interface. A maximum slip of 53 m reaching the trench and an insignificant stress drop (< 3 MPa) at the shallowest portion of the fault were estimated. Based on our modeling results and the past experimental studies, it is suggested that friction at the shallow near-trench portion should be inherently low both before and during the earthquake. This result provides perspectives on the shallow slip behavior along the plate boundary, in which the strain energy accumulation at the deep portion of the fault accounts for the anomalous large shallow slip, but shallow mechanical coupling does not. A large shallow slip has been considered as a result of the release of sufficiently large strain energy at the shallow portion of the plate interface, but we suggest that shallow slips similar to that during the 2011 Tohoku-Oki earthquake may occur in any subduction zones where the energy sufficiently accumulates only in the deeper portion. [ABSTRACT FROM AUTHOR]

This study proposes a new method for the uncertainty estimation of coseismic slip distribution on the plate interface deduced from real-time global navigation satellite system (GNSS) data and explores its application for tsunami inundation prediction. Jointly developed by the Geospatial Information Authority of Japan and Tohoku University, REGARD (REal-time GEONET Analysis system for Rapid Deformation monitoring) estimates coseismic fault models (a single rectangular fault model and slip distribution model) in real time to support tsunami prediction. The estimated results are adopted as part of the Disaster Information System, which is used by the Cabinet Office of the Government of Japan to assess tsunami inundation and damage. However, the REGARD system currently struggles to estimate the quantitative uncertainty of the estimated result, although the obtained result should contain both observation and modeling errors caused by the model settings. Understanding such quantitative uncertainties based on the input data is essential for utilizing this resource for disaster response. We developed an algorithm that estimates the coseismic slip distribution and its uncertainties using Markov chain Monte Carlo methods. We focused on the Nankai Trough of southwest Japan, where megathrust earthquakes have repeatedly occurred, and used simulation data to assume a Hoei-type earthquake. We divided the 2951 rectangular subfaults on the plate interface and designed a multistage sampling flow with stepwise perturbation groups. As a result, we successfully estimated the slip distribution and its uncertainty at the 95% confidence interval of the posterior probability density function. Furthermore, we developed a new visualization procedure that shows the risk of tsunami inundation and the probability on a map. Under the algorithm, we regarded the Markov chain Monte Carlo samples as individual fault models and clustered them using the k-means approach to obtain different tsunami source scenarios. We then calculated the parallel tsunami inundations and integrated the results on the map. This map, which expresses the uncertainties of tsunami inundation caused by uncertainties in the coseismic fault estimation, offers quantitative and real time insights into possible worst-case scenarios. [ABSTRACT FROM AUTHOR]

An earthquake with a magnitude of 6.7 occurred in the Japan Sea off Yamagata on June 18, 2019. The mainshock had a source mechanism of reverse-fault type with a compression axis of WNW–ESE direction. Since the source area is positioned in a marine area, seafloor seismic observation is indispensable for obtaining the precise distribution of the aftershocks. The source area has a water depth of less than 100 m, and fishing activity is high. It is difficult to perform aftershock observation using ordinary free-fall pop-up type ocean bottom seismometers (OBSs). We developed a simple anchored-buoy type OBS for shallow water depths and performed the seafloor observation using this. The seafloor seismic unit had three-component seismometers and a hydrophone. Two orthogonal tiltmeters and an azimuth meter monitored the attitude of the package. For seismic observation at shallow water depth, we concluded that an anchored-buoy system would have the advantage of avoiding accidents. Our anchored-buoy OBS was based on a system used in fisheries. We deployed three anchored-buoy OBSs in the source region where the water depth was approximately 80 m on July 5, 2019, and two of the OBSs were recovered on July 13, 2019. Temporary land seismic stations with a three-component seismometer were also installed. The arrival times of P- and S-waves were read from the records of the OBSs and land stations, and we located hypocenters with correction for travel time. A preliminary location was performed using absolute travel time and final hypocenters were obtained using the double-difference method. The aftershocks were distributed at a depth range of 2.5 km to 10 km and along a plane dipping to the southeast. The plane formed by the aftershocks is consistent with the focal mechanism of the mainshock. The activity region of the aftershocks was positioned in the upper part of the upper crust. Focal mechanisms were estimated using the polarity of the first arrivals. Although many aftershocks had a reverse-fault focal mechanism similar to the focal solution of the mainshock, normal-fault type and strike–slip fault type focal mechanisms were also estimated. [ABSTRACT FROM AUTHOR]

Rapid estimation of the coseismic fault model for medium-to-large-sized earthquakes is key for disaster response. To estimate the coseismic fault model for large earthquakes, the Geospatial Information Authority of Japan and Tohoku University have jointly developed a real-time GEONET analysis system for rapid deformation monitoring (REGARD). REGARD can estimate the single rectangular fault model and slip distribution along the assumed plate interface. The single rectangular fault model is useful as a first-order approximation of a medium-to-large earthquake. However, in its estimation, it is difficult to obtain accurate results for model parameters due to the strong effect of initial values. To solve this problem, this study proposes a new method to estimate the coseismic fault model and model uncertainties in real time based on the Bayesian inversion approach using the Markov Chain Monte Carlo (MCMC) method. The MCMC approach is computationally expensive and hyperparameters should be defined in advance via trial and error. The sampling efficiency was improved using a parallel tempering method, and an automatic definition method for hyperparameters was developed for real-time use. The calculation time was within 30 s for 1 × 106 samples using a typical single LINUX server, which can implement real-time analysis, similar to REGARD. The reliability of the developed method was evaluated using data from recent earthquakes (2016 Kumamoto and 2019 Yamagata-Oki earthquakes). Simulations of the earthquakes in the Sea of Japan were also conducted exhaustively. The results showed an advantage over the maximum likelihood approach with a priori information, which has initial value dependence in nonlinear problems. In terms of application to data with a small signal-to-noise ratio, the results suggest the possibility of using several conjugate fault models. There is a tradeoff between the fault area and slip amount, especially for offshore earthquakes, which means that quantification of the uncertainty enables us to evaluate the reliability of the fault model estimation results in real time. [ABSTRACT FROM AUTHOR]

There are significant differences between the middle and southern segments of the Japan Trench in terms of the seismic and aseismic slips on the plate interface and seismic velocity structures. Although the large coseismic slip of the 2011 Tohoku-Oki earthquake was limited to the middle segment, the observed negative residual gravity anomaly area in the southern segment corresponds to the postseismic slip area of the Tohoku-Oki earthquake. A density distribution model can explain the different slip behaviours of the two segments by considering their structural differences. The model indicates that the plate interface in the south was covered with a thick channel layer, as indicated by seismic survey imaging, and this layer resulted in a residual gravity anomaly. Numerical simulations which assumed evident frictional heterogeneity caused by the layer in the south efficiently reproduced M9 earthquakes recurring only in the middle, followed by evident postseismic slips in the south. This study proposes that although the layer makes the megathrust less compliant to seismic slip, it promotes aseismic slips following the growth of seismic slips on the fault in an adjacent region. [ABSTRACT FROM AUTHOR]

Obvious crustal deformation is observed during a postseismic period as well as a coseismic period associated with a large earthquake. Major mechanisms of transient postseismic deformation are known as afterslip and viscoelastic relaxation. Since the viscoelastic relaxation occurs as a response to a coseismic slip, postseismic deformation provides information on coseismic deformation through the viscoelastic response. However, most previous studies have not thoroughly utilized postseismic geodetic observational data for revealing coseismic slip behaviors. In this study, we developed a slip inversion method that simultaneously estimates coseismic slip and postseismic slip distributions from coseismic and postseismic geodetic observational data using viscoelastic Green's function (viscoelastic inversion method). We investigated the performance of the viscoelastic inversion method via two synthetic tests: one assumed a strike–slip event along an inland fault, while the other assumed a dip–slip event along a plate interface in a subduction zone. Both synthetic tests demonstrated that when extensive postseismic observational data were given, the viscoelastic inversion method provided a superior spatial resolution of coseismic slip distributions compared to conventional elastic inversion distributions. We also applied the viscoelastic inversion method to co- and post-seismic deformations associated with the 2011 Tohoku-oki earthquake. The seafloor geodetic observational network of the off-Tohoku region has been widely extended after the occurrence of the mainshock. Using this extended seafloor geodetic observational data, we successfully improved the spatial resolution of the coseismic slip distribution through the viscoelastic inversion method. Furthermore, using the seafloor observational data during the postseismic period, our inversion method enables us to obtain high spatial resolution of the coseismic slip in the offshore area and a reasonable coseismic slip distribution even if seafloor observational data during the coseismic period are unavailable. These results clarify the importance of deploying a geodetic observational network even after large coseismic events to assess past coseismic slip behaviors by considering the viscoelasticity of the Earth. [ABSTRACT FROM AUTHOR]

On December 7, 2012, an earthquake occurred within the Pacific Plate near the Japan Trench, which was composed of deep reverse- and shallow normal-faulting subevents (Mw 7.2 and 7.1, respectively) with a time interval of ~10 s. It had been known that the stress state within the plate was characterized by shallow tensile and deep horizontal compressional stresses due to the bending of the plate (bending stress). This study estimates the fault model of the doublet earthquake utilizing tsunami, teleseismic, and aftershock data and discusses the stress state within the incoming plate and spatiotemporal changes seen in it after the 2011 Tohoku-Oki earthquake. We obtained the vertical extents of the fault planes of deep and shallow subevents as ~45–70 km and ~5 (the seafloor)–35 km, respectively. The down-dip edge of the shallow normal-faulting seismic zone (~30–35 km) deepened significantly compared to what it was in 2007 (~25 km). However, a quantitative comparison of the brittle strength and bending stress suggested that the change in stress after the Tohoku-Oki earthquake was too small to deepen the down-dip end of the seismicity by ~10 km. To explain the seismicity that occurred at a depth of ~30–35 km, the frictional coefficient in the normal-faulting depth range required would have had to be ~0.07 ≤ μ ≤ ~0.2, which is significantly smaller than the typical friction coefficient. This suggests the infiltration of pore fluid along the bending faults, down to ~30–35 km. It is considered that the plate had already yielded to a depth of ~35 km before 2011 and that the seismicity of the area was reactivated by the increase in stress from the Tohoku-Oki earthquake. [ABSTRACT FROM AUTHOR]

The tsunami disasters that occurred in Indonesia, Chile, and Japan have inflicted serious casualties and damaged social infrastructures. Tsunami forecasting systems are thus urgently required worldwide. We have developed a real-time tsunami inundation forecast system that can complete a tsunami inundation and damage forecast for coastal cities at the level of 10-m grid size in less than 20 min. As the tsunami inundation and damage simulation is a vectorizable memory-intensive program, we incorporate NEC’s vector supercomputer SX-ACE. In this paper, we present an overview of our system. In addition, we describe an implementation of the program on SX-ACE and evaluate its performance of SX-ACE in comparison with the cases using an Intel Xeon-based system and the K computer. Then, we clarify that the fulfillment of a real-time tsunami inundation forecast system requires a system with high-performance cores connected to the memory subsystem at a high memory bandwidth such as SX-ACE. [ABSTRACT FROM AUTHOR]

We determined the three-dimensional structure of the crust and upper mantle off the Boso Peninsula, Japan, by analyzing seismograms recorded by ocean-bottom seismometers and land stations between 2011 and 2013. We employed seismic tomography to determine the P- and S-wave velocity structures and earthquake locations simultaneously. The tomographic image shows that the mantle parts of the Pacific and the Philippine Sea plates have high-velocity anomalies. The upper boundary of the Philippine Sea plate is delineated as approximately 2-6 km shallower than that previously estimated from land-based data for the area 140.5°E-141.5°E and 35°N-35.5°N. A pronounced low-velocity anomaly in P- and S-waves with low- V / V ratio (1.5-1.6) was observed at depths shallower than 20 km in the overriding North American plate. This anomaly may be caused by the presence of rocks with a low- V / V ratio, such as quartzite, and the water expelled from the subducted Pacific and Philippine Sea plates. [ABSTRACT FROM AUTHOR]

We determined hypocenter distributions and focal mechanisms of earthquakes off the Boso Peninsula, Japan, based on seafloor observations using ocean-bottom seismometers (OBSs). From OBS data, we detected and located earthquakes that were closer to the trenches than had previously been reported. The determined focal depths of relocated hypocenters were systematically shallower than those in the Japan Meteorological Agency catalog determined from land-based data. The hypocenter distribution was separated into two groups: one associated with the Pacific plate (PAC) and the other with the North American and/or Philippine Sea plates (PHS). For the former group, both the focal depths and mechanisms of low-angle thrust-faulting earthquakes were consistent with the geometry of the PAC determined by previous studies. Among the latter group, we selected low-angle thrust-faulting earthquakes with a slip direction parallel to the direction of PHS subduction or with a dip direction parallel to the PHS to delineate the upper boundary of the PHS. The depth of the plate upper boundary off the Boso Peninsula was found to be approximately 6 km shallower than previously reported estimates. [ABSTRACT FROM AUTHOR]

Background: A GPS/acoustic (GPS/A) geodetic observation technique allows us to determine far offshore plate motion in order to understand the mechanism of megathrust earthquakes. In this technique, the distance between a sea-surface platform and seafloor transponders is estimated using the two-way traveltimes (TWT) of acoustic signals. TWTs are determined by maximizing the cross-correlation coefficient between the transmitted and returned signals. However, this analysis caused significantly wrong detection of TWT when the correlogram has an enlarged secondary envelope due to the enlarged amplitude of multiple signals depending on the relative spatial geometry between the ship and the transponder. The handled manual rereading of thousands of correlograms to obtain correct TWTs needs enormous time, and human errors may cause. To prevent these difficulties, an automated TWT determination procedure is valid to process numerous GPS/A data efficiently not only without human errors but also with high precision. Proposed methods: We developed automated methods for precisely analyzing GPS/A data. Method 1: The maximum peak in the observed correlogram is read, and a synthetic correlogram is then subtracted from the observation. Then, the same operation is applied to the subtracted waveform. This procedure is iterated until the correlation coefficient lowers than a pre-defined threshold. A true traveltime is defined as the fastest traveltime during the iterations. Method 2: The observed correlograms are divided into several groups based on their similarity through cluster analysis, and a master waveform in each group is selected. Then, the traveltime residual between the maximum and true peaks in the master waveform is manually evaluated. The obtained residual is employed as the correction value for each slave waveform. Further, we employed a seismic data projection to visually inspect the reliability of obtained results. Results: We confirmed that both new methods accurately correct misreadings in the current method, which amount to 0.4-0.8 ms roughly corresponding to 30-60 cm difference in the slant range. Conclusions: Thus, the proposed algorithms significantly improve the estimation of the transponder location. Further analyses are required to determine the arbitrary threshold values and to construct fully automated algorithms. [ABSTRACT FROM AUTHOR]

After a large subduction earthquake, crustal deformation continues to occur, with a complex pattern of evolution. This postseismic deformation is due primarily to viscoelastic relaxation of stresses induced by the earthquake rupture and continuing slip (afterslip) or relocking of different parts of the fault. When postseismic geodetic observations are used to study Earth's rheology and fault behaviour, it is commonly assumed that short-term (a few years) deformation near the rupture zone is caused mainly by afterslip, and that viscoelasticity is important only for longer-term deformation. However, it is difficult to test the validity of this assumption against conventional geodetic data. Here we show that new seafloor GPS (Global Positioning System) observations immediately after the great Tohoku-oki earthquake provide unambiguous evidence for the dominant role of viscoelastic relaxation in short-term postseismic deformation. These data reveal fast landward motion of the trench area, opposing the seaward motion of GPS sites on land. Using numerical models of transient viscoelastic mantle rheology, we demonstrate that the landward motion is a consequence of relaxation of stresses induced by the asymmetric rupture of the thrust earthquake, a process previously unknown because of the lack of near-field observations. Our findings indicate that previous models assuming an elastic Earth will have substantially overestimated afterslip downdip of the rupture zone, and underestimated afterslip updip of the rupture zone; our knowledge of fault friction based on these estimates therefore needs to be revised. [ABSTRACT FROM AUTHOR]

In order to understand the characteristics of shallow very low-frequency (VLF) events as revealed by recent ocean-floor observation studies, we perform a trial simulation of earthquake cycles in the Tonankai district by taking the characteristics of the 1944 Tonankai earthquake and assuming that slow earthquakes occur on numerous small asperities. Our simulation results show that the increase of moment release rate of shallower VLF events in the pre-seismic stage of a megathrust earthquake is higher than that of deeper VLF events. This increase may make leveling change due to VLF swarms detectable at Dense Oceanfloor Network system for Earthquakes and Tsunamis (DONET). We also introduce the time series of hydraulic pressure data at DONET, comparing with the leveling change expected from our numerical simulation. Since leveling change due to shallower VLF swarms is so local as to be incoherent, removal of the moving-averaged data from the data stacked by four nearby observation points in the same node may be useful to detect the short-term local leveling change. [ABSTRACT FROM AUTHOR]

A global barotropic ocean model forced by atmospheric disturbances is developed for the detection of seafloor vertical displacements from in situ ocean bottom pressure (OBP) data. The model accuracy is validated by deep-sea OBP data at more than 100 sites obtained over the global ocean. Parameters and boundary conditions including the horizontal resolution incorporated in the ocean model are tested in order to accurately simulate the nontidal (>2 days) OBP variations. The horizontal resolution is found to the factor that most significantly affects the simulated result. The finer the horizontal resolution applied, the smaller the model variability is. The model accuracy is highest when the horizontal resolution is 1/12°, but deteriorates when the horizontal resolution is finer than 1/12°. This may indicate a failure of the energy dissipation parameterization in the barotropic ocean model. Using the developed 1/12° model, the root-mean-square of the observed nontidal OBP component can be reduced by 18 % as an average of all the OBP data used. It is found that the 1/12° model is useful for the detection of a slow seafloor vertical displacement of centimeters related to the 2011 Tohoku-Oki earthquake from in situ OBP records near the hypocenter of the earthquake. [ABSTRACT FROM AUTHOR]

Shear-wave anisotropy in Earth's mantle helps constrain the lattice-preferred orientation of anisotropic minerals due to viscous flow. Previous studies at the Japan Trench subduction zone using land-based seismic networks identified strong anisotropy in the mantle wedge, reflecting viscous flow induced by the subducting slab. Here we map anisotropy in the previously uninvestigated offshore region by analyzing shear waves from interplate earthquakes that are recorded by a new seafloor network (the S-net). The newly detected anisotropy is not in the mantle wedge but only in the overlying crust (∼0.1 s time delay and trench-parallel fast direction). The distinct lack of anisotropy indicates that the forearc mantle wedge offshore is decoupled from the slab and does not participate in the viscous flow, in sharp contrast with the rest of the mantle wedge. A stagnant forearc mantle wedge provides a stable and cold tectonic environment that is important for the petrological evolution and earthquake processes of subduction zones. Knowledge of shear-wave anisotropy is important to understanding the structure and dynamics of the subduction zone mantle wedge. Here, the authors find unambiguous evidence that forearc anisotropy resides in the upper-plate crust, while weak anisotropy in the most seaward part of the mantle wedge indicates decoupling from the slab [ABSTRACT FROM AUTHOR]

We applied tomographic inversion and receiver function analysis to seismic data from ocean-bottom seismometers and land-based stations to understand the structure and its relationship with slow slip events off Boso, Japan. First, we delineated the upper boundary of the Philippine Sea Plate based on both the velocity structure and the locations of the low-angle thrust-faulting earthquakes. The upper boundary of the Philippine Sea Plate is distorted upward by a few kilometers between 140.5 and 141.0°E. We also determined the eastern edge of the Philippine Sea Plate based on the delineated upper boundary and the results of the receiver function analysis. The eastern edge has a northwest–southeast trend between the triple junction and 141.6°E, which changes to a north–south trend north of 34.7°N. The change in the subduction direction at 1–3 Ma might have resulted in the inflection of the eastern edge of the subducted Philippine Sea Plate. Second, we compared the subduction zone structure and hypocenter locations and the area of the Boso slow slip events. Most of the low-angle thrust-faulting earthquakes identified in this study occurred outside the areas of recurrent Boso slow slip events, which indicates that the slow slip area and regular low-angle thrust earthquakes are spatially separated in the offshore area. In addition, the slow slip areas are located only at the contact zone between the crustal parts of the North American Plate and the subducting Philippine Sea Plate. The localization of the slow slip events in the crust–crust contact zone off Boso is examined for the first time in this study. Finally, we detected a relatively low-velocity region in the mantle of the Philippine Sea Plate. The low-velocity mantle can be interpreted as serpentinized peridotite, which is also found in the Philippine Sea Plate prior to subduction. The serpentinized peridotite zone remains after the subduction of the Philippine Sea Plate and is likely distributed over a wide area along the subducted slab. [ABSTRACT FROM AUTHOR]

We assessed the positioning accuracy of GNSS-Acoustic (GNSS-A) measurement using a slackly moored buoy. A key feature of real-time buoy-based GNSS-A measurement is that positioning must be performed via single ranging from an arbitrary observation position and encompassing the drifting range of the buoy. In this study, slack-line mooring was employed to resist strong Kuroshio current of up to 5.5 knots in the Nankai Trough. During a year-long sea trial, the buoy drifted within a circle of ~ 4000 m radius, which is much larger than the dimension of a seafloor transponder array. In the sea trial, more than 500 successful regular and on-demand ranging measurements were obtained at various observation positions. The results show that the horizontal positioning accuracy (2σ) of the array for single ranging was 46 cm when the buoy was inside the array and 97 cm when it was outside the array. These accuracies are comparable to those achieved through single ranging in traditional ship-based surveys at the same site. In addition to the distance dependency, we observed directivity in the accuracy, depending on the geometry between an observation point and seafloor transponders. To interpret the accuracy degradation and directivity as a function of observation position, we calculated error ellipses using a dilution-of-precision (DOP) analysis procedure for GNSS-A positioning. The error ellipses clearly illustrate that the variance in array positions is large in the line-of-sight direction from an observation point to the array center. We translated the error ellipses to positioning accuracies using the standard deviation of travel time residuals when an array position is estimated using all pings. The positioning accuracy resulting from the translation corresponded to that obtained using array positions in a buoy observation. The DOP analysis results and their translation into positioning accuracies enabled a detailed assessment of the directional accuracy of GNSS-A positioning at an arbitrary observation position, which is important for real-time measurement of seafloor movement in the event of a huge earthquake. [ABSTRACT FROM AUTHOR]