Isotopic Fingerprints of Paleoclimates during the Last 30,000 Years in Deep Confined Groundwaters of Southern India

Isotopic and geochemical evidence of paleoclimates, especially for the last glaciation, has been obtained from deep confined groundwaters of southern India. The δ¹³C, δ¹⁸O, chloride, and deuterium analyses of groundwaters show distinct excursions inferred to be related to climatic variations. The arid climatic episode associated with the last glaciation (18,000 ± 2000 yr B.P.) is conspicuously identified by signatures of relatively enriched δ¹³C (−10 to −12‰ PDB) and δ¹⁸O (−5.3 to −4.8‰ SMOW) values, and high chloride concentration (80 to 160 mg/l). The transition from an arid to humid period ca. 12,000–8000 yr B.P. is shown by a decreasing trend in the δ¹³C (−9.5 to −17‰) and δ¹⁸O (−4.5 to −6.3‰) contents of groundwaters. The late Holocene (since 4000 yr B.P.), marked by a more humid but unstable climate, is identified by further depletion of δ¹³C (−13 to −20‰) and δ¹⁸O (−5.2 to −6.3‰). Similar variation between δ¹⁸O and chloride values in confined groundwaters further demonstrates two distinct climatic excursions (arid and humid) governed by the “amount effect.” This is the first time that isotopic and geochemical signatures related to changing paleoclimates have been identified in the confined groundwaters of the southern Indian landmass.     

Optical dating of liquefaction features to constrain prehistoric earthquakes in Upper Assam, NE India—some preliminary results

Optical dating of liquefied sand structures formed during major earthquakes in Upper Assam, northeast of India, has been carried out to constrain the timing of prehistoric earthquakes in this seismically active region. The bleaching of source material and of the same material during the creation of a liquefaction feature was tested using quartz extracted from 21 samples associated with two different liquefaction structures in Upper Assam. Due to the poor sensitivity of quartz from this region, a sensitisation procedure was used in our SAR protocol to reduce the scatter in optical ages. Various internal consistency tests of the measurement protocol and the excellent reproducibility of the OSL ages indicate that the dose estimates from the quartz are accurate and the optical ages reliable. The preliminary OSL ages indicate that the liquefaction features were formed between 1430 AD and 1630 AD. This study demonstrates that using OSL, ‘direct dating’ of prehistoric earthquakes may be possible, if sand blows from liquefied dykes are preserved.     

Stable isotope systematics and geochemical signatures constraining groundwater hydraulics in the mining environment of the Korba Coalfield,Central India

The present study explores the use of stable isotopes and major ion chemistry of various water sources such as mine water, groundwater, river water, and rainwater to identify the dominant hydrogeochemical process that controls the water quality in the active opencast coal mining area of the Korba Coalfield, India. Different hydrochemical models have revealed that the study area is mostly characterized by two facies alongside a dominance by ion exchange: i.e., a slightly mineralized Ca²⁺–Mg²⁺–Cl^? and a meteorically derived Ca²⁺–HCO₃^? fresh water facies. In the isotope bivariate plot, the samples plot to the right of the LMWL and slopes observed in pre (4.94) and post-monsoon (5.85) seasons are not as steep as the LMWL (7.95) and this warrants the suggestion that meteoric water is the major source that replenishes the dynamic groundwater resource in the study area after being subjected to evaporation. The negative d-excess values (<?0) and enrichment of the δ¹⁸O ratio observed at some locations suggests a non-equilibrium process and a “mechanism of mixing”. A noteworthy fact is that isotopically river water is in tandem with the mine water samples for both seasons and testifies to a hydraulic connectivity between Hasdeo River water and mine water through a major fault. It is substantiated by mine inflows assessed by the sump test for Kusmunda and Gevra, which is two times higher than the predicted values by Darcy’s law. The wide variation in the isotopic composition is attributed to different vapor sources viz., southwest monsoon (SW) that originates in the Arabian Sea and locally driven air moisture from surface water bodies.     

A study on the Wilkins and Forchheimer equations used in coarse granular media flow

Complexity of the pore geometry and the random nature of flow velocity make it difficult to predict and represent post laminar flow through porous media. Present study experimentally investigates the applicability of Forchheimer and Wilkins equations for post laminar flow where Darcy’s law is invalid due to predominant inertial effect. It is observed that both porosity and media size have significant influence over the coefficients of the Forchheimer coefficients. To incorporate the effect of porosity and media size, behaviour of Forchheimer coefficients are investigated with hydraulic radius as characteristic length. An inversely proportional variation trend is found for all the present and earlier reported data. A new empirical relation between Forchheimer coefficients and hydraulic radius is obtained which can be universally applicable for all media size and porosity. Coefficients of the Wilkins equation are found to be non-deviating for different hydraulic radius in the present study and in the reported literature validating its applicability in predicting the non laminar flow through porous media. Further the Wilkins equation is modified after incorporating the correction factors for better applicability on the field.     

Tomographic approach for gas hydrate investigation in Kerala-Konkan region,India

Seismic tomography is an effective means of estimating velocity and structure from multichannel seismic (MCS) reflection data. In this study we have followed a 2D approach to arrive at the probable velocity field configuration from multichannel seismic data and infer the presence of gas hydrates/free-gas in the offshore Kerala-Konkan region, along the eastern part of a seismic line on which a bottom simulating reflector (BSR) has previously been identified. Tomographic modeling consists of the identification of reflection phases and picking of respective travel times for various source-receiver positions. These picks were then utilized to arrive at a 2D velocity field following a forward and inversion approach using a ray tracing technique. The modeling for the first time brought out the finer scale velocity structure under the region of investigation. Modeling through the 2D approach shows lateral variation in velocity field along the studied segment of the seismic line. The results indicate a thin (∼50–60 m) sedimentary cover with velocity ranging from 1,770 to 1,850 m/s. A sedimentary layer with high P-wave velocity 1,980–2,100 m/s below the sea floor was interpreted as the hydrate layer. The thickness of this layer varies between 110 and 140 m. The hydrate layer is underlain by a low-velocity layer having velocities in the range 1,660–1,720 m/s. This low velocity may represent a free gas layer, whose thickness varies between 50 and 100 m located below the hydrated layer. The investigation suggests the occurrence of gas hydrate underlain by free gas in some parts of the Kerala-Konkan offshore region.     

Study of forward–backward solute dispersion profiles in a semi-infinite groundwater system

ABSTRACT

Forward–backward solute dispersion with an intermediate point source in one-dimensional semi-infinite homogeneous porous media is studied in this paper. Solute transport under sorption conditions, first-order decay and zero-order production terms are included. The first type of boundary condition is taken as a constant point source at an intermediate point from where forward and backward solute dispersion is examined. The Laplace transform method is adopted to solve the governing equation analytically. All the analytical results are obtained in graphical form to investigate the forward–backward solute transport in porous media for various hydrological input data. The graphical nature of the analytical solution is compared with numerical data taken from existing literature and similar results are obtained. Also, numerical solution of the governing equation is obtained by the Crank-Nicolson finite difference scheme and validated with the analytical solution, which demonstrates good agreement between them. Accuracy of the solution is also observed by using RMSE.