Normal faulting mechanisms observed in the northern foothills of the Central–Western Pyrenees are remarkable, since one expects thrust faults at a convergent plate boundary. To understand the mechanisms involved, we used numerical modeling and investigated the impact of the following processes: gravitational potential energy associated with topography and dense crustal blocks; isostatic compensation in response to denudation and/or sedimentation. To decipher the effects of each process, we designed three model geometries and added a pre-existing weak fault where most of the seismicity occurs. To evaluate our model results, we derived the fault slip rate from the focal mechanisms in the region where we have the fault in our model. We found a slip rate of ∼15m/Ma, which is in agreement with our modelling results. We conclude that flexural rebound induced by surface processes is able to explain the seismicity in Central–Western Pyrenees.