Cumulate nodules as evidence for convective fractionation in a phonolite magma chamber

The physical mechanism by which chemical zonation develops in magma chambers has been controversial partly because unambiguous geological constraints have been lacking. The 11,000 years B.P. eruption of Laacher See Volcano produced a zoned tephra deposit and also ejected crystal-rich nodules which provide a snapshot of the materials crystallising at the magma chamber margins. New data on petrography and chemical compositions of nodules, their cumulate minerals and interstitial glasses are used to deduce the chemical evolution of the phonolite melt due to fractional crystallisation of the mineral assemblages. These data, together with those on the vertical zonation of the melt in the bulk of the chamber, are shown to be consistent with a model of stratification of the chamber by convective fractionation, in which a thin boundary layer of residual melt from fractional crystallisation ascends at the chamber side and accumulates at the roof. Crystallisation could have provided buoyancy to drive convection by enriching incompatible volatile components (mainly water) in the residual melt. Available fluid dynamic studies of single- and double-diffusive boundary layers are used to assess convection in the Laacher See chamber. The boundary layer is likely to have been: (1) laminar, which implies that the density gradient in the chamber steepened upwards; (2) in the counterflow regime, in which compositional and thermal layers flow in opposite directions; and (3) thin ( 10 cm), estimated from theory for a flat wall, suggesting that wall morphology could be important in determining boundary layer characteristics. Estimates of mass transfer rates due to this mechanism suggest that the chamber could have become stratified in a time of the order of 10³ years.     

Co-existing calcic amphiboles in calc-alkaline andesites: Possible evidence of a zoned magma chamber

Hornblende-biotite andesites erupted from Mount Price and Clinker Peak volcanoes, southwestern British Columbia, contain two texturally and compositionally distinct calcic amphiboles: pargasitic hornblende xenocrysts and magnesio-hornblende microphenocrysts. Disequilibrium relationships exhibited by these amphiboles and associated minerals suggest that the magnesio-hornblendes precipitated under chemical and thermal conditions that were intermediate between those under which pargasitic hornblende and biotite, respectively, crystallized. Experimental studies of crystallization in double-diffusive systems (Chen and Turner, 1980; Turner, 1980; McBirney, 1980) suggest that these varied magmatic environments can be explained as a consequence of progressive crystallization within a zoned magma chamber. Although gravitational settling may have played a role, the observed mineral assemblages probably developed by convective mixing of crystals precipitated at the cooling margins with those crystallized in the interior of the compositionally stratified magma column.     

Replenishment of a mafic magma in a zoned felsic magma chamber beneath Rishiri Volcano, Japan

The trachytic Tanetomi lava from Rishiri Volcano, northern Japan, provides useful information concerning how a replenished mafic magma mixes with a compositionally zoned felsic magma in a magma chamber. The Tanetomi lava was erupted in the order of Lower lava 1 (LL1, 59.2-59.8 wt.% in SiO₂), Lower lava 2 (LL2, 58.4-59.1 wt.%), and Upper lava (UL, 59.9-65.1 wt.%). Evidence for mixing with a mafic magma is observed only in the LL2, in which a greater amount of crystals derived from the mafic magma occurs in rocks with higher SiO₂ content. The whole-rock compositional trend of the Tanetomi lavas is fairly smooth except for the LL2 lava composition, which scatter along the main composition trend. There is no reasonable composition of basaltic magma on the extrapolation of the LL2 composition trend, and the trend cannot be explained by a simple two-component magma mixing. Before the replenishment, the felsic magma was zoned in composition (58-65 wt.% in SiO₂) and temperature (1030-920°C) in the magma chamber located at the pressure of ~2 kbar. The compositional variation of the main felsic magma was produced by extraction of a fractionated interstitial melt from mush zones along the chamber walls and its subsequent mixing with the main magma (boundary layer fractionation). The LL1 magma tapped the magma chamber soon after the replenishment, before the mafic magma mixed with the overall felsic magma. Then the basalt magma mixed heterogeneously with the upper part of the felsic magma by forced convection as a fountain during injection. The mixing of the basalt magma with compositionally zoned felsic magma resulted in the characteristic composition trend of the LL2. The fraction of basaltic magma in the LL2 magma is estimated to be at most 10%. Despite such a small proportion, the basalt magma was mixed completely with the felsic magma, probably because the crystallinity of undercooled basalt magma was low enough to behave as a liquid.     

