Correlation analyses of effects of temperature on physical and mechanical properties of clay

In determining the physical and mechanical parameters of clay, it is sometimes necessary to determine them indirectly from other parameters since they cannot be measured directly from laboratory or field tests. In order to determine the effect of temperature on the behavior of clay, an indirect approach is used here by analyzing the changes of mass (\(\Delta m\)), density (\(\rho\)), porosity (\(\phi\)), P-wave velocity (\({v_p}\)), thermal conductivity (\(\lambda\)), specific heat capacity (c), resistivity (R) and uniaxial compressive strength (f) of clay from eastern China for a temperature range between 20 and 800 °C. The results indicate that temperature has a significant effect on these parameters. Comparisons between \(\Delta m\) and \(\rho\), \(\Delta m\) and \({v_p}\), \(\rho\) and \({v_p}\), \(\phi\) and \(\lambda\), \({v_p}\) and f, R and f show a linear change among these parameters,whereas the relationships among \(\Delta m\) and \(\phi\), \(\phi\) and \({v_p}\), \(\phi\) and R, \({v_p}\) and \(\lambda\), \(\phi\) and f are exponential. It is difficult to obtain these relationships by using regression analysis with high levels of accuracy. Further refinement is therefore required.     

Electrical conductivity of NaCl-bearing aqueous fluids to 600 °C and 1 GPa

The electrical conductivity of aqueous fluids containing 0.01, 0.1, and 1 M NaCl was measured in an externally heated diamond cell to 600 °C and 1 GPa. These measurements therefore more than double the pressure range of previous data and extend it to higher NaCl concentrations relevant for crustal and mantle fluids. Electrical conductivity was generally found to increase with pressure and fluid salinity. The conductivity increase observed upon variation of NaCl concentration from 0.1 to 1 M was smaller than from 0.01 to 0.1 M, which reflects the reduced degree of dissociation at high NaCl concentration. Measured conductivities can be reproduced (R ² = 0.96) by a numerical model with log \(\sigma\) = ?1.7060– 93.78/T + 0.8075 log c + 3.0781 log \(\rho\) + log \(\varLambda\) ₀(T, \(\rho\)), where \(\sigma\) is the conductivity in S m^?1, T is temperature in K, c is NaCl concentration in wt%, \(\rho\) is the density of pure water (in g/cm³) at given pressure and temperature, and \(\varLambda\) ₀ (T, \(\rho\)) is the molar conductivity of NaCl in water at infinite dilution (in S cm² mol^?1), \(\varLambda\) ₀ = 1573–1212 \(\rho\) + 537 062/T–208 122 721/T ². This model allows accurate predictions of the conductivity of saline fluids throughout most of the crust and upper mantle; it should not be used at temperatures below 100 °C. In general, the data show that already a very small fraction of NaCl-bearing aqueous fluid in the deep crust is sufficient to enhance bulk conductivities to values that would be expected for a high degree of partial melting. Accordingly, aqueous fluids may be distinguished from hydrous melts by comparing magnetotelluric and seismic data. H₂O–NaCl fluids may enhance electrical conductivities in the deep crust with little disturbance of v _p or v _p/v _s ratios. However, at the high temperatures in the mantle wedge above subduction zones, the conductivity of hydrous basaltic melts and saline aqueous fluids is rather similar, so that distinguishing these two phases from conductivity data alone is difficult. Observed conductivities in forearc regions, where temperatures are too low to allow melting, may be accounted for by not more than 1 wt% of an aqueous fluid with 5 wt% NaCl, if this fluid forms a continuous film or fills interconnected tubes.     

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

Middendorfite,K<Subscript>3</Subscript>Na<Subscript>2</Subscript>Mn<Subscript>5</Subscript>Si<Subscript>12</Subscript>(O,OH)<Subscript>36</Subscript> · 2H<Subscript>2</Subscript>O,a new mineral species from the Khibiny pluton,Kola Peninsula

