Petrogenesis of a voluminous Quaternary adakitic volcano: the case of Baru volcano

The origin of adakite magmas remains controversial because initially the term adakite had petrogenetic significance implying an origin by direct melting of the eclogitized subducting oceanic crust. Many models have been produced for their origin, and until now there has not been a straightforward method to discriminate between these models in a given adakite suite. Here, we use detailed chronological and geochemical studies of selected adakitic edifices that allows for the determination of the magmatic output rate parameter (Qe), which has been correlated with the rates of magma generation deep within subduction zones. By providing temporal and eruption rate estimates, we provide constraints on the possible petrogenetic processes involved in the generation of adakite-like signatures. Adakite magmas derived from the melting of the subducting slab should be volumetrically insignificant when compared to the adakite-like magmas produced by typical arc magma generation processes. In this study, we use this observation and the extraordinary stratigraphic exposure from Miocene to present in an adakitic volcano in Panama and to study the temporal and chemical variation in erupted magmas to estimate rates of magma generation. Detailed chemical and geochronological analyses of Baru volcano indicate that the volcanic edifice was constructed in its entirety during the Quaternary and magmas display adakite-like features such as steep rare earth elements patterns, pronounced depletions in the heavy rare earth elements, low Y, high Sr, and high Sr/Y. The magmatic output rates (Qe) that we have calculated show that compared to other typical adakitic volcanoes, most of the volcanic edifice of Baru volcano was constructed extremely rapidly (<~213 k.a.) and in time frames that are similar to typical arc volcanoes. The observed chemical and mineralogical variation, coupled with the high magma production rates, indicate that Baru volcano is more representative of a typical arc volcano than a small-volume melt of the subducting oceanic crust. The technique we outline may have broader application in determining the petrogenetic conditions of other adakite suites.     

Origin and significance of poikilitic and mosaic peridotite xenoliths in the western Pannonian Basin: geochemical and petrological evidences

Peridotite xenoliths erupted by alkali basaltic volcanoes in the western Pannonian Basin can be divided into two fundamentally contrasting groups. Geochemical characteristics of the abundant protogranular, porphyroclastic and equigranular nodules suggest that these samples originate from an old consolidated and moderately depleted lithospheric mantle domain. In contrast, the geochemical features of the worldwide rare, but in the Pannonian Basin relatively abundant, poikilitic xenoliths attest to a more complex evolution. It has been argued that the origin of the peculiar texture and chemistry may be intimately linked to melt/rock reactions at successively decreasing liquid volumes in a porous melt flow system. The most likely site where such reactions can take place is the asthenosphere–lithosphere boundary. In this context, poikilitic xenoliths may provide petrological and geochemical evidence for reactions between magmatic liquids issued from the uprising asthenosphere and the solid mantle rocks of the lithosphere. These reactions are important agents of the thermal erosion of the lithosphere; thus, they could have considerably contributed to the thinning of the lithosphere in the Pannonian region. We suggest that in the Pannonian Basin, there could be a strong relation between the unusual abundance of poikilitic mantle xenoliths and the strongly eroded lithosphere.     

Overstepping the garnet isograd: a comparison of QuiG barometry and thermodynamic modeling

The consequences of overstepping the garnet isograd reaction have been investigated by comparing the composition of garnet formed at overstepped P–T conditions (the overstep or “OS” model) with the P–T conditions that would be inferred by assuming garnet nucleated in equilibrium with the matrix assemblage at the isograd (the equilibrium or “EQ” model). The garnet nucleus composition formed at overstepped conditions is calculated as the composition that produces the maximum decrease in Gibbs free energy from the equilibrated, garnet-absent, matrix assemblage for the bulk composition under study. Isopleths were then calculated for this garnet nucleus composition assuming equilibrium with the matrix assemblage (the EQ model). Comparison of the actual P–T conditions of nucleation (the OS model) with those inferred from the EQ model reveals considerable discrepancy between the two. In general, the inferred garnet nucleation P–T conditions (the EQ model) are at a lower temperature and higher or lower pressure (depending on the coexisting calcic phase(s)) than the actual (OS model) nucleation conditions. Moreover, the degree of discrepancy increases with the degree of overstepping. Independent estimates of the pressure of nucleation of garnet were made using the Raman shift of quartz inclusions in garnet (quartz-in-garnet or QuiG barometry). To test the validity of this method, an experimental synthesis of garnet containing quartz inclusions was made at 800 °C, 20 kbar, and the measured Raman shift reproduced the synthesis conditions to within 120 bars. Raman band shifts from three natural samples were then used to calculate an isochore along which garnet was presumed to have nucleated. Model calculations were made at several temperatures along this isochore (the OS model), and these P–T conditions were compared to those computed assuming equilibrium nucleation (the EQ model) to estimate the degree of overstepping displayed by these samples. A sample from the garnet isograd in eastern Vermont is consistent with overstepping of around 10 degrees and 0.6 kbar (affinities of around 2 kJ/mole garnet). A sample from the staurolite–kyanite zone in the same terrane requires overstepping of around 50 °C and 2–5 kbar (affinities of around 10–18 kJ/mole garnet). A similar amount of overstepping was inferred for a blueschist sample from Sifnos, Greece. These results indicate that overstepping of garnet nucleation reactions may be common and pronounced in regionally metamorphosed terranes, and that the P–T conditions and paths inferred from garnet zoning studies may be egregiously in error.     

