PGE and isotopic characteristics of Shergol and Suru Valley Ophiolites,Western Ladakh: Implications for supra-subduction tectonics along Indus Suture Zone

Present study reports the PGE-geochemistry of mantle peridotites and Nd-isotope geochemistry of arc related mafic rocks from the Indus Suture Zone (ISZ), western Ladakh. The total PGE concentration of the Shergol and Suru Valley peridotites (∑PGE = 96–180 ppb) is much higher than that of the primitive mantle and global ophiolitic mantle peridotites. The studied peridotites show concave upward PGE-distribution patterns with higher palladium-group PGE/Iridium-group PGE ratios (i.e., 0.8–2.9) suggesting that the partial melting is not the sole factor responsible for the evolution of these peridotites. The observed PGE-distribution patterns are distinct from residual/refractory mantle peridotites, which have concave downward or flat PGE-distribution patterns. Relative enrichment of palladium-group PGE as well as other whole-rock incompatible elements (e.g., LILE and LREE) and higher Pd/Ir ratio (1.1–5.9) reflects that these peridotites have experienced fluid/melt interaction in a supra-subduction zone (SSZ) tectonic setting. Also, the Shergol mafic intrusives and Dras mafic volcanics, associated with the studied peridotites, have high ¹⁴³Nd/¹⁴⁴Nd ratios (i.e., 0.512908–0.513078 and 0.512901–0.512977, respectively) and positive ε_Nd(t) (calculated for t = 140 Ma) values (i.e., +5.3 to +8.6 and + 5.1 to +6.6, respectively), indicating derivation from depleted mantle sources within an intra-oceanic arc setting, similar to Spongtang and Nidar ophiolites from other parts of Ladakh Himalaya. The transition from SSZ-type Shergol and Suru Valley peridotites to Early Cretaceous tholeiitic Shergol mafic intrusives followed by tholeiitic to calc-alkaline Dras mafic volcanics within the Neo-Tethys Ocean exhibit characteristics of subduction initiation mechanism analogous to the Izu-Bonin-Mariana arc system within western Pacific.     

Multilevel Graph Partitioning for Three-Dimensional Discrete Fracture Network Flow Simulations

Mathematical Geosciences - We present a topology-based method for mesh-partitioning in three-dimensional discrete fracture network (DFN) simulations that takes advantage of the intrinsic...     

Contamination of sediments in the floodplain wetlands of the lower uMngeni River,Kwa-Zulu Natal,South Africa

This study describes the geochemical and physico-chemical characterization of contaminant status of six wetlands of the lower uMngeni River floodplain, KwaZulu-Natal Province, South Africa. At all sites, organic matter content, calcium carbonate and metal concentrations were highest within the fine-sediment dominated surface laminae, decreasing with depth. Exceptions were displayed by arsenic (As), nickel (Ni) and vanadium (V), presumed to be derived from normal geochemical processes. Geochemical indicies indicated high contamination factors for lead (Pb) and V, pointing to moderate anthropogenic metal pollution contribution. Most conspicuous are the contamination high factors for Pb across most sites. Fe and Pb contamination is ubiquitous in all surface laminae, with contamination factors being either border-line, moderate or of high contamination. Enrichment factors for As, Mn and P at site S6 indicate possible anthropogenic sources. The data also shows a possibility of anthropogenic input of P and Zn at site S3. Apart from the enrichment factors established for Ni across all sites, all other metals indicate some degree of enrichment. A range of variable results for the mid to deeper laminae are indicative of natural processes with some human influence. Results of ANOVA confirm the above. Principal component analysis reveals a 46.4% variance from component 1 elements (Al, Ca, Cr, Cu, Fe, Mg, Mn and Zn) and a 17.8% variance from component 2 (Al, Cu and Fe). Given the rapid rate of development in this region, the protection of these floodplain wetlands must receive high conservation priority from the local municipality.     

Robust system size reduction of discrete fracture networks: a multi-fidelity method that preserves transport characteristics

