Patchy distribution of magma that fed the Bishop Tuff supereruption: Evidence from matrix glass major and trace-element compositions

NASA Astrophysics Data System (ADS)

Gualda, G. A. R.; Ghiorso, M. S.; Hurst, A. A.; Allen, M. C.; Bradshaw, R. W.

2017-12-01

For more than 40 years, the Bishop Tuff has been the archetypical example of a singular, zoned magma body that fed a supereruption. Early-erupted material is pyroxene-free and crystal poor (<20 wt. %), presumably erupted from the upper parts of the magma body; late-erupted material is orthopyroxene and clinopyroxene-bearing, commonly more crystal rich (up to 30 wt. % crystals), and presumably tapped magma from the lower portions of the magma body. Fe-Ti oxide compositions suggest higher crystallization temperatures for late-erupted magmas (as high as 820 °C) than for early-erupted magmas (as low as 700 °C). Pressures and temperatures derived from major element compositions of glass inclusions led Gualda & Ghiorso (2013, CMP) to suggest an alternative model of lateral juxtaposition of two main magma bodies - each one feeding early-erupted and late-erupted units. Chamberlain et al. (2015, JPet) and Evans et al. (2016, AmMin) recently disputed this interpretation. We present a large dataset of matrix glass compositions for 161 pumice clasts that span the stratigraphy of the deposit. We calculate crystallization pressures based on major-element glass compositions using rhyolite-MELTS geobarometry, and crystallization temperatures based on Zr in glass using zircon saturation geothermometry. We apply the same methods to 1538 major-element and 615 trace-element analyses from Chamberlain et al. The results overwhelmingly demonstrate that there is no difference in crystallization temperature or pressure between early and late-erupted magmas. Crystallization pressures and temperatures are unimodal, with modes of 150 MPa and 730 °C (calibration of Watson & Harrison). Our results strongly support lateral juxtaposition of two main magma bodies. Smaller units recognized by Chamberlain et al. crystallized at the same pressures as the main bodies - this suggests the coexistence of larger and smaller magma bodies at the time of the Bishop Tuff supereruption. We compare our findings for the Bishop Tuff with results for very large and supereruptions elsewhere in the world. We argue that supereruptions typically mobilize a complex patchwork of magma bodies that reside within specific levels of the crust. They reveal moments of high-melt productivity in the crust, unlike what we observe in the Earth today.

Metamorphism, graphite crystallinity, and sulfide anatexis of the Rampura-Agucha massive sulfide deposit, northwestern India

NASA Astrophysics Data System (ADS)

Mishra, Biswajit; Bernhardt, Heinz-Jurgen

2009-02-01

Located adjacent to the Banded Gneissic Complex, Rampura-Agucha is the only sulfide ore deposit discovered to date within the Precambrian basement gneisses of Rajasthan. The massive Zn-(Pb) sulfide orebody occurs within graphite-biotite-sillimanite schist along with garnet-biotite-sillimanite gneiss, calc-silicate gneisses, amphibolites, and garnet-bearing leucosomes. Plagioclase-hornblende thermometry in amphibolites yielded a peak metamorphic temperature of 720-780°C, whereas temperatures obtained from Fe-Mg exchange between garnet and biotite (580-610°C) in the pelites correspond to postpeak resetting. Thermodynamic considerations of pertinent silicate equilibria, coupled with sphalerite geobarometry, furnished part of a clockwise P- T- t path with peak P- T of ˜6.2 kbar and 780°C, attained during granulite grade metamorphism of the major Zn-rich stratiform sedimentary exhalative deposits orebody and its host rocks. Arsenopyrite composition in the metamorphosed ore yielded a temperature [and log f( S 2)] range of 352°C (-8.2) to 490°C (-4.64), thus indicating its retrograde nature. Contrary to earlier research on the retrogressed nature of graphite, Raman spectroscopic studies on graphite in the metamorphosed ore reveal variable degree of preservation of prograde graphite crystals (490 ± 43°C with a maximum at 593°C). The main orebody is mineralogically simple (sphalerite, pyrite, pyrrhotite, arsenopyrite, galena), deformed and metamorphosed while the Pb-Ag-rich sulfosalt-bearing veins and pods that are irregularly distributed within the hanging wall calc-silicate gneisses show no evidence of deformation and metamorphism. The sulfosalt minerals identified include freibergite, boulangerite, pyrargyrite, stephanite, diaphorite, Mn-jamesonite, Cu-free meneghinite, and semseyite; the last three are reported from Agucha for the first time. Stability relations of Cu-free meneghinite and semseyite in the Pb-Ag-rich ores constrain temperatures at >550°C and <300°C, respectively. Features such as (1) low galena-sphalerite interfacial angles, (2) presence of multiphase sulfide-sulfosalt inclusions, (3) microcracks filled with galena (±pyrargyrite) without any hydrothermal alteration, and (4) high contents of Zn, Ag (and Sb) in galena, indicate partial melting in the PbS-Fe0.96S-ZnS-(1% Ag2S ± CuFeS2) system, which was critical for metamorphic remobilization of the Rampura-Agucha deposit.

A geodynamic constraint on Archean continental geotherms

NASA Astrophysics Data System (ADS)

Bailey, R. C.

2003-04-01

Dewey (1988) observed that gravitational collapse appears to currently limit the altitudes of large plateaus on Earth to about 3 to 5 km above sea level. Arndt (1999) summarized the evidence for the failure of large parts of the continental crust to reach even sea-level during the Archean. If this property of Archean continental elevations was also enforced by gravitational collapse, it permits an estimation of the geothermal gradient in Archean continental crust. If extensional (collapse) tectonics is primarily a balance between gravitational power and the power consumed by extensional (normal) faulting in the upper brittle crust, as analysed by Bailey (1999), then it occurs when continental elevations above ocean bottoms exceed about 0.4 times the thickness of the brittle crust (Bailey, 2000). Assuming an Archean oceanic depth of about 5 km, it follows that that the typical thickness of Archean continental brittle crustal must have been less than about 12 km. Assuming the brittle-ductile transition to occur at about 350 degrees Celsius, this suggests a steep geothermal gradient of at least 30 degrees Celsius per kilometer for Archean continents, during that part of the Archean when continents were primarily submarine. This result does not help resolve the Archean thermal paradox (England and Bickle, 1984) whereby the high global heat flow of the Archean conflicts with the rather shallow crustal Archean geotherms inferred from geobarometry. In fact, the low elevation of Archean continental platforms raises another paradox, a barometric one: that continents were significantly below sea-level implies, by isostasy, that continental crustal thicknesses were significantly less than 30 km, yet the geobarometric data utilized by England and Bickle indicated burial pressures of Archean continental material of up to 10 kb. One resolution of both paradoxes (as discussed by England and Bickle) would be to interpret such deep burials as transient crustal thickening events of duration less than the crustal thermal equilibriation time (about 10 to 30 Ma). Temporary entrainment in the wake of basal eclogite ``sinkers'' might provide such transient burial. Vlaar's (1994) modelling of this eclogite delamination process (tectonically elaborated by Zegers and van Keken (2001)) indicates such sinker events would be significantly shorter than 10 Ma. The topographic re-equilibriation of a hot moho above such a process would be similarly short (Kaufmann and Royden, 1994).

