We present a novel inverse method for discriminating regional deformation and long-term fault creep by inversion of GNSS velocities observed at the spatial scale of intraplate faults by exploiting the different spatial signatures of these two mechanisms. In doing so our method provides a refined estimate of the upper bound of the strain accumulation process. As case study, we apply this method to a six year GNSS campaign (2003–2008) set up in the southern portion of the Pollino Range over the Castrovillari and Pollino faults. We show that regional deformation alone cannot explain the observed deformation pattern and implies high geodetic strain rate, with a WSW-ENE extension of 86±41×10−9/yr. Allowing for the possibility of fault creep, the modelling of GNSS velocities is consistent with their uncertainties and they are mainly explained by a shallow creep over the Pollino fault, with a normal/strike-slip mechanism up to 5 mm/yr. The regional strain rate decrease by about 70 percent and is characterized by WNW-ESE extension of 24±28×10−9/yr. The large uncertainties affecting our estimate of regional strain rate do not allow infering whether the tectonic regime of the area is extensional or strike-slip, although the latter is slightly more likely.