We present a new automated earthquake detection and location method based on beamforming (or back projection) and template matching and apply it to study the seismicity of the Southwestern Alps. We use beamforming with prior knowledge of the 3-D variations of seismic velocities as a first detection run to search for earthquakes that are used as templates in a subsequent matched-filter search. Template matching allows us to detect low signal-to-noise ratio events and thus to obtain a high spatiotemporal resolution of the seismicity in the Southwestern Alps. We describe how we address the problem of false positives in energy-based earthquake detection with supervised machine learning and how to best leverage template matching to iteratively r...