SAR- and gravity change-based characterization of the distribution pattern of pyroclastic flow deposits at Mt. Merapi during the past 10 years

Mt. Merapi, Indonesia, is one of the most active and dangerous volcanoes in the Torrid Zone. This volcano has erupted frequently and has produced pyroclastic flows following the collapse of the summit lava dome. We used Synthetic Aperture Radar (SAR) data acquired by JERS-1 and RADARSAT-1 satellites from April 1996 to July 2006 to clarify the distribution patterns of the pyroclastic flow deposits. The extent of the deposits, termed P-zones, was accurately extracted by ratio operation and low-level feature extraction from SAR intensity images. These images highlighted temporal changes of the distribution area, perimeter, flow distance, included angle, and collapse direction. To validate the image-processing results, reflectance spectra of the rock samples collected after the eruption in June 2006 were measured in a laboratory. The reflectance spectra of all samples showed similar characteristics to the reference spectra, which were derived from atmospheric correction of Hyperion sensor image data covering the lava dome at the summit. Therefore, P-zones were confirmed to be the pyroclastic flow deposits originating from destruction of the lava dome at the summit. The image-processing results clarified that the extent of the distribution areas, perimeter, flow distances, and included angle of the P-zones were variable among the eruptions, while the collapse direction had a constant pattern. The collapse pattern followed a clockwise change from the south toward the west. By comparing the ratio maps of Bouguer gravity anomaly data in two periods, the change was interpreted to originate from the inclination of the conduit and the formation of shallow and deep magma reservoirs.     

A hypothesis of the Senoumi submarine megaslide in Suruga Bay in Japan—based on the undrained dynamic-loading ring shear tests and computer simulation

The distinctive bathymetric feature exists in the Suruga Bay, Japan. It has been called as Senoumi (Stone flower sea) from old times. Senoumi is a 30?km wide and 20?km long concave feature. Its origin has not been explained yet; however, the feature might be a combined consequence of intensive tectonic activity in the plate border, landslides, and a submarine flow coming from the Oi River. If the Senoumi was caused by a landslide, the latter would be larger than any on-land landslide in Japan. The downshelf “exit” from this feature is much narrower than its central part. This is not usual shape of landslides, but it is similar to the liquefied landslides such as those in quick clays which mobilize great strength reduction after failure. To study Senoumi as a landslide, the shear behaviors of the following three soil samples were investigated by the cyclic and seismic undrained stress control ring shear tests. One sample is volcanic ash taken from the base of landslide deposits (mass transport deposits), from 130 to 190?m deep layer below the submarine floor which was drilled and cored by the Integrated Ocean Drilling Program Expedition 333. Another two samples are the Neogene silty–sand and silt taken from the Omaezaki hill adjacent to the Senoumi, because the shear zone might have been formed in Neogene layers extending from on-land to the continental shelf. The largest strength reduction from peak to steady-state shear resistance in the undrained cyclic loading test was found in volcanic ash. The strength reduction in Neogene silty–sand was smaller than volcanic ash, while the Neogene silt mobilized the least post-failure strength reduction. An integrated model simulating the initiation and motion of earthquake-induced rapid landslides (landslide simulation (LS)-RAPID, Sassa et al. Landslides 7–3:219–236, 2010) was applied to this study. The steady-state shear resistance and other geotechnical parameters measured by the undrained ring shear tests and the greatest strong motion record in the 2011 off-the-Pacific Coast of Tohoku earthquake (M _w 9.0), also known as “2011 Tohoku Earthquake” at the observation point MYG004 (2,933?gal) were input to this model. As the result, it was found that landslides would be triggered by 0.30–1.0 times of MYG004 in volcanic ash, 0.4–1.0 times of MYG004 in Neogene silty–sand and Neogene silt, though the depth and area of triggered landslides were different in soils and intensity of shaking. Feature, created by LS-RAPID using the parameters of volcanic ash, was most similar to the Senoumi in depth and extent. The result obtained from this study includes a hypothesis to be proved, but presents the strong need to investigate the risk of the large-scale submarine landslides which could enhance tsunami wave and possibly enlarge the submarine landslide retrogressively into the adjacent coastal plain by the upcoming mega earthquake in the Nankai Trough.     

Characterization of polychlorinated biphenyls and brominated flame retardants in sediments from riverine and coastal waters of Surabaya, Indonesia

