New high-resolution 2D seismic imaging of fluid escape structures in the Makran subduction zone,Arabian Sea

Seabed fluid escape is active in the Makran subduction zone, Arabian Sea. Based on the new high-resolution 2D seismic data, acoustic blanking zones and seafloor mounds are identified. Acoustic blanking zones include three kinds of geometries: Bell-shaped, vertically columnar and tilted zones. The bell-shaped blanking zone is characterized by weak and discontinuous reflections in the interior and up-bending reflections on the top, interpreted as gas chimneys. Vertically columnar blanking zone is interpreted as side-imaged gas chimneys associated with focused fluid flow and topped by a seafloor anomaly expressed as a localized reflection discontinuity, which may together serve as a vent structure. Tilted acoustic blanking zone could be induced by accretionary thrust activity and rapid sedimentation surrounding slope. Seafloor mounds occur at the sites of bell-shaped acoustic blanking zone and may be associated with the material intrusion. Bottom simulating refectors (BSRs) are widely distributed and exhibit a series of characteristics including diminished amplitude, low continuity as well as local shoaling overlapping with these acoustic blanking zones. The large amount of gases dissociated from the gas hydrates migrated upwards and then arrived at the near-seafloor sediments, followed by the formation of the gas hydrates and hence the seafloor mound.     

Pore-water geochemistry in methane-seep sediments of the Makran accretionary wedge off Pakistan: Possible link to subsurface methane hydrate

Cold seeps are pervasive along the continental margin worldwide,and are recognized as hotspots for elemental cycling pathway on Earth.In this study,analyses of pore water geochemical compositions of one-400 cm piston core(S3) and the application of a mass balance model are conducted to assess methane-associated biogeochemical reactions and uncover the relationship of methane in shallow sediment with gas hydrate reservoir at the Makran accretionary wedge off Pakistan.The results revealed that approximately 77% of sulfate is consumed by the predominant biogeochemical process of anaerobic oxidation of methane.However,the estimated sulfate-methane interface depth is-400 cm below sea floor with the methane diffusive flux of 0.039 mol/(m~2·a),suggesting the activity of methane seepage.Based on the δ~(13)C_(DIC) mass balance model combined with the contribution proportion of different dissolved inorganic carbon sources,this study calculated the δ~(13)C of the exogenous methane to be-57.9‰,indicating that the exogenous methane may be a mixture source,including thermogenic and biogenic methane.The study of pore water geochemistry at Makran accretionary wedge off Pakistan may have considerable implications for understanding the specific details on the dynamics of methane in cold seeps and provide important evidence for the potential occurrence of subsurface gas hydrate in this area.     

Historical tsunami in the Makran Subduction Zone off the southern coasts of Iran and Pakistan and results of numerical modeling

Tsunami hazard in the Makran Subduction Zone (MSZ), off the southern coasts of Iran and Pakistan, was studied by numerical modeling of historical tsunami in this region. Although the MSZ triggered the second deadliest tsunami in the Indian Ocean, among those known, the tsunami hazard in this region has yet to be analyzed in detail. This paper reports the results of a risk analysis using five scenario events based on the historic records, and identifies a seismic gap area in western Makran off the southern coast of Iran. This is a possible site for a future large earthquake and tsunami. In addition, we performed numerical modeling to explain some ambiguities in the historical reports. Based on the modeling results, we conclude that either the extreme run-up of 12–15 m assigned for the 1945 Makran tsunami in the historical record was produced by a submarine landslide triggered by the parent earthquake, or that these reports are exaggerated. The other possibility could be the generation of the huge run-up heights by large displacements on splay faults. The results of run-up modeling reveal that a large earthquake and tsunami in the MSZ is capable of producing considerable run-up heights in the far field. Therefore, it is possible that the MSZ was the source of the tsunami encountered by a Portuguese fleet in Dabhul in 1524.     

Effective medium modeling to assess gas hydrate and free-gas evident from the velocity structure in the Makran accretionary prism, offshore Pakistan

Seismic properties of sediments are strongly influenced by pore fluids. Stiffness of unconsolidated marine sediment increases with the presence of gas hydrate and decreases with the presence of gas. A strong bottom-simulating reflector (BSR) observed on a seismic profile in the Makran accretionary prism, offshore Pakistan, indicates the presence of gas hydrate and free-gas across the BSR. Elastic properties of gas depend largely on pressure and temperature. We, therefore, first determine the elastic modulus of gas at pressure and temperature calculated at the BSR depth in the study region. The interval velocities derived from the seismic data are interpreted by the effective medium theory, which is a combination of self-consistent approximation and differential effective medium theories, together with a smoothing approximation, for assessment of gas hydrate and free-gas. The results show the saturations of gas hydrate and free-gas as 22 and 2.4% of pore space, respectively, across the BSR.     

冷泉系统甲烷气体渗漏的声学特征及潜在环境效应

