The Magmatic-hydromagmatic Transition During Explosive Eruptions: Highlights From Mwi (magma- Water Interaction ) Laboratory Experiments

NASA Astrophysics Data System (ADS)

Mincione, V.; Trigila, R.

The transition from a mere magmatic to a hydromagmatic regime has been frequently recognized for basaltic explosive eruptions. Indeed, due to great efficiency of MWI to transform the magma thermal energy in mechanical energy, the detection of this transition is particularly relevant for the evaluation of volcanic hazard. The study of this transition is being realized using the MAQUA pressure vessel built on purpose for MWI experiments. With this setup the distribution of the produced mechanical energy among magma fragmentation, system expansion and transport of fragmented material has been measured. By means of two different experimental procedures, the fragmentation is achieved either because of water exsolution (magmatic regime) and as the result of magma-water interaction (hydromagmatic regime). In both cases, acoustic microphones on input and output pressure tubings have been connected with our data aquisition system allowing the record of acoustic signals. During the magmatic fragmentation some acoustic waves were emitted, whereas after a few seconds following the water injection a shockwave was recorded. In the latter case, special particles are produced as the result of the magma-water interaction ("interactive particles" Zimanowski et al., JGR, 102, b1, 803-814, 1997). More experiments are in progress in order to better constrain the occurrence of the shock wave associated with the explosive MWI.

MAGMA: analysis of two-channel microarrays made easy.

PubMed

Rehrauer, Hubert; Zoller, Stefan; Schlapbach, Ralph

2007-07-01

The web application MAGMA provides a simple and intuitive interface to identify differentially expressed genes from two-channel microarray data. While the underlying algorithms are not superior to those of similar web applications, MAGMA is particularly user friendly and can be used without prior training. The user interface guides the novice user through the most typical microarray analysis workflow consisting of data upload, annotation, normalization and statistical analysis. It automatically generates R-scripts that document MAGMA's entire data processing steps, thereby allowing the user to regenerate all results in his local R installation. The implementation of MAGMA follows the model-view-controller design pattern that strictly separates the R-based statistical data processing, the web-representation and the application logic. This modular design makes the application flexible and easily extendible by experts in one of the fields: statistical microarray analysis, web design or software development. State-of-the-art Java Server Faces technology was used to generate the web interface and to perform user input processing. MAGMA's object-oriented modular framework makes it easily extendible and applicable to other fields and demonstrates that modern Java technology is also suitable for rather small and concise academic projects. MAGMA is freely available at www.magma-fgcz.uzh.ch.

Hypersolidus geothermal energy from the moving freeze-fracture-flow boundary

NASA Astrophysics Data System (ADS)

Carrigan, Charles; Eichelberger, John; Sigmundsson, Freysteinn; Papale, Paolo; Sun, Yunwei

2014-05-01

Rhyolitic magmas at low pressure undergo much of their crystallization over a small temperature interval just above the solidus. This hypersolidus material has a high energy density and effective heat capacity because of stored heat of crystallization, yet may sustain fractures and therefore admit heat exchange with fluids because of its interlocking crystal framework. Rhyolitic magmas emplaced near the liquidus should at first cool rapidly, owing to internal convection, modest crystallization with declining temperature, and extreme temperature gradients at their boundaries. However, once the solidus is approached the rapid rise in effective heat capacity should result in low temperature gradients and rates of heat flow within the bodies. They are suspended for a time in the hypersolidus state. Prodigious quantities of heat can be released from these thermal masses by hydrothermal systems, natural or perhaps stimulated, fracturing their way inward from the margins. The fracture front drives the solidus isotherm ahead of it. Heat of crystallization in front of the advancing solidus is transferred across the thin, moving boundary zone to the external fluid, which advects it away. Once the material is below (outboard of) the solidus, it behaves as normal rock and cools rapidly, having a heat capacity only about 20% that of water. Variations on this theme were published by Lister (1974) for mid-ocean ridges, Hardee (1980) for lava lakes, and Bjornsson et al (1982) for Grimsvotn and Heimaey, who cited possible geothermal energy exploitiation. This scenario is consistent with a number of observations: 1. The geophysical rarity of imaging mostly liquid magma in the shallow crust, despite common petrologic evidence that silicic magma has undergone shallow storage. 2. More common imaging of "partial melt" volumes, whose inferred properties suggest some, but not dominant proportion of melt. 3. Evidence that pure-melt rhyolitic eruptions may have drained relatively shallow hypersolidus plutons. 4. Downward propagating thin conductive boundary zone observed in repeated coring of Kilauea Iki lava lake, Hawaii 5. Record enthalpy flow and temperature during flow-testing of Iceland Deep Drilling Project (IDDP)-1 in Krafla Caldera by Landsvirkjun Co. Production came from a 2.1-km-deep 500oC "magma" contact zone, from the vicinity of which fresh rhyolite glass-bearing felsite and crystal-poor rhyolite glass fragments were recovered. The hypothesis of a moving freeze-fracture-flow boundary raises the possibility of ultra-high-temperature, natural or engineered geothermal systems in volcanic areas. We believe that this prospect, as well as the benefit to understanding volcanic hazards at restless calderas, gives merit to further exploration of the hypersolidus regime beneath Krafla Caldera.

