Petrography and Geochemistry of Dolomites of Samanasuk Formation, Dara Adam Khel Section, Kohat Ranges, Pakistan

Replacement dolomite occurs in Jurassic Samanasuk Formation in Dara Adam khel area of Kohat ranges, North-Western Himalayas, Pakistan. This study, for the first time, document the process of dolomitization and evolution of strata bound dolomitic bodies. Field investigation, petrography and geochemistry helped in unraveling the formation of several dolomitic bodies. Petrographically dolomites comprises of: (1) medium grain crystalline planer subhedral dolomite (Dol-I); (2) fine grained crystalline anhedral non-planer dolomite rhombs (Dol-II); (3) medium to coarse grained crystalline subhedral-anhedral non-planer dolomite (Dol-III) and coarse to very coarse grained crystalline saddle dolomite cements (SD). The saddle dolomites (SD) postdate the replacement dolomites and precede telogenetic calcite (TC) cements. Stable O and C isotope analysis shows that these dolomites have δ18Ovpdb ranging from -4.09% to -10.4 whereas the δ13Cvpdb ranges from +0.8 to +2.51. Major and trace elements data show that Sr concentrations of 145.5 to 173 ppm; Fe contents of 2198 to 8215 ppm; and Mn contents of 93.5 to 411 ppm. Petrographically replacive dolomites, saddle dolomite, and δ18Ovpdb values depicts neomorphism of replacement dolomites that were formed earlier were exposed to late dolomitizing fluids. As a result of basin uplift during the Himalayan orogeny in Eocene time, dolomitization event was stopped through occurrence of meteoric water. The Main Boundary Thrust (MBT) and its splays were most likely essential conduits that channelized dolomitizing fluids from siliciclastic rocks that were buried deeply into the Jurassic carbonates rocks, leading to more extreme dolomitization.

First record of stable isotopes (δ13C, δ18O) and element ratios (Mg/Ca, Sr/Ca) of Middle to Late Jurassic belemnites from the Indian Himalayas and their potential for palaeoenvironmental reconstructions

Middle to Late Jurassic belemnites from the Spiti and Zanskar valleys in the Indian Himalayas were used for stable isotope (δ13C, δ18O) and element (Mg/Ca, Sr/Ca) analyses. Although the Himalayan orogeny deformed and altered a large portion of the collected fossils, cathodoluminescence and scanning electron microscopy in combination with analyses of iron and manganese contents allowed the identification of belemnites believed to still retain their original chemical composition. Results indicate a long-term temperature decrease from the Middle Callovian–Oxfordian to the Tithonian, which is proposed to have been caused by a concomitant drift of eastern Gondwana into higher palaeolatitudes. Reconstructed absolute temperatures depend on the used equation and assumed δ18O value of seawater, but most likely varied between 17.6 °C to 27.6 °C in the Kimmeridgian and Tithonian with average values between 22 °C to 24 °C. This way, temperatures were similar to slightly warmer than today at comparable latitudes. The reconstruction of absolute temperatures for the Middle Callovian–Oxfordian was hindered by a larger number of poorly preserved belemnites representing this time interval.

Where is the missing felsic crust of Greater Indian continent?

According to the plate tectonics theory, continental crust (CC), especially the felsic upper and middle continental crust (UCC and MCC), cannot subduct due to its buoyancy. Therefore, most, if not all, of the felsic crustal mass will be preserved in continental collision zones or eroded by the surface process. Consequently, the continent-continent convergence is generally slower and more short-lived than oceanic plate subduction. However, the long-duration, fast convergence, and imbalance of crustal mass in the India-Asia collisional system challenge the classical rules of plate tectonics. Systematic compilation and calculations indicate ~20-47% of the felsic crust in Greater India is missing during collision. Based on the phase equilibria modeling and density calculations, we explore the pressure-temperature-dependent density evolutions of UCC and MCC and demonstrate they are denser than the surrounding mantle at P >7-8 GPa when the phase transition from coesite to stishovite occurs. The phase equilibria induced density evolution is further integrated into the thermo-mechanical model, which confirm the deep subduction of Greater Indian continent with its felsic UCC and MCC. Analytical studies of the slab-pull forces in the subduction zone indicate the Greater Indian continent can subduct spontaneously under its own negative buoyancy when it is dragged to a depth of ~170 km by the preceding oceanic slab. The great slab-pull force, induced by the negative buoyancy of UCC and MCC below 170 km, not only contributes to the long-lasting fast convergence between India and Asia, but also explains the crustal mass imbalance during the Himalayan orogeny.

