Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.

Tracing pre-eruptive magma degassing using (Pb-210/Ra-226) disequilibria in the volcanic deposits of the 1980-1986 eruption of Mount St. Helens

Disequilibria between Pb-210 and Ra-226 can be used to trace magma degassing, because the intermediate nuclides, particularly Rn-222, are volatile. Products of the 1980-1986 eruptions of Mount St. Helens have been analysed for (Pb-210/Ra-226). Both excesses and deficits of Pb-210 are encountered suggesting rapid gas transfer. The time scale of diffuse, non-eruptive gas escape prior to 1980 as documented by Pb-210 deficits is on the order of a decade using the model developed by Gauthier and Condomines (Earth Planet. Sci. Lett. 172 (1999) 111-126) for a non-renewed magma chamber and efficient Rn removal. The time required to build-up Pb-210 excess is much shorter (months) as can be observed from steady increases of (Pb-210/Ra-226) with time during 1980-1982. The formation of Pb-210 excess requires both rapid gas transport through the magma and periodic blocking of gas escape routes. Superposed on this time trend is the natural variability of (Pb-210/Ra-226) in a single eruption caused by tapping magma from various depths. The two time scales of gas transport, to create both Pb-210 deficits and Pb-210 excesses, cannot be reconciled in a single event. Rather Pb-210 deficits are associated with pre-eruptive diffuse degassing, while Pb-210 excesses document the more vigorous degassing associated with eruption and recharge of the system. (c) 2006 Elsevier B.V. All rights reserved.

Life cycle of Deccan trap magma chambers: A crystal scale elemental and strontium isotopic investigation

The Deccan Trap basalts are the remnants of a massive series of lava flows that erupted at the K/T boundary and covered 1-2 million km2 of west-central India. This eruptive event is of global interest because of its possible link to the major mass extinction event, and there is much debate about the duration of this massive volcanic event. In contrast to isotopic or paleomagnetic dating methods, I explore an alternative approach to determine the lifecycle of the magma chambers that supplied the lavas, and extend the concept to obtain a tighter constraint on Deccan’s duration. My method relies on extracting time information from elemental and isotopic diffusion across zone boundaries in individual crystals. I determined elemental and Sr-isotopic variations across abnormally large (2-5 cm) plagioclase crystals from the Thalghat and Kashele “Giant Plagioclase Basalts” from the lowermost Jawhar and Igatpuri Formations respectively in the thickest Western Ghats section near Mumbai. I also obtained bulk rock major, trace and rare earth element chemistry of each lava flow from the two formations. Thalghat flows contain only 12% zoned crystals, with 87 Sr/86Sr ratios of 0.7096 in the core and 0.7106 in the rim, separated by a sharp boundary. In contrast, all Kashele crystals have a wider range of 87Sr/86Sr values, with multiple zones. Geochemical modeling of the data suggests that the two types of crystals grew in distinct magmatic environments. Modeling intracrystalline diffusive equilibration between the core and rim of Thalghat crystals led me to obtain a crystal growth rate of 2.03x10-10 cm/s and a residence time of 780 years for the crystals in the magma chamber(s). Employing some assumptions based on field and geochronologic evidence, I extrapolated this residence time to the entire Western Ghats and obtained an estimate of 25,000–35,000 years for the duration of Western Ghats volcanism. This gave an eruptive rate of 30–40 km3/yr, which is much higher than any presently erupting volcano. This result will remain speculative until a similarly detailed analytical-modeling study is performed for the rest of the Western Ghats formations.

Life Cycle of Deccan Trap Magma Chambers: A Crystal Scale Elemental and Strontium Isotopic Investigation

The Deccan Trap basalts are the remnants of a massive series of lava flows that erupted at the K/T boundary and covered 1-2 million km2 of west-central India. This eruptive event is of global interest because of its possible link to the major mass extinction event, and there is much debate about the duration of this massive volcanic event. In contrast to isotopic or paleomagnetic dating methods, I explore an alternative approach to determine the lifecycle of the magma chambers that supplied the lavas, and extend the concept to obtain a tighter constraint on Deccan’s duration. My method relies on extracting time information from elemental and isotopic diffusion across zone boundary in an individual crystal. I determined elemental and Sr-isotopic variations across abnormally large (2-5 cm) plagioclase crystals from the Thalghat and Kashele “Giant Plagioclase Basalts” from the lowermost Jawhar and Igatpuri Formations respectively in the thickest Western Ghats section near Mumbai. I also obtained bulk rock major, trace and rare earth element chemistry of each lava flow from the two formations. Thalghat flows contain only 12% zoned crystals, with 87Sr/86Sr ratios of 0.7096 in the core and 0.7106 in the rim, separated by a sharp boundary. In contrast, all Kashele crystals have a wider range of 87Sr/86Sr values, with multiple zones. Geochemical modeling of the data suggests that the two types of crystals grew in distinct magmatic environments. Modeling intracrystalline diffusive equilibration between the core and rim of Thalghat crystals led me to obtain a crystal growth rate of 2.03x10-10 cm/s and a residence time of 780 years for the crystals in the magma chamber(s). Employing some assumptions based on field and geochronologic evidence, I extrapolated this residence time to the entire Western Ghats and obtained an estimate of 25,000 – 35,000 years for the duration of Western Ghats volcanism. This gave an eruptive rate of 30 – 40 km3/yr, which is much higher than any presently erupting volcano. This result will remain speculative until a similarly detailed analytical-modeling study is performed for the rest of the Western Ghats formations.

