This review synthesizes existing research to elucidate the factors driving the distinct tectonic behaviors in the western and eastern Makran subduction zone, focusing on seismic activity, uplift rate, convergence rate, coupling, and subduction angle. The literature identifies the asymmetry in pressure and the variation in subduction angles between the western and eastern parts of the Makran as key factors in defining the region’s tectonic patterns. The western region has a steeper subduction angle, resulting in lower pressure, reduced coupling, and decreased seismic activity. This disparity arises from different interactions between the subducted and overriding plates. This article offers an overview of the Makran subduction zone, identifies some knowledge gaps, and suggests directions for future research to improve our understanding of this complex geological region. The review highlights the need for more comprehensive GPS stations and targeted studies on subduction dip angles to better understand the region’s tectonic dynamics.

On 14 November 2021, a doublet earthquake, each event of which had an Mw of 6.1, struck near Fin in the Simply Folded Belt (SFB) in southern Iran. The first quake occurred at 12:07:04 UTC, followed by a second one just a minute and a half later. The SFB is known for its blind thrust faults, typically not associated with surface ruptures. These earthquakes are usually linked to the middle and lower layers of the sedimentary cover. Identifying the faults that trigger earthquakes in the region remains a significant challenge and is subject to high uncertainty. This study aims to identify and determine the fault(s) that may have caused the doublet earthquake. To achieve this goal, we utilized the DInSAR method using Sentinel-1 to detect deformation, followed by finite-fault inversion and magnetic interpretation to determine the location, geometry, and slip distribution of the fault(s). Bayesian probabilistic joint inversion was used to model the earthquake sources and derive the geometric parameters of potential fault planes. The study presents two potential fault solutions—one dipping to the north and the other to the south. Both solutions showed no significant difference in strike and fault location, suggesting a single fault. Based on the results of the seismic inversion, it appears that a north-dipping fault with a strike, dip, and rake of 257°, 74°, and 77°, respectively, is more consistent with the geological setting of the area. The fault plane has a width of roughly 3.6 km, a length of 13.4 km, and a depth of 5.6 km. Our results revealed maximum displacements along the radar line of sight reaching values of up to −360 mm in the ascending orbit, indicating an unknown fault with horizontal displacements at the surface ranging from −144 to 170 mm and maximum vertical displacements between −204 and 415 mm. Aeromagnetic data for Iran were utilized with an average flight-line spacing of 7.5 km. The middle of the data observation period was considered to apply the RTP filter, and the DRTP method was used. We calculated the gradient of the residual anomaly in the N-S direction due to the direction of the existing faults and folds. The gradient map identified the fault and potential extension of the observed anomalies related to a fault with an ENE-WSW strike, which could extend to the ~ E-W. We suggest that earthquakes occur in the sedimentary cover of the SFB where subsurface faulting is involved, with Hormuz salt acting as an important barrier to rupture. The multidisciplinary approach used in this study, including InSAR and magnetic data, underscores the importance of accurate fault characterization. These findings provide valuable insights into the seismic hazard of the area.

We have developed a multiscale approach for solving 2D and 3D nonlinear inverse problems of gravity data in estimating the basement topography. The inversion is carried out in two stages in which the long-wavelength features of the basement are first estimated from smoothed gravity data via a stochastic optimization algorithm. The solution of this stage is used as the starting model for a deterministic optimization algorithm to reconstruct the short-wavelength features from the full-spectrum gravity data. The forward problem is capable of handling lateral and vertical variations in the density of sediments. Two cases are considered regarding prior knowledge about the density: (1) The density contrast between sediments at the surface and the underlying basement and its vertical variations are a priori known, and (2) only the density contrast at the surface is known with its vertical gradient to be recovered in the inversion. In the former case, the unknowns of the problem are the depths, whereas in the latter case, they are the depths and density gradients defined individually for each prism. Therefore, the inverse problem is ill-posed and has many local minima. The stochastic optimization algorithm uses a random initial model and estimates a coarse model of the basement topography. By repeating the stochastic inversion, an ensemble of solutions is formed defining an equivalent domain in the model space supposed to be within the neighborhood of the global minimum of which several starting solutions are extracted for the secondary deterministic inversion. The presented methodology has been tested successfully in converging to the global minima in 2D and 3D cases with 50 and 2352 total number of prisms, respectively. Finally, the inversion algorithm is used to calculate the thickness of the sediments in the South Caspian Basin using the EIGEN-6c4 global gravity model.