Gravity evidence for an extinct magma chamber beneath Syrtis Major, Mars: a look at the magmatic plumbing system

Syrtis Major is an ancient basaltic shield volcano on Mars with a basal diameter of 1100 km. The free-air gravity anomaly is 126 mGal at spherical harmonic degree 50 and reaches its maximum amplitude over the 2 km deep topographic caldera. The observed gravity anomaly cannot be explained by flexurally supported surface topography and requires the presence of a buried, high-density load. The geologically most reasonable interpretation of this high-density load is that it represents the magma chamber of Syrtis Major, now solidified and filled at least in part by dense igneous cumulates. Pyroxene is likely to be the dominant cumulate mineral in this system, although olivine may also be present. Gravity models presented here define the structure of the buried load and in essence provide a look at the magmatic plumbing system of this volcano. The preferred model involves a buried load that is approximately 300×600 km across, roughly twice as large as the topographic caldera. Both the buried load and the caldera are elongated in the north-south direction. In the center of the buried load, the minimum thickness is 2.8 km for an olivine-dominated cumulate system or 3.9 km for a pyroxene-dominated system. The best terrestrial analog for this structure is the Bushveld Complex, an igneous cumulate body that is similar in size and thickness to the Syrtis Major structure. Assuming that the mean crustal density is 2600 kg m⁻³ due to impact brecciation, the elastic lithosphere at Syrtis Major was 10-15 km thick at the time when the topographic load was emplaced. This corresponds to a lithospheric thermal gradient of 28-52 K/km and a surface heat flux of 70-130 mW m⁻². Higher resolution gravity data, such as that which is planned for the 2005 Mars Reconnaissance Orbiter, will permit further refinement of the dimensions of this structure.     

Caldera formation by magma withdrawal from a reservoir beneath a volcanic edifice

Caldera formation has been explained by magma withdrawal from a crustal reservoir, but little is known about the conditions that lead to the critical reservoir pressure for collapse. During an eruption, the reservoir pressure is constrained to lie within a finite range: it cannot exceed the threshold value for eruption, and cannot decrease below another threshold value such that feeder dykes get shut by the confining pressure, which stops the eruption. For caldera collapse to occur, the critical reservoir pressure for roof failure must therefore be within this operating range. We use an analytical elastic model to evaluate the changes of reservoir pressure that are required for failure of roof rocks above the reservoir with and without a volcanic edifice at Earth's surface. With no edifice at Earth's surface, faulting in the roof region can only occur in the initial phase of reservoir inflation and affects a very small part of the focal area. Such conditions do not allow caldera collapse. With a volcanic edifice, large tensile stresses develop in the roof region, whose magnitude increase as the reservoir deflates during an eruption. The edifice size must exceed a threshold value for failure of the roof region before the end of eruption. The largest tensile stresses are reached at Earth's surface, indicating that faulting starts there. Failure affects an area whose horizontal dimensions depend on edifice and chamber dimensions. For small and deep reservoirs, failure conditions cannot be achieved even if the edifice is very large. Quantitative predictions are consistent with observations on a number of volcanoes.     

日本上地壳S波衰减成像