Petrology and geochemistry of the Quaternary Caldera-forming, Phonolitic Granadilla Eruption, Tenerife (Canary Islands)

The Granadilla eruption at 600 ka was one of the largest phonolitic explosive eruptions from the Las Cañadas volcano on Tenerife, producing a classical plinian eruptive sequence of a widespread pumice fall deposit overlain by an ignimbrite. The eruption resulted in a major phase of caldera collapse that probably destroyed the shallow-level magma chamber system. Granadilla pumices contain a diverse phenocryst assemblage of alkali feldspar + biotite + sodian diopside to aegirine–augite + titanomagnetite + ilmenite + nosean/haüyne + titanite + apatite; alkali feldspar is the dominant phenocryst and biotite is the main ferromagnesian phase. Kaersutite and partially resorbed plagioclase (oligoclase to sodic andesine) are present in some eruptive units, particularly in pumice erupted during the early plinian phase, and in the Granadilla ignimbrite at the top of the sequence. Associated with the kaersutite and plagioclase are small clots of microlitic plagioclase and kaersutite interpreted as quenched blebs of tephriphonolitic magma within the phonolite pumice. The Granadilla Member has previously been recognized as an example of reverse-then-normal compositional zonation, where the zonation is primarily expressed in terms of substantial variations in trace element abundances with limited major element variation (cryptic zonation). Evidence for cryptic zonation is also provided by the chemistry of the phenocryst phases, and corresponding changes in intensive parameters (e.g. T, f O2, f  H2O). Geothermometry estimates indicate that the main body of phonolite magma had a temperature gradient from 860 °C to ∼790 °C, with hotter magma (≥900 °C) tapped at the onset and terminal phases of the eruption. The reverse-then-normal chemical and thermal zonation reflects the initial tapping of a partially hybridized magma (mixing of phonolite and tephriphonolite), followed by the more sequential tapping of a zoned and relatively large body of highly evolved phonolite at a new vent and during the main plinian phase. This suggests that the different magma types within the main holding chamber could have been laterally juxtaposed, as well as in a density-stratified arrangement. Correlations between the presence of mixed phenocryst populations (i.e. presence of plagioclase and kaersutite) and coarser pumice fall layers suggest that increased eruption vigour led to the tapping of hybridized and/or less evolved magma probably from greater depths in the chamber. New oxygen isotope data for glass and mineral separates preclude syn-eruptive interaction between the vesiculating magma and hydrothermal fluids as the cause of the Sr isotope disequilibrium identified previously for the deposit. Enrichment in radiogenic Sr in the pumice glass has more likely been due to low-temperature exchange with meteoric water that was enriched in 87Sr by sea spray, which may be a common process affecting porous and glassy pyroclastic deposits on oceanic islands.

Timing and source constraints on the relationship between mafic and felsic intrusions in the Emeiashan large igneous province

