Diapir structure and its constraint on gas hydrate accumulation in the Makran accretionary prism,offshore Pakistan

The Makran accretionary prism is located at the junction of the Eurasian Plate, Arabian Plate and Indian Plate and is rich in natural gas hydrate (NGH) resources. It consists of a narrow continental shelf, a broad continental slope, and a deformation front. The continental slope can be further divided into the upper slope, middle slope, and lower slope. There are three types of diapir structure in the accretionary prism, namely mud diapir, mud volcano, and gas chimney. (1) The mud diapirs can be grouped into two types, namely the ones with low arching amplitude and weak-medium activity energy and the ones with high arching amplitude and medium-strong activity energy. The mud diapirs increase from offshore areas towards onshore areas in general, while the ones favorable for the formation of NGH are mainly distributed on the middle slope in the central and western parts of the accretionary prism. (2) The mud volcanoes are mainly concentrated along the anticline ridges in the southern part of the lower slope and the deformation front. (3) The gas chimneys can be grouped into three types, which are located in piggyback basins, active anticline ridges, and inactive anticline ridges, respectively. They are mainly distributed on the middle slope in the central and western parts of the accretionary prism and most of them are accompanied with thrust faults. The gas chimneys located at different tectonic locations started to be active at different time and pierced different horizons. The mud diapirs, mud volcanoes, and gas chimneys and thrust faults serve as the main pathways of gas migration, and thus are the important factors that control the formation, accumulation, and distribution of NGH in the Makran accretionary prism. Mud diapir/gas chimney type hydrate develop in the middle slope, mud volcano type hydrate develop in the southern lower slope and the deformation front, and stepped accretionary prism type hydrate develop on the central and northern lower slope. The middle slope, lower slope and deformation front in the central and western parts of the Makran accretionary prism jointly constitute the NGH prospect area.     

Reservoir Potential Assessment of the Tortonian–Calabrian Sediments from Offshore Makran(Southwest Pakistan) through a Multidisciplinary Approach

The Tortonian–Calabrian strata of Offshore Makran(Pakistan) is investigated for the purpose of reservoir assessment. The stratigraphy and petrophysics indicate that the Neogene sediments have tight reservoir potential with porosities in the range of 3%–9% and 40%–50% water saturation. The mineralogical cross plots indicate a mixed lithology with an abundance of dolomite and calcite, together with minor quartz content and shale. The seismic interpretation demonstrates medium–high amplitude reflec...     

Seismic quality factors across a bottom simulating reflector in the Makran Accretionary Prism, Arabian Sea

The hydrate-bearing sediments above the bottom simulating reflector (BSR) are associated with low attenuation or high quality factor (Q), whereas underlying gas-bearing sediments exhibit high attenuation. Hence, estimation of Q can be important for qualifying whether a BSR is related to gas hydrates and free-gas. This property is also useful for identifying gas hydrates where detection of BSR is dubious. Here, we calculate the interval Q for three submarine sedimentary layers bounded by seafloor, BSR, one reflector above and another reflector below the BSR at three locations with moderate, strong and no BSR along a seismic line in the Makran accretionary prism, Arabian Sea for studying attenuation (Q⁻¹) characteristics of sediments. Interval Q for hydrate-bearing sediments (layer B) above the BSR are estimated as 191 ± 11, 223 ± 12, and 117 ± 5, whereas interval Q for the underlying gas-bearing sediments (layer C) are calculated as 112 ± 7, 107 ± 8 and 124 ± 11 at moderate, strong and no BSR locations, respectively. The large variation in Q is observed at strong BSR. Thus Q can be used for ascertaining whether the observed BSR is due to gas hydrates, and for identifying gas hydrates at places where detection of BSR is rather doubtful. Interval Q of 98 ± 4, 108 ± 5, and 102 ± 5, respectively, at moderate, strong and no BSR locations for the layer immediately beneath the seafloor (layer A) show almost uniform attenuation.     

P-T evolution of metapelites from the Bajgan complex in the Makran accretionary prism,south eastern Iran

The Bajgan Complex, one of the basement constituents of the arc massif in Iranian Makran forms a rugged, deeply incised terrain. The complex consists of pelitic schists with minor psammitic and basic schists, calc silicate rocks, amphibolites, marbles, metavolcanosediments, mafic and felsic intrusives as well as ultramafic rocks. Metapelitic rocks show an amphibolite facies regional metamorphism and contain garnet, biotite, white mica, quartz, albite ± rutile ± apatite. Thermobarometry of garnet schist yields pressure of more than 9 kbar and temperatures between 560 and 675 °C. The geothermal gradient obtained for the peak of regional metamorphism is 19 °C/km, corresponding to a depth of ca. 31 km. Replacement of garnet by chlorite and epidote suggest greenschist facies metamorphism due to a decrease in temperature and pressure through exhumation and retrograde metamorphism (370–450 °C and 3–6 kbar). The metapelitic rocks followed a ‘clockwise’ P–T path during metamorphism, consistent with thermal decline following tectonic thickening. The formation of medium-pressure metamorphic rocks is related to presence of active subduction of the Neotethys Oceanic lithosphere beneath Eurasia in the Makran.     

Characteristics of major and trace elements in surface sediments of the Makran Accretionary Prism,Pakistan and their implications for natural gas hydrates