Filling Provenance in Fracture Cavity Formation within Aksu Area, Tarim Basin, NW China: Indicators from Major and Trace Element, Carbon-Oxygen, and Strontium Isotope Compositions

With an aim to increase the understanding about sedimentary environment and isotopic and chemical characteristics of fillings in fracture cavities with multiple compositions, we conducted scanning electron microscope (SEM), fluid inclusion testing (FIT), common and trace element chemistry, full analysis testing, isotopic compositions (δ13C, δ18O, 87Sr/86Sr), and apatite fission track testing to study the formation environment of Aksu area, Tarim Basin. According to outfield and microscope observations, combined with SEM results, three textural and compositional type fractures and cavities were distinguished. Through fine analysis of geochemistry characteristics on fractures, cavities, multiple filling periods, and environments were interpreted. Constrained by rare earth element (REE) pattern diagram, relationships between carbon and oxygen isotopes, strontium isotope, the compositional patterns, and generation environment of the fracture and cavities were determined. The results show that (1) cavity, fracture filling, and wall rock primarily consist of calcite, with a proportion of 56.85%, 80.48%, and 81.00%, respectively. (2) Four fracture sets have been distinguished in the Ordovician limestone of the karst cave, Middle-Late Caledonian (Set 1), Early Hercynian (Set 2), Indo-Yanshanian (Set 3), and Himalayan orogeny (Set 4). Two stages of cave filling deposition are distinguished. Stage I was coeval with the Middle-Late Caledonian Set 1 fractures and is attributable to the circulation of freshwater fluid. Stage II was coeval with the Early Hercynian Set 2 fractures and is attributable to deep hydrothermal fluid circulation. (3) Cavity, fracture filling, and wall rock in Ordovician strata are slightly influenced by diagenesis alteration and territorial supply. Three significant filling stages were distinguished, freshwater fluid with strong oxidizing environment (Middle-Late Caledonian), hydrothermal fluid with authigenic abnormal enrichment (indicating obvious hypoxic sedimentary water, Early Hercynian), and high-temperature hydrothermal fluid from deep earth (primarily influenced by magmatism, Indo-Yanshanian, and Himalayan).

A Study on the Heterogeneity Characteristics of Geological Controls on Coalbed Methane Accumulation in Gujiao Coalbed Methane Field, Xishan Coalfield, China

Commercial exploration and exploitation of coalbed methane (CBM) in Gujiao coalbed methane (CBM) field, Xishan coalfield, have rapidly increased in recent decades. The Gujiao CBM field has shown strong gas distribution heterogeneity, low gas content, and wide distribution of wells with low production. To better understand the geological controlling mechanism on gas distribution heterogeneity, the coal reservoir evolution history and CBM accumulation process have been studied on the base of numerical simulation work. The burial history of coal reservoir can be classified into six stages: shallowly buried stage; deeply buried stage; uplifting stage; short-term tectonic subsidence stage; large-scale uplifting stage; and sustaining uplifting and structural inversion stage. Mostly, coal seams have experienced two-time thermal metamorphisms with twice hydrocarbon-generation processes in this area, whereas in the southwest part, the coal seams in there suffered three-time thermal metamorphisms and hydrocarbon-generation processes. The critical tectonic events of the Indosinian, Yanshanian, and Himalayan orogenies affect different stages of the CBM reservoir accumulation evolution process. The Indosinian orogeny mainly controls the primary CBM generation. The Yanshanian orogeny dominates the second and third gas generation and migration processes. The Himalayan orogeny mainly affects the gas dissipation process and current CBM distribution heterogeneity.