We propose a multi-fidelity system reduction technique that uses weighted graphs paired with three-dimensional discrete fracture network (DFN) modelling for efficient simulation of subsurface flow and transport in fractured media. DFN models are used to simulate flow and transport in subsurface fractured rock with low-permeability. One method to alleviate the heavy computational overhead associated with these simulations is to reduce the size of the DFN using a graph representation of it to identify the primary flow sub-network and only simulate flow and transport thereon. The first of these methods used unweighted graphs constructed solely on DFN topology and could be used for accurate predictions of first-passage times. However, these techniques perform poorly when predicting later stages of the mass breakthrough. We utilize a weighted-graph representation of the DFN where edge weights are based on hydrological parameters in the DFN that allows us to exploit the kinematic quantities derivable a posteriori from the flow solution obtained on the graph representation of the DFN to perform system reduction and predict the later stages of the breakthrough curve with high fidelity. We also propose and demonstrate the use of an adaptive pruning algorithm with error control that produces a pruned DFN sub-network whose predicted mass breakthrough agrees with the original DFN within a user-specified tolerance. The method allows for the level of accuracy to be a user-controlled parameter.     

Molecular composition and paleobotanical origin of Eocene resin from northeast India

The molecular composition of fossil resins from early to middle Eocene coal from northeast India, has been analyzed for the first time to infer their paleobotanical source. The soluble component of fossil resin was analyzed using gas chromatography–mass spectrometry (GC–MS). The resin extracts are composed of cadalene-based C₁₅ sesquiterpenoids and diagenetically altered triterpenoids. The macromolecular composition was investigated using pyrolysis gas chromatography-mass spectrometry (Py-GC–MS) and Fourier transform infrared (FTIR) spectroscopy. The major pyrolysis products are C₁₅ bicyclic sesquiterpenoids, alkylated naphthalenes, benzenes and a series of C₁₇–C₃₄ n-alkene-n-alkane pairs. Spectroscopic analysis revealed the dominance of aliphatic components. The presence of cadalene-based sequiterpenoids confirms the resin to be Class II or dammar resin, derived from angiosperms of Dipterocarpaceae family. These sesquiterpenoids are often detected in many SE Asian fluvio-deltaic oils. Dipterocarpaceae are characteristic of warm tropical climate suggesting the prevalence of such climate during early Eocene in northeast India.     

Pointing of HAGAR telescope mirrors

An array of seven atmospheric Cherenkov telescopes was commissioned at a high altitude site in Hanle in the Ladakh region of the Himalayas. The array called HAGAR has been designed to observe celestial γ-rays of energy >100 GeV. Each telescope is altitude-azimuth mounted and carries seven parabolic mirrors whose optic axes are co-aligned with the telescope axis. The telescopes point and track a celestial source using a PC-based drive control system. Two important issues in positioning of each HAGAR telescope are pointing accuracy of telescope axis and co-alignment of mirrors’ optic axes with the telescope axis. We have adopted a three pronged strategy to address these issues, namely use of pointing models to improve pointing accuracy of the telescopes, RA-DEC scan technique to measure the pointing offsets of the mirrors and mechanical fine-tuning of off-axis mirrors by sighting a distant stationary light source. This paper discusses our efforts in this regard as well as the current status of pointing and monitoring of HAGAR telescopes.     

Convolution-based particle tracking method for transient flow

A convolution-based particle tracking (CBPT) method was recently developed for calculating solute concentrations (Robinson et al., Comput Geosci 14(4): 779–792, 2010). This method is highly efficient but limited to steady-state flow conditions. Here, we present an extension of this method to transient flow conditions. This extension requires a single-particle tracking process model run, with a pulse of particles introduced at a sequence of times for each source location. The number and interval of particle releases depends upon the transients in the flow. Numerical convolution of particle paths obtained at each release time and location with a time-varying source term is performed to yield the shape of the plume. Many factors controlling transport such as variation in source terms, radioactive decay, and in some cases linear processes such as sorption and diffusion into dead-end pores can be simulated in the convolution step for Monte Carlo-based analysis of transport uncertainty. We demonstrate the efficiency of the transient CBPT method, by showing that it requires fewer particles than traditional random walk particle tracking methods to achieve the same levels of accuracy, especially as the source term increases in duration or is uncertain. Since flow calculations under transient conditions are often very expensive, this is a computationally efficient yet accurate method.     

Delineating floodplain and upland areas for hydrologic models: a comparison of methods

