Semi-empirical model for retrieval of soil moisture using RISAT-1 C-Band SAR data over a sub-tropical semi-arid area of Rewari district,Haryana (India)

We have estimated soil moisture (SM) by using circular horizontal polarization backscattering coefficient (\(\sigma ^{\mathrm{o}}_{\mathrm{RH}}\)), differences of circular vertical and horizontal \(\sigma ^{\mathrm{o}} \, (\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}})\) from FRS-1 data of Radar Imaging Satellite (RISAT-1) and surface roughness in terms of RMS height (\({\hbox {RMS}}_{\mathrm{height}}\)). We examined the performance of FRS-1 in retrieving SM under wheat crop at tillering stage. Results revealed that it is possible to develop a good semi-empirical model (SEM) to estimate SM of the upper soil layer using RISAT-1 SAR data rather than using existing empirical model based on only single parameter, i.e., \(\sigma ^{\mathrm{o}}\). Near surface SM measurements were related to \(\sigma ^{\mathrm{o}}_{\mathrm{RH}}\), \(\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}}\) derived using 5.35 GHz (C-band) image of RISAT-1 and \({\hbox {RMS}}_{\mathrm{height}}\). The roughness component derived in terms of \({\hbox {RMS}}_{\mathrm{height}}\) showed a good positive correlation with \(\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}} \, (R^{2} = 0.65)\). By considering all the major influencing factors (\(\sigma ^{\mathrm{o}}_{\mathrm{RH}}\), \(\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}}\), and \({\hbox {RMS}}_{\mathrm{height}}\)), an SEM was developed where SM (volumetric) predicted values depend on \(\sigma ^{\mathrm{o}}_{\mathrm{RH}}\), \(\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}}\), and \({\hbox {RMS}}_{\mathrm{height}}\). This SEM showed \(R^{2}\) of 0.87 and adjusted \(R^{2}\) of 0.85, multiple R=0.94 and with standard error of 0.05 at 95% confidence level. Validation of the SM derived from semi-empirical model with observed measurement (\({\hbox {SM}}_{\mathrm{Observed}}\)) showed root mean square error (RMSE) = 0.06, relative-RMSE (R-RMSE) = 0.18, mean absolute error (MAE) = 0.04, normalized RMSE (NRMSE) = 0.17, Nash–Sutcliffe efficiency (NSE) = 0.91 (\({\approx } 1\)), index of agreement (d) = 1, coefficient of determination \((R^{2}) = 0.87\), mean bias error (MBE) = 0.04, standard error of estimate (SEE) = 0.10, volume error (VE) = 0.15, variance of the distribution of differences \(({\hbox {S}}_{\mathrm{d}}^{2}) = 0.004\). The developed SEM showed better performance in estimating SM than Topp empirical model which is based only on \(\sigma ^{\mathrm{o}}\). By using the developed SEM, top soil SM can be estimated with low mean absolute percent error (MAPE) = 1.39 and can be used for operational applications.     

Effect of uncertainty in $$V_{\mathrm{S}}{-}N$$ correlations on seismic site response analysis

This paper studies the reliability of the calculated shear-wave velocity (\(V_{\mathrm{S}}\)) from different available \(V_{\mathrm{S}}\) and SPT-N correlations in terms of seismic site response analysis. In the present study, various \(V_{\mathrm{S}}\)–N correlations developed for different regions around the globe have been used to calculate the bound of \(V_{\mathrm{S}}\) variations with depth at three different sites in Kolkata city. This bound has later been used to generate the random \(V_{\mathrm{S}}\) profiles using the Monte Carlo simulation. Equivalent linear site response analysis has been performed to study the response of those generated profiles under different input motion excitations. Strong-to-weak ground motion records have been used for this purpose. The amplification spectra of the generated \(V_{\mathrm{S}}\) profiles using all soil types and specific soil-type \(V_{\mathrm{S}}\)–N correlations show significant variations. The study also shows that the \(V_{\mathrm{S}}\)–N correlation may result in quite different \(V_{\mathrm{S},30}\) values and subsequently it may lead to the different site classes according to the NEHRP 2003 classification. So, the random choice of the \(V_{\mathrm{S}}\)–N correlation, where the direct measurement of \(V_{\mathrm{S}}\) is not available, may affect the outcome of seismic hazard analysis significantly. The study points out the need for accurate estimation of the \(V_{\mathrm{S}}\) profile either from in-situ determination or using site-specific correlation.     

Characteristics of surface ozone in Agra,a sub-urban site in Indo-Gangetic Plain

