Seismic waves generated by rockfalls contain valuable information on the properties of these events. However, as rockfalls mainly occur in mountainous regions, the generated seismic waves can be affected by strong surface topography variations. We present a methodology for investigating the influence of topography using a Spectral-Element-based simulation of 3D wave propagation in various geological media. This methodology is applied here to Dolomieu crater on the Piton de la Fournaise volcano, Reunion Island, but it can be used for other sites, taking into account local topography and medium properties. The complexity of wave fields generated by single-point forces is analyzed for different velocity models and topographies. Ground-motion amplification is studied relative to flat reference models, showing that Peak Ground Velocity (PGV) and total kinetic energy can be amplified by factors of up to 10 and 20, respectively. Simulations with Dolomieu-like crater shapes suggest that curvature variations are more influential than depth variations. Topographic effects on seismic signals from rockfalls at Dolomieu crater are revealed by inter-station spectral ratios. Results suggest that propagation along the topography rather than source direction dominates the spectral ratios and that resulting radiation patterns can be neglected. The seismic signature of single rockfall impacts is studied. Using Hertz contact theory, impact force and duration are estimated and then used to scale simulations, achieving order-of-magnitude agreement with observed signal amplitudes and frequency thresholds. Our study shows that combining Hertz theory with high-frequency seismic wave simulations on real topography improves the quantitative analysis of rockfall seismic signals.

Rockfalls generate seismic signals that can be used to detect and monitor rockfall activity. Event locations can be estimated on the basis of arrival times, amplitudes or polarization of these seismic signals. However, surface topography variations can significantly influence seismic wave propagation and hence compromise results. Here we specifically use the signature of topography on the seismic signal to better constrain the source location. Seismic impulse responses are predicted using Spectral Element based simulation of 3D wave propagation in realistic geological media. Subsequently, rockfalls are located by minimizing the misfit between simulated and observed inter-station energy ratios. The method is tested on rockfalls at Dolomieu crater, Piton de la Fournaise volcano, Reunion Island. Both single boulder impacts and distributed granular flows are successfully located. The complete rockfall trajectories are determined by analyzing the signals in sliding time windows. Results from the highest frequency band (here 13-17\,Hz) yield the best spatial resolution, making it possible to distinguish detachment positions less than 100\,m apart. By taking into account surface topography, both vertical and horizontal signal components can be used. Limitations and the noise robustness of the location method are assessed using synthetic signals. Precise representation of the topography controls the location resolution, which is not significantly affected by the assumed impact direction. Tests on the network geometry reveal best resolution when the seismometers triangulate the source. We conclude that this method can improve the monitoring of rockfall activity in real time once a simulated database for the region of interest is created.