Investigating the Response of Tectonic Loading in the Himalayas on the Peripheral Forebulge through Drainage Network Analysis

<p>Drainage divide migration is a conspicuous natural process through which a landscape evolves. In response to a forced climatic and tectonic disturbance, susceptible river networks transfer the transient signals to the entire river basin, which results in an incision or aggradation. The Himalayan orogeny and subduction of the Indian plate have resulted in an upward flexure in the Indian lithosphere known as a peripheral forebulge. A forebulge can flexurally uplift and migrate following the variation in tectonic load. The emergence of the central Indian plateau is a consequence of the upwarping of the Indian lithosphere (Bilham et al. 2003).  In this work, we are trying to assess the drainage network dynamics between the Narmada and Ganga river systems, which drain the uplifted central Indian plateau. We have calculated the Chi(χ) metrics, steepness index (Ksn), knickpoints for the channels in the study area. We have generated Topographic swath profiles to analyze the topographic variations on the plateau. It has been observed from the results that the rivers in the study area lack dynamic equilibrium, and river capturing is an evident response to the perturbations. Our analysis shows that the Narmada River tributaries are gaining drainage area and aggressing Northwards by capturing adjacent Ganga river tributaries. The field observations show a variation in the surface slope and presence of knickpoints (waterfalls) along the "aggressor" drainages. We propose a model to show a correlation between the tectonic loading of Himalayas, movement of forebulge, and its feedback to the river systems present on the forebulge.</p>

Paleomagnetic results from the western Himalaya indicate multi-stage India-Eurasia collision

<p>The classical model for the collision between India and Eurasia, which resulted in the formation of the Himalayan orogeny, is a single-stage continent-continent collision event at around 55 – 50 Ma. However, it has also been proposed that the India-Eurasia collision was a multi-stage process involving an intra-oceanic Trans-Tethyan subduction zone south of the Eurasian margin. We present paleomagnetic data constraining the location the Kohistan-Ladakh arc, a remnant of this intra-oceanic subduction zone, to a paleolatitude of 8.1 ± 5.6 °N between 66 – 62 Ma. Comparing this result with new paleomagnetic data from the Eurasian Karakoram terrane, and previous paleomagnetic reconstructions of the Lhasa terrane reveals that the Trans-Tethyan Subduction zone was situated 600 – 2,300 km south of the contemporaneous Eurasian margin at the same time as the first ophiolite obduction event onto the northern Indian margin. Our results confirm that the collision was a multistage process involving at least two subduction systems. Collision began with docking between India and the Trans-Tethyan subduction zone in the Late Cretaceous and Early Paleocene, followed by the India-Eurasia collision in the mid-Eocene. The final stage of India-Eurasia collision occurred along the Shyok-Tsangpo suture zone, rather than the Indus-Tsangpo. The addition of the Kshiroda oceanic plate, north of India after the Paleocene reconciles the amount of convergence between India and Eurasia with the observed shortening across the India–Eurasia collision system. Our results constrain the total post-collisional convergence accommodated by crustal deformation in the Himalaya to 1,350 – 2,150 km, and the north-south extent of the northwestern part of Greater India to < 900 km.</p>

Crustal structure across Himalayan Frontal Thrust (HFT) using high-resolution seismic datasets

