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.     

Laboratory investigations of inert gas flow behaviors in compact sandstone

The seepage evolution behavior of compact rock is significant for the stability and safety of many engineering applications. In this research, both hydrostatic and triaxial compression tests were conducted on compact sandstone using an inert gas, namely argon. A triaxial compression test with a water permeability measurement was carried out to study the difference between the gas permeability and water permeability evolutions during the complete stress–strain process. Based on the experimental data, the hydrostatic stress-dependent gas permeability was discussed firstly. A second-order function was proposed to predict and explain the gas slippage effect. The mechanical properties and crack development of the sandstone samples were discussed to better understand the permeability evolution with crack growth during the complete stress–strain process. The results show that the gas permeability evolution can be divided into five stages according to the different crack growth stages. Then, the permeability changes in the crack closure stress \( \sigma_{\text{cc}} \), crack initiation stress \( \sigma_{\text{ci}} \), crack damage stress \( \sigma_{\text{cd}} \) and peak stress \( \sigma_{\text{p}} \) with confining pressures were analyzed. Finally, we found that the difference between the corrected gas permeability and water permeability can be attributed to the interaction between the water and sandstone grains.     

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.     

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.     

Arsenic in shallow groundwater of Bangladesh: investigations from three different physiographic settings

Occurrences of arsenic (As) in the Bengal Basin of Bangladesh show close relationships with depositional environments and sediment textures. Hydrochemical data from three sites with varying physiography and sedimentation history show marked variations in redox status and dissolved As concentrations. Arsenic concentration in groundwater of the Ganges Flood Plain (GFP) is characteristically low, where high Mn concentrations indicate redox buffering by reduction of Mn(IV)-oxyhydroxides. Low DOC, \( {\text{HCO}}^{ - }_{3} \), \( {\text{NH}}^{ + }_{4} \) and high \( {\text{NO}}^{ - }_{3} \) and \( {\text{SO}}^{{2 - }}_{4} \) concentrations reflect an elevated redox status in GFP aquifers. In contrast, As concentration in the Ganges Delta Plain (GDP) is very high along with high Fe and low Mn. In the Meghna Flood Plain (MFP), moderate to high As and Fe concentrations and low Mn are detected. Degradation of organic matter probably drives redox reactions in the aquifers, particularly in MFP and GDP, thereby mobilising dissolved As. Speciation calculations indicate supersaturation with respect to siderite and vivianite in the groundwater samples at MFP and GDP, but groundwater in the GFP wells is generally supersaturated with respect to rhodochrosite. Values of log P_CO2 at MFP and GDP sites are generally higher than at the GFP site. This is consistent with Mn(IV)-redox buffering suggested at the GFP site compared to Fe(III)-redox buffering at MFP and GDP sites.     

The preference of nitrate uptake in Chinese prickly ash estimated by δ<Superscript>15</Superscript>N values and cation concentrations

The nitrogen isotopic compositions of plant tissue could reflect its uptake of and preference for ammonium or nitrate. However, various factors may influence the field-collected δ¹⁵N values under field condition, which causes the interpretation problematic. The spatial variation of nitrogen (N) concentrations and the isotopic compositions were investigated in the soils and tissues of Chinese prickly ash from the southwest China to the east China. The objectives were to investigate the variation in soil and tissue δ¹⁵N values and N forms taken up by the plant. The leaf and root δ¹⁵N values varied significantly in response to the pattern of soil δ¹⁵N values. The difference in δ¹⁵N values between the leaves and roots was 2.57‰ and may be caused by an increase in the transport of unassimilated \( {\text{NO}}_{3}^{ - } \) and \( {\text{NH}}_{4}^{ + } \) to the leaves. Leaf nitrogen was significantly and positively correlated with leaf potassium and negatively related to leaf calcium. Because potassium is the favoured counter-cation for nitrate transport in the xylem, the enrichment of ¹⁵N in leaf relative to root induced by preferenced uptake of nitrate should be accompanied by significant and positive relationship of leaf nitrogen with leaf potassium concentrations. These results suggest that Chinese prickly ash prefers \( {\text{NO}}_{3}^{ - } \) over \( {\text{NH}}_{4}^{ + } \).     

