Testing a model for predicting the timing and location of shallow landslide initiation in soil‐mantled landscapes

The growing availability of digital topographic data and the increased reliability of precipitation forecasts invite modelling efforts to predict the timing and location of shallow landslides in hilly and mountainous areas in order to reduce risk to an ever‐expanding human population. Here, we exploit a rare data set to develop and test such a model. In a 1·7 km² catchment a near‐annual aerial photographic coverage records just three single storm events over a 45 year period that produced multiple landslides. Such data enable us to test model performance by running the entire rainfall time series and determine whether just those three storms are correctly detected. To do this, we link a dynamic and spatially distributed shallow subsurface runoff model (similar to TOPMODEL) to an in?nite slope model to predict the spatial distribution of shallow landsliding. The spatial distribution of soil depth, a strong control on local landsliding, is predicted from a process‐based model. Because of its common availability, daily rainfall data were used to drive the model. Topographic data were derived from digitized 1 : 24 000 US Geological Survey contour maps. Analysis of the landslides shows that 97 occurred in 1955, 37 in 1982 and ?ve in 1998, although the heaviest rainfall was in 1982. Furthermore, intensity–duration analysis of available daily and hourly rainfall from the closest raingauges does not discriminate those three storms from others that did not generate failures. We explore the question of whether a mechanistic modelling approach is better able to identify landslide‐producing storms. Landslide and soil production parameters were ?xed from studies elsewhere. Four hydrologic parameters characterizing the saturated hydraulic conductivity of the soil and underlying bedrock and its decline with depth were ?rst calibrated on the 1955 landslide record. Success was characterized as the most number of actual landslides predicted with the least amount of total area predicted to be unstable. Because landslide area was consistently overpredicted, a threshold catchment area of predicted slope instability was used to de?ne whether a rainstorm was a signi?cant landslide producer. Many combinations of the four hydrological parameters performed equally well for the 1955 event, but only one combination successfully identi?ed the 1982 storm as the only landslide‐producing storm during the period 1980–86. Application of this parameter combination to the entire 45 year record successfully identi?ed the three events, but also predicted that two other landslide‐producing events should have occurred. This performance is signi?cantly better than the empirical intensity–duration threshold approach, but requires considerable calibration effort. Overprediction of instability, both for storms that produced landslides and for non‐producing storms, appears to arise from at least four causes: (1) coarse rainfall data time scale and inability to document short rainfall bursts and predict pressure wave response; (2) absence of local rainfall data; (3) legacy effect of previous landslides; and (4) inaccurate topographic and soil property data. Greater resolution of spatial and rainfall data, as well as topographic data, coupled with systematic documentation of landslides to create time series to test models, should lead to signi?cant improvements in shallow landslides forecasting. Copyright © 2003 John Wiley & Sons, Ltd.     

Accretion and exhumation at a Variscan active margin, recorded in the Saxothuringian flysch

The Saxothuringian flysch basin, on the north flank of the Central European Variscides, was fed and eventually overthrust by the northwestern, active margin of the Tepla-Barrandian terrane. Clast spectra, mineral composition and isotopic ages of detrital mica and zircon have been analyzed in order to constrain accretion and exhumation of rocks in the orogenic wedge. The earliest clastic sediments preserved are of early Famennian age (ca. 370?Ma). They are exposed immediately to the NW of the suture, and belong to the par-autochthon of the foreland. Besides ultramafic (?ophiolite) material, these rocks contain clasts derived from Early Paleozoic continental slope sediments, originally deposited at the NW margin of the Saxothuringian basin. These findings, together with the paleogeographic position of the Famennian clastics debris on the northwestern passive margin, indicate that the Saxothuringian narrow ocean had been closed by that time. Microprobe analyses of detrital hornblendes suggest derivation from the “Randamphibolit” unit, now present in the middle part of the Saxothuringian allochthon (Münchberg nappes). Detrital zircons of metamorphic rocks formed a little earlier (ca. 380?Ma) indicate rapid recycling at the tectonic front. The middle part of the flysch sequence (ca. early to middle Viséan), both in the par-autochthon and in the allochthon, contains abundant clasts of Paleozoic rocks derived from the northwestern slope and rise, together with debris of Cadomian basement, 500-Ma granitoids and 380?Ma (early Variscan) crystalline rocks. All of these source rocks were still available in the youngest part of the flysch (c. middle to late Viséan), but some clasts record, in addition, accretion of the northwestern shelf. Our findings permit deduction of minimum rates of tectonic shortening well in excess of 10–30?mm per year, and rates of exhumation of ca. 3?mm/a, and possibly more.     