Hydrogeochemical investigation and groundwater quality assessment of Pratapgarh district,Uttar Pradesh

Hydogrochemical investigation of groundwater resources of Paragraph district has been carried out to assess the solute acquisition processes and water quality for domestic and irrigation uses. Fifty-five groundwater samples were collected and analyzed for pH, electrical conductivity, total dissolved solids, hardness, major anions (F^?, Cl^?, NO₃, HCO₃ ^?, SO₄ ^2?) and cations (Ca²⁺, Mg²⁺, Na⁺, K⁺). Study results reveal that groundwater of the area is alkaline in nature and HCO₃ ^?, Cl^?, Mg²⁺, Na⁺ and Ca²⁺ are the major contributing ions to the dissolved solids. The hydrogeochemical data suggest that weathering of rock forming minerals along with secondary contributions from agricultural and anthropogenic sources are mainly controlling the groundwater composition of Pratapgarh district. Alkaline earth metals (Ca²⁺+Mg²⁺) exceed alkalis (Na⁺+K⁺) and weak acid (HCO₃ ^?) dominate over strong acids (Cl^?+SO₄ ^2?) in majority of the groundwater samples. Ca-Mg-HCO₃ and Ca-Mg-Cl-HCO₃ are the dominant hydrogeochemical facies in the groundwater of the area. The computed saturation indices demonstrate oversaturated condition with respect to dolomite and calcite and undersaturated with gypsum and fluorite. A comparison of groundwater quality parameters in relation to specified limits for drinking water shows that concentrations of TDS, F^?, NO₃ ^? and total hardness exceed the desirable limits in many water samples. Quality assessment for irrigation uses reveal that the groundwater is good for irrigation. However, values of salinity, sodium adsorption ratio (SAR), residual sodium carbonate (RSC), %Na and Kelley index are exceeding the prescribed limit at some sites, demanding adequate drainage and water management plan for the area.     

Remote Sensing based assessment of glacial lake growth on Milam glacier,Goriganga basin,Kumaon Himalaya

The climate change of the twentieth century had an evident effect on glacier environments of the Himalaya. Temporal images of Indian Remote Sensing satellites provide an opportunity to monitor the recession of glacier and development of glacial lakes in the Himalayan cryosphere with a cost to time benefit ratio. The recession of Milam glacier and subsequent growth of a proglacial lake near the snout was analysed using Resourcesat-1 and Resourcesat-2 data. The recession of 480 m during 2004 to 2011 and growth of 47 epiglacial ponds over Milam glacier shows the glacier is in a state of imbalance and losing the ice by downwasting.     

Lithospheric structure below transantarctic mountain using receiver function analysis of TAMSEIS data