地球介质的衰减特性可以从地震波中的频谱信息中获取,并且能够提供地球内部结构非均匀性的重要信息.与大量的速度结构研究相比,对衰减结构的研究却较少,且大多仅对日本群岛局部地区做研究.本文运用ML振幅层析成像技术反演获得了整个日本群岛上地壳的品质因子Q0值（1 Hz时的Q值）.计算中共利用了日本Hi-net台网971个台站记录到的5559个地震事件,超过60000 条Sg波最大振幅-频率数据,并且满足震中距小于2°、震源深度小于10 km.计算结果表明,低Q0值区域主要存在于日本群岛的中央火山带,与日本新生代以来的火山分布非常一致;同时,位于从神户到京都东北地区、以及四国岛东部等地的非火山区,断层密集,对应的Q0值也很低;另外,西太平洋俯冲海沟与日本群岛东海岸之间的区域也是低Q0值区,这说明高温、岩石强烈破碎、深厚的沉积物是造成地震波强烈衰减的主要原因.高Q0值主要分布在中央火山带前缘与日本东海岸之间.     

A laboratory model of a replenished magma chamber

It has recently been suggested that periodic influxes of hot but heavy magma into the base of a basaltic magma chamber can remain isolated from the rest of the chamber while the new magma cools and crystallization proceeds. When thermal equilibrium is almost complete, the suspended crystals settle out and the residual, less dense liquid can then mix with the fluid above. In the present paper the basic fluid-dynamical processes underlying this model have been investigated in laboratory experiments using aqueous solutions. The lower layer was hot KNO₃ solution, for which saturated solutions become less dense as the temperature decreases. With a cold, deeper layer of less dense NaNO₃ or K₂CO₃ above the lower layer, there was strong convective transfer of heat through a sharp interface separating the layers, at a rate which is predicted here drawing on previous studies carried out with oceanographic applications in mind. Once crystallization began, non-equilibrium effects became important and the observed temperatures differ somewhat from those predicted. In the experiments crystals grew mainly from the bottom rather than while in suspension, but this is not an essential aspect of the model. The important fact is that the density of the residual liquid in the lower layer decreased until it became equal to that of the upper layer, and then the interface broke down so that the two layers mixed thoroughly together, leaving a layer of KNO₃ crystals at the base. No crystallization at all occurred when the hot input liquid was forced to mix initially with the cold solution already in the chamber.     

Constraints on melting and magma production in the crust

Major intrusions of granitic rocks are found in several tectonic settings and, in all cases, crustal melts may contribute to the volumes of granitic magma. High-grade metamorphism and partial melting of the crust take place predominantly under fluid-absent conditions. We present a model for calculating the amounts of melt that may be formed by fluid-absent breakdown of micas and amphiboles in common crustal rock types (pelitic, quartzofeldspathic, intermediate and mafic). Melt proportions depend mainly on the kind of source rock, the pressure at which melting takes place, the temperature and the hydrous mineral (H₂O) content of the source. As a consequence of the pressure dependence of water solubility in silicate melts, any given source rock will produce more melt, by a given fluid-absent reaction, at lower pressure. At a given pressure, higher-temperature reactions can produce more melt from a given source rock. Based on a survey of the compositions of common rock types, we show that the amounts of melt can vary from < 10to> 50vol.%. Thus, crustal rocks vary widely in their “fertility” as magma sources, depending on the types and amounts of hydrous phases they contain. In general, muscovite breakdown in pelites will yield only small quantities of melt and lead to migmatite formation. Biotite breakdown in pelites occurs at higher temperature and, because most high-grade pelites (below granulite grade) are biotite-rich, can yield up to about 50 vol.% of granitoid melt. Rocks of intermediate composition and hornblende-rich mafic rocks are potentially highly fertile magma sources also, provided that the high temperatures necessary for biotite and hornblende breakdown are realized. Pyroxene-rich mafic rocks and quartzofeldspathic rocks are much less potentially fertile. Data suggest that mechanisms exist for the efficient segregation of melt and restite in systems with < 30and probably< 20vol.% melt. The pressure-temperature history of a region can greatly influence crustal source fertility through its control over the occurrence of subsolidus dehydration and melting equilibria.     