Middendorfite, a new mineral species, has been found in a hydrothermal assemblage in Hilairite hyperperalkaline pegmatite at the Kirovsky Mine, Mount Kukisvumchorr apatite deposit, Khibiny alkaline pluton, Kola Peninsula, Russia. Microcline, sodalite, cancrisilite, aegirine, calcite, natrolite, fluorite, narsarsukite, labuntsovite-Mn, mangan-neptunite, and donnayite are associated minerals. Middendorfite occurs as rhombshaped lamellar and tabular crystals up to 0.1 × 0.2 × 0.4 mm in size, which are combined in worm-and fanlike segregations up to 1 mm in size. The color is dark to bright orange, with a yellowish streak and vitreous luster. The mineral is transparent. The cleavage (001) is perfect, micalike; the fracture is scaly; flakes are flexible but not elastic. The Mohs hardness is 3 to 3.5. Density is 2.60 g/cm³ (meas.) and 2.65 g/cm³ (calc.). Middendorfite is biaxial (?), α = 1.534, β = 1.562, and γ = 1.563; 2V (meas.) = 10°. The mineral is pleochroic strongly from yellowish to colorless on X through brown on Y and to deep brown on Z. Optical orientation: X = c. The chemical composition (electron microprobe, H₂O determined with Penfield method) is as follows (wt %): 4.55 Na₂O, 10.16 K₂O, 0.11 CaO, 0.18 MgO, 24.88 MnO, 0.68 FeO, 0.15 ZnO, 0.20 Al₂O₃, 50.87 SiO₂, 0.17 TiO₂, 0.23 F, 7.73 H₂O; ?O=F₂?0.10, total is 99.81. The empirical formula calculated on the basis of (Si,Al)₁₂(O,OH,F)₃₆ is K_3.04(Na_2.07Ca_0.03)_Σ2.10(Mn_4.95Fe_0.13Mg_0.06Ti_0.03Zn_0.03)_Σ5.20(Si_11.94Al_0.06)_Σ12O_27.57(OH)_8.26F_0.17 · 1.92H₂O. The simplified formula is K₃Na₂Mn₅Si₁₂(O,OH)₃₆ · 2H₂O. Middenforite is monoclinic, space group: P2₁/m or P2₁. The unit cell dimensions are a = 12.55, b = 5.721, c = 26.86 Å; β = 114.04°, V = 1761 ų, Z = 2. The strongest lines in the X-ray powder pattern [d, Å, (I)(hkl)] are: 12.28(100)(002), 4.31(81)(11\(\overline 4 \)), 3.555(62)(301, 212), 3.063(52)(008, 31\(\overline 6 \)), 2.840(90)(312, 021, 30\(\overline 9 \)), 2.634(88)(21\(\overline 9 \), 1.0.\(\overline 1 \)0, 12\(\overline 4 \)), 2.366(76)(22\(\overline 6 \), 3.1.\(\overline 1 \)0, 32\(\overline 3 \)), 2.109(54)(42–33, 42–44, 51\(\overline 9 \), 414), 1.669(64)(2.2.\(\overline 1 \)3, 3.2.\(\overline 1 \)3, 62\(\overline 3 \), 6.1.\(\overline 1 \)3), 1.614(56)(5.0.\(\overline 1 \)6, 137, 333, 71\(\overline 1 \)). The infrared spectrum is given. Middendorfite is a phyllosilicate related to bannisterite, parsenttensite, and the minerals of the ganophyllite and stilpnomelane groups. The new mineral is named in memory of A.F. von Middendorff (1815–1894), an outstanding scientist, who carried out the first mineralogical investigations in the Khibiny pluton. The type material of middenforite has been deposited at the Fersman Mineralogical Museum, Russian Academy of Sciences, Moscow.     

Characteristics of global strong earthquakes and their implications for the present-day stress pattern