Geochemical and isotopic evidence for the petrogenesis and emplacement tectonics of the Serra dos Órgãos batholith in the Ribeira belt, Rio de Janeiro, Brazil

NASA Astrophysics Data System (ADS)

Machado, Rômulo; Philipp, Ruy Paulo; McReath, Ian; Peucat, Jean Jacques

2016-07-01

The Serra dos Órgãos batholith in the State of Rio de Janeiro (Brazil) is a NE-SW-trending elongated body that occupies ca. 5000 km2 in plan view. It is a foliated intrusion, especially at its borders and is crosscut by syn-magmatic shear zones, with foliations that are moderately-to steeply-dipping to the northwest and moderately-to shallow-dipping in the center and to the southeast, in a configuration of a large laccolith. It was emplaced between 560 and 570 Ma, during an extensional episode that was part of a series of events that comprise the Brasiliano Orogeny in SE Brazil, and which includes deformation, metamorphism and granite intrusion during the interval between 630 and 480 Ma. The two main rock types in the batholith are biotite-hornblende monzogranite, and biotite leucogranite, with subordinate tonalite, granodiorite, diorite, quartz diorite (enclaves), aplite and pegmatite. Harker-type diagrams help show two rock groups with similar trends of evolution: a dioritic and a granitic. The first one is tholeiitic, whereas the second is calc-alkaline, with medium-to high-K calc-alkaline affinity and metaluminous to slightly peraluminous character. In both groups strong decrease in Al2O3, MgO, FeOT and CaO relative to silica contents are observed, which is compatible with trends of fractional crystallization involving clinopyroxene and/or hornblende, plagioclase, opaque minerals, apatite, microcline and biotite. The Sr and Nd isotopic data suggest recycling of a Paleoproterozoic crust as an important petrological process to generate the batholith rocks. Geothermometry (amphibole composition) and geobarometry (saturation in zircon and apatite) indicate that most of the batholith solidified at mid to lower crustal levels at about 750 °C and between 5 and 5.5 kbar. We consider that Serra dos Órgãos crustal protoliths underwent melting caused by the interaction with hotter mafic magma at the base of the crust. These two magmas, with distinct initial compositions and rheology, probably underwent mixing and mingling. This process continued during the rise of the magma through the crust, which was accompanied by magmatic differentiation. The main feature that characterizes the post-collisional Serra dos Órgãos granite magmatism is the connection with high angle ductile shear zones of continental scale and presence to a greater or lesser extent of mafic magmas.

Mantle compositions below petit-spot volcanoes of the NW Pacific Plate

NASA Astrophysics Data System (ADS)

Hirano, N.

2017-12-01