Three-dimensional numerical analysis of magma transport through a pre-existing fracture in the crust

Magmas are transported through pre-existing fractures in many repeatedly erupting volcanoes. The study of this special process of magma transport is fundamentally important to understand the mechanisms and conditions of volcanic eruptions. In this paper, we numerically simulate the magma propagation process through a pre-existing vertical fracture in the crust by using the combined finite difference method (FDM), finite element method (FEM) and discontinuous deformation analysis (DDA) approach. FDM is used to analyze magma flow in the pre-existing fracture, FEM is used to calculate the opening of the fracture during magma intrusion, and DDA is used to deal with the contact of the closed fracture surfaces. Both two-dimensional (2D) and three-dimensional (3D) examples are presented. Parametric studies are carried out to investigate the influence of various physical and geometric parameters on the magma transport in the pre-existing fracture. We have considered magma chamber depth ranging from 7 km to 10 km under the crust surface, magma viscosity ranging from 2 × 10⁻² to 2 × 10⁻⁷ MPa s, and the density difference between the magma and host rock ranging from 300 to 700 kg/m³. The numerical results indicate that (1) the fluid pressure p varies gradually along the depth, (2) the shape of the magma body during propagation is like a torch bar and its width ranges from 2 m to 4 m approximately in the 3D case and 10 m to 50 m in the 2D case for the same physical parameters used, (3) the crust surface around the pre-existing fracture begins to increase on both sides of the fracture, forms a trough between them, then gradually uplifts during the transport of the magma, and finally takes the shape of a crater when the magma reaches the surface. We have also examined the influence of physical and geometric parameters on the minimum overpressure for magma transport in the 3D case. The numerical results show that our numerical technique presented in this paper is an effective tool for simulating magma transport process through pre-existing fractures in the crust.     

Evolution of double diffusion convection in a felsic magma chamber

The onset of double diffusion convection (DDC) is modeled in a two-dimensional case in respect to magma chambers. The viscosity model for the melt takes into account the effects of temperature and concentration of the dissolved component (H₂O). The upper boundary of the convecting magma chamber is assumed to be anhydrous and at constant temperature, whereas the lower boundary is treated as being hydrous permeable with a temperature greater than that within the upper boundary. The case of positive compositional and thermal buoyancy of melt is studied assuming a H₂O diffusion coefficient small in comparison with thermal diffusivity. The DDC has been modeled using a system of equations solved by the finite difference method on a square grid. The convective pattern evolution has been studied for fixed boundary conditions as well as for cooling and degassing. Due to the higher viscosity in the upper zone, the upper boundary layer is thicker than the lower one. The variation of water concentration in this zone of the convective cell can be significant. In nature, the high gradient of water concentration can be responsible for the observed variations of water content in minerals crystallized from a granite melt (e.g., biotite). Because of a high Lewis number (= 100), temperature variations in the magma chamber decay much faster than the water concentration. In this case the intensive convection can continue at a constant temperature due to the non-zero water content in the chamber. In principle, the effect can be applied to the formation of magmatic bodies. If the cooling and degassing system reaches a uniform temperature distribution prior to the crystallization temperature, water content throughout the body may still remain variable.     

The seismic detection of a magma chamber beneath Mount Etna,sicily

Macroseismic studies, linear refraction profiles, and a two-dimensional seismic array study on Mount Etna, have all detected anomalous low velocity zones beneath the volcano. Seismic travel time delays together with high frequency attenuation observations confirm the presence of a large volume of partial melt beneath the volcano. A simplified three-dimensional model of the main storage system of Mount Etna is presented.Paper presented at the Symposium Volcanoes of the Earth and Planets, held at the University of Lancaster, March 17, 1981.     

The roles of magma and groundwater in the phreatic eurptions at Inyo Craters,Long Valley Caldera,California