Several I- and A-type granite, syenite plutons and spatially associated, giant Fe–Ti–V deposit-bearing mafic ultramafic layered intrusions occur in the Pan–Xi(Panzhihua–Xichang) area within the inner zone of the Emeishan large igneous province (ELIP). These complexes are interpreted to be related to the Emeishan mantle plume. We present LA-ICP-MS and SIMS zircon U–Pb ages and Hf–Nd isotopic compositions for the gabbros, syenites and granites from these complexes. The dating shows that the age of the felsic intrusive magmatism (256.2 ± 3.0–259.8 ± 1.6 Ma) is indistinguishable from that of the mafic intrusive magmatism (255.4 ± 3.1–259.5 ± 2.7 Ma) and represents the final phase of a continuous magmatic episode that lasted no more than 10 Myr. The upper gabbros in the mafic–ultramafic intrusions are generally more isotopically enriched (lower eNd and eHf) than the middle and lower gabbros, suggesting that the upper gabbros have experienced a higher level of crustal contamination than the lower gabbros. The significantly positive eHf(t) values of the A-type granites and syenites (+4.9 to +10.8) are higher than those of the upper gabbros of the associated mafic intrusion, which shows that they cannot be derived by fractional crystallization of these bodies. They are however identical to those of the mafic enclaves (+7.0 to +11.4) and middle and lower gabbros, implying that they are cogenetic. We suggest that they were generated by fractionation of large-volume, plume-related basaltic magmas that ponded deep in the crust. The deep-seated magma chamber erupted in two stages: the first near a density minimum in the basaltic fractionation trend and the second during the final stage of fractionation when the magma was a low density Fe-poor, Si-rich felsic magma. The basaltic magmas emplaced in the shallowlevel magma chambers differentiated to form mafic–ultramafic layered intrusions accompanied by a small amount of crustal assimilation through roof melting. Evolved A-type granites (synenites and syenodiorites) were produced dominantly by crystallization in the deep crustal magma chamber. In contrast, the I-type granites have negative eNd(t) [-6.3 to -7.5] and eHf(t) [-1.3 to -6.7] values, with the Nd model ages (T Nd DM2) of 1.63-1.67 Ga and Hf model ages (T Hf DM2) of 1.56-1.58 Ga, suggesting that they were mainly derived from partial melting of Mesoproterozoic crust. In combination with previous studies, this study also shows that plume activity not only gave rise to reworking of ancient crust, but also significant growth of juvenile crust in the center of the ELIP.

The SEA-CALIPSO volcano imaging experiment at Montserrat : plans, campaigns at sea and on land, scientific results, and lessons learned

Since 1995 the eruption of the andesitic Soufrière Hills Volcano (SHV), Montserrat, has been studied in substantial detail. As an important contribution to this effort, the Seismic Experiment with Airgunsource-Caribbean Andesitic Lava Island Precision Seismo-geodetic Observatory (SEA-CALIPSO) experiment was devised to image the arc crust underlying Montserrat, and, if possible, the magma system at SHV using tomography and reflection seismology. Field operations were carried out in October–December 2007, with deployment of 238 seismometers on land supplementing seven volcano observatory stations, and with an array of 10 ocean-bottom seismometers deployed offshore. The RRS James Cook on NERC cruise JC19 towed a tuned airgun array plus a digital 48-channel streamer on encircling and radial tracks for 77 h about Montserrat during December 2007, firing 4414 airgun shots and yielding about 47 Gb of data. The main objecctives of the experiment were achieved. Preliminary analyses of these data published in 2010 generated images of heterogeneous high-velocity bodies representing the cores of volcanoes and subjacent intrusions, and shallow areas of low velocity on the flanks of the island that reflect volcaniclastic deposits and hydrothermal alteration. The resolution of this preliminary work did not extend beyond 5 km depth. An improved three-dimensional (3D) seismic velocity model was then obtained by inversion of 181 665 first-arrival travel times from a more-complete sampling of the dataset, yielding clear images to 7.5 km depth of a low-velocity volume that was interpreted as the magma chamber which feeds the current eruption, with an estimated volume 13 km3. Coupled thermal and seismic modelling revealed properties of the partly crystallized magma. Seismic reflection analyses aimed at imaging structures under southern Montserrat had limited success, and suggest subhorizontal layering interpreted as sills at a depth of between 6 and 19 km. Seismic reflection profiles collected offshore reveal deep fans of volcaniclastic debris and fault offsets, leading to new tectonic interpretations. This chapter presents the project goals and planning concepts, describes in detail the campaigns at sea and on land, summarizes the major results, and identifies the key lessons learned.

Magmatic evolution of topaz-bearing granite stocks within the Wiborg rapakivi granite batholith