Monogenetic petit-spot volcanoes of a few kilometers in diameter and <300 m in height form volcanic clusters on the subducting NW Pacific plate offshore from NE Japan. Three of these petit-spot provinces form clusters with extents of 1,000-10,000 km2, containing between 15 to 90 monogenetic volcanoes, respectively (Hirano et al., 2008). The magmas that form these volcanoes originate below the lithosphere and ascend along the concavely flexed zone of the outer-rise prior to plate subduction at the trench (Hirano et al., 2006). This forms a unique opportunity to geochemically examine the mantle beneath the oceanic crust in a region outside of the well-examined but spatially restricted areas of mid-oceanic ridges and hotspots, indicating that these petit-spot lavas and associated xenoliths can directly provide the information on the asthenospheric and lithospheric material within and beneath old and subducting plates. Recent research into the geochemistry of petit-spot lavas and the petrography of xenoliths within these lavas indicates that the conventional subducting lithospheric theories require some revision in terms of the nature of subducting lithospheric and asthenospheric materials (e.g., heterogeneous asthenosphere and the presence of a higher geothermal gradient than the conventional GDH1 model; Machida et al., 2015; Yamamoto et al., 2014). The fact that the majority of the petit-spot lava samples do not contain olivine phenocrysts and have differentiated compositions (45-52 wt% SiO2, Mg# values of 50-65) indicates that these magmas have undergone differentiation in a magma chamber. However, geobarometry indicates that the deepest-sourced associated peridotitic xenoliths were derived from a depth of 42 km (Yamamoto et al., 2014). This indicates that melt fractionation must have occurred at depths greater than the middle lithosphere, a situation where the depth of fractionation could correlate with the rotation of the σ3 stress axis from the extensionally lower to the compressional upper part of the lithosphere. This rotation is the result of concave flexure prior to the outer rise of the subduction zone (Valentine & Hirano, 2010). Pilet et al. (2016) and Yamamoto et al. (2009) reported that these xenoliths were derived from a metasomatized region of the mantle, with this region metasomatized by prior melts of petit-spot magmas in the province. The strategic analysis of xenocrystic olivines from several petit-spot volcanoes also indicates that more depleted compositions are located in areas more proximal to the trench. This indicates that the lithospheric mantle in this region must have been significantly metasomatized prior to the onset of trench subduction.

Depth of Formation of Ferropericlase Included in Super-Deep Diamonds

NASA Astrophysics Data System (ADS)

Anzolini, C.; Nestola, F.; Gianese, A.; Nimis, P.; Harris, J. W.

2017-12-01

Super-deep diamonds are believed to have formed at depths of at least 300 km depth (Harte, 2010). A common mineral inclusion in these diamonds is ferropericlase, (Mg,Fe)O (see Kaminsky, 2012 and references therein). Ferropericlase (fPer) is the second most abundant mineral in the lower mantle, comprising approximately 16-20 wt% (660 to 2900 km depth), and inclusions of fPer in diamond are often considered to indicate a lower-mantle origin (Harte et al., 1999). Samples from São Luiz/Juina, Brazil, are noteworthy for containing nanometer-sized magnesioferrite (Harte et al., 1999; Wirth et al., 2014; Kaminsky et al., 2015; Palot et al., 2016). Based upon a phase diagram valid for 1 atm, such exsolutions would place the origin of this assemblage in the uppermost part of the lower mantle. However, a newly reported phase diagram for magnesioferrite demonstrates that the latter is not stable at such pressures and, thus, it cannot exsolve directly from fPer at lower-mantle conditions (Uenver-Thiele et al., 2017). Here we report the investigation of two fPer inclusions, extracted from a single São Luiz diamond, by single-crystal X-ray diffraction and field emission scanning electron microscopy. Both techniques showed micrometer-sized exsolutions of magnesioferrite within the two fPers. We also completed elastic geobarometry (see Angel et al., 2015), which determined an estimate for the depth of entrapment of the two ferropericlase - diamond pairs. In the temperature range between 1273 and 1773 K, pressures varied between 9.88 and 12.34 GPa (325-410 km depth) for one inclusion and between 10.69 and 13.16 GPa (350-440 km depth) for the other one. These results strengthen the hypothesis that solitary fPer inclusions might not be reliable markers for a lower-mantle provenance. This work was supported by Fondazione CaRiPaRo and ERC-2012-StG 307322 to FN. Angel, R.J., et al. (2015) Russ Geol Geophys, 56, 211-220; Harte, B. (2010) Mineral Mag, 74, 189-215; Harte, B., et al. (1999) The Geochemical Society, Special Publication, 6, 125-153; Kaminsky, F. (2012) Earth-Sci Rev, 110, 127-147; Kaminsky, F., et al. (2015) Earth Planet Sci Lett, 417, 49-56; Palot, M., et al. (2016) Lithos, 265, 237-243; Uenver-Thiele, L., et al. (2017) Am Mineral, 102, 632-642; Wirth, R., et al. (2014) Earth Planet Sci Lett, 404, 365-375.

Crustal structure in the Elko-Carlin Region, Nevada, during Eocene gold mineralization: Ruby-East Humboldt metamorphic core complex as a guide to the deep crust

USGS Publications Warehouse

Howard, K.A.

2003-01-01

The deep crustal rocks exposed in the Ruby-East Humboldt metamorphic core complex, northeastern Nevada, provide a guide for reconstructing Eocene crustal structure ~50 km to the west near the Carlin trend of gold deposits. The deep crustal rocks, in the footwall of a west-dipping normal-sense shear system, may have underlain the Pinon and Adobe Ranges about 50 km to the west before Tertiary extension, close to or under part of the Carlin trend. Eocene lakes formed on the hanging wall of the fault system during an early phase of extension and may have been linked to a fluid reservoir for hydrothermal circulation. The magnitude and timing of Paleogene extension remain indistinct, but dikes and tilt axes in the upper crust indicate that spreading was east-west to northwest-southeast, perpendicular to a Paleozoic and Mesozoic orogen that the spreading overprinted. High geothermal gradients associated with Eocene or older crustal thinning may have contributed to hydrothermal circulation in the upper crust. Late Eocene eruptions, upper crustal dike intrusion, and gold mineralization approximately coincided temporally with deep intrusion of Eocene sills of granite and quartz diorite and shallower intrusion of the Harrison Pass pluton into the core-complex rocks. Stacked Mesozoic nappes of metamorphosed Paleozoic and Precambrian rocks in the core complex lay at least 13 to 20 km deep in Eocene time, on the basis of geobarometry studies. In the northern part of the complex, the presently exposed rocks had been even deeper in the late Mesozoic, to >30 km depths, before losing part of their cover by Eocene time. Nappes in the core plunge northward beneath the originally thicker Mesozoic tectonic cover in the north part of the core complex. Mesozoic nappes and tectonic wedging likely occupied the thickened midlevel crustal section between the deep crustal core-complex intrusions and nappes and the overlying upper crust. These structures, as well as the subsequent large-displacement Cenozoic extensional faulting and flow in the deep crust, would be expected to blur the expression of any regional structural roots that could correlate with mineral belts. Structural mismatch of the mineralized upper crust and the tectonically complex middle crust suggests that the Carlin trend relates not to subjacent deeply penetrating rooted structures but to favorable upper crustal host rocks aligned within a relatively coherent regional block of upper crust.

High-pressure behaviour of Cr-Fe-Mg-Al spinels: applications to diamond geobarometry

NASA Astrophysics Data System (ADS)

Periotto, Benedetta; Bruschini, Enrico; Nestola, Fabrizio; Lenaz, Davide; Princivalle, Francesco; Andreozzi, Giovanni B.; Bosi, Ferdinando

2014-05-01

Spinels belonging to the chromite - magnesiochromite - hercynite (FeCr2O4-MgCr2O4-FeAl2O4) system are among the most common inclusions found in diamonds (Stachel and Harris 2008). In particular, although FeCr2O4 and MgCr2O4 components sum to between 85 and 88% of spinels found in diamonds, hercynite FeAl2O4 plays a not negligible role in determining their thermo-elastic properties with concentrations reaching 7-9 % (other minor end-members like MgAl2O4, MgFe2O4 and Fe2O3 rarely reach 2-3% in total, see Lenaz et al. 2009). Recent studies were focused on the determination of the diamond formation pressure by the so-called "elastic method" (see for example Nestola et al. 2011 and references therein). It was demonstrated that accurate and precise thermo-elastic parameters are fundamental to minimize the uncertainty of formation pressure. In this work we have determined the equations of state at room temperature of three synthetic spinel end-members chromite - magnesiochromite - hercynite and one natural spinel crystal extracted from a diamond (from Udachnaya mine, Siberia, Russia) by single-crystal X-ray diffraction in situ at high-pressure. A diamond-anvil cell was mounted on a STADI IV diffractometer equipped with a point detector and motorized by SINGLE software (Angel and Finger 2011). The natural crystal was investigated to test (and possibly validate) the "empirical prediction model", capable to provide bulk modulus and its first pressure derivative only knowing the composition of the spinels found in diamonds. Such prediction model could be used to obtain pressure of formation for the diamond-spinel pair through the elastic method. Details and results will be discussed. The research was funded by the ERC Starting Grant to FN (grant agreement n° 307322). References Angel R.J., Finger L.W. (2011) SINGLE A program to control single-crystal diffractometers. Journal of Applied Crystallography, 44, 247-251. Lenaz D., Logvinova A.M., Princivalle F., Sobolev N. (2009) Structural parameters of chromite included in diamond and kimberlites from Siberia: a new tool for discriminating source. American Mineralogist, 94, 1067-1070. Nestola F., Nimis P., Ziberna L., Longo M., Marzoli A., Harris J.W., Manghnani M.H., Fedortchouk Y. (2011) First crystal-structure determination of olivine in diamond: composition and implications for provenance in the Earth's mantle. Earth and Planetary Science Letters, 305, 249-255. Stachel, T., and Harris, J.W. (2008) The origin of cratonic diamonds - constraints from mineral inclusions. Ore Geology Reviews, 34, 5-32.

Single inclusion piezobarometry confirms high-temperature decompression path for Variscan granulites

NASA Astrophysics Data System (ADS)

Angel, Ross; Alvaro, Matteo; Mazzucchelli, Mattia; Nimis, Paolo; Nestola, Fabrizio

2016-04-01

The identification and chemistry of inclusions trapped in host minerals during growth of the host phase have long been used to infer P-T points on metamorphic paths. The determination of the remnant pressure on the inclusion, e.g., using data from X-ray diffractometry, birefringence analysis or Raman spectroscopy, provides an alternative method of barometry using elasticity theory. A remnant pressure in an inclusion is developed because the inclusion and the host have different thermal expansion and compressibilities, and the inclusion does not expand in response to P and T as would a free crystal. Instead it is restricted to expand only as much as the host mineral, and this constriction in volume can result in inclusions exhibiting over-pressures when the host is studied at room conditions. This concept has been known for a long time, but satisfactory quantitative modelling of inclusion-host systems based on non-linear elasticity theory and precise thermal-pressure euqations of state has only recently come available (Angel et al., 2014, 2015), even though it is still restricted to elastically isotropic minerals. No mineral is elastically isotropic, but garnets and diamond are almost so. Calculations show that diamonds trapped as inclusions in host silicates at P and T within the stability field of diamond should exhibit zero pressure when the samples are recovered to room conditions. However, some diamond inclusions in garnets in granulites are reported to exhibit significant residual overpressures (e.g., Kotková et al., 2011). This indicates that the inclusion was elastically re-equilibrated (e.g., by plastic flow in the garnet host) at high temperatures and lower pressures in the stability field of graphite, consistent also with the observed partial inversion of diamond to graphite. In this case, the elastic analysis of the diamond-in-garnet inclusions provides qualitative independent evidence that the Variscan granulites underwent pressure reduction at high temperatures. The extension of single inclusion piezobarometry to elastically anisotropic minerals will allow quantitative analysis of diamonds trapped in other minerals such as kyanite. This work was supported by ERC starting grant 307322 to Fabrizio Nestola and by the MIUR-SIR grant "MILE DEEp" (RBSI140351) to M. Alvaro. Angel R.J., Mazzucchelli M.L., Alvaro M., Nimis P. & Nestola F. (2014) Geobarometry from host-inclusion systems: the role of elastic relaxation. Am. Mineral., 99, 2146-2149. Angel R.J., Nimis P., Mazzucchelli M.L., Alvaro M. & Nestola F. (2015) How large are departures from lithostatic pressure? Constraints from host-inclusion elasticity. J. Metamorphic Geol., 33, 801-813. Kotková J, O'Brien P.J & Ziemann M.A. (2011)Diamond and coesite discovered in Saxony-type granulite: Solution to the Variscan garnet peridotite enigma. Geology, 39, 667-670.

Carbon-saturated monosulfide melting in the shallow mantle: solubility and effect on solidus

NASA Astrophysics Data System (ADS)

Zhang, Zhou; Lentsch, Nathan; Hirschmann, Marc M.

2015-12-01

We present high-pressure experiments from 0.8 to 7.95 GPa to determine the effect of carbon on the solidus of mantle monosulfide. The graphite-saturated solidus of monosulfide (Fe0.69Ni0.23Cu0.01S1.00) is described by a Simon and Glatzel (Z Anorg Allg Chem 178:309-316, 1929) equation T (°C) = 969.0[ P (GPa)/5.92 + 1]0.39 (1 ≤ P ≤ 8) and is 80 ± 25 °C below the melting temperature found for carbon-free conditions. A series of comparison experiments using different capsule configurations and preparations document that the observed solidus-lowering is owing to graphite saturation and not an artifact of different capsules or hydrogen contamination. Concentrations of carbon in quenched graphite-saturated monosulfide melt measured by electron microprobe are 0.1-0.3 wt% in monosulfide melt and below the detection limit (<0.2 wt%) in crystalline monosulfide solid solution. Although there is only a small amount of carbon dissolved in monosulfide melts, the substantial effect on monosulfide solidus temperature means that the carbon-saturated monosulfide (Fe0.69Ni0.23Cu0.01S1.00) solidus intersects continental mantle geotherms inferred from diamond inclusion geobarometry at 6-7 GPa ( 200 km), whereas carbon-free monosulfide (Fe0.69Ni0.23Cu0.01S1.00) solidus does not. The composition investigated (Fe0.69Ni0.23Cu0.01S1.00) has a comparatively low metal/sulfur (M/S) ratio and low Ni/(Fe + Ni), but sulfides with higher (M/S) and with greater Ni/(Fe + Ni) should melt at lower temperatures and these should have a broader melt stability field in the diamond formation environment and in the continental lithosphere. Low carbon solubility in monosulfide melt excludes the possibility that diamonds are crystallized from sulfide melt. Although monosulfide melt can store no more than 2 ppm C in a bulk mantle with 225 ppm S, melts with higher M/S could be a primary host of carbon in the deeper part of the upper mantle. For example, the storage capacity of C in sulfide melts in the deep upper mantle ( 400 km) for a depleted mantle domain (MORB source, 120 ± 30 ppm S) is estimated to be 57 ±_{30}^{63} ppm, and so all the C could be in a sulfide melt. In an enriched (OIB source, 225 ± 25 ppm S) mantle domain, the C stored in sulfide melt in the deep upper mantle is estimated to be 86 ±_{44}^{92} ppm, which would amount to about half the available carbon.

Crystallisation condition of the Quaternary basanites of volcanic centre Black Rock, monogenetic field Lunar Crater

NASA Astrophysics Data System (ADS)

Turova, Mariia; Plechov, Pavel; Scherbakov, Vasily; Larin, Nikolay

2017-04-01

The Lunar Crater volcanic field is located in a tension zone Basin and Range Province (USA). This tension is connected with dives oceanic plate under the continental plate [1]. Lunar Crater consists of flows basalt, basanite, trachybasalt has a different age [2]. In this work we investigate the youngest rock - basanite. The basanite is highly crystalline consisting of about megacrysts (3-10 cm) 30-60 wt% phenocrysts ( 800-1500 µm) and microphenocrysts (100-800 µm) and 40-60% microlites (<100 µm). This type of crystal allocated on the basis of size and different chemical composition. The basanite contains about 40 wt % of olivine phenocrysts and microphenocrysts; 35 % clinopyroxene phenocrysts and microphenocrysts. The other phenocrysts and microphenocrysts are feldspar and spinel. Phenocrysts of olivine plagioclase and clinopyroxene are the features of dissolution. The groundmass (<100 µm) consist of microlites olivine, clinopyroxene sanidine Ti-magnetite. Megacrysts are crystals range from I to l0 cm, are free of inclusions, and are unzoned. Basanite also bearing homeogenic enclaves and amphibole-feldspar-clinopyroxene cumulates. This size 4 mm-1.5 cm. Also in some cumulat indentified mineral renit. We determined pressure of the formation of clinopyroxene assemblage using the clinopyroxene barometer based on the relationship between pressure and the volumes of the unit cell and polyhedron M1 in the mineral structure [3]. The pressure is 18-20 kbar for megacrysts, for phenocrysts 15-18 kbar, for microphenocryst 6-8 kbar, for microlites 1,5-3 kbar. Moreover megacrysts are depleted of REE, compared with phenocrysts. Possibly, megacrysts are formed from the same basanite magma during earlier stage of crystallization [4]. Oxygen barometer data shows that the grains were formed in Redox conditions about FMQ+0.2. Temperature and oxygen fugacity conditions were estimated for microphenocrysts and groundmass crystallization only. Bibliography 1. Zoback M. L., Anderson R. E., Thompson G. A. Cainozoic evolution of the state of stress and style of tectonism of the Basin and Range province of the western United States //Philosophical Transactions of the Royal Society of London A: Mathematical, Physical and Engineering Sciences. - 1981. - T. 300. - №. 1454. - C. 407-434. 2. Wood, X., and Keinle, Y., 1990, Volcanoes of North America: Cambridge,United Kingdom, Cambridge University Press, 354 p. 3. Nimis P. Clinopyroxene geobarometry of magmatic rocks. Part 2. Structural geobarometers for basic to acid, tholeiitic and mildly alkaline magmatic systems //Contributions to Mineralogy and Petrology. - 1999. - T. 135. - №. 1. - C. 62-74. 4. Ballhaus C., Berry R. F., Green D. H. High pressure experimental calibration of the olivine-orthopyroxene-spinel oxygen geobarometer: implications for the oxidation state of the upper mantle //Contributions to Mineralogy and Petrology. - 1991. - T. 107. - №. 1. - C. 27-40.

EosFit-Pinc: a GUI program to calculate pressures in host-inclusion systems

NASA Astrophysics Data System (ADS)

Angel, Ross; Alvaro, Matteo; Mazzucchelli, Mattia; Nestola, Fabrizio

2017-04-01

A remnant pressure in an inclusion trapped inside a host mineral is developed because the inclusion and the host have different thermal expansion and compressibilities, and the inclusion does not expand in response to P and T as would a free crystal. Instead it is restricted to expand only as much as the cavity of the host mineral, and this constriction in volume can result in inclusions exhibiting over-pressures when the host is studied at room conditions. The remnant pressure of the inclusion, measured by X-ray diffractometry, birefringence analysis or Raman spectroscopy, can then be used with the equations of state (EoS) of the host and inclusion to constrain the P and T at entrapment. This concept has been known for a long time, but satisfactory quantitative modelling of inclusion-host systems based on non-linear elasticity theory and precise EoS has only recently come available (Angel et al., 2014, 2015), even though calculations still assume isotropic elastic properties. The elasticity calculations to determine entrapment conditions involving the EoSs for both the host and the inclusion are complex if thermodynamically-realistic EoS are employed. We have therefore developed a simple GUI program, EosFit-Pinc that performs all of the necessary calculations under the assumptions of isotropic elasticity. Equations of state of the host and the inclusion can be loaded as files created by other software in the EosFit7 program suite, or imported directly from thermodynamic databases such as Thermocalc. The complete range of EoS types supported by EosFit-7 are available in EosFit-Pinc. Fluid EoS can be provided in the form of PVT tables, which allows fluid inclusions to be modelled. Once loaded, the EoS of the host and inclusion can be used to calculate the entrapment isomeke from the measured remnant pressure of the inclusion. Or the final pressure can be calculated if the entrapment conditions are known or estimated. Calculations of the isochors of both the host and inclusion phases, and their mutual isomekes, can be performed, and output is provided in a format suitable for external plotting programs. The program EosFit-Pinc and the EosFit7 program suite are available at www.rossangel.net This work was supported by ERC starting grant "INDIMEDEA" (307322) to F. Nestola and by the MIUR-SIR grant "MILE DEEp" (RBSI140351) to M. Alvaro. Angel R.J., Mazzucchelli M.L., Alvaro M., Nimis P. & Nestola F. (2014) Geobarometry from host-inclusion systems: the role of elastic relaxation. Am. Mineral., 99, 2146-2149. Angel R.J., Nimis P., Mazzucchelli M.L., Alvaro M. & Nestola F. (2015) How large are departures from lithostatic pressure? Constraints from host-inclusion elasticity. J. Metamorphic Geol., 33, 801-813.

Clinopyroxene geobarometry of magmatic rocks. Part 2. Structural geobarometers for basic to acid, tholeiitic and mildly alkaline magmatic systems

NASA Astrophysics Data System (ADS)

Nimis, Paolo

The crystal structures of 212 experimentally synthesized, igneous clinopyroxenes were modeled from electronprobe chemical data. The coexisting melts span a wide range of petrologically relevant, dry and hydrous compositions, characterized by variable enrichment in silica and alkalis. Experimental conditions pertain to Earth's crust and uppermost mantle (P=0-24kbar garnet absent) and a variety of fO2 values (from CCO-buffered to air-buffered) and mineral assemblages (Cpx+/-Opx+/-Pig+/-Ol+/-Plag+/-Spl +/-Mt+/-Amp+/-Ilm). Unit-cell volume (Vcell) versus M1-polyhedron volume (VM1) relations were investigated over a range of pressures and temperatures using data derived from structure modeling and corrected for thermal expansivity and compressibility. The relationships between pressure and clinopyroxene structural parameters were found to be dependent on the nature of the coexisting melt. To reduce compositional effects, only clinopyroxenes belonging to mildly alkaline (MA) and tholeiitic (TH) series were considered. Pressure was modeled as a linear function of Vcell, VM1, and Mg/(Mg+Fe2+)Cpx ratio. A calibration based on the whole data set (MA+TH) reproduced the experimental pressures within 1.4kbar at the 1-σ level. The maximum residuals were 3.5kbar and 3.9kbar for MA- and TH-clinopyroxenes, respectively. Better statistics were obtained by considering MA- and TH-clinopyroxenes separately. A calibration based on the 69 MA-clinopyroxenes reproduced the experimental pressures within 1.1kbar (1σ) and with a maximum residual of 2.7kbar. A calibration based on the 143 TH-clinopyroxenes reproduced the experimental pressures within 1.0kbar (1σ) and with a maximum residual of 3.4kbar. When these geobarometers are applied to natural samples for which P is unknown, the correction for compressibility is necessarily made through a trial-and-error procedure. This expedient propagates an additional error that increases the above uncertainties and residuals by a factor of about 2. Applications to natural, igneous rocks for which the pressures of crystallization could be constrained based on experimental, petrological or geological evidence yielded pressure estimates that reproduced the expected values to within ca. 2kbar. Compared to the MA-formulation, the TH-formulation appears to be less robust to variations in magma composition. When applied to high-pressure (>10kbar) clinopyroxenes synthesized from very low Na (Na2O<1.5%) melts, the latter geobarometer can underestimate P by as much as 6kbar. Calculation of P through the present geobarometers requires clinopyroxene major-element composition and an independent, accurate estimate of crystallization T. Underestimating T by 20°C propagates into a 1-kbar increase in calculated P. The proposed geobarometers are incorporated in the CpxBar software program, which is designed to retrieve the pressure of crystallization from a clinopyroxene chemical analysis.

Thermobarometry for spinel lherzolite xenoliths in alkali basalts

NASA Astrophysics Data System (ADS)

Ozawa, Kazuhito; Youbi, Nasrrddine; Boumehdi, Moulay Ahmed; Nagahara, Hiroko

2016-04-01

Application of geothermobarometers to peridotite xenoliths has been providing very useful information on thermal and chemical structure of lithospheric or asthenospheric mantle at the time of almost instantaneous sampling by the host magmas, based on which various thermal (e.g., McKenzie et al., 2005), chemical (e.g., Griffin et al., 2003), and rheological (e.g., Ave Lallemant et al., 1980) models of lithosphere have been constructed. Geothermobarometry for garnet or plagioclase-bearing lithologies provide accurate pressure estimation, but this is not the case for the spinel peridotites, which are frequently sampled from Phanerozoic provinces in various tectonic environments (Nixon and Davies, 1987). There are several geobarometers proposed for spinel lherzolite, such as single pyroxene geothermobarometer (Mercier, 1980) and geothermobarometer based on Ca exchange between olivine and clinopyroxene (Köhler and Brey, 1990), but they have essential problems and it is usually believed that appropriated barometers do not exist for spinel lherzolites (O'Reilly et al., 1997; Medaris et al., 1999). It is thus imperative to develop reliable barometry for spinel peridotite xenoliths. We have developed barometry for spinel peridotite xenoliths by exploiting small differences in pressure dependence in relevant reactions, whose calibration was made through careful evaluation of volume changes of the reactions. This is augmented with higher levels of care in application of barometer by choosing mineral domains and their chemical components that are in equilibrium as close as possible. This is necessary because such barometry is very sensitive to changes in chemical composition induced by transient state of the system possibly owing to pressure and temperature changes as well as chemical modification, forming chemical heterogeneity or zoning frequently reported from various mantle xenoliths (Smith, 1999). Thus very carful treatment of heterogeneity, which might be trivial for geothermobarometry based on reactions with large and distinct volume changes, is necessary. Specification of mineral domains and their components representing the thermal state of the mantle just before xenolith extraction is one of the major tasks for the establishment of reliable geothermobarometry for spinel lherzolite xenoliths. Systematic variations of such mineralogical information among xenoliths transported by a single volcanic eruption guarantees proper estimation of a mantle geotherm. For the development of such geobarometry, it is important to choose appropriate xenolith locality, where previous studies provide enough information and where many xenolith samples are available for extending a range of derivation depth. Spinel lherzolite xenoliths in alkali basalts from Bou Ibalhatene maars in the Middle Atlas in Morocco are suitable study target. Geochemical, geochronological, petrological, and rheological aspects of the spinel lherzolite xenoliths have been studied (Raffone et al. 2009; El Messbahi et al., 2015; Witting et al., 2010; El Azzouzi et al., 2010), which show that they represent fragments of the lithospheric mantle formed and modified since 1.7Ga before their extraction from Miocene to recent. We have pinpointed portions of minerals in the xenolith samples and their components representing condition just before their entrapment in magmas, on which appropriate geothermobarometers are applied and detected ~0.5GPa pressure difference (1.5-2.0GPa) for ~100°C variation in temperatures (950-1050°C).

Kulanaokuaiki 3: Product of an Energetic, Diatreme-Like Eruption at Kilauea

NASA Astrophysics Data System (ADS)

Fiske, R. S.; Rose, T. R.; Swanson, D. A.

2006-12-01

Kulanaokuaiki 3 (K-3), one of five units of the Kulanaokuaiki tephra, was erupted at ~AD 850 and blanketed large near-summit areas. Most complete remnants today are found in the Koa`e fault system and on the volcano`s south flank, S and SE of the summit. There, K-3 consists mostly of crystal-rich scoria lapilli contained in two sub-units, generally 1-8 cm thick, separated by a <1 cm "parting" of coarse ash and/or reticulite lapilli. Fine ash (<0.5 mm) makes up <3% of the two scoria units, increasing upward to ~10%. Dense lithic clasts are contained in both sub-units; ~85% of these consist of a wide variety of basalt (some enclosed in cored bombs), and ~12% are fine-coarse gabbro (some containing interstitial glass w/vesicles). The lithics are typically fresh, suggesting that the eruptive conduit pierced pristine parts of the volcano`s edifice rather than long-established, hydrothermally altered conduit systems. Erosion has stripped most K-3 from the south flank, leaving its lithics as scattered lags. Dense clasts, >4 kg and 18 cm across, are found as far as 7 km from the summit; progressively smaller clasts (~3-4 cm) fell at the coastline, 17 km away. The K-3 scoria deposits are unremarkable to the eye, but this belies cryptic vertical zonation that characterizes these units at widespread south-flank localities. The specific gravity of scoria lapilli (7-10 mm dia.) decreases upward in the lower sub-unit, accompanied by decreasing whole-rock MgO values. The pattern is reversed in the upper sub-unit, where specific gravity and MgO values increase upward. Available information suggests the specific gravity and MgO variations correlate with percentages of phenocrystic olivine. Preliminary geobarometry of pyroxene-glass pairs suggests that some gabbro was crystallizing at 5-7 km depth before exploding from the volcano-- far deeper than expected in a phreatomagmatic eruption. We interpret that CO2, known to be released in huge volumes from Kilauea`s summit, and which initially exsolves from basaltic magma at ~10 km depth, was the likely propellant for the diatreme-like K-3 eruption. While reaming a conduit to the surface, the streaming CO2, knicked the upper part of a magma body (likely dike-shaped), initiating its disintegration. The first pulse of the eruption released scoria that, along with spalled conduit wall rocks, erupted to form the lower K-3 sub-unit. Following a brief pause, when the air partly cleared to form the mid-K-3 parting, a second pulse entrained scoria originating from progressively deeper and more olivine-rich parts of the magma body. As a result, scoria containing greater percentages of phenocrystic olivine was erupted, and these were showered over the south flank to produce the observed upside-down "magma-chamber grading" in the upper K-3 sub-unit. Multi-mach exit velocities are visualized, and entrained lithic clasts may have been carried to heights of 15-20 km. These clasts were carried to the southeast as they fell through high-level northwesterly winds.

Petrologic evolution of the Caetano magmatic system: What can we learn from a dissected, 34 Ma caldera in the northern Great Basin, western U.S.A.?

NASA Astrophysics Data System (ADS)

Watts, K. E.; Colgan, J. P.; John, D. A.; Henry, C. D.

2012-12-01

Eruption of the >1,100 km3 Caetano Tuff and formation of the Caetano caldera occurred during the mid-Tertiary ignimbrite flare-up in the Great Basin. Post-collapse extension and faulting created a series of tilted fault blocks that expose >4 km thick intracaldera tuff, two generations of resurgent granitic plutons, silicic ring-fracture intrusions, a tuff dike that fed the early eruption, and pre- and post-caldera andesites. We integrate new petrologic data for extrusive and intrusive Caetano units with geologic mapping and geochronology to provide an exceptional view into the inner workings of a large caldera center. The Caetano Tuff is a phenocryst-rich (~30-50%) ignimbrite with a mineralogy of plagioclase + sanidine + quartz + biotite + orthopyroxene + Fe-Ti oxides ± hornblende + accessory zircon and allanite. Plagioclase crystals in the Caetano Tuff and cogenetic intrusive units span a wide compositional range (>30 mol% An) and have diverse petrographic textures ranging from euhedral phenocrysts to anhedral, sieved crystals with melt-rich cores. Plagioclase compositions measured by electron microprobe for whole rock thin sections are consistent with compositional zoning of the intracaldera tuff shown by XRF whole rock analyses, oligoclase (~10-30 mol% An) and andesine (~30-50 mol% An) in the most evolved (75-77% SiO2) and least evolved (72-74% SiO2) tuff units, respectively. However, orthopyroxene compositions are apparently decoupled from the host tuff composition, with the highest Mg#s (~60-70%) occurring in the most evolved tuff samples. In the Caetano Tuff, equilibrium pairs of Fe-Ti oxides yield an average eruption temperature of 745°C, which is consistent with the average Ti-in-zircon temperature of 750±70°C (1 stdev, n=90 spots) obtained from Ti concentrations measured by SHRIMP for single zircons. Application of Al-in-hornblende geobarometry indicates an average equilibration pressure of 4.5±0.1 kbar, corresponding to mid-crustal magma storage depths of ~14-15 km. In light of our new petrologic data, we highlight the following key points: (1) Diverse crystal cargoes, disequilibrium textures, and wide compositional oscillations in single phenocrysts and among discrete mineral populations indicate prolonged and complex episodes of magma assembly and growth. Based on zircon U-Pb SHRIMP ages that range from ~34-37 Ma, assembly and growth may have spanned ~2-3 Ma, or a 34 Ma Caetano magma chamber may have assimilated older igneous rocks in and around the caldera. (2) Mineral chemistry, U-Pb and Ar-Ar geochronology, O isotope geochemistry, and whole rock major and trace element geochemistry indicate a genetic connection between the Caetano Tuff and resurgent granitic plutons, supporting the role of linked volcanic-plutonic components in caldera settings. (3) Generation and eruption of crystal-rich "monotonous" rhyolite calls into question the prevailing paradigms of crystal-poor rhyolites derived from crystal mushes, or crystal-rich "monotonous intermediates" derived from homogeneous dacitic magma reservoirs. The Caetano Tuff may be a representative end member of caldera-forming eruptions that is important for understanding large-volume rhyolite genesis in the shallow-middle crust.

Archaean and Palaeoproterozoic metamorphic events in the Orekhov-Pavlograd compressional zone, Ukrainian Shield

NASA Astrophysics Data System (ADS)

Yurchenko, A. V.

2012-04-01

The Orekhov-Pavlograd zone (OPZ) is located between the Mesoarchaean-Neoarchaean Middle Dnieper Province and the Mesoarchaean-Palaeoproterozoic Azov Province in the eastern Ukrainian Shield. The OPZ consists of Archaean and Palaeoproterozoic high-grade metamorphic rocks. According U-Pb isotope analyses Archaean methaigneous rocks have age of 3.5-3.3 Ga, and latest AR events dated form both individual grains and metamorphic rims in the tonalite and the granitic vein occurred at about 2.88 Ga ego. Paleoproterozoic zircons from a hornblende granulite have a concordia age of 2.08 Ga [1]. P-T conditions of the 3.5-3.3 Ga processes calculated from the Ti content in zircon are of 730-760°C. Metamorphic event dated as 2.88 Ga is more preserved and detected in some amphibolites after mafic dykes. According to different methods of hornblende-plagioclase geothermometry along with Al- and Ti-geobarometry of hornblende, the amphibolites have formed at temperature of 735-749 °C and pressure of 5.2 to 7.8 kbar. P-T conditions of Paleoproterozoic metamorphic processes have been calculated for a Paleoproterozoic high-Al paragneiss and mafic rocks. On the base of the computer software THERIAK-DOMINO [2], near-isothermal decompression from ca. 8.5 to 6.0 kbar at 650 °C and then to 5.8 kbar at 740 °C has been determined for small irregular garnet grains (grs 4-7% and XMg 0.36-0.37) associated with the same biotite and plagioclase. P-T conditions obtained by means of the P-T pseudosection calculation are identical within errors to those defined by the Grt + Bt + Pl + Ozt geothermometer by [3] and the geobarometer by [4], T = 675 °C and P = 5.6 kbar. Temperature and pressure calculated for assemblage Grt-Pl-Opx-Amph-Ilm-Ru (mafic rock) by using the TWEEQU method shows: 1) high values of pressure and temperature (ca. 7 kbar and 800 °C) are linked with the first metamorphic event with Opx-Cpx assemblage, 2) moderate values (ca. 5 kbar and ca. 600 °C) are referred to the second metamorphic event when amphibole was crystallized instead of orthopyroxene. The latest metamorphic reworking took place at P = 3.3-4 kbar and T = ca. 600 °C. The resulting Paleoproterozoic P-T-t path suggests a clockwise P-T evolution of the OPZ area. Preferences: 1. S.B.Lobach-Zhuchenko, Yu.S.Egorova, A.V.Yurchenko, V.V.Balagansky, G.V.Artemenko, V.P.Chekulaev, N.A.Arestova, 2009. Mineralogical Journal (Kiev), 31 (1): 3-10. 2. de Capitani, C., 2005. THERIAK - DOMINO User's Guide, Version 1402005, http://titan. minpet.unibas.ch/minpet/theriak/theruser.html (04/2007). 3. Ferry, J.M., Spear, F.S., 1978. Experimental calibration of the partitioning of Fe and Mg between biotite and garnet. Contrib. Mineral. Petrol. 66, 113-117. 4. Koziol, A.M., Newton, R.C., 1989. Grossular activity-composition relationship in ternary garnets determined by reversed displaced-equilibrium experiments. Contrib. Mineral. Petrol. 103, 423-433.