Earthquakes occurred on the surface of the Earth contain comprehensive and abundant geodynamic connotations, and can serve as important sources for describing the present-day stress field and regime. An important advantage of the earthquake focal mechanism solution is the ability to obtain the stress pattern information at depth in the lithosphere. During the past several decades, an increasing number of focal mechanisms were available for estimating the present-day stress field and regime. In the present study, altogether 553 focal mechanism data ranging from the year 1976 to 2017 with Mw \(\ge \)7.0 were compiled in the Global/Harvard centroid moment tensor (CMT) catalogue, the characteristics of global strong earthquakes and the present-day stress pattern were analyzed based on these data. The majority of global strong earthquakes are located around the plate boundaries, shallow-focus, and thrust faulting (TF) regime. We grouped 518 of them into 12 regions (Boxes) based on their geographical proximity and tectonic setting. For each box, the present-day stress field and regime were obtained by formal stress inversion. The results indicated that the maximum horizontal principal stress direction was \(\sim \)N–S-trending in western North America continent and southwestern Indonesia, \(\sim \)NNE–SSW-trending in western Middle America and central Asia, \(\sim \)NE–SW in southeastern South America continent and northeastern Australia, \(\sim \)NEE–SWW-trending in western South America continent and southeastern Asia, \(\sim \)E–W-trending in southeastern Australia, and \(\sim \)NW–SE-trending in eastern Asia. The results can provide additional constraints to the driving forces and geodynamic models, allowing them to explain the current plate interactions and crustal tectonic complexities better.     

Geochemical assessment of fluoride enrichment and nitrate contamination in groundwater in hard-rock aquifer by using graphical and statistical methods

This systematic study was carried out with objective to delineate the various sources responsible for \(\hbox {NO}_{3}^{-}\) contamination and \(\hbox {F}^{-}\) enrichment by utilizing statistical and graphical methods. Since Central Ground Water Board, India, indicated susceptibility of \(\hbox {NO}_{3}^{-}\) contamination and \(\hbox {F}^{-}\) enrichment, in most of the groundwater, \(\hbox {NO}_{3}^{-}\) and \(\hbox {F}^{-}\) concentration primarily observed \({>}45\) and \({>}1.5~\hbox {mg/L}\), respectively, i.e., higher than the permissible limit for drinking water. Water Quality Index (WQI) indicates \({\sim }22.81\%\) groundwater are good-water, \({\sim }71.14\%\) groundwater poor-water, \({\sim }5.37\%\) very poor-water and 0.67% unsuitable for drinking purpose. Piper diagram indicates \({\sim }59.73\%\) groundwater hydrogeochemical facies are Ca–Mg–\(\hbox {HCO}_{3 }\) water-types, \({\sim }28.19\%\) Ca–Mg–\(\hbox {SO}_{4}\)–Cl water-types, \({\sim }8.72\%\) Na–K–\(\hbox {SO}_{4}\)–Cl water-types and 3.36% Na–K–\(\hbox {HCO}_{3 }\) water-types. This classification indicates dissolution and mixing are mainly controlling groundwater chemistry. Salinity diagram indicate \({\sim }44.30\%\) groundwater under in low sodium and medium salinity hazard, \({\sim }49.66\%\) groundwater fall under low sodium and high salinity hazard, \({\sim }3.36\%\) groundwater fall under very-high salinity hazard. Sodium adsorption ratio indicates \({\sim }97\%\) groundwater are in excellent condition for irrigation. The spatial distribution of \(\hbox {NO}_{3}^{-}\) indicates significant contribution of fertilizer from agriculture lands. Fluoride enrichment occurs in groundwater through the dissolution of fluoride-rich minerals. By reducing the consumption of fertilizer and stress over groundwater, the water quality can be improved.     

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.     

Providing optimum models for predicting the water loss after harvesting and extracting trees aboveground biomass from the Hyrcanian mixed-beech forests

