Coseismic and postseismic slip of the 2011 magnitude-9 Tohoku-Oki earthquake

Volcanic arc rigidity variations illuminated by coseismic deformation of the 2011 Tohoku-oki M9

Rock strength has long been linked to lithospheric deformation and seismicity. However, independent constraints on the related elastic heterogeneity are missing, yet could provide key information for solid Earth dynamics. Using coseismic Global Navigation Satellite Systems (GNSS) data for the 2011 M9 Tohoku-oki earthquake in Japan, we apply an inverse method to infer elastic structure and fault slip simultaneously. We find compliant material beneath the volcanic arc and in the mantle wedge within the partial melt generation zone inferred to lie above ~100 km slab depth. We also identify low-rigidity material closer to the trench matching seismicity patterns, likely associated with accretionary wedge structure. Along with traditional seismic and electromagnetic methods, our approach opens up avenues for multiphysics inversions. Those have the potential to advance earthquake and volcano science, and in particular once expanded to InSAR type constraints, may lead to a better understanding of transient lithospheric deformation across scales.

Complex tsunamigenic near-trench seafloor deformation during the 2011 Tohoku-Oki earthquake

The near-trench coseismic rupture behaviour of the 2011 Tohoku-Oki earthquake remains poorly understood due to the scarcity of near-field observations. Differential bathymetry offers a unique approach to studying offshore coseismic seafloor deformation but has a limited horizontal resolution. Here we use differential bathymetry estimates with improved horizontal resolutions to investigate near-trench coseismic slip behaviours in the 2011 Tohoku-Oki earthquake. In the main rupture region, a velocity-strengthening behaviour in the shallow fault is observed. By contrast, the seafloor uplift decreases towards the trench, but the trend inverts near the backstop interface outcrop, revealing significant off-fault deformation features. Amongst various competing off-fault effects observed, we suggest that inelastic deformation plays a predominant role in near-trench tsunami excitation. Large trench-bleaching rupture is also observed immediately north of 39°, delimiting the northern extent of the main rupture region. Overall, striking spatial heterogeneity of the shallow rupture behaviour is revealed for the region.

Scour ponds from unusually large tsunamis on a beach-ridge plain in eastern Hokkaido, Japan

Scour ponds from unusually large tsunamis cut across the crest of a beach ridge in Kiritappu marsh, eastern Hokkaido. No fewer than ten of these ponds were imaged by photogrammetry as elongate topographic depressions as large as 5 m by 30 m. Sediments in these ponds are underlain by unconformities that were detected with ground-penetrating radar and observed directly in cores and a slice sample. Sediment deposits in the ponds contain peat and volcanic ash layers, the ages of which suggest that the scouring occurred during tsunamis generated by spatially extensive thrust ruptures along the southern Kuril trench, most recently during the early seventeenth century and its predecessor during the thirteenth-fourteenth century. Some of the ponds appear to have been formed during one tsunami and refreshed during later successors. This evidence of recurrent erosion suggests that the shoreline may retreat as part of earthquake-related cycles of coastal uplift and subsidence.

Integrated mechanical environment of pre- and post-rupture fault and asperity origin of the 2011 giant Tohoku-Oki earthquake

It is a key to know mechanical environment (ME) of pre- and post-rupture fault of giant earthquakes at subduction zones for predicting earthquake and tsunami disaster. However, we know little about its details till now. In this paper, using the inverted stress change three hours before and three hours after the mainshock in the rupture zone of the 2011 Tohoku-Oki M_w 9.0 earthquake, we show a quantitative integrated ME in the rupture zone, including principal stress, pore-fluid pressure and friction strength. We discover from this environment a large asperity composed of two asperities induced by relatively high friction coefficients and relatively lower pore-fluid pressures. The integrate ME quantitatively explained the reasons of the overshoot and relatively lower shear strength of the trench, which caused huge displacement and tsunami at the trench. We suggest that the asperities favor the horst and graben structure system which provides a geology environment for interseismic stress accumulation and thus for breeding the megathrust tsunami earthquake.

Palaeo-Tsunami Events on the Coasts of Cyprus

Cyprus has a long history of tsunami activity, as described in archaeological and geological records. Although the study area has experienced tsunamis in the past and constitutes an area threatened by this hazard both from the Cyprean arc and from the neighboring Hellenic arc, field research on tsunami evidence on the coastal zone of Cyprus still remains scarce. It is clear from the literature that large boulder accumulations are an important feature along the coasts of Cyprus, testifying to extreme events. A detailed field survey revealed that at various locations cited in the literature as hosting geomorphological evidence of past tsunamis, no such evidence was identified. It is likely that the high touristic activity that has been occurring on the coasts of Cyprus during the last 20 years may have affected tsunami indicators such as boulder accumulations. Tsunamis are unpredictable and infrequent but potentially large-impact natural disasters. The latest strong tsunami that caused damage to the Cypriot coast was centuries ago, when the population and economic growth and development at the Cypriot shoreline did not exist. Today, the coastal zone hosts a higher population as well as increasing touristic activity, highlighting the need for better preparedness, awareness raising and for tsunami-related risk reduction.