The lithospheric structure of Antarctica has been investigated from P- (PRF) and S- receiver functions (SRF) using the seismological data from Trans-Antarctic Mountain Seismic Experiment (TAMSEIS). For the stations deployed on the thick ice sheet, estimation of crustal parameters from PRF may be erroneous as the Moho conversions may interfere with the reverberations within the thick ice sheet. However, the free surface multiples are well observed in PRF. On the other hand, in SRFs, the primary conversions of interest and multiples are separated by the mother S-phase. Therefore, it is advantageous to interpret PRF and SRF jointly for the regions where we have thick low velocity layer at the top such as ice or sediments. The crustal structure and corresponding parameters have already been estimated by various workers, but here we interpret the PRF and SRF jointly to minimize the ambiguity and map the lithospheric architecture below TAM. Our analysis reveals that the average crustal thickness beneath the east Antarctica craton is ～44 km with Vp/Vs ranging between ～1.7 and 1.9. Below Trans-Antarctic Mountain (TAM), the average crustal thickness is ～36 km with higher Vp/Vs of ～1.8–2.0. The rift and the volcanic affected coastal region show erratic depths and Vp/Vs, primarily due to the absence of either primary conversion or multiples in the receiver functions. A small number of stations far from the volcano show that the crust is thinnest (～26 to 34 km thick) in the coastal part. The contribution of this study is the mapping of the lithospheric configuration, not done so far using SRF. The SRF section along a profile spanning E-, W- Antarctica and TAM reveals that the lithospheric thickness in the coast is ～80 km and below TAM it is ～120 km. In the central thick ice cover region, the lithosphere thickens upto ～150 km towards Vostok highlands. The most intriguing feature in our SRF section is that the crust and lithosphere are shallow below TAM compared to the E- Antarctica. Further, we observe a mid-lithospheric low velocity layer confined mostly below TAM, suggesting that the thermal buoyancy could be the prime cause for the upliftment of TAM.     

Nine years of GPS measurements of crustal deformation at Imphal,Indo-Burmese wedge

We report results of 9 years of GPS measurements of crustal deformation at Imphal, Manipur, a site located in the Indo-Burmese wedge of northwest Sunda arc. The analysis of these measurements suggests that the site moves at a rate of about 36.3±0.5 mm/year towards N55° in the ITRF2008. With respect to the Indian plate it moves at a rate of 16.7 mm/year towards N222°, i.e., predominantly towards southwest. The site is located about 15 km east of the Churachandpur Mao fault (CMF), which is reported to accommodate part of the India-Sunda motion. The site motion is not significantly affected by the earthquakes that occurred in the nearby region. However, the 2004 Sumatra-Andaman earthquake caused a coseismic displacement of ～ 3–5 mm predominantly towards southwest. The site motion is almost linear, with some seasonal variation, and does not show any evidence of accelerated slip or slow earthquake on the CMF or along the plate boundary.     

Thermal and rheological controls on the formation of mafic enclaves or banded pumice

Magma mixing can occur in a fluid manner to produce banded pumice or in a brittle manner to form enclaves. We propose that the critical control on mixing style is a competition between developing networks of crystals in the intruding magma that impart a strength to the magma and melting and disrupting those networks in the host. X-ray computed tomography analysis demonstrates that banded pumice from the 1915 Mt. Lassen eruption lacks crystal networks. In contrast, rhyodacite hosts with mafic enclaves from Chaos Crags contain well-developed networks of large crystals. We present a one-dimensional conductive cooling model that predicts mixing style, either ductile or brittle, as a function of magma compositions, temperatures, and the size of the intruding dike. Our model relies on three assumptions: (1) Mixing is initiated by the injection of a hot dike into a cooler magma body with a yield strength; (2) when magma crystallinity exceeds a critical value, 13 vol% plagioclase, the magma develops a yield strength; and (3) when total crystallinity exceeds 40 vol%, the magma has a penetrative crystal network and is effectively solid. Importantly, because the two magmas are of different compositions, their crystallinities and viscosities do not have the same variations with temperature. As the intruding magma cools, it crystallizes from the outside in, while simultaneously, host magma temperature near the intruder rises. Mixing of the two magmas begins when the host magma is heated sufficiently to (1) disrupt the crystal network and (2) initiate convection. If the shear stress exerted by the convecting host magma on the dike is greater than the yield strength of the dike margin (and dike crystallinity does not exceed 40 %), then fluid mixing occurs, otherwise enclaves form by brittle deformation of the dike. Application of the model to magma compositions representative of Lassen and Chaos Crags shows that emplacement of dikes <1 m thick should produce enclaves, whereas thicker dikes should generate fluid mixing and form banded pumice within days to weeks of emplacement. Similar relationships apply to other modeled magmatic systems, including Pinatubo, Unzen, and Ksudach/Shtuybel’ volcanoes. For all studied systems, the absolute size of the intruding dike, not just its proportion relative to the host, influences mixing style.     

Primary melt from Sannome-gata volcano,NE Japan arc: constraints on generation conditions of rear-arc magmas