The Inyo Craters (North Inyo Crater and South Inyo Crater), and a third crater, Summit Crater, are the largest of more than a dozen 650- to 550-yr-B.p. phreatic craters that lie in a 1-km-square area at the south end of the Inyo Volcanic Chain, on the west side of the Long Valley Caldera in eastern California. The three craters are aligned within a 1-km-long northsouth system of fissures and normal faults, and coincide in age with aligned magmatic vents farther north in the Inyo Volcanic Chain, suggesting that they were all produced by intrusion of one or more dikes. To study the sequence and mechanisms of the eruptions, the deposits were mapped, sampled, and compared with subsurface stratigraphy obtained from the core of a slant hole drilled directly below the center of South Inyo Crater from the southwest. The deposits from the two Inyo Craters are fine-grained (median diameter less than 1 mm), are several meters thick at the crater walls, and cover at most a few km² of ground surface. Stratigraphic relationships between the Inyo Craters and Summit Crater indicate that the eruptions proceeded from north to south, overlapped slightly in time, and produced indistinctly plane-parallel bedded, poorly sorted deposits, containing debris derived primarily from within 450 m of the surface. Debris from the deepest identifiable unit (whose top is at 450 m depth) is present at the very base of both Inyo Craters deposits, suggesting that the eruptive vents were open and tapping debris from at least that depth, probably along preexisting fractures, even at their inception. According to ballistic studies, the greatest velocity of ejected blocks was of the order of 100 m/s. All eruptions, particularly the least powerful, selectively removed debris from the finest-grained, most easily eroded subsurface units. Although juvenile fragments have been previously identified in these deposits, they are confined primarily to the grain-size fraction smaller than 0.25 mm dia. and probably did not constitute more than several percent of the deposit. It is therefore suggested that these juvenile fragments were not the main source of heat for the eruptions, and that the eruptions were caused either by: (1) heating of water by fragmented magma that was not ejected before the eruption shut off; (2) slow heating (over months to years) of groundwater under confined conditions without fragmentation of magma, followed by a second process (pressure buildup, seismic faulting, or intrusions) that breached the confinement; or (3) breach of a pre-existing confined geothermal aquifer.     

Crystal retention,fractionation and crustal assimilation in a convecting magma chamber,Nisyros Volcano,Greece

Nisyros island is a calc-alkaline volcano, built up during the last 100 ka. The first cycle of its subaerial history includes the cone-building activity with three phases, each characterized by a similar sequence: (1) effusive and explosive activity fed by basaltic andesitic and andesitic magmas; and (2) effusive andextrusive activity fed by dacitic and rhyolitic magmas. The second eruptive cycle includes the caldera-forming explosive activity with two phases, each consisting of the sequence: (1) rhyolitic phreatomagmatic eruptions triggering a central caldera collapse; and (2) extrusion of dacitic-rhyolitic domes and lava flows. The rocks of this cycle are characteized by the presence of mafic enclaves with different petrographic and chemical features which testify to mixing-mingling processes between variously evolved magmas. Jumps in the degree of evolution are present in the stratigraphic series, accompanied by changes in the porphyritic index. This index ranges from 60% to about 5% and correlates with several teochemical parameters, including a negative correlation with Sr isotope ratios (0.703384–0.705120). The latter increase from basaltic andesites to intermediate rocks, but then slightly decrease in the most evolved volcanic rocks. The petrographic, geochemical and isotopic characteristics can be largely explained by processes occurring in a convecting, crystallizing and assimilating magma chamber, where crystal sorting, retention, resorption and accumulation take place. A group of crystal-rich basaltic andesites with high Sr and compatible element contents and low incompatible elements and Sr isotope ratios probably resulted from the accumulation of plagioclase and pyroxene in an andesitic liquid. Re-entrainment of plagioclase crystals in the crystallizing magma may have been responsible for the lower ⁸⁷Sr/⁸⁶Sr in the most evolved rocks. The gaps in the degree of evolution with time are interpreted as due to liquid segregation from a crystal mush once critical crystallinity was reached. At that stage convection halted, and a less dense, less porphyritic, more evolved magma separated from a denser crystal-rich magma portion. The differences in incompatible element enrichment of pre-and post-caldera dacites and the chemical variation in the post-caldera dome sequence are the result of hybridization of post-caldera dome magmas with more mafic magmas, as represented by the enclave compositions. The occurrence of the quenched, more mafic magmas in the two post-caldera units suggests that renewed intrusion of mafic magma took place after each collapse event.     