The high-pressure stability of chlorite and other hydrates in subduction mélanges: experiments in the system Cr<Subscript>2</Subscript>O<Subscript>3</Subscript>–MgO–Al<Subscript>2</Subscript>O<Subscript>3</Subscript>–SiO<Subscript>2</Subscript>–H<Subscript>2</Subscript>O

The solubility of chromium in chlorite as a function of pressure, temperature, and bulk composition was investigated in the system Cr₂O₃–MgO–Al₂O₃–SiO₂–H₂O, and its effect on phase relations evaluated. Three different compositions with X _Cr = Cr/(Cr + Al) = 0.075, 0.25, and 0.5 respectively, were investigated at 1.5–6.5 GPa, 650–900 °C. Cr-chlorite only occurs in the bulk composition with X _Cr = 0.075; otherwise, spinel and garnet are the major aluminous phases. In the experiments, Cr-chlorite coexists with enstatite up to 3.5 GPa, 800–850 °C, and with forsterite, pyrope, and spinel at higher pressure. At P > 5 GPa other hydrates occur: a Cr-bearing phase-HAPY (Mg_2.2Al_1.5Cr_0.1Si_1.1O₆(OH)₂) is stable in assemblage with pyrope, forsterite, and spinel; Mg-sursassite coexists at 6.0 GPa, 650 °C with forsterite and spinel and a new Cr-bearing phase, named 11.5 Å phase (Mg:Al:Si = 6.3:1.2:2.4) after the first diffraction peak observed in high-resolution X-ray diffraction pattern. Cr affects the stability of chlorite by shifting its breakdown reactions toward higher temperature, but Cr solubility at high pressure is reduced compared with the solubility observed in low-pressure occurrences in hydrothermal environments. Chromium partitions generally according to \(X_{\text{Cr}}^{\text{spinel}}\) ? \(X_{\text{Cr}}^{\text{opx}}\) > \(X_{\text{Cr}}^{\text{chlorite}}\) ≥ \(X_{\text{Cr}}^{\text{HAPY}}\) > \(X_{\text{Cr}}^{\text{garnet}}\). At 5 GPa, 750 °C (bulk with X _Cr = 0.075) equilibrium values are \(X_{\text{Cr}}^{\text{spinel}}\) = 0.27, \(X_{\text{Cr}}^{\text{chlorite}}\) = 0.08, \(X_{\text{Cr}}^{\text{garnet}}\) = 0.05; at 5.4 GPa, 720 °C \(X_{\text{Cr}}^{\text{spinel}}\) = 0.33, \(X_{\text{Cr}}^{\text{HAPY}}\) = 0.06, and \(X_{\text{Cr}}^{\text{garnet}}\) = 0.04; and at 3.5 GPa, 850 °C \(X_{\text{Cr}}^{\text{opx}}\) = 0.12 and \(X_{\text{Cr}}^{\text{chlorite}}\) = 0.07. Results on Cr–Al partitioning between spinel and garnet suggest that at low temperature the spinel- to garnet-peridotite transition has a negative slope of 0.5 GPa/100 °C. The formation of phase-HAPY, in assemblage with garnet and spinel, at pressures above chlorite breakdown, provides a viable mechanism to promote H₂O transport in metasomatized ultramafic mélanges of subduction channels.     

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.     

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.     

Effect of solar flare on the equatorial electrojet in eastern Brazil region

