Large-scale mass wasting in ancient volcanic materials

Giant landslides are significant hazards associated with many active volcanic edifices. We describe a similar feature on ancient (>4 Ma) volcanic deposits subject to active tectonism. The landslide is approximately 3 km long by 1 km wide, with an estimated depth of 400 m. Side margins are straight and parallel, mimicking regional structure; narrow valleys incised down these margins provide low-strength side-release surfaces. Between these is a giant slump consisting of at least four, largely intact, discrete blocks that have moved down-dip a distance of >500 m. A series of flows with areal extents ranging from 0.01 to 0.5 km² extends from the front of the failure. The materials represent an eroded sequence of andesite flows on the flanks of a stratovolcano. These have undergone two phases of hydrothermal alteration, and are deeply weathered to low-density (1040±80 kg m⁻³) silt (59%) and clay (35%) materials with strength properties typical of weathered silts (c=26±3 kN m⁻²; φ=42±8°). The size and location of this landslide preclude detailed geotechnical investigation of the failure. The worth of numerical stability analysis as an alternative technique in assessing the nature of the failure and hence the risk it poses to nearby communities is investigated. Sensitivity analysis identified likely conditions under which initial failure may have occurred: analyses for sensitivity to strength and earthquake acceleration needed conversion to critical combinations (F=1.0) of water table and strength/acceleration to remove the overriding influence of water table fluctuations. Failure was likely initiated either by a high water table level (83-84%), or some combination of intensity VII-IX earthquake waves together with water table heights of 40-80%. A general hazard assessment indicates that the risk associated with creep and catastrophic failure of the main mass is small, whereas the risk from flow failures near the toe of the landslide may be high. Important parameters (hydrological regime, flow failure morphology, age of initiation, and rates of movement) requiring closer investigation are identified. Development of a model is crucial to assessing the hazard associated with a feature such as that described here. With limited resources, a detailed stability analysis is a powerful tool as an initial stage in hazard analysis.     

Determining Q of near-surface materials from Rayleigh waves

High-frequency (≥2 Hz) Rayleigh wave phase velocities can be inverted to shear (S)-wave velocities for a layered earth model up to 30 m below the ground surface in many settings. Given S-wave velocity (V_S), compressional (P)-wave velocity (V_P), and Rayleigh wave phase velocities, it is feasible to solve for P-wave quality factor Q_P and S-wave quality factor Q_S in a layered earth model by inverting Rayleigh wave attenuation coefficients. Model results demonstrate the plausibility of inverting Q_S from Rayleigh wave attenuation coefficients. Contributions to the Rayleigh wave attenuation coefficients from Q_P cannot be ignored when Vs/V_P reaches 0.45, which is not uncommon in near-surface settings. It is possible to invert Q_P from Rayleigh wave attenuation coefficients in some geological setting, a concept that differs from the common perception that Rayleigh wave attenuation coefficients are always far less sensitive to Q_P than to Q_S. Sixty-channel surface wave data were acquired in an Arizona desert. For a 10-layer model with a thickness of over 20 m, the data were first inverted to obtain S-wave velocities by the multichannel analysis of surface waves (MASW) method and then quality factors were determined by inverting attenuation coefficients.     

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.     

Chemical alteration and REE mobilization in meteorites from hot and cold deserts

The effects of terrestrial weathering on REE mobilization are evaluated for a variety of uncommon meteorites found in Antarctica and in hot deserts. The meteorites analyzed include 7 non-cumulate eucrites, 10 shergottites, 3 nakhlites, 2 lunar meteorites, 4 angrites, 10 acapulcoites, 1 winonaite, and 1 brachinite. In-situ concentration measurements of lanthanides and selected other minor and trace elements were made on individual grains by secondary ion mass spectrometry (SIMS). In Antarctic meteorites, oxidation converts Ce³⁺ to Ce⁴⁺, which is less soluble than the trivalent REE, resulting in Ce anomalies. The mineral most affected is low-Ca pyroxene. However, not all grains of a given mineral are, and distinct analyses of a single grain can even yield REE patterns with and without Ce anomalies. The effect is most pronounced for Antarctic eucrites in which Ce anomalies are observed not only in individual minerals but also in whole rock samples. Although Ce anomalies are observed in meteorites from hot deserts as well, the most characteristic signs of chemical alteration in this environment are a LREE enrichment with a typical crustal signature, as well as Sr, Ba and U contaminations. These can modify the whole rock REE patterns and disturb the isotope systematics used to date these objects. The LREE contamination is highly heterogeneous, affecting some grains and not others of a given mineral (mainly olivine and low-Ca pyroxene, the two minerals with the lowest REE concentrations). The major conduit for REE movement is through shock-induced cracks and defects, and the highest levels of contamination are found in altered material filling such veins and cracks. Meteorites that experienced low shock levels and those that are highly recrystallized are the least altered.