Timing of magma extraction during the Campanian Ignimbrite eruption (Campi Flegrei Caldera)

A core drilled within the northern part of the city of Napoli has offered the unique opportunity to observe in one single sequence the superposition of the four pyroclastic flow units emplaced during the Campanian Ignimbrite (CI) eruption. Such a stratigraphic succession has never been encountered before in natural or in man made exposures. Therefore the CI sequence was reconstructed only on the basis of stratigraphic correlations and compositional data (in literature). The occurrence of four superposed CI flows, together with all the data available (in literature) allowed us to better constrain the chemical stratigraphy of the deposit and the compositional structure of the CI magma chamber. The CI magma chamber includes two cogenetic magma layers, separated by a compositional gap. The upper magma layer was contaminated by interaction with radiogenic fluids. The two magma layers were extruded either individually or simultaneously during the course of the eruption. In the latter case they produced a hybrid magma. But no evidence of input of new geochemically and isotopically distinct magma batches just prior or during the eruption has been found. Comparison with the exposed CI deposits has permitted reconstruction of variable eruption phases and related magma withdrawal and caldera collapse episodes. The eruption was likely to have began with phreatomagmatic explosions followed by the formation of a sustained plinian eruption column fed by the simultaneous extraction from both magma layers. Towards the end of this phase the upward migration of the fragmentation surface and the decrease in magma eruption rate and/or activation of fractures formed an unstable pulsating column that was fed only by the most-evolved magma layer. This plinian phase was followed by the collapse of the eruption column and the beginning of caldera formation. At this stage expanded pyroclastic flows fed by the upper magma layer in the chamber generated. During the following major caldera collapse episode, the maximum mass discharge rate was reached and both magma layers were tapped, generating expanded pyroclastic flows. Towards the end of the eruption, only the deeper and less differentiated magma layer was tapped producing more concentrated pyroclastic flows that traveled short distances.     

Incipient magma chamber formation as a result of repetitive intrusions

An analytical solution for periodic magma intrusions in conduits was developed to study the onset of shallow magma chamber formation. The solution is based on determining when a repetitive series of intrusions can cause the wall rock of a conduit to reach its melt temperature. The results show that magma chamber formation in conduits is a strong function of the volume rate of intrusion and that magma chamber formation is likely when the intrusion rate exceeds 10^?3 km³/ yr. which agrees with observations by other investigators. Once this critical value of intrusion rate is reached, magma chambers are likely to begin forming after only a few intrusive pulses (less than ten). Results for both cylindrical conduits and dikes show cylindrical conduits are more favourable for the formation of shallow magma chambers.     

Effusive eruptions from a large silicic magma chamber: the Bearhead Rhyolite,Jemez volcanic field,NM