A total of 33 surficial sediments from riverine and coastal waters from Surabaya, Indonesia were analyzed for PCBs and BFRs. Concentrations of PCBs (62 congeners), PBDEs (14 congeners) and HBCDs (3 isomers) varied from <DL-420, <DL-35 and <DL-5.4 ng g⁻¹ dw, respectively. Higher concentrations of these compounds were found in riverine than coastal sediments. Their levels and distribution were influenced by proximity to the point sources and TOC. The predominant congeners were CB-153, -28, -138, -149, -180, -33 and BDE-209, -207, -206, -197, -196, -183, -99, -47 for PCBs and PBDEs, respectively, and γ-isomer for HBCDs. Debromination of BDE-209 might be taking place producing lower toxic congeners in sediment. Levels of PCBs in riverine sediments were comparable with some polluted areas worldwide, but PBDEs and HBCDs were lower. Hazard assessment of PCBs indicated possible toxic potential, particularly in areas close to point sources.     

Formation process of two massive dams following rainfall-induced deep-seated rapid landslide failures in the Kii Peninsula of Japan

Extreme heavy rainfall due to Typhoon Talas on September 2–4, 2011 in the Kii Peninsula, Japan, triggered numerous floods and landslides. This study investigates the mechanism and the entire process of rainfall-induced deep-seated landslides forming two massive dams in the Kuridaira and Akatani valleys, respectively. The mechanism of the rapid deep-seated landslides is examined through a series of laboratory experiments on samples from sliding surfaces by using undrained high-stress dynamic-loading ring-shear apparatus. The test results indicate that the failure of samples is triggered by excess pore water pressure generation under a shear displacement from 2 to 7 mm with a pore pressure ratio ranging from 0.33 to 0.37. The rapid movement of landslides is mainly attributed to high mobility due to the liquefaction behavior of both sandstone-rich and shale samples. Geomorphic settings and landslide mobility are major contributing factors to the dam formation. Additionally, shear displacement control tests show that a certain amount of shear displacement between 2 and 7 mm along the sliding surfaces of the gravitationally deformed slopes might have led to the failures. Importantly, computer simulation with LS-RAPID software using input parameters obtained from physical experiments is employed to interpret the entire formation process of the abovementioned two landslide dams. The simulation results are examined in accordance with the observed on-site geomorphic features and recorded data to explain the possibility of sliding processes. The results further point out that local failures are initiated from the lower middle part of the landslide bodies where the geological boundary exists. This condition most probably influences the landslide initiation in the two case studies. This research is therefore helpful for hazard assessment of slopes that are susceptible to deep-seated landslides and other sequential processes in areas with geology and geomorphology similar to that of the Kii Peninsula.     

Reconstruction of Tsunami Inland Propagation on December 26, 2004 in Banda Aceh, Indonesia, through Field Investigations

This paper presents the results from an extensive field data collection effort following the December 26, 2004 earthquake and tsunami in Banda Aceh, Sumatra. The data were collected under the auspices of TSUNARISQUE, a joint French-Indonesian program dedicated to tsunami research and hazard mitigation, which has been active since before the 2004 event. In total, data from three months of field investigations are presented, which detail important aspects of the tsunami inundation dynamics in Banda Aceh. These include measurements of runup, tsunami wave heights, flow depths, flow directions, event chronology and building damage patterns. The result is a series of detailed inundation maps of the northern and western coasts of Sumatra including Banda Aceh and Lhok Nga. Among the more important findings, we obtained consistent accounts that approximately ten separate waves affected the region after the earthquake; this indicates a high-frequency component of the tsunami wave energy in the extreme near-field. The largest tsunami wave heights were on the order of 35 m with a maximum runup height of 51 m. This value is the highest runup value measured in human history for a seismically generated tsunami. In addition, our field investigations show a significant discontinuity in the tsunami wave heights and flow depths along a line approximately 3 km inland, which the authors interpret to be the location of the collapse of the main tsunami bore caused by sudden energy dissipation. The propagating bore looked like a breaking wave from the landward side although it has distinct characteristics. Patterns of building damage are related to the location of the propagating bore with overall less damage to buildings beyond the line where the bore collapsed. This data set was built to be of use to the tsunami community for the purposes of calibrating and improving existing tsunami inundation models, especially in the analysis of extreme near-field events.     

Assessment of the effects of vegetation on soil erosion risk by water: a case of study of the Batta watershed in Tunisia