Influence of normal tide and the Great Tsunami as recorded through hourly-resolution micro-analysis of a mussel shell

We report here hourly variations of Mg/Ca, Sr/Ca, and Ba/Ca ratios in a Mediterranean mussel shell (Mytilus galloprovincialis) collected at the Otsuchi bay, on the Pacific coast of northeastern Japan. This bivalve was living in the intertidal zone, where such organisms are known to form a daily or bidaily growth line comprised of abundant organic matter. Mg/Ca ratios of the inner surface of the outer shell layer, corresponding to the most recent date, show cyclic changes at 25-90 μm intervals, while no interpretable variations are observed in Sr/Ca and Ba/Ca ratios. High Mg/Ca ratios were probably established by (1) cessation of the external supply of Ca and organic layer forming when the shell is closed at low tide, and (2) the strong binding of Mg to the organic layer, but not of Sr and Ba. Immediately following the great tsunami induced by the 2011 Tohoku earthquake, Mg/Ca enrichment occurred, up to 10 times that of normal low tide, while apparent Ba/Ca enrichment was observed for only a few days following the event, therefore serving a proxy of the past tsunami. Following the tsunami, periodic peaks and troughs in Mg/Ca continued, perhaps due to a biological memory effect as an endogenous clock.

Decaying Post-Seismic Deformation Observed on the Korean Peninsula Following the 2011 Tohoku-Oki Earthquake

We investigated decaying post-seismic deformation observed on the Korean Peninsula associated with the 2011 Mw 9.0 Tohoku-Oki earthquake using Global Navigation Satellite System (GNSS). The GNSS velocity vectors were estimated in five periods from 2005 to 2019. A co-seismic offset of the Korean Peninsula caused by the 2011 earthquake was inversely proportional to epicentral distances. According to the temporal variations of two components (magnitude and direction) of the GNSS velocity vector with the epicentral distance, the difference between the eastern and western regions for the two components becomes smaller over time. For approximately nine years after the 2011 event, the direction for the crustal movement in South Korea showed a recovery pattern returning to the pre-earthquake motion. In addition, the recovery patterns of the crustal movement were observed differently with the regional geologic structure (e.g., the crustal thickness) and each period. Our estimates of the decay in post-seismic deformation of the Korean Peninsula suggest that post-seismic relaxation will be complete within 5–20 years after the 2011 earthquake. The results suggest that the crustal movement on the Korean Peninsula is gradually recovering to its pre-earthquake motion.

Investigating a tsunamigenic megathrust earthquake in the Japan Trench

The 2011 Tohoku-oki earthquake occurred in the Japan Trench 10 years ago, where devastating earthquakes and tsunamis have repeatedly resulted from subduction of the Pacific plate. Densely instrumented seismic, geodetic, and tsunami observation networks precisely recorded the event, including seafloor observations. A large coseismic fault slip that unexpectedly extended to a shallow part of megathrust fault was documented. Strong lateral variations of the coseismic slip near the trench were recorded from marine geophysical studies, along with a possible cause of these variations. The seismic activities in east Japan are still higher than those before the earthquake, and crustal deformation is still occurring. Although the recurrence probability of a great earthquake (magnitude = ~9) in the Japan Trench in the near future is very low, a large normal fault earthquake seaward of the Japan Trench is a concerning possibility.

Near-surface softening and healing in eastern Honshu associated with the 2011 magnitude-9 Tohoku-Oki Earthquake

The near-surface part of the crust, also called the skin of the earth, is the arena of human activity of which the stiffness is of great concern to engineers in infrastructure construction. The stiffness reduction of near-surface geomaterials also plays a vital role in geohazards triggering. However, the physical mechanism behind the material softening is still not fully understood. Here, we report a coseismic shear-wave velocity reduction in the near surface by up to a few tens of percent during the strongest shaking from the 11 March 2011 Tohoku-Oki Earthquake and a subsequent two-stage healing process including a rapid recovery within a few minutes and a slow recovery over many years. We also present a theoretical contact model between mineral grains in geomaterials containing multiple metastable contacts at small separations due to the oscillatory hydration interaction, which can explain the emergence of different stages in the healing process.

Quantifying the inundation impacts of earthquake-induced surface elevation change by hydrological and hydraulic modeling

Current estimates of flood hazards are often based on the assumption that topography is static. When tectonic and/or anthropogenic processes change the land surface elevation, the spatial patterns of floods might also change. Here, we employ the hydrological and hydraulic modeling to simulate floods in the Kujukuri Plain, Japan, in the years 2004 and 2013, when two severe floods occurred. In between the two floods, land surface elevations were changed by the 2011 Tohoku-Oki earthquake. The effects of land surface elevation changes on inundation areas were quantified by changing input topographies. Our results showed that, without taking into account land surface elevation changes, around 10% of inundation areas were underestimated at the time of flood events in the year 2013. The spatial distribution of inundation locations varied with local topographical features, for example, the areas with backmarsh and valley fill deposits were sensitive to the extent of inundation by land surface elevation changes. The sub-watershed near the coastal shoreline having below-zero meter elevation areas showed that the earthquake-induced land surface elevation changes exacerbated an additional 22% inundation area. This study suggests that the inundation areas will increase in catchments suffering severe settlements, which highlights the necessity of taking into account the spatio-temporal changes of land surface elevations on the assessment of flood hazards.