Ground deformation of Nisyros Volcano (Greece) for the period 1995–2002: Results from DInSAR and DGPS observations

Differential GPS (DGPS) and Differential Interferometric Synthetic Aperture Radar (DInSAR) analyses were applied to the Kos-Yali-Nisyros Volcanic Field (SE Hellenic Volcanic Arc) to quantify the ground deformation of Nisyros Volcano. After intense seismic activity in 1996, a GPS network was installed in June 1997 and re-occupied annually up to 2002. A general uplift ranging from 14 to 140 mm was determined at all stations of the network. The corresponding horizontal displacements ranged from 13 to 53 mm. The displacement vectors indicate that the island is undergoing extension towards the East, West and South. A two-source “Mogi” model combined with assumed motion along the Mandraki Fault was constructed to fit the observed deformation. The best-fit model assumes sources at a depth of 5500 m NW of the centre of the island and at 6500 m offshore ESE of Yali Island. DInSAR analysis using four pairs of images taken between May 1995 and September 2000 suggests that deformation was occurring during 1995 before the start of the seismic crisis. An amplitude of at least 56 mm along the slant range appeared for the period 1996 through 1999. This deformation is consistent with the two-source model invoked in DGPS modelling. Surface evidence of ground deformation is expressed in the contemporaneous reactivation of the Mandraki Fault. In addition, a 600 m long N-S trending irregular rupture in the caldera floor was formed between 2001 and 2002. This rupture is interpreted as the release of surface stress in the consolidated epiclastic and hydrothermal sediments of the caldera floor.     

Gondwana breakup: no evidence for a Davie Fracture Zone offshore northern Mozambique,Tanzania and Kenya

Plate tectonic reconstructions assume a major inactive transform fault, the Davie Fracture Zone, in the West Somali Basin, along which Madagascar is thought to have migrated southwards following Gondwana breakup in the Mesozoic. Based on the interpretation of reflection seismic data, we show that the Walu Ridge offshore Kenya and the Kerimbas Basin offshore northern Mozambique are tectonically unrelated to the southward motion of Madagascar and correlate with Late Cretaceous volcanism and inversion in Kenya and the evolution of the East African Rift System respectively. Offshore Tanzania, geophysical data do not show basement structures indicating the presence of a major transform fault. These results challenge the commonly supported transform margin concept and imply a more southerly pre‐breakup position of Madagascar within Gondwana. Opening of the West Somali Basin by SW‐propagating oblique rifting and seafloor spreading is proposed.     

Assessment of small-diameter shallow wells for managed aquifer recharge at a site in southern Styria,Austria

An approach to establish the recharge component of managed aquifer recharge (MAR) has recently been proposed that uses small-diameter shallow wells installed using relatively inexpensive drilling methods such as direct push. As part of further development of that approach, a generalized procedure is presented for a technical and economic assessment of the approach’s potential in comparison to other systems. Following this procedure, the use of small-diameter wells was evaluated both experimentally and numerically for a site located in southern Styria, Austria. MAR is currently done at the site using a horizontal pipe infiltration system, and system expansion has been proposed with a target rate of 12 l/s using small-diameter wells as one possible option. A short-duration single-well field recharge experiment (recharge rate 1.3–3.5 l/s) was performed (recharge by gravity only). Numerical modeling of the injection test was used to estimate hydraulic conductivity (K). Quasi-steady-state, single-well recharge simulations for different locations, as well as a long-term transient simulation, were performed using the K value calibrated from the field injection test. Results indicate that a recharge capacity of 4.1 l/s was achievable with a maximum head rise of 0.2 m at the injection well. Finally, simulations were performed for three different well fields (4, 6 and 8 wells, respectively) designed to infiltrate a target rate of 12 l/s. The experimental and numerical assessments, supported by a cost analysis of the small-diameter wells, indicate that the small-diameter wells are a viable, cost-effective recharge approach at this and other similar sites.     