The amount of methane leaked from deep sea cold seeps is enormous and potentially affects the global warming,ocean acidification and global carbon cycle. It is of great significance to study the methane bubble movement and dissolution process in the water column and its output to the atmosphere. Methane bubbles produce strong acoustic impedance in water bodies, and bubble strings released from deep sea cold seeps are called "gas flares"which expressed as flame-like strong backscatter in the water column. We characterized the morphology and movement of methane bubbles released into the water using multibeam water column data at two cold seeps. The result shows that methane at site I reached 920 m water depth without passing through the top of the gas hydrate stability zone(GHSZ, 850 m), while methane bubbles at site II passed through the top of the GHSZ(597 m) and entered the non-GHSZ(above 550 m). By applying two methods on the multibeam data, the bubble rising velocity in the water column at sites I and II were estimated to be 9.6 cm/s and 24 cm/s, respectively. Bubble velocity is positively associated with water depth which is inferred to be resulted from decrease of bubble size during methane ascending in the water. Combined with numerical simulation, we concluded that formation of gas hydrate shells plays an important role in helping methane bubbles entering the upper water bodies, while other factors, including water depth, bubble velocity, initial kinetic energy and bubble size, also influence the bubble residence time in the water and the possibility of methane entering the atmosphere. We estimate that methane gas flux at these two sites is 0.4×10~6–87.6×10~6 mol/a which is extremely small compared to the total amount of methane in the ocean body, however, methane leakage might exert significant impact on the ocean acidification considering the widespread distributed cold seeps. In addition, although methane entering the atmosphere is not observed, further research is still needed to understand its potential impact on increasing methane concentration in the surface seawater and gas-water interface methane exchange rate, which consequently increase the greenhouse effect.     

Structural elements of the Makran region, Oman sea and their potential relevance to tsunamigenisis

The character of convergence along the Arabian–Iranian plate boundary changes radically eastward from the Zagros ranges to the Makran region. This appears to be due to collision of continental crust in the west, in contrast to subduction of oceanic crust in the east. The Makran subduction zone with a length of about 900 km display progressively older and highly deformed sedimentary units northward from the coast, together with an increase in elevation of the ranges. North of the Makran ranges are large subsiding basins, flanked to the north by active volcanoes. Based on 2D seismic reflection data obtained in this study, the main structural provinces and elements in the Gulf of Oman include: (i) the structural elements on the northeastern part of the Arabian Plate and, (ii) the Offshore Makran Accretionary Complex. Based on detailed analysis of these data on the northeastern part of the Arabian Plate five structural provinces and elements—the Musendam High, the Musendam Peneplain, the Musendam Slope, the Dibba Zone, and the Abyssal Plain have been identified. Further, the Offshore Makran Accretionary Complex shown is to consist Accretionary Prism and the For-Arc Basin, while the Accretionary Prism has been subdivided into the Accretionary Wedge and the Accreted/Colored Mélange. Lastly, it is important to note that the Makran subduction zone lacks the trench. The identification of these structural elements should help in better understanding the seismicity of the Makran region in general and the subduction zone in particular. The 1945 magnitude 8.1 tsunamigenic earthquake of the Makran and some other historical events are illustrative of the coastal region’s vulnerability to future tsunami in the area, and such data should be of value to the developing Indian Ocean Tsunami Warning System.     

Deformation history and processes during accretion of seamounts in subduction zones: The example of the Durkan Complex (Makran,SE Iran)

The Durkan Complex is a tectonic element of the Makran Accretionary Prism (SE Iran) that includes fragments of Late Cretaceous seamounts. In this paper, the results of map- to micro-scale structural studies of the western Durkan Complex are presented with the aim to describe its structural and tectono-metamorphic evolution. The Durkan Complex consists of several tectonic units bordered by mainly NNW-striking thrusts. Three main deformation phases (D₁, D₂, and D₃) are distinguished and likely occurred from the Late Cretaceous to the Miocene–Pliocene. D₁ is characterized by sub-isoclinal to close and W-verging folds associated with an axial plane foliation and shear zone along the fold limbs. This phase records the accretion of fragments of the seamount within the Makran at blueschist facies metamorphic conditions (160–300 °C and 0.6 – 1.2 GPa). D₂ is characterized by open to close folds with sub-horizontal axial plane that likely developed during the exhumation of previously accreted seamount fragments. An upper Paleocene – Eocene siliciclastic succession unconformably sealed the D₁ and D₂ structures and is, in turn, deformed by W-verging thrust faults typical of D₃. The latter likely testifies for a Miocene – Pliocene tectonic reworking of the accreted seamount fragments with the activation of out of sequence thrusts. Our results shed light on the mechanism of accretion of seamount materials in the accretionary prisms, suggesting that seamount slope successions favour the localization and propagation of the basal décollement. This study further confirms that the physiography of the subducting plates plays a significant role in the tectonic evolution of the subduction complexes.     

Preliminary estimation of the tsunami hazards associated with the Makran subduction zone at the northwestern Indian Ocean

