We characterize and compare creep preceding and following the complex 2011 Pampahasi landslide (∼40 Mm³±50 %) in the city of La Paz, Bolivia, using spaceborne radar interferometry (InSAR) that combines displacement records from both distributed and point scatterers. The failure remobilized deposits of an ancient complex landslide in weakly cemented, predominantly fine-grained sediments and affected ∼1.5 km² of suburban development. During the 30 months preceding failure, about half of the toe area was creeping at 3–8 cm a⁻¹ and localized parts of the scarp area showed displacements of up to 14 cm a⁻¹. Changes in deformation in the 10 months following the landslide demonstrate an increase in slope activity and indicate that stress redistribution resulting from the discrete failure decreased stability of parts of the slope. During that period, most of the landslide toe and areas near the head scarp accelerated, respectively, to 4–14 and 14 cm a⁻¹. The extent of deformation increased to cover most, or probably all, of the 2011 landslide as well as adjacent parts of the slope and plateau above. The InSAR-measured displacement patterns, supplemented by field observations and optical satellite images, reveal complex slope activity; kinematically complex, steady-state creep along pre-existing sliding surfaces accelerated in response to heavy rainfall, after which slightly faster and expanded steady creeping was re-established. This case study demonstrates that high-quality ground-surface motion fields derived using spaceborne InSAR can help to characterize creep mechanisms, quantify spatial and temporal patterns of slope activity, and identify isolated small-scale instabilities; such details are especially useful where knowledge of landslide extent and activity is limited. Characterizing slope activity before, during, and after the 2011 Pampahasi landslide is particularly important for understanding landslide hazard in La Paz, half of which is underlain by similar large paleolandslides.

The extent and behaviour of glaciers during the mid-Piacenzian warm period illustrate the sensitivity of the cryosphere to atmospheric CO₂concentrations above pre-industrial levels. Knowledge of glaciation during this period is restricted to globally or regionally averaged records from marine sediments and to sparse terrestrial glacial deposits in mid-to-high latitudes. Here we expand the Pliocene glacial record to the tropics by reporting recurrent large-scale glaciation in the Bolivian Andes based on stratigraphic and paleomagnetic analysis of a 95-m sequence of glacial sediments underlying the 2.74-Ma Chijini Tuff. Paleosols and polarity reversals separate eight glacial diamictons, which we link to cold periods in the benthic oxygen isotope record. The glaciations appear to coincide with the earliest glacial activity at high northern latitudes and with events in Antarctica, including the strong M2 cold peak and terminal Pliocene climate deterioration. This concordance suggests inter-hemispheric climate linkages during the late Pliocene and requires that the Central Andes were at least as high in the late Pliocene as today. Our record fills a critical gap in knowledge of Earth systems during the globally warm mid-Piacenzian and suggests a possible driver of faunal migration preceding the large-scale biotic interchange in the Americas during the earliest Pleistocene.

In this work, a simple methodology for preliminarily assessing the magnitude of potential landslide-induced impulse waves’ attenuation in mountain lakes is presented. A set of metrics is used to define the geometries of theoretical mountain lakes of different sizes and shapes and to simulate impulse waves in them using the hydrodynamic software Flow-3D. The modeling results provide the ‘wave decay potential’, a ratio between the maximum wave amplitude and the flow depth at the shoreline. Wave decay potential is highly correlated with what is defined as the ‘shape product’, a metric that represents lake geometry. The relation between these two parameters can be used to evaluate wave dissipation in a natural lake given its geometric properties, and thus estimate expected flow depth at the shoreline. This novel approach is tested by applying it to a real-world event, the 2007 landslide-generated wave in Chehalis Lake (Canada), where the results match well with those obtained using the empirical equation provided by ETH Zurich (2019 Edition). This work represents the initial stage in the development of this method, and it encourages additional research and modeling in which the influence of the impacting characteristics on the resulting waves and flow depths is investigated.