<p>The Himalayan fold-thrust belt has been developing due to the northward convergence of the Indian plate against the Eurasian plate since ~55 Ma. Three major thrust systems: Main Central Thrust (MCT), Main Boundary Thrust (MBT), and Himalayan Frontal Thrust (HFT) are distinctly observed in the Himalayan orogeny from north to south indicating southward propagation of active deformation. These active thrust systems produced several devastating earthquakes in the past such as 1905 Kangra (Mw 7.8), 1934 Nepal-Bihar (Mw 8), and 1950 Assam (Mw 8.6) earthquakes. Presently HFT is found to be the tectonically very active zone that accommodates a strain rate of ~10-15 mm/year and is a zone for great threats in near future to the societies residing over the Himalayan foothills. The present study carried out in the lower Siwalik Himalaya near Pawalgarh in Nainital District of Uttarakhand, India with an objective to estimate the velocity model across HFT in the locality. To accomplish the objective, seismic data were acquired along three profiles of a cumulative length of ~13 km using a seismic thumper as a source and 96 vertical component geophones with the natural frequency of 5 Hz and Remote Acquisition Unites (RAUs) as sensors and data loggers, respectively, and with a group and shot interval of 20 m and near offset of 100 m. Highly uneven Himalayan terrain causes large static errors. In order to overcome this challenge, we used Real Time Kinematics (RTK) to estimate more precise source and receiver surface elevation. In the pre-processing phase of acquired seismic data, three different shots taken at the same location are vertically stacked to eliminate random non-coherent noises and improve the SNR of the data. We then applied a low-frequency array filter (LFAF) to suppress the ground roll using velocity estimates from the ambient noise tomography (ANT). We process the data by implementing conventional seismic processing techniques including normal move-out (NMO) correction, velocity analysis followed by stacking. In the stack section, we observe a northward dipping reflector extending from the surface to ~ 1- 1.25 s TWT indicating evidence of HFT. Another reflector observed at ~3-4 s TWT demarcating the extent of overlying sedimentary deposits on the top of the under-thrusting lithosphere. Rocks of the Siwalik Himalaya mainly composed of sedimentary deposits of sandstone mudstone, and alluvial deposits. Average velocity obtained from the refraction tomography ~ 2900 m/s matches well with rock type in the region. Thus, the high-resolution crustal structure across the highly active HFT can be crucial to understand the earthquake mechanism in the locality and for a better hazard assessment.</p>

Contemporary Earthquake Hazards in the West-Northwest Himalaya: A Statistical Perspective through Natural Times

Himalayan earthquakes have deep societal and economic impact. In this article, we implement a surrogate method of nowcasting (Rundle et al., 2016) to determine the current state of seismic hazard from large earthquakes in a dozen populous cities from India and Pakistan that belong to the west-northwest part of Himalayan orogeny. For this, we (1) perform statistical inference of natural times, intersperse counts of small-magnitude events between pairs of succeeding large events, based on a set of eight probability distributions; (2) compute earthquake potential score (EPS) of 14 cities from the best-fit cumulative distribution of natural times; and (3) carry out a sensitivity testing of parameters—threshold magnitude and area of city region. Formulation of natural time (Varostos et al., 2005) based on frequency–magnitude power-law statistics essentially avoids the daunting need of seismicity declustering in hazard estimation. A retrospective analysis of natural time counts corresponding to M≥6 events for the Indian cities provides an EPS (%) as New Delhi (56), Chandigarh (86), Dehradun (83), Jammu (99), Ludhiana (89), Moradabad (84), and Shimla (87), whereas the cities in Pakistan observe an EPS (%) as Islamabad (99), Faisalabad (88), Gujranwala (99), Lahore (89), Multan (98), Peshawar (38), and Rawalpindi (99). The estimated nowcast values that range from 38% to as high as 99% lead to a rapid yet useful ranking of cities in terms of their present progression to the regional earthquake cycle of magnitude ≥6.0 events. The analysis inevitably encourages scientists and engineers from governments and industry to join hands for better policymaking toward land-use planning, insurance, and disaster preparation in the west-northwest part of active Himalayan belt.

A new Himalayan species of Alfieriella (Coleoptera: Cryptophagidae)

A new species of Alfieriella Wittmer, 1935 (Coleoptera, Cryptophagidae), Alfieriella senguptai sp. n. from China and India, is described. This is the first formal record of the genus Alfieriella and the tribe Hypocoprini from the Himalayan region. Alfieriella senguptai is the largest member of Alfieriella, and its presence in a cold, high-altitude environment conforms to Bergmann’s rule. The distribution of the genus Alfieriella may be associated with the breakup of the Tethys Ocean and the origin of A. senguptai influenced by the Himalayan orogeny. A distribution map and a key to species of Alfieriella are also provided.