Increasing benthic vent formation: a threat to Japan's ancient lake

An autonomous underwater vehicle (AUV) was deployed in Lake Biwa from 2000 to 2012. In December 2009, ebullition of turbid water was first found in the deepest area (> 90 m) of the North Basin. Follow-up investigations in April and December 2010 and January 2012 confirmed the existence of benthic vents similar to the vents observed in other deep lakes. Importantly, vent numbers per unit travel distance in Lake Biwa dramatically increased from only two vents (0.37 vents km⁻¹) in December 2009 to 54 vents (5.28 vents km⁻¹) in January 2012, which could be related to recent tectonic activity in Japan, e.g., the M9.1 Tohoku earthquake in March 2011 and slow earthquakes along the Nankai Trough from 2006 to 2018. Continuous back-up investigations from 2014 to 2019 revealed additional benthic vents in the same area. The sudden increase in benthic vent activity (liquid and gaseous ebullitions) have significant potential to alter lake biogeochemistry and, ultimately, degrade Japan’s major drinking water source and may be a harbinger of major crustal change in the near future.

Rheology of the Zagros Lithosphere from Post-Seismic Deformation of the 2017 Mw7.3 Kermanshah, Iraq, Earthquake

We use 2018–2020 Sentinel-1 InSAR time series data to study post-seismic deformation processes following the 2017 Mw 7.3 Kermanshah, Iraq earthquake. We remove displacements caused by two large aftershock sequences from the displacement field. We find that for a six month period the response is dominated by afterslip along the up-dip extension of the coseismic rupture zone, producing up to 6 cm of radar line-of-sight displacements. The moment magnitude of afterslip is Mw 5.9 or 12% of the mainshock moment. After that period, the displacement field is best explained by viscoelastic relaxation and a lower crustal viscosity of η l c = 1 − 0.4 + 0.8 × 10 19 Pas . The viscosity of the uppermost mantle is not constrained by the data, except that it is larger than 0.6 × 10 19 Pas . The relatively high lower crustal and uppermost mantle viscosities are consistent with a cold and dry lithosphere of the Zagros region.

Tomography of the source zone of the great 2011 Tohoku earthquake

The mechanism and rupture process of the giant 2011 Tohoku-oki earthquake (Mw 9.0) are still poorly understood due to lack of permanent near-field observations. Using seismic arrival times recorded by dense seismograph networks on land and at ocean floor, we determine a detailed seismic tomography model of the megathrust zone beneath the Tohoku forearc. Our results show that the coseismic slip of the 2011 Tohoku-oki earthquake initiated at a boundary between a down-dip high-velocity anomaly and an up-dip low-velocity anomaly. The slow anomaly at shallow depths near the Japan trench may reflect low-rigidity materials that are close to the free surface, resulting in large slip and weak high-frequency radiation. Our new tomographic model can account for not only large slip near the trench but also weak high-frequency radiation from the shallow rupture areas.

Using data recorded by a new seafloor seismic network, the authors reveal the detailed 3D structure of the source zone of the great 2011 Tohoku-oki earthquake, which sheds new light on the mechanism of the great earthquake and tsunami.

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.

Advances in gravity analyses for studying volcanoes and earthquakes

This report highlights the usefulness and applicability of various gravimetric methods for studying earthquakes and volcanic activities. A high-resolution gravity anomaly map of Japan reveals areas with very steep horizontal gradients, where potential seismic faults are likely to be buried. Such traditional geoprospecting is coupled with novel cosmic-ray radiography to produce a fine-resolution (<100 m) three-dimensional density structure of a volcano. On the other hand, temporal gravity changes provide invaluable information about the process of earthquake faulting, volcanic eruptions, caldera formation, etc. Specifically, in this report we present our previous work on gravity research for solid earth science: (1) the first detection of coseismic gravity changes, (2) the virtual visualization of the rising and falling of magma in a conduit of Asama volcano, and (3) the large-scale lateral movement of magma during the Miyake-jima eruption in 2000.

Rapid mantle flow with power-law creep explains deformation after the 2011 Tohoku mega-quake

The deformation transient following large subduction zone earthquakes is thought to originate from the interaction of viscoelastic flow in the asthenospheric mantle and slip on the megathrust that are both accelerated by the sudden coseismic stress change. Here, we show that combining insight from laboratory solid-state creep and friction experiments can successfully explain the spatial distribution of surface deformation in the first few years after the 2011 M_w 9.0 Tohoku-Oki earthquake. The transient reduction of effective viscosity resulting from dislocation creep in the asthenosphere explains the peculiar retrograde displacement revealed by seafloor geodesy, while the slip acceleration on the megathrust accounts for surface displacements on land and offshore outside the rupture area. Our results suggest that a rapid mantle flow takes place in the asthenosphere with temporarily decreased viscosity in response to large coseismic stress, presumably due to the activation of power-law creep during the post-earthquake period.

