Field Survey of Non-tectonic Surface Displacements Caused by the 2016 Kumamoto Earthquake Around Aso Valley

Similar documents
Semi-automated landform classification for hazard mapping of soil liquefaction by earthquake

DETECTION OF CRUSTAL DEFORMATION OF THE NORTHERN PAKISTAN EARTHQUAKE BY SATELLITE DATA. Submitted by Japan **

GROUND SURFACE VISUALIZATION USING RED RELIEF IMAGE MAP FOR A VARIETY OF MAP SCALES

NGIAs' roles in successful disaster response

Case Study of Japan: Crustal deformation monitoring with GNSS and InSAR

Mamoru KOARAI and Takayuki NAKANO. (Published online: 18 December 2013) Abstract

Cenozoic Extensional Basin Development and Sedimentation in SW Montana

InSAR-derived Crustal Deformation and Reverse Fault Motion of the 2017 Iran-Iraq Earthquake in the Northwestern Part of the Zagros Orogenic Belt

S. Toda, S. Okada, D. Ishimura, and Y. Niwa International Research Institute of Disaster Science, Tohoku University, Japan

Neotectonic Implications between Kaotai and Peinanshan

Tatsuo Sekiguchi* and Hiroshi Sato*

GEOMORPHOLOGIC RESEARCH USING LIDAR DATA

RESEARCH OF DAMAGED CONDITION BY THE 2016 KUMAMOTO EARTHQUAKE AND GROUND INVESTIGATION ON STONE WALLS AND ERATH STRUCTURES IN KUMAMOTO CASTLE

Fault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction

Continuous Caldera Changes in Miyakejima Volcano after Hiroyuki HASEGAWA, Hiroshi P. SATO and Junko IWAHASHI

68. Izu-Torishima. Summary. Latitude: 30 29'02" N, Longitude: '11" E, Elevation: 394 m (Ioyama) (Elevation Point) (68.

Monthly Volcanic Activity Report (March, 2011)

Monthly Volcanic Activity Report (February 2016)

EXTRACTION OF FLOODED AREAS DUE THE 2015 KANTO-TOHOKU HEAVY RAINFALL IN JAPAN USING PALSAR-2 IMAGES

Latitude: 43 25'03" N, Longitude: '52" E, Elevation: 1,692 m (Maruyama) (Triangulation Point)

Structural Modelling of Inversion Structures: A case study on South Cambay Basin

An analysis on the relationship between land subsidence and floods at the Kujukuri Plain in Chiba Prefecture, Japan

PREDICTION OF AVERAGE SHEAR-WAVE VELOCITY FOR GROUND SHAKING MAPPING USING THE DIGITAL NATIONAL LAND INFORMATION OF JAPAN

NATIONWIDE SITE AMPLIFICATION ZONATION STUDY USING JAPAN ENGINEERING GEOMORPHOLOGIC CLASSIFICATION MAP

32. Hijiori. Summary. (32. Hijiori) Latitude: 38 35'57" N, Longitude: '42" E, Elevation: 552 m (Sankakuyama) (Elevation Point - measured by JMA)

Comparison between predicted liquefaction induced settlement and ground damage observed from the Canterbury earthquake sequence

Case study of Japan: Reference Frames in Practice

Continuously Monitored by JMA. Latitude: 24 45'02" N, Longitude: '21" E, Elevation: 169 m (Suribachiyama) (GSI Measuring Point)

The Effects of Hydraulic Structures on Streams Prone to Bank Erosion in an Intense Flood Event: A Case Study from Eastern Hokkaido

The Geographical features of the Hiroshima landslide disaster triggered by heavy rainfall on August 20, 2014

DETAILED NEAR SURFACE GEOPHYSICAL SURVEY AT AN AREA STRICKEN BY THE 2016 KUMAMOTO EARTHQUAKES. Abstract

Chapter 15 Structures

NCC roles in successful disaster response I.R.IRAN

Applying Hazard Maps to Urban Planning

Challenge to the tunnel construction that supports national resilience

Roles of natural levees in the Ara River alluvial fan on flood management

