This study delves into slip distribution on the North Tabriz Fault (NTF), a critical aspect of seismic hazard analysis due to its proximity to the Tabriz metropolis. The study operates within a uniform elastic half-space, maintaining constant values for fault geometry and regional rheological parameters throughout the research. To calculate strain boundary conditions, permanent and periodic global positioning system (GPS) data from the northwest region were utilized. The fault was constrained perpendicularly while allowing tangential movement, facilitating the determination of its annual slip rate using the boundary element method, with the Okada analytical model serving as the fundamental solution. The findings underscore the intricate relationship between the fault’s slip rate and boundary conditions, revealing a predominant right-lateral strike-slip motion. The study offers two slip rate assessments, obtained through earthquake focal mechanisms and GPS velocity data, yielding values of 5 and 5.5
mm
/
year
{\rm{mm}}/{\rm{year}}
, respectively. Importantly, the alignment of these calculated slip rates with paleo-seismological data underscores the credibility of the results generated via the boundary element method, distinguishing it as a reliable approach when compared to other numerical and analytical techniques. This research provides valuable insights into the behavior and slip dynamics of the NTF, which is pivotal for assessing seismic risks.