This study brings new insights into the magmatic evolution of natural F-enriched peraluminous granitic systems. The Artjärvi, Sääskjärvi and Kymi granite stocks within the 1.64 Ga Wiborg rapakivi granite batholith have been investigated by petrographic, geochemical, experimental and melt inclusion methods. These stocks represent late-stage leucocratic and weakly peraluminous intrusive phases typical of rapakivi granites worldwide. The Artjärvi and Sääskjärvi stocks are multiphase intrusions in which the most evolved phase is topaz granite. The Kymi stock contains topaz throughout and has a well-developed zoned structure, from the rim to the center: stockscheider pegmatite equigranular topaz granite porphyritic topaz granite. Geochemically the topaz granites are enriched in F, Li, Be, Ga, Rb, Sn and Nb and depleted in Mg, Fe, Ti, Ba, Sr, Zr and Eu. The anomalous geochemistry and mineralogy of the topaz granites are essentially magmatic in origin; postmagmatic reactions have only slightly modified the compositions. The Kymi equigranular topaz granite shows the most evolved character, and the topaz granites at Artjärvi and Sääskjärvi resemble the less evolved porphyritic topaz granite of the Kymi stock. Stockscheiders are found at the roof contacts of the Artjärvi and Kymi stocks. The stockscheider at Artjärvi is composed of biotite-rich schlieren and pegmatite layers parallel to the contact. The schlieren layering is considered to have formed by velocity-gradient sorting mechanism parallel to the flow, which led to the accumulation of mafic minerals along the upper contact of the topaz granite. Cooling and contraction of the topaz granite formed fractures parallel to the roof contact and residual pegmatite magmas were injected along the fractures and formed the pegmatite layers. The zoned structure of the Kymi stock is the result of intrusion of highly evolved residual melt from deeper parts of the magma chamber along the fractured contact between the porphyritic granite crystal mush and country rock. The equigranular topaz granite and marginal pegmatite (stockscheider) crystallized from this evolved melt. Phase relations of the Kymi equigranular topaz granite have been investigated utilizing crystallization experiments at 100 to 500 MPa as a function of water activity and F content. Fluorite and topaz can crystallize as liquidus phases in F-rich peraluminous systems, but the F content of the melt should exceed 2.5 - 3.0 wt % to facilitate crystallization of topaz. In peraluminous F-bearing melts containing more than 1 wt % F, topaz and muscovite are expected to be the first F-bearing phases to crystallize at high pressure, whereas fluorite and topaz should crystallize first at low pressure. Overall, the saturation of fluorite and topaz follows the reaction: CaAl2Si2O8 (plagioclase) + 2[AlF3]melt = CaF2 (fluorite) + 2Al2SiO4F2 (topaz). The obtained partition coefficient for F between biotite and glass D(F)Bt/glass is 1.89 to 0.80 (average 1.29) and can be used as an empirical fluormeter to determine the F content of coexisting melts. In order to study the magmatic evolution of the Kymi stock, crystallized melt inclusions in quartz and topaz grains in the porphyritic and the equigranular topaz granites and the marginal pegmatite were rehomogenized and analyzed. The homogenization conditions for the melt inclusions from the granites were 700 °C, 300 MPa, and 24 h, and for melt inclusions from the pegmatite, 700 °C, 100 MPa, and 24/96 h. The majority of the melt inclusions is chemically similar to the bulk rocks (excluding H2O content), but a few melt inclusions in the equigranular granite show clearly higher F and low K2O contents (on average 11.6 wt % F, 0.65 wt % K2O). The melt inclusion compositions indicate coexistence of two melt fractions, a prevailing peraluminous and a very volatile-rich, possibly peralkaline. Combined petrological, experimental and melt inclusion studies of the Kymi equigranular topaz granite indicate that plagioclase was the liquidus phase at nearly water-saturated (fluid-saturated) conditions and that the F content of the melt was at least 2 wt %. The early crystallization of biotite and the presence of muscovite in crystallization experiments at 200 MPa contrasts with the late-stage crystallization of biotite and the absence of muscovite in the equigranular granite, indicating that crystallization pressure may have been lower than 200 MPa for the granite.

Chromite-silicate chemistry of the Neoarchean Sittampundi Complex, southern India: Implications for subduction-related arc magmatism

The Neoarchean layered anorthositic complex at Sittampundi in southern India is known for its chromitite layers that are mostly associated with anorthosite (An(90-100)). The chromitites contain FeAl-rich chromites concentrated in layers between amphibole-rich layers with a dominant mineralogy of amphibole-spinel-plagiocase+/-sapphirine. The chromite-rich layers contain only amphibole and plagioclase. Mineral compositions illustrated by X-ray composition maps and profiles show subtle chemical differences. The chrome spinels are of refractory grade with Cr2O3 and Al2O3 contents varying between 34-40 wt.% and 23-28 wt.%. The chromite compositions are noticeably different from those in layered igneous intrusions of the Bushveld-Stillwater type. The existence of original highly calcic plagioclase, FeAl-rich chromite, and magmatic amphibole is consistent with derivation from a parental magma of hydrous tholeiitic composition that was most likely generated in a supra-subduction zone arc setting. In terms of mineralogy and field relations, the Sittampundi chromitites are remarkably similar to anorthosite-hosted chromitites in the Neoarchean Fiskensset anorthositic complex, Greenland. We propose that the Sittampundi chromitites formed by partial melting of unusually aluminous harzburgite in a hydrated mantle wedge above a subduction zone. This melting process produced hydrous, aluminous basalt, which fractionated at depth to give rise to a variety of high-alumina basalt compositions from which the anorthositic complex with its cumulate chromite-rich and amphibole-rich layers formed within the magma chamber of a supra-subduction zone arc. (C) 2011 Elsevier B.V. All rights reserved.