Supporting the response to the 2018 lower East Rift Zone and summit collapse at Kīlauea Volcano, Hawaiʻi

4. The map below shows a meandering stream. Points A, B, C, and D represent locations along the stream bottom.

Radar Remote Sensing: Monitoring Ground Deformations and Geohazards from Space

The March 11, 2011, Tohoku-oki earthquake (Japan): surface displacement and source modelling

Blind fault Configuration in Osaka, Japan based on 2D gravity inversion

Latitude: 34 31'13" N, Longitude: '45" E, Elevation: 508 m (Miyatsukayama) (Triangulation Point - Toshima)

COSMORPHOLOGY - May 2009

19. Esan Continuously Monitored by JMA

CPO Science Middle School Earth Science Learning System Correlated to Ohio Science Academic Content Standards for Earth Science, grades 6-8

Open Access FULL PAPER. Yukitoshi Fukahata * and Manabu Hashimoto

Date: / Page # Version A Page 1

DAMAGE DETECTION OF THE 2008 SICHUAN, CHINA EARTHQUAKE FROM ALOS OPTICAL IMAGES

4.Mashu. Summary. Latitude: 43 34'20" N, Longitude: '39" E, Elevation: 857 m (Kamuinupuri) (Elevation Point) (4. Mashu)

Monthly Volcanic Activity Report (July, 2012)

RADAR Remote Sensing Application Examples

PALEOGEOGRAPHY of NYS. Definitions GEOLOGIC PROCESSES. Faulting. Folding 9/6/2012. TOPOGRAPHIC RELIEF MAP of NYS GRADATIONAL TECTONIC

89. Wakamiko. Summary. Latitude: ' N, Longitude: ' E, Depth: -77 m (Central cone) (89. Wakamiko)

Structural Styles and Geotectonic Elements in Northwestern Mississippi: Interpreted from Gravity, Magnetic, and Proprietary 2D Seismic Data

Down-stream process transition (f (q s ) = 1)

A physical feature of the Earth s surface

Study on the feature of surface rupture zone of the west of Kunlunshan pass earthquake ( M S 811) with high spatial resolution satellite images

RELATIONSHIP BETWEEN AGE OF GROUND AND LIQUEFACTION OCCURRENCE IN THE 2011 GREAT EAST JAPAN EARTHQUAKE

Map Mysteries Activities

ENGINEER S CERTIFICATION OF FAULT AREA DEMONSTRATION (40 CFR )

WHAT IS THE EARTH MADE OF? LITHOSPHERE AND HYDROSPHERE

Monthly Volcanic Activity Report (August 2015)

22 Bulletin of the Geospatial Information Authority of Japan, Vol.59 December, 2011 Fig. 1 Crustal deformation (horizontal) associated with the 2011 o

Geology 300, Physical Geology Spring 2019 Quiz Ch 19, Plate Tectonics Name

DAMAGE ANALYSIS OF WATER SUPPLY PIPES DUE TO THE 2004 NIIGATA-KEN CHUETSU, JAPAN EARTHQUAKE ABSTRACT

Monthly Volcanic Activity Report (February, 2011)

9. SLOPE AND EMBANKMENT FAILURES

Identification of Lateral Spread Features in the Western New Madrid Seismic Zone J. David Rogers and Briget C. Doyle

PHYSICAL GEOGRAPHY: An overview. Definitions. Faulting. Folding GEOLOGIC PROCESSES 9/17/2013 TOPOGRAPHIC RELIEF

General Geologic Setting and Seismicity of the FHWA Project Site in the New Madrid Seismic Zone

water erosion lithosphere Describe the process of erosion and deposition. chemical weathering Dissolving limestone is an example of.

ESTIMATION OF LANDFORM CLASSIFICATION BASED ON LAND USE AND ITS CHANGE - Use of Object-based Classification and Altitude Data -

Structural Style and Tectonic Evolution of the Nakhon Basin, Gulf of Thailand

Unit 7.2 W.E.D. & Topography Test

Continuously Monitored by JMA. Latitude: 34 23'49" N, Longitude: '13" E, Elevation: 432 m (Miyatsukayama) (Spot elevation measured by JMA)