In the present study, measurements of surface ozone (\(\hbox {O}_{3}\)) and its precursors (NO and \(\hbox {NO}_{2}\)) were carried out at a sub-urban site of Agra (\(27{^{\circ }}10'\hbox {N}\), \(78{^{\circ }}05'\hbox {E}\)), India during May 2012–May 2013. During the study period, average concentrations of \(\hbox {O}_{3}\), NO, and \(\hbox {NO}_{2}\) were \(39.6 \pm 25.3\), \(0.8 \pm 0.8\) and \(9.1 \pm 6.6 \, \hbox {ppb}\), respectively. \(\hbox {O}_{3}\) showed distinct seasonal variation in peak value of diurnal variation: summer \({>}\) post-monsoon \({>}\) winter \({>}\) monsoon. However, \(\hbox {NO}_{2}\) showed highest levels in winter and lowest in monsoon. The average positive rate of change of \(\hbox {O}_{3}\) (08:00–11:00 hr) was highest in April (16.3 ppb/hr) and lowest in August (1.1 ppb/hr), while average negative rate of change of \(\hbox {O}_{3}\) (17:00–19:00 hr) was highest in December (–13.2 ppb/hr) and lowest in July (–1.1 ppb/hr). An attempt was made to identify the \(\hbox {VOC--NO}_{\mathrm{x}}\) sensitivity of the site using \(\hbox {O}_{3}/\hbox {HNO}_{3}\) ratio as photochemical indicator. Most of the days this ratio was above the threshold value (12–16), which suggests \(\hbox {NO}_{\mathrm{x}}\) sensitivity of the site. The episodic event of ozone was characterized through meteorological parameters and precursors concentration. Fine particles (\(\hbox {PM}_{2.5}\)) cause loss of ozone through heterogeneous reactions on their surface and reduction in solar radiation. In the study, statistical analyses were used to estimate the amount of ozone loss.     

New paleomagnetic results on $$\sim $$2367 Ma Dharwar giant dyke swarm,Dharwar craton,southern India: implications for Paleoproterozoic continental reconstruction

Here we report new paleomagnetic results and precise paleopole position of the extensional study on \(\sim \)2367 Ma mafic giant radiating dyke swarm in the Dharwar craton, southern India. We have sampled 29 sites on 12 dykes from NE–SW Karimnagar–Hyderabad dykes and Dhone–Gooty sector dykes, eastern Dharwar craton to provide unambiguous paleomagnetism evidence on the spectacular radiating dyke swarm and thereby strengthening the presence of single magmatic event at \(\sim \)2367 Ma. A total of 158 samples were subjected to detailed alternating field and thermal demagnetization techniques and the results are presented here along with previously reported data on the same dyke swarm. The remanent magnetic directions are showing two components, viz., seven sites representing four dykes show component (A) with mean declination of \(94{{}^{\circ }}\) and mean inclination of \(-\,70{{}^{\circ }}\) (\(\hbox {k}=87\), \(\upalpha _{95}=10{{}^{\circ }}\)) and corresponding paleopole at \(16{{}^{\circ }}\hbox {N}\), \(41{{}^{\circ }}\hbox {E}\) (\(\hbox {dp}=15{{}^{\circ }}\) and \(\hbox {dm}=17{{}^{\circ }}\)) and 22 sites representing 8 dykes yielded a component (B) with mean declination of \(41{{}^{\circ }}\) and mean inclination of \(-\,21{{}^{\circ }}\) (\(\hbox {k}=41\), \(\upalpha _{95}=9{{}^{\circ }}\)) with a paleopole at \(41{{}^{\circ }}\hbox {N}\), \(200{{}^{\circ }}\hbox {E}\) (\(\hbox {dp}=5{{}^{\circ }}\) and \(\hbox {dm}=10{{}^{\circ }}\)). Component (A) results are similar to the previously reported directions from the \(\sim \)2367 Ma dyke swarm, which have been confirmed fairly reliably to be of primary origin. The component (B) directions appear to be strongly overprinted by the 2080 Ma event. The grand mean for the primary component (A) combined with earlier reported studies gives mean declination of \(97{{}^{\circ }}\) and mean inclination of \(-\,79{{}^{\circ }}\) (\(\hbox {k}=55\), \(\upalpha _{95}=3{{}^{\circ }}\)) with a paleopole at \(15{{}^{\circ }}\hbox {N}\), \(57{{}^{\circ }}\hbox {E}\) (\(\hbox {dp}=5{{}^{\circ }}\), \(\hbox {dm}=6{{}^{\circ }}\)). Paleogeographical position for the Dharwar craton at \(\sim \)2367 Ma suggests that there may be a chance to possible spatial link between Dharwar dykes of Dharwar craton (India), Widgemooltha and Erayinia dykes of Yilgarn craton (Australia), Sebanga Poort Dykes of Zimbabwe craton (Africa) and Karelian dykes of Kola-Karelia craton (Baltica Shield).     

