Origin and dynamics of groundwater salinity in the alluvial plains of western Delhi and adjacent territories of Haryana State,India

Groundwater salinity is a widespread problem and a challenge to water resources management. It is an increasing concern in the alluvial plains of Delhi and neighbouring Haryana state as well as a risk for agricultural production water supply and sustainable development. This study aims to identify potential sources of dissolved salts and the driving mechanisms of salinity ingress in the shallow aquifer. It combines a comprehensive review of environmental conditions and the analysis of groundwater samples from 25 sampling points. Major ions are analysed to describe the composition and distribution of saline groundwater and dissolution/precipitation dynamics. Density stratification and local upconing of saline waters were identified by multilevel monitoring and temperature logging. Bromide–chloride ratios hold information on the formation of saline waters, and nitrate is used as an indicator for anthropogenic influences. In addition, stable isotope analysis helps to identify evaporation and to better understand recharge processes and mixing dynamics in the study region. The results lead to the conclusion that surface water and groundwater influx into the poorly drained semiarid basin naturally results in the accumulation of salts in soil, sediments and groundwater. Human‐induced changes of environmental conditions, especially the implementation of traditional canal and modern groundwater irrigation, have augmented evapotranspiration and led to waterlogging in large areas. In addition, water‐level fluctuations and perturbation of the natural hydraulic equilibrium favour the mobilisation of salts from salt stores in the unsaturated zone and deeper aquifer sections. The holistic approach of this study demonstrates the importance of various salinity mechanisms and provides new insights into the interference of natural and anthropogenic influences. Copyright © 2011 John Wiley & Sons, Ltd.     

Encounters with an unearthly mudstone: Understanding the first mudstone found on Mars

The Sheepbed mudstone forms the base of the strata examined by the Curiosity rover in Gale Crater on Mars, and is the first bona fide mudstone known on another planet. From images and associated data, this contribution proposes a holistic interpretation of depositional regime, diagenesis and burial history. A lake basin probably received sediment pulses from alluvial fans. Bed cross‐sections show millimetre to centimetre‐scale layering due to distal pulses of fluvial sediment injections (fine‐grained hyperpycnites), fall‐out from river plumes, and some aeolian supply. Diagenetic features include mineralized synaeresis cracks and millimetre‐scale nodules, as well as stratiform cementation. Clay minerals were initially considered due to in situ alteration, but bulk rock chemistry and mineralogy suggests that sediments were derived from variably weathered source rocks that probably contained pre‐existing clay minerals. X‐ray diffraction analyses show contrasting clay mineralogy in closely spaced samples, consistent with at least partial detrital supply of clay minerals. A significant (ca 30 wt%) amorphous component is consistent with little post‐depositional alteration. Theoretical modelling of diagenetic reactions, as well as kinetic considerations, suggest that the bulk of diagenetic clay mineral formation occurred comparatively late in diagenesis. Diagenetic features (synaeresis cracks and nodules) were previously thought to reflect early diagenetic gas formation, but an alternative scenario of synaeresis crack formation via fabric collapse of flocculated clays appears more likely. The observed diagenetic features, such as solid nodules, hollow nodules, matrix cement and ‘raised ridges’ (synaeresis cracks) can be explained with progressive alteration of olivine/glass in conjunction with centrifugal and counter diffusion of reactive species. Anhydrite‐filled fractures in the Sheepbed mudstone occurred late in diagenesis when fluid pressures built up to exceed lithostatic pressure. Generating fluid overpressure by burial to facilitate hydraulic fracturing suggests a burial depth of at least 1000 m for the underlying strata that supplied these fluids.     

Modeling of source parameters of the 15 December 2015 Deogarh earthquake of M<Subscript>w</Subscript> 4.0

We herein present source parameters and focal mechanism of a rare cratonic upper crustal earthquake of M_w4.0, which occurred at 8 km depth (centroid depth) below a region near Deogarh, Jharkhand. For our study, we used broadband waveform data from a seismic network of 15 three-component seismographs in the eastern Indian craton. The average seismic moment, moment magnitude and source radius are estimated to be 1.1 × 10¹⁵ N-m, 4.0 and 180.6 m, respectively. The high average stress drop of 14.27 MPa could be attributed to its lower-crustal origin. The mean corner frequency is calculated to be 4.1 Hz. To study the source mechanism, we perform a deviatoric constrained full waveform moment tensor inversion of multiple point sources on the band-passed (0.06 – 0.14 Hz) broadband displacement data of the Deogarh event, using ISOLA software. The best fit is obtained for the source at 8 km centroid depth, with a moment magnitude 3.7, and a right-lateral strike-slip mechanism with strike 162°, dip 72° and rake 169°. The P-axis orients N24°E, which is parallel to the direction of the absolute plate motion direction of Indian plate, while T-axis orients E-W, which is parallel to the strike of the pre-existing Damodar Graben (DG) of Gondwana age. The occurrence of this earthquake is attributed to the neotectonic reactivation of a fault associated with the E-W trending DG shear zone.     