A. rift valley B. deep trench C. volcanic cone D. uplifted mountain. A. a volcanic island arc B. a large watershed area

GEOL104: Exploring the Planets LAB 5: PLANETARY TECTONICS

Answers: Internal Processes and Structures (Isostasy)

Geologic Landforms Seen on Aerial Photos Instructor Notes

Interpretive Map Series 24

Estimation of Damage Ratios of Wooden Houses and Water Distribution Pipelines in an Earthquake Scenario for the Tokyo Metropolitan Region

The Little Colorado River

Gullfoss. 8.1 Why Has Gullfoss Two Oblique Steps?

Episodic growth of fault-related fold in northern Japan observed by SAR interferometry

Correspondence to: O. Kazaoka Published: 12 November 2015

Crustal deformation by the Southeast-off Kii Peninsula Earthquake

Keywords: Tohoku earthquake, Liquefaction, Hot spot, Pipeline damage, Topographic history

Latitude: 42 49'36" N, Longitude: '41" E, Elevation: 1,898 m (Ezo-Fuji) (Elevation Point)

The Tohoku Earthquake 2011 Case.

Crustal Boundaries. As they move across the asthenosphere and form plate boundaries they interact in various ways. Convergent Transform Divergent

Earth Science, (Tarbuck/Lutgens) Chapter 10: Mountain Building

Aira/Sakura-jima. Kyushu, Japan N, E; summit elev m. All times are local (= UTC + 9 hours)

Lecture # 6. Geological Structures

Urgent Report of the Landslide Disasters by the Mid Niigata Prefecture Earthquake in 2004, Japan

Lecture 9 faults, folds and mountain building

TAKE HOME EXAM 8R - Geology

Special edition paper

4. Geotechnical and Geological Aspects. 4.1 Geotechnical Aspects

Structural Geology Lab. The Objectives are to gain experience

Seismic Activity and Crustal Deformation after the 2011 Off the Pacific Coast of Tohoku Earthquake

Transcription:

Field Survey of Non-tectonic Surface Displacements Caused by the 2016 Kumamoto Earthquake Around Aso Valley 47 Field Survey of Non-tectonic Surface Displacements Caused by the 2016 Kumamoto Earthquake Around Aso Valley Takayuki NAKANO, Tomokazu KOBAYASHI, Kazuki YOSHIDA and Satoshi FUJIWARA (Published online: 28 December 2016) Abstract The 2016 Kumamoto earthquake caused a large number of surface ruptures (cracks, depressions) and surface earthquake faults emerged along known active faults. Depressional ruptures occurred in Aso Valley. Remarkable horizontal displacements were observed by satellite SAR analysis around the depressional (graben-like) ruptures area. As a result of field surveys around the area, shortening deformations of constructions were identified along the northern edge of the deformed areas as detected by SAR interferogram as well as opening ruptures around the southern edge of the deformed area. It is presumed that these horizontal displacements with depressional ruptures and shortening deformations were caused by lateral movement of shallow underground layers. 1. Preface The 2016 Kumamoto Earthquake caused a large number of surface ruptures such as open cracks, shortening cracks, depressions, and surface earthquake faults emerged along known active faults. These surface ruptures were associated with surface earthquake faults, gravity deformation of slope, ground liquefaction, etc. The Geospatial Information Authority of Japan (GSI) was able to clearly identify these ruptures in the aerial photos taken by GSI in the earthquake disaster area and then analyzed and interpreted these ruptures over a wide area (Nakano et al., 2016; Yoshida et al., 2016). Among these surface ruptures, the surface ruptures that occurred in Aso Valley which is located to the north of Aso Volcano, were distinctive. Also, topical surface displacements were detected by satellite SAR analysis using ALOS-2 data around this area (Fujiwara et al., 2016a). Therefore, field surveys were conducted in the area in order to confirm the actual condition of the surface displacements. Hase, 2011) and flatland (flood plain) extended into Aso Valley. This flood plain was formed to reclaim the Caldera Lake, and marsh zones are distributed in northern part of Aso Valley where reclamation was delayed (GSI, 1994). Also, natural levees and former river beds formed by Kurokawa River which was flowing east to west in Aso Valley are distributed on the flood plain. Alluvial fans formed by erosion debris from central cones of Aso Volcano and taluses formed by depositing of debris from caldera wall are distributed at the periphery of the flood plain (GSI, 1994). Highest seismic intensity during the Kumamoto Earthquake in this area (Uchinomaki, Aso City) was 6- (Japan Meteorological Agency, 2016). 2. Outline of Study Area Figure 1 shows the area of interest in this study. This area is located in the northwest of the Aso Valley which is formed at the north side of Aso Volcano within Aso Caldera. Between ca. 21,000 years ago and ca. 8,000 years ago, there was a lake in Aso Caldera where lake sediments deposited (Ono and Watanabe, 1985; Uchikoshiyama and Fig. 1 Study area (GSI Maps is used as the background map. Shadedrelief map is created by using 10 m DEM data of GSI)

48 Bulletin of the Geospatial Information Authority of Japan, Vol.64 December, 2016 GEONET, which is the continuous GNSS observation network operated by the GSI, recorded approximately 13 cm horizontal displacement in the northwestward direction and vertical subsidence of approximately 8 cm at the Aso site (GSI, 2016a). 3. Observed Phenomena in the Study Area 3.1 Surface Ruptures Identified From Aerial Photos Figure 2 shows the distribution map of surface ruptures identified from aerial photos taken after the 2016 Kumamoto Earthquake in the study area. The background map shows landform classification of the 1:30,000 scale land condition map of volcano Mt. Aso. In this area, clear surface ruptures were identified from the Kurumagaeri area to Uchinomaki area, Aso City. The surface ruptures are arranged linearly or in a meandering pattern, and as a whole are arranged in north-east to southwest direction. Stereoscopic vision interpretation of aerial photos made it possible to confirm that the ground between some major parallel surface ruptures collapsed in graben-like shapes around the Matoishi area and the Shimonohara area, Aso City. Fig. 2 Distribution map of surface ruptures as identified from aerial photos taken in the study area. GSI Maps and the 1:30,000 scale land condition map of volcano Mt. Aso are used as the background map. Also, as the shape of the meandering surface ruptures in the Kurumagaeri area matches the shape of the banked-up land of the former river bed, it is presumed that these meandering surface ruptures were caused by ground subsidence associated with ground liquefaction of former river bed. However, the relationship between location of surface ruptures and landform does not always match in other areas. 3.2 Surface Displacements Detected by SAR Analysis Figure 3 shows the surface ruptures as overlaid on a representative SAR interferogram using ALOS-2 data in the study area. A non-interference zone approximately 1.3km in diameter can be identified in the Matoishi area. Around Shimonohara area, there is a mixture of non-interference zones and interference fringes in the range of approximately 2 km from east-to-west and approximately 1.3 km from south-to-north. The diameter of this interference fringe zone is approximately 2.5 km around the Uchinomaki area. According to Fujiwara et al. (2016a), surface deformations can be observed in these areas having a horizontal displacement of more than 2m in a north-northwest direction. Many surface ruptures are distributed around the southeastern periphery of these surface displacements. Also, GSI (2016b) suggests that these surface displacements have almost no vertical component and displace mostly horizontally toward the caldera wall (Fig. 4). These features are also clarified by the differential analysis between pre- and post-earthquake digital elevation model (DEM) produced by airborne laser survey (JSEG 2016 Kumamoto-Oita Earthquake disaster investigation team, 2016a). In comparison with landform classification of Figure 2, the surface displacements in the Shimonohara area and Uchinomaki area seem to move toward the marsh zone. The SAR interferograms indicate that these surface displacements only occur within certain areas and do not occur in other areas. This phenomenon suggests that largescale surface displacements were located in extremely shallow zones (Fujiwara et al., 2016a).

Field Survey of Non-tectonic Surface Displacements Caused by the 2016 Kumamoto Earthquake Around Aso Valley 49 The field surveys were carried out on May 11 and August 15, 2016. Surveys were carried out at the following places in each area: 1) places of clear graben-like ruptures located along the southern edge of horizontal displacement as detected by SAR interferogram, 2) places located around the northern edge of horizontal displacement as detected by SAR interferogram; although clear ruptures could not be identified from aerial photos. It was expected that open cracks associated with extension field would be found at the points as described at above item 1) and shortening cracks associated with compression field would be found at the point as described at above item 2). Fig. 3 Surface ruptures distribution and representative SAR interferogram by ALOS-2 data (Apr. 15, 2016-Apr. 29, 2016) of study area 5. Findings by Field Survey 5.1 Kurumagaeri Area Figure 5 shows a close up figure of the Kurumagaeri area shown in the Figure 2. In this area, graben-like ruptures were confirmed along former river bed. Major ruptures cross a paddy field surface and a road, which matches the result of aerial photos interpretation (Photo 1), and the portion of the road between two ruptures subsided approximately 30 cm (Photo 2). Fig. 4 Surface displacement map around the study area as revealed by MAI (Multiple Aperture Interferometry) analysis (modified from GSI (2016b)). MAI is a technique for measuring ground displacement along flight direction of the satellite (azimuth offset) using two SAR images. 4. Outline of Field Survey Field surveys were conducted at the distinctive interference fringe areas in order to understand the actual condition of the surface displacements. Our main survey areas were Kurumagaeri, Matoishi and Uchinomaki (Fig. 2, 3). Fig. 5 Distribution map of surface ruptures around Kurumagaeri (The 1:30,000 scale land condition map of volcano Mt. Aso is used as the background map)