The ecological and biological attributes of trees stand as well as the water cycle in forests are substantially related to variations of water storage capacity in forest ecosystems. This study aimed to figure out a protocol for monitoring the water storage capacity variations in the Hyrcanian mixed-beech stands after harvesting and extracting trees form the forest. A total of 174 trees were felled and weighed, and destructive sampling following lines of exploitation was carried out for measuring the water content in aboveground biomass of trees. Curve estimation regression analyses including the tree biophysical variables (breast height diameter DBH, total height H, basic wood density \(\rho\)) were used for examining the prediction accuracy. Nonlinear models were log-transformed, and systematic bias was corrected by correction factor depending on standard error of estimate when back-transforming to the originally dependent value. The findings showed that the power-law models were the best functional form for predicting the dependent variables. Using only DBH in the simple power model explained 76% of total variance (Adj.R ² = 0.76) with a low Akaike information criterion (AIC) and normalized root-mean-square error RMSE%, indicating a high accuracy of prediction (\({\text{AIC}} \approx - 238\); RMSE% = 11.7). Adding H and \(\rho\) in the linear log-transformed power models with different interaction terms increased the certainty of prediction with the highest accuracy (\({\text{Adj}}.R^{2} = 0.86,\;{\text{AIC}} = - 329,\;{\text{RMSE}}\% = 9\)). Considering diverse conditions for natural forest sites, the optimum models including the biophysical variables may have associated parameters in the other forests having different stand types and compositions.     

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

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

Moisture rainout fraction over the Indian Ocean during austral summer based on $$^{18}\hbox {O}/{}^{16}\hbox {O}$$ ratios of surface seawater,rainwater at latitude range of 10°N–60°S

Oxygen isotope ratios (\(^{18}\hbox {O}/^{16}\hbox {O}\)) of surface seawater and rainwater samples from the Indian Ocean region (10°N–60°S) during austral summer collected onboard ORV Sagar Nidhi during 2011–2013 have been measured along with salinity, sea surface temperature and relative humidity. The rainwater is isotopically lighter (by \(4.6\pm 2.7\permille )\) compared to the equilibrium condensation of the vapour arising from the seawater at the ambient condition. The isotopic composition of the vapour at high altitude responsible for the rain formation at the sampling location is estimated from a global atmospheric water isotope model (IsoGSM2). The apparent deficit of \(\sim \)5\(\permille \) can be explained by invoking a high degree of rainout (on average, about 70% of the overhead atmospheric moisture) during transport of the source vapour to the sampling location undergoing a Rayleigh fractionation. The required rainout fraction is higher (\(\sim \)80%) in the latitude belt 40°–60°S compared to the equatorial belt (\(\sim \)60%). The pattern of variation in the rainout fraction with latitude is consistent with the well-known evaporation/precipitation processes in the Indian Ocean.     

Micromechanical Modeling of Yield in Isotropic Non-Cohesive Particulate Materials

We perform a micromechanical analysis of general isotropic non-cohesive particulate materials idealized as three-dimensional random assemblies of uniform spheres with a simple linear elastic inter-particle contact force law and inter-particle Coulomb friction law. We obtain analytical relationships between the inter-particle friction coefficient \(\mu\) (or inter-particle friction angle \(\phi _\mu = \tan ^{-1} \mu\)) on the microscale and the material friction angle \(\phi\) on the macroscale. Our micromechanical analysis directly employs force and moment equilibrium (together with compatibility and the contact constitutive assumptions noted) rather than energy methods, and thus can account for the effects of particle rotation, and in particular the effects of mechanisms or zero-energy modes due to particle rotation. To explore the effects of particle rotation, we perform analyses with particle rotation either allowed or prohibited. To validate the analytical results obtained here, we compare the \(\phi\) versus \(\phi _\mu\) curves determined theoretically to those obtained by the discrete element method (DEM) for six randomly packed specimens of 3430–29, 660 uniform spherical elements with uniform inter-element Coulomb friction in Fleischmann et al. in Geotech Geol Eng 32(4):1081–1100, (2014). The \(\phi\) versus \(\phi _\mu\) curves derived here show remarkable agreement with those obtained via DEM simulations in Fleischmann et al. in Geotech Geol Eng 32(4):1081–1100, (2014), especially for the case in which particle rotation is not artificially restrained.     

