Mechanisms that control seismic activity in low strain rate areas such as western Europe remain poorly understood. For example, in spite of low shortening rates of <0.5 mm/yr, the Western Alps and the Pyrenees are underlain by moderate but frequent seismicity detectable by instruments. Beneath the elevated part of these mountain ranges, analysis of earthquake focal mechanisms indicates extension, which is commonly interpreted as the result of gravitational collapse. Here we show that erosional processes are the predominant control on present-day deformation and seismicity. We demonstrate, using finite element modeling, that erosion induces extension and rock uplift of the elevated region of mountain ranges accommodating relatively low overall convergence. Our results suggest that an erosion rate of ∼1 mm/yr can lead to extension in mountain ranges accommodating significant shortening of <3 mm/yr. Based on this study, the seismotectonic framework and seismic hazard assessment for low strain rate areas need to be revisited, because erosion-related earthquakes could increase seismic hazard.

The Pyrenean mountain range is a slowly deforming belt with continuous and moderate seismic activity. To quantify its deformation field, we present the velocity field estimated from a GPS survey of the Pyrenees spanning 18 yr. The PotSis and ResPyr networks, including a total of 85 GPS sites, were installed and first measured in 1992 and 1995–1997, respectively, and remeasured in 2008 and 2010. We obtain a deformation field with velocities less than 1 mm yr⁻¹ across the range. The estimated velocities for individual stations do not differ significantly from zero with 95 per cent confidence. Even so, we estimate a maximum extensional horizontal strain rate of 2.0 ± 1.7 nanostrain per year in a N–S direction in the western part of the range. We do not interpret the vertical displacements due to their large uncertainties. In order to compare the horizontal strain rates with the seismic activity, we analyse a set of 194 focal mechanisms using three methods: (i) the ‘r’ factor relating their P and T axes, (ii) the stress tensors obtained by fault slip inversion and (iii) the strain-rate tensors. Stress and strain-rate tensors are estimated for: (i) the whole data set, (ii) the eastern and western parts of the range separately, and (iii) eight zones, which are defined based on the seismicity and the tectonic patterns of the Pyrenees. Each of these analyses reveals a lateral variation of the deformation style from compression and extension in the east to extension and strike-slip in the west of the range. Although the horizontal components of the strain-rate tensors estimated from the seismic data are slightly smaller in magnitude than those computed from the GPS velocity field, they are consistent within the 2σ uncertainties. Furthermore, the orientations of their principal axes agree with the mapped active faults.