50 Bulletin of the Geospatial Information Authority of Japan, Vol.64 December, 2016 Photo 1 Graben-like ruptures along former river bed in Kurumagaeri area (taken towards the northeast on May 11, 2016) cracks was approximately 130 cm. Also, at the Taisho Bridge over the Kurokawa River, the boundary between the bridge and embankment was greatly shortened with approximately 280 cm relatively towards the north (Photo 6). A bulge-like up-lift of the river bed can be identified on the aerial photo taken after the earthquake (Fig. 7). These damages were also reported by other organizations (e.g. JSEG 2016 Kumamoto-Oita Earthquake disaster investigation team (2016b)). Photo 2 Western edge of graben-like ruptures along former river bed in the Kurumagaeri area (taken towards the south on May 11, 2016) Fig. 6 Overlay of ruptures distribution and field surveyed-cracks on the SAR interferogram (Apr. 15, 2016-Apr. 29, 2016) around the Matoishi area 5.2 Matoishi Area Figure 6 shows a close up image of the Matoishi area as shown in the Figure 3. At spots with graben-like ruptures located along the southern edge of the horizontal displacement areas as detected by SAR interferogram, there were open cracks with vertical offsets of up to 70 cm (Photos 3, 4). The total width of open cracks was approximately 95 cm. Around the northern edge of the horizontal displacement areas (Adogase, Aso City), a large number of shortening and open cracks perpendicular to the road were found on the road running north-northwest to south-southeast (Photo 5). The total shortening length of multiple shortening Photo 3 Graben-like rupture in the Matoishi area (taken toward the northwest on May 11, 2016)

Field Survey of Non-tectonic Surface Displacements Caused by the 2016 Kumamoto Earthquake Around Aso Valley 51 Photo 4 Graben-like rupture in the Matoishi area (taken toward the northeast on May 11, 2016) Fig. 7 Bulge-like uplift of Kurokawa River bed identified by aerial photo (taken on Apr. 16, 2016) 5.3 Uchinomaki area Figure 8 shows a close up image of the Uchinomaki area shown in Figure 3. At spots with graben-like ruptures located along the eastern edge of the interference fringe (Ozato/Yuyama, Aso City), a depression was identified between two major ruptures extending in the north-northeast to south-southwest direction (Photo 7). The vertical offset of the ruptures was approximately 50 cm with a side slip of approximately 5 cm. Photo 5 Shortening deformation of the road and lateral gutter (taken toward the northeast on May 11, 2016). The lateral gutter was shortened approximately 20 cm. Photo 6 Shortening deformation of Taisho Bridge (taken toward the west-southwest on Aug. 15, 2016). Embankment was displaced approximately 280 cm relatively towards the north (right side of photo). Fig. 8 Overlay of the surface ruptures and field surveyed-cracks on the SAR interferogram (Apr. 15, 2016-Apr. 29, 2016) around the Uchinomaki area

