Earthquakes trigger widespread landsliding in tectonically active landscapes. The effects of strong ground shaking on hillslope stability persist into the post-seismic stage; rates of landsliding remain elevated in the years following an earthquake. The mechanisms that control the spatial pattern and rate of ongoing landsliding are poorly constrained, hindering our ability to reliably forecast how landscapes and landslide hazard evolve. To address this, we undertook a detailed geotechnical investigation in which we subjected representative rock samples to dynamic loading, simulating the effects of earthquake ground shaking on hillslopes of different configuration. Our results indicate that post-seismic hillslope strength is not an intrinsic rock property; rather, it responds to the amplitude of imposed dynamic loads and the degree of pre-existing shear surface formation within the rock. This path-dependent behavior results from differences in the character of fractures generated by dynamic loads of different amplitude, and the ways in which apertures are mobilized or degraded in subsequent (post-seismic) shearing. Sensitivity to dynamic loading amplitude is greater in shallow landslides in which shear surfaces are yet to fully form; such hillslopes can be strengthened or weakened by earthquake events, depending on their characteristics. In contrast, deeper landslides on steeper hillslopes in which shear surfaces have largely developed are less likely to display differences in behavior in response to dynamic loading because strain accumulation along pre-existing fractures is dominant. Our results demonstrate the need to consider path-dependent hillslope stability in numerical models used to forecast how landscapes respond to earthquakes and how post-seismic hazard evolves