Simulating climate change impact on soil erosion using RUSLE model <Emphasis Type="BoldItalic">−</Emphasis> A case study in a watershed of mid-Himalayan landscape

Climate change, particularly due to the changed precipitation trend, can have a severe impact on soil erosion. The effect is more pronounced on the higher slopes of the Himalayan region. The goal of this study was to estimate the impact of climate change on soil erosion in a watershed of the Himalayan region using RUSLE model. The GCM (general circulation model) derived emission scenarios (HadCM3 A2a and B2a SRES) were used for climate projection. The statistical downscaling model (SDSM) was used to downscale the precipitation for three future periods, 2011–2040, 2041–2070, and 2071–2099, at large scale. Rainfall erosivity (R) was calculated for future periods using the SDSM downscaled precipitation data. ASTER digital elevation model (DEM) and Indian Remote Sensing data – IRS LISS IV satellite data were used to generate the spatial input parameters required by RUSLE model. A digital soil-landscape map was prepared to generate spatially distributed soil erodibility (K) factor map of the watershed. Topographic factors, slope length (L) and steepness (S) were derived from DEM. Normalised difference vegetation index (NDVI) derived from the satellite data was used to represent spatial variation vegetation density and condition under various land use/land cover. This variation was used to represent spatial vegetation cover factor. Analysis revealed that the average annual soil loss may increase by 28.38, 25.64 and 20.33% in the 2020s, 2050s and 2080s, respectively under A2 scenario, while under B2 scenario, it may increase by 27.06, 25.31 and 23.38% in the 2020s, 2050s and 2080s, respectively, from the base period (1985–2013). The study provides a comprehensive understanding of the possible future scenario of soil erosion in the mid-Himalaya for scientists and policy makers.     

Reconstruction of specific mass balance for glaciers in Western Himalaya using seasonal sensitivity characteristic(s)

Seasonal sensitivity characteristics (SSCs) were developed for Naradu, Shaune Garang, Gor Garang and Gara glaciers, Western Himalaya to quantify the changes in mean specific mass balance using monthly temperature and precipitation perturbations. The temperature sensitivities were observed high during summer (April–October) and precipitation sensitivities during winter months (November–March), respectively. The reconstructed mass balance correlates well with the field and remote sensing measurements, available between 1980 and 2014. Further, SSCs were used with the monthly mean temperatures and precipitation estimates of ERA 20CM ensemble climate reanalysis datasets to reconstruct the specific mass balance for a period of 110 years, between 1900 and 2010. Mass balance estimates suggest that the Shaune Garang, Gor-Garang and Gara glaciers have experienced both positive and negative mass balance, whereas the Naradu glacier has experienced only negative mass balance since 1900 AD. Further, a cumulative loss of \(-133 \pm 21.5\) m.w.e was estimated for four glaciers during the observation period. This study is the first record from Indian Himalaya in evaluating the mass balance characteristics over a century scale.     

Stable carbon and oxygen isotope study on benthic foraminifera: Implication for microhabitat preferences and interspecies correlation

Stable isotopes of benthic foraminifera have widely been applied in micropalaeontological research to understand vital effects in foraminifera. Isotopic fractionations are mainly controlled by ontogeny, bottom/pore water chemistry, habitat preference, kinetic effect and respiration. Discontinuous abundance of a species for isotopic analysis has forced us to select multiple species from down-core samples. Thus standardisation factors are required to convert isotopic values of one species with respect to other species. The present study is pursued on isotopic values of different pairs of benthic foraminifera from the Krishna–Godavari basin and Peru offshore to understand habitat-wise isotopic variation and estimation of isotopic correction factors for the paired species (Cibicides wuellerstorfi–Bulimina marginata, Ammonia spp.–Loxostomum amygdalaeformis and Bolivina seminuda–Nonionella auris). Infaunal species (B. marginata, Ammonia spp. and N. auris) show a lighter carbon isotopic excursion with respect to the epifaunal to shallow infaunal forms (C. wuellerstorfi, L. amygdalaeformis and B. seminuda). These lighter \(\updelta ^{13}\) \(\hbox {C}\) values are related to utilisation of \(\hbox {CO}_{2}\) produced by anaerobic remineralisation of organic matter. However, enrichment of \(\updelta ^{18}\) \(\hbox {O}\) for the deeper microhabitat (bearing lower pH and decreased \({\hbox {CO}_{3}}^{2-})\) is only recorded in case of B. marginata. It is reverse in case of N. auris and related to utilisation of respiratory \(\hbox {CO}_{2}\) and internal dissolve inorganic carbon pool. Estimation of interspecies isotopic correction factors for the species pairs (\(\updelta ^{13}\) \(\hbox {C}\) of C. wuellerstorfi–B. marginata, L. amygdalaeformis–Ammonia spp., N. auris–B. seminuda) and \(\updelta ^{18}\) \(\hbox {O}\) of C. wuellerstorfi–B. marginata are statistically reliable and may be used in palaeoecological studies.     

Scaling of response spectrum and duration for aftershocks

