  In this paper, we investigate earthquake-induced landslides using a
geostatistical model that includes a latent spatial effect (LSE). The LSE
represents the spatially structured residuals in the data, which are
complementary to the information carried by the covariates. To determine
whether the LSE can capture the residual signal from a given trigger, we test
whether the LSE is able to capture the pattern of seismic shaking caused by an
earthquake from the distribution of seismically induced landslides, without
prior knowledge of the earthquake being included in the statistical model. We
assess the landslide intensity, i.e., the expected number of landslide
activations per mapping unit, for the area in which landslides triggered by the
Wenchuan (M 7.9, May 12, 2008) and Lushan (M 6.6, April 20, 2013) earthquakes
overlap. We chose an area of overlapping landslides in order to test our method
on landslide inventories located in the near and far fields of the earthquake.
We generated three different models for each earthquake-induced landslide
scenario: i) seismic parameters only (as a proxy for the trigger); ii) the LSE
only; and iii) both seismic parameters and the LSE. The three configurations
share the same set of morphometric covariates. This allowed us to study the
pattern in the LSE and assess whether it adequately approximated the effects of
seismic wave propagation. Moreover, it allowed us to check whether the LSE
captured effects that are not explained by the shaking levels, such as
topographic amplification. Our results show that the LSE reproduced the shaking
patterns in space for both earthquakes with a level of spatial detail even
greater than the seismic parameters. In addition, the models including the LSE
perform better than conventional models featuring seismic parameters only.