The effect of solar flare, sudden commencement of magnetic storm and of the disturbances ring current on the equatorial electrojet in the Eastern Brazil region, where the ground magnetic declination is as large as \(20^{^{\circ }}\hbox {W}\) is studied based on geomagnetic data with one minute resolution from Bacabal during November–December 1990. It is shown that the mean diurnal vector of the horizontal field was aligned along \(2{^{\circ }}\hbox {E}\) of north at Huancayo and \(30{^{\circ }}\hbox {W}\) of north at Bacabal during the month of December 1990. Number of solar flares that occurred on 30 December 1990 indicated the direction of solar flare related \(\Delta H\) vector to be aligned along \(5{^{\circ }}\hbox {E}\) of north at Huancayo and \(28{^{\circ }}\hbox {W}\) of north at Bacabal. This is expected as the solar flare effects are due to the enhanced conductivity in the ionosphere. The SC at 2230 UT on 26 November 1990 produced a positive impulse in \(\Delta X\) and negative impulse in \(\Delta Y\) at Bacabal with \(\Delta H\) vector aligned along \(27{^{\circ }}\hbox {W}\) of north. At Huancayo the \(\Delta H\) vector associated with SC is aligned along \(8{^{\circ }}\hbox {E}\) of north, few degrees east to the alignment of the diurnal vector of H. The magnetic storm that followed the SC had a minimum Dst index of –150 nT. The corresponding storm time disturbance in \(\Delta X\) at Huancayo as well as at Bacabal were about –250 nT but \(\Delta Y\) at Bacabal was about +70 nT and very small at Huancayo, that give the alignment of the H vector due to ring current about \(16{^{\circ }}\hbox {W}\) of north at Bacabal and almost along N–S at Huancayo. Thus alignment of the \(\Delta H\) vector due to ring current at Bacabal is \(14{^{\circ }}\hbox {E}\) of the mean direction of \(\Delta H\) vector during December 1990. This is consistent with the direction of ring current dependent on the dipole declination at the ring current altitude which is about \(5{^{\circ }}\hbox {W}\) of north over Bacabal and the deviation of declination due to the ring current during disturbed period given by the angle (\(\psi \)-D).     

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}}}}}\).     

Comparison of earthquake source characteristics in the Kachchh Rift Basin and Saurashtra horst,Deccan Volcanic Province,western India

Seismic source parameters of small to moderate sized intraplate earthquakes that occurred during 2002–2009 in the tectonic blocks of Kachchh Rift Basin (KRB) and the Saurashtra Horst (SH), in the stable continental region of western peninsular India, are studied through spectral analysis of shear waves. The data of aftershock sequence of the 2001 Bhuj earthquake (\(M_{w}\) 7.7) in the KRB and the 2007 Talala earthquake (\(M_{w}\) 5.0) in the SH are used for this study. In the SH, the seismic moment (\(M_{o})\), corner frequency \((f_{c})\), stress drop (\(\varDelta \sigma \)) and source radius (r) vary from \(7.8\times 10^{11}\) to \(4.0\times \)10\(^{16}\) N-m, 1.0–8.9 Hz, 4.8–10.2 MPa and 195–1480 m, respectively. While in the KRB, these parameters vary from \(M_{o} \sim 1.24 \,\times \, 10^{11}\) to \(4.1 \times 10^{16}\) N-m, \(f_{c }\sim \) 1.6 to 13.1 Hz, \(\varDelta \sigma \sim 0.06\) to 16.62 MPa and \(r \sim 100\) to 840 m. The kappa (K) value in the KRB (0.025–0.03) is slightly larger than that in the SH region (0.02), probably due to thick sedimentary layers. The estimated stress drops of earthquakes in the KRB are relatively higher than those in SH, due to large crustal stress concentration associated with mafic/ultramafic rocks at the hypocentral depths. The results also suggest that the stress drop value of intraplate earthquakes is larger than the interplate earthquakes. In addition, it is observed that the strike-slip events in the SH have lower stress drops, compared to the thrust and strike-slip events.     

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).     

Synthesis of monticellite–forsterite and merwinite–forsterite symplectites in the CaO–MgO–SiO<Subscript>2</Subscript> model system: influence of temperature and water content on microstructure evolution

