GIS-based weights-of-evidence modelling of rainfall-induced landslides in small catchments for landslide susceptibility mapping

Landslide susceptibility mapping is a vital tool for disaster management and planning development activities in mountainous terrains of tropical and subtropical environments. In this paper, the weights-of-evidence modelling was applied, within a geographical information system (GIS), to derive landslide susceptibility map of two small catchments of Shikoku, Japan. The objective of this paper is to evaluate the importance of weights-of-evidence modelling in the generation of landslide susceptibility maps in relatively small catchments having an area less than 4 sq km. For the study area in Moriyuki and Monnyu catchments, northeast Shikoku Island in west Japan, a data set was generated at scale 1:5,000. Relevant thematic maps representing various factors (e.g. slope, aspect, relief, flow accumulation, soil depth, soil type, land use and distance to road) that are related to landslide activity were generated using field data and GIS techniques. Both catchments have homogeneous geology and only consist of Cretaceous granitic rock. Thus, bedrock geology was not considered in data layering during GIS analysis. Success rates were also estimated to evaluate the accuracy of landslide susceptibility maps and the weights-of-evidence modelling was found useful in landslide susceptibility mapping of small catchments.     

Spatio-temporal evolution of Coulomb stress in the Pacific slab inverted from the seismicity rate change and its tectonic interpretation in Hokkaido, Northern Japan

We estimated spatio-temporal evolution of Coulomb stress within the subducted Pacific slab in Hokkaido from the analysis of seismicity rate change. For this purpose we used earthquake catalog from the Institute of Seismology and Volcanology (ISV), Hokkaido University for the period 1993/4/1–2006/12/31 after relocating to compensate location errors due to the heterogeneous P- and S-wave structure beneath Hokkaido. We found that spatial pattern of Coulomb stress change inverted from the seismicity rate change is comparable with static change in Coulomb stress estimated from dislocation models. Our results and analyses reveal important insights on spatio-temporal pattern of deformation of the subducted Pacific slab in terms of Coulomb stress change. We found that the 2003 Tokachi Oki earthquake (Mw = 8.0) pervasively perturbed Coulomb stress in a regional scale with a significant impact to trigger the 2004 Kushiro Oki earthquake. The 2004 Kushiro Oki earthquake (Mw = 7.0) is another significant stressing event that changed the pattern of Coulomb stress in the area. We found that stressing events with magnitude smaller than 7.0 has minimal impact on Coulomb stress change in the Pacific slab. Similarly, comparatively deep focused large earthquakes could not change Coulomb stress significantly. Further the pattern of Coulomb stress change after the 2003 Tokachi Oki earthquake correlates the pattern of afterslip distribution in Hokkaido.     

Geochemistry of Mesozoic intracontinental basalts from Yunnan, southern China: Implications for geochemical evolution of the subcontinental lithosphere

Summary Southwestern Yunnan, comprising the Yangtze and Shan-Thai microcontinents and the Simao block, has successively undergone subduction of an oceanic plate, followed by a collision of the microcontinents and intracontinental rifting associated with basaltic volcanism during Late Paleozoic to Mesozoic.The Triassic Nanjian basalts, erupted on the Yangtze microcontinent, have more enriched isotopic ratios and higher LREE/HFSE and LREE/HREE ratios. This suggests the existence of an enriched subcontinental lithosphere under the Yangtze microcontinent which stabilized over long periods of the earth's history (> 2Ga).The Middle Jurassic Simao basalts have more depleted geochemical features and also have element enrichments characteristic of a subduction zone environment, although the basalts were erupted in an intracontinental graben. It may be inferred that the lithospheric mantle of the Simao block was modified by subduction processes during Latest Carboniferous to Late Triassic prior to the onset of the Middle Jurassic continental rifting. The lack of correlation between depletion of HFSE, Y and HREE, and relative enriched Nd isotopic ratios suggests that the source depletion of the Simao basalts is not an old feature and has been contemporaneous with the subduction-related enrichment through mantle metasomatism shortly before the basalts were produced.The Middle Jurassic Baoshan basalts which erupted during the continental rifting on the Shan-Thai microcontinent have an Sr-Nd isotopic composition similar to the bulk earth and higher concentrations of incompatible trace elements. These features suggest that the subcontinental lithosphere under the Shan-Thai microcontinent underwent mantle metasomatism just prior to eruption of the Baoshan basalt.