We present a preliminary estimation of tsunami hazard associated with the Makran subduction zone (MSZ) at the northwestern Indian Ocean. Makran is one of the two main tsunamigenic zones in the Indian Ocean, which has produced some tsunamis in the past. Northwestern Indian Ocean remains one of the least studied regions in the world in terms of tsunami hazard assessment. Hence, a scenario-based method is employed to provide an estimation of tsunami hazard in this region for the first time. The numerical modeling of tsunami is verified using historical observations of the 1945 Makran tsunami. Then, a number of tsunamis each resulting from a 1945-type earthquake (M _w 8.1) and spaced evenly along the MSZ are simulated. The results indicate that by moving a 1945-type earthquake along the MSZ, the southern coasts of Iran and Pakistan will experience the largest waves with heights of between 5 and 7 m, depending on the location of the source. The tsunami will reach a height of about 5 m and 2 m in northern coast of Oman and eastern coast of the United Arab Emirates, respectively.     

Geochemical characteristics, Ar–Ar ages and original tectonic setting of the Band-e-Zeyarat/Dar Anar ophiolite, Makran accretionary prism, S.E. Iran

The Makran accretionary prism in southeastern Iran contains extensive Mesozoic zones of melange and large intact ophiolites, representing remnants of the Tethys oceanic crust that was subducted beneath Eurasia. To the north of the Makran accretionary prism lies the Jaz Murian depression which is a subduction-related back-arc basin. The Band-e-Zeyarat/Dar Anar ophiolite is one of the ophiolite complexes; it is located on the west side of the Makran accretionary prism and Jaz Murian depression, and is bounded by two major fault systems. The principal rock units of this complex are a gabbro sequence which includes low- and high-level gabbros, an extensive sheeted diabase dike sequence, late intrusive rocks which consist largely of trondhjemites and diorites, and volcanic rocks which are largely pillow basalts interbedded with pelagic sedimentary rocks, including radiolarian chert. Chondrite- and primitive-mantle-normalized incompatible trace element data and age-corrected Nd, Pb, and Sr isotopic data indicate that the Band-e-Zeyarat/Dar Anar ophiolite was derived from a midocean ridge basalt-like mantle source. The isotopic data also reveal that the source for basalts was Indian-Ocean-type mantle. Based on the rare earth element (REE) data and small isotopic range, all the rocks from the Band-e-Zeyarat/Dar Anar ophiolite are cogenetic and were derived by fractionation from melts with a composition similar to average E-MORB; fractionation was controlled by the removal of clinopyroxene, hornblende and plagioclase. Three ⁴⁰Ar–³⁹Ar plateau ages of 140.7±2.2, 142.9±3.5 and 141.7±1.0 Ma, and five previously published K–Ar ages ranging from 121±4 to 146±5 Ma for the hornblende gabbros suggest that rocks from this ophiolite were formed during the Late Jurassic–Early Cretaceous. Plate reconstructions suggest that the rocks of this complex appear to be approximately contemporaneous with the Masirah ophiolite which has crystallization age of (150 Ma). Like Masirah, the rocks from the Band-e-Zeyarat/Dar Anar ophiolite complex represent southern Tethyan ocean crust that was formed distinctly earlier than crust preserved in the 90–100 Ma Bitlis-Zagros ophiolites (including the Samail ophiolite).     

Makran ophiolitic basalts (SE Iran) record Late Cretaceous Neotethys plume-ridge interaction

ABSTRACT

The Makran complex in southeast Iran provides a spectacular subduction-related accretionary complex to understand the mechanism of oceanic accretion and the evolution of subduction zones. In this paper, we present new major and trace element data as well as isotopic compositions of mafic volcanic blocks from the Makran ophiolitic mélange complex (OMC). Our aim is to assess the genesis of these rocks and discuss their implications on the evolution of Neotethys Ocean. These volcanic blocks are composed mainly of basalts with minor trachytes. The Makran lavas are occasionally interlayered with tuff layers. Zircons from these tuffs give U-Pb ages of 95 Ma, which is well in accordance with the reposted microfossil data for the interlayered pelagic limestones with pillow lavas. Makran basalts can be geochemically subdivided into four groups; normal to transitional MORB, enriched-MORB, Plume-type MORB and alkaline (-OIB-like) basalts. The OIB-like pillow lavas are represented by high values of Th/Tb (6.3–7.4) which are higher than other basalts (group 1 = 0.3–0.8; groups 2 = 0.7–1.6; group 3 = 1.58–1.36).¹⁴³Nd/¹⁴⁴Nd_(t) ratios for basalts ranges from 0.51247 to 0.51292, whereas ⁸⁷Sr/⁸⁶Sr_(t) isotopic composition of the OMC lavas varies from 0.704433 to 0.709466. The Pb isotopic composition of the lavas are quite high, ranging from 15.49–15.66 for ²⁰⁷Pb/²⁰⁴Pb_(t), 18.09–19.12 for ²⁰⁶Pb/²⁰⁴Pb_(t) and 37.80–39.23 for ²⁰⁸Pb/²⁰⁴Pb_(t). The chemistry of these rocks suggests that they were formed most likely in an oceanic setting with clear plume-ridge interaction. These rocks can form from partial melting of a highly heterogeneous mantle source, which is extensively metasomatized with deep mantle OIB-type components. We suggest these rocks have been generated in an oceanic ridge with plume-ridge interaction, similar to the Iceland-Reykjanes Ridge, before being fragmented and accreted into the Makran accretionary complex.