Episodic crustal growth in the Bundelkhand craton of central India shield: Constraints from petrogenesis of the tonalite–trondhjemite–granodiorite gneisses and K-rich granites of Bundelkhand tectonic zone

Tonalite–trondhjemite–granodiorite gneisses (TTG) and K-rich granites are extensively exposed in the Mesoarchean to Paleoproterozoic Bundelkhand craton of central India. The TTGs rocks are coarse- grained with biotite, plagioclase feldspar, K-feldspar and amphibole as major constituent phases. The major minerals constituting the K-rich granites are K-feldspar, plagioclase feldspar and biotite. They are also medium to coarse grained. Mineral chemical studies show that the amphiboles of TTG are calcic amphibole hastingsite, plagioclase feldspars are mostly of oligoclase composition, K-feldspars are near pure end members and biotites are solid solutions between annite and siderophyllite components. The K-rich granites have biotites of siderophyllite–annite composition similar to those of TTGs, plagioclase feldspars are oligoclase in composition, potassic feldspars have \(\hbox {X}_{\mathrm{K}}\) ranging from 0.97 to 0.99 and are devoid of any amphibole. The tonalite–trondhjemite–granodiorite gneiss samples have high \(\hbox {SiO}_{2}\) (64.17–74.52 wt%), \(\hbox {Na}_{2}\hbox {O}\) (3.11–5.90 wt%), low Mg# (30–47) and HREE contents, with moderate \((\hbox {La/Yb})_{\mathrm{CN}}\) values (14.7–33.50) and Sr/Y ratios (4.85–98.7). These geochemical characteristics suggest formation of the TTG by partial melting of the hydrous basaltic crust at pressures and depths where garnet and amphibole were stable phases in the Paleo-Mesoarchean. The K-rich granite samples show high \(\hbox {SiO}_{2}\) (64.72–76.73 wt%), \(\hbox {K}_{2}\hbox {O}\) (4.31–5.42), low \(\hbox {Na}_{2}\hbox {O}\) (2.75–3.31 wt%), Mg# (24–40) and HREE contents, with moderate to high \((\hbox {La/Yb})_{\mathrm{CN}}\) values (9.26–29.75) and Sr/Y ratios (1.52–24). They differ from their TTG in having elevated concentrations of incompatible elements like K, Zr, Th, and REE. These geochemical features indicate formation of the K-granites by anhydrous partial melting of the Paleo-Mesoarchean TTG or mafic crustal materials in an extensional regime. Combined with previous studies it is interpreted that two stages of continental accretion (at 3.59–3.33 and 3.2–3.0 Ga) and reworking (at 2.5–1.9 Ga) occurred in the Bundelkhand craton from Archaean to Paleoproterozoic.     

Competing effects of crystal chemistry and silicate melt composition on trace element behavior in magmatic systems: insights from crystal/silicate melt partitioning of the REE,HFSE, Sn,In, Ga,Ba, Pt and Rh

We present new partition coefficients for the REE, HFSE, Sn, In, Ga, Ba, Pt and Rh between clinopyroxene, olivine and basaltic melt as a function of crystal chemistry and melt composition at temperatures of 1190–1300 °C and 1-bar pressure. Two components, namely \(\mathrm {Al_2O_3}\) and \(\mathrm {Na_2O}\), were chosen to be investigated since they are known to affect the structure of silicate melts and especially clinopyroxene crystal chemistry. The amount of \(^{[4]}\mathrm{Al}\) in clinopyroxene will result in an increase of \(D_i^\mathrm{{cpx/melt}}\) even after applying a correction factor to account for the effect of melt polymerization. Moreover, the positive correlation between \(^{[4]}\mathrm{Al}\) and \(D_i^\mathrm{{cpx/melt}}\) is not restricted to the REE, but also applies for Sn, Ga, In, and Ba. The addition of up to 2.6 wt% \(\mathrm {Na_2O}\) to the silicate melt universally increases the \(D_i^\mathrm{{cpx/melt}}\) without any concomitant change in crystal chemistry or a significant effect in melt polymerization. This compositional effect is likely due to the ability of Na to break REE–Al complexes in the melt. Our results emphasize the importance of considering all variables that affect the behavior of trace elements in magmatic systems before applying the lattice strain model and derive meaningful results for the changes in the parameters of the crystallographic sites.     

Mineral shock signatures in rocks from Dhala (Mohar) impact structure,Shivpuri district,Madhya Pradesh,India