---

Geochemie mesozoischer interkontinentaler Basalte aus Jünnan, Südchina: Hinweise auf die geochemische Entwicklung der subkontinentalen Lithosphäre

Zusammenfassung Südwestjünnan umfaßt den Jangtse und den Shan-Thai Mikrokontinent und den Simao Block. Das Gebiet wurde von aufeinander folgenden Subduktionsphasen einer ozeanischen Platte betroffen, auf die Kollision der Mikrokontinente und interkontinentales Rifting folgte. Dieses war mit basaltischem Vulkanismus während des späten Paläozoikums bis ins Mesozoikum assoziiert.Die triassischen Nanjian-Basalte, die auf dem Jangtse Mikrokontinent eruptierten, haben mehr angereicherte Isotopenverhältnisse und höhere LREE/HFSE und LREE/HREE Verhältnisse. Dieses weist auf eine angereicherte subkontinentale Lithosphäre unter dem Jangtse Mikrokontinent hin, die sich während langer Perioden der Erdgeschichte stabilisierte (>2Ga).Die mittel jurassischen Simao-Basalte haben eine mehr verarmte geochemische Signatur aber auch Elementanreicherungen, die für ein Subduktionszonen-Milieu charakteristisch sind, obwohl die Basalte in einem interkontinentalen Graben ausgetreten sind. Man kann daraus schließen, daß der lithosphärische Mantel des Simao-Blockes durch Subduktionsprozesse während des jüngsten Karbons bis in die späte Trias vor dem Beginn des mittel-jurassischen kontinentalen Riftings modifiziert worden war. Das Fehlen einer Korrelation zwischen der Anreicherung von HFSE, Y und HREE und relativ angereicherter Nd-Isotopenverhältnisse weist darauf hin, daß die Verarmung der Quelle der Simaobasalte nicht weit zurückreicht. Sie dürfte viel eher gleichaltrig mit der subduktions-bedingten Anreicherung durch Mantel-Metasomatose kurz vor der Entstehung der Basalte sein.Die mittel-jurassischen Baoshan-Basalte, die während des kontinentalen Riftings auf den Shan-Thai Mikrokontintent eruptierten, haben eine Sr-Nd-Isotopensignatur, die ähnlich der Gesamterde ist, jedoch höhere Konzentrationen inkompatibler Spurenelemente zeigt. All dies legt nahe, daß die subkontinentale Lithosphäre unter dem Shan-Thai-Mikrokontinent kurz vor der Eruption der Baoshan-Basalte von Mantel-Metasomatose betroffen worden ist.

---

---

With 8 Figures

---

Visiting Fellow, Geology Department, Australian National University, Australia     

Volcanism in the Sumisu Rift, I. Major element, volatile, and stable isotope geochemistry

A bimodal volcanic suite with KAr ages of 0.05–1.40 Ma was collected from the Sumisu Rift using alvin. These rocks are contemporaneous with island arc tholeiite lavas of the Izu-Ogasawara arc 20 km to the east, and provide a present day example of volcanism associated with arc rifting and back-arc basin initiation. Major element geochemistry of the basalts is most similar to that of basalts found in other, more mature back-arc basins, which indicates that back-arc basins need not begin their magmatic evolution with lavas bearing strong arc signatures.Volatile concentrations distinguish Sumisu Rift basalts from island arc basalts and MORB. H₂O contents, which are at least four times greater than in MORB, suppress plagioclase crystallization. This suppression results in a more mafic fractionating assemblage, which prevents Al₂O₃ depletion and delays the initiation of Fe₂O_3^(tot) and TiO₂ enrichment. However, unlike arc basalts,Fe^3+/ΣFe ratios are only slightly higher than in MORB and are insufficient to cause magnetite saturation early enough to suppress Fe₂O_3^(tot) and TiO₂ enrichment. Thus, major element trends are more similar to those of MORB than arcs.H₂O, CO₂ and S are undersaturated relative to pure phase solubility curves, indicating exsolution of an H₂O-rich mixed gas phase. HighH₂O/S, highδD, and low (MORB-like)δ³⁴S ratios are considered primary and distinctive of the back-arc basin setting.     

