Deccan plume,lithosphere rifting,and volcanism in Kutch,India

Kutch (northwest India) experienced lithospheric thinning due to rifting and tholeiitic and alkalic volcanism related to the Deccan Traps K/T boundary event. Alkalic lavas, containing mantle xenoliths, form plug-like bodies that are aligned along broadly east–west rift faults. The mantle xenoliths are dominantly spinel wehrlite with fewer spinel lherzolite. Wehrlites are inferred to have formed by reaction between transient carbonatite melts and lherzolite forming the lithosphere. The alkalic lavas are primitive (Mg# = 64–72) relative to the tholeiites (Mg# = 38–54), and are enriched in incompatible trace elements. Isotope and trace element compositions of the tholeiites are similar to what are believed to be the crustally contaminated Deccan tholeiites from elsewhere in India. In terms of Hf, Nd, Sr, and Pb isotope ratios, all except two alkalic basalts plot in a tight cluster that largely overlap the Indian Ridge basalts and only slightly overlap the field of Reunion lavas. This suggests that the alkalic magmas came largely from the asthenosphere mixed with Reunion-like source that welled up beneath the rifted lithosphere. The two alkalic outliers have an affinity toward Group I kimberlites and may have come from an old enriched (metasomatized) asthenosphere. We present a new model for the metasomatism and rifting of the Kutch lithosphere, and magma generation from a CO₂-rich lherzolite mantle. In this model the earliest melts are carbonatite, which locally metasomatized the lithosphere. Further partial melting of CO₂-rich lherzolite at about 2–2.5 GPa from a mixed source of asthenosphere and Reunion-like plume material produced the alkalic melts. Such melts ascended along deep lithospheric rift faults, while devolatilizing and exploding their way up through the lithosphere. Tholeiites may have been generated from the main plume head further south of Kutch.     

Frequency–Magnitude Distribution of −3.7 ≤ M W ≤ 1 Mining-Induced Earthquakes Around a Mining Front and b Value Invariance with Post-Blast Time

We investigated frequency-magnitude distribution (FMD) of acoustic emissions (AE) occurring near an active mining front in a South African gold mine, using a catalog developed from an AE network, which is capable of detecting AEs down to M _W  ?5. When records of blasts were removed, FMDs of AEs obeyed a Gutenberg?Richter law with similar b values, not depending on post-blasting time from the initial 1-min interval through more than 30 h. This result denies a suggestion in a previous study (Richardson and Jordan Bull Seismol Soc Am, 92:1766–1782, 2002) that new fractures generated by blasting disturb the size distribution of background events, which they interpreted as slip events on existing weak planes. Our AE catalog showed that the GR law with b ～ 1.2 was valid between M _W  ?3.7 and 0 for AEs around the mining front. Further, using the mine’s seismic catalog, which covers a longer time period of the same area, we could extend the validity range of the GR law with the same b value up to M _W 1.     

Recording-based identification of site liquefaction

Reconnaissance reports and pertinent research on seismic hazards show that liquefaction is one of the key sources of damage to geotechnical and structural engineering systems. Therefore, identifying site liquefaction conditions plays an important role in seismic hazard mitigation. One of the widely used approaches for detecting liquefaction is based on the time-frequency analysis of ground motion recordings, in which short-time Fourier transform is typically used. It is known that recordings at a site with liquefaction are the result of nonlinear responses of seismic waves propagating in the liquefied layers underneath the site. Moreover, Fourier transform is not effective in characterizing such dynamic features as time-dependent frequency of the recordings rooted in nonlinear responses. Therefore, the aforementioned approach may not be intrinsically effective in detecting liquefaction. An alternative to the Fourier-based approach is presented in this study, which proposes time-frequency analysis of earthquake ground motion recordings with the aid of the Hilbert-Huang transform （HHT）, and offers justification for the HHT in addressing the liquefaction features shown in the recordings. The paper then defines the predominant instantaneous frequency （PIF） and introduces the PiF-related motion features to identify liquefaction conditions at a given site. Analysis of 29 recorded data sets at different site conditions shows that the proposed approach is effective in detecting site liquefaction in comparison with other methods.     

Effect of water on the fluorine and chlorine partitioning behavior between olivine and silicate melt

Halogens show a range from moderate (F) to highly (Cl, Br, I) volatile and incompatible behavior, which makes them excellent tracers for volatile transport processes in the Earth’s mantle. Experimentally determined fluorine and chlorine partitioning data between mantle minerals and silicate melt enable us to estimate Mid Ocean Ridge Basalt (MORB) and Ocean Island Basalt (OIB) source region concentrations for these elements. This study investigates the effect of varying small amounts of water on the fluorine and chlorine partitioning behavior at 1280?°C and 0.3 GPa between olivine and silicate melt in the Fe-free CMAS+F–Cl–Br–I–H₂O model system. Results show that, within the uncertainty of the analyses, water has no effect on the chlorine partitioning behavior for bulk water contents ranging from 0.03 (2) wt% H₂O (D_Cl ^ol/melt = 1.6?±?0.9 × 10^?4) to 0.33 (6) wt% H₂O (D_Cl ^ol/melt = 2.2?±?1.1 × 10^?4). Consequently, with the effect of pressure being negligible in the uppermost mantle (Joachim et al. Chem Geol 416:65–78, 2015), temperature is the only parameter that needs to be considered for the determination of chlorine partition coefficients between olivine and melt at least in the simplified iron-free CMAS+F–Cl–Br–I–H₂O system. In contrast, the fluorine partition coefficient increases linearly in this range and may be described at 1280?°C and 0.3 GPa with (R ²?=?0.99): \(D_{F}^{\text{ol/melt}}\ =\ 3.6\pm 0.4\ \times \ {{10}^{-3}}\ \times \ {{X}_{{{\text{H}}_{\text{2}}}\text{O}}}\left( \text{wt }\!\!\%\!\!\text{ } \right)\ +\ 6\ \pm \ 0.4\times \,{{10}^{-4}}\). The observed fluorine partitioning behavior supports the theory suggested by Crépisson et al. (Earth Planet Sci Lett 390:287–295, 2014) that fluorine and water are incorporated as clumped OH/F defects in the olivine structure. Results of this study further suggest that fluorine concentration estimates in OIB source regions are at least 10% lower than previously expected (Joachim et al. Chem Geol 416:65–78, 2015), implying that consideration of the effect of water on the fluorine partitioning behavior between Earth’s mantle minerals and silicate melt is vital for a correct estimation of fluorine abundances in OIB source regions. Estimates for MORB source fluorine concentrations as well as chlorine abundances in both mantle source regions are within uncertainty not affected by the presence of water.     

Differing responses of zircon,chevkinite-(Ce), monazite-(Ce) and fergusonite-(Y) to hydrothermal alteration: Evidence from the Keivy alkaline province,Kola Peninsula,Russia

Mineralogy and Petrology - A quartzolite from the Rova occurrence, Keivy alkali granite province, Kola Peninsula, Russia, is used to examine the differing responses of certain rare-metal minerals...     

Recharge mechanism and processes controlling groundwater chemistry in a Precambrian sedimentary terrain: a case study from Central India

Unplanned and unsustainable extraction has created stress on groundwater resources in many parts of India. The stress symptoms are more pronounced in hard rock areas, where the aquifer potentials are comparatively low. Present research targeted an area of 3000 km² in the interstream region between the Kharun and Seonath rivers, which is one such region in Central India. In spite of being a water-stressed area, so far, little is understood about processes of recharge, amount of recharge and processes controlling chemical quality, which are key inputs for groundwater management in the area. This study presents an appraisal of recharge mechanism, recharge rate and prevailing water–rock interactions in the study area. Stable isotope composition of groundwater when compared to that of rainfall indicates monsoon rainfall as the primary source of groundwater recharge. Winter rains, which are characteristically enriched in heavier isotopes, do not contribute notably to groundwater recharge. Recharge is rapid with minor or no evaporative enrichment before recharge. Further, analysis of stable isotopes show that ‘macropore recharge’ is dominant in limestone or calcareous shale, covering more than 70% of the study area. Also apparent is the vertical connectivity amongst the aquifers. However, active intermixing of surface water and groundwater is not a predominant process. Annual groundwater recharge from rainfall, as derived from chloride mass balance, is 105.26 million cubic metre. Groundwater is predominantly of bicarbonate type, irrespective of its hydrostratigraphic (lithology) setting. Dissolution of carbonates and gypsum (occurring as veins), weathering of feldspar and ion exchange of clay minerals are amongst the most likely processes controlling the regional groundwater chemistry.     

Simulations of energy balance components at snow-dominated montane watershed by land surface models

The quantification of energy interactions among land surface, atmosphere, and surface vegetation is significant to comprehend the hydrological cycle in montane watersheds. Moreover, elevation change is an essential in causing variations in energy fluxes. Thus, estimating the major components of energy interactions is essential for better understanding of the hydrological process. The advanced land surface models (LSMs); the common land model (CLM) and variables infiltration capacity (VIC) are used to estimate accurate hydrometeorological variables. These hydrometeorological variables such as net radiation and sensible, latent, and ground heat fluxes were estimated using CLM and VIC at upper and lower meteorological stations in Sierra Nevada Mountain, California, USA. The estimated fluxes were compared with observations at each site. The estimated daily and monthly net radiation and sensible heat flux from both models showed good agreement with the observations (R ≥ 0.84). The CLM-modeled estimates showed lower trends during the rainfall periods, which occurred mainly during winter at both sites. In comparison, the estimated daily and monthly latent heat flux from CLM at both sites showed better results with lower RMSE and bias than that from VIC, which underestimated latent heat flux. Both models overestimated ground heat flux, and the variation trend was similar to observation. For sensitivity analysis, according to elevation change, all the estimated energy fluxes had slightly different values at the upper and lower met stations. In future studies, parameterization for the LSMs will be conducted for more robust estimations of hydrometeorological variables in montane watersheds.     

Assimilation of altimetry data for nonlinear shallow-water tides: Quarter-diurnal tides of the Northwest European Shelf

Non-linear tidal constituents, such as the overtide M₄ or the compound tide MS₄, are generated by interaction in shallow seas of the much larger astronomically forced “primary” tidal constituents (e.g., M₂, S₂). As such, errors in modeling these “secondary” shallow-water tides might be expected to be caused first of all by errors in modeling the primary constituents. Thus, in the context of data assimilation, observations of primary-constituent harmonic constants can indirectly constrain shallow-water constituents. Here we consider variational data assimilation for primary and secondary tidal constituents as a coupled problem, using a simple linearized perturbation theory for weak interactions of the dominant primary constituents. Variation of the resulting penalty functional leads to weakly non-linear Euler–Lagrange equations, which we show can be solved approximately with a simple two-stage scheme. In the first stage, data for the primary constituents are assimilated into the linear shallow water equations (SWE), and the resulting inverse solutions are used to compute the quadratic interactions in the non-linear SWE that constitute the forcing for the secondary constituents. In the second stage, data for the compound or overtide constituent are assimilated into the linear SWE, using a prior forced by the results of the first stage. We apply this scheme to assimilation of TOPEX/Poseidon and Jason altimetry data on the Northwest European Shelf, comparing results to a large set of shelf and coastal tide gauges. Prior solutions for M₄, MS₄ and MN₄ computed using inverse solutions for M₂, S₂, and N₂ dramatically improve fits to validation tide gauges relative to unconstrained forward solutions. Further assimilation of along-track harmonic constants for these shallow-water constituents reduces RMS differences to below 1 cm on the shelf, approaching the accuracy of the validation tide gauge harmonic constants.     

Rhyolites at Kerlingarfjöll, Iceland: the evolution and lifespan of silicic central volcanoes

Kerlingarfjöll central volcano is Iceland’s second largest outcrop of Quaternary rhyolite and is part of the Icelandic Western Rift Zone. Geochemical and Ar/Ar age data show that at least 21 different rhyolite eruptions have taken place at Kerlingarfjöll over the last 350 ka. Ar/Ar dating was carried out on samples of obsidian which showed variable reproducibility, illustrating the difficulty in dating young Icelandic volcanics. Nevertheless, reasonable estimates of eruption age have been derived for a number of eruptive units that are consistent with observed stratigraphy, enabling an understanding of the temporal evolution of Kerlingarfjöll. Two rhyolite magma types are present. The first is an older, low-Nb rhyolite that was erupted episodically along a cryptic curved fracture system, to form a discontinuous ring of rhyolite mountains, between 350 and 250 ka. This discontinuous ring is similar to structures observed at other volcanoes in Iceland, suggesting that the development of a curved fracture that acts as a pathway for episodic silicic eruptions is a feature of central volcano development. The second magma is a younger, high-Nb rhyolite that was erupted episodically between 250 and 68 ka in the northern part of Kerlingarfjöll, forming two clusters, both of which have areas of intense hydrothermal activity. Repose periods for rhyolite volcanism are thought to be on the order of tens of thousands of years, and it is possible that Kerlingarfjöll will erupt rhyolite again in the future.     

Predicting phosphorus dynamics in complex terrains using a variable source area hydrology model

Phosphorus (P) loss from agricultural watersheds has long been a critical water quality problem, the control of which has been the focus of considerable research and investment. Preventing P loss depends on accurately representing the hydrological and chemical processes governing P mobilization and transport. The Soil and Water Assessment Tool (SWAT) is a watershed model commonly used to predict run‐off and non‐point source pollution transport. SWAT simulates run‐off employing either the curve number (CN) or the Green and Ampt methods, both assume infiltration‐excess run‐off, although shallow soils underlain by a restricting layer commonly generate saturation‐excess run‐off from variable source areas (VSA). In this study, we compared traditional SWAT with a re‐conceptualized version, SWAT‐VSA, that represents VSA hydrology, in a complex agricultural watershed in east central Pennsylvania. The objectives of this research were to provide further evidence of SWAT‐VSA's integrated and distributed predictive capabilities against measured surface run‐off and stream P loads and to highlight the model's ability to drive sub‐field management of P. Thus, we relied on a detailed field management database to parameterize the models. SWAT and SWAT‐VSA predicted discharge similarly well (daily Nash–Sutcliffe efficiencies of 0.61 and 0.66, respectively), but SWAT‐VSA outperformed SWAT in predicting P export from the watershed. SWAT estimated lower P loss (0.0–0.25 kg ha^?1) from agricultural fields than SWAT‐VSA (0.0–1.0+ kg ha^?1), which also identified critical source areas – those areas generating large run‐off and P losses at the sub‐field level. These results support the use of SWAT‐VSA in predicting watershed‐scale P losses and identifying critical source areas of P loss in landscapes with VSA hydrology. Copyright © 2014 John Wiley & Sons, Ltd.