Coupling Thermal and Chemical Signatures of Crustal Magma Bodies: Energy-Constrained Eruption, Recharge, Assimilation, and Fractional Crystallization (E'RAχFC)

NASA Astrophysics Data System (ADS)

Bohrson, W. A.; Spera, F. J.

2004-12-01

Energy-Constrained Eruption, Recharge, Assimilation and Fractional Crystallization (E'RAχFC) tracks the evolution of an open-system magmatic system by coupling conservation equations governing energy, mass and species (isotopes and trace elements). By linking the compositional characteristics of a composite magmatic system (host magma, recharge magma, wallrock, eruptive reservoir) to its mass and energy fluxes, predictions can be made about the chemical evolution of systems characterized by distinct compositional and thermal characteristics. An interesting application of E'RAχFC involves documenting the influence distinct thermal regimes have on the chemical evolution of magmatic systems. Heat transfer between a magma-country rock system at epizonal depths can be viewed as a conjugate heat transfer problem in which the average country rock-magma boundary temperature, Tb, is governed by the relative vigor of hydrothermal convection in the country rock vs. magma convection. For cases where hydrothermal circulation is vigorous and magmatic heat is efficiently transported away from the boundary, contact aureole temperatures (~Tb) are low. In cases where magmatic heat can not be efficiently transported away from the boundary and hydrothermal cells are absent or poorly developed, Tb is relatively high. Simultaneous solution of the differential equations governing momentum and energy conservation and continuity for the coupled hydrothermal-magmatic conjugate heat transfer system enables calculation of the characteristic timescale for EC-RAFC evolution and development of hydrothermal deposits as a function of material and medium properties, sizes of systems and relative efficiency of hydrothermal vs. magmatic heat transfer. Characteristic timescales lie in the range 102-106 yr depending on system size, magma properties and permeability among other parameters. In E'RAχFC, Tb is approximated by the user-defined equilibration temperature, Teq, which is the temperature at which all parts of the composite magmatic system achieve thermal equilibrium. Comparison of the results of three EC-AFC simulations at different Teq (1150° C, 1050° C, 1000° C) for a mafic magma intruding middle-upper crust of mafic-intermediate composition illustrate the distinctions that can be imparted by a range of thermal regimes. Model parameters relevant to the following results include: initial Sr concentration, isotope composition and bulk D for host magma are 700 ppm, 0.7035, and 1.5, respectively; those for wallrock are 230 ppm, 0.7100, 0.05. The 1150° C case (i.e., high Tb) yields the least crust-like Sr isotope signatures. The mass of wallrock that reaches thermal equilibrium is relatively small (0.26, normalized to the mass of initial host magma), although the degree of melting is high (97%). In contrast, the 1000° C case (i.e., low Tb) yields the most crust-like Sr isotope signatures. This case is also characterized by the largest mass of wallrock (0.98, normalized to the mass of initial host magma) that achieves thermal equilibrium, but the degree to which this wallrock melts is small (10%). A fundamental issue that derives from these results is the relationship between the chemical evolution of the hydrothermal system and the chemical evolution of associated melt and cumulates. In particular, to what extent can predictions be made from the thermal interactions between magma and wallrock on the chemical signatures of the associated magmatic rocks and hydrothermal deposits?

Lunar magma transport phenomena

NASA Technical Reports Server (NTRS)

Spera, Frank J.

1992-01-01

An outline of magma transport theory relevant to the evolution of a possible Lunar Magma Ocean and the origin and transport history of the later phase of mare basaltic volcanism is presented. A simple model is proposed to evaluate the extent of fractionation as magma traverses the cold lunar lithosphere. If Apollo green glasses are primitive and have not undergone significant fractionation en route to the surface, then mean ascent rates of 10 m/s and cracks of widths greater than 40 m are indicated. Lunar tephra and vesiculated basalts suggest that a volatile component plays a role in eruption dynamics. The predominant vapor species appear to be CO CO2, and COS. Near the lunar surface, the vapor fraction expands enormously and vapor internal energy is converted to mixture kinetic energy with the concomitant high-speed ejection of vapor and pyroclasts to form lunary fire fountain deposits such as the Apollo 17 orange and black glasses and Apollo 15 green glass.

Thermohydrodynamic model: Hydrothermal system, shallowly seated magma chamber

NASA Astrophysics Data System (ADS)

Kiryukhin, A. V.

1985-02-01

The results of numerical modeling of heat exchange in the Hawaiian geothermal reservoir demonstrate the possibility of appearance of a hydrothermal system over a magma chamber. This matter was investigated in hydrothermal system. The equations for the conservation of mass and energy are discussed. Two possible variants of interaction between the magma chamber and the hydrothermal system were computated stationary dry magma chamber and dry magma chamber changing volume in dependence on the discharge of magma and taking into account heat exchange with the surrounding rocks. It is shown that the thermal supplying of the hydrothermal system can be ensured by the extraction of heat from a magma chamber which lies at a depth of 3 km and is melted out due to receipt of 40 cubic km of basalt melt with a temperature of 1,300 C. The initial data correspond with computations made with the model to the temperature values in the geothermal reservoir and a natural heat transfer comparable with the actually observed values.

Petrogenesis of Mount Rainier andesite: magma flux and geologic controls on the contrasting differentiation styles at stratovolcanoes of the southern Washington Cascades

USGS Publications Warehouse

Sisson, Thomas W.; Salters, V.J.M.; Larson, P.B.

2013-01-01

The dominant cause of magmatic evolution at Mount Rainier, however, is inferred to be a version of in situ crystallization-differentiation and mixing (Langmuir, 1989) wherein small magma batches stall as crustal intrusions and solidify extensively, yielding silicic residual liquids with trace element concentrations influenced by accessory mineral saturation. Subsequent magmas ascending through the intrusive plexus entrain and mix with the residual liquids and low-degree re-melts of those antecedent intrusions, producing hybrid andesites and dacites. Mount St. Helens volcanic rocks have geochemical similarities to those at Mount Rainier, and may also result from in situ differentiation and mixing due to low and intermittent long-term magma supply, accompanied by modest crustal assimilation. Andesites and dacites of Mount Adams isotopically overlap the least contaminated Mount Rainier magmas and derive from similar parental magma types, but have trace element variations more consistent with progressive crystallization-differentiation, probably due to higher magma fluxes leading to slower crystallization of large magma batches, allowing time for progressive separation of minerals from melt. Mount Adams also sits atop the southern projection of a regional anticlinorium, so Eocene sediments are absent, or are at shallow crustal levels, and so are cold and difficult to assimilate. Differences between southwest Washington stratovolcanoes highlight some ways that crustal geology and magma flux are primary factors in andesite generation.

Magmagenesis at Soufriere volcano St Vincent, Lesser Antilles Arc

USGS Publications Warehouse

Heath, E.; Macdonald, R.; Belkin, H.; Hawkesworth, C.; Sigurdsson, Haraldur

1998-01-01

Soufriere volcano of St Vincent (3 wt %, whereas various projections onto phase diagrams are more consistent with relatively anhydrous magmas. Primary magmas at Soufriere were generated by around 15% melting of mid-ocean ridge basalt type mantle sources which had been modified by addition of fluids released from the slab containing contributions from subducted sediments and mafic crust.

Lithium deposits hosted in intracontinental rhyolite calderas

NASA Astrophysics Data System (ADS)

Benson, T. R.; Coble, M. A.; Mahood, G. A.

2017-12-01

Lithium (Li) is classified as a technology-critical element due to the increasing demand for Li-ion batteries, which have a high power density and a relatively low cost that make them optimal for energy storage in mobile electronics, the electrical power grid, and hybrid and electric vehicles. Given that many projections for Li demand exceed current economic reserves and the market is dominated by Australia and Chile, discovery of new domestic Li resources will help diversify the supply chain and keep future technology costs down. Here we show that lake sediments preserved within intracontinental rhyolite calderas have the potential to host Li deposits on par with some of the largest Li brine deposits in the world. We compare Li concentrations of rhyolite magmas formed in a variety of tectonic settings using in situ SHRIMP-RG measurements of homogenized quartz-hosted melt inclusions. Rhyolite magmas that formed within thick, felsic continental crust (e.g., Yellowstone and Hideaway Park, United States) display moderate to extreme Li enrichment (1,500 - 9,000 ppm), whereas magmas formed in thin crust or crust comprised of accreted arc terranes (e.g., Pantelleria, Italy and High Rock, Nevada) contain Li concentrations less than 500 ppm. When the Li-enriched magmas erupt to form calderas, the cauldron depression serves as an ideal catchment within which meteoric water that leached Li from intracaldera ignimbrite, nearby outflow ignimbrite, and caldera-related lavas can accumulate. Additional Li is concentrated in the system through near-neutral, low-temperature hydrothermal fluids circulated along ring fractures as remnant magma solidifies and degasses. Li-bearing hectorite and illite clays form in this alteration zone, and when preserved in the geological record, can lead to a large Li deposit like the 2 Mt Kings Valley Li deposit in the McDermitt Caldera, Nevada. Because more than 100 large Cenozoic calderas occur in the western United States that formed on eruption of moderately to extremely Li-rich magma, previously unidentified Li resources may occur in the largest calderas young enough to preserve hydrothermally-altered caldera lake sediments.

Probing the Source of Explosive Volcanic Eruptions (Sergey Soloviev Medal Lecture)

NASA Astrophysics Data System (ADS)

Eichelberger, John C.

2015-04-01

What if we knew where magma is located under a volcano and its current state? Such information would transform volcanology. For extreme events, we typically know where the vulnerabilities are: people, lifelines, and critical infrastructure, but seldom do we know the 'source term' beforehand. For restless calderas such as Campi Flegrei, Italy and Yellowstone, USA, the threat is silicic magma within the caldera itself. Great effort has gone into finding such bodies through surface measurements. 'Discovery' is declared when consensus is achieved. But there is a difference between consensus and knowledge. By following certain conventions in finding magma bodies (aseismic volume, seismic attenuation, Mogi source location, water and CO2 content of melt inclusions) and depicting them in accepted ways (oblate spheroids or lenses with an impossible solid/liquid boundary discontinuity), we perpetuate myths that mislead even ourselves. The consensus view of the Long Valley Caldera, USA, magma reservoir has evolved over 40 years from a 104 km3 balloon to two tiny pockets of magma, in part because drilling revealed a temperature of 100°C at 3 km depth over the 'balloon'. Oil and gas exploration is free of fanciful reservoirs because there is ground truth. Geophysics and geology define a possible reservoir and a well is drilled. If oil is not there, the model needs revision. The situation is worse for conditions of magma storage. The heretofore-unknowable roof zone of magma chambers has been invoked for separating melt from crystals and/or for accumulating vapor and evolved magma leading to eruption. Anything is possible when there are no data. The accidental (but technically remarkable) drilling discovery of rhyolite magma at 2,100 m depth under Krafla Caldera, Iceland by Landsvirkjun Co. and the Iceland Deep Drilling Project opens the door to properly detect magma and to understand how magma evolves, energizes hydrothermal systems, and erupts. A new project before the International Continental Scientific Drilling Program (ICDP) would continuously core through the margin of the magma body, accompanied by state-of-the-art geophysics, geochemical analyses and 3-D mass/heat transport modeling. Coring of molten rock has been conducted with success in lava lakes. Gradients in phase assemblage and composition will provide definitive tests of models of mass/heat transfer and magma evolution. By knowing 'the answer', techniques for finding magma will likewise be tested, making Krafla an international magma laboratory. In fact, Krafla may resemble the state of neighboring Askja Caldera system prior its 1875 eruption, with hidden rhyolite being brewed in a basalt-fired caldera crucible. Additionally, the observed high permeability and sustained power output from the magma body's margin implies self-sustained thermal fracturing, i.e. an 'Enhanced Geothermal System' an order of magnitude more powerful than conventional geothermal. The cost should be balanced against the higher cost of ignorance. For tsunamis, Sergey Soloviev showed there is no substitute for direct measurements at depth, despite technical and economic obstacles. He also led the way in Russian - American cooperation on natural hazards, thereby mitigating the risk of the ultimate hazard, of humans to each other.

Jurassic magmatism in Dronning Maud Land: synthesis of results of the MAMOG project

USGS Publications Warehouse

Leat, P.T.; Curtis, M.L.; Riley, T.R.; Ferraccioli, F.

2007-01-01

The Jurassic Karoo large igneous province (LIP) of Antarctica, and its conjugate margin in southern Africa, is critical for investigating important questions about the relationship of basaltic LIPs to mantle plumes. Detailed aerogeophysical, structural, anisotropy of magnetic susceptibility (AMS), geochronological and geochemical investigations completed under the British Antarctic Survey’s MAMOG project have provided some of the answers. Across most of the area, magma volumes were small compared to those in southern Africa. Jurassic dikes intruding the Archean craton are sparse and the Jutulstraumen trough, a Jurassic rift, is interpreted, from aerogeophysical data, as largely amagmatic. The largest volumes of magma were emplaced along the margin of the craton and close to the Africa-Antarctica rift. Although dikes were emplaced by both vertical and horizontal flow, overwhelmingly magmas in Dronning Maud Land were locally derived, and not emplaced laterally from distant sources. Basaltic magmatism was protracted in Dronning Maud Land (several dike emplacement episodes between ~206 and 175 Ma), and the small magma volumes resulted in highly diverse magma compositions, including picrites and ferropicrites interpreted to have been derived from hot mantle in a mantle plume. The protracted magmatism before the locally ~177 Ma flood lava eruptions, and evidence for a radiating dike swarm, favor a model of mantle plume incubation for 20-30 million years before flood lava eruption.

Magma Energy Research, 79-1. Semiannual report, October 1, 1978-March 31, 1979

DOE Office of Scientific and Technical Information (OSTI.GOV)

Traeger, R.K.; Colp, J.L.; Neel, R.R.

1979-07-01

A major effort in evaluating Kilauea Iki lava lake has been completed. The physical model based on FY 76 geophysical experiments is not correct in that a low viscosity, liquid lens of appreciable thickness does not exist. Mathematical models of the cooling of the lava lake and the state of solidification of the liquid lens were verified by thermal profile and permeability measurements. New jet-augmented drilling concepts successfully penetrated the viscous, multi-phase molten rock region in some locations where conventional drilling failed. Heat transfer studies in the lake suggest injection of fluids to enhance convection may be useful to extractmore » energy from magma chamber margins. Other activities resulted in the completion and successful testing of a 800 cc simulation facility for evaluating simulated magma properties at temperatures to 1500/sup 0/C and pressures to 4 kbar. In materials compatibility studies, thermodynamic stability diagrams were developed for 15 pure metals in basaltic magma systems and compatibility tests completed. Results are being used to define simple alloy systems which may be compatible with magmas and to identify other superalloy materials candidates.« less

Quantitative models for magma degassing and ground deformation (bradyseism) at Campi Flegrei, Italy: Implications for future eruptions

USGS Publications Warehouse

Bodnar, R.J.; Cannatelli, C.; de Vivo, B.; Lima, A.; Belkin, H.E.; Milia, A.

2007-01-01

Campi Flegrei (Phlegrean Fields) is an active volcanic center near Naples, Italy. Numerous eruptions have occurred here during the Quaternary, and repeated episodes of slow vertical ground movement (bradyseism) have been documented since Roman times. Here, we present a quantitative model that relates deformation episodes to magma degassing and fracturing at the brittle-ductile transition in a magmatic-hydrothermal enviromnent. The model is consistent with field and laboratory observations and predicts that uplift between 1982 and 1984 was associated with crystallization of ???0.83 km3 of H2O-saturated magma at 6 km depth. During crystallization, ???6.2 ?? 1010 kg of H2O and 7.5 ?? 108 kg of CO2, exsolved from the magma and generated ???7 ?? 1015 J of mechanical (P??V) energy to drive the observed uplift. For comparison, ???1017 J of thermal energy was released during the 18 May 1980 lateral blast at Mount St. Helens. ?? 2007 The Geological Society of America.

Residual glasses and melt inclusions in basalts from DSDP Legs 45 and 46 - Evidence for magma mixing. [Deep Sea Drilling Project

NASA Technical Reports Server (NTRS)

Dungan, M. A.; Rhodes, J. M.

1978-01-01

Microprobe analyses of natural glasses in basalts recovered by Legs 45 and 46 of the Deep Sea Drilling Project are reported and interpreted in the context of other geochemical, petrographic and experimental data on the same rocks (Rhodes et al., 1978). Residual glass compositions in the moderately evolved aphyritic and abundantly phyric basalts within each site indicate that none of the units is related to any other or to a common parent by simple fractional crystallization. The compositional trends, extensive disequilibrium textures in the plagioclase phenocrysts and the presence in evolved lavas of refractory plagioclase and olivine phenocrysts bearing primitive melt inclusions provide evidence that magma mixing had a major role in the genesis of the Leg 45 and 46 basalts. The magma parental to these basalts was most likely characterized by high Mg/(Mg + Fe/+2/), CaO/Al2O3, CaO/Na2O and low lithophile concentrations. A mixing model involving incremental enrichment of magmaphile elements by repeated episodes of mixing of relatively primitive and moderately evolved magmas, followed by a small amount of fractionation is consistent with the characteristics of the basalts studied.

The Campi Flegrei Deep Drilling Project: using borehole measurements to discriminate magmatic and geothermal effects in caldera unrest

NASA Astrophysics Data System (ADS)

De Natale, Giuseppe; Troise, Claudia; Carlino, Stefano; Troiano, Antonio; Giulia Di Giuseppe, Maria; Piochi, Monica; Somma, Renato; Tramelli, Anna; Kilburn, Christopher

2015-04-01

Large calderas are potentially the most risky volcanic areas in the world since they are capable of producing huge eruptions whose major effects can involve human life and activities from regional to global scale. Calderas worldwide are characterized by frequent episodes of unrest which, only in few cases, culminate with eruptions. This ambiguous behavior is generally explained in terms of magma intrusion or disturbance of geothermal fluids in the shallow crust, which are both source of ground deformations and seismicity. A major goal is to determine the relative contribution of each process, because the potential for eruptions significantly enhanced if magma movements emerge as the primary component. A very important case study is the active Campi Flegrei caldera, hosting part of the large city of Naples (Southern Italy). In the framework of the Campi Flegrei Deep Drilling Project new filed data from pilot borehole have been recorded (permeability and in situ stress) by using a novel procedure of Leak Off Test. These new data, particularly the actual permeability, are fundamental to calibrate the caldera unrest models at Campi Flegrei and, , to put constrains to forecast the maximum future eruptive scenario. We show here that these new data, integrated by fluid-dynamical modeling, allow to assess that only about a third of the maximum uplift recorded in 1982-1984 may be due to shallow aquifer perturbation, so that the remaining part should be due to magma inflow, corresponding to about 0.05 Km3 of new magma if we assume a sill-like reservoir located at 4 km of depth. Considering an almost equivalent magma inflow for the 1969-1972 unrest, which showed a similar uplift, we got a total magma inflow of 0.1 Km3. It is then very important to assess the times for cooling of such accumulated magma, in order to assess the eruption hazard.

Energy Saving Melting and Revert Reduction Technology (Energy-SMARRT): Clean Steel Casting Production

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kuyucak, Selcuk; Li, Delin

2013-12-31

Inclusions in steel castings can cause rework, scrap, poor machining, and reduced casting performance, which can obviously result in excess energy consumption. Significant progress in understanding inclusion source, formation and control has been made. Inclusions can be defined as non-metallic materials such as refractory, sand, slag, or coatings, embedded in a metallic matrix. This research project has focused on the mold filling aspects to examine the effects of pouring methods and gating designs on the steel casting cleanliness through water modeling, computer modeling, and melting/casting experiments. Early in the research project, comprehensive studies of bottom-pouring water modeling and low-alloy steelmore » casting experiments were completed. The extent of air entrainment in bottom-poured large castings was demonstrated by water modeling. Current gating systems are designed to prevent air aspiration. However, air entrainment is equally harmful and no prevention measures are in current practice. In this study, new basin designs included a basin dam, submerged nozzle, and nozzle extension. The entrained air and inclusions from the gating system were significantly reduced using the new basin method. Near the end of the project, there has been close collaboration with Wescast Industries Inc., a company manufacturing automotive exhaust components. Both computer modeling using Magma software and melting/casting experiments on thin wall turbo-housing stainless steel castings were completed in this short period of time. Six gating designs were created, including the current gating on the pattern, non-pressurized, partially pressurized, naturally pressurized, naturally pressurized without filter, and radial choke gating without filter, for Magma modeling. The melt filling velocity and temperature were determined from the modeling. Based on the simulation results, three gating designs were chosen for further melting and casting experiments on the same casting pattern using the lip pouring method. It was observed again that gating designs greatly influenced the melt filling velocity and the number of inclusion defects. The radial choked gating showed improvements in casting cleanliness and yield over the other gatings, even though no mold filters were used in the gating system.« less

Numerical simulation of magma chamber dynamics.

NASA Astrophysics Data System (ADS)

Longo, Antonella; Papale, Paolo; Montagna, Chiara Paola; Vassalli, Melissa; Giudice, Salvatore; Cassioli, Andrea

2010-05-01

Magma chambers are characterized by periodic arrivals of deep magma batches that give origin to complex patterns of magma convection and mixing, and modify the distribution of physical quantities inside the chamber. We simulate the transient, 2D, multi-component homogeneous dynamics in geometrically complex dyke+chamber systems, by means of GALES, a finite element parallel C++ code solving mass, momentum and energy equations for multi-component homogeneous gas-liquid (± crystals) mixtures in compressible-to-incompressible flow conditions. Code validation analysis includes several cases from the classical engineering literature, corresponding to a variety of subsonic to supersonic gas-liquid flow regimes (see http://www.pi.ingv.it/~longo/gales/gales.html). The model allows specification of the composition of the different magmas in the domain, in terms of ten major oxides plus the two volatile species H2O and CO2. Gas-liquid thermodynamics are modeled by using the compositional dependent, non-ideal model in Papale et al. (Chem.. Geol., 2006). Magma properties are defined in terms of local pressure, temperature, and composition including volatiles. Several applications are performed within domains characterized by the presence of one or more magma chambers and one or more dykes, with different geometries and characteristic size from hundreds of m to several km. In most simulations an initial compositional interface is placed at the top of a feeding dyke, or at larger depth, with the deeper magma having a lower density as a consequence of larger volatile content. The numerical results show complex patterns of magma refilling in the chamber, with alternating phases of magma ingression and magma sinking from the chamber into the feeding dyke. Intense mixing takes place in feeding dykes, so that the new magma entering the chamber is always a mixture of the deep and the initially resident magma. Buoyant plume rise occurs through the formation of complex convective patterns, giving origin to a density-stratified magma chamber.

Amplification of seismic waves beneath active volcanoes

NASA Astrophysics Data System (ADS)

Navon, O.; Lensky, N. G.; Collier, L.; Neuberg, J.; Lyakhovsky, V.

2003-04-01

Long-period (LP) seismic events are typical for many volcanoes and are attributed to energy leaking from waves traveling along the conduit - country-rock interface. While the wave propagation is well understood, their actual trigger mechanism and their energy source are not. Here we test the hypothesis that energy may be supplied by volatile-release from a supersaturated melt. If bubbles are initially in equilibrium with the melt in the conduit, and the melt is suddenly decompressed, the transfer of volatiles from the supersaturated melt into the bubbles transforms stored potential energy into expansion work. For example, small dome collapse, opening of a crack or a displacement along the brittle part of the conduit may decompress the magma by a few bars and create the needed supersaturation. This energy is released over the timescale of accelerated expansion, which is longer than a typical LP event. Following decompression, when the transfer of volatiles into bubbles is fast enough, expansion accelerates and the bulk viscosity of the bubbly magma is negative (Lensky et al., 2002). New calculations show that under such conditions a sinusoidal P-wave is amplified. We note that seismic waves created by tectonic earthquakes that are not associated with net decompression, do not lead to net release of volatiles or to net expansion. In this case, the bulk viscosity is positive and waves traveling through the magma should attenuate. The proposed model explains how weak seismic signals may be amplified as they travel through a conduit that contains supersaturated bubbly magma. It provides the general framework for amplifying volcanic seismicity such as long-period events.

Applying Fractal Dimensions and Energy-Budget Analysis to Characterize Fracturing Processes During Magma Migration and Eruption: 2011-2012 El Hierro (Canary Islands) Submarine Eruption

NASA Astrophysics Data System (ADS)

López, Carmen; Martí, Joan; Abella, Rafael; Tarraga, Marta

2014-07-01

The impossibility of observing magma migration inside the crust obliges us to rely on geophysical data and mathematical modelling to interpret precursors and to forecast volcanic eruptions. Of the geophysical signals that may be recorded before and during an eruption, deformation and seismicity are two of the most relevant as they are directly related to its dynamic. The final phase of the unrest episode that preceded the 2011-2012 eruption on El Hierro (Canary Islands) was characterized by local and accelerated deformation and seismic energy release indicating an increasing fracturing and a migration of the magma. Application of time varying fractal analysis to the seismic data and the characterization of the seismicity pattern and the strain and the stress rates allow us to identify different stages in the source mechanism and to infer the geometry of the path used by the magma and associated fluids to reach the Earth's surface. The results obtained illustrate the relevance of such studies to understanding volcanic unrest and the causes that govern the initiation of volcanic eruptions.

Superhot fluids circulating close to magma intrusions: a contribution from analogue modelling

NASA Astrophysics Data System (ADS)

Montanari, Domenico; Agostini, Andrea; Bonini, Marco; Corti, Giacomo

2017-04-01

Magma overpressure at the time of the emplacement at shallow crustal levels may lead to deformation (i.e. forced folding, fracturing and faulting) in the country rock, both at local and regional scale. To get insights into this process, we reproduced and analysed in the laboratory the fracture/fault network associated with the emplacement of magma at shallow crustal levels. We used a mixture of quartz sand and K-feldspar fine sand as an analogue for the brittle crust, and polyglycerols for the magma. The models were able to reproduce complex 3D architectures of deformation resulting from magma emplacement, with different deformation patterns -invariably dominated by forced folding and associated brittle faulting/fracturing- resulting from variable parameters. These results provide useful hints into geothermal researches. Fractures and faults associated with magma emplacement are indeed expected to significantly influence the distribution and migration of superhot geothermal fluids near the edge of the magma intrusion. These structures can therefore be considered as potential targets for geothermal or mineral deposits exploration. In this perspective, the results of analogue models may provide useful geometric and conceptual constraints for field work, numerical modeling, and particularly seismic interpretation for achieving a better understanding and tuning of the integrated conceptual model concerning the circulation of supercritical fluids. The research leading to these results has received funding from the European Community's Seventh Framework Programme under grant agreement No. 608553 (Project IMAGE).

The trigger mechanism of low-frequency earthquakes on Montserrat

NASA Astrophysics Data System (ADS)

Neuberg, J. W.; Tuffen, H.; Collier, L.; Green, D.; Powell, T.; Dingwell, D.

2006-05-01

A careful analysis of low-frequency seismic events on Soufrièere Hills volcano, Montserrat, points to a source mechanism that is non-destructive, repetitive, and has a stationary source location. By combining these seismological clues with new field evidence and numerical magma flow modelling, we propose a seismic trigger model which is based on brittle failure of magma in the glass transition. Loss of heat and gas from the magma results in a strong viscosity gradient across a dyke or conduit. This leads to a build-up of shear stress near the conduit wall where magma can rupture in a brittle manner, as field evidence from a rhyolitic dyke demonstrates. This brittle failure provides seismic energy, the majority of which is trapped in the conduit or dyke forming the low-frequency coda of the observed seismic signal. The trigger source location marks the transition from ductile conduit flow to friction-controlled magma ascent. As the trigger mechanism is governed by the depth-dependent magma parameters, the source location remains fixed at a depth where the conditions allow brittle failure. This is reflected in the fixed seismic source locations.

Lifetime and size of shallow magma bodies controlled by crustal-scale magmatism

NASA Astrophysics Data System (ADS)

Karakas, Ozge; Degruyter, Wim; Bachmann, Olivier; Dufek, Josef

2017-06-01

Magmatic processes on Earth govern the mass, energy and chemical transfer between the mantle, crust and atmosphere. To understand magma storage conditions in the crust that ultimately control volcanic activity and growth of continents, an evaluation of the mass and heat budget of the entire crustal column during magmatic episodes is essential. Here we use a numerical model to constrain the physical conditions under which both lower and upper crustal magma bodies form. We find that over long durations of intrusions (greater than 105 to 106 yr), extensive lower crustal mush zones develop, which modify the thermal budget of the upper crust and reduce the flux of magma required to sustain upper crustal magma reservoirs. Our results reconcile physical models of magma reservoir construction and field-based estimates of intrusion rates in numerous volcanic and plutonic localities. Young igneous provinces (less than a few hundred thousand years old) are unlikely to support large upper crustal reservoirs, whereas longer-lived systems (active for longer than 1 million years) can accumulate magma and build reservoirs capable of producing super-eruptions, even with intrusion rates smaller than 10-3 to 10-2 km3 yr-1. Hence, total duration of magmatism should be combined with the magma intrusion rates to assess the capability of volcanic systems to form the largest explosive eruptions on Earth.

Geomorphic Mapping of Lava Flows on Mars, Earth, and Mercury

NASA Astrophysics Data System (ADS)

Golder, K. B.; Burr, D. M.

2018-06-01

To advance understanding of flood basalts, we have mapped lava flows on three planets, Mars, Earth, and Mercury, as part of three projects. The common purpose of each project is to investigate potential magma sources and/or emplacement conditions.

Experimental Study of Lunar and SNC Magmas

NASA Astrophysics Data System (ADS)

Rutherford, Malcolm J.

2000-08-01

The research described in this progress report involved the study of petrological, geochemical and volcanic processes that occur on the Moon and the SNC parent body, generally accepted to be Mars. The link between these studies is that they focus on two terrestrial-type parent bodies somewhat smaller than earth, and the fact that they focus on the role of volatiles in magmatic processes and on processes of magma evolution on these planets. The work on the lunar volcanic glasses has resulted in some exciting new discoveries over the years of this grant. During the tenure of the present grant, we discovered a variety of metal blebs in the A17 orange glass. Some of these Fe-Ni metal blebs occur in the glass; others were found in olivine phenocrysts which we find make up about 2 vol % of the orange glass magma. The importance of these metal spheres is that they fix the oxidation state of the parent magma during the eruption, and also indicate changes during the eruption. They also yield important information about the composition of the gas phase present, the gas which drove the lunar fire-fountaining. In an Undergraduate senior thesis project, Nora Klein discovered a melt inclusion that remained in a glassy state in one of the olivine phenocrysts. Analyses of this inclusion gave additional information on the CO2, CO and S contents of the orange glass magma prior to its reaching the lunar surface. The composition of lunar volcanic gases has long been one of the puzzles of lunar magmatic processes. One of the more exciting findings in our research over the past year has been the study of magmatic processes linking the SNC meteorite source magma composition with the andesitic composition rocks found at the Pathfinder site. In this project, graduate student Michelle Minitti showed that there was a clear petrologic link between these two magma types via fractional removal of crystals from the SNC parent melt, but the process only worked if there was at least 1 wt % dissolved water in the melt.

Experimental Study of Lunar and SNC Magmas

NASA Technical Reports Server (NTRS)

Rutherford, Malcolm J.

2000-01-01

The research described in this progress report involved the study of petrological, geochemical and volcanic processes that occur on the Moon and the SNC parent body, generally accepted to be Mars. The link between these studies is that they focus on two terrestrial-type parent bodies somewhat smaller than earth, and the fact that they focus on the role of volatiles in magmatic processes and on processes of magma evolution on these planets. The work on the lunar volcanic glasses has resulted in some exciting new discoveries over the years of this grant. During the tenure of the present grant, we discovered a variety of metal blebs in the A17 orange glass. Some of these Fe-Ni metal blebs occur in the glass; others were found in olivine phenocrysts which we find make up about 2 vol % of the orange glass magma. The importance of these metal spheres is that they fix the oxidation state of the parent magma during the eruption, and also indicate changes during the eruption. They also yield important information about the composition of the gas phase present, the gas which drove the lunar fire-fountaining. In an Undergraduate senior thesis project, Nora Klein discovered a melt inclusion that remained in a glassy state in one of the olivine phenocrysts. Analyses of this inclusion gave additional information on the CO2, CO and S contents of the orange glass magma prior to its reaching the lunar surface. The composition of lunar volcanic gases has long been one of the puzzles of lunar magmatic processes. One of the more exciting findings in our research over the past year has been the study of magmatic processes linking the SNC meteorite source magma composition with the andesitic composition rocks found at the Pathfinder site. In this project, graduate student Michelle Minitti showed that there was a clear petrologic link between these two magma types via fractional removal of crystals from the SNC parent melt, but the process only worked if there was at least 1 wt % dissolved water in the melt.

Crustal-scale magmatism and its control on the longevity of magmatic systems

NASA Astrophysics Data System (ADS)

Karakas, Ozge; Degruyter, Wim; Bachmann, Olivier; Dufek, Josef

2017-04-01

Constraining the duration and evolution of crustal magma reservoirs is crucial to our understanding of the eruptive potential of magmatic systems, as well as the volcanic:plutonic ratios in the crust, but estimates of such parameters vary widely in the current literature. Although no consensus has been reached on the lifetime of magma reservoirs, recent studies have revealed about the presence, location, and melt fraction of multi-level (polybaric) storage zones in the crust. If magma accumulates at different crustal levels, it must redistribute significant enthalpy within the crustal column and therefore must influence the lifetime of magma plumbing systems. However, an evaluation of the mass and heat budget of the entire crustal column is lacking. Here, we use a two-dimensional thermal model to determine the thermal conditions under which both lower and upper crustal magma bodies form. We find that large lower crustal mush zones supply heat to the upper crust and reduce the amount of thermal energy necessary to form subvolcanic reservoirs. This indicates that the crust is thermally viable to sustain partially molten magma reservoirs over long timescales (>10^5-106 yr) for a range of magma fluxes (10^-4 to 10^-2 km^3/yr). Our results reconcile physical models of crustal magma evolution and field-based estimates of intrusion rates in numerous magmatic provinces (which include both volcanic and plutonic lithologies). We also show that young magmatic provinces (< 105 yr old) are unlikely to support large upper crustal reservoirs, whereas longer-lived systems (> 106 yr) can accumulate magma and build reservoirs capable of triggering supereruptions, even with intrusion rates as low as ≤10^-2 km^3/yr. Hence, the total duration of magmatism is critical in determining the size of the magma reservoirs, and should be combined with the magma intrusions rates to assess the capability of volcanic systems to form the largest eruptions on Earth.

Timescales and conditions of crystallization in the Pokai and Chimpanzee Ignimbrites, Taupo Volcanic Zone, New Zealand

NASA Astrophysics Data System (ADS)

Connor, M.; Gualda, G. A.; Gravley, D. M.

2013-12-01

Silicic magmas give rise to explosive eruptions that are both of scientific and societal interest. The central Taupo Volcanic Zone in New Zealand has been volcanically active for 2 Ma and represents the most active volcanic area in the world today. Particularly intense volcanic activity took place as part of a flare-up event that occurred from ~340 to ~240 ka, when 7 large ignimbrite eruptions took place, as well as many smaller eruptions, which erupted a total of at least 3000 km3 of magma. This project seeks to identify the conditions and timescales over which magma bodies that gave rise to these ignimbrite eruptions evolved. We aim to understand how much of the tens of thousands of years between successive eruptions were characterized by the presence of large bodies of silicic magma within the crust, as well the magma distribution within the crust during those times. We focus on the Chimpanzee and Pokai ignimbrites, which together erupted ~150 km3 of magma. The Pokai ignimbrite erupted at ~275 ka, while the Chimpanzee ignimbrite (undated) erupted between ~320 and 275 ka. Pumice clasts from the Chimpanzee and Pokai ignimbrite were collected in the field. Pumice bulk densities were measured using a submersion technique. Quartz and plagioclase crystals were extracted through a crushing, sieving, and winnowing procedure. Whole crystals were hand-picked under a conventional microscope, mounted on epoxy, and polished to expose grain interiors. Grain mounts were analyzed under an SEM using back-scattered electron, cathodoluminescence (CL), and energy-dispersive x-ray (EDX) imaging. Bulk-densities vary from 0.42 to 0.81 g/cm3 for Pokai and between 0.52 and 0.64 g/cm3 for Chimpanzee pumice clasts. Plagioclase is the dominant crystal phase in both units. Several plagioclase crystals have inclusions of orthopyroxene, ilmenite, magnetite, and zircon, which in some cases form clusters. Quartz is rare but is present in pumice from both deposits. Both plagioclase and quartz have glass inclusions. Quantitative EDX analyses under the SEM reveal high-silica rhyolite compositions (~78 wt. % SiO2), consistent with crystallization at shallow crustal levels (~100 MPa). Rhyolite-MELTS barometry confirms this suggestion, with estimated crystallization pressures of 70-110 MPa. Such crystallization depths are characteristic of magma bodies that give rise to ignimbrites in the Taupo Volcanic Zone. CL images of a small number of quartz crystals reveal zoning that is variable in character, from unzoned to distinctly zoned. Zoned crystals have several discrete zones, each one with its own characteristic homogeneous brightness; bright CL rims are observed in these crystals. Contacts between zones can be used to calculate crystallization times using diffusional relaxation assuming that CL intensity variations are due to Ti zonation. We estimate that quartz crystallization took place over decades to a few centuries, with rims crystallizing over months to years at most. These short timescales suggest that magma bodies with eruptible magma are very short-lived, much shorter-lived than the 10s of ka repose periods between eruptions that make up the flare-up event.

Petrology and Geochemistry of New Paired Martian Meteorites Larkman Nunatak 12240 and Larkman Nunatak 12095

NASA Technical Reports Server (NTRS)

Funk, R. C.; Peslier, A. H.; Brandon, A. D.; Humayun, M.

2016-01-01

Two of the latest Martian meteorites found in Antarctica, paired olivine-phyric shergottites LAR 12240 and LAR 12095, are described in order to decipher their petrological context, and place constraints on the geological history of Mars. This project identifies all phases found in LAR 12240 and 12095 and analyzes them for major and trace elements. The textural relationships among these phases are examined in order to develop a crystallization history of the magma(s) that formed these basalts.

The Krafla International Testbed (KMT): Ground Truth for the New Magma Geophysics

NASA Astrophysics Data System (ADS)

Brown, L. D.; Kim, D.; Malin, P. E.; Eichelberger, J. C.

2017-12-01

Recent developments in geophysics such as large N seismic arrays , 4D (time lapse) subsurface imaging and joint inversion algorithms represent fresh approaches to delineating and monitoring magma in the subsurface. Drilling at Krafla, both past and proposed, are unique opportunities to quantitatively corroborate and calibrate these new technologies. For example, dense seismic arrays are capable of passive imaging of magma systems with resolutions comparable to that achieved by more expensive (and often logistically impractical) controlled source surveys such as those used in oil exploration. Fine details of the geometry of magma lenses, feeders and associated fluid bearing fracture systems on the scale of meters to tens of meters are now realistic targets for surface seismic surveys using ambient energy sources, as are detection of their temporal variations. Joint inversions, for example of seismic and MT measurements, offer the promise of tighter quantitative constraints on the physical properties of the various components of magma and related geothermal systems imaged by geophysics. However, the accuracy of such techniques will remain captive to academic debate without testing against real world targets that have been directly sampled. Thus application of these new techniques to both guide future drilling at Krafla and to be calibrated against the resulting borehole observations of magma are an important step forward in validating geophysics for magma studies in general.

Crystalline heterogeneities and instabilities in thermally convecting magma chamber

NASA Astrophysics Data System (ADS)

Culha, C.; Suckale, J.; Qin, Z.

2016-12-01

A volcanic vent can supply different densities of crystals over an eruption time period. This has been seen in Hawai'i's Kilauea Iki 1959 eruption; however it is not common for all Kilauea or basaltic eruptions. We ask the question: Under what conditions can homogenous magma chamber cultivate crystalline heterogeneities? In some laboratory experiments and numerical simulations, a horizontal variation is observed. The region where crystals reside is identified as a retention zone: convection velocity balances settling velocity. Simulations and experiments that observe retention zones assume crystals do not alter the convection in the fluid. However, a comparison of experiments and simulations of convecting magma with crystals suggest that large crystal volume densities and crystal sizes alter fluid flow considerably. We introduce a computational method that fully resolves the crystalline phase. To simulate basaltic magma chambers in thermal convection, we built a numerical solver of the Navier-Stoke's equation, continuity equation, and energy equation. The modeled magma is assumed to be a viscous, incompressible fluid with a liquid and solid phase. Crystals are spherical, rigid bodies. We create Rayleigh-Taylor instability through a cool top layer and hot bottom layer and update magma density while keeping crystal temperature and size constant. Our method provides a detailed picture of magma chambers, which we compare to other models and experiments to identify when and how crystals alter magma chamber convection. Alterations include stratification, differential settling and instabilities. These characteristics are dependent on viscosity, convection vigor, crystal volume density and crystal characteristics. We reveal that a volumetric crystal density variation may occur over an eruption time period, if right conditions are met to form stratifications and instabilities in magma chambers. These conditions are realistic for Kilauea Iki's 1959 eruption.

The variation of magma discharge during basaltic eruptions

NASA Technical Reports Server (NTRS)

Wadge, G.

1981-01-01

The rate at which basaltic magma is discharged during many eruptions varies substantially. An individual eruption has an eruption rate, which is the volumetric rate of discharge averaged over the whole or a major part of an eruption, and an effusion rate, which is the volumetric flux rate at any given time. In many cases, the effusion rate soon reaches a maximum after a short period of waxing flow (partly because of magmatic expansion); it then falls more slowly in the later parts of the eruption. The release of elastic strain energy from stored magma and the subvolcanic reservoir during eruption can give a waning flow of this type an exponential form. A comparison of the eruption rates of eruptions of Mauna Loa, Kilauea and Etna shows that for each volcano there is a trend of decreasing effusion rate with increasing duration of eruption. It is noted that this relationship is not predicted by a simple elastic model of magma release. Two other processes are invoked to explain the eruptive histories of these volcanoes: modification of the eruptive conduits and the continued supply of magma from depth during eruption.

Episodic inflation events at Akutan Volcano, Alaska, during 2005-2017

NASA Astrophysics Data System (ADS)

Ji, Kang Hyeun; Yun, Sang-Ho; Rim, Hyoungrea

2017-08-01

Detection of weak volcano deformation helps constrain characteristics of eruption cycles. We have developed a signal detection technique, called the Targeted Projection Operator (TPO), to monitor surface deformation with Global Positioning System (GPS) data. We have applied the TPO to GPS data collected at Akutan Volcano from June 2005 to March 2017 and detected four inflation events that occurred in 2008, 2011, 2014, and 2016 with inflation rates of about 8-22 mm/yr above the background trend at a near-source site AV13. Numerical modeling suggests that the events should be driven by closely located sources or a single source in a shallow magma chamber at a depth of about 4 km. The inflation events suggest that magma has episodically accumulated in a shallow magma chamber.

Assimilation by lunar mare basalts: Melting of crustal material and dissolution of anorthite

NASA Astrophysics Data System (ADS)

Finnila, A. B.; Hess, P. C.; Rutherford, M. J.

1994-07-01

We discuss techniques for calculating the amount of crustal assimilation possible in lunar magma chambers and dikes based on thermal energy balances, kinetic rates, and simple fluid mechanical constraints. Assuming parent magmas of picritic compositions, we demonstrate the limits on the capacity of such magmas to melt and dissolve wall rock of anorthitic, troctolitic, noritic, and KREEP (quartz monzodiorite) compositions. Significant melting of the plagioclase-rich crustal lithologies requires turbulent convection in the assimilating magma and an efficient method of mixing in the relatively buoyant and viscous new melt. Even when this occurs, the major element chemistry of the picritic magmas will change by less than 1-2 wt %. Diffusion coefficients measured for Al2O3 from an iron-free basalt and an orange glass composition are 10-12 sq m/s at 1340 C and 10-11 sq m/s at 1390 C. These rates are too slow to allow dissolution of plagioclase to significantly affect magma compositions. Picritic magmas can melt significant quantities of KREEP, which suggests that their trace element chemistry may still be affected by assimilation processes; however, mixing viscous melts of KREEP composition with the fluid picritic magmas could be prohibitively difficult. We conclude that only a small part of the total major element chemical variation in the mare basalt and volcanic glass collection is due to assimilation/fractional crystallization processes near the lunar surface. Instead, most of the chemical variation in the lunar basalts and volcanic glasses must result from assimilation at deeper levels or from having distinct source regions in a heterogeneous lunar mantle.

Assimilation by lunar mare basalts: Melting of crustal material and dissolution of anorthite

NASA Technical Reports Server (NTRS)

Finnila, A. B.; Hess, P. C.; Rutherford, M. J.

1994-01-01

We discuss techniques for calculating the amount of crustal assimilation possible in lunar magma chambers and dikes based on thermal energy balances, kinetic rates, and simple fluid mechanical constraints. Assuming parent magmas of picritic compositions, we demonstrate the limits on the capacity of such magmas to melt and dissolve wall rock of anorthitic, troctolitic, noritic, and KREEP (quartz monzodiorite) compositions. Significant melting of the plagioclase-rich crustal lithologies requires turbulent convection in the assimilating magma and an efficient method of mixing in the relatively buoyant and viscous new melt. Even when this occurs, the major element chemistry of the picritic magmas will change by less than 1-2 wt %. Diffusion coefficients measured for Al2O3 from an iron-free basalt and an orange glass composition are 10(exp -12) sq m/s at 1340 C and 10(exp -11) sq m/s at 1390 C. These rates are too slow to allow dissolution of plagioclase to significantly affect magma compositions. Picritic magmas can melt significant quantities of KREEP, which suggests that their trace element chemistry may still be affected by assimilation processes; however, mixing viscous melts of KREEP composition with the fluid picritic magmas could be prohibitively difficult. We conclude that only a small part of the total major element chemical variation in the mare basalt and volcanic glass collection is due to assimilation/fractional crystallization processes near the lunar surface. Instead, most of the chemical variation in the lunar basalts and volcanic glasses must result from assimilation at deeper levels or from having distinct source regions in a heterogeneous lunar mantle.

Reporting from the Iceland Deep Drilling Project

NASA Astrophysics Data System (ADS)

Urban, Karl

2017-04-01

Geoscience-related topics are in many cases difficult to communicate to the public: Often they include dead soil which not easily tells lively stories. And it is hard to sell those topics to editors of public media. In addition the topics might also be politically supercharged if they are resource-related with a visible environmental impact. Therefore any researcher involved might be overcautious while talking to journalists. With a grant from the EGU Science Journalist Fellowship I travelled to Iceland in autumn 2016 to report about the Iceland Deep Drilling Project (IDDP). The project which started just weeks prior to my arrival aimed to drill the deepest borehole in a volcanically active region. During earlier trials the borehole collapsed or the drill string unintentionally hit magma. If successful the IDDP promises a much higher level of geothermal energy harvested. The IDDP was therefore ideally suited to be sold to public media outlets since Iceland's volcanic legacy easily tells a lively story. But the drilling's potential environmental impact makes it a political topic in Iceland - even though geothermal energy has a positive public perception. Therefore the IDDP included some pitfalls I observed several times before while reporting about geoscience research. Those could be circumvented if researchers and journalists knew better about their expectations before any interview takes place.

Thermal evolution of plutons: a parameterized approach.

PubMed

Spera, F

1980-01-18

A conservation-of-energy equation has been derived for the spatially averaged magma temperature in a spherical pluton undergoing simultaneous crystallization and both internal (magma) and external (hydrothermal fluid) thermal convection. The model accounts for the dependence of magma viscosity on crystallinity, temperature, and bulk composition; it includes latent heat effects and the effects of different initial water concentrations in the melt and quantitatively considers the role that large volumes of circulatory hydrothermal fluids play in dissipating heat. The nonlinear ordinary differential equation describing these processes has been solved for a variety of magma compositions, initial termperatures, initial crystallinities, volume ratios of hydrothermal fluid to magma, and pluton sizes. These calculations are graphically summarized in plots of the average magma temperature versus time after emplacement. Solidification times, defined as the time necessary for magma to cool from the initial emplacement temperature to the solidus temperature vary as R(1,3), where R is the pluton radius. The solidification time of a pluton with a radius of 1 kilometer is 5 x 10(4) years; for an otherwise identical pluton with a radius of 10 kilometers, the solidification time is approximately 10(6) years. The water content has a marked effect on the solidification time. A granodiorite pluton with a radius of 5 kilometers and either 0.5 or 4 percent (by weight) water cools in 3.3 x 10(5) or 5 x 10(4) years, respectively. Convection solidification times are usually but not always less than conduction cooling times.

Experimental study of lunar and SNC (Mars) magmas

NASA Technical Reports Server (NTRS)

Rutherford, Malcolm J.

1994-01-01

The overall objectives of this research were to evaluate the role of C-O-S-Cl degassing processes in explaining vesiculation, oxidation state and fire-fountaining of lunar magmas by analysis of individual lunar glass spherules, and by experimental determination of equilibrium abundances and diffusion rates of C, S and Cl melt species in lunar glass compositions; and to determine possible primitive SNC magma compositions and the mineralogy of the mantle from which they were derived, and to evaluate P, T, XH2O etc. conditions at which they crystallize to form the SNC meteorites. After funding for one year, a project on the A15 volcanic green glass has been completed to the point of writing a first manuscript. Carbon-oxygen species C-O and CO2 are below detection limits (20 ppm) in these glasses, but there is up to 500 ppm S with concentrations both increasing and decreasing toward the spherule margins. Calculations and modeling indicate that C species could have been present in the volcanic gases, however. In a second project, experiments with low PH2O have resulted in refined estimates of the early intercumulus melt composition in the Chassigny meteorite which is generally accepted as a sample from Mars.

Chemical consequences of compaction within the freezing front of a crystallizing magma ocean

NASA Astrophysics Data System (ADS)

Hier-Majumder, S.; Hirschmann, M. M.

2013-12-01

The thermal and compositional evolution of planetary magma oceans have profound influences on the early development and differentiation of terrestrial planets. During crystallization, rejection of elements incompatible in precipitating solids leads to petrologic and geochemical planetary differentiation, including potentially development of a compositionally stratified early mantle and evolution of thick overlying atmospheres. In cases of extremely efficient segregation of melt and crystals, solidified early mantles can be nearly devoid of key incompatible species including heat-producing (U, Th, K) and volatile (H,C,N,& noble gas) elements. A key structural component of a crystallizing magma ocean is the partially molten freezing front. The dynamics of this region influences the distribution of incompatible elements between the earliest mantle and the initial surficial reservoirs. It also can be the locus of heating owing to the dissipation of large amounts of tidal energy potentially available from the early Moon. The dynamics are influenced by the solidification rate, which is coupled to the liberation of volatiles owing to the modulating greenhouse effects in the overlying thick atmosphere. Compaction and melt retention in the freezing front of a magma ocean has received little previous attention. While the front advances during the course of crystallization, coupled conservation of mass, momentum, and energy within the front controls distribution and retention of melt within this layer. Due to compaction within this layer, melt distribution is far from uniform, and the fraction of melt trapped within this front depends on the rate of freezing of the magma ocean. During phases of rapid freezing, high amount of trapped melt within the freezing front retains a larger quantity of dissolved volatiles and the reverse is true during slow periods of crystallization. Similar effects are known from inferred trapped liquid fractions in layered mafic intrusions. Here we develop a simple 1-D model of melt retention in the freezing front of a crystallizing magma ocean, and apply it to the thermal and chemical evolution of the early Earth.

Assimilation by Lunar Mare Basalts: Melting of Crustal Material and Dissolution of Anorthite

NASA Technical Reports Server (NTRS)

Finnila, A. B.; Hess, P. C.; Rutherford, M. J.

1994-01-01

We discuss techniques for calculating the amount of crustal assimilation possible in lunar magma chambers and dikes based on thermal energy balances, kinetic rates, and simple fluid mechanical constraints. Assuming parent magmas of picritic compositions, we demonstrate the limits on the capacity of such magmas to melt and dissolve wall rock of anorthitic, troctolitic, noritic, and KREEP (quartz monzodiorite) compositions. Significant melting of the plagioclase-rich crustal lithologies requires turbulent convection in the assimilating magma and an efficient method of mixing in the relatively buoyant and viscous new melt. Even when this occurs, the major element chemistry of the picritic magmas will change by less than 1-2 wt %. Diffusion coefficients measured for Al2O3 from an iron-free basalt and an orange glass composition are 10(exp -12) m(exp 2) s(exp -1) at 1340 C and 10(exp -11) m(exp 2) s(exp -1) at 1390 C. These rates are too slow to allow dissolution of plagioclase to significantly affect magma compositions. Picritic magmas can melt significant quantities of KREEP, which suggests that their trace element chemistry may still be affected by assimilation processes; however, mixing viscous melts of KREEP composition with the fluid picritic magmas could be prohibitively difficult. We conclude that only a small part of the total major element chemical variation in the mare basalt and volcanic glass collection is due to assimilation/fractional crystallization processes near the lunar surface. Instead, most of the chemical variation in the lunar basalts and volcanic glasses must result from assimilation at deeper levels or from having distinct source regions in a heterogeneous lunar mantle.

Creep, dike intrusion, and magma chamber deflation model for the 2000 Miyake eruption and the Izu islands earthquakes

NASA Astrophysics Data System (ADS)

Ozawa, S.; Miyazaki, S.; Nishimura, T.; Murakami, M.; Kaidzu, M.; Imakiire, T.; Ji, X.

2004-02-01

Analysis of Global Positioning System data shows shrinkage of Miyake Island and the widening between Nijima and Kozu Islands during the period of the Miyake island volcanic activity and the ensuing Izu islands earthquakes in 2000. The estimated time evolution of a model consisting of a dike, creeping faults, and a Mogi source suggests that a crack opening on Miyake Island occurred immediately after the start of the seismic activities on 26 June and ended within several days at the west coast of Miyake Island. After the seabed eruption on 27 June, magma migrated from Miyake Island to Kozu and Nijima Islands within several days. The estimated volume of intruded magma totals around 1.2 × 109 m3. Associated with the magma intrusion between Miyake and Kozu Islands, left-lateral and right-lateral creep motions occurred in regions off the west coast of Miyake Island and near Kozu Island. The accumulated moment energy is equivalent to an earthquake of Mw 6.6 and 6.6 for right-lateral and left-lateral creeping faults, respectively. The estimated magma chamber continued deflation beneath the southwestern part of Miyake Island from 26 June, totaling around 0.12 × 109 m3 in volume change, in addition to the collapse volume of 0.6 × 109 m3 at the summit of Mount Oyama on Miyake Island. The volume change on Miyake Island can be compensated by the migrated magma toward Kozu Island from the deflation source beneath Miyake Island. The deflation speed of the magma chamber beneath Miyake Island decreased and increased before and after the eruption of 14-15 July and 18 August, suggesting a change in balance of mass influx and draining out rate of the magma chamber.

Exploration of Geothermal Natural Resources from Menengai Caldera at Naruku, Kenya

NASA Astrophysics Data System (ADS)

Patlan, E.; Wamalwa, A.; Thompson, L. E.; Kaip, G.; Velasco, A. A.

2011-12-01

The Menengai Caldera, a large, dormant volcano, lies near the city of Naruku, Kenya (0.20°S, 36.07°E) and presents a significant natural geothermal energy resource that will benefit local communities. Kenya continues to explore and exploit its only major energy resource: geothermal energy. The Geothermal Development Company (GDC) of Kenya and University of Texas at El Paso (UTEP) have initially deployed seven seismic stations to address the volcanic hazards and associated processes that occurs through the analysis of data collection from seismic sensors that record ground motion. Seven more sensors are planned to be deployed in Aug. 2011. In general, the internal state and activity of the caldera is an important component to the understanding of porosity of the fault system, which is derived from the magma movement of the hot spot, and for the exploitation of geothermal energy. We analyze data from March to May 2011 to investigate the role of earthquakes and faults in controlling the caldera processes, and we find 15 events occurred within the caldera. We will utilize the double difference earthquake location algorithm (HypoDD) to analyze the local events in order to find active faulting of the caldera and the possible location of the magma chamber. For future work, we will combine the exiting data with the new seismic station to image the location of the caldera magma chamber.

European collaboration for improved monitoring of Icelandic volcanoes: Status of the FUTUREVOLC project after the initial 18 months

NASA Astrophysics Data System (ADS)

Dumont, Stéphanie; Parks, Michelle; Sigmundsson, Freysteinn; Vogfjörð, Kristín; Einarsdóttir, Heiðveig Maria; Tumi Gudmundsson, Magnús; Kristinsson, Ingvar; Loughlin, Sue; Ilyinskaya, Evgenia; Hooper, Andrew; Kylling, Arve; Witham, Claire; Bean, Chris; Braiden, Aoife; Ripepe, Maurizio; Prata, Fred; Pétur Heiðarsson, Einar; Other Members Of The Futurevolc Team

2014-05-01

The FUTUREVOLC project funded by the European Union (FP7) is devoted to volcanic hazard assessment and establishing an integrated volcanological monitoring procedure through a European collaboration. To reach these objectives the project combines broad expertise from 26 partners from 10 countries, focusing on the four most active volcanoes of Iceland: Grímsvötn, Katla, Hekla and Bárdarbunga. The geological setting of Iceland, the high rate of eruptions and the various eruption styles make this country an optimal natural laboratory to study volcanic processes from crustal depths to the atmosphere. The project, which began on 1 October 2012, integrates advanced monitoring and analytical techniques in an innovative way, focusing on (i) detailed monitoring to improve our understanding of the seismic/magmatic unrest, in order to estimate the amount of magma available for an eruption and to provide early warnings (ii) the dynamics of magma in the conduit and a near real time estimation of the mass eruption rate and (iii) observing and modelling the plume dynamics. The project design considers effective collaboration between partners and aims for efficient cross-disciplinary workflows. A major step during the first 18 months of the project was the installation of additional equipment in the volcanic regions of Iceland to reinforce and complement the existing monitoring. The instruments include: seismometers, GPS stations, MultigGAS detectors, DOAS, infrasonic arrays, electric field sensors, radars, and optical particle sizers. Data streaming is designed to withstand extreme weather conditions. The FUTUREVOLC project has an open data policy for real and near-time data. Implementation of a data hub is currently under way, based on open access to data from the 2010 Eyjafjallajökull eruption. Access to volcano monitoring data through a common interface will allow timely information on magma movements facilitated through combined analysis. A key part of the project is to combine ground based measurements with satellite observations for both monitoring ground deformation (e.g. GPS, tilt and InSAR) and quantifying the plume parameters (using ground and space based IR sensors). An important complementary effort is the establishment of an "Icelandic Volcanoes" permanent Geohazard Supersite by the Committee on Earth Observation Satellites (CEOS). The support by CEOS plays a vital role in providing greater access to past and future SAR acquisitions over the coming years, therefore promoting a better understanding of magma migration over long time periods. New methods will be developed within FUTUREVOLC, taking into consideration all aspects of a volcanic crisis - not only to improve our understanding of magmatic processes, but also for the delivery of relevant information to civil protection and local authorities. A report on the lessons learned from Eyjafjallajökull 2010 is in preparation in this regard.

Magma Differentiation and Storage Inferred from Crystal Textures at Harrat Rahat Volcanic Field, Kingdom of Saudi Arabia

NASA Astrophysics Data System (ADS)

Witter, M. R.; Mahood, G. A.; Stelten, M. E.; Downs, D. T.; Zahran, H. M.

2015-12-01

We present results of a petrographic study of Harrat Rahat volcanic field in western Saudi Arabia as part of a collaborative project between the U.S.G.S. and the Saudi Geological Survey. Lavas range in composition from alkali basalt to trachyphonolite. Basalts have <2-10 vol.% phenocrysts of euhedral olivine and plagioclase (± minor clinopyroxene). In intermediate lavas, phenocrysts (<5 vol.%) of olivine and plagioclase are resorbed, and plagioclase also exhibits sieve textures and strong zoning, indicative of complex magmatic histories. Trachyphonolite lavas have 0-35 vol.% large phenocrysts of anorthoclase and trace fayalitic olivine but are characterized by a size distribution of crystals that is seriate in hand specimen, so that most exceeded 45% crystals at the time of eruption. Some contain groundmass alkali amphibole. Crystal size distributions (CSD) of crystal-rich trachyphonolites produce simple linear trends (see below), which are interpreted as signifying that all the crystals are related through a common nucleation and growth history, at more or less constant pressure. Linear CSDs indicate no loss of small crystals due to reheating of magmas by recharge, no gain of small crystals due to late-stage nucleation on ascent or degassing, and no addition of large phenocrysts by crystal accumulation or magma mixing. Experimental studies demonstrate that silica-undersaturated evolved magmas like those erupted at Harrat Rahat can form by fractionation of alkali basalts at crustal depths greater than ~25 km. The observed phenocryst assemblage in the trachyphonolites, however, forms at shallow depths, ~2-4 km, according to MELTS modeling. Coupled with CSD data, this suggests that deep extraction events yield crystal-poor trachyphonolite magmas that rise to the upper crust where they undergo crystallization. Extensive shallow crystallization of trachyphonolites may have triggered eruptions by causing vapor saturation, which lowers magma density via vesiculation and has the potential to explosively disrupt wallrocks. Based on the complex crystal textures in the intermediate lavas and their similarity in age to the trachyphonolites, ~120 ka, we suggest that most of the intermediate magmas form by magma mixing when rising basalts intercept and entrain shallow trachyphonolite magma.

Volcanic Centers in the East Africa Rift: Volcanic Processes with Seismic Stresses to Identify Potential Hydrothermal Vents

NASA Astrophysics Data System (ADS)

Patlan, E.; Wamalwa, A. M.; Kaip, G.; Velasco, A. A.

2015-12-01

The Geothermal Development Company (GDC) in Kenya actively seeks to produce geothermal energy, which lies within the East African Rift System (EARS). The EARS, an active continental rift zone, appears to be a developing tectonic plate boundary and thus, has a number of active as well as dormant volcanoes throughout its extent. These volcanic centers can be used as potential sources for geothermal energy. The University of Texas at El Paso (UTEP) and the GDC deployed seismic sensors to monitor several volcanic centers: Menengai, Silali, and Paka, and Korosi. We identify microseismic, local events, and tilt like events using automatic detection algorithms and manual review to identify potential local earthquakes within our seismic network. We then perform the double-difference location method of local magnitude less than two to image the boundary of the magma chamber and the conduit feeding the volcanoes. In the process of locating local seismicity, we also identify long-period, explosion, and tremor signals that we interpret as magma passing through conduits of the magma chamber and/or fluid being transported as a function of magma movement or hydrothermal activity. We used waveform inversion and S-wave shear wave splitting to approximate the orientation of the local stresses from the vent or fissure-like conduit of the volcano. The microseismic events and long period events will help us interpret the activity of the volcanoes. Our goal is to investigate basement structures beneath the volcanoes and identify the extent of magmatic modifications of the crust. Overall, these seismic techniques will help us understand magma movement and volcanic processes in the region.

Physical mechanisms leading to two-dimensional gas content evolution within a volcanic conduit

NASA Astrophysics Data System (ADS)

Collombet, M.; Burgisser, A.; Chevalier, L. A. C.

2017-12-01

The eruption of viscous magma at the Earth's surface often gives rise to abrupt regime changes. The transition from the gentle effusion of a lava dome to brief but powerful explosions is a common regime change. This transition is often preceded by the sealing of the shallow part of the volcanic conduit and the accumulation of volatile-rich magma underneath, a situation that collects the energy to be brutally released during the subsequent explosion. While conduit sealing is well-documented, volatile accumulation has proven harder to characterize. In this study, we use a 2D conduit flow numerical model including gas loss within the magma and into the wallrock to follow the evolution of gas content during a regime transition. Using various initial porosity distributions, permeability laws and boundary conditions, we track the physical parameters that prevent or enhance gas escape from the magma. Our approach aims to identify the physical processes controlling eruptive transitions and to highlight the importance of using field data observations to constrain numerical models.

a Test to Prove Cloud Whitening THEORY!

NASA Astrophysics Data System (ADS)

Buttram, J. W.

2011-12-01

Climate science researchers believe our planet can possibly tolerate twice the present carbon dioxide levels with no upwards temperature change, IF we could increase the amount of energy reflected back out into space by about 2.0%. (c)Cloudtec basically alters a blend of seawater and applies heat derived from magma to it at a temperature exceeding 2,000 degrees F. The interaction of seawater and magma displaces the oxygen, causing the volume of water to vaporize and expand over 4,000 times - transforming billions of tons of seawater into thousands of cubic miles of white, maritime, stratocumulus clouds to reflect the incident Sun's rays back out into space. A 6 month test to prove Cloud Whitening Theory will cost 6 million dollars. (No profit added.) This study will enable everyone on the planet with a computer the transparency to use satellite imagery and check out for themselves - if and when Cloud Whitening is occurring. If Cloud Whitening Theory is validated, (c)Cloudtec's innovation can strategically create the clouds we need to reflect the Sun's rays back out into space and help neutralize the projected 3.6 degrees F rise in temperature. Based on reasonable calculations of anthropogenic global warming: this one move alone would be comparable to slashing global carbon dioxide emissions by over 60% over the next 40 years.

Magma replenishment and volcanic unrest inferred from the analysis of VT micro-seismicity and seismic velocity changes at Piton de la Fournaise Volcano

NASA Astrophysics Data System (ADS)

Brenguier, F.; Rivemale, E.; Clarke, D. S.; Schmid, A.; Got, J.; Battaglia, J.; Taisne, B.; Staudacher, T.; Peltier, A.; Shapiro, N. M.; Tait, S.; Ferrazzini, V.; Di Muro, A.

2011-12-01

Piton de la Fournaise volcano (PdF) is among the most active basaltic volcanoes worldwide with more than one eruption per year on average. Also, PdF is densely instrumented with short-period and broad-band seismometers as well as with GPS receivers. Continuous seismic waveforms are available from 1999. Piton de la Fournaise volcano has a moderate inter-eruptive seismic activity with an average of five detected Volcano-Tectonic (VT) earthquakes per day with magnitudes ranging from 0.5 to 3.5. These earthquakes are shallow and located about 2.5 kilometers beneath the edifice surface. Volcanic unrest is captured on average a few weeks before eruptions by measurements of increased VT seismicity rate, inflation of the edifice summit, and decreased seismic velocities from correlations of seismic noise. Eruptions are usually preceded by seismic swarms of VT earthquakes. Recently, almost 50 % of seismic swarms were not followed by eruptions. Within this work, we aim to gather results from different groups of the UnderVolc research project in order to better understand the processes of deep magma transfer, volcanic unrest, and pre-eruptive magma transport initiation. Among our results, we show that the period 1999-2003 was characterized by a long-term increase of VT seismicity rate coupled with a long-term decrease of seismic velocities. These observations could indicate a long-term replenishment of the magma storage area. The relocation of ten years of inter-eruptive micro-seismicity shows a narrow (~300 m long) sub-vertical fault zone thus indicating a conduit rather than an extended magma reservoir as the shallow magma feeder system. Also, we focus on the processes of short-term volcanic unrest and prove that magma intrusions within the edifice leading to eruptions activate specific VT earthquakes that are distinct from magma intrusions that do not lead to eruptions. We thus propose that, among the different pathways of magma transport within the edifice, only one will allow magma to reach the edifice summit. Moreover, we have identified transient seismic velocity changes lasting a few weeks that could be associated with unreported lateral magma intrusions not leading to eruptions. The clustering of pre-eruptive micro-seismicity between mid 1999-2003 shows that seismic events repeat over successive seismic swarms and suggests that the magma pathway is spatially separated from the seismic faults. Also, the inversion for focal mechanisms shows dominant sub-horizontal P-axes indicating that part of the pre-eruptive micro-seismicity is due to the horizontal compressive stress induced by magma injection. Finally, the analysis of long-term GPS data recorded on the edifice flank shows a constant lateral displacement rate of 3.5 cm/year. More work will be needed in order to infer the possible mutual interactions between magma unrest and transport and the large-scale deformation of the edifice flank.

Energy-Constrained Recharge, Assimilation, and Fractional Crystallization (EC-RAχFC): A Visual Basic computer code for calculating trace element and isotope variations of open-system magmatic systems

NASA Astrophysics Data System (ADS)

Bohrson, Wendy A.; Spera, Frank J.

2007-11-01

Volcanic and plutonic rocks provide abundant evidence for complex processes that occur in magma storage and transport systems. The fingerprint of these processes, which include fractional crystallization, assimilation, and magma recharge, is captured in petrologic and geochemical characteristics of suites of cogenetic rocks. Quantitatively evaluating the relative contributions of each process requires integration of mass, species, and energy constraints, applied in a self-consistent way. The energy-constrained model Energy-Constrained Recharge, Assimilation, and Fractional Crystallization (EC-RaχFC) tracks the trace element and isotopic evolution of a magmatic system (melt + solids) undergoing simultaneous fractional crystallization, recharge, and assimilation. Mass, thermal, and compositional (trace element and isotope) output is provided for melt in the magma body, cumulates, enclaves, and anatectic (i.e., country rock) melt. Theory of the EC computational method has been presented by Spera and Bohrson (2001, 2002, 2004), and applications to natural systems have been elucidated by Bohrson and Spera (2001, 2003) and Fowler et al. (2004). The purpose of this contribution is to make the final version of the EC-RAχFC computer code available and to provide instructions for code implementation, description of input and output parameters, and estimates of typical values for some input parameters. A brief discussion highlights measures by which the user may evaluate the quality of the output and also provides some guidelines for implementing nonlinear productivity functions. The EC-RAχFC computer code is written in Visual Basic, the programming language of Excel. The code therefore launches in Excel and is compatible with both PC and MAC platforms. The code is available on the authors' Web sites http://magma.geol.ucsb.edu/and http://www.geology.cwu.edu/ecrafc) as well as in the auxiliary material.

Ardnamurchan 3D cone-sheet architecture explained by a single elongate magma chamber

PubMed Central

Burchardt, Steffi; Troll, Valentin R.; Mathieu, Lucie; Emeleus, Henry C.; Donaldson, Colin H.

2013-01-01

The Palaeogene Ardnamurchan central igneous complex, NW Scotland, was a defining place for the development of the classic concepts of cone-sheet and ring-dyke emplacement and has thus fundamentally influenced our thinking on subvolcanic structures. We have used the available structural information on Ardnamurchan to project the underlying three-dimensional (3D) cone-sheet structure. Here we show that a single elongate magma chamber likely acted as the source of the cone-sheet swarm(s) instead of the traditionally accepted model of three successive centres. This proposal is supported by the ridge-like morphology of the Ardnamurchan volcano and is consistent with the depth and elongation of the gravity anomaly underlying the peninsula. Our model challenges the traditional model of cone-sheet emplacement at Ardnamurchan that involves successive but independent centres in favour of a more dynamical one that involves a single, but elongate and progressively evolving magma chamber system. PMID:24100542

Ardnamurchan 3D cone-sheet architecture explained by a single elongate magma chamber.

PubMed

Burchardt, Steffi; Troll, Valentin R; Mathieu, Lucie; Emeleus, Henry C; Donaldson, Colin H

2013-10-08

The Palaeogene Ardnamurchan central igneous complex, NW Scotland, was a defining place for the development of the classic concepts of cone-sheet and ring-dyke emplacement and has thus fundamentally influenced our thinking on subvolcanic structures. We have used the available structural information on Ardnamurchan to project the underlying three-dimensional (3D) cone-sheet structure. Here we show that a single elongate magma chamber likely acted as the source of the cone-sheet swarm(s) instead of the traditionally accepted model of three successive centres. This proposal is supported by the ridge-like morphology of the Ardnamurchan volcano and is consistent with the depth and elongation of the gravity anomaly underlying the peninsula. Our model challenges the traditional model of cone-sheet emplacement at Ardnamurchan that involves successive but independent centres in favour of a more dynamical one that involves a single, but elongate and progressively evolving magma chamber system.

Low-frequency seismic events in a wider volcanological context

NASA Astrophysics Data System (ADS)

Neuberg, J. W.; Collombet, M.

2006-12-01

Low-frequency seismic events have been in the centre of attention for several years, particularly on volcanoes with highly viscous magmas. The ultimate aim is to detect changes in volcanic activity by identifying changes in the seismic behaviour in order to forecast an eruption, or in case of an ongoing eruption, forecast the short and longterm behaviour of the volcanic system. A major boost in recent years arose through several attempts of multi-parameter volcanic monitoring and modelling programs, which allowed multi-disciplinary groups of volcanologists to interpret seismic signals together with, e.g. ground deformation, stress field analysis and petrological information. This talk will give several examples of such multi-disciplinary projects, focussing on the joint modelling of seismic source processes for low-frequency events together with advanced magma flow models, and the signs of magma movement in the deformation and stress field at the surface.

Evolution of the Campanian Ignimbrite Magmatic System II: Trace Element and Th Isotopic Evidence for Open-System Processes

NASA Astrophysics Data System (ADS)

Bohrson, W. A.; Spera, F. J.; Fowler, S.; Belkin, H.; de Vivo, B.

2005-12-01

The Campanian Ignimbrite, a large volume (~200 km3 DRE) trachytic to phonolitic ignimbrite was deposited at ~39.3 ka and represents the largest of a number of highly explosive volcanic events in the region near Naples, Italy. Thermodynamic modeling of the major element evolution using the MELTS algorithm (see companion contribution by Fowler et al.) provides detailed information about the identity of and changes in proportions of solids along the liquid line of descent during isobaric fractional crystallization. We have derived trace element mass balance equations that explicitly accommodate changing mineral-melt bulk distribution coefficients during crystallization and also simultaneously satisfy energy and major element mass conservation. Although major element patterns are reasonably modeled assuming closed system fractional crystallization, modeling of trace elements that represent a range of behaviors (e.g. Zr, Nb, Th, U, Rb, Sm, Sr) yields trends for closed system fractionation that are distinct from those observed. These results suggest open-system processes were also important in the evolution of the Campanian magmatic system. Th isotope data yield an apparent isochron that is ~20 kyr younger than the age of the deposit, and age-corrected Th isotope data indicate that the magma body was an open-system at the time of eruption. Because open-system processes can profoundly change isotopic characteristics of a magma body, these results illustrate that it is critical to understand the contribution that open-system processes make to silicic magma bodies prior to assigning relevance to age or timescale information derived from isotope systematics. Fluid-magma interaction has been proposed as a mechanism to change isotopic and elemental characteristics of magma bodies, but an evaluation of the mass and thermal constraints on such a process suggest large-scale fluid-melt interaction at liquidus temperatures is unlikely. In the case of the magma body associated with the Campanian Ignimbrite, the most likely source of open-system signatures is assimilation of partial melts of compositionally heterogeneous basement composed of older cumulates and intrusive equivalents of volcanic activity within the Campanian region. Additional trace element modeling, explicitly evaluating the mass and energy balance effects that fluid, solids, and melt have on trace element evolution, will further elucidate the contributions of open vs. closed system processes within the Campanian magma body.

Three-dimensional elastic wave analysis of the October 2, 2004 tremor at Mount St. Helens, Washington.

NASA Astrophysics Data System (ADS)

Denlinger, R. P.

2006-12-01

On October 2, 2004, three component, broad-frequency band seismometers as far as 250 km away from Mount St. Helens volcano detected a prominent low-frequency resonance associated with the onset of eruptive volcanic activity. The energy is dominantly in the 0.5 to 10 Hz range. Since gas emissions were low, I investigate a gas-free tremor mechanism generated by sudden extension of a pressurized, cylindrical, visco- elastic magma conduit beneath the mountain as it forces its way through a brittle crust. The concept is that rapid faulting of the crust around and above the magma results in a sudden drop in resistance and, consequently, a concurrent extension of the magma in the magma column at a rate greater than magma can resupply the column. The result is a rapid pressure drop in the magma, causing the column to oscillate and radiate elastic energy into the surrounding crust. The full wavefield generated by this physical mechanism is analyzed using the finite volume method of Leveque (2002), with the approximation outlined in Langseth and Leveque (2000) for hyperbolic conservation laws. In this method, the three-dimensional equations of elasticity are written as a first order set of conservation equations, with a solution vector composed of 3 velocity components and 6 stress components. The eigenvectors of the jacobian matrices of these conservation equations are used in a fourth-order Taylors series expansion of the solution vector around a small increment in time. This method is robust, allows waves to cleanly propagate off of a finite computational grid, and includes surface and interface waves. The method also allows for extremely large contrasts in elastic moduli across internal boundaries in the grid, necessary to accommodate the large variations in rigidity between a hot, visco-elastic magma at depth and the Earth's crust. Analyzing the October 2, 2004 tremor observed at Johnson Ridge Observatory, 9 km north of the mountain, for pressure drop in a 10 km long conduit, I obtain 0.3 +/- .05 MPa, which is consistent with elastic analysis of crustal deformation around the mountain during this time period. Langseth, J.O., and R.J. LeVeque, 2000, A wave propagation method for 3D hyperbolic conservation laws, J. Comp. Phys., 165, 126-166. Leveque, R.J., 2002, Finite volume methods for hyperbolic problems, Cambridge U. Press, 558 p.

Composition of basaltic lavas sampled by phase-2 of the Hawaii Scientific Drilling Project: Geochemical stratigraphy and magma types

NASA Astrophysics Data System (ADS)

Rhodes, J. M.; Vollinger, M. J.

2004-03-01

This paper presents major and trace element compositions of lavas from the entire 3098 m stratigraphic section sampled by phase-2 of the Hawaii Scientific Drilling Project. The upper 245 m are lavas from Mauna Loa volcano, and the lower 2853 m are lavas and volcanoclastic rocks from Mauna Kea volcano. These intervals are inferred to represent about 100 ka and 400 ka respectively of the eruptive history of the two volcanoes. The Mauna Loa tholeiites tend to be higher in SiO2 and lower in total iron, TiO2, alkalis, and incompatible elements at a given MgO content than Mauna Kea lavas. The transition from Mauna Loa to Mauna Kea lavas is all the more pronounced because the Mauna Loa tholeiites overlie a thin sequence of postshield Mauna Kea alkalic to transitional tholeiitic lavas. The Mauna Loa tholeiites display well-developed coherent trends with MgO that are indistinguishable in most respects from modern lavas. With depth, however, there is a slight decline in incompatible element abundances, and small shifts to depleted isotopic ratios. These characteristics suggest small changes in melt production and source components over time, superimposed on shallow melt segregation. The Mauna Kea section is subdivided into a thin, upper 107 m sequence of postshield tholeiites, transitional tholeiites and alkali basalts of the Hamakua volcanics, overlying four tholeiitic magma types that are intercalated throughout the rest of the core. These four magma types are recognized on the basis of MgO-normalized SiO2 and Zr/Nb values. Type-1 lavas (high SiO2 and Zr/Nb) are ubiquitous below the postshield lavas and are the dominant magma type on Mauna Kea. They are inter-layered with the other three lava types. Type-2 lavas (low SiO2 but high Zr/Nb) are found only in the upper core, and especially above 850 m. Type-3 lavas (low SiO2 and Zr/Nb) are very similar to tholeiites from Loihi volcano and are present only below 1974 m. There are only 3 discrete samples of type-4 lavas (high SiO2 and low Zr/Nb), which are present in the upper and lower core. The differences between these magma types are inferred to reflect changes in melt production, depth of melt segregation, and differences in plume source components over about 400 ka of Mauna Kea's eruptive history. At the start of this record, eruption rates were high, and two distinct tholeiitic magmas (type-1 and 3) were erupting concurrently. These two magmas require two distinct source components, one similar to that of modern Loihi tholeiites and the other close to that of Kilauea magmas. Subsequently, the Loihi-like source of the type-3 magmas was exhausted, and these lavas are absent from the remainder of the core. For the next 200 ka or so, the eruptive sequence consists of inter-layered type-1 and -2 lavas that are derived from a common Mauna Kea source, the major difference between the two being the depth at which the melts segregated from the source. At around 440 ka (corresponding with the transition in the core from submarine to subaerial lavas) eruption rates began to decline and low-MgO lavas are suddenly much more abundant in the record. Continuing gradual decline in melting and eruption rates was accompanied by a decline in normalized SiO2 content of the type-1 magmas, and the eventual onset of postshield magmatism.

Geological Model of Supercritical Geothermal Reservoir on the Top of the Magma Chamber

NASA Astrophysics Data System (ADS)

Tsuchiya, N.

2017-12-01

We are conducting supercritical geothermal project, and deep drilling project named as "JBBP: Japan Beyond Brittle Project" The temperatures of geothermal fields operating in Japan range from 200 to 300 °C (average 250 °C), and the depths range from 1000 to 2000 m (average 1500 m). In conventional geothermal reservoirs, the mechanical behavior of the rocks is presumed to be brittle, and convection of the hydrothermal fluid through existing network is the main method of circulation in the reservoir. In order to minimize induced seismicity, a rock mass that is "beyond brittle" is one possible candidate, because the rock mechanics of "beyond brittle" material is one of plastic deformation rather than brittle failure. To understand the geological model of a supercritical geothermal reservoir, granite-porphyry system, which had been formed in subduction zone, was investigated as a natural analog of the supercritical geothermal energy system. Quartz veins, hydrothermal breccia veins, and glassy veins are observed in a granitic body. The glassy veins formed at 500-550 °C under lithostatic pressures, and then pressures dropped drastically. The solubility of silica also dropped, resulting in formation of quartz veins under a hydrostatic pressure regime. Connections between the lithostatic and hydrostatic pressure regimes were key to the formation of the hydrothermal breccia veins, and the granite-porphyry system provides useful information for creation of fracture clouds in supercritical geothermal reservoirs. A granite-porphyry system, associated with hydrothermal activity and mineralization, provides a suitable natural analog for studying a deep-seated geothermal reservoir where stockwork fracture systems are created in the presence of supercritical geothermal fluids. I describe fracture networks and their formation mechanisms using petrology and fluid inclusion studies in order to understand this "beyond brittle" supercritical geothermal reservoir, and a geological model for "Beyond Brittle" and "Supercritical" geothermal reservoir, which is located at the top of magma chamber of granite-porphyry system, will be revealed.

Deep Basalt Aquifers in Orcus Patera, Elysium Basin Mars: Perspectives for Exobiology Exploration

NASA Technical Reports Server (NTRS)

Grin, E. A.; Cabrol, N. A.

1998-01-01

Direct indicators of shorelines, spillways, and terraces allowed to determine the extent of the Elysium Paleolake between the contour-lines 1000 and 500 m below the Martian datum. The Elysium Paleolake is bordered north by Orcus Patera (14N/181W), which lies west of the Tartarus Montes and Tartarus Colles. The Orcus Patera displays an ellipse-shaped collapsed caldera of 360-km long and 100-km wide. Viking topographic data show that the bottom of the caldera is located at 2500 below the Martian datum, and surrounded by a steep-walled ram art which crest is located at about 0 m elevation. Considering the localization of Orcus Patera in the Elysium paleolake, its altimetry, and the magmatic origin of this caldera, we propose the existence of a paleolake in Orcus Patera generated (a) by juvenile water from magma during the Noachian period, and (b) by intermittent influx of the Elysium Basin from Hesperian to Amazonian. Results are encouraging to consider this site as a potential high-energy source environment for microbial communities. are circumscribed by a 50-km wide lava field mapped as Noachian material. The structure of Orcus Patera represents the record of material erupted from a magmatic reservoir. The caldera is enclosed by steep inner walls (25% measured from topographic data), values which could be in agreement with the presence of a deep magmatic reservoir, as suggested by the typology of Crumpler et.al. The depth of the caldera might be due to the collapse of the magma reservoir, and the release of gases accompanying the magma thermal evolution. Origins of water for the paleolake(s): The water that generated a paleolake in Orcus Patera may have come from two origins: (1) Juvenile water: Plescia and Crips estimated a magma H20 content by weight between 0.5% and 1.5% using for the first value a comparison with terrestrial basalt, and for the second values from a Martian meteorite. The amount of H20 can be estimated by the volume of erupted lava, and the lava content of the caldera. In this study, we adopt a water content of 1%. The total volume of magma that has been contained in the caldera, and the volume of lava contained in the observed lava field is about 110 x 10(exp 6) cubic km, that gives a total volume of 1.10 x 10(exp 6) cubic km of water. The juvenile water expelled by the overpressure within the magma chamber charged with desolved water-vapor may have moved into the crust. The decrease in overburden pressure led to bubble formation. The ascent of these bubbles generated a pressurization of the magma, which was sufficient to fracture the overlaying magma layer, (2) Water from Elysium paleolake. During the Amazonian, the rise of the Elysium paleolake level generated an overspilling that supplied the caldera with water. The southern portion of the crest shows a deep gap 12-km wide at -1500 m elevation, locating the gap between 500 to 1000 in below the assumed water of Elysium paleolake, thus facilitating the influx of Elysium paleolake water into Orcus Patera. Bathymetric calculations give a floor area of 25,500 sq km at -2000 m elevation, and a water volume of 42,000 cubic km, with a lake-level at -1500 m. A substantial amount of water may have percolated through the fractured lava, and part of the volume may have overspilled the northern crest of Orcus Patera to debouch in the Tartarus Montes region. We envision the formation of a subsurface aqueous environment in basaltic rocks at the contact of the two water-source origins, possibly the percolating surface lake water, and more likely the juvenile water. Similarly to terrestrial calderas, Orcus Patera might be surrounded by ring-fractures caused by the collapse of the magma chamber that followed the release of gases. These ring-fractures may have been covered later by sedimentation in the caldera (lacustrine, aeolian, and volcanic), and by mass wasting. The detumescence of the magma in the caldera, and the vesiculation of the juvenile water may have operated simultaneously. Comparatively to terrestrial melts, Martian iron-rich melts are denser. This greater density implies greater effusion rates (eight-times terrestrial values), and larger fissuration widths (two-times terrestrial ones). With increasing vesiculation of magma, the bubbles interact with one-another because there are of similar pressure. They make a magma froth at the contact with the caldera surface, and on the walls of the fractures. In the saturated magma, froth, where the volume ratio of gases-to-liquid is about 4:1, the bubbles form a huge surface area of interconnected spaces. Bubbles near the caldera surface disrupt the magma, and fragmentation takes place, which moves downward through the magma column. On Earth, the bubbles are likely to grow between 1 and 50 mm in diameter due to the difference between the magma surface tension, and the bubble supersaturation pressure. The Martian low-pressure at surface level is likely to accelerate the expansion of the bubbles, and increase their final diameter and number, creating more voids in the magma. The strong magma froth with enclosed juvenile water bubbles interconnected with exsolved gas bubbles constitute a potential geothermal environment for geochemical energy production from basalt and water that does not require excessive temperatures. This process can start at +20C. Similar types of environments have been shown on Earth as potential energy sources for microbial metabolism, and could have provided deep aqueous basaltic niches for possible Martian microorganisms, even geologically recently. During the Amazonian, combination of volcanism and water activity still existed on Mars. Moreover, this type of potential niches open ways for investigation of possible oases of extinct or extant life, not only on paleolakes, and surface hydrothermalism spring areas, but also all large systems of fossae, which combine hydrologic and volcanic activities, and which provide an energy source, and an underground shelter to prevent surface UV bombardment. Additional information contained in the original.

Ultrafast syn-eruptive degassing and ascent trigger high-energy basic eruptions.

PubMed

Giuffrida, Marisa; Viccaro, Marco; Ottolini, Luisa

2018-01-09

Lithium gradients in plagioclase are capable of recording extremely short-lived processes associated with gas loss from magmas prior to extrusion at the surface. We present SIMS profiles of the 7 Li/ 30 Si ion ratio in plagioclase crystals from products of the paroxysmal sequence that occurred in the period 2011-2013 at Mt. Etna (Italy) in an attempt to constrain the final ascent and degassing processes leading to these powerful eruptions involving basic magma. The observed Li concentrations reflect cycles of Li addition to the melt through gas flushing, and a syn-eruptive stage of magma degassing driven by decompression that finally produce significant Li depletion from the melt. Modeling the decreases in Li concentration in plagioclase by diffusion allowed determination of magma ascent timescales that are on the order of minutes or less. Knowledge of the storage depth beneath the volcano has led to the quantification of a mean magma ascent velocity of ~43 m/s for paroxysmal eruptions at Etna. The importance of these results relies on the application of methods, recently used exclusively for closed-system volcanoes producing violent eruptions, to open-conduit systems that have generally quiet eruptive periods of activity sometimes interrupted by sudden re-awakening and the production of anomalously energetic eruptions.

Bulk Viscosity of Bubbly Magmas and the Amplification of Pressure Waves

NASA Astrophysics Data System (ADS)

Navon, O.; Lensky, N. G.; Neuberg, J. W.; Lyakhovsky, V.

2001-12-01

The bulk viscosity of magma is needed in order to describe the dynamics of a compressible bubbly magma flowing in conduits and to follow the attenuation of pressure waves travelling through a compressible magma. We developed a model for the bulk viscosity of a suspension of gas bubbles in an incompressible Newtonian liquid that exsolves volatiles (e.g. magma). The suspension is modeled as a close pack of spherical cells, consisting of gas bubbles centered in spherical shells of a volatile-bearing liquid. Following a drop in the ambient pressure the resulting dilatational motion and driving pressure are obtained in terms of the two-phase cell parameters, i.e. bubble radius and gas pressure. By definition, the bulk viscosity of a fluid is the relation between changes of the driving pressure with respect to changes in the resulted expansion strain-rate. Thus, we can use the two-phase solution to define the bulk viscosity of a hypothetical cell, composed of a homogeneously compressible, one-phase, continuous fluid. The resulted bulk viscosity is highly non-linear. At the beginning of the expansion process, when gas exsolution is efficient, the expansion rate grows exponentially while the driving pressure decreases slightly. That means that bulk viscosity is formally negative. The negative value reflects the release of the energy stored in the supersaturated liquid (melt) and its conversion to mechanical work during exsolution. Later, when bubbles are large enough and the gas influx decreases significantly, the strain rate decelerates and the bulk viscosity becomes positive as expected in a dissipative system. We demonstrate that amplification of seismic wave travelling through a volcanic conduit filled with a volatile saturated magma may be attributed to the negative bulk viscosity of the compressible magma. Amplification of an expansion wave may, at some level in the conduit, damage the conduit walls and initiate opening of new pathways for magma to erupt.

The dynamics of magma chamber refilling at the Campi Flegrei caldera.

NASA Astrophysics Data System (ADS)

Montagna, Chiara Paola; Vassalli, Melissa; Longo, Antonella; Papale, Paolo; Giudice, Salvatore; Saccorotti, Gilberto

2010-05-01

The volcanologic and petrologic reconstructions of several eruptions during the last tens of thousand years of volcanism at the Campi Flegrei caldera show that in most cases a small, chemically evolved, partially degassed magma chamber was refilled by magma of deeper origin shortly before the eruption. New magma input in a shallow chamber is revealed from a variety of indicators, well described in the literature, that include major-trace element and isotope heterogeneities, and crystal-liquid disequilibria (e.g., Arienzo et al., Bull. Volcanol., 2009). In the case of the 4100 BP Agnano Monte Spina eruption, representing the highest intensity and magnitude event of the last epoch of activity, it has been suggested that the refilling occurred within a few tens of hours from the start of the eruption. Notably, in such a case the two end-member magmas that mixed shortly before eruption onset are not recognized as individual members in the deposits, rather, their composition and characteristics are reconstructed from small scale disequilibria, revealing that a relatively short time was sufficient for efficient mixing of the liquid components. In order to investigate the dynamics of magma chamber refilling and mixing at Campi Flegrei we have applied the GALES code (Longo et al., Geophys. Res. Lett., 2006) in a series of numerical simulations. The initial and boundary conditions have been defined in the frame of two subsequent projects coordinated by INGV and funded by the Italian Civil Protection Department, that gather a large number of experts on Campi Flegrei, and are consistent with the bulk of knowledge on the deep magmatic system. In all cases an initial compositional interface is placed at a certain depth, with non-degassed, buoyant magma placed below. The simulations investigate both the dynamics in a very large, 8 km deep reservoir revealed by seismic tomography (Zollo et al., Geophys. Res. Lett., 2008), and those in shallower and smaller chamber systems connected by dykes and representative of pre-eruptive conditions. The numerical results reveal the complex dynamics of magma mixing, dominated by the interplay between buoyant magma rise and dense magma sinking. In all simulated cases efficient mixing takes place at dyke levels, the buoyant magma entering the chamber is already a mixture of the two initial end-members, and the initial deep magma is never found as an individual component in the chamber. Over the time scale of our longest simulation (about 8 hours of real time), and with reference to the spatial resolution of our simulations (max 1 m), the magma chamber is occupied by a nearly homogeneous mixture of the two initial end-members, with minor but still visible density stratification continuously perturbed by the rise of small buoyant plumes. Consistent with the observations, an eruption occurring a few tens of hours after new magma ingression would be characterized by a magmatic composition intermediate between the two initial end-members, that can therefore be revealed only from small-scale heterogeneities and possibly from crystal-liquid disequilibria.

Physical state of the very early Earth

NASA Astrophysics Data System (ADS)

Abe, Yutaka

1993-09-01

The earliest surface environment of the Earth is reconstructed in accordance with the planetary formation theory. Formation of an atmosphere is an inevitable consequence of Earth's formation. The atmosphere near the close of accretion is composed of 200 ˜ 300 bars of H 2 and H 2O, and several tens of bars of CO and CO 2. Either by the blanketing effect of the proto-atmosphere or heating by large planetesimal impacts a magma ocean is formed during accretion. We can distinguish three stages for the thermal evolution of the magma ocean and proto-crust. Stage 0 is characterized by a super-liquidus (or completely molten) regime near the surface. At this stage the surface of the Earth is covered by a super-liquidus magma ocean. No chemical differentiation is expected during this stage. Once the energy flux released by planet formation decreases to the 200 W/m 2 level the super-liquidus magma ocean then disappears within a time interval of 1 m.y. This is the transition from stage 0 to 1. Stage 1 is characterized by a partially molten magma ocean. In the magma ocean consisting of 20 ˜ 30% partial melt, heat transport is controlled by melt-solid separation (a type of compositional convection) rather than thermal convection. Chemical differentiation of the mantle mainly occurs in this stage. Once the energy flux drops to the 160 W/m 2 level, more than 90% of water vapor in the proto-atmosphere condense to form the proto-oceans. Several tens of bars of CO and CO 2 remain in the atmosphere just after formation of the oceans. Water oceans are occasionally evaporated by large impacts. After each such event, recondensation of the ocean takes several hundred years. Although the surface is covered by a chilled proto-crust, it is short-lived because of extensive volcanic resurfacing activity as well as meteorite impacts resurfacing. This stage ends when the energy flux drops to 0.1 ˜ 1 W/m 2 level. The duration time of stage 1 is estimated to be several hundred million years (the best estimate is about 400 m.y.). Stage 2 is characterized by solid state convection. This stage continues to the present day. One of the most important change on the proto-Earth is the transition from stage 1 to 2, which occurs several hundred million years after the Earth formation. Long-lived crust is formed only after this transition.

Development of innovative biopolymers and related composites for bone tissue regeneration: study of their interaction with human osteoprogenitor cells.

PubMed

Basile, Maria Assunta; d'Ayala, Giovanna Gomez; Laurienzo, Paola; Malinconico, Mario; Della Ragione, Fulvio; Oliva, Adriana

2012-01-01

In the framework of a project aiming to improve the properties of poly(ε-caprolactone) (PCL)-based devices, we prepared novel composites and tested their in vitro biocompatibility and osteogenic capacity on human mesenchymal stromal cells (MSC) from bone marrow. We prepared two functionalized derivatives, PCL-g-MAGMA and PCL-g-DMAEA, by insertion of anhydride groups by radical grafting of maleic anhydride (MA) and glycidyl-methacrylate (GMA) molecules, and by insertion of N-(dimethylamino)ethylacrylate (DMAEA) of tertiary amines groups, respectively. In addition, in order to improve the osteoconductive properties of the materials, we also prepared the corresponding composites containing the mineral component of bone, namely hydroxyapatite (HA). Mesenchymal stromal cells (MSC) derived from bone marrow were prepared, plated onto a number of discs obtained from these functionalized derivatives and tested in terms of adhesion and vitality (by MTT test and SEM observation), and the expression of alkaline phosphatase, the early marker of osteoblastic phenotype. The biological in vitro assessment of the functionalized materials, PCL-g-MAGMA and PCL-g-DMAEA, appeared promising only in part, in particular the cells exhibited very poor adhesion to PCL-g-MAGMA. On the contrary, the related composites, PCL-g-MAGMA-HA and PCL-g-DMAEA-HA clearly showed that the addition of HA greatly ameliorated the cell-material interaction. In particular, a surprisingly increased response characterized PCL-g-MAGMA-HA, either in terms of adhesion and vitality or in terms of alkaline phosphatase activity. Altogether these studies showed that the addition of HA nanowhiskers resulted for all basic materials, in particular PCL-g-MAGMA, in improved cell adhesion and performance.

Fault-Magma Interactions during Early Continental Rifting: Seismicity of the Magadi-Natron-Manyara basins, Africa

NASA Astrophysics Data System (ADS)

Weinstein, A.; Oliva, S. J.; Ebinger, C.; Aman, M.; Lambert, C.; Roecker, S. W.; Tiberi, C.; Muirhead, J.

2017-12-01

Although magmatism may occur during the earliest stages of continental rifting, its role in strain accommodation remains weakly constrained by largely 2D studies. We analyze seismicity data from a 13-month, 39-station broadband seismic array to determine the role of magma intrusion on state-of-stress and strain localization, and their along-strike variations. Precise earthquake locations using cluster analyses and a new 3D velocity model reveal lower crustal earthquakes along projections of steep border faults that degas CO2. Seismicity forms several disks interpreted as sills at 6-10 km below a monogenetic cone field. The sills overlie a lower crustal magma chamber that may feed eruptions at Oldoinyo Lengai volcano. After determining a new ML scaling relation, we determine a b-value of 0.87 ± 0.03. Focal mechanisms for 66 earthquakes, and a longer time period of relocated earthquakes from global arrays reveal an along-axis stress rotation of 50 o ( N150 oE) in the magmatically active zone. Using Kostrov summation of local and teleseismic mechanisms, we find opening directions of N122ºE and N92ºE north and south of the magmatically active zone. The stress rotation facilitates strain transfer from border fault systems, the locus of early stage deformation, to the zone of magma intrusion in the central rift. Our seismic, structural, and geochemistry results indicate that frequent lower crustal earthquakes are promoted by elevated pore pressures from volatile degassing along border faults, and hydraulic fracture around the margins of magma bodies. Earthquakes are largely driven by stress state around inflating magma bodies, and more dike intrusions with surface faulting, eruptions, and earthquakes are expected.

Determination of Magma Ascent Rates From D/H Fractionation in Olivine-Hosted Melt Inclusions

NASA Astrophysics Data System (ADS)

Gaetani, G. A.; Bucholz, C. E.; Le Roux, V.; Klein, F.; Ghiorso, M. S.; Wallace, P. J.; Sims, K. W. W.

2016-12-01

The depths at which magmas are stored and the rates at which they ascend to Earth's surface are important controls on the dynamics of volcanic eruptions. Eruptive style is influenced by the rate at which magma ascends from the reservoir to the surface through its effect on vapor bubble nucleation, growth, and coalescence. However, ascent rates are difficult to quantify because few accurate geospeedometers are appropriate for a process occurring on such short timescales. We developed a new approach to determining ascent rates on the basis of D/H fraction associated with diffusive H2O loss from olivine-hosted melt inclusions. The utility of this approach was demonstrated on olivine-hosted melt inclusions in a hyaloclastite recovered from within Dry Valley Drilling Project core 3 from Hut Point Peninsula, Antarctica. All of the melt inclusions are glassy and contain vapor bubbles. The volumes of melt inclusions and vapor bubbles were determined by X-ray microtomography, and the density of CO2 within each bubble was determined using Raman spectroscopy. Olivines were then polished to expose individual inclusions and analyzed for volatiles and dDVSMOW by secondary ion mass spectrometry. Total CO2 was reconstructed by summing CO2 in the included glass and vapor bubble. Entrapment pressures calculated on the basis of reconstructed CO2 and maximum H2O concentrations using the MagmaSat solubility model [1] indicate a depth of origin of 24 km - in good agreement with the seismically determined depth to the Moho beneath Ross Island [2]. Magma ascent rates were determined using a finite difference model for melt inclusion dehydration during magma ascent. The positive correlation between H2O and CO2 is consistent with diffusive loss during ascent, but does not provide direct information on magma ascent rate. In contrast, the slope of the negative correlation between H2O and dDVSMOW is a reflection of transport time and, therefore, ascent rate. If it is assumed that magmas did not stall between the Moho and the surface, our results indicate an ascent rate of 0.1 m/s. Our new approach has broad applicability to determining magma ascent rates for both active and extinct volcanic centers in all tectonic environments. References: [1] Ghiorso and Gualda (2015) Cont Miner Pet 169; [2] Finotello et al. (2011) Geophys J Int 185:85-92.

Origin of reverse compositional and textural zoning in granite plutons by localized thermal overturn of stratified magma chambers

NASA Astrophysics Data System (ADS)

Trubač, Jakub; Janoušek, Vojtěch; Žák, Jiří; Somr, Michael; Kabele, Petr; Švancara, Jan; Gerdes, Axel; Žáčková, Eliška

2017-04-01

This study integrates gravimetry and thermal modelling with petrology, U-Th-Pb monazite and zircon geochronology and whole-rock geochemistry of the early Carboniferous Říčany Pluton, Bohemian Massif, in order to discuss the origin of compositional and textural zoning in granitic plutons and complex histories of horizontally stratified, multiply replenished magma chambers. The pluton consists of two coeval, nested biotite (-muscovite) granite facies: outer one, strongly porphyritic (SPm) and inner one, weakly porphyritic (WPc). Their contact is concealed but is likely gradational over several hundreds of meters. The two facies have nearly identical modal composition, are subaluminous to slightly peraluminous and geochemically evolved. Mafic microgranular enclaves, commonly associated with K-feldspar phenocryst patches, are abundant in the pluton center and indicate a repeated basic magma injection and its multistage interactions with the granitic magma and nearly solidified cumulates. Furthermore, the gravimetric data show that the nested pluton is only a small outcrop of a large anvil-like body reaching the depth of at least 14 km, where the pluton root is expected. Trace-element compositions reveal that the pluton is doubly reversely zoned. On the pluton scale, the outer SRG is geochemically more evolved than the inner WPc. On the scale of individual units, outward whole-rock geochemical variations within each facies (SPm, WPc) are compatible with fractional crystallization dominated by feldspars. The proposed genetic model invokes vertical overturn of a deeper, horizontally stratified anvil-shaped magma chamber. The overturn was driven by reactivation of resident felsic magma from the K-feldspar-rich crystal mush. The energy for the melt remobilization, extraction and subsequent ascent is thought to be provided by a long-lived thermal anomaly above the pluton feeding zone, enhanced by the multiple injections of hot basic magmas. In general, it is concluded that the three-dimensional shape of the granitic bodies exerts a first-order control on their cooling histories and thus also on their physico-chemical evolution. Thicker and longer lived portions of magma chambers are the favourable sites for extensive fractionation and/or, potentially vigorous interaction with the basic magmas. These hot domains are then particularly prone to rejuvenation and subsequent extraction of highly mobile magma leading potentially to volcanic eruptions.

Multiple magma emplacement and its effect on the superficial deformation: hints from analogue models

NASA Astrophysics Data System (ADS)

Montanari, Domenico; Bonini, Marco; Corti, Giacomo; del Ventisette, Chiara

2017-04-01

To test the effect exerted by multiple magma emplacement on the deformation pattern, we have run analogue models with synchronous, as well as diachronous magma injection from different, aligned inlets. The distance between injection points, as well as the activation in time of injection points was varied for each model. Our model results show how the position and activation in time of injection points (which reproduce multiple magma batches in nature) strongly influence model evolution. In the case of synchronous injection at different inlets, the intrusions and associated surface deformation were elongated. Forced folds and annular bounding reverse faults were quite elliptical, and with the main axis of the elongated dome trending sub-parallel to the direction of the magma input points. Model results also indicate that the injection from multiple aligned sources could reproduce the same features of systems associated with planar feeder dikes, thereby suggesting that caution should be taken when trying to infer the feeding areas on the basis of the deformation features observed at the surface or in seismic profiles. Diachronous injection from different injection points showed that the deformation observed at surface does not necessarily reflect the location and/or geometry of their feeders. Most notably, these experiments suggest that coeval magma injection from different sources favor the lateral migration of magma rather than the vertical growth, promoting the development of laterally interconnected intrusions. Recently, some authors (Magee et al., 2014, 2016; Schofield et al., 2015) have suggested that, based on seismic reflection data analysis, interconnected sills and inclined sheets can facilitate the transport of magma over great vertical distances and laterally for large distances. Intrusions and volcanoes fed by sill complexes may thus be laterally offset significantly from the melt source. Our model results strongly support these findings, by reproducing in the laboratory a strong lateral magma migration, and suggesting a possible mechanism. The models also confirmed that lateral magma migration could take place with little or no accompanying surface deformation. The research leading to these results has received funding from the European Community's Seventh Framework Programme under grant agreement No. 608553 (Project IMAGE). References: Magee et al., 2014. Basin Research, v. 26, p. 85-105, doi: 10 .1111 /bre.12044. Magee et al., 2016. Geosphere, v. 12, p. 809-841, ISSN: 1553-040X. Schofield et al., 2015. Basin Research, v. 29, p. 41-63, doi:10.1111/bre.12164.

Amplification of seismic waves beneath active volcanoes

NASA Astrophysics Data System (ADS)

Navon, O.; Lensky, N. G.; Collier, L.; Neuberg, J.; Lyakhovsky, V.

2003-04-01

Long-period (LP) seismic events are typical of many volcanoes and are attributed to energy leaking from waves traveling through the volcanic conduit or along the conduit - country-rock interface. The LP events are triggered locally, at the volcanic edifice, but the source of energy for the formation of tens of events per day is not clear. Energy may be supplied by volatile-release from a supersaturated melt. If bubbles are present in equilibrium with the melt in the conduit, and the melt is suddenly decompressed, transfer of volatiles from the supersaturated melt into the bubbles transforms stored potential energy into expansion work. For example, small dome collapses may decompress the conduit by a few bars and lead to solubility decrease, exsolution of volatiles and, consequently, to work done by the expansion of the bubbles under pressure. This energy is released over a timescale that is similar to that of LP events and may amplify the original weak seismic signals associated with the collapse. Using the formulation of Lensky et al. (2002), following the decompression, when the transfer of volatiles into bubbles is fast enough, expansion accelerates and the bulk viscosity of the bubbly magma is negative. New calculations show that under such conditions a sinusoidal P-wave is amplified. We note that seismic waves created by tectonic earthquakes that are not associated with net decompression, do not lead to net release of volatiles or to net expansion. In this case, the bulk viscosity is positive and waves traveling through the magma should attenuate. The proposed model explains how weak seismic signals may be amplified as they travel through a conduit that contains supersaturated bubbly magma. It provides the general framework for amplifying volcanic seismicity such as the signals associated with long-period events.

Physical rock properties in and around a conduit zone by well-logging in the Unzen Scientific Drilling Project, Japan

USGS Publications Warehouse

Ikeda, R.; Kajiwara, T.; Omura, K.; Hickman, S.

2008-01-01

The objective of the Unzen Scientific Drilling Project (USDP) is not only to reveal the structure and eruption history of the Unzen volcano but also to clarify the ascent and degassing mechanisms of the magma conduit. Conduit drilling (USDP-4) was conducted in 2004, which targeted the magma conduit for the 1990-95 eruption. The total drilled length of USDP-4 was 1995.75??m. Geophysical well logging, including resistivity, gamma-ray, spontaneous potential, sonic-wave velocity, density, neutron porosity, and Fullbore Formation MicroImager (FMI), was conducted at each drilling stage. Variations in the physical properties of the rocks were revealed by the well-log data, which correlated with not only large-scale formation boundaries but also small-scale changes in lithology. Such variations were evident in the lava dike, pyroclastic rocks, and breccias over depth intervals ranging from 1 to 40??m. These data support previous models for structure of the lava conduit, in that they indicate the existence of alternating layers of high-resistivity and high P-wave velocity rocks corresponding to the lava dikes, in proximity to narrower zones exhibiting high porosity, low resistivity, and low P-wave velocity. These narrow, low-porosity zones are presumably higher in permeability than the adjacent rocks and may form preferential conduits for degassing during magma ascent. ?? 2008 Elsevier B.V.

Project Hotspot: Temporal Compositional Variation in Basalts of the Kimama Core and Implications for Magma Source Evolution, Snake River Scientific Drilling Project, Idaho

NASA Astrophysics Data System (ADS)

Potter, K. E.; Shervais, J. W.; Champion, D.; Duncan, R. A.; Christiansen, E. H.

2012-12-01

Project Hotspot produced continuous core from three drill sites in the Snake River plain, including 1912 m of core from the Kimama drill site on the axis of the plain. Ongoing major and trace element chemical characterization of the Kimama core and new 40Ar/39Ar and paleomagnetic age data demonstrate temporal variations in the evolution of Snake River Plain volcanism. Cyclic fluctuations in magma chemistry identify over a hundred chemically distinct basalt flow groups (comprising 550 individual lava flows) within 54 periods of volcanic activity, separated by hiatuses of decades to many millennia. From a surface age of 700 ka to a bottom-hole age of 6.5 Ma, the Kimama core records the presence of several nearly coeval but compositionally different lava flows, ranging from highly evolved lavas to non-evolved tholeiites. Determining whether Kimama lavas are genetically unrelated or extreme differentiates of a single magma batch relies upon a combination of detailed chemostratigraphy and absolute and relative age data. Age and geochemical data introduce new ideas on the role of multiple magma sources and/or differentiation processes in the development of central Snake River Plain volcanic systems. The relatively short gestation of evolved liquids is demonstrated throughout the Kimama core, with evidence for cyclic fractionation of mafic lavas at depths of 318 m, 350 m, 547 m, and 1078 m. Here, highly evolved lava flows (FeOT 16.0-18.4 wt %; TiO2 3.43-4.62 wt %) are stratigraphically bounded by more primitive tholeiitic basalts (FeOT 9.9-14.8 wt%; TiO2 1.22-3.56 wt%) within the same inclination range, suggesting that cyclic fractionation is a regular feature of shield volcano development on the central Snake River Plain. Between 1.60 ± 0.13 Ma (453.5 m depth) and 1.54 ± 0.15 Ma (320.0 m depth), Kimama lavas ranged in composition from primitive tholeiite (FeOT 11.7 wt %; TiO2 1.76 wt %) to evolved basalt (FeOT 16.0 wt %; TiO2 4.00 wt %). At depths of 1119 m and 1138 m, evolved lava flows (FeOT 17.2 and 17.0 wt %; TiO2 4.20 and 4.09 wt %, respectively) of negative polarity are stratigraphically bounded by more primitive tholeittic lava flows (FeOT 13.6 and 14.5 wt %; TiO2 2.92 and 3.24 wt %, respectively) of positive polarity, a chronological transition that may represent many millennia and magma source variability. Kimama core stratigraphy as well as paleomagnetic, and radiometric age data demonstrate that mafic volcanism on the central Snake River Plain has been relatively continuous for the last 6.5 Ma. The compositional variability in Kimama basalts introduces broader implications for the timing of cyclic fractionation processes and the development of regional magma sources.

Petrology and Physics of Magma Ocean Crystallization

NASA Technical Reports Server (NTRS)

Elkins-Tanton, Linda T.; Parmentier, E. M.; Hess, P. C.

2003-01-01

Early Mars is thought to have been melted significantly by the conversion of kinetic energy to heat during accretion of planetesimals. The processes of solidification of a magma ocean determine initial planetary compositional differentiation and the stability of the resulting mantle density profile. The stability and compositional heterogeneity of the mantle have significance for magmatic source regions, convective instability, and magnetic field generation. Significant progress on the dynamical problem of magma ocean crystallization has been made by a number of workers. The work done under the 2003 MFRP grant further explored the implications of early physical processes on compositional heterogeneity in Mars. Our goals were to connect early physical processes in Mars evolution with the present planet's most ancient observable characteristics, including the early, strong magnetic field, the crustal dichotomy, and the compositional characteristics of the SNC meteorite's source regions as well as their formation as isotopically distinct compositions early in Mars's evolution. We had already established a possible relationship between the major element compositions of SNC meteorite sources and processes of Martian magma ocean crystallization and overturn, and under this grant extended the analysis to the crucial trace element and isotopic SNC signatures. This study then demonstrated the ability to create and end the magnetic field through magma ocean cumulate overturn and subsequent cooling, as well as the feasibility of creating a compositionally- and volumetrically-consistent crustal dichotomy through mode-1 overturn and simultaneous adiabatic melting.

Dikes, joints, and faults in the upper mantle

NASA Astrophysics Data System (ADS)

Wilshire, H. G.; Kirby, S. H.

1989-04-01

Three different types of macroscopic fractures are recognized in upper-mantle and lower-crustal xenoliths in volcanic rocks from around the world: (1) joints that are tensile fractures not occupied by crystallized magma products (2) dikes that are tensile fractures occupied by mafic magmas crystallized to pyroxenites, gabbros or hydrous-mineral-rich rocks, (3) faults that are unfilled shear fractures with surface markings indicative of shear displacement. In addition to intra-xenolith fractures, xenoliths commonly have polygonal or faceted shapes that represent fractures exploited during incorporation of the xenoliths into the host magma that brought them to the surface. The various types of fractures are considered to have formed in response to the pressures associated with magmatic fluids and to the ambient tectonic stress field. The presence of fracture sets and crosscutting relations indicate that both magma-filled and unfilled fractures can be contemporaneous and that the local stress field can change with time, leading to repeated episodes of fracture. These observations give insight into the nature of deep fracture processes and the importance of fluid-peridotite interactions in the mantle. We suggest that unfilled fractures were opened by volatile fluids exsolved from ascending magmas to the tops of growing dikes. These volatile fluids are important because they are of low viscosity and can rapidly transmit fluid pressure to dike and fault tips and because they lower the energy and tectonic stresses required to extend macroscopic cracks and to allow sliding on pre-existing fractures. Mantle seismicity at depths of 20-65 km beneath active volcanic centers in Hawaii corresponds to the depth interval where CO 2-rich fluids are expected to be liberated from ascending basaltic magmas, suggesting that such fluids play an important role in facilitating earthquake instabilities in the presence of tectonic stresses. Other phenomena related to the fractures include permeation of peridotite by fluid inclusions derived by degassing of magmas, partial melting of peridotite and dike rocks, and metasomatic alteration of peridotite host rock by magmas emplaced in fractures. These effects of magmatism generally reduce the bulk density of peridotite and might also reduce seismic velocities. The velocity contrasts between fractured and unfractured peridotite might be detected by seismic-velocity profiling techniques.

Integrating field, textural, and geochemical monitoring to track eruption triggers and dynamics: a case study from Piton de la Fournaise

NASA Astrophysics Data System (ADS)

Gurioli, Lucia; Di Muro, Andrea; Vlastélic, Ivan; Moune, Séverine; Thivet, Simon; Valer, Marina; Villeneuve, Nicolas; Boudoire, Guillaume; Peltier, Aline; Bachèlery, Patrick; Ferrazzini, Valérie; Métrich, Nicole; Benbakkar, Mhammed; Cluzel, Nicolas; Constantin, Christophe; Devidal, Jean-Luc; Fonquernie, Claire; Hénot, Jean-Marc

2018-04-01

The 2014 eruption at Piton de la Fournaise (PdF), La Réunion, which occurred after 41 months of quiescence, began with surprisingly little precursory activity and was one of the smallest so far observed at PdF in terms of duration (less than 2 days) and volume (less than 0.4 × 106 m3). The pyroclastic material was composed of golden basaltic pumice along with fluidal, spiny iridescent and spiny opaque basaltic scoria. Density analyses performed on 200 lapilli reveal that while the spiny opaque clasts are the densest (1600 kg m-3) and most crystalline (55 vol. %), the golden pumices are the least dense (400 kg m-3) and crystalline (8 vol. %). The connectivity data indicate that the fluidal and golden (Hawaiian-like) clasts have more isolated vesicles (up to 40 vol. %) than the spiny (Strombolian-like) clasts (0-5 vol. %). These textural variations are linked to primary pre-eruptive magma storage conditions. The golden and fluidal fragments track the hotter portion of the melt, in contrast to the spiny fragments and lava that mirror the cooler portion of the shallow reservoir. Exponential decay of the magma ascent and output rates through time revealed depressurization of the source during which a stratified storage system was progressively tapped. Increasing syn-eruptive degassing and melt-gas decoupling led to a decrease in the explosive intensity from early fountaining to Strombolian activity. The geochemical results confirm the absence of new input of hot magma into the 2014 reservoir and confirm the emission of a single shallow, differentiated magma source, possibly related to residual magma from the November 2009 eruption. Fast volatile exsolution and crystal-melt separation (second boiling) were triggered by deep pre-eruptive magma transfer and stress field change. Our study highlights the possibility that shallow magma pockets can be quickly reactivated by deep processes without mass or energy (heat) transfer and produce hazardous eruptions with only short-term elusive precursors.

Impact-induced melting and heating of planetary interiors - implications for the thermo-chemical evolution of planets and crystallization of magma oceans

NASA Astrophysics Data System (ADS)

Wuennemann, K.; Manske, L.; Zhu, M.; Nakajima, M.; Breuer, D.; Schwinger, S.; Plesa, A. C.

2017-12-01

Large collisions and giant impact events play an important role in the thermo-chemical evolution of planets during their early and late accretion phases. Besides material that is delivered by differentiated and primitive projectiles a significant amount of the kinetic impact energy is transferred to the planets interior resulting in heating and widespread melting of matter. As a consequence, giant impacts are thought to form global magma oceans. The amount and distribution of impact-induced heating and melting has been previously estimated by scaling laws derived from small-scale impact simulations and experiments, simple theoretical considerations, and observations at terrestrial craters. We carried out a suite of numerical models using the iSALE shock physics code and an SPH code combined with the ANEOS package to investigate the melt production in giant impacts and planetary collision events as a function of impactor size and velocity, and the target temperature. Our results are consistent with previously derived scaling laws only for smaller impactors (<10 km in diameter), but significantly deviate for larger impactors: (1) for hot planets, where the temperature below the lithosphere lies close to the solidus temperature, the melt production is significantly increased for impactors comparable in the size to the depth of the lithosphere. The resulting crater structures would drown in their own melt and only large igneous provinces (local magma oceans) would remain visible at the surface;(2) even bigger impacts (planetary collisions) generate global magma oceans; (3) impacts into a completely solidified (cold) target result in more localized heating in comparison to impacts into a magma ocean, where the impact-induced heating is distributed over a larger volume. In addition, we investigate the influence of impacts on a cooling and crystallization of magma oceans and use the lunar magma ocean as an example.

Cooling of a magmatic system under thermal chaotic mixing

NASA Astrophysics Data System (ADS)

Petrelli, Maurizio; El Omari, Kamal; Le Guer, Yves; Perugini, Diego

2015-04-01

The cooling of a melt undergoing chaotic advection is studied numerically for a magma with a temperature-dependent viscosity in a 2D cavity with moving boundary. Different statistical mixing and energy indicators are used to characterize the efficiency of cooling by thermal chaotic mixing. We show that different cooling rates can be obtained during the thermal mixing even of a single basaltic magmatic batch undergoing chaotic advection. This process can induce complex temperature patterns inside the magma chamber. The emergence of chaotic dynamics strongly affects the temperature field during time and greatly increases the cooling rates. This mechanism has implications for the lifetime of a magmatic body and may favor the appearance of chemical heterogeneities in igneous systems as a result of different crystallization rates. Results from this study also highlight that even a single magma batch can develop, under chaotic thermal advection, complex thermal and therefore compositional patterns resulting from different cooling rates, which can account for some natural features that, to date, have received unsatisfactory explanations. Among them, the production of magmatic enclaves showing completely different cooling histories compared with the host magma, compositional zoning in mineral phases, and the generation of large-scale compositionally zoning observed in many plutons worldwide.

Survey of possibility for volcanic energy development

NASA Astrophysics Data System (ADS)

1990-03-01

Volcanic areas, clarification of heat source structure, evaluation of resources and problems on utilization techniques were arranged to search the possibility of future volcanic heat source. It is necessary to improve the exploration accuracy by combining geophysical exploration with geological and geochemical surveys in order to explorate a magma reservoir. Especially, seismic exploration is effective. The surveying procedure is as follows: confirmation of magma existence and grasping the whole image, evaluation of resources, clarification of three-dimensional distribution of magma in a promising area, and heat structure survey by heat flow measurement and others to construct more accurate model for resources. This model is verified finally by practical drilling. Promising areas which are worthy of development, are active volcanic areas in Kyushu, Hakkoda nad Hokkaido. It is desirable to make drilling to the depth of 3 km or magma reservoir to develop the future heat source. It is also required to improve the thermal resistance and corrosion resistance of materials to be used. Heat extraction by a single well is most realistic and the closed coaxial double pipe heat exchanger or open heat exchanger in the well will be used to improve the extraction.

Crustal forensics in arc magmas

NASA Astrophysics Data System (ADS)

Davidson, Jon P.; Hora, John M.; Garrison, Jennifer M.; Dungan, Michael A.

2005-01-01

The geochemical characteristics of continental crust are present in nearly all arc magmas. These characteristics may reflect a specific source process, such as fluid fluxing, common to both arc magmas and the continental crust, and/or may reflect the incorporation of continental crust into arc magmas either at source via subducted sediment, or via contamination during differentiation. Resolving the relative mass contributions of juvenile, mantle-derived material, versus that derived from pre-existing crust of the upper plate, and providing these estimates on an element-by-element basis, is important because: (1) we want to constrain crustal growth rates; (2) we want to quantitatively track element cycling at convergent margins; and (3) we want to determine the origin of economically important elements and compounds. Traditional geochemical approaches for determining the contributions of various components to arc magmas are particularly successful when applied on a comparative basis. Studies of suites from multiple magmatic systems along arcs, for which differentiation effects can be individually constrained, can be used to extrapolate to potential source compositions. In the Lesser Antilles Arc, for example, differentiation trends from individual volcanoes are consistent with open-system evolution. However, such trends do not project back to a common primitive magma composition, suggesting that differentiation modifies magmas that were derived from distinct mantle sources. We propose that such approaches should now be complemented by petrographically constrained mineral-scale isotope and trace element analysis to unravel the contributing components to arc magmas. This innovative approach can: (1) better constrain true end-member compositions by returning wider ranges in geochemical compositions among constituent minerals than is found in whole rocks; (2) better determine magmatic evolution processes from core-rim isotopic or trace element profiles from the phases contained in magmas; and (3) constrain rates of differentiation by applying diffusion-controlled timescales to element profiles. An example from Nguaruhoe Volcano, New Zealand, underscores the importance of such a microsampling approach, showing that mineral isotopic compositions encompass wide ranges, that whole-rock isotopic compositions are consequently simply element-weighted averages of the heterogeneous crystal cargo, and that open-system evolution is proved by core-rim variations in Sr isotope ratios. Nguaruhoe is just one of many systems examined through microanalytical approaches. The overwhelming conclusion of these studies is that crystal cargoes are not truly phenocrystic, but are inherited from various sources. The implication of this realization is that the interpretation of whole-rock isotopic data, including the currently popular U-series, needs careful evaluation in the context of petrographic observations.

Heat-transfer measurements of the 1983 Kilauea lava flow

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hardee, H.C.

1983-10-07

Convective heat flow measurements of a basaltic lava flow were made during the 1983 eruption of Kilauea volcano in Hawaii. Eight field measurements of induced natural convection were made, giving heat flux values that ranged from 1.78 to 8.09 kilowatts per square meter at lava temperatures of 1088 and 1128 degrees Celsius, respectively. These field measurements of convective heat flux at subliquidus temperatures agree with previous laboratory measurements in furnace-melted samples of molten lava, and are useful for predicting heat transfer in magma bodies and for estimating heat extraction rates for magma energy.

Heat transfer measurements of the 1983 kilauea lava flow.

PubMed

Hardee, H C

1983-10-07

Convective heat flow measurements of a basaltic lava flow were made during the 1983 eruption of Kilauea volcano in Hawaii. Eight field measurements of induced natural convection were made, giving heat flux values that ranged from 1.78 to 8.09 kilowatts per square meter at lava temperatures of 1088 and 1128 degrees Celsius, respectively. These field measurements of convective heat flux at subliquidus temperatures agree with previous laboratory measurements in furnace-melted samples of molten lava, and are useful for predicting heat transfer in magma bodies and for estimating heat extraction rates for magma energy.

Miocene silicic volcanism in southwestern Idaho: Geochronology, geochemistry, and evolution of the central Snake River Plain

USGS Publications Warehouse

Bonnichsen, B.; Leeman, W.P.; Honjo, N.; McIntosh, W.C.; Godchaux, M.M.

2008-01-01

New 40Ar-39Ar geochronology, bulk rock geochemical data, and physical characteristics for representative stratigraphic sections of rhyolite ignimbrites and lavas from the west-central Snake River Plain (SRP) are combined to develop a coherent stratigraphic framework for Miocene silicic magmatism in this part of the Yellowstone 'hotspot track'. The magmatic record differs from that in areas to the west and east with regard to its unusually large extrusive volume, broad lateral scale, and extended duration. We infer that the magmatic systems developed in response to large-scale and repeated injections of basaltic magma into the crust, resulting in significant reconstitution of large volumes of the crust, wide distribution of crustal melt zones, and complex feeder systems for individual eruptive events. Some eruptive episodes or 'events' appear to be contemporaneous with major normal faulting, and perhaps catastrophic crustal foundering, that may have triggered concurrent evacuations of separate silicic magma reservoirs. This behavior and cumulative time-composition relations are difficult to relate to simple caldera-style single-source feeder systems and imply complex temporal-spatial development of the silicic magma systems. Inferred volumes and timing of mafic magma inputs, as the driving energy source, require a significant component of lithospheric extension on NNW-trending Basin and Range style faults (i.e., roughly parallel to the SW-NE orientation of the eastern SRP). This is needed to accommodate basaltic inputs at crustal levels, and is likely to play a role in generation of those magmas. Anomalously high magma production in the SRP compared to that in adjacent areas (e.g., northern Basin and Range Province) may require additional sub-lithospheric processes. ?? Springer-Verlag 2007.

Fluid-driven fracture and melt transport through lithosphere on earth and terrestrial planets

NASA Astrophysics Data System (ADS)

Fialko, Yuri Alex

Fluid-driven fracture is a fundamental geophysical phenomenon operating in planetary interiors on many scales. A few examples of geological processes involving fluid transport via self-induced fractures include melt segregation in the mantle, magma ascent through the lithosphere, crustal accretion at mid-ocean ridges and volcanic "hot spots", migration of metamorphic and sedimentary fluids in the crust, etc. Overall, fluid-driven (in particular, magma-driven) fracture plays a major role in chemical differentiation of the upper mantle. Because our ability to make direct observations of the dynamics and styles of fluid-driven fracture is quite limited, our understanding of this phenomenon relies on theoretical models that use fundamental physical principles and available field data to constrain the behavior of fluid-driven cracks at depth. This thesis proposes new and more accurate ways of theoretical and experimental description of magma transport in self-induced fractures, or dikes. Dike propagation is a complex process that involves elastic and inelastic deformation of the host rocks, rock fracture, viscous flow of magma, heat transfer, and phase transitions (e.g., rock crystallization and fusion, volatile exolution etc.). We consider relationships between different physical processes associated with magma transport in dikes by solving appropriate boundary value problems of continuum mechanics and heat and mass transfer. The first chapter of this thesis revises existing interpretations of available experimental data bearing on the role of fracture resistance in the overall energy balance during dike propagation. It is shown for the first time that the experimental data indicate that the rock tensile fracture energy, which is not a material property at elevated confining pressures, may substantially increase under in-situ stress conditions. The second chapter concentrates on the interaction between magma flow, heat transfer and phase changes associated with dike emplacement, and discusses some important implications of our results for the generation of the Earth's crust at mid-ocean ridges. In particular, we find that the thermal arrest lengths of typical mid-ocean ridge dikes are of the order of the wavelength of crustal thickness variations and transform fault spacing along slow spreading ridges. This suggests that thermal controls on the crustal melt delivery system could be an important factor in modulating these variations. The third chapter deals with fluid-mechanical aspects of lateral dike propagation in volcanic rift zones. We demonstrate the existence of a feedback between viscous pressure losses during magma transport at depth and the along-strike surface topography of a rift zone. Our estimated values of the along-strike slopes resulting from such a feedback are in general agreement with observations in Hawaiian rift zones. The fourth chapter explores mechanisms of emplacement of giant dike swarms that might have played a role in splitting continents and producing mass extinctions. We reconcile field observations of chilled margins, low crustal contamination, and large dike thicknesses with the theoretically inferred turbulent mode of magma flow in such dikes.

A Simulated Chlorine-Saturated Lunar Magmatic System at the Surface and At Depth

NASA Astrophysics Data System (ADS)

DiFrancesco, N.; Nekvasil, H.; Lindsley, D. H.

2016-12-01

Analysis of igneous minerals present in lunar rocks has provided evidence that volatiles such as water, chlorine and fluorine were concentrated in melts present at or near the lunar surface. While at depth, pressure on a magma allows these gases to remain dissolved in a silicate liquid, however as the magma ascends and depressurizes, these components become saturated and begin exsolving. While at pressure, it's possible for these components, specifically Cl, to form complexes in the melt with major cations such as Na, K, and Fe as well as trace elements such as Zn and Li. While dissolved in the melt, it may be possible for the Cl to inhibit the ability for these cations to enter into crystalline phases such as olivine, plagioclase, or pyroxene, potentially altering the composition of minerals associated with the melt. As the magma rises, these compounds are able to boil off from the magma, changing its bulk composition by effectively removing these cations as halides in a vapor phase. The goals of this project are to experimentally ascertain the nature of minerals sublimated by this degassing, and the effects that this process may have on the evolution and liquid line of decent for a cooling lunar magma. This is accomplished by crystallizing volatile-rich synthetic lunar basalts both at high and zero pressure and analyzing both vapor deposits and solidified liquids. Experimental data simulating volatile-rich magma degassing and crystallization at the lunar surface, and within the lunar crust has demonstrated that typical KREEP basalts (potentially rich in Cl) will crystallize more magnesian and calcic phases at high pressure, and subsequently lose alkalis and iron to a vapor phase at low pressure. We see evidence of vapor deposits and volatile element enrichment in returned Apollo samples such as "Rusty Rock", and on the surface of orange glass beads.

Seismic Forecasting of Eruptions at Dormant StratoVolcanoes

NASA Astrophysics Data System (ADS)

White, R. A.

2015-12-01

Seismic monitoring data provide important constraints on tracking magmatic ascent and eruption. Based on direct experience with over 25 and review of over 10 additional eruption sequences at 24 volcanoes, we have identified 4 phases of precursory seismicity. 1) Deep (>20 km) low frequency (DLF) earthquakes occur near the base of the crust as magma rises toward crustal reservoirs. This seismicity is the most difficult to observe, owing to generally small magnitudes (M<2.5) the significant depth. 2) Distal volcano-tectonic (DVT) earthquakes occur on tectonic faults from a 2 to 30+ km distance laterally from (not beneath) the eventual eruption site as magma intrudes into and rises out of upper crustal reservoirs to depths of 2-3 km. A survey of 111 eruptions of 83 previously dormant volcanoes, (including all eruptions of VEI >4 since 1955) shows they were all preceded by significant DVT seismicity, usually felt. This DVT seismicity is easily observed owing to magnitudes generally reaching M>3.5. The cumulative DVT energy correlates to the intruding magma volume. 3) Low frequency (LF) earthquakes, LF tremor and contained explosions occur as magma interacts with the shallow hydrothermal system (<2 km depth), while the distal seismicity dies off.4) Shortly after this, repetitive self-similar proximal seismicity may occur and may dominate the seismic records as magma rises to the surface. We present some examples of this seismic progression to demonstrate that data from a single short-period vertical station are often sufficient to forecast eruption onsets.

Seismic signature of crustal magma and fluid from deep seismic sounding data across Tengchong volcanic area

NASA Astrophysics Data System (ADS)

Bai, Z. M.; Zhang, Z. Z.; Wang, C. Y.; Klemperer, S. L.

2012-04-01

The weakened lithosphere around eastern syntax of Tibet plateau has been revealed by the Average Pn and Sn velocities, the 3D upper mantle velocity variations of P wave and S wave, and the iimaging results of magnetotelluric data. Tengchong volcanic area is neighboring to core of eastern syntax and famous for its springs, volcanic-geothermal activities and remarkable seismicity in mainland China. To probe the deep environment for the Tengchong volcanic-geothermal activity a deep seismic sounding (DSS) project was carried out across the this area in 1999. In this paper the seismic signature of crustal magma and fluid is explored from the DSS data with the seismic attribute fusion (SAF) technique, hence four possible positions for magma generation together with some locations for porous and fractured fluid beneath the Tengchong volcanic area were disclosed from the final fusion image of multi seismic attributes. The adopted attributes include the Vp, Vs and Vp/Vs results derived from a new inversion method based on the No-Ray-Tomography technique, and the migrated instantaneous attributes of central frequency, bandwidth and high frequency energy of pressure wave. Moreover, the back-projected ones which are mainly consisted by the attenuation factor Qp , the delay-time of shear wave splitting, and the amplitude ratio between S wave and P wave + S wave were also considered in this fusion process. Our fusion image indicates such a mechanism for the surface springs: a large amount of heat and the fluid released by the crystallization of magma were transmitted upward into the fluid-filled rock, and the fluid upwells along some pipeline since the high pressure in deep, thus the widespread springs of Tengchong volcanic area were developed. Moreover, the fusion image, regional volcanic and geothermal activities, and the seismicity suggest that the main risk of volcanic eruption was concentrated to the south of Tengchong city, especially around the shot point (SP) Tuantian. There are typical tectonic and deep origin mechanisms for the moderate-strong earthquakes nearby SP Tuantian, and precaution should be added on this area in case of the potential earthquake. Our fusion image also clearly revealed that there exist two remarkable positions on the Moho discontinuity through which the heat from the upper mantle was transmitted upward, and this is attributed to the widely distributed hot material within the crust and upper mantle. We acknowledge the financial support of the Ministry of Land and Resources of China (SinoProbe-02-02), and the National Nature Science Foundation of China (No. 41074033 and No. 40830315). Key Words: Seismic Signature, Magma, Tengchong Volcanic Area, Deep Seismic Sounding, Seismic Attribute Fusion Li, Chang, van der Hilst, D., Meltzer, A.S., Engdahl, E.R., 2008. Subduction of the Indian lithosphere beneath the Tibetan Plateau and Burma. Earth Planet. Sci. Lett. 274. doi:10.1016/j.epsl.2008.07.016. Lebedev, S., van der Hilst, R.D., 2008. Global upper-mantle tomography with the automated multi-mode surface and S waveforms. Geophys. J. Int. 173 (2), 505-518. Wang C.Y. and Huangfu G.. 2004. Crustal structure in Tengchong Volcano-Geothermal Area, western Yunnan, China. Tectonophysics, 380: 69-87.

STRAWBERRY CRATER ROADLESS AREAS, ARIZONA.

USGS Publications Warehouse

Wolfe, Edward W.; Light, Thomas D.

1984-01-01

The results of a mineral survey conducted in the Strawberry Crater Roadless Areas, Arizona, indicate little promise for the occurrence of metallic mineral or fossil fuel resources in the area. The area contains deposits of cinder, useful for the production of aggregate block, and for deposits of decorative stone; however, similar deposits occur in great abundance throughout the San Francisco volcanic field outside the roadless areas. There is a possibility that the Strawberry Crater Roadless Areas may overlie part of a crustal magma chamber or still warm pluton related to the San Francisco Mountain stratovolcano or to basaltic vents of late Pleistocene or Holocene age. Such a magma chamber or pluton beneath the Strawberry Crater Roadless Areas might be an energy source from which a hot-, dry-rock geothermal energy system could be developed, and a probable geothermal resource potential is therefore assigned to these areas. 9 refs.

Strawberry Crater Roadless Areas, Arizona

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wolfe, E.W.; Light, T.D.

1984-01-01

The results of a mineral survey conducted in 1980 in the Strawberry Crater Roadless Areas, Arizona, indicate little promise for the occurrence of metallic mineral or fossil fuel resources in the area. The area contains deposits of cinder, useful for the production of aggregate block, and for deposits of decorative stone; however, similar deposits occur in great abundance throughout the San Francisco volcanic field outside the roadless areas. There is a possibility that the Strawberry Crater Roadless Areas may overlie part of a crustal magma chamber or still warm pluton related to the San Francisco Mountain stratovolcano or to basalticmore » vents of late Pleistocene or Holocene age. Such a magma chamber or pluton beneath the Strawberry Crater Roadless Areas might be an energy source from which a hot-, dry-rock geothermal energy system could be developed, and a probable geothermal resource potential is therefore assigned to these areas.« less

Hydrovolcanism

NASA Technical Reports Server (NTRS)

Sheridan, M. F.; Wohletz, K. H.

1985-01-01

Hydrovolcanism is a common phenomena produced by the interaction of magma or magmatic heat with an external source of water, such as a surface body, an aquifer, or a glacier. The effects include hydrofracture of existing rock units in the subsurface and the formation of hyaloclastites in a subaqueous environment. Hydroexplosions originate within a few kilometers of the surface. They may be relatively small, phreatic events or devastating complex blasts. Large-scale experiments determined that the optimal mixing ratio of water to basaltic melt (thermite plus silicates) for efficient conversion of thermal energy into mechanical energy is in the range of 0.1 to 0.3. Based on experimental results, eruptions can be classified as dominantly magmatic if the ratio of external water to magma is less than 0.2. Eruptions with water/melt ratios in the range of 0.2 to 1.0 are highly explosive and carry tephra in a hot vapor that contains dominantly superheated (dry) steam.

Dike propagation energy balance from deformation modeling and seismic release

NASA Astrophysics Data System (ADS)

Bonaccorso, Alessandro; Aoki, Yosuke; Rivalta, Eleonora

2017-06-01

Magma is transported in the crust mainly by dike intrusions. In volcanic areas, dikes can ascend toward the free surface and also move by lateral propagation, eventually feeding flank eruptions. Understanding dike mechanics is a key to forecasting the expected propagation and associated hazard. Several studies have been conducted on dike mechanisms and propagation; however, a less in-depth investigated aspect is the relation between measured dike-induced deformation and the seismicity released during its propagation. We individuated a simple x that can be used as a proxy of the expected mechanical energy released by a propagating dike and is related to its average thickness. For several intrusions around the world (Afar, Japan, and Mount Etna), we correlate such mechanical energy to the seismic moment released by the induced earthquakes. We obtain an empirical law that quantifies the expected seismic energy released before arrest. The proposed approach may be helpful to predict the total seismic moment that will be released by an intrusion and thus to control the energy status during its propagation and the time of dike arrest.