Variations of trace gases over the Bay of Bengal during the summer monsoon

In situ measurements of near-surface ozone (\(\hbox {O}_{3})\), carbon monoxide (CO), and methane (\(\hbox {CH}_{4})\) were carried out over the Bay of Bengal (BoB) as a part of the Continental Tropical Convergence Zone (CTCZ) campaign during the summer monsoon season of 2009. \(\hbox {O}_{3}\), CO and \(\hbox {CH}_{4}\) mixing ratios varied in the ranges of 8–54 ppbv, 50–200 ppbv and 1.57–2.15 ppmv, respectively during 16 July–17 August 2009. The spatial distribution of mean tropospheric \(\hbox {O}_{3}\) from satellite retrievals is found to be similar to that in surface \(\hbox {O}_{3}\) observations, with higher levels over coastal and northern BoB as compared to central BoB. The comparison of in situ measurements with the Monitoring Atmospheric Composition & Climate (MACC) global reanalysis shows that MACC simulations reproduce the observations with small mean biases of 1.6 ppbv, –2.6 ppbv and 0.07 ppmv for \(\hbox {O}_{3}\), CO and \(\hbox {CH}_{4}\), respectively. The analysis of diurnal variation of \(\hbox {O}_{3}\) based on observations and the simulations from Weather Research and Forecasting coupled with Chemistry (WRF-Chem) at a stationary point over the BoB did not show a net photochemical build up during daytime. Satellite retrievals show limitations in capturing \(\hbox {CH}_{4}\) variations as measured by in situ sample analysis highlighting the need of more shipborne in situ measurements of trace gases over this region during monsoon.     

Indian summer monsoon forcing on the deglacial polar cold reversals

The deglacial transition from the last glacial maximum at \(\sim \)20 kiloyears before present (ka) to the Holocene (11.7 ka to Present) was interrupted by millennial-scale cold reversals, viz., Antarctic Cold Reversal (\(\sim \)14.5–12.8 ka) and Greenland Younger Dryas (\(\sim \)12.8–11.8 ka) which had different timings and extent of cooling in each hemisphere. The cause of this synchronously initiated, but different hemispheric cooling during these cold reversals (Antarctic Cold Reversal \(\sim \)3\(^{\circ }\hbox {C}\) and Younger Dryas \(\sim \)10\(^{\circ }\hbox {C}\)) is elusive because \(\hbox {CO}_{2}\), the fundamental forcing for deglaciation, and Atlantic meridional overturning circulation, the driver of antiphased bipolar climate response, both fail to explain this asymmetry. We use centennial-resolution records of the local surface water \(\delta ^{18}\hbox {O}\) of the Eastern Arabian Sea, which constitutes a proxy for the precipitation associated with the Indian Summer Monsoon, and other tropical precipitation records to deduce the role of tropical forcing in the polar cold reversals. We hypothesize a mechanism for tropical forcing, via the Indian Summer Monsoons, of the polar cold reversals by migration of the Inter-Tropical Convergence Zone and the associated cross-equatorial heat transport.     

Pressuremeter Modulus and Limit Pressure of Clayey Soils Using <Emphasis Type="Italic">GMDH</Emphasis>-Type Neural Network and Genetic Algorithms