A new restoration of the NFP20-East cross section and possible tectonic overpressure in the Penninic Adula Nappe (Central Alps)

NASA Astrophysics Data System (ADS)

Pleuger, J.; Podladchikov, Y.

2012-04-01

The Adula Nappe in the eastern Central Alps is one of the four units in the Alps from which ultrahigh-pressure rocks have been reported. Several very different models for its tectonic history have been published but none of these models is fully satisfactory. In the models of Schmid et al. (1996) and Engi et al. (2001), the main mechanism of exhumation is assumed to be extrusion. The extrusion models require top-to-the-hinterland, i.e. top-to-the-south faulting in the hanging wall of the exhuming nappe for which there is no evidence. Froitzheim et al. (2003) proposed a scenario with two different subduction zones, an internal one in which the South Penninic and Briançonnais domains were subducted, and an external one in which the North Penninc domain and the European margin, including the Adula nappe, were subducted. In this model, the exhumation of the Adula nappe results from the subduction of the overlying sub-Briançonnais and sub-South-Penninic mantle in the internal subduction zone. The Adula nappe would then have been exhumed from below into a top-to-the-north shear zone also affecting the overriding Briançonnais units. The main shortcoming of this model is that otherwise there is little evidence for two Alpine subduction zones. All the models cited above are based on the conversion of peak pressures obtained from geobarometry to depth by assuming lithostatic pressures. This results in a much greater burial depth of the Adula Nappe with respect to the surrounding units which poses major problems when trying to reconcile maximum burial depths of the Penninic nappes with their structural record. We performed a new restoration of the NFP20-East cross section (Schmid et al. 1996) without applying a lithostatic pressure-to-depth conversion but a purely geometrical restoration of deformation events in the Penninic nappe stack. The major constraints on these reconstructions are given by strain estimates for the major deformation phases in the units overlying the Adula Nappe (Mayerat Demarne 1994) and zircon fission track ages (Flisch 1986) indicating that the Austroalpine units have not been more than 10 km below surface after the Palaeocene. The maximum pressures of eclogites from the Adula nappe reported in the literature are about 1.8 times as high as the lithostatic pressures derived from our cross section restoration. Given that tectonic overpressure in an orogen may be as high as lithostatic pressure (Petrini and Podladchikov 2000), the results of our cross section restoration suggest that the exceptionally high pressures recorded by the Adula Nappe may not be due to exceptionally deep burial but, at least partly, to tectonic overpressure. Engi, M., Berger, A. & Roselle, G.T. 2001: Geology 29, 1143-1146. Flisch, M. 1986: Bull. Ver. Schweiz. Pet.-Geol.-Ing. 53, 23- 49. Froitzheim, N., Pleuger, J., Roller, S. & Nagel, T. 2003: Geology 31, 925-928. Mayerat Demarne, A.M. 1994: Beitr. Geol. Karte Schweiz, 165. Petrini, K. & Podladchikov, Yu. 2000: J. metamorphic Geol.18, 67-77. Schmid, S.M., Pfiffner, O.A., Froitzheim, N., Schönborn, G. & Kissling, E. 1996: Tectonics 15, 1036-1064.

