Suito H.,Geospatial Information Authority of Japan
Earth, Planets and Space | Year: 2017
This study develops a three-dimensional viscoelastic model using the finite element method to understand the postseismic deformation that followed the 2011 Tohoku-Oki earthquake. The question of understanding which elements of the viscoelastic media affect the surface deformation is of particular importance. We first examined the individual effects of two different viscoelastic media, the mantle wedge and the oceanic mantle, which produce almost opposite deformation patterns. The mantle wedge controls eastward motion, uplift of the Pacific coastal and offshore regions, and extension across a broad area. In contrast, the oceanic mantle controls dominantly offshore westward motion, subsidence across a broad area, minor uplift of the surrounding areas, and contraction offshore. These differences are the most important issues for understanding the viscoelastic relaxation caused by subduction earthquakes. We then developed four different models to clarify which elements of the viscoelastic media affect the observed surface deformation. The simplest model, with uniform viscosity for all viscoelastic media, could explain the horizontal deformation but not the vertical deformation. The second model, with different viscosities for the mantle wedge and the oceanic mantle, could explain the onshore observations but could not explain the seafloor observations. The third model, which includes a thin weak layer beneath the subducting slab, could essentially explain the near-field onshore and seafloor observations but could not explain the far-field data. The final depth-dependent model was able to explain the far-field data as well as the near-field data. In these typical models, it is of particular importance to consider the different viscosities between the mantle wedge and the oceanic mantle and to include a thin weak layer beneath the slab, which has a dramatic impact on the seafloor deformation. Far-field data as well as near-field data are also important for constraining the viscoelastic structure; the former is sensitive to viscoelastic relaxation at greater depths. Clearly, viscoelastic relaxation alone cannot explain the observed deformation. A combined viscoelastic and afterslip model is necessary for constructing a complete postseismic deformation model.[Figure not available: see fulltext.] © 2017 The Author(s).
Ozawa S.,Geospatial Information Authority of Japan
Geophysical Research Letters | Year: 2014
A slow slip event occurred off the coast of the Boso peninsula, Japan, from approximately 28 December 2013 to 10 January 2014. The estimated aseismic slip expanded slightly southward and westward over time with a moment magnitude of 6.5, which was the smallest value since 1996. The recurrence interval has decreased from approximately 6.4 to 2.2years from 1996 to 2014. One explanation of this shortening is the change in Coulomb failure stress due to the 2011 Mw9.0 Tohoku earthquake and its afterslip. Another interpretation is related to a scenario observed in several numerical simulation studies, in which the recurrence interval of slow slip becomes shorter as the time nears a large earthquake. This case will constrain the physical processes of a slip cycle. The Boso slow slip events together with the Tohoku earthquake and its afterslip changed the stress state for the anticipated interplate earthquake near the Sagami trough. Key Points Spatial and temporal slow slip process in Boso slow slip, Japan Shortening of the recurrence interval DCFS is increasing after the Tohoku earthquake ©2014. American Geophysical Union. All Rights Reserved.
Kobayashi T.,Geospatial Information Authority of Japan |
Morishita Y.,Geospatial Information Authority of Japan |
Yarai H.,Geospatial Information Authority of Japan
Earth, Planets and Space | Year: 2015
We have successfully detected widely distributed ground displacements for the 2015 Gorkha earthquake by applying a ScanSAR-based interferometry analysis of Advanced Land Observing Satellite 2 (ALOS-2) L-band data. A major displacement area extends with a length of about 160 km in the east-west direction, and the most concentrated crustal deformation with ground displacement exceeding 1 m is located 20-30 km east from Kathmandu. A quasi-vertical displacement estimated by combining the ascending and the descending data indicates upheaval of about 1.4 m at maximum. We inverted the synthetic aperture radar interferometry (InSAR) data including both of the main shock (moment magnitude (Mw) 7.8) and the largest aftershock (Mw 7.3) to construct a slip distribution model. Our model shows a nearly pure reverse fault motion with a slip amount of approximately 6 m at maximum, and the spatial extent is zonally distributed within a distance of 50 to 100 km from the surface along downdip direction. The downdip end of the slip is quite consistent with that of the interseismic coupling area geodetically inferred in previous studies. On the other hand, there is no significant slip at shallow depth in spite of the fact that the plate interface is thought to be fully locked there, may be suggesting that there still remains a potential of fault slip. The slip distribution unnaturally bifurcates in the east, and we can identify a clear-cut slip deficit area with a radius of ~10 km just west side of the Mw 7.3 event, where the slip amount reaches only 20 cm at most. This area is presumably subjected to a strong shear stress which should promote a reverse fault slip. There is a possibility to produce a fault slip equivalent to Mw ~7.0 in the future although we do not know if the slip heterogeneity would be smoothed out by a seismic event or an aseismic event. © 2015 Kobayashi et al.
Munekane H.,Geospatial Information Authority of Japan
Journal of Geodesy | Year: 2013
The sub-daily noise in horizontal global positioning system (GPS) kinematic time series arising from monument tilts is quantitatively evaluated using tiltmeter data at GPS stations from the Japanese nationwide global navigation satellite system network. The estimated tilt-induced monument displacements show characteristics that are typical of those caused by thermal tilts of the monuments. The root mean square of the displacements is typically a few millimetres, with notable inter-seasonal variations. The stacked amplitude spectra of the monument displacements have peaks at the tidal bands S1 and S2, and their higher tones. The peaks at the S1 and S2 bands in the amplitude spectra are reduced by 41 and 43 % for the north-south component and 36 and 53 % for the east-west component, respectively, after correcting for the monument displacements. The monument displacements due to the thermal tilts of the monuments may also be a favourable candidate for sub-daily noise at the S1 and S2 bands found in other GPS networks. © 2013 Springer-Verlag Berlin Heidelberg.
Yarai H.,Geospatial Information Authority of Japan |
Ozawa S.,Geospatial Information Authority of Japan
Journal of Geophysical Research: Solid Earth | Year: 2013
The time evolution of afterslip on a plate boundary experiencing interplate earthquakes is expected to show logarithmic decay. The global positioning system network in Japan has been monitoring transient deformation since the occurrence of two large interplate earthquakes with moment magnitudes of 6.8 and 6.7 in the Hyuga-nada area, southwest Japan, in 1996. The spatial and temporal evolution of aseismic interplate slip based on crustal deformation data indicates that afterslip followed the two earthquakes and gradually declined to background rates by around 2004 with total moment magnitude of 7.3. However, quasi-periodic slow slip events suddenly began within the afterslip area in 2005 with approximately one year duration and two-year recurrence interval. The moment magnitudes of the three slow slip events since January 2005 range from 6.7 to 6.8. This differs greatly from the expected behavior of logarithmic decay over time. Both velocity-strengthening and velocity-weakening rate-and-state modes have been implicated as the cause of afterslip, whose location is complementary to the main shock area of velocity-weakening, while a slow slip event occurs in the velocity-weakening area with different frictional properties from those of an afterslip area. In light of the seemingly different frictional properties, the coexistence of afterslip and slow slip events in the same area would provide additional information about precisely how the plate interface is behaving. The monitoring of these slow slip events should give the clues to understanding the coexistence of long-term afterslip and slow slip events and the increasing risk of earthquakes in neighboring areas. © 2013. American Geophysical Union. All Rights Reserved.
Imakiire T.,Geospatial Information Authority of Japan |
Koarai M.,Geospatial Information Authority of Japan
Soils and Foundations | Year: 2012
This paper gives an overview of the crustal deformation caused by "the 2011 off the Pacific Coast of Tohoku Earthquake" (hereinafter: "the Tohoku Earthquake") as detected by the GPS Earth Observation Network (GEONET), a GPS continuous observation system operated by the Geospatial Information Authority of Japan (GSI), and of ground deformation resulting from liquefaction triggered by the earthquake motion as determined by a leveling survey. A very large part of the Japanese archipelago was affected by the crustal deformation caused by the mainshock. The subsidence in the Pacific coastal area of the Tohoku region is especially remarkable. Based on a resurvey after the earthquake, the GSI revised the coordinates and heights of survey marks, including the Origin of the Japanese Horizontal Control Network and the Origin of the Japanese Vertical Control Network. We estimated the geometry of the seismogenic fault of the Tohoku Earthquake, as well as a slip model for the boundary between the Pacific Plate and the North American Plate, based on crustal deformation data. Postseismic crustal deformation was also observed by GEONET. Even though the area that coseismicaly subsided is undergoing partial uplift from postseismic deformation, such uplift is not rapid enough to make up for the coseismic subsidence within the next several years. The earthquake caused liquefaction in large parts of the Kanto Region, especially on the Tokyo Bay side and at the lower reaches of the Tone River. Liquefaction was concentrated in areas of land reclamation and former river channels, where it is easy to recognize using time-series geospatial information such as olden topographic maps and olden aerial photographs. To determine the amount of land subsidence (cm) due to liquefaction in Urayasu City, the GSI carried out a leveling survey and compared the differences between pre-and post-earthquake LIDAR survey data with the relative subsidence obtained by leveling. © 2012 The Japanese Geotechnical Society.
Munekane H.,Geospatial Information Authority of Japan
Earth, Planets and Space | Year: 2012
GPS kinematic time series are analyzed to estimate slip distributions for the Mw 9.0 2011 off the Pacific coast of Tohoku Earthquake sequence. Empirical orthogonal function (EOF) analysis is employed to enhance the signal-to-noise ratio of the original time series. The coseismic, and subsequent postseismic, deformations of the foreshock can be described by a single mode, suggesting that the extent of the source for both events must be similar. The total moment magnitude of the afterslip following the foreshock is estimated to be Mw 7.1 with a decay time of 0.63 days. The magnitude of the afterslip was larger for its duration than was anticipated by the scaling law for a typical slow earthquake, although two previous earthquakes in adjacent regions showed the same tendency as that in the present case. The pattern of slip of the mainshock and the subsequent afterslips and aftershocks indicates that each slip occurs in a region adjacent to that of the previous slips in a complementary manner. Finally, in the course of the EOF analysis, the modes representing the thermal expansion of the GPS pillars are clearly identified. © The Society of Geomagnetism and Earth, Planetary and Space Sciences (SGEPSS); The Seismological Society of Japan; The Volcanological Society of Japan; The Geodetic Society of Japan; The Japanese Society for Planetary Sciences; TERRAPUB.
Kuroishi Y.,Geospatial Information Authority of Japan
Nature Photonics | Year: 2016
According to Einstein's theory of relativity, the passage of time changes in a gravitational field. On Earth, raising a clock by 1 cm increases its apparent tick rate by 1.1 parts in 1018, allowing chronometric levelling through comparison of optical clocks. Here, we demonstrate such geopotential measurements by determining the height difference of master and slave clocks separated by 15 km with an uncertainty of 5 cm. A subharmonic of the master clock laser is delivered through a telecom fibre to synchronously operate the distant clocks. Clocks operated under such phase coherence reject clock laser noise and facilitate proposals for linking clocks and interferometers. Taken over half a year, 11 measurements determine the fractional frequency difference between the two clocks to be 1,652.9(5.9) × 10-18, consistent with an independent measurement by levelling and gravimetry. Our system demonstrates a building block for an internet of clocks, which may constitute ‘quantum benchmarks’, serving as height references with dynamic responses. © 2016 Nature Publishing Group
Tobita M.,Geospatial Information Authority of Japan
Earth, Planets and Space | Year: 2016
The time series of a postseismic deformation is commonly fitted by a logarithmic or exponential decay function. However, the high-quality postseismic Global Navigation Satellite System (GNSS) time series of the 2011 Mw 9 Tohoku-Oki earthquake indicates that a single decay function cannot be used to represent the postseismic behaviour. We therefore combined the logarithmic (log) and exponential (exp) decay functions and developed methods for obtaining global solutions using nonlinear least squares calculations for such complex functions. Our models significantly improved the fitting performance of the postseismic time series and the prediction performance of the evolution of postseismic deformation. The solutions obtained by the proposed models and methods enabled distinction between the contributions of the log and exp functions, and explanation of characteristic phenomena such as the subsidence that occurs immediately after an earthquake is reversed to an uplift. The analysis of the solutions may suggest that there has been a continuous increase in the contribution of viscoelastic relaxation to postseismic deformation in eastern Japan, whereas the contribution of afterslip has rapidly decreased. The short-term prediction performance and the universal applicability of the proposed models to the Tohoku-Oki earthquake have contributed to the detection of a slow-slip event in the Tokai region. Rather than the existence of a unique single relaxation time for each surface site, our results suggest a unique single relaxation time for each postseismic deformation mechanism at a given subsurface location. Although the predictions were highly dependent on the assigned steady velocities and the long-term relaxation time constants, they indicate that the coseismic subsidence of the Yamoto station in Miyagi prefecture will recover around the year 2020. The estimated relaxation time constants of the present models appeared to be uniform throughout eastern Japan. © 2016 Tobita.
Morishita Y.,Geospatial Information Authority of Japan |
Hanssen R.F.,Technical University of Delft
IEEE Transactions on Geoscience and Remote Sensing | Year: 2015
Temporal decorrelation is one of the main limitations of synthetic aperture radar (SAR) interferometry. For nonurban areas, its mechanism is very complex, as it is very dependent of vegetation types and their temporal dynamics, actual land use, soil types, and climatological circumstances. Yet, an a priori assessment and comprehension of the expected coherence levels of interferograms are required for designing new satellite missions (in terms of frequency, resolution, and repeat orbits), for choosing the optimal data sets for a specific application, and for feasibility studies for new interferometric applications. Although generic models for temporal decorrelation have been proposed, their parameters depend heavily on the land use in the area of interest. Here, we report the behavior of temporal decorrelation for a specific class of land use: pasture on drained peat soils. We use L-, C-, and X-band SAR observations from the Advanced Land Observation Satellite (ALOS), European Remote Sensing Satellite, Envisat, RADARSAT-2, and TerraSAR-X missions. We present a dedicated temporal decorrelation model using three parameters and demonstrate how coherent information can be retrieved as a function of frequency, repeat intervals, and coherence estimation window sizes. New satellites such as Sentinel-1 and ALOS-2, with shorter repeat intervals than their predecessors, would enhance the possibility to obtain a coherent signal over pasture. © 2014 IEEE.