Homogeneous single crystals of synthetic monticellite with the composition \({\text{Ca}}_{0.88}{\text{Mg}}_{1.12}{\text{SiO}}_4\) (Mtc I) were annealed in a piston-cylinder apparatus at temperatures between 1000 and \(1200\,^{\circ }\hbox {C}\), pressures of 1.0–1.4 GPa, for run durations from 10 min to 24 h and applying bulk water contents ranging from 0.0 to 0.5 wt% of the total charge. At these conditions, Mtc I breaks down to a fine-grained, symplectic intergrowth. Thereby, two types of symplectites are produced: a first symplectite type (Sy I) is represented by an aggregate of rod-shaped forsterite immersed in a matrix of monticellite with end-member composition (Mtc II), and a second symplectite type (Sy II) takes the form of a lamellar merwinite–forsterite intergrowth. Both symplectites may form simultaneously, where the formation of Sy I is favoured by the presence of water. Sy I is metastable with respect to Sy II and is successively replaced by the latter. For both symplectite types, the characteristic spacing of the symplectite phases is independent of run duration and is only weeakly influenced by the water content, but it is strongly temperature dependent. It varies from about 400 nm at \(1000\,^{\circ }\hbox {C}\) to 1200 nm at \(1100\,^{\circ }\hbox {C}\) in Sy I, and from 300 nm at \(1000\,^{\circ }\hbox {C}\) to 700 nm at \(1200\,^{\circ }\hbox {C}\) in Sy II. A thermodynamic analysis reveals that the temperature dependence of the characteristic spacing of the symplectite phases is due to a relatively high activation energy for chemical segregation by diffusion within the reaction front as compared to the activation energy for interface reactions at the reaction front. The temperature dependence of the characteristic lamellar spacing and the temperature-time dependence of overall reaction progress have potential for applications in geo-thermometry and geo-speedometry.     

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.     

Ammonium Silicate Stability Relations

A method has been developed to control ammonium fugacity, \(f_{{\text{NH}}_{3}}\), at elevated temperatures and pressures. The method uses an internal nitrogen buffer, the assemblage Cr + CrN, in conjunction with a traditional external hydrogen buffer. In this manner, all gas fugacities in the system N-O-H can be calculated.The Cr + CrN buffer has been applied to study equilibria between buddingtonite (ammonium feldspar), ammonium muscovite, sillimanite, and quartz at a constant gas pressure of 2,000 bars. Two of the five relevant reactions were measured experimentally; from these data, it is possible to calculate isothermal sections at 500, 600, and 700° C.Below 600° C, ammonium muscovite is stable even at extremely low levels of \(f_{{\text{NH}}_{3}}\), while buddingtonite requires \(f_{{\text{NH}}_{3}}\;\geqq\;10^4\) bars. Release of NH₃ during progressive metamorphism can be achieved by three processes: thermal decomposition, dehydration, and cation exchange. Within the crust, \(f_{{\text{NH}}_{3}}\) predominates over \(f_{{\text{N}}_{2}}\) by several orders of magnitude; but on the surface, nitrogen released as NH₃ by metamorphism will be oxidized to N₂. Biological materials provide important intermediate storage for nitrogen compounds during the nitrogen cycle.     

Assessment of spatial distribution of contaminants and their levels in soil and water resources of Calabar,Nigeria using geophysical and geological data

The contamination levels of soils and water resources in Calabar, Nigeria have been investigated using resistivity (vertical electrical sounding and electrical resistivity tomography), geochemical analyses of soil and water resources and textural analysis. Sixty randomly sited VES sites were investigated in two seasons while ERT investigations were performed along four profiles. The geochemical investigations were spread across seasons in order to track seasonal changes in physico-chemical parameters: hydrogen ion concentration (pH), electrical conductivity, total dissolved solids, chloride ion (Cl^?), nitrate ion (\( {\text{NO}}_{ 3}^{ - } \)), bicarbonate (\( {\text{HCO}}_{ 3}^{ - } \)), sulphate ion (\( {\text{SO}}_{ 4}^{2 - } \)), calcium ion (Ca²⁺), sodium ion (Na⁺), potassium ion (K⁺) and magnesium ion (Mg²⁺). Additionally, concentrations of ammonium, aluminium and nitrite ions in soils were determined. Results show that ionic concentrations in the sand-dominated soils and water are within permissible limits and baseline standards. The resistivities follow known trends in the area. However, at the central waste disposal site, a localised thin (< 5 m), low resistivity (< 15 Ωm) anomaly suspected to be due to contamination by leachates was observed. Comparatively, the contaminated area is also characterised by marginal increase in ionic concentrations. Strong attenuation capacities of overlying and adjoining clay/lateritic sediments and optimal design of the waste dump site probably reduced the spread of contaminants. The contaminated zone need to be closely monitored so that it does not extend to the aquifers. Hence, all strategies presently being used in managing wastes in Calabar should be sustained.     