A concrete study combining optical microscopy, Raman spectroscopy and X-ray diffractometry, was carried out on subsurface samples of basement granite and melt breccia from Mohar (Dhala) impact structure, Shivpuri district, Madhya Pradesh, India. Optical microscopy reveals aberrations in the optical properties of quartz and feldspar in the form of planar deformation feature-like structures, lowered birefringence and mosaics in quartz, toasting, planar fractures and ladder texture in alkali feldspar and near-isotropism in bytownite. It also brings to light incidence of parisite, a radioactive rare mineral in shocked granite. Raman spectral pattern, peak positions, peak widths and multiplicity of peak groups of all minerals, suggest subtle structural/crystallographic deviations. XRD data further reveals minute deviations of unit cell parameters of quartz, alkali feldspar and plagioclase, with respect to standard \({\upalpha }\)-quartz, high- and low albite and microcline. Reduced cell volumes in these minerals indicate compression due to pressure. The \(\hbox {c}_{0}/\hbox {a}_{0}\) values indicate an inter-tetrahedral angle roughly between \(120^{\mathrm{o}}\) and \(144^{\mathrm{o}}\), further pointing to a possible pressure maxima of around 12 GPa. The observed unit cell aberration of minerals may indicate an intermediate stage between crystalline and amorphous stages, thereby, signifying possible overprinting of decompression signatures over shock compression effects, from a shock recovery process.     

Appraisal of groundwater quality in a crystalline aquifer: a chemometric approach

The purpose of this study is to assess the groundwater quality and identify the processes that control the groundwater chemistry in a crystalline aquifer. A total of 72 groundwater samples were collected during pre- and post-monsoon seasons in the year 2014 in a semi-arid region of Gooty Mandal, Anantapur district, Andhra Pradesh, India. The study utilized chemometric analysis like basic statistics, Pearson’s correlation coefficient (r), principal component analysis (PCA), Gibbs ratio, and index of base exchange to understand the mechanism of controlling the groundwater chemistry in the study area. The results reveal that groundwater in the study area is neutral to slightly alkaline in nature. The order of dominance of cations is Na⁺ > Ca²⁺ > Mg²⁺ > K⁺ while for anions, it is \( {\mathrm{HCO}}_3^{-}>{\mathrm{Cl}}^{-} \)>\( {\mathrm{NO}}_3^{-} \)>\( {\mathrm{SO}}_4^{2-} \)>\( {\mathrm{CO}}_3^{2-}>{\mathrm{F}}^{-} \) in both seasons. Based on the Piper classification, most of the groundwater samples are identified as of sodium bicarbonate (\( {\mathrm{Na}}^{+}-{\mathrm{HCO}}_3^{-}\Big) \) type. According to the results of the principal component analysis (PCA), three factors and two factors were identified pre and post monsoon, respectively. The present study indicates that the groundwater chemistry is mostly controlled by geogenic processes (weathering, dissolution, and ion exchange) and some extent of anthropogenic activities.     

The effect of liquid composition on the partitioning of Ni between olivine and silicate melt