Soil erosion by water is a serious environmental problem which affects particularly the agriculture of developing countries. Due to specific factors, such as high rainfall intensity, steep slopes and vegetation scarcity, Tunisia is prone to soil erosion. Taking this into account, the main objective of this study was to estimate the soil erosion risk in the Batta watershed in Tunisia using qualitative and quantitative erosion model with remote sensing data and geographic information system (GIS). Moreover, a developed method that deals with evaluating the impact of vegetation on soil erosion by water is also applied. This method used multi-temporal satellite images with seasonal variability and the transformed soil adjusted vegetation index (TSAVI) which is appropriate in arid and semi-arid areas. For both erosion models, the results show that a large area of the Batta watershed is seriously affected by erosion. This potentially high risk is due especially to severe slopes, poor vegetation coverage and high soil erodibility in this watershed. Furthermore, the use of multi-temporal satellite images and vegetation index show that the effect of vegetation is a significant factor to protect the soil against erosion. The risk is especially serious in the summer season, but it decreases with the growth of vegetation cover in spring. Erosion model, combined with a GIS and remote sensing, is an adequate method to evaluate the soil erosion risk by water. The findings can be used by decision makers as a guideline to plan appropriate strategies for soil and water conservation practices.     

Groundwater exploration using integrated geophysics method in hard rock terrains in Mount Betung Western Bandar Lampung,Indonesia

The presence of hard rock in Mount Betung has caused the misalignment of the groundwater aquifers,and resulted in many drilling failures for groundwater.An integrated geophysics method using gravity survey and Geoelectric Vertical Electrical Soundings(VES)were conducted to study the effect of basement and hard rock on groundwater prospects.From the gravity method,38 mapping points were carried out randomly,with a distance of 1-2 km in-between.Meanwhile,from the geoelectric method,51 VES points were acquired at the foot of Mount Betung.The acquisition was conducted with a Schlumberger configuration with AB/2=1 m to 250 m.The results show the Bouguer Anomaly in the west is 50-68 mgal due to the presence of hard rock in Mount Betung.This anomaly responds to relatively shallow hard rocks near surface.Hard rocks composed of andesite and breccia normally present at the depth of 5-180 m during well construction.Resistivity isopach mapping from VES data(at AB/2=50 m,100 m,and 150 m)shows the dominant constituents of hard rock.Fractures in hard rock contribute to secondary porosity,which could be a prospect zone that transmit groundwater.This finding shows that the fractures are randomly scattered,causing several well failures that have been worked.Furthermore,the fractures in the hard rock at the foot of Mount Betung acts as conduits between recharge at Mount Betung and the aquifer in the Bandar Lampung Basin.     

Observations of intraseasonal variability in the Sunda Strait throughflow

Using velocity profiles observed by bottom-mounted ADCPs, we identified strong intraseasonal variability in the Sunda Strait throughflow. This intraseasonal variability, with typical periods of 20–40 days and the strongest energy occurring in the boreal spring, can reverse the Sunda Strait throughflow. Further analysis showed this intraseasonal variability to be closely related to local zonal wind and the sea level gradient along the strait. These observations confirm for the first time the existence of Kelvin-wave-like signals in the Sunda Strait, propagating from the equatorial Indian Ocean. This study also provides new insights into the effects of Kelvin waves on the Sunda Strait throughflow.     

Hourly sea surface temperature retrieval using the Japanese Geostationary Satellite,Multi-Functional Transport Satellite (MTSAT)

An algorithm has been developed for retrieving sea surface temperature (SST) from hourly data transmitted from the Japanese Advanced Meteorological Imager (JAMI) aboard a Japanese geostationary satellite, Multi-functional Transport Satellite (MTSAT)-1R. Threshold tests screening cloudy pixels are empirically adjusted to cases of daytime with/without sun glitter, and nighttime. The Non-Linear SST (NLSST) equation, including several new additional terms, is used to calculate MTSAT SST. The estimated SST is compared with drifting and moored buoy measurements, with the result that the bias of the MTSAT SST is nearly 0.0°K. The root mean square (rms) error is about 0.8°K, and it is 0.7°K under the condition that the satellite zenith angle is less than 50°. It is demonstrated that the hourly MTSAT SST produced by the algorithm developed here captures diurnal SST variations in the equatorial sea in mid-November 2006.