Pressuremeter modulus (\(E_{M}\)) and limit pressure (\(P_{L}\)) are used for the calculation of the settlement and bearing capacity of foundation respectively. As the determination of these parameters from pressuremeter test (PMT) is relatively time-consuming and expensive, various empirical correlations have been proposed to correlate the \(E_{M}\) and \(P_{L}\) to other soil parameters. For the existing equations are incapable of estimating these PMT parameters well, in present research group method of data handling type neural network is used to estimate the \(E_{M}\) and \(P_{L}\) of clayey soils. The \(E_{M}\) and \(P_{L}\) were modeled as a function of three variables including the moisture content (\(\omega\)), plasticity index and corrected SPT blow counts (\(N_{60}\)). A database containing 51 data sets have been used for training and testing of the models. The performances of proposed models are compared with those of existing empirical equations. The results demonstrate that appreciable improvement with respect to the other correlations has been achieved. At the end, sensitivity analysis of the obtained models has been performed to study the influence of input parameters on model outputs and shows that the \(N_{60}\) is the most influential parameter on the PMT parameters.     

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

Lithologic boundaries from gravity and magnetic anomalies over Proterozoic Dalma volcanics

Dalma volcanics (DVs) has intruded the older Singhbhum Group of Metapelites. Despite DVs being rich in mineralisation, its boundaries are not clearly demarcated. Gravity and magnetic surveys have been attempted for mapping the boundaries in DVs. These surveys were made in the northern fringe of the DVs over an area of \(\sim \)0.70 \(\hbox {km}^{2}\) along 13 parallel lines at 50 m spacing. The data was acquired at \(\sim \)25 \(\hbox {m}\) spacing. The surveys were taken for determination of lithological boundaries, depths and nature of causative source using Euler depth solutions and radially averaged power spectrum (RAPS). Residual anomaly maps of gravity and magnetic intensity show the same trend as that of Bouguer gravity anomaly and total magnetic intensity anomaly map indicating towards shallow sources. The magnetic map in general follows the same pattern as that of gravity anomaly maps. The map shows coincident high gravity and magnetic anomalies. These anomalies together with resistivity signatures confirm that the northern fringe of DVs hosts volcanogenic massive sulphide settings. The Euler depth solution delineated the lateral boundaries and nature of the source. It seems that the source is of spherical nature lying within a depth range of 25–40 m. The obtained lithological (vertical) units from RAPS are between Lower DVs, Upper DVs and Singhbhum Group Metapelites at depths of \(\sim \)15, \(~\sim \)25 and \(\sim \)40 \(\hbox {m}\), respectively. The metallogeny is associated with the Upper DVs and the corresponding delineated lithological (vertical) unit is indicative of the top of the ore body. Good agreement is observed with the geological succession from the drilling data and resistivity data. The findings suggest that the northern fringe of DVs could be a preferred target for drilling.     

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.     

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.     

Phase diagram and thermodynamic properties of AIPO<Subscript>4</Subscript> based on first-principles calculations and the quasiharmonic approximation

We calculated the phase diagram of \(\hbox {AlPO}_{4}\) up to 15 GPa and 2,000 K and investigated the thermodynamic properties of the high-pressure phases. The investigated phases include the berlinite, moganite-like, \(\hbox {AlVO}_{4},\, P2_1/c\), and \(\hbox {CrVO}_{4}\) phases. The computational methods used include density functional theory, density functional perturbation theory, and the quasiharmonic approximation. The investigated thermodynamic properties include the thermal equation of state, isothermal bulk modulus, thermal expansivity, and heat capacity. With increasing pressure, the ambient phase berlinite transforms to the moganite-like phase, and then to the \(\hbox {AlVO}_{4}\) and \(P2_1/c\) phases, and further to the \(\hbox {CrVO}_{4}\) phase. The stability fields of the \(\hbox {AlVO}_{4}\) and \(P2_1/c\) phases are similar in pressure but different in temperature, as the \(\hbox {AlVO}_{4}\) phase is stable at low temperatures, whereas the \(P2_1/c\) phase is stable at high temperatures. All of the phase relationships agree well with those obtained by quench experiments, and they support the stabilities of the moganite-like, \(\hbox {AlVO}_{4}\), and \(P2_1/c\) phases, which were not observed in room-temperature compression experiments.