Large continental silicic magma systems commonly produce voluminous ignimbrites and associated caldera collapse events. Less conspicuous and relatively poorly documented are cases in which silicic magma chambers of similar size to those associated with caldera-forming events produce dominantly effusive eruptions of small-volume rhyolite domes and flows. The Bearhead Rhyolite and associated Peralta Tuff Member in the Jemez volcanic field, New Mexico, represent small-volume eruptions from a large silicic magma system in which no caldera-forming event occurred, and thus may have implications for the genesis and eruption of large volumes of silicic magma and the long-term evolution of continental silicic magma systems.⁴⁰Ar/³⁹Ar dating reveals that most units mapped as Bearhead Rhyolite and Peralta Tuff (the Main Group) were erupted during an ∼540 ka interval between 7.06 and 6.52 Ma. These rocks define a chemically coherent group of high-silica rhyolites that can be related by simple fractional crystallization models. Preceding the Main Group, minor amounts of unrelated trachydacite and low silica rhyolite were erupted at ∼11–9 and ∼8 Ma, respectively, whereas subsequent to the Main Group minor amounts of unrelated rhyolites were erupted at ∼6.1 and ∼1.5 Ma.The chemical coherency, apparent fractional crystallization-derived geochemical trends, large areal distribution of rhyolite domes (∼200 km²), and presence of a major hydrothermal system support the hypothesis that Main Group magmas were derived from a single, large, shallow magma chamber. The ∼540 ka eruptive interval demands input of heat into the system by replenishment with silicic melts, or basaltic underplating to maintain the Bearhead Rhyolite magma chamber.Although the volatile content of Main Group magmas was within the range of rhyolites from major caldera-forming eruptions such as the Bandelier and Bishop Tuffs, eruptions were smaller volume and dominantly effusive. Bearhead Rhyolite domes occur at the intersection of faults, and are cut by faults, suggesting that the magma chamber was structurally vented preventing volatiles from accumulating to levels high enough to trigger a caldera-forming eruption.     

Experimental evidence for the role of accessory phases in magma genesis

Recent experimental studies have established petrogenetic models based on melting processes involving major phases. The possible residual character of trace-element-enriched accessory phases is not considered for temperatures well above the solidus in these models. In contrast, geochemists, applying trace element data to independently test the experimentally-based models, have concluded that residual (or fractionating) accessory phases may have an essential role in controlling the trace element (especially REE) distributions in magmas.Some recent experimental work provides data on the stability of potentially significant accessories such as sphene, rutile, apatite, zoisite and mica in basaltic compositions at elevated P and T. Sphene is stable to 1000°C with 60% melting of a hydrous tholeiite at 15 kbar. At higher pressure, rutile is the only Ti-rich accessory phase, and is present to at least 1000°C and high degrees of melting. Published REE data on sphene and rutile suggest that these phases may be important in controlling REE distribution in some magmas. For example, island are high-Mg, low-Ca-Ti tholeiites with low REE abundances and U-shaped patterns (Hickey and Frey, 1979) may reflect the role of sphene. In addition to rutile, similar close-packed Ti-rich accessory phases such as priderite, perovskite, crichtonite and loveringite may occur in mantle-derived magmas. These phases readily accommodate the REE but their possible role needs experimental confirmation.Apatite is recorded in hawaiite (1.16% P₂O_s) with 2% H₂O added at 5–6 kbar and 1050°C within 30°C of the liquidus, but at present no other experimental data are available on its high P, T stability, although thermodynamic calculations indicate that F may increase its stability markedly. Apatite is well known in high-pressure inclusions and as a phenocryst phase in rocks of the alkaline and calc-alkaline series.Ilmenite is known as a near-liquidus phase in some mafic magmas at 5–10 kbar, but its stability decreases to near-solidus at 25–30 kbar. Zoisite occurs in hydrous mafic compositions at mantle pressures, but it is confined to temperatures < 780°C. Finally, mica has a wide temperature range of stability at mantle pressures, especially in potassic magmas, and phlogopitic mica is stable to 1040°C at 20–25 kbar in a hydrous, K-rich “tholeiite” (1.6% K₂O).     

Silicic magma entering a basaltic magma chamber: eruptive dynamics and magma mixing — an example from Salina (Aeolian islands,Southern Tyrrhenian Sea)