Seasonal changes in the mesozooplankton biomass and community structure in subarctic and subtropical time-series stations in the western North Pacific

Seasonal changes in mesozooplankton biomass and their community structures were observed at time-series stations K2 (subarctic) and S1 (subtropical) in the western North Pacific Ocean. At K2, the maximum biomass was observed during the spring when primary productivity was still low. The annual mean biomasses in the euphotic and 200- to 1000-m layers were 1.39 (day) and 2.49 (night) g C m^?2 and 4.00 (day) and 3.63 (night) g C m^?2, respectively. Mesozooplankton vertical distribution was bimodal and mesopelagic peak was observed in a 200- to 300-m layer; it mainly comprised dormant copepods. Copepods predominated in most sampling layers, but euphausiids were dominant at the surface during the night. At S1, the maximum biomass was observed during the spring and the peak timing of biomass followed those of chlorophyll a and primary productivity. The annual mean biomasses in the euphotic and 200- to 1000-m layers were 0.10 (day) and 0.21 (night) g C m^?2 and 0.47 (day) and 0.26 (night) g C m^?2, respectively. Copepods were dominant in most sampling layers, but their mean proportion was lower than that in K2. Mesozooplankton community characteristics at both sites were compared with those at other time-series stations in the North Pacific and with each other. The annual mean primary productivities and sinking POC fluxes were equivalent at both sites; however, mesozooplankton biomasses were higher at K2 than at S1. The difference of biomasses was probably caused by differences of individual carbon losses, population turnover rates, and trophic structures of communities between the two sites.     

Effects of a continental slope along the western boundary on the abyssal circulation

To investigate effects of a continental slope along the western boundary on the abyssal circulation, numerical experiments using multi-level models were carried out. An ocean which extends over the northern and southern hemispheres is forced by cooling inside the ocean at the southwest corner of the basin and uniform heating through the sea surface. When the reference density for the cooling is vertically uniform, effects of the slope emerge clearly for the slope with considerably broad width. The deep western boundary current flowing over the slope feeds no bottom flows in the southern hemisphere, and carries the warmed deep water into the northern hemisphere. This leads to the increased meridional density gradient, which results in the modification of deep flow patterns. When the reference density is vertically distributed, the upper and lower northward flowing western boundary currents form in the deep layer. As the density stratification relaxes the topographic control, the westward intensification of the upper boundary current is achieved over the slope. The intensified flow is accompanied by the countercurrent and they form the horizontal recirculation over the slope. However, the effects are confined around the slope region and the interior flow patterns do not change. The lower boundary current is not significantly affected by the slope and has the large width with no countercurrent. It is found that the actual continental slope does not have significant effects on the gross feature of the thermohaline circulation.     

Global distribution and atmospheric transport of chlorinated hydrocarbons: HCH (BHC) isomers and DDT compounds in the Western Pacific,Eastern Indian and Antarctic Oceans

Concentrations of chlorinated hydrocarbons such as HCH isomers and DDT compounds were determined in air and surface water samples taken from the Western Pacific, Eastern Indian and Antarctic Oceans. The most interesting finding was their presence in measurable concentrations in the Antarctic Ocean. Chlorinated hydrocarbon pesticides are widely distributed in the open ocean environment over both the northern and southern hemispheres, and some characteristic distribution patterns of pesticide species in different oceanic regions were observed both in air and water samples. HCH residues found in the northern hemisphere were much higher in concentration than those in the southern hemisphere. On the other hand, higher concentrations of DDT residues were found in the tropical regions, but their levels were not so different between both the northern and southern hemispheres. HCH isomers found in the northern hemisphere had the following order of concentrations:-HCH> HCH>-HCH, while in the southern hemisphere-HCH was apparently dominant. DDT compound compositions were rather uniform in all the oceans surveyed, and more than 50% wasp,p-DDT. These facts can be explained by the world wide situation regarding pesticide use and the physicochemical properties of the pesticides such as their vapor pressures and water solubilities. In addition, the meridional circulation of the atmosphere, particularly the mass flows of the Hadley and Ferrel cells in the troposphere, also contributes to the atmospheric transport and global distribution of these pesticides.     