Magnitude of long-term non-lithostatic pressure variations in lithospheric processes: insight from thermo-mechanical subduction/collision models

NASA Astrophysics Data System (ADS)

Gerya, Taras

2014-05-01

On the one hand, the principle of lithostatic pressure is habitually used in metamorphic geology to calculate paleo-depths of metamorphism from mineralogical pressure estimates given by geobarometry. On the other hand, it is obvious that this lithostatic (hydrostatic) pressure principle should only be valid for an ideal case of negligible deviatoric stresses during the long-term development of the entire tectono-metamorphic system - the situation, which newer comes to existence in natural lithospheric processes. The question is therefore not "Do non-lithostatic pressure variations exist?" but " What is the magnitude of long-term non-lithostatic pressure variations in various lithospheric processes, which can be recorded by mineral equilibria of respective metamorphic rocks?". The later question is, in particular, relevant for various types of high-pressure (HP) and ultrahigh-pressure (UHP) rocks, which are often produced in convergent plate boundary settings (e.g., Hacker and Gerya, 2013). This question, can, in particular, be answered with the use of thermo-mechanical models of subduction/collision processes employing realistic P-T-stress-dependent visco-elasto-brittle/plastic rheology of rocks. These models suggest that magnitudes of pressure deviations from lithostatic values can range >50% underpressure to >100% overpressure, mainly in the regions of bending of rheologically strong mantle lithosphere (Burg and Gerya, 2005; Li et al., 2010). In particular, strong undepresures along normal faults forming within outer rise regions of subducting plates can be responsible for downward water suction and deep hydration of oceanic slabs (Faccenda et al., 2009). Weaker HP and UHP rocks of subduction/collision channels are typically subjected to lesser non-lithostatic pressure variations with characteristic magnitudes ranging within 10-20% from the lithostatic values (Burg and Gerya, 2005; Li et al., 2010). The strength of subducted crustal rocks and the degree of confinement of the subduction/collision channel are the key factors controlling this magnitude (Burg and Gerya, 2005; Li et al., 2010). High-temperature (>700 C) UHP rocks formed by continental crust subduction typically demonstrate negligible non-lithostatic pressure variations at peak metamorphic conditions, although these variations can be larger at the prograde stage (Gerya et al., 2008; Li et al., 2010). However, the variability of tectonic mechanisms by which UHP rocks can form (e.g., Sizova et al., 2012; Hacker and Gerya, 2013) precludes generalization of this result for all types of UHP-complexes. References Burg, J.-P., Gerya, T.V. (2005) Viscous heating and thermal doming in orogenic metamorphism: numerical modeling and geological implications. J. Metamorph. Geol., 23, 75-95. Faccenda, M., Gerya, T.V., Burlini, L. (2009) Deep slab hydration induced by bending related variations in tectonic pressure. Nature Geoscience, 2, 790-793. Gerya T.V., Perchuk, L.L., Burg J.-P. (2008) Transient hot channels: perpetrating and regurgitating ultrahigh-pressure, high temperature crust-mantle associations in collision belts. Lithos, 103, 236-256. Hacker, B., Gerya, T.V. (2013) Paradigms, new and old, for ultrahigh-pressure tectonism. Tectonophysics, 603, 79-88. Li, Z., Gerya, T.V., Burg, J.P. (2010) Influence of tectonic overpressure on P-T paths of HP-UHP rocks in continental collision zones: Thermomechanical modelling. J. Metamorphic Geol., 28, 227-247. Sizova, E., Gerya, T., Brown M. (2012) Exhumation mechanisms of melt-bearing ultrahigh pressure crustal rocks during collision of spontaneously moving plates. Journal of Metamorphic Geology, 30, 927-955.