52 Bulletin of the Geospatial Information Authority of Japan, Vol.64 December, 2016 Around the northern and northwestern edge of horizontal displacement (Ozato/Nishikozono, Aso City), a large number of shortening cracks were identified on the roads and bridge around the junction of Yuura River and Kabaru River and roads around paddy fields located in marsh zones (Photos 8, 9). Shin-kabaru Bridge over the Yuura River was lifted up due to a shortening deformation (Photo 10). Photo 7 Graben-like rupture at the Ozato/Yuyama area (taken toward the northwest on May 11, 2016) Photo 9 Shortening cracks of road and curbstone around paddy fields located in marsh zone (taken toward the east on May 11, 2016) Photo 10 Uplifting of Shin-kabaru Bridge (taken toward the north on May 11, 2016) Photo 8 Shortening cracks at the junction of Yuura River and Kabaru River (taken towards the west on May 11, 2016) 6. Discussions Here, the generating factors of graben-like ruptures and shortening deformations identified around Matoishi and Uchinomaki area are discussed based on the result of field survey and interpretation of SAR interferograms. The SAR interferograms and the result of MAI analysis indicated the block shapes with over 2m horizontal displacement without clear vertical displacement in the Aso Valley (Fujiwara et al., 2016a; GSI, 2016b), and it was indicated that the horizontal displacement occurred in the shallow zone relatively (Fujiwara et al., 2016a).

Field Survey of Non-tectonic Surface Displacements Caused by the 2016 Kumamoto Earthquake Around Aso Valley 53 In the field survey, the graben-like opening ruptures were identified along the southern or eastern edge of the horizontal displacement as detected by SAR analysis located on the extensional field, and a large number of shortening deformations were identified around the northern or northwestern edge of the horizontal displacement located on the compressional field. The amount of horizontal displacements measured on-site coincided with the observation amount by SAR interferometry in order level. These results show that the horizontal displacements were non-tectonic displacements caused by lateral movement of blocks nearby shallow underground layers. In the Uchinomaki area, there are sediment layers on the bottom of former lake near the depth of 50 m (Hase et al, 2010). Also, in some of the thermal wells of Uchinomaki hot springs, breakage of the well was confirmed near the depth of 50 m (Kyushu Environmental Office, 2016). From these facts, it is presumed that the former lake bottom sediments with high water content liquefied with the strong motion of the Kumamoto Earthquake, which became a sliding zone and the block of each area uniformly moved laterally (Fujiwara et al., 2016b). However, the detailed mechanism and driving force are unclear at present and remain open question to future study. 7. Conclusions The field surveys were conducted in and around surface ruptures areas as identified from aerial photos and remarkable surface displacement areas as detected by satellite SAR analysis, and the actual ground conditions of the surface displacements were investigated. As a result, a large number of graben-like open ruptures were identified along the southern or eastern edge and a large number of shortening deformations were identified around the northern or northwestern edge of the horizontal displacements using SAR analysis. Taking into account these results comprehensively, while the formation mechanism is still unclear, it can be speculated that the large horizontal displacements are nontectonic displacements such as lateral movement. Acknowledgements PALSAR-2 data are provided from JAXA (Japan Aerospace Exploration Agency) through joint cooperative agreement between GSI and JAXA and the Earthquake Working Group under a cooperative research contract with JAXA. The ownership of PALSAR-2 data belongs to JAXA. The products of the numerical weather model were provided by JMA (Japan Meteorological Agency) under the agreement between GSI and JMA. References Fujiwara, S., H. Yarai, T. Kobayashi, Y. Morishita, T. Nakano, B. Miyahara, H. Nakai, Y. Miura, H. Ueshiba, Y. Kakiage, and H. Une (2016a): Small-Displacement Linear Surface Ruptures Of The 2016 Kumamoto Earthquake Sequence Detected By ALOS-2 SAR Interferometry, Earth, Planets and Space, 68:160, doi: 10.1186/s40623-016-0534-x. Fujiwara, S., H. Yarai, T. Kobayashi, Y. Morishita, T. Nakano, H. Une, B. Miyahara, H. Nakai, Y. Miura, H. Ueshiba and Y. Kakiage (2016b): Small-displacement Linear Surface Ruptures of the Kumamoto Earthquake Detected by ALOS-2 SAR Interferometry, Journal of Geospatial Information Authority of Japan, 128, 147-162. (in Japanese) Geospatial Information Authority of Japan (1994): 1:30,000 scale land condition map of volcano Mt. Aso. (in Japanese) Geospatial Information Authority of Japan (2016a): Crustal Deformation Associated With The 2016 Kumamoto Earthquake Detected By GEONET, http://www.gsi. go.jp/chibankansi/chikakukansi_kumamoto20160414.ht ml (accessed 5 Oct. 2016). (in Japanese) Geospatial Information Authority of Japan (2016b): Crustal Deformation around the Faults Revealed by MAI Analysis, http://www.gsi.go.jp/cais/topic160428-index -e.html (accessed 5 Oct. 2016). Hase, Y., Y. Miyabuchi, N. Haruta, N. Sasaki and T. Yumoto (2010): Change Of Sedimentary Facies And Topographic Process After The Late Period Of The Last Glacial Age Of Asodani In The Northern Area Of Aso Caldera In Centeral Kyushu, Japan, Bulletin of Goshoura Cretaceous Museum, 11, 1-10. (in Japanese with English abstract) Japan Meteorological Agency (2016): Seismic intensity

