In January and February 2001, El Salvador was hit by two strong earthquakes that triggered thousands of landslides, causing 1,259 fatalities and extensive damage. The analysis of aerial and SPOT-4 satellite images taken a few days after the events allowed us to map 6,491 coesismic landslides, which occurred in 14 study areas extending for about 400 km2. Four different Multivariate Adaptive Regression Splines (MARS) models were produced by using different covariate sets and landslide inventories, the latter containing the slope failures triggered by an extreme rainfall event of November 2009 and those induced by the earthquakes of 2001. Moreover, two validation scenarios were employed, including the 25% and 95% of the mapped landslides, respectively. The results of our experiment revealed that: (i) MARS algorithm provides reliable predictions of coesismic landslides; (ii) models calibrated with rainfall-induced landslides predict with acceptable accuracy landslides caused by deep earthquakes and distributed over vast areas; (iii) the best accuracy is achieved by models trained with both preparatory and trigger variables; (iv) a small portion of the landslides produced by an earthquake can be used to calibrate MARS predictive models that help to identify slopes where yet unreported landslides may have occurred.