The continuous limit of the multiple lens effect and the optical scalar equation

We study the continuous limit of the multiple gravitational lensing theory based on the thin lens approximation. We define a new, light-path dependent angular diameter distance     and show that it satisfies the optical scalar equation. The distance provides relations between quantities used in gravitational lensing theory (the convergence, the shear and the twist terms) and those used in scalar optics theory (the rates of expansion, shear and rotation).     

Phase and group velocities and Q of mantle Love and Rayleigh waves of the first two modes and their azimuthal dependences for the 1963 Kurile Islands earthquake

Phase and group velocities and Q of mantle Love and Rayleigh waves from the 1963 Kurile Islands earthquake (M_w = 8.5) were determined over 37 great circle paths by a time variable filtering technique, in a period range 100–500 s for the fundamental modes and 100–275 s for the first higher modes. The preliminary reference Earth model (PREM) explains reasonably well the average dispersion results for the fundamental Love and Rayleigh waves. There exists a small, but significant inconsistency between the observation and the model for the first higher Love and Rayleigh waves. The Q structure of PREM is inconsistent with the observation for the fundamental Love waves, but explains other observations reasonably well. The dispersion of each mode shows a clear azimuthal dependence from which the four azimuthal windows were established. The phase and group velocity measurements for each window were, in general, shown to be mutually consistent. The azimuthal variations are largest for the first higher Rayleigh waves, indicating strong lateral heterogeneity in the structure of the low velocity zone. The first of the four windows is characterized by the largest fraction of Precambrian shields and the second window by the largest fraction of normal oceans. A comparison of these two windows may give some insight into deep lateral heterogeneity between continents and oceans. The observed phase and group velocities of the first window are systematically higher than those of the second window for the fundamental Love and Rayleigh waves at periods up to 400 s, and for the first higher Love and Rayleigh waves up to 175 s. Their differences are greatest for the first higher Rayleigh waves and least for the fundamental Rayleigh waves. Although the fundamental Rayleigh waves show the least velocity differences, their persistence up to a period of longer than 300 s is in striking contrast with some of the pure path phase velocities derived earlier for continents and oceans. A set of models for continents and oceans. PEM-C and PEM-O are not consistent with our observation. The third azimuthal window is characterized by trench-marginal seas and the fourth window by mountainous areas, typically the Asian high plateaus from northern China to the Middle East through Tibet. A comparison of these two windows gives some information about deep structural differences between subduction zones and continental collision zones, both belonging to plate convergence zones. For the fundamental and the first higher Love waves, the phase and group velocities for the third window are markedly low, whereas those for the fourth window are somewhat comparable to those for the second window. Slow Rayleigh waves are evident for two windows, with the fourth window apparently being the slowest for the fundamental Rayleigh above 200 s and for the first higher Rayleigh. For the fundamental Rayleigh waves, the third window is very slow below 200 s, but becomes progressively fast as the period increases and tends to be the fastest window around 400 s, suggesting a deep seated high velocity anomaly beneath trench-marginal seas. The dispersion characteristics of the fourth window indicate a thick high velocity lid with an extensive low velocity zone beneath it. The shield-like lithosphere, coupled with an extensive low velocity zone, may be a characteristic feature of continental collision zones. The particle motion of the fundamental Love waves was found not to be purely transverse to a great-circle connecting the epicenter to a station. The departure from the purely transverse motion is systematic among different periods, different G arrivals (G₂, G₃,…) and different stations, which may be interpreted as being due to lateral refraction.