We report the results of experiments designed to separate the effects of temperature and pressure from liquid composition on the partitioning of Ni between olivine and liquid, \(D_{\text{Ni}}^{\text{ol/liq}}\). Experiments were performed from 1300 to 1600 °C and 1 atm to 3.0 GPa, using mid-ocean ridge basalt (MORB) glass surrounded by powdered olivine in graphite–Pt double capsules at high pressure and powdered MORB in crucibles fabricated from single crystals of San Carlos olivine at one atmosphere. In these experiments, pressure and temperature were varied in such a way that we produced a series of liquids, each with an approximately constant composition (~12, ~15, and ~21 wt% MgO). Previously, we used a similar approach to show that \(D_{\text{Ni}}^{\text{ol/liq}}\) for a liquid with ~18 wt% MgO is a strong function of temperature. Combining the new data presented here with our previous results allows us to separate the effects of temperature from composition. We fit our data based on a Ni–Mg exchange reaction, which yields \(\ln \left( {D_{\text{Ni}}^{\text{molar}} } \right) = \frac{{ -\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } }}{RT} + \frac{{\Delta _{r(1)} S_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } }}{R} - \ln \left( {\frac{{X_{\text{MgO}}^{\text{liq}} }}{{X_{{{\text{MgSi}}_{ 0. 5} {\text{O}}_{ 2} }}^{\text{ol}} }}} \right).\) Each subset of constant composition experiments displays roughly the same temperature dependence of \(D_{\text{Ni}}^{\text{ol/liq}}\) (i.e.,\(-\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\)) as previously reported for liquids with ~18 wt% MgO. Fitting new data presented here (15 experiments) in conjunction with our 13 previously published experiments (those with ~18 wt% MgO in the silicate liquid) to the above expression gives \(-\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = 3641 ± 396 (K) and \(\Delta _{r(1)} S_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = ? 1.597 ± 0.229. Adding data from the literature yields \(-\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = 4505 ± 196 (K) and \(\Delta _{r(1)} S_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = ? 2.075 ± 0.120, a set of coefficients that leads to a predictive equation for \(D_{\text{Ni}}^{\text{ol/liq}}\) applicable to a wide range of melt compositions. We use the results of our work to model the melting of peridotite beneath lithosphere of varying thickness and show that: (1) a positive correlation between NiO in magnesian olivine phenocrysts and lithospheric thickness is expected given a temperature-dependent \(D_{\text{Ni}}^{\text{ol/liq}} ,\) and (2) the magnitude of the slope for natural samples is consistent with our experimentally determined temperature dependence. Alternative processes to generate the positive correlation between NiO in magnesian olivines and lithospheric thickness, such as the melting of olivine-free pyroxenite, are possible, but they are not required to explain the observed correlation of NiO concentration in initially crystallizing olivine with lithospheric thickness.     

Thermodynamics and phonon dispersion of pyrope and grossular silicate garnets from ab initio simulations

The phonon dispersion and thermodynamic properties of pyrope (\(\hbox {Mg}_3\hbox {Al}_2\hbox {Si}_3\hbox {O}_{12}\)) and grossular (\(\hbox {Ca}_3\hbox {Al}_2\hbox {Si}_3\hbox {O}_{12}\) ) have been computed by using an ab initio quantum mechanical approach, an all-electron variational Gaussian-type basis set and the B3LYP hybrid functional, as implemented in the Crystal program. Dispersion effects in the phonon bands have been simulated by using supercells of increasing size, containing 80, 160, 320, 640, 1280 and 2160 atoms, corresponding to 1, 2, 4, 8, 16 and 27 \(\mathbf {k}\) points in the first Brillouin zone. Phonon band structures, density of states and corresponding inelastic neutron scattering spectra are reported. Full convergence of the various thermodynamic properties, in particular entropy (S) and specific heat at constant volume (\(C_\mathrm{{V}}\)), with the number of \(\mathbf {k}\) points is achieved with 27 \(\mathbf {k}\) points. The very regular behavior of the S(T) and \(C_\mathrm{{V}}(T)\) curves as a function of the number of \(\mathbf {k}\) points, determined by high numerical stability of the code, permits extrapolation to an infinite number of \(\mathbf {k}\) points. The limiting value differs from the 27-\(\mathbf {k}\) case by only 0.40 % at 100 K for S (the difference decreasing to 0.11 % at 1000 K) and by 0.29 % (0.05 % at 1000 K) for \(C_\mathrm{{V}}\). The agreement with the experimental data is rather satisfactory. We also address the problem of the relative entropy of pyrope and grossular, a still debated question. Our lattice dynamical calculations correctly describe the larger entropy of pyrope than grossular by taking into account merely vibrational contributions and without invoking “static disorder” of the Mg ions in dodecahedral sites. However, as the computed entropy difference is found to be larger than the experimental one by a factor of 2–3, present calculations cannot exclude possible thermally induced structural changes, which could lead to further conformational contributions to the entropy.     

Strain responses of frozen clay with thermal gradient under triaxial creep

Thermal gradient is one of the main features for the temperature distribution in artificial frozen shaft lining (FSL). The time-dependent strain responses and the corresponding heterogeneity characteristics of frozen soils with thermal gradient are of potential significance for stability assessment and prediction of FSL, especially of the FSL embedded in thick alluvium. A series of triaxial creep tests were carried out on frozen saturated clay under various thermal gradients and creep stresses. The experimental results indicated that the triaxial creep curves for frozen clay with thermal gradient exhibit viscous characteristics, and the creep rate \(\Delta \varepsilon_{\text{a}} /\Delta t\) decreases with the increase in creep time \(t\) and decrease in thermal gradient. The stress–strain curve under different \(t\) showed that the creep stress has a marked growth when axial strain \(\varepsilon_{\text{a}} \le 1\,\%\). However, when \(\varepsilon_{\text{a}} \ge 1\,\%\), the growth rate decreases gradually. The deviation between measured radial strain \(\varepsilon_{\text{r}}^{\text{m}}\) under the middle specimen section height SSH and the calculated radial strain \(\varepsilon_{\text{r}}^{\text{c}}\) from the volumetric strain increases following a unified equation with the increase in axial strain. The radial strain \(\varepsilon_{\text{r}}^{\text{f}}\) for frozen clay with thermal gradient after experiment increases with the increase in SSH, and the slope of \( \varepsilon_{\text{r}}^{\text{f}} - {\text{SSH}} \) curve is significantly dependent on the thermal gradient and creep stress. The variation of \(\varepsilon_{\text{r}}^{\text{m}} - \varepsilon_{\text{r}}^{\text{c }}\) during experiment and \(\varepsilon_{\text{r}}^{\text{f}}\) distribution after experiment are the macro-responses of internal micro-heterogeneities in frozen soils induced from thermal gradient, and are closely related to strain rate and its variation. These observations and findings provide an insight into the creep mechanism and the estimation method of creep deformation for frozen soils with thermal gradient.     

Physico-chemical conditions of four calc-alkaline granitoid plutons of Chhotanagpur Gneissic Complex,eastern India: Tectonic implications

Petrography and mineralogy of four calc-alkaline granitoid plutons Agarpur, Sindurpur, Raghunathpur and Sarpahari located from west to east of northern Purulia of Chhotanagpur Gneissic Complex, eastern India, are investigated. The plutons, as a whole, are composed of varying proportions of Qtz–Pl–Kfs–Bt–Hbl±Px–Ttn–Mag–Ap–Zrn±Ep. The composition of biotite is consistent with those of calc-alkaline granitoids. Hornblende–plagioclase thermometry, aluminium-in-hornblende barometry and the assemblage sphene–magnetite–quartz were used to determine the P, T and \(f_{\mathrm{O}_2}\) during the crystallisation of the parent magmas in different plutons. The plutons are crystallised under varying pressures (6.2–2.4 kbar) and a wide range of temperatures (896–\(718{^{\circ }}\hbox {C}\)) from highly oxidised magmas (log \(f_{\mathrm{O}_2}\) \(-11.2\) to \(-15.4\) bar). The water content of the magma of different plutons varied from 5.0 to 6.5 wt%, consistent with the calc-alkaline nature of the magma. Calc-alkaline nature, high oxygen fugacity and high \(\hbox {H}_{2}\hbox {O}_{{\mathrm{melt}}}\) suggest that these plutons were emplaced in subduction zone environment. The depths of emplacement of these plutons seem to increase from west to east. Petrologic compositions of these granitoids continuously change from enderbite (opx-tonalite: Sarpahari) in the east to monzogranite (Raghunathpur) to syenogranite (Sindurpur) to alkali feldspar granite (Agarpur) in the west. The water contents of the parental magmas of different plutons also increase systematically from east to west. No substantial increase in the depth of emplacement is found in these plutons lying south and north of the major shear zone passing through the study area suggesting the strike-slip nature of the east–west shear zone.     

Ionic conductivity in gem-quality single-crystal alkali feldspar from the Eifel: temperature,orientation and composition dependence

We measured the ion conductivity of single-crystal alkali feldspar originating from two different locations in the Eifel/Germany, named Volkesfeld and Rockeskyller sanidine and having potassium site fractions \(C_\mathrm{K}\) of 0.83 and 0.71, respectively. The dc conductivities resulting from electrochemical impedance spectroscopy over the temperature range of 300–900 \(^{\circ }\hbox {C}\) show a weak composition dependence but pronounced differences between the b-direction [\(\perp (010)\)] and \(c^{*}\)-direction [\(\perp (001)\)] of the monoclinic feldspar structure. Conductivity activation energies obtained from the observed linear Arrhenius plots are close to 1.2 eV in all cases, which is closely similar to the activation energies of the \(^{22}\mathrm{Na}\) tracer diffusivity in the same crystals. Taking into account literature data on K tracer diffusion and diffusion correlation effects, the present results point to a predominance of the interstitialcy mechanism over the vacancy mechanism in mass and charge transport on the alkali sublattice in potassium-rich alkali feldspar.     

Data-driven nonlinear modeling studies on removal of Acid Yellow 59 using Si-doped multi-walled carbon nanotubes

Data-driven modeling of removal of color index name of Acid Yellow 59 from aqueous solutions using multi-walled carbon nanotubes by multiple (non)linear regression and artificial neural networks (ANN) models based on leave-one-out cross-validation to predict the adsorbed dye amount per unit mass of adsorbent (mg g^?1) and performance evaluation of the proposed multiple (non)linear regression and ANN models is the main novel contributor of the present study. Initial dye concentration, adsorbent concentration, reaction time, and temperature were determined as explanatory variables and input neurons for multiple (non)linear regression and ANN models, respectively. The total number of experiments was determined as 1280 statistically. The results showed that multilayer perception ANN model (\(R^{2}_{\text{training}}\) = 0.9997, \(R^{2}_{\text{testing}}\) = 0.9993, RMSE = 0.7678, MAE of 0.0007) predicted q _t better than multiple (non)linear regression model (\(R^{2}_{\text{adj}}\) = 0.9645, \(R^{2}_{\text{pred}}\) = 0.9633, SE = 9.55) and MLR (R ² = 0.9543, SE = 10.87) models. The results justified the accuracy of ANN in prediction of q _t , significantly.     

The Numba ductile deformation zone (northwest Cameroon): A geometric analysis of folds based on the Fold Profiler method

The Numba ductile deformation zone (NDDZ) is characterised by folds recorded during the three deformation phases that affected the banded amphibole gneiss. Fold-shape analyses using the program Fold Profiler with the aim to show the importance of folding events in the structural analysis of the NDDZ and its contribution to the Pan-African orogeny in central Africa have been made. Classical field method, conic sections method and Ramsay’s fold classification method were applied to (i) have the general orientation of folds, (ii) analyze the fold shapes and (iii) classify the geometry of the folded bands. Fold axes in banded amphibole gneiss plunge moderately (\(<\)15\(^{\circ }\)) towards the NNE or SSW. The morphology of F\(_{1}\), F\(_{2}\) and F\(_{3}\) folds in the study area clearly points to (i) Z-shape folds with SE vergence and (ii) a dextral sense of shear motion. Conic section method reveals two dominant families: F\(_{1}\) and F\(_{3}\) folds belong to parabolic shape folds, while F\(_{2}\) folds belong to parabolic shape and hyperbolic shape folds. Ramsay’s scheme emphasizes class 1C (for F\(_{1}\), F\(_{2}\) and F\(_{3}\) folds) and class 3 (for F\(_{2}\) folds) as main fold classes. The co-existence of the various fold shapes can be explained by (i) the structuration of the banded gneiss, (ii) the folding mechanisms that associate shear with a non-least compressive or flattening component in a ductile shear zone and (iii) the change in rheological properties of the band during the period of fold formation. These data allow us to conclude that the Numba region underwent ductile dextral shear and can be integrated (i) in a correlation model with the Central Cameroon Shear Zone (CCSZ) and associated syn-kinematic intrusions and (ii) into the tectonic model of Pan-African belt of central Africa in Cameroon.     

Principle of Formation of the Earth’s Crust under Natural Stress Conditions

This article formulates the experimentally substantiated physical principle that the natural stress condition of the Earth’s crust is formed due to the superposition of stress fields, which is caused by the gravitational forces of the Earth and by tectonic and astrophysical forces that are produced by physical processes in space. The natural stress field is represented by the stress tensor regulatory components: \(\sigma _{z}^{{\text{N}}}\) = \( - \gamma H + {{\sigma }_{{zT}}} + {{\sigma }_{{z{\text{AF}}}}}\), \(\sigma _{x}^{{\text{N}}}\) = \( - \lambda \gamma H + {{\sigma }_{{xT}}} + {{\sigma }_{{x{\text{AF}}}}}\), \(\sigma _{y}^{{\text{N}}}\) = \( - \lambda \gamma H + {{\sigma }_{{yT}}} + {{\sigma }_{{y{\text{AF}}}}}\).     

Rethinking communication in risk interpretation and action

We have studied the attenuation characteristics of eastern Himalaya and southern Tibet by using local earthquake data set that consists of 123 well-located events, recorded by the Himalayan Nepal Tibet Seismic Experiment operated during 2001–2003. We have used single backscattering model to calculate frequency-dependent values of coda Q (\(Q_\mathrm{c}\)). The estimation of \(Q_\mathrm{c}\) is made at central frequencies 2, 4, 8 and 12 Hz through five lapse time windows from 10 to 50 s starting at double the travel time of the S-wave. The observed \(Q_\mathrm{c}\) is found to be strongly frequency-dependent and follows a similar trend as observed in other tectonically active parts of the Himalaya. The trend of variation of \(Q_\mathrm{c}\) with lapse time and the corresponding apparent depths is also studied. Increase in \(Q_\mathrm{c}\) values with the lapse time suggests that the deeper part of the study region is less heterogeneous than the shallower part. The observed values of \(Q_0\) (\(Q_\mathrm{c}\) at 1 Hz) and frequency parameter n indicate that the medium beneath the study area is highly heterogeneous and tectonically very active. A regionalization of the estimated \(Q_0\) is carried out, and a contour map is prepared for the whole region. Some segments of Lesser Himalaya and Sub-Himalaya exhibit very low \(Q_0\) , while the whole Tethyan Himalaya and some parts of Greater Himalaya are characterized by low \(Q_0\) values. Our results are comparable with those obtained from tectonically active regions in the world.     

Petrography and geochemistry of the Proterozoic sandstones of Somanpalli Group from Pomburna area,Eastern Belt of Pranhita–Godavari Valley,central India: Implications for provenance,weathering and tectonic setting

In this paper, we, for the first time, report geochemistry of sandstone from Somanpalli Group from Pomburna area in the Eastern Belt of Pranhita–Godavari (PG) Valley, central India and studied to infer their provenance, intensity of paleo-weathering and depositional tectonic setting. Petrographic study of sandstones show QFL modal composition of arenite. Chemical results show high \(\hbox {SiO}_{2}\) and CIA but lower \(\hbox {Al}_{2}\hbox {O}_{3}, \hbox {TiO}_{2}\), Rb, Sr, \(\hbox {K}_{2}\hbox {O}\) indicating mixed sources. Major elements chemistry parameters such as, \(\hbox {K}_{2}\hbox {O/Al}_{2}\hbox {O}_{3}\) ratio and positive correlation of Rb with \(\hbox {K}_{2}\hbox {O}\), reflects a warm and humid climate for study area. The tectonic discrimination plots (\(\hbox {SiO}_{2}/20\)–\(\hbox {K}_{2}\hbox {O} + \hbox {Na}_{2}\hbox {O}\)–\(\hbox {TiO}_{2} + \hbox {Fe}_{2}\hbox {O}_{3} + \hbox {MgO};\,\hbox {K}_{2}\hbox {O}/\hbox {Na}_{2}\hbox {O}\) vs. \(\hbox {SiO}_{2}\); Th–Sc–Zr/20) indicate dominantly passive margin and slight active tectonic setting. Concentrations of Zr, Nb, Y, and Th are higher compared to the UCC values and the trends in Th/Cr, Th/Co, La/Sc and Cr/Zr ratios support a felsic and mafic source for these sandstones and deposition in passive margin basin. Chondrite normalized REE pattern reflects LREE depletion, negative Eu anomaly and flat HREE similar to UCC, felsic components. ICV value (0.95) also supports tectonically quiescent passive margin settings. CIA values (74) indicate high degree of chemical weathering and warm and humid paleoclimatic condition.     

Near InfraRed Imaging Spectrograph (NIRIS) for ground-based mesospheric OH(6-2) and $$\hbox {O}_{2}$$(0-1) intensity and temperature measurements

This paper describes the development of a new Near InfraRed Imaging Spectrograph (NIRIS) which is capable of simultaneous measurements of OH(6-2) Meinel and \(\hbox {O}_{2}\)(0-1) atmospheric band nightglow emission intensities. In this spectrographic technique, rotational line ratios are obtained to derive temperatures corresponding to the emission altitudes of 87 and 94 km. NIRIS has been commissioned for continuous operation from optical aeronomy observatory, Gurushikhar, Mount Abu (\(24.6^{\circ }\hbox {N}\), \(72.8^{\circ }\hbox {E}\)) since January 2013. NIRIS uses a diffraction grating of 1200 lines \(\hbox {mm}^{-1}\) and 1024\(\times \)1024 pixels thermoelectrically cooled CCD camera and has a large field-of-view (FOV) of \(80^{\circ }\) along the slit orientation. The data analysis methodology adopted for the derivation of mesospheric temperatures is also described in detail. The observed NIRIS temperatures show good correspondence with satellite (SABER) derived temperatures and exhibit both tidal and gravity waves (GW) like features. From the time taken for phase propagation in the emission intensities between these two altitudes, vertical phase speed of gravity waves, \(c_{z}\), is calculated and along with the coherent GW time period ‘\(\tau \)’, the vertical wavelength, \(\lambda _{z}\), is obtained. Using large FOV observations from NIRIS, the meridional wavelengths, \(\lambda _{y}\), are also calculated. We have used one year of data to study the possible cause(s) for the occurrences of mesospheric temperature inversions (MTIs). From the statistics obtained for 234 nights, it appears that in situ chemical heating is mainly responsible for the observed MTIs than the vertical propagation of the waves. Thus, this paper describes a novel near infrared imaging spectrograph, its working principle, data analysis method for deriving OH and \(\hbox {O}_{2}\) emission intensities and the corresponding rotational temperatures at these altitudes, derivation of gravity wave parameters (\(\tau \), \(c_{z}\), \(\lambda _{z}\), and \(\lambda _{y})\), and results on the statistical study of MTIs that exist in the earth’s mesospheric altitudes.