Modeling competing hydraulic fracture propagation with the extended finite element method

We present an extended finite element framework to numerically study competing hydraulic fracture propagation. The framework is capable of modeling fully coupled hydraulic fracturing processes including fracture propagation, elastoplastic bulk deformation and fluid flow inside both fractures and the wellbore. In particular, the framework incorporates the classical orifice equation to capture fluid pressure loss across perforation clusters linking the wellbore with fractures. Dynamic fluid partitioning among fractures during propagation is solved together with other coupled factors, such as wellbore pressure loss (\(\Delta p_w\)), perforation pressure loss (\(\Delta p\)), interaction stress (\(\sigma _\mathrm{int}\)) and fracture propagation. By numerical examples, we study the effects of perforation pressure loss and wellbore pressure loss on competing fracture propagation under plane-strain conditions. Two dimensionless parameters \(\Gamma = \sigma _\mathrm{int}/\Delta p\) and \(\Lambda = \Delta p_w/\Delta p\) are used to describe the transition from uniform fracture propagation to preferential fracture propagation. The numerical examples demonstrate the dimensionless parameter \(\Gamma \) also works in the elastoplastic media.     

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.     

Hydrocarbon source rock potential evaluation of the Late Paleocene Patala Formation,Salt Range,Pakistan: Organic geochemical and palynofacies approach

Organic geochemical and palynofacies analyses were carried out on shale intervals of the Late Paleocene Patala Formation at Nammal Gorge Section, western Salt Range, Pakistan. The total organic carbon content and Rock-Eval pyrolysis results indicated that the formation is dominated by type II and type III kerogens. Rock-Eval \({T}_{\mathrm{max}}\) vs. hydrogen index (HI) and thermal alteration index indicated that the analysed shale intervals present in the formation are thermally mature. \(S_{1}\) and \(S_{2}\) yields showed poor source rock potential for the formation. Three palynofacies assemblages including palynofacies-1, palynofacies-2 and palynofacies-3 were identified, which are prone to dry gas, wet gas and oil generation, respectively. The palynofacies assessment revealed the presence of oil/gas and gas prone type II and type III kerogens in the formation and their deposition on proximal shelf with suboxic to anoxic conditions. The kerogen macerals are dominated by vitrinite and amorphinite with minor inertinite and liptinite. The kerogen macerals are of both marine and terrestrial origin, deposited on a shallow shelf. Overall, the dark black carbonaceous shales present within the formation act as a source rock for hydrocarbons with poor-to-moderate source rock quality, while the grey shales act as a poor source rock for hydrocarbon generation.     

Short-temporal variation of soil organic carbon in different land use systems in the Ramsar site 2027 ‘Presa Manuel Ávila Camacho’ Puebla

The soil organic carbon (SOC) was determined in 40 sites at two depths (0–10 and 10–20 cm) for different uses of soil during one year (February 2014–February 2015). The total SOC stored in the analysed soil from the Ramsar site was \(9{,}67 \times 10^{5}\) t, from which 40% was stored in induced pasture, followed by the red oak forest with shrubbery secondary vegetation, rain-fed agriculture and human settlements (24%, 23%, and 13%, respectively); the last was evaluated to determine how the proximity of the city impacts the SOC. The SOC concentrations present significant differences with respect to soil depth (\(p=0.0\)) and land use (\(p=0.0\)). The temporal distribution maps showed that SOC did not present significant variations in the short-term. A relation between SOC and bulk density was found (\(r = -0.654\), \(p=0.00\)), with respect to other physicochemical properties. Moreover, a significant relation between SOC and stored total nitrogen (\(r = 0.585\); \(p = 0.00\)) was found. This work represents the first study that analyses the current condition of the soils in the Ramsar site ‘Presa Manuel Ávila Camacho’.     

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

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

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.     

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.     

