We present a new global model of plate motions and strain rates in plate boundary zones constrained by horizontal geodetic velocities. This Global Strain Rate Model (GSRM v.2.1) is a vast improvement over its predecessor both in terms of amount of data input as in an increase in spatial model resolution by factor of ∼2.5 in areas with dense data coverage. We determined 6739 velocities from time series of (mostly) continuous GPS measurements; i.e., by far the largest global velocity solution to date. We transformed 15,772 velocities from 233 (mostly) published studies onto our core solution to obtain 22,511 velocities in the same reference frame. Care is taken to not use velocities from stations (or time periods) that are affected by transient phenomena; i.e., this data set consists of velocities best representing the interseismic plate velocity. About 14% of the Earth is allowed to deform in 145,086 deforming grid cells (0.25° longitude by 0.2° latitude in dimension). The remainder of the Earth's surface is modeled as rigid spherical caps representing 50 tectonic plates. For 36 plates we present new GPS-derived angular velocities. For all the plates that can be compared with the most recent geologic plate motion model, we find that the difference in angular velocity is significant. The rigid-body rotations are used as boundary conditions in the strain rate calculations. The strain rate field is modeled using the Haines and Holt method, which uses splines to obtain an self-consistent interpolated velocity gradient tensor field, from which strain rates, vorticity rates, and expected velocities are derived. We also present expected faulting orientations in areas with significant vorticity, and update the no-net rotation reference frame associated with our global velocity gradient field. Finally, we present a global map of recurrence times for Mw=7.5 characteristic earthquakes.

Triple junctions involving non-subductable plates extend beyond local implications, crucial for studying the geology of convergent plate boundary zones. However, kinematic models overlook Cyprus-Anatolia motion due to limited geodetic constraints. Our study area comprises Cyprus, southern Turkey, and the Levant coast, focusing on the Kahramanmaraş triple junction where a destructive earthquake sequence occurred on February 6, 2023. We present precise positioning data merged with published velocities, constructing an up-to-date velocity field for the interseismic period. Employing two kinematic approaches, we analyze its tectonic implications. In Cyprus, we find the relative motion of Africa (Sinai Plate) and Anatolia is partitioned between convergence in the Cyprus subduction, with rate 3.5-6.2 mm/yr, progressively decreasing from west to east and left-lateral transpressive Kyrenia fault, situated along the northern coast of Cyprus, with rate 3.3-4.2 mm/yr. The relative strike-slip motion between Arabia and Anatolia is partitioned between the East Anatolian Fault (slip rates 5.2-6.2 mm/yr) and some secondary faults such as Çardak and Malatya faults (slip rates 2.0-1.7 mm/yr respectively) and causes distributed deformation for a 50-60 km wide region. The largest second invariant strain rate tensors from the continuum kinematic model also coincide with the same region which can be defined as the East Anatolian shear zone. A shear partitioning system exists around the Kahramanmaraş triple junction, from Cyprus to southeast Turkey. The Levant Fault has a 3.5-4.7 mm/yr left-lateral slip rate, decreasing northward as part of it is transferred to offshore faults. Strain rates appear to be relatively small in the Taurus range and Adana/Cilicia basin, transitioning from extensional/transtensional to compressional from east to west. The up component of GNSS velocities on the Taurus karstic plateau reaches 1.2 mm/yr, lower than Quaternary uplift rates on its southern edge. We examine Taurus uplift and Adana/Cilicia basin subsidence, considering slow convergence through earthquake distribution, seismic tomography, and modeling.