Earthquake segmentation in northern Chile correlates with curved plate geometry

We performed an integrated analysis of the coseismic slip, afterslip and aftershock activity of the 2014 M_w 8.1 Pisagua earthquake. This earthquake seems to be spatially located between two major historical earthquakes, the 1868 M_w 8.8 earthquake in southern Peru and the 1877 M_w 8.5 earthquake in northern Chile. Continuous GPS data were used to model the coseismic slip of the mainshock and the largest aftershock (M_w 7.6). The afterslip was modeled for 273 days (end of year 2014) after the largest aftershock, revealing two patches of afterslip: a southern patch between the mainshock and the largest aftershock and a patch to the north of the mainshock. Observations from the seismic network indicate that aftershocks were concentrated near the southern patch. Conversely, the northern patch contained hardly any aftershocks, indicating a dominant aseismic slip. The Pisagua earthquake occurred within a prominent, curved section of the Andean subduction zone. This section may have acted as a barrier for the largest historical earthquakes and as an isolated segment during the Pisagua earthquake.

Effects of Spatiotemporal Filtering on the Periodic Signals and Noise in the GPS Position Time Series of the Crustal Movement Observation Network of China

Analysis of Global Positioning System (GPS) position time series and its common mode components (CMC) is very important for the investigation of GPS technique error, the evaluation of environmental loading effects, and the estimation of a realistic and unbiased GPS velocity field for geodynamic applications. In this paper, we homogeneously processed the daily observations of 231 Crustal Movement Observation Network of China (CMONOC) Continuous GPS stations to obtain their position time series. Then, we filtered out the CMC and evaluated its effects on the periodic signals and noise for the CMONOC time series. Results show that, with CMC filtering, peaks in the stacked power spectra can be reduced at draconitic harmonics up to the 14th, supporting the point that the draconitic signal is spatially correlated. With the colored noise suppressed by CMC filtering, the velocity uncertainty estimates for both of the two subnetworks, CMONOC-I (≈16.5 years) and CMONOC-II (≈4.6 years), are reduced significantly. However, the CMONOC-II stations obtain greater reduction ratios in velocity uncertainty estimates with average values of 33%, 38%, and 54% for the north, east, and up components. These results indicate that CMC filtering can suppress the colored noise amplitudes and improve the precision of velocity estimates. Therefore, a unified, realistic, and three-dimensional CMONOC GPS velocity field estimated with the consideration of colored noise is given. Furthermore, contributions of environmental loading to the vertical CMC are also investigated and discussed. We find that the vertical CMC are reduced at 224 of the 231 CMONOC stations and 170 of them are with a root mean square (RMS) reduction ratio of CMC larger than 10%, confirming that environmental loading is one of the sources of CMC for the CMONOC height time series.

Acoustic Emission/Seismicity at Depth Beneath an Artificial Lake after the 2011 Tohoku Earthquake

Acoustic emission (AE)/seismicity activity increased near the city of Sendai, Japan, after the 11 March 2011 Tohoku earthquake in a newly seismically active region near the Nagamachi-Rifu fault, which caused a magnitude 5.0 earthquake in 1998. The source of this activity was around 12 km beneath an artificial lake. At the same time, activity on the Nagamachi-Rifu fault nearly ceased. More than 1550 micro-earthquakes were observed between 11 March 2011 and 1 August 2012, of which 63% exhibited similar waveforms and defined 64 multiplets. It appears that crustal extension of about 2 m during the Tohoku earthquake and additional extension of about 1 m during the following year changed the stress field in this region, thus generating micro-earthquakes and controlling their frequency. However, it has been presumed that crustal movement during the Tohoku earthquake did not affect the direction of principal stress, and that these events induced repeated quasi-static slips at asperities and the resultant micro-earthquakes.

Groundwater oxygen isotope anomaly before the M6.6 Tottori earthquake in Southwest Japan

Geochemical monitoring of groundwater in seismically-active regions has been carried out since 1970s. Precursors were well documented, but often criticized for anecdotal or fragmentary signals, and for lacking a clear physico-chemical explanation for these anomalies. Here we report - as potential seismic precursor - oxygen isotopic ratio anomalies of +0.24‰ relative to the local background measured in groundwater, a few months before the Tottori earthquake (M 6.6) in Southwest Japan. Samples were deep groundwater located 5 km west of the epicenter, packed in bottles and distributed as drinking water between September 2015 and July 2017, a time frame which covers the pre- and post-event. Small but substantial increase of 0.07‰ was observed soon after the earthquake. Laboratory crushing experiments of aquifer rock aimed to simulating rock deformation under strain and tensile stresses were carried out. Measured helium degassing from the rock and ¹⁸O-shift suggest that the co-seismic oxygen anomalies are directly related to volumetric strain changes. The findings provide a plausible physico-chemical basis to explain geochemical anomalies in water and may be useful in future earthquake prediction research.