Aftershocks have the potential to cause collapse of a structure that has been already damaged by the preceding main shock. Seismic safety of a structure should therefore be ascertained through a damage analysis using the anticipated main shock and few larger-aftershock motions. Simulation of aftershock motions needs characterization of the seismic hazard due to aftershocks, and therefore it will be useful to develop a conditional scaling model that can predict the response spectrum of an anticipated aftershock motion consistent with the design spectrum of the main shock motion anticipated at the same station. In this study an attempt is made to develop a conditional scaling model for the pseudo spectral velocity spectrum via linear regression analysis on the aftershock and main shock recordings for the 1999 Chi–Chi earthquake. It is shown that it may be possible to obtain a simpler and approximate version of the conditional model from an unconditional model. Damage-causing potential of a ground motion also depends on its strong motion duration (SMD) and therefore a conditional scaling model is developed for SMD of the aftershock motion in several narrow frequency-bands. The model is developed for the larger-aftershock motions and it is shown that a reasonable replacement of such a model may be obtainable directly from an unconditional model. Finally, a simple weighted averaging scheme is proposed to obtain the composite SMD from the SMDs for different frequency bands by using the pseudo spectral acceleration spectrum of the motion.     

How Bayesian data assimilation can be used to estimate the mathematical structure of a model

In previous work, we presented a method for estimation and correction of non-linear mathematical model structures, within a Bayesian framework, by merging uncertain knowledge about process physics with uncertain and incomplete observations of dynamical input-state-output behavior. The resulting uncertainty in the model input-state-output mapping is expressed as a weighted combination of an uncertain conceptual model prior and a data-derived probability density function, with weights depending on the conditional data density. Our algorithm is based on the use of iterative data assimilation to update a conceptual model prior using observed system data, and thereby construct a posterior estimate of the model structure (the mathematical form of the equation itself, not just its parameters) that is consistent with both physically based prior knowledge and with the information in the data. An important aspect of the approach is that it facilitates a clear differentiation between the influences of different types of uncertainties (initial condition, input, and mapping structure) on the model prediction. Further, if some prior assumptions regarding the structural (mathematical) forms of the model equations exist, the procedure can help reveal errors in those forms and how they should be corrected. This paper examines the properties of the approach by investigating two case studies in considerable detail. The results show how, and to what degree, the structure of a dynamical hydrological model can be estimated without little or no prior knowledge (or under conditions of incorrect prior information) regarding the functional forms of the storage–streamflow and storage–evapotranspiration relationships. The importance and implications of careful specification of the model prior are illustrated and discussed.     

On characterizing the temporal dominance patterns of model parameters and processes

Diagnostic analyses of hydrological models intend to improve the understanding of how processes and their dynamics are represented in models. Temporal patterns of parameter dominance could be precisely characterized with a temporally resolved parameter sensitivity analysis. In this way, the discharge conditions are characterized, that lead to a parameter dominance in the model. To achieve this, the analysis of temporal dynamics in parameter sensitivity is enhanced by including additional information in a three‐tiered framework on different aggregation levels. Firstly, temporal dynamics of parameter sensitivity provide daily time series of their sensitivities to detect variations in the dominance of model parameters. Secondly, the daily sensitivities are related to the flow duration curve (FDC) to emphasize high sensitivities of model parameters in relation to specific discharge magnitudes. Thirdly, parameter sensitivities are monthly averaged separately for five segments of the FDC to detect typical patterns of parameter dominances for different discharge magnitudes. The three methodical steps are applied on two contrasting catchments (upland and lowland catchment) to demonstrate how the temporal patterns of parameter dynamics represent different hydrological regimes. The discharge dynamic in the lowland catchment is controlled by groundwater parameters for all discharge magnitudes. In contrast, different processes are relevant in the upland catchment, because the dominances of parameters from fast and slow runoff components in the upland catchment are changing over the year for the different discharge magnitudes. The joined interpretation of these three diagnostic steps provides deeper insights of how model parameters represent hydrological dynamics in models for different discharge magnitudes. Thus, this diagnostic framework leads to a better characterization of model parameters and their temporal dynamics and helps to understand the process behaviour in hydrological models. Copyright © 2015 John Wiley & Sons, Ltd.     

Temporal variation of b value associated with M ��4 earthquakes in the reservoir-triggered seismic environment of the Koyna�CWarna region, Western India

It is generally found that the b values associated with reservoir-triggered seismicity (RTS) are higher than the regional b values in the frequency magnitude relation of earthquakes. In the present study, temporal and spatial variation of b value is investigated using a catalog of 3,000 earthquakes from August 2005 through December 2010 for the Koyna?CWarna region in Western India, which is a classical site of RTS globally. It is an isolated (30?×?20?km²) zone of seismicity where earthquakes of up to M ??5 are found to occur during phases of loading and unloading of the Koyna and Warna reservoirs situated 25?km apart. For the Warna region, it is found that low b values of 0.6?C0.9 are associated with earthquakes of M ??4 during the loading phase. The percentage correlation of the occurrence of an M????4 earthquake with a low b value outside the 1?? or 2?? level is as high as 78?%. A drastic drop in the b value of about 50?% being reported for an RTS site may be an important precursory parameter for short-term earthquake forecast in the future.