54 Bulletin of the Geospatial Information Authority of Japan, Vol.64 December, 2016 database, http://www.data.jma.go.jp/svd/eqdb/data/shin do/event.php?id=9901415&q=7401133 (accessed 5 Oct. 2016). (in Japanese) JSEG 2016 Kumamoto-Oita Earthquake disaster investigation team (2016a): Estimation Of Ground Displacements Around Aso-Caldera Caused By The 2016 Kumamoto Earthquake, From The Geomorphic Image Analysis Of Temporal Lidar Dems, http://www.jseg.or.jp/00- main/pdf/20160822_kumamoto_rep.pdf (accessed 1 Dec. 2016). (in Japanese) JSEG 2016 Kumamoto-Oita Earthquake disaster investigation team (2016b): Horst Graben And Its Surrounding Terrain Displacement That Appeared In The Aso Caldera, http://www.jseg.or.jp/00-main/pdf/ 20160716_Aso.pdf (accessed 1 Dec. 2016). (in Japanese) Kyushu Environmental Office (2016): Survey Results On The Damage Situation To Pasture, Grassland, Springs And Hot Springs In The Aso Area Due To The 2016 Kumamoto Earthquake, http://kyushu.env.go.jp/ pre_20 16/28_1.html (accessed 1 Dec. 2016). (in Japanese) Nakano, T., H. Une, T. Sekiguchi, H. Sakai and K. Yoshida (2016): Distribution Of Surface Cracks Caused By The 2016 Kumamoto Earthquake Interpreted From Aerial Photos, Paper presented at the Japan Geoscience Union meeting 2016, MIS34-P55. (in Japanese with English abstract) Ono, K. and K. Watanabe (1985): Geological Map of Aso Volcano, Geological Survey of Japan. (in Japanese) Uchikoshiyama, U. and Y. Hase (2011): Aquatic Environments Reconstructed Since The Last Glacial Based On Fossil Diatom Assemblages From The Northern Aso Caldera, Central Kyushu, Japan, Earth Science (ChikyuKagaku), 65, 63-79. (in Japanese with English abstract) Yoshida, K., T. Sekiguchi and T. Nakano (2016): Mapping Of Surface Cracks Derived From The Kumamoto Earthquake, Journal of Geospatial Information Authority of Japan, 128, 201-206. (in Japanese)