A b map implying the first eastern rupture of the Nankai Trough earthquakes

The Nankai Trough megathrust earthquakes inflicted catastrophic damage on Japanese society and more widely. Most research is aimed at identifying strongly coupled regions that are considered as a major source of future disastrous earthquakes. Here we present a b-value map for the entire Nankai Trough zone. The b value, which represents the rate of occurrence of small earthquakes relative to larger ones, is inversely dependent on differential stresses, and has been used to detect highly stressed areas on fault planes in various tectonic situations. A remarkable finding is that the b value is inversely correlated with the slip-deficit rate (SDR). Moreover, the b value for the areas of high SDR in the eastern part is lower than that in the western part, indicating that differential stress on asperities in the eastern part is higher than that in the western part. This may explain the history of the Nankai Trough earthquakes, in which the eastern part tends to rupture first.

Earthquakes generated from the Nankai Trough have caused much devastation over the years. Here, the authors present a b-value map for the Nankai Trough zone, where the Eastern part of the trough has lower b-values than the West, which may help to explain why the Eastern part tends to rupture first.

Intermittent Granular Dynamics at a Seismogenic Plate Boundary

Earthquakes at seismogenic plate boundaries are a response to the differential motions of tectonic blocks embedded within a geometrically complex network of branching and coalescing faults. Elastic strain is accumulated at a slow strain rate on the order of 10^{-15}  s^{-1}, and released intermittently at intervals >100  yr, in the form of rapid (seconds to minutes) coseismic ruptures. The development of macroscopic models of quasistatic planar tectonic dynamics at these plate boundaries has remained challenging due to uncertainty with regard to the spatial and kinematic complexity of fault system behaviors. The characteristic length scale of kinematically distinct tectonic structures is particularly poorly constrained. Here, we analyze fluctuations in Global Positioning System observations of interseismic motion from the southern California plate boundary, identifying heavy-tailed scaling behavior. Namely, we show that, consistent with findings for slowly sheared granular media, the distribution of velocity fluctuations deviates from a Gaussian, exhibiting broad tails, and the correlation function decays as a stretched exponential. This suggests that the plate boundary can be understood as a densely packed granular medium, predicting a characteristic tectonic length scale of 91±20  km, here representing the characteristic size of tectonic blocks in the southern California fault network, and relating the characteristic duration and recurrence interval of earthquakes, with the observed sheared strain rate, and the nanosecond value for the crack tip evolution time scale. Within a granular description, fault and blocks systems may rapidly rearrange the distribution of forces within them, driving a mixture of transient and intermittent fault slip behaviors over tectonic time scales.

Along-trench variation in seafloor displacements after the 2011 Tohoku earthquake

The 2011 Tohoku-oki earthquake was the largest earthquake ever observed with seafloor geodetic techniques in and around its source region. Large crustal deformation associated with both the coseismic rupture and the rapid postseismic deformation has been reported. However, these observations are insufficient to describe the postseismic deformation processes occurring around the broad rupture area. We report the first results of seafloor Global Positioning System and acoustic ranging (GPS-A) observations based on repeated campaign surveys conducted over nearly 4 years using the extended GPS-A network deployed along the Japan Trench in September 2012. The observed postseismic displacement rates (DRs) show evident spatial variation along the trench: (i) distinct landward DRs in the large coseismic slip area [primary rupture area (PRA)], evidencing the predominance of viscoelastic relaxation; (ii) remarkable trenchward DRs in the south of the PRA, indicating rapid afterslip; and (iii) slight trenchward DRs in the north of the PRA. These features provide great insights into constructing a more complete model of viscoelastic relaxation, and they also indicate spatial variation of afterslip and fault locking along the plate interface with clear spatial resolution, providing invaluable information for the improvement of seismic hazard assessment.

Coseismic seafloor deformation in the trench region during the Mw8.8 Maule megathrust earthquake

The M_w 8.8 megathrust earthquake that occurred on 27 February 2010 offshore the Maule region of central Chile triggered a destructive tsunami. Whether the earthquake rupture extended to the shallow part of the plate boundary near the trench remains controversial. The up-dip limit of rupture during large subduction zone earthquakes has important implications for tsunami generation and for the rheological behavior of the sedimentary prism in accretionary margins. However, in general, the slip models derived from tsunami wave modeling and seismological data are poorly constrained by direct seafloor geodetic observations. We difference swath bathymetric data acquired across the trench in 2008, 2011 and 2012 and find ~3–5 m of uplift of the seafloor landward of the deformation front, at the eastern edge of the trench. Modeling suggests this is compatible with slip extending seaward, at least, to within ~6 km of the deformation front. After the M_w 9.0 Tohoku-oki earthquake, this result for the Maule earthquake represents only the second time that repeated bathymetric data has been used to detect the deformation following megathrust earthquakes, providing methodological guidelines for this relatively inexpensive way of obtaining seafloor geodetic data across subduction zone.