To accurately identify the natural gas hydrates (NGH) in the sea area of the Makran Accretionary Prism, Pakistan, this paper presents the testing and analysis of major and trace elements in sediment samples taken from two stations (S₂ and S₃) in the area by the China Geological Survey. As shown by testing results, all major elements are slightly different in content between the two stations except SiO₂ and CaO. This also applies to the trace elements that include Sr and Ba primarily and Cr, Ni and Zn secondarily. It can be concluded in this study that the tectonic setting of the Makran Accretionary Prism is dominated by oceanic island arc and that provenance of the Makran Accretionary Prism is dominated by felsic igneous provenance, which is at the initial weathering stage and mainly consists of granodiorite. Besides terrigenous detritus, there are sediments possibly originating from Makran-Bela Ophiolite from the northwestern part and Murray Ridge igneous rocks from the southeastern part. The V/Cr, Ni/Co, and V/(V+Ni) ratios indicate that sediments of the two stations are in an oxidation-suboxidation environment. However, the authors infer that the sedimentary environment of the sediments 3.0 m below the seafloor tends to be gradually transformed into a reduction environment by comparison with the Qiongdongnan Basin in the South China Sea where NGH has been discovered. The sediments in the Makran Accretionary Prism are rich in organic matter, with total organic carbon (TOC) content greater than 1%. According to comprehensive research, the organic matter in the sediments mainly originates from marine algae and has high TOC content, which is favorable for the formation of NGH.     

Iranshahr blueschist: subduction of the inner Makran oceanic crust

The Makran accretionary prism in SE Iran and SW Pakistan is one of the most extensive subduction accretions on Earth. It is characterized by intense folding, thrust faulting and dislocation of the Cenozoic units that consist of sedimentary, igneous and metamorphic rocks. Rock units forming the northern Makran ophiolites are amalgamated as a mélange. Metamorphic rocks, including greenschist, amphibolite and blueschist, resulted from metamorphism of mafic rocks and serpentinites. In spite of the geodynamic significance of blueschist in this area, it has been rarely studied. Peak metamorphic phases of the northern Makran mafic blueschist in the Iranshahr area are glaucophane, phengite, quartz±omphacite+epidote. Post peak minerals are chlorite, albite and calcic amphibole. Blueschist facies metasedimentary rocks contain garnet, phengite, albite and epidote in the matrix and as inclusions in glaucophane. The calculated P–T pseudosection for a representative metabasic glaucophane schist yields peak pressure and temperature of 11.5–15 kbar at 400–510 °C. These rocks experienced retrograde metamorphism from blueschist to greenschist facies (350–450 °C and 7–8 kbar) during exhumation. A back arc basin was formed due to northward subduction of Neotethys under Eurasia (Lut block). Exhumation of the high‐pressure metamorphic rocks in northern Makran occurred contemporarily with subduction. Several reverse faults played an important role in exhumation of the ophiolitic and HP‐LT rocks. The presence of serpentinite shows the possible role of a serpentinite diapir for exhumation of the blueschist. A tectonic model is proposed here for metamorphism and exhumation of oceanic crust and accretionary sedimentary rocks of the Makran area. Vast accretion of subducted materials caused southward migration of the shore.     

Modeling the near-field effects of the worst-case tsunami in the Makran subduction zone

As a first step towards the development of inundation maps for the northwestern Indian Ocean, we simulated the near-field inundation of two large tsunami in the Makran subduction zone (MSZ). The tsunami scenarios were based on large historical earthquakes in the region. The first scenario included the rupture of about 500 km of the plate boundary in the eastern MSZ, featuring a moment magnitude of M_w 8.6. The second scenario involved the full rupture of the plate boundary resulting from a M_w 9 earthquake. For each scenario, the distribution of tsunami wave height along the coastlines of the region is presented. Also, detailed runup modeling was performed at four main coastal cities in the region for the second scenario. To investigate the possible effect of splay fault branching on tsunami wave height, a hypothetical splay fault was modeled which showed that it can locally increase the maximum wave height by a factor of 2. Our results showed that the two tsunami scenarios produce a runup height of 12-18 m and 24-30 m, respectively. For the second scenario, the modeled inundation distance was between 1 and 5 km.     

Misattributed tsunami 3: the Mw 7.7 2013.9.24 SE Pakistan earthquake

The M_w 7.7 earthquake that struck SE Pakistan on 24 September 2013 at 11.29.48 UTC was a sinistral strike-slip event on a branch of the Ornach-Nal-Chaman fault system which hereabouts separates the Eurasian Plate from the Indian Plate. Although the focus was at a depth of 15 km and 400 km inland the earthquake was accompanied by the emergence of an island off the Makran coast and the generation of a tsunami with a peak amplitude of 27 cm at Muscat (Oman) and 20 cm at Chah Bahar (Iran). At DART marine buoy 23228 in the Indian Ocean 500 km to the south a series of seismic Rayleigh waves about 4 min after the main shock was followed 54 min later by a tsunami with a peak amplitude of 1 cm. The Rayleigh series is here attributed to seafloor vibration during accelerated subduction of the Arabian Plate beneath the Eurasian Plate, and the tsunami to the development or reactivation of one or more reverse faults on the seaward portion of the Makran imbricate fan. As in the 2010.2.27 M_w 8.8 Maule (Chile), the 2004.12.26 M_w 9.2 Sumatra–Andaman, the 2005.3.28 M_w 8.7 Nias (Indonesia) and the 2011.3.11 M_w 9.0 Tohoku (Japan) earthquakes, the link between tsunami generation and slip on the megathrust is thus very indirect, to the detriment of attempts to mitigate coastal hazards using teleseismic data when nearshore seafloor monitoring would probably prove more effective.