Experimental determination of liquidus H<Subscript>2</Subscript>O contents of haplogranite at deep-crustal conditions

The liquidus water content of a haplogranite melt at high pressure (P) and temperature (T) is important, because it is a key parameter for constraining the volume of granite that could be produced by melting of the deep crust. Previous estimates based on melting experiments at low P (≤0.5 GPa) show substantial scatter when extrapolated to deep crustal P and T (700–1000 °C, 0.6–1.5 GPa). To improve the high-P constraints on H₂O concentration at the granite liquidus, we performed experiments in a piston–cylinder apparatus at 1.0 GPa using a range of haplogranite compositions in the albite (Ab: NaAlSi₃O₈)—orthoclase (Or: KAlSi₃O₈)—quartz (Qz: SiO₂)—H₂O system. We used equal weight fractions of the feldspar components and varied the Qz between 20 and 30 wt%. In each experiment, synthetic granitic composition glass + H₂O was homogenized well above the liquidus T, and T was lowered by increments until quartz and alkali feldspar crystalized from the liquid. To establish reversed equilibrium, we crystallized the homogenized melt at the lower T and then raised T until we found that the crystalline phases were completely resorbed into the liquid. The reversed liquidus minimum temperatures at 3.0, 4.1, 5.8, 8.0, and 12.0 wt% H₂O are 935–985, 875–900, 775–800, 725–775, and 650–675 °C, respectively. Quenched charges were analyzed by petrographic microscope, scanning electron microscope (SEM), X-ray diffraction (XRD), and electron microprobe analysis (EMPA). The equation for the reversed haplogranite liquidus minimum curve for Ab_36.25Or_36.25Qz_27.5 (wt% basis) at 1.0 GPa is \(T = - 0.0995 w_{{{\text{H}}_{ 2} {\text{O}}}}^{ 3} + 5.0242w_{{{\text{H}}_{ 2} {\text{O}}}}^{ 2} - 88.183 w_{{{\text{H}}_{ 2} {\text{O}}}} + 1171.0\) for \(0 \le w_{{{\text{H}}_{ 2} {\text{O}}}} \le 17\) wt% and \(T\) is in °C. We present a revised \(P - T\) diagram of liquidus minimum H₂O isopleths which integrates data from previous determinations of vapor-saturated melting and the lower pressure vapor-undersaturated melting studies conducted by other workers on the haplogranite system. For lower H₂O (<5.8 wt%) and higher temperature, our results plot on the high end of the extrapolated water contents at liquidus minima when compared to the previous estimates. As a consequence, amounts of metaluminous granites that can be produced from lower crustal biotite–amphibole gneisses by dehydration melting are more restricted than previously thought.     

Trend analysis of evapotranspiration over India: Observed from long-term satellite measurements