Possibility of short-term probabilistic forecasts for large earthquakes making good use of the limitations of existing catalogs

Earthquakes are quite hard to predict. One of the possible reasons can be the fact that the existing catalogs of past earthquakes are limited at most to the order of 100 years, while their characteristic time scale is sometimes greater than that time span. Here we rather use these limitations positively and characterize some large earthquake events as abnormal events that are not included there. When we constructed probabilistic forecasts for large earthquakes in Japan based on similarity and difference to their past patterns-which we call known and unknown abnormalities, respectively-our forecast achieved probabilistic gains of 5.7 and 2.4 against a time-independent model for main shocks with the magnitudes of 7 or above. Moreover, the two abnormal conditions covered 70% of days whose maximum magnitude was 7 or above.

Physical applications of GPS geodesy: a review

Geodesy, the oldest science, has become an important discipline in the geosciences, in large part by enhancing Global Positioning System (GPS) capabilities over the last 35 years well beyond the satellite constellation's original design. The ability of GPS geodesy to estimate 3D positions with millimeter-level precision with respect to a global terrestrial reference frame has contributed to significant advances in geophysics, seismology, atmospheric science, hydrology, and natural hazard science. Monitoring the changes in the positions or trajectories of GPS instruments on the Earth's land and water surfaces, in the atmosphere, or in space, is important for both theory and applications, from an improved understanding of tectonic and magmatic processes to developing systems for mitigating the impact of natural hazards on society and the environment. Besides accurate positioning, all disturbances in the propagation of the transmitted GPS radio signals from satellite to receiver are mined for information, from troposphere and ionosphere delays for weather, climate, and natural hazard applications, to disturbances in the signals due to multipath reflections from the solid ground, water, and ice for environmental applications. We review the relevant concepts of geodetic theory, data analysis, and physical modeling for a myriad of processes at multiple spatial and temporal scales, and discuss the extensive global infrastructure that has been built to support GPS geodesy consisting of thousands of continuously operating stations. We also discuss the integration of heterogeneous and complementary data sets from geodesy, seismology, and geology, focusing on crustal deformation applications and early warning systems for natural hazards.

Impacts of the 2011 tsunami on the subtidal polychaete assemblage and the following recolonization in Onagawa Bay, northeastern Japan

The ecological impacts of the 2011 Great East Japan Earthquake and tsunami and the following recolonization of the subtidal benthic polychaete community were examined by monthly pre- and post-quake field surveys that were conducted in Onagawa Bay from 2007 to 2013. Before the tsunami, the species composition in this benthic community was constant and was dominated by cirratulid and magelonid polychaetes. The density and biomass of benthic polychaetes drastically decreased after the tsunami, and the polychaete community fluctuated during the 2 years after the natural disaster. Spionid and capitellid polychaetes were dominant at this period. In June 2013, the community entered a new constant stage dominated by maldanids, which is different from the pre-quake community. Ecological impacts due to chemical pollution were suggested in addition to the tsunami disturbance. These overlapping effects and physical, chemical and biological factors affected the recovery and recolonization of the polychaete community after the natural disaster.

Aftereffects of Subduction-Zone Earthquakes: Potential Tsunami Hazards along the Japan Sea Coast

The 2011 Tohoku-Oki Earthquake is a typical subduction-zone earthquake and is the 4th largest earthquake after the beginning of instrumental observation of earthquakes in the 19th century. In fact, the 2011 Tohoku-Oki Earthquake displaced the northeast Japan island arc horizontally and vertically. The displacement largely changed the tectonic situation of the arc from compressive to tensile. The 9th century in Japan was a period of natural hazards caused by frequent large-scale earthquakes. The aseismic tsunamis that inflicted damage on the Japan Sea coast in the 11th century were related to the occurrence of massive earthquakes that represented the final stage of a period of high seismic activity. Anti-compressive tectonics triggered by the subduction-zone earthquakes induced gravitational instability, which resulted in the generation of tsunamis caused by slope failing at the arc-back-arc boundary. The crustal displacement after the 2011 earthquake infers an increased risk of unexpected local tsunami flooding in the Japan Sea coastal areas.

Helium anomalies suggest a fluid pathway from mantle to trench during the 2011 Tohoku-Oki earthquake

Geophysical evidence suggests that fluids along fault planes have an important role in generating earthquakes; however, the nature of these fluids has not been well defined. The 2011 magnitude 9.0 Tohoku-Oki earthquake ruptured the interface between the subducting Pacific plate and the overlying Okhotsk plate. Here we report a sharp increase in mantle-derived helium in bottom seawater near the rupture zone 1 month after the earthquake. The timing and location indicate that fluids were released from the mantle on the seafloor along the plate interface. The movement of the fluids was rapid, with a velocity of ~4 km per day and an uncertainty factor of four. This rate is much faster than what would be expected from pressure-gradient propagation, suggesting that over-pressurized fluid is discharged along the plate interface.