The conditions under which rear-arc magmas are generated were estimated using primary basalts from the Sannome-gata volcano, located in the rear of the NE Japan arc. Scoriae from the volcano occur with abundant crustal and mantle xenoliths, suggesting that the magma ascended rapidly from the upper mantle. The scoriae show significant variations in their whole-rock compositions (7.9–11.1 wt% MgO). High-MgO scoriae (MgO > ~9.5 wt%) have mostly homogeneous ⁸⁷Sr/⁸⁶Sr ratios (0.70318–0.70320), whereas low-MgO scoriae (MgO < ~9 wt%) have higher ⁸⁷Sr/⁸⁶Sr ratios (>0.70327); ratios tend to increase with decreasing MgO content. The high-MgO scoriae are aphyric, containing ~5 vol% olivine microphenocrysts with Mg# [100 × Mg/(Mg + Fe²⁺)] of up to 90. In contrast, the low-MgO scoriae have crustal xenocrysts of plagioclase, alkali feldspar, and quartz, and the mineralogic modes correlate negatively with whole-rock MgO content. On the basis of these observations, it is inferred that the high-MgO scoriae represent primary or near-primary melts, while the low-MgO scoriae underwent considerable interaction with the crust. Using thermodynamic analysis of the observed petrological features of the high-MgO scoriae, the eruption temperature of the magmas was constrained to 1,160–1,220 °C. Given that the source mantle was depleted MORB-source mantle, the primary magma was plausibly generated by ~7 % melting of a garnet-bearing spinel peridotite; taking this into consideration, and considering the constraints of multi-component thermodynamics, we estimated that the primary Sannome-gata magma was generated in the source mantle with 0.5–0.6 wt% H₂O at 1,220–1,230 °C and at ~1.8 GPa, and that the H₂O content of the primary magma was 6–7 wt%. The rear-arc Sannome-gata magma was generated by a lower degree of melting of the mantle at greater depths and lower temperatures than the frontal-arc magma from the Iwate volcano, which was also estimated to be generated by ~15 % melting of the source mantle with 0.6–0.7 wt% H₂O at ~1,250 °C and at ~1.3 GPa.     

Settling and compaction of chromite cumulates employing a centrifuging piston cylinder and application to layered mafic intrusions

The time scales and mechanics of gravitationally driven crystal settling and compaction is investigated through high temperature (1,280–1,500 °C) centrifuge-assisted experiments on a chromite-basalt melt system at 100–1,500g (0.5 GPa). Subsequently, the feasibility of this process for the formation of dense chromite cumulate layers in large layered mafic intrusions (LMIs) is assessed. Centrifugation leads to a single cumulate layer formed at the gravitational bottom of the capsule. The experimentally observed mechanical settling velocity of a suspension of ~24 vol% chromite is calculated to be about half (~0.53) of the Stokes settling velocity, with a sedimentation exponent n of 2.35 (3). Gravitational settling leads to an orthocumulate layer with a porosity of 0.52 (all porosities as fraction). Formation times for such a layer from a magma with initial chromite contents of 0.1–1 vol% are 140–3.5 days, equal to a growth rate of 0.007–0.3 m/day for grain sizes of 1–2 mm. More compacted chromite layers form with increasing centrifugation time and acceleration through chemical compaction: An increase of grain contact areas and grain sizes together with a decrease in porosity is best explained by pressure dissolution at grain contacts, reprecipitation and grain growth into the intergranular space and a concomitant expulsion of intergranular melt. The relation between the porosity in the cumulate pile and effective pressure integrated over time (Δρ · h · a · t) is best fit with a logarithmic function, in fact confirming that a (pressure) dissolution–reprecipitation process is the dominant mechanism of compaction. The experimentally derived equation allows calculating compaction times: 70–80 % chromite at the bottom of a 1-m-thick chromite layer are reached after 9–250 years, whereas equivalent compaction times are 0.2–0.9 years for olivine (both for 2 mm grain size). The experiments allow to determine the bulk viscosities of chromite and olivine cumulates to be of magnitude 10⁹ Pa s, much lower than previously reported. As long as melt escape from the compacting cumulate remains homogeneous, fluidization does not play any role; however, channelized melt flow may lead to suspension and upward movement of cumulate crystals. In LMIs, chromitite layers are typically part of a sequence with layers of mafic minerals, compaction occurs under the additional weight of the overlying layers and can be achieved in a few years to decades.