Owing to the lack of consistent spatial time series data on actual evapotranspiration (ET), very few studies have been conducted on the long-term trend and variability in ET at a national scale over the Indian subcontinent. The present study uses biome specific ET data derived from NOAA satellite’s advanced very high resolution radiometer to investigate the trends and variability in ET over India from 1983 to 2006. Trend analysis using the non-parametric Mann–Kendall test showed that the domain average ET decreased during the period at a rate of \(0.22\,\hbox {mm year}^{-1}\). A strong decreasing trend (\(m = -1.75\, \hbox {mm year}^{-1}\), \(F = 17.41\), \(P\) 0.01) was observed in forest regions. Seasonal analyses indicated a decreasing trend during southwest summer monsoon (\(m= -0.320\, \hbox {mm season}^{-1}\,\hbox {year}^{-1})\) and post-monsoon period (\(m= -0.188\, \hbox {mm season}^{-1 }\,\hbox {year}^{-1})\). In contrast, an increasing trend was observed during northeast winter monsoon (\(m = 0.156 \,\hbox {mm season}^{-1 }\,\hbox {year}^{-1})\) and pre-monsoon (\(m = 0.068\, \hbox {mm season}^{-1 }\,\hbox {year}^{-1})\) periods. Despite an overall net decline in the country, a considerable increase ( \(4 \,\hbox {mm year}^{-1}\)) was observed over arid and semi-arid regions. Grid level correlation with various climatic parameters exhibited a strong positive correlation (\(r \!>\!0.5\)) of ET with soil moisture and precipitation over semi-arid and arid regions, whereas a negative correlation (\(r\) \(-0.5\)) occurred with temperature and insolation in dry regions of western India. The results of this analysis are useful for understanding regional ET dynamics and its relationship with various climatic parameters over India. Future studies on the effects of ET changes on the hydrological cycle, carbon cycle, and energy partitioning are needed to account for the feedbacks to the climate.     

Heat capacity of hydrous trachybasalt from Mt Etna: comparison with CaAl<Subscript>2</Subscript>Si<Subscript>2</Subscript>O<Subscript>8</Subscript> (An)–CaMgSi<Subscript>2</Subscript>O<Subscript>6</Subscript> (Di) as basaltic proxy compositions

The specific heat capacity (C _p) of six variably hydrated (~3.5 wt% H₂O) iron-bearing Etna trachybasaltic glasses and liquids has been measured using differential scanning calorimetry from room temperature across the glass transition region. These data are compared to heat capacity measurements on thirteen melt compositions in the iron-free anorthite (An)–diopside (Di) system over a similar range of H₂O contents. These data extend considerably the published C _p measurements for hydrous melts and glasses. The results for the Etna trachybasalts show nonlinear variations in, both, the heat capacity of the glass at the onset of the glass transition (i.e., C _p ^g ) and the fully relaxed liquid (i.e., C _p ^l ) with increasing H₂O content. Similarly, the “configurational heat capacity” (i.e., C _p ^c  = C _p ^l  ? C _p ^g ) varies nonlinearly with H₂O content. The An–Di hydrous compositions investigated show similar trends, with C _p values varying as a function of melt composition and H₂O content. The results show that values in hydrous C _p ^g , C _p ^l and C _p ^c in the depolymerized glasses and liquids are substantially different from those observed for more polymerized hydrous albitic, leucogranitic, trachytic and phonolitic multicomponent compositions previously investigated. Polymerized melts have lower C _p ^l and C _p ^c and higher C _p ^g with respect to more depolymerized compositions. The covariation between C _p values and the degree of polymerization in glasses and melts is well described in terms of SM_hydrous and NBO/T _hydrous. Values of C _p ^c increase sharply with increasing depolymerization up to SM_hydrous ~ 30–35 mol% (NBO/T _hydrous ~ 0.5) and then stabilize to an almost constant value. The partial molar heat capacity of H₂O for both glasses (\( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{g}} \)) and liquids (\( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{l}} \)) appears to be independent of composition and, assuming ideal mixing, we obtain a value for \( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{l}} \) of 79 J mol^?1 K^?1. However, we note that a range of values for \( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{l}} \) (i.e., ~78–87 J mol^?1 K^?1) proposed by previous workers will reproduce the extended data to within experimental uncertainty. Our analysis suggests that more data are required in order to ascribe a compositional dependence (i.e., nonideal mixing) to \( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{l}} \).