Broadscale postseismic gravity change following the 2011 Tohoku-Oki earthquake and implication for deformation by viscoelastic relaxation and afterslip

The analysis of GRACE gravity data revealed postseismic gravity increase by 6 μGal over a 500 km scale within a couple of years after the 2011 Tohoku-Oki earthquake, which is nearly 40-50% of the coseismic gravity change. It originates mostly from changes in the isotropic component corresponding to the M_rr moment tensor element. The exponential decay with rapid change in a year and gradual change afterward is a characteristic temporal pattern. Both viscoelastic relaxation and afterslip models produce reasonable agreement with the GRACE free-air gravity observation, while their Bouguer gravity patterns and seafloor vertical deformations are distinctly different. The postseismic gravity variation is best modeled by the biviscous relaxation with a transient and steady state viscosity of 10¹⁸ and 10¹⁹ Pa s, respectively, for the asthenosphere. Our calculated higher-resolution viscoelastic relaxation model, underlying the partially ruptured elastic lithosphere, yields the localized postseismic subsidence above the hypocenter reported from the GPS-acoustic seafloor surveying.

Disturbance of shallow marine soft-bottom environments and megabenthos assemblages by a huge tsunami induced by the 2011 M9.0 Tohoku-Oki earthquake

Huge tsunami waves associated with megathrust earthquakes have a severe impact on shallow marine ecosystems. We investigated the impact of a tsunami generated by the 2011 M9.0 Tohoku-Oki earthquake on the seafloor and large benthic animals in muddy and sandy ria coasts (Otsuchi and Funakoshi bays) in northeastern Japan. We conducted underwater field surveys using scuba equipment in water depths of <20 m before the tsunami (September 2010) and after the tsunami (September 2011 and September 2012). During the study period, episodic changes in topography and grain-size composition occurred on the seafloor of the study area. Megabenthos sampling revealed a distinct pattern of distribution succession for each benthic species. For example, the protobranch bivalve Yoldia notabilis (Bivalvia: Nuculanidae) and the heterodont bivalve Felaniella usta (Bivalvia: Ungulinidae) disappeared after the tsunami event, whereas the distribution of the venus clam Gomphina melanaegis (Bivalvia: Veneridae) remained unchanged. In addition, the patterns of succession for a single species, such as the giant button top shell Umbonium costatum (Gastropoda: Trochidae) and the heart urchin Echinocardium cordatum (Echinoidea: Loveniidae), varied between the two bays studied. Our data also show that reestablishment of some benthic animal populations began within 18 months of the tsunami disturbance.

Radioactivity analysis following the Fukushima Dai-ichi nuclear accident

A total of 118 samples were analyzed using HPGe γ-spectrometry. (131)I, (134)Cs, (137)Cs and (136)Cs were detected in aerosol air samples that were collected 22 days after the accident with values of 1720 µBq m(-)³, 247 µBq m(-)³, 289 µBq m(-)³ and 23 µBq m(-)³, respectively. (131)I was detected in rainwater and soil samples and was also measurable in vegetables collected between April 2 and 13, 2011, with values ranging from 0.55 Bq kg(-1) to 2.68 Bq kg(-1). No (131)I was detected in milk, drinking water, seawater or marine biota samples.

Compressive sensing of frequency-dependent seismic radiation from subduction zone megathrust ruptures

Megathrust earthquakes rupture a broad zone of the subducting plate interface in both along-strike and along-dip directions. The along-dip rupture characteristics of megathrust events, e.g., their slip and energy radiation distribution, reflect depth-varying frictional properties of the slab interface. Here, we report high-resolution frequency-dependent seismic radiation of the four largest megathrust earthquakes in the past 10 y using a compressive-sensing (sparse source recovery) technique, resolving generally low-frequency radiation closer to the trench at shallower depths and high-frequency radiation farther from the trench at greater depths. Together with coseismic slip models and early aftershock locations, our results suggest depth-varying frictional properties at the subducting plate interfaces. The shallower portion of the slab interface (above ∼15 km) is frictionally stable or conditionally stable and is the source region for tsunami earthquakes with large coseismic slip, deficient high-frequency radiation, and few early aftershocks. The slab interface at intermediate depths (∼15–35 km) is the main unstable seismogenic zone for the nucleation of megathrust quakes, typically with large coseismic slip, abundant early aftershocks, and intermediate- to high-frequency radiation. The deeper portion of the slab interface (∼35–45 km) is seismically unstable, however with small coseismic slip, dominant high-frequency radiation, and relatively fewer aftershocks.

Observation of Earth's free oscillation by dense GPS array: after the 2011 Tohoku megathrust earthquake