Effect of $$\hbox {CO}_2$$ solubility on dissolution rate of calcite in saline aquifers for temperature range of 50–100 $$^{\circ } \hbox {C}$$∘C and pressures up to 600 bar: alterations of fractures geometry in carbonate rocks by $$\hbox {CO}_2$$-acidified brines

An empirical model is developed to predict the dissolution rate of calcite in saline solutions that are saturated with respect to dissolved \(\hbox {CO}_2\) over a broad range of both subcritical and supercritical conditions. The focus is on determining the rate of calcite dissolution within a temperature range of 50–100 \(^\circ \hbox {C}\) and pressures up to 600 bar, relevant for \(\hbox {CO}_2\) sequestration in saline aquifers. A general reaction kinetic model is used that is based on the extension of the standard Arrhenius equation with an added, solubility-dependent, pH term to account for the saturated concentration of dissolved \(\hbox {CO}_2\). On the basis of this kinetic model, a new rate equation is obtained using multi-parameter, nonlinear regression of experimental data to determine the dissolution of calcite as a function of temperature, pressure and salinity. Different models for the activity coefficient of \(\hbox {CO}_2\) dissolved in saline solutions are accounted for. The new rate equation helps us obtain good agreement with experimental data, and it is applied to study the geochemically induced alterations of fracture geometry due to calcite dissolution.     

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.     

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.     

Thermal equation of state of CaIrO<Subscript>3</Subscript> post-perovskite

The pressure–volume–temperature (P–V–T) relation of CaIrO₃ post-perovskite (ppv) was measured at pressures and temperatures up to 8.6 GPa and 1,273 K, respectively, with energy-dispersive synchrotron X-ray diffraction using a DIA-type, cubic-anvil apparatus (SAM85). Unit-cell dimensions were derived from the Le Bail full profile refinement technique, and the results were fitted using the third-order Birth-Murnaghan equation of state. The derived bulk modulus \( K_{T0} \) at ambient pressure and temperature is 168.3 ± 7.1 GPa with a pressure derivative \( K_{T0}^{\prime } \) = 5.4 ± 0.7. All of the high temperature data, combined with previous experimental data, are fitted using the high-temperature Birch-Murnaghan equation of state, the thermal pressure approach, and the Mie-Grüneisen-Debye formalism. The refined thermoelastic parameters for CaIrO₃ ppv are: temperature derivative of bulk modulus \( (\partial K_{T} /\partial T)_{P} \) = ?0.038 ± 0.011 GPa K^?1, \( \alpha K_{T} \) = 0.0039 ± 0.0001 GPa K^?1, \( \left( {\partial K_{T} /\partial T} \right)_{V} \) = ?0.012 ± 0.002 GPa K^?1, and \( \left( {\partial^{2} P/\partial T^{2} } \right)_{V} \) = 1.9 ± 0.3 × 10^?6 GPa² K^?2. Using the Mie-Grüneisen-Debye formalism, we obtain Grüneisen parameter \( \gamma_{0} \) = 0.92 ± 0.01 and its volume dependence q = 3.4 ± 0.6. The systematic variation of bulk moduli for several oxide post-perovskites can be described approximately by the relationship K _T0  = 5406.0/V(molar) + 5.9 GPa.     

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.     

Shear Behaviors of the Binary Structural Plane