Attenuation tomography. Modelling regional love waves: Imperial Valley to Pasadena

Abstract to Part I

The inverse problem of seismic wave attenuation is solved by an iterative back-projection method. The seismic wave quality factor, Q, can be estimated approximately by inverting the S-to-P amplitude ratios. Effects of various uncertain ties in the method are tested and the attenuation tomography is shown to be useful in solving for the spatial variations in attenuation structure and in estimating the effective seismic quality factor of attenuating anomalies.

Back-projection attenuation tomography is applied to two cases in southern California: Imperial Valley and the Coso-Indian Wells region. In the Coso-Indian Wells region, a highly attenuating body (S-wave quality factor (Q_β ≈ 30) coincides with a slow P-wave anomaly mapped by Walck and Clayton (1987). This coincidence suggests the presence of a magmatic or hydrothermal body 3 to 5 km deep in the Indian Wells region. In the Imperial Valley, slow P-wave travel-time anomalies and highly attenuating S-wave anomalies were found in the Brawley seismic zone at a depth of 8 to 12 km. The effective S-wave quality factor is very low (Q_β ≈ 20) and the P-wave velocity is 10% slower than the surrounding areas. These results suggest either magmatic or hydrothermal intrusions, or fractures at depth, possibly related to active shear in the Brawley seismic zone.

No-block inversion is a generalized tomographic method utilizing the continuous form of an inverse problem. The inverse problem of attenuation can be posed in a continuous form , and the no-block inversion technique is applied to the same data set used in the back-projection tomography. A relatively small data set with little redundancy enables us to apply both techniques to a similar degree of resolution. The results obtained by the two methods are very similar. By applying the two methods to the same data set, formal errors and resolution can be directly computed for the final model, and the objectivity of the final result can be enhanced.

Both methods of attenuation tomography are applied to a data set of local earthquakes in Kilauea, Hawaii, to solve for the attenuation structure under Kilauea and the East Rift Zone. The shallow Kilauea magma chamber, East Rift Zone and the Mauna Loa magma chamber are delineated as attenuating anomalies. Detailed inversion reveals shallow secondary magma reservoirs at Mauna Ulu and Puu Oo, the present sites of volcanic eruptions. The Hilina Fault zone is highly attenuating, dominating the attenuating anomalies at shallow depths. The magma conduit system along the summit and the East Rift Zone of Kilauea shows up as a continuous supply channel extending down to a depth of approximately 6 km. The Southwest Rift Zone, on the other hand, is not delineated by attenuating anomalies, except at a depth of 8-12 km, where an attenuating anomaly is imaged west of Puu Kou. The Ylauna Loa chamber is seated at a deeper level (about 6-10 km) than the Kilauea magma chamber. Resolution in the Mauna Loa area is not as good as in the Kilauea area, and there is a trade-off between the depth extent of the magma chamber imaged under Mauna Loa and the error that is due to poor ray coverage. Kilauea magma chamber, on the other hand, is well resolved, according to a resolution test done at the location of the magma chamber.

Abstract to Part II

Long period seismograms recorded at Pasadena of earthquakes occurring along a profile to Imperial Valley are studied in terms of source phenomena (e.g., source mechanisms and depths) versus path effects. Some of the events have known source parameters, determined by teleseismic or near-field studies, and are used as master events in a forward modeling exercise to derive the Green's functions (SH displacements at Pasadena that are due to a pure strike-slip or dip-slip mechanism) that describe the propagation effects along the profile. Both timing and waveforms of records are matched by synthetics calculated from 2-dimensional velocity models. The best 2-dimensional section begins at Imperial Valley with a thin crust containing the basin structure and thickens towards Pasadena. The detailed nature of the transition zone at the base of the crust controls the early arriving shorter periods (strong motions), while the edge of the basin controls the scattered longer period surface waves. From the waveform characteristics alone, shallow events in the basin are easily distinguished from deep events, and the amount of strike-slip versus dip-slip motion is also easily determined. Those events rupturing the sediments, such as the 1979 Imperial Valley earthquake, can be recognized easily by a late-arriving scattered Love wave that has been delayed by the very slow path across the shallow valley structure.

Deep Fractionation of Clinopyroxene in the East Pacific Rise 13 degrees N: Evidence from High MgO MORB and Melt Inclusions

Mid-ocean ridge basalts (MORBs) from East Pacific Rise (EPR) 13 degrees N are analysed for major and trace elements, both of which show a continuous evolving trend. Positive MgO-Al2O3 and negative MgO-Sc relationships manifest the cotectic crystallization of plagioclase and olivine, which exist with the presence of plagioclase and olivine phenocrysts and the absence of clinopyroxene phenocrysts. However, the fractionation of clinopyroxene is proven by the positive correlation of MgO and CaO. Thus, MORB samples are believed to show a "clinopyroxene paradox". The highest magnesium-bearing MORB sample E13-3B (MgO = 9.52%) is modelled for isobaric crystallization with COMAGMAT at different pressures. Observed CaO/Al2O3 ratios can be derived from E13-3B only by fractional crystallization at pressure > 4 +/- 1 kbar, which necessitates clinopyroxene crystallization and is not consistent with cotectic crystallization of olivine plus plagioclase in the magma chamber (at pressure similar to 1 kbar). The initial compositions of the melt inclusions, which could represent potential parental magmas, are reconstructed by correcting for post-entrapment crystallization (PEC). The simulated crystallization of initial melt inclusions also produce observed CaO/Al2O3 ratios only at > 4 +/- 1 kbar, in which clinopyroxene takes part in crystallization. It is suggested that MORB magmas have experienced clinopyroxene fractionation in the lower crust, in and below the Moho transition zone. The MORB magmas have experienced transition from clinopyroxene+plagioclase+olivine crystallization at > 4 +/- 1 kbar to mainly olivine+plagioclase crystallization at < 1 kbar, which contributes to the explanation of the "clinopyroxene paradox".

Major element, trace element, and Sr, Nd and Pb isotope studies of Cenozoic basalts from the South China Sea

The whole rock K-Ar ages of basalts from the South China Sea basin vary from 3.8 to 7.9 Ma, which suggest that intra-plate volcanism after the cessation of spreading of the South China Sea (SCS) is comparable to that in adjacent regions around the SCS, i.e., Leiqiong Peninsula, northern margin of the SCS, Indochina block, and so on. Based on detailed petrographic studies, we selected many fresh basaltic rocks and measured their major element, trace element, and Sr-Nd-Pb isotope compositions. Geochemical characteristics of major element and trace element show that these basaltic rocks belong to alkali basalt magma series, and are similar to OIB-type basalt. The extent of partial melting of mantle rock in source region is very low, and magma may experience crystallization differentiation and cumulation during the ascent to or storing in the high-level magma chamber. Sr-Nd-Pb isotopic data of these basaltic rocks imply an inhomogeneous mantle below the South China Sea. The nature of magma origin has a two end-member mixing model, one is EM2 (Enriched Mantle 2) which may be originated from mantle plume, the other is DMM (Depleted MORB Mantle). Pb isotopic characteristics show the Dupal anomaly in the South China Sea, and combined with newly found Dupal anomaly at Gakkel ridge in Arctic Ocean, this implies that Dupal anomaly is not only limited to South Hemisphere. In variation diagrams among Sr, Nd and Pb, the origin nature of mantle below the SCS is similar to those below Leiqiong peninsula, northern margin of the SCS and Indochina peninsula, and is different from those below north and northeast China. This study provides geochemical constraints on Hainan mantle plume.

Periodical mixing of MORB magmas near East Pacific Rise 13N: evidence from modeling and zoned plagioclase phenocrysts

Thirty-six basalt samples from near East Pacific Rise 13N are analyzed for major and trace elements. Different types of zoned plagioclase phenocrysts in basalts are also backscatter imaged, and major element profiles scanned and analyzed for microprobe. Basalts dredged from a restricted area have evolved to different extents (MgO=9.38wt%-6.76wt%). High MgO basalts are modeled for crystallization to MgO of about 7wt%, and resulted in the Ni contents (28 ppm) that are generally lower than that in observed basalts (> 60 ppm). It suggests that low MgO basalts may have experienced more intensive magma mixing. High MgO (9.38wt%) basalt is modeled for self-"mixing-crystallization", and the high Ni contents in low MgO basalts can be generated in small scale and periodical self-mixing of new magma (high MgO). "Mixing-crystallization" processes that low MgO magmas experienced accord with recent 226Ra/230Th disequilibria studies for magma residence time, in which low MgO magmas have experienced more circles of "mixing-crystallization" in relatively longer residence time. Magma mixing is not homogeneous in magma chamber, however, low MgO magmas are closer to stable composition produced by periodical "mixing-crystallization", which is also an important reason for magma diversity in East Pacific Rise. Zoned plagioclase phenocrysts can be divided into two types: with and without high An# cores, both of which have multiple reversed An# zones, suggesting periodical mixing of their host magmas. Cores of zoned plagioclase in low MgO (7.45wt%) basalt differ significantly with their mantle in An#, but are similar in An# with microlite cores (products of equilibrium crystallization) in high MgO (9.38wt%) basalt, which further shows that plagioclase phenocryst cores in low MgO basalts may have formed in their parental magmas before entering into the magma chamber.

东天山镁铁-超镁铁岩铜镍硫化物矿床通道式成矿机制与岩体含矿性评价研究

A mafic-ultramafic complex belt well developed in Eastern Tianshan, Xinjiang, NW China, which contains a series of Cu-Ni sulfide deposits. This area is the important production basis for Cu-Ni deposits, including Tulargen deposit, Hulu deposit, Huangshan-Huangshandong deposit, Hulu deposit, Xiangshan deposit, Tianyu deposit, Chuanzhu deposit. In China, especially Eastern Tianshan, it is prevalent that large Cu-Ni deposits occurred in small intrusions, typically including Jinchuan, Kalatongke, et al., so the ore-forming mechanism and evaluation rule for those small intrusions are very meaningful and of universal significance. On the basis of the research to typical Cu-Ni deposits, ore-forming conditions and processes are summarized through which to evaluate the ore-bearing potential for barren intrusions and unexplored mafic-ultramafic intrusions. By the contrast, metallogenic rule and mechanism of ore genesis are concluded, and evaluation system is preliminarily set up on the basis of these conclusions. Quantitatively simulation for the composition of olivine is introduced for the first time in China to discuss the interaction between magma and sulfide, and a new method to calculate the Mg-Fe composition of primitive magma is developed. Interaction between magma and sulfide liquid is used to get the Ni content in sulfide liquid. Sulfur isotopic characteristics in sulfide minerals in country rocks and ores are used to judge crustal sulfur introduction, which is applied for the first time in China. Re-Os isotopic characteristics are related to the ore-forming process, to interpret the process of enrichment of chalcophile elements. On the basis of the evaluation system, Mati, Chuanzhu, Luodong, Xiadong, those intrusions are evaluated to their ore-bearing potential. According to the studies to typical Cu-Ni deposits, conduit-type ore-forming model is set up, and the characteristics of the model are concluded systematically. The evaluation system and conduit-type ore-forming model can be helpful to the evaluation of mafic-ultramafic intrusions in this and similar mafic-ultramafic intrusion belts. The studied typical deposits and mafic-ultramafic intrusion include Tulargen deposit, Hulu deposit, Huangshandong deposit, Chuanzhu deposit, Mati intrusion，Luodong intrusion, Xiadong intrusion, and others. Through studies, there are similar characteristics for Tulargen and Hulu deposits in magma origin, composition of primitive magma（MgO=12.5%, FeO=12% and MgO=11%, FeO=10.5% respectively）, magma evolution, mechanism of sulfide segregation and conduit-type ore-forming process. By Re-Os isotopic system, the ore forming date of Tulargen deposit is 265.6±9.2Ma, which is consistent to regional metallogenic event, but little younger. The Mg-Fe composition of primitive magma of Baishiquan, Huangshandong area, Kalatongke is lower than that of Tulargen and Hulu deposit, showing common basalt composition. The Mg# value（Mg#=(Mg/Mg+Fe）increases gradually from Kalatongke to Baishiquan to Huangshan-Huangshandong East. Baishiquan intrusions show relatively higher crustal contamination by evidence of trace element, which indicates the lower magma original source, from depleted mantle to crust. One break is the discovery of komatiitic intrusion, Xiadong intrusion, which shows characteristics of highly magnesium (Max Fo=96). The primitive magma is calculated of MgO=28%，FeO=9%, belonging to komatiitic magma. Tectonic evolution of Eastern Tianshan is discussed. By the statistics of ore-forming data of porphyry copper deposits, magmatic sulfide Cu-Ni deposits, orogenic hydrothermal gold deposits, we believe that those deposits are the successive products of oceanic subduction, are and back-arc basin collision and post-orogenic extention. And Cu-Ni sulfide deposits and orogenic gold deposits occurred in the stage of post-orogenic extention. According to the conclusions, the conduit-type ore-forming mechanism of magmatic sulfide deposit is set up, and its characteristics and conditions are concluded as well. The conduit-type ore-forming system includes magma generation, sulfide segregation, enrichment of chalcophile elements, interaction of sulfide and magma, sulfide collection in limited space in magma conduit and bottom of the chamber, which make a whole ore-forming system.The ore-forming process of Cu-Ni sulfide deposits is concluded as three steps: 1. mantle derived magma rises upward to the middle-upper crust; 2. magma suffers crustal contamination of different degrees and assimilates crustal sulfur, which leads to sulfur saturation and sulfide segregation. Sulfide liquid interacts with magma and concentrates chalcophile elements; 3. enriched sulfide located in the conduit(Tulargen) or bottom of the chamber (Hulu). Depleted magma rises upward continuously to form barren complexes. For the practical cases, Tulargen deposit represents the feeding conduit, and Hulu deposit represents the bottom of the staging magma chamber. So the deeper of west of Tulargen and southwest of Hulu are the favorite locate for ore location. The evaluation for ore potential can be summarized as follows: (1) Olivine can be served as indicator for magma evolution and events of sulfide segregation; (2) Sulfur isotopic characteristics is an efficient method to judge sulfur origin for magmatic sulfide deposit; (3) Re-Os content of the ores can indicate interaction between sulfide and silicate magma and crustal contamination; (4) PGE mineralization is effected by degree of partial melting of mantle; (5) Cu/Zr is efficient parameter to judge sulfide segregation; (6) The effects of multiple magma fractionation and emplacement are important, for inverse order shows the destruction to previous solid lithofacies and orebodies. Mati, Chuanzhu, Xiadong, Luodong, mafic-ultramafic intrusions are evaluated using evaluation system above. Remarkable Ni depletion is found in olivine of Mati, and southwest of the intrusion can be hopeful location for ore location. Chuanzhu intrusion has remarkable evidence of sulfide segregation, but the intrusion represents the narrow feeder conduit, so the wide part of the conduit maybe the favorite location for sulfide to deposit. The ore potential of Luodong and Xiadong is not good. Both the intrusions show no Ni depletion in olivine, and there is no sulfide in country rocks, so no crustal sulfur is added into the magmatic system. For Sidingheishan, a very large intrusion, the phenomenon of sulfide segregation is found, but there are no favorite places for sulfide to deposit. So the Cu-Ni ore potential maybe not good, but PGE mineralization should be evaluated further.

新增生岩石圈与软流圈的相互作用：来自地幔橄榄岩捕虏体的证据

The petrology and geochemistry of peridotites entrained in Beiyan Cenozoic alkaline basalts within the middle segment of Tan-Lu fault zone and clinopyroxene megacrysts in the late Mesozoic and Cenozoic alkaline basaltic rocks from the North China Craton, have been systematically investigated. The main conclusions are obtained as follows. The peridotites entrained in alkaline basalts at Beiyan, Shandong Province, China are comprised of dominantly spinel lherzolites and spinel wehrlites with porphyroclastic, granuloblastic textures to resorption textures. The xenoliths are fertile in major element compositions (High CaO, TiO2, Low MgO, Cr2O3). The olivine Fo (= 100×Mg / (Mg+Fe) possesses a low and very large range of 81.0 to 91.0. The peridotites contain high percentages (Lherzolites: 10 - 19% in volume; Wehrlites: 24 - 28% in volume) of clinopyroxene with spongy textures. The Sr and Nd isotopic ratios of clinopyroxene separates from peridotites and pyroxenite xenoliths have a depleted and small range fall within the area of MORB, similar to newly-accreted lithospheric mantle. However, the appearance of many wehrlites and highly enriched LREE pattern suggest that this newly-accreted lithospheric mantle was considerably modified and reconstructed recently through the peridotite-asthenospheric melt interaction. The upwelling of asthenosphere from late Cretaceous to Eogene and upper mantle shearing of the Tan-Lu fault played an important role in the modification and reconstruction of the newly-accreted lithospheric mantle. The clinopyroxene megacrysts in the late Mesozoic and Cenozoic alkaline basaltic rocks from the eatern North China Craton are different in aspects of major elements, trace elements and isotopic composition. The characteristics of texture, mineral compositions and geochemistry as well as the Fe-Mg partitioning between the crystal and the melt indicates that the Al-augites in the Cenozoic basalts represent high-pressure crystallization products of alkaline basaltic melts. Thus, both of clinopyroxene megacrysts and host basalts could be derived from a same primitive magma. However, the Al-augites in the late Mesozoic basaltic rocks represent accidentally-included xenocrysts of basaltic components which had crystallized in the depth from a previously melting episode. The more depleted Sr-Nd isotopic compositions of Cenozoic megacrysts compared with those of host alkaline basalts and tholeiites demonstrate that even the alkali basalts could not completely represent primitive magma initiating in asthenosphere. That is to say, the Cenozoic alkaline basalts were more or less modified by some enriched Sr-Nd isotopic components during their eruption. Meanwhile, the tholeiites were not the products formed only by fractional crystallization of alkaline basaltic magma or different degrees of partial melting. It may result from the contribution of lithospheric mantle materials or crust contamination in magma chamber to alkali basaltic magmas.