Day-scale Earth's free oscillation after large earthquakes has been detected by underground instruments such as strainmeters, gravimeters and seismometers, to investigate Earth's internal structure, geodynamics, and source properties of earthquakes. Here we show that Global Positioning System (GPS) can also detect the signals of the Earth's free oscillation. A dense GPS array in Japan (GEONET) recorded the surface deformation following the 2011 Tohoku megathrust earthquake. A simple array analysis over 300 stations reduces local noise in GPS time series. We find that the dense GPS array truly detected both spheroidal and toroidal fundamental modes in three-direction displacement. This new tool has a strong potential to investigate the free oscillations particularly in low-frequency bands.

Patient handling system for carbon ion beam scanning therapy

Our institution established a new treatment facility for carbon ion beam scanning therapy in 2010. The major advantages of scanning beam treatment compared to the passive beam treatment are the following: high dose conformation with less excessive dose to the normal tissues, no bolus compensator and patient collimator/ multi‐leaf collimator, better dose efficiency by reducing the number of scatters. The new facility was designed to solve several problems encountered in the existing facility, at which several thousand patients were treated over more than 15 years. Here, we introduce the patient handling system in the new treatment facility. The new facility incorporates three main systems, a scanning irradiation system (S‐IR), treatment planning system (TPS), and patient handling system (PTH). The PTH covers a wide range of functions including imaging, geometrical/position accuracy including motion management (immobilization, robotic arm treatment bed), layout of the treatment room, treatment workflow, software, and others. The first clinical trials without respiratory gating have been successfully started. The PTH allows a reduction in patient stay in the treatment room to as few as 7 min. The PTH plays an important role in carbon ion beam scanning therapy at the new institution, particularly in the management of patient handling, application of image‐guided therapy, and improvement of treatment workflow, and thereby allows substantially better treatment at minimum cost.

PACS numbers: 87.56.‐v; 87.57.‐s; 87.55.‐x

An overview of respiratory medicine during the Tsunami Disaster at Tohoku, Japan, on March 11, 2011

A 9.0-magnitude earthquake, with an epicenter 150 km east of the Tohoku area of Japan in the deep Pacific Ocean, triggered an unexpectedly huge tsunami and caused 19,000 casualties along 500 km of the Pacific coast in northern Japan, as documented by Ozawa et al. [1]. This brief review provides an overview of the tsunami disaster and ensuing respiratory medical conditions, from emergency rescue conditions to chronic stage diseases, leading to a series of reviews, original articles, and case reports. Successive mega-earthquakes in the previous decade around the world may suggest a shift in the Earth's geological state from the stable to the active stage. The "tsunami lung" term is not limited to the near-drowning lung; inhalation of the sandy wave or of the sludge and slime near industrialized areas can cause inflammation and rare semi-acute phase fungal infections. Although the long-term outcomes of tsunami-related respiratory medicine need further analysis, determining how to reduce the extent of damage is a critical and central issue.

Recurrent slow slip event likely hastened by the 2011 Tohoku earthquake

Slow slip events (SSEs) are another mode of fault deformation than the fast faulting of regular earthquakes. Such transient episodes have been observed at plate boundaries in a number of subduction zones around the globe. The SSEs near the Boso Peninsula, central Japan, are among the most documented SSEs, with the longest repeating history, of almost 30 y, and have a recurrence interval of 5 to 7 y. A remarkable characteristic of the slow slip episodes is the accompanying earthquake swarm activity. Our stable, long-term seismic observations enable us to detect SSEs using the recorded earthquake catalog, by considering an earthquake swarm as a proxy for a slow slip episode. Six recurrent episodes are identified in this way since 1982. The average duration of the SSE interoccurrence interval is 68 mo; however, there are significant fluctuations from this mean. While a regular cycle can be explained using a simple physical model, the mechanisms that are responsible for the observed fluctuations are poorly known. Here we show that the latest SSE in the Boso Peninsula was likely hastened by the stress transfer from the March 11, 2011 great Tohoku earthquake. Moreover, a similar mechanism accounts for the delay of an SSE in 1990 by a nearby earthquake. The low stress buildups and drops during the SSE cycle can explain the strong sensitivity of these SSEs to stress transfer from external sources.

Clues from joint inversion of tsunami and geodetic data of the 2011 Tohoku-oki earthquake

The 2011 Tohoku-oki (Mw 9.1) earthquake is so far the best-observed megathrust rupture, which allowed the collection of unprecedented offshore data. The joint inversion of tsunami waveforms (DART buoys, bottom pressure sensors, coastal wave gauges, and GPS-buoys) and static geodetic data (onshore GPS, seafloor displacements obtained by a GPS/acoustic combination technique), allows us to retrieve the slip distribution on a non-planar fault. We show that the inclusion of near-source data is necessary to image the details of slip pattern (maximum slip ~48 m, up to ~35 m close to the Japan trench), which generated the large and shallow seafloor coseismic deformations and the devastating inundation of the Japanese coast. We investigate the relation between the spatial distribution of previously inferred interseismic coupling and coseismic slip and we highlight the importance of seafloor geodetic measurements to constrain the interseismic coupling, which is one of the key-elements for long-term earthquake and tsunami hazard assessment.