The Pollara tuff-ring resulted from two explosive eruptions whose deposits are separated by a paleosol 13 Ka old. The oldest deposits (LPP, about 0.2 km³) consist of three main fall units (A, B, C) deposited from a subplinian column whose height (7–14 km) increased with time from A to C, as a consequence of the increased magma discharge rate during the eruption (1–8x10⁶ kg/s). A highly variable juvenile population characterizes the eruption. Black, dense, highly porphyritic, mafic ejecta (SiO₂=50–55%) almost exclusively form A deposits, whereas grey, mildly vesiculated, mildly porphyritic pumice (SiO₂=56–67%) and white, highly vesiculated, nearly aphyric pumice (SiO₂=66–71%) predominate in B and C respectively. Mafic cumulates are abundant in A, while crystalline lithic ejecta first appear in B and increase upward. The LPP result from the emptying of an unusual and unstable, compositionally zoned, shallow magma chamber in which high density mafic melts capped low density salic ones. Evidence of the existence of a short crystal fractionation series is found in the mafic rocks; the andesitic pumice results from complete blending between rhyolitic and variously fractionated mafic melts (salic component up to 60 wt%), whereas bulk dacitic compositions mainly result from the presence of mafic xenocrysts within rhyolitic glasses. Viscosity and composition-mixing diagrams show that blended liquids formed when the visosities of the two end members had close values. The following model is suggested: 1. A rhyolitic magma rising through the metamorphic basement enterrd a mafic magma chamber whose souter portions were occupied by a highly viscous, mafic crystal mush. 2. Under the pressure of the rhyolitic body the nearly rigid mush was pushed upwards and mafic melts were squeezed against the walls of the chamber, beginning roof fracturing and mingling with silicic melts. 3. When the equilibrium temperature was reached between mafic and silicic melts, blended liquids rapidly formed. 4. When fractures reached the surface, the eruption began by the ejection of the mafic melts and crystal mush (A), followed by the emission of variously mingled and blended magmas (B) and ended by the ejection of nearly unmixed rhyolitic magma (C).     

Life-time of a stratified magma chamber recorded in ultramafic xenoliths from Ichinomegata volcano,northeastern Japan

Variations in major-element chemistry and modal compositions of the mantle xenolith-bearing calc-alkalic ejecta from Ichinomegata volcano are inferred to be due to mixing of three magmatic end members: Basalt I (SiO₂ 51 wt% , MgO 8.5 wt%), Basalt II (SiO₂ 54 wt%, MgO 5 wt%), and Dacite (SiO₂ 65 wt%, MgO₂ wt%). Ultramafic xenoliths are found in mafic mixtures of Dacite-Basalt I and Dacite-Basalt II. The thermal histories of the xenoliths in both mixtures are compared with each other. Chemical compositions of olivine and orthopyroxene in xenoliths suggest that xenoliths in Basalt I were equilibrated at about 800 °C, while those in Basalt II were also equilibrated originally at about 800 °C but were subsequently annealed at about 1000 °C for more than 10^2–3 years prior to the eruption.The chemical composition of Basalt I indicates that it can coexist with upper mantle peridotite and it is an appropriate candidate for a carrier of ultramafic xenoliths from the upper mantle. On the other hand, Basalt II is fractionated and it cannot be directly derived from the upper mantle. Two pulses of xenolith-bearing basalt injection into a dacite magma chamber are inferred to have occurred. The first injection did not lead to eruption and subsequently formed a dacite/basalt stratified magma chamber. In the lower layer, the basalt was slightly differentiated to become Basalt II and ultramafic xenoliths carried by the first pulse were annealed at the bottom of the layer. The duration of the annealing of the xenoliths implies a minimum life-time of the Dacite-Basalt II stratification in the magma chamber beneath Ichinomegata of 10^2–3 years. The second injection of the xenolith-bearing basalt (Basalt I) was immediately followed by eruption, and all the magmas were effused with mixing in a conduit. Consequently, the ultramafic xenoliths carried by the second pulse are not annealed.