2D probabilistic prediction of sparsely measured earth properties constrained by geophysical imaging fully accounting for tomographic reconstruction ambiguity

Many hydrological, environmental, or engineering exploration tasks require predicting spatially continuous scenarios of sparsely measured borehole logging data. We present a methodology to probabilistically predict such scenarios constrained by ill-posed geophysical tomography. Our approach allows for transducing tomographic reconstruction ambiguity into the probabilistic prediction of spatially continuous target parameter scenarios. It is even applicable to data sets where petrophysical relations in the survey area are non-unique, i.e., different facies related petrophysical relations may be present. We employ static two-layer artificial neural networks (ANNs) for prediction and additionally evaluate, whether the training performance of the ANNs can be used to rank geophysical tomograms, which are mathematically equal reconstructions of physical parameter distributions in the ground. We illustrate our methodology using a realistic synthetic database for maximal control about the prediction performance and ranking potential of the approach. For doing so, we try to link geophysical radar and seismic tomography as input parameters to porosity of the ground as target parameter of ANN. However, the approach is flexible and can cope with any combination of geophysical tomograms and hydrologic, environmental or engineering target parameters. Ranking of equivalent geophysical tomograms based on additional borehole logging data is found to be generally possible, but risks remain that the ranking based on the ANN training performance does not fully coincide with the closeness of geophysical tomograms to ground truth. Since geophysical field data sets do usually not offer control options similar to those used in our synthetic database, we do not recommend the utilization of recurrent ANNs to learn weights for the individual geophysical tomograms used in the prediction procedure.     

Seismic images of a collision zone offshore NW Sabah/Borneo

Multichannel reflection seismic data from the southern South China Sea, refraction and gravity modelling were used to investigate the compressional sedimentary structures of the collision-prone continental margin off NW Borneo. An elongated imbricate deepwater fan, the toe Thrust Zone bounds the Northwest Borneo Trough to the southeast. The faults separating the individual imbricates cut through post-Early Miocene sediments and curve down to a carbonate platform at the top of the subsiding continental Dangerous Grounds platform that forms the major detachment surface. The age of deformation migrates outward toward the front of the wedge. We propose crustal shortening mechanisms as the main reason for the formation of the imbricate fan. At the location of the in the past defined Lower Tertiary Thrust Sheet tectonostratigraphic province a high velocity body was found but with a much smaller extend than the previously defined structure. The high velocity structure may be interpreted either as carbonates that limit the transfer of seismic energy into the sedimentary layers beneath or as Paleogene Crocker sediments dissected by remnants of a proto-South China Sea oceanic crust that were overthrust onto a southward migrating attenuated continental block of the Dangerous Grounds during plate convergence.     

Trends,drivers and impacts of changes in swidden cultivation in tropical forest-agriculture frontiers: A global assessment

This meta-analysis of land-cover transformations of the past 10–15 years in tropical forest-agriculture frontiers world-wide shows that swidden agriculture decreases in landscapes with access to local, national and international markets that encourage cattle production and cash cropping, including biofuels. Conservation policies and practices also accelerate changes in swidden by restricting forest clearing and encouraging commercial agriculture. However, swidden remains important in many frontier areas where farmers have unequal or insecure access to investment and market opportunities, or where multi-functionality of land uses has been preserved as a strategy to adapt to current ecological, economic and political circumstances. In some areas swidden remains important simply because intensification is not a viable choice, for example when population densities and/or food market demands are low. The transformation of swidden landscapes into more intensive land uses has generally increased household incomes, but has also led to negative effects on the social and human capital of local communities to varying degrees. From an environmental perspective, the transition from swidden to other land uses often contributes to permanent deforestation, loss of biodiversity, increased weed pressure, declines in soil fertility, and accelerated soil erosion. Our prognosis is that, despite the global trend towards land use intensification, in many areas swidden will remain part of rural landscapes as the safety component of diversified systems, particularly in response to risks and uncertainties associated with more intensive land use systems.