Binary-medium contact interfaces widely exist in rock engineering. They have significant impacts on the safety of rock engineering due to their poor shear behavior. A material of different strength is produced by pouring mortar of a different sand-to-cement ratio on the top of a rock-like mortar material (with the ratio of 1:1), thereby forming a binary-medium structural plane. Then, direct shear test is performed on the structural plane by applying different normal stresses. The shear strength parameters of the structural plane (cohesion \(c\) and friction angle \(\varphi\)) are obtained from the Mohr–Coulomb criterion. Moreover, the mechanical behaviors of the structural plane are compared with the unitary-medium specimen in the shearing process. A similar shear stress–shear displacement rule is observed in the shearing process. However, the peak and residual shear strengths of the binary structural plane are far lower than those of the unitary ones. The difference between the unitary and binary planes at cohesion \(c\) decreases with the increasing sand–cement ratio, whereas a up-down trend is observed in the friction angle \(\varphi\) with the increasing sand-to-cement ratio. When the upper and lower parts of the structural plane are different in sand-to-cement ratio, the cohesion \(c\) of the structural plane slightly increases with the increasing ratio of the upper specimen. However, when the two parts are identical in ratio, the cohesion of the structural plane reaches the peak, and its friction angle \(\varphi\) substantially increases with the increasing ratio of the upper part.     

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.     

Characterization of organic carbon in black shales of the Kachchh basin,Gujarat, India

Thirty-three black shale samples from four locations on the onland Kachchh basin, western India were analyzed to characterize organic carbon (OC), thermal maturity and to determine the hydrocarbon potential of the basin. Upper Jurassic black shales from the Jhuran Formation (Dhonsa and Kodki areas) are characterized by the presence of chlorite, halloysite, high \(T_{\mathrm{max}}\), low OC, low hydrogen index and high oxygen index. These parameters indicate the OC as type IV kerogen, formed in a marine environment. The rocks attained thermal maturity possibly during Deccan volcanism. Early Eocene samples of the Naredi Formation (Naliya-Narayan Sarovar Road (NNSR) and the Matanomadh areas) are rich in TOC, smectite, chlorite and framboidal pyrite, but have low \(T_{\mathrm{max}}\). These indicate deposition of sediments in a reducing condition, probably in a lagoonal/marsh/swamp environment. Organic carbon of the Naredi Formation of NNSR may be considered as immature type III to IV kerogen, prone to generate coal. Core samples from the Naredi Formation of the Matanomadh area show two fold distribution in terms of kerogen. Organic carbon of the upper section is immature type III to IV kerogen, but the lower section has type II to III kerogen having potential to generate oil and gas after attaining appropriate thermal maturity.     

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.     

The geochemical and petrological characteristics of prenatal caldera volcano: a case of the newly formed small dacitic caldera,Hijiori, Northeast Japan

Evaluating the magma depth and its physical properties is critical to conduct a better geophysical assessment of magma chambers of caldera volcanoes that may potentially cause future volcanic hazards. To understand pre-eruptive conditions of a magma chamber before its first appearance at the surface, this paper describes the case of Hijiori caldera volcano in northeastern Japan, which emerged approximately 12,000 years ago at a place where no volcano ever existed. We estimated the depth, density, bulk modulus, vesicularity, crystal content, and bulk H\(_2\)O content of the magma chamber using petrographic interpretations, bulk and microchemical compositions, and thermodynamic calculations. The chemical mass balance calculations and thermodynamic modeling of the erupted magmas indicate that the upper portion of the Hijiori magmatic plumbing system was located at depths between 2 and 4 km, and had the following characteristics: (1) pre-eruptive temperature: about 780 \(^{\circ }\)C; (2) bulk magma composition: 66 ± 1.5 wt% SiO\(_{2}\); (3) bulk magmatic H\(_2\)O: approximately 2.5 wt%, and variable characteristics that depend on depth; (4) crystal content: \(\le\)57 vol%; (5) bulk modulus of magma: 0.1–0.8 GPa; (6) magma density: 1.8–2.3 g/cm³; and (7) amount of excess magmatic H\(_2\)O: 11–32 vol% or 48–81 mol%. The range of melt water contents found in quartz-hosted melt inclusions (2–9 wt%) suggests the range of depth phenocrysts growth to be wide (2\(\sim\)13 km). Our data suggest the presence of a vertically elongated magma chamber whose top is nearly solidified but highly vesiculated; this chamber has probably grown and re-mobilized by repeated injections of a small amount of hot dacitic magma originated from the depth.     

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.