A spatially distributed representation of basin hydrology and transport processes in hydrologic models facilitates the identification of critical source areas and the placement of management and conservation measures. Floodplains are critical landscape features that differ from neighbouring uplands in terms of their hydrological processes and functions. Accordingly, an important step in watershed modelling is the representation of floodplain and upland areas within a watershed. The aim of this study is (1) to evaluate four floodplain–upland delineation methods that use readily available topographic data (topographic wetness index, slope position, uniform flood stage, and variable flood stage) with regard to their suitability for hydrological models and (2) to introduce an evaluation scheme for the delineated landscape units. The methods are tested in three U.S. watersheds ranging in size from 334 to 629 km² with different climatic, hydrological, and geomorphological characteristics. Evaluation of the landscape delineation methods includes visual comparisons, error matrices (i.e. cross‐tabulations of delineated vs reference data), and geometric accuracy metrics. Reference data were obtained from the Soil Survey Geographic (SSURGO) database and Federal Emergency Management Agency (FEMA) flood maps. Results suggest that the slope position and the variable flood stage method work very well in all three watersheds. Overall percentages of floodplain and upland areas allocated correctly were obtained by comparing delineated and reference data. Values range from 83 to 93% for the slope position and from 80 to 95% for the variable flood stage method. Future studies will incorporate these two floodplain–upland delineation methods into the subwatershed‐based hydrologic model Soil and Water Assessment Tool (SWAT) to improve the representation of hydrological processes within floodplain and upland areas. Copyright © 2016 John Wiley & Sons, Ltd.     

Monsoon circulation interaction with Western Ghats orography under changing climate

In this study, the authors have investigated the likely future changes in the summer monsoon over the Western Ghats (WG) orographic region of India in response to global warming, using time-slice simulations of an ultra high-resolution global climate model and climate datasets of recent past. The model with approximately 20-km mesh horizontal resolution resolves orographic features on finer spatial scales leading to a quasi-realistic simulation of the spatial distribution of the present-day summer monsoon rainfall over India and trends in monsoon rainfall over the west coast of India. As a result, a higher degree of confidence appears to emerge in many aspects of the 20-km model simulation, and therefore, we can have better confidence in the validity of the model prediction of future changes in the climate over WG mountains. Our analysis suggests that the summer mean rainfall and the vertical velocities over the orographic regions of Western Ghats have significantly weakened during the recent past and the model simulates these features realistically in the present-day climate simulation. Under future climate scenario, by the end of the twenty-first century, the model projects reduced orographic precipitation over the narrow Western Ghats south of 16°N that is found to be associated with drastic reduction in the southwesterly winds and moisture transport into the region, weakening of the summer mean meridional circulation and diminished vertical velocities. We show that this is due to larger upper tropospheric warming relative to the surface and lower levels, which decreases the lapse rate causing an increase in vertical moist static stability (which in turn inhibits vertical ascent) in response to global warming. Increased stability that weakens vertical velocities leads to reduction in large-scale precipitation which is found to be the major contributor to summer mean rainfall over WG orographic region. This is further corroborated by a significant decrease in the frequency of moderate-to-heavy rainfall days over WG which is a typical manifestation of the decrease in large-scale precipitation over this region. Thus, the drastic reduction of vertical ascent and weakening of circulation due to ??upper tropospheric warming effect?? predominates over the ??moisture build-up effect?? in reducing the rainfall over this narrow orographic region. This analysis illustrates that monsoon rainfall over mountainous regions is strongly controlled by processes and parameterized physics which need to be resolved with adequately high resolution for accurate assessment of local and regional-scale climate change.     

Archaean crustal thickness of greenstone granite belts of South India

Archaean crustal thickness for the Dharwar craton is estimated using potash index and Rb?Sr crustal thickness grid. The volcanics of the Dharwar greenstone belts appear to have evolved in a less than 20 km thick crust. Whereas the tonalite-trondhjemite pebbles of the Dharwar conglomerates (3250±150 m.y.) were derived from gneisses that evolved in a crust less than 20 km thick, the bulk of the peninsular gneisses and associated granitoids were emplaced in a crust 25 to 35 km thick. The 2000 m.y. old Closepet granite suite was emplaced in a crust thicker than 30 km. It is deduced that the continental crust in the region thickened from 15 to 35 km during a span of about 1000 m.y. between 3250±150 to 2000 m.y. ago. Calculations show that Archaean gecthermal gradients in Dharwar craton were three to four times steeper when compared to the present 10.5°C/km. The thin crust and the steep geothermal gradients are reflected by the emplacement of high magnesia basalts, layered igneous complexes and the strong iron enrichment trend shown by Dharwar metavolcanics.