Mixing and differentiation in the Oruanui rhyolitic magma, Taupo, New Zealand: evidence from volatiles and trace elements in melt inclusions

Large pyroclastic rhyolites are snapshots of evolving magma bodies, and preserved in their eruptive pyroclasts is a record of evolution up to the time of eruption. Here we focus on the conditions and processes in the Oruanui magma that erupted at 26.5 ka from Taupo Volcano, New Zealand. The 530 km³ (void-free) of material erupted in the Oruanui event is comparable in size to the Bishop Tuff in California, but differs in that rhyolitic pumice and glass compositions, although variable, did not change systematically with eruption order. We measured the concentrations of H₂O, CO₂ and major and trace elements in zoned phenocrysts and melt inclusions from individual pumice clasts covering the range from early to late erupted units. We also used cathodoluminescence imaging to infer growth histories of quartz phenocrysts. For quartz-hosted inclusions, we studied both fully enclosed melt inclusions and reentrants (connecting to host melt through a small opening). The textures and compositions of inclusions and phenocrysts reflect complex pre-eruptive processes of incomplete assimilation/partial melting, crystallization differentiation, magma mixing and gas saturation. ‘Restitic’ quartz occurs in seven of eight pumice clasts studied. Variations in dissolved H₂O and CO₂ in quartz-hosted melt inclusions reflect gas saturation in the Oruanui magma and crystallization depths of ∼3.5–7 km. Based on variations of dissolved H₂O and CO₂ in reentrants, the amount of exsolved gas at the beginning of eruption increased with depth, corresponding to decreasing density with depth. Pre-eruptive mixing of magma with varying gas content implies variations in magma bulk density that would have driven convective mixing. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Modes of raising northeastern Tibet probed by explosion seismology

New wide-angle reflection and refraction seismic data provide constraints on the structure of the upper lithosphere, and test models of its evolution to raise the northeastern part of the Tibetan Plateau. Amplitudes observed for reflections from the crust-mantle boundary are sufficiently large to suggest that there is no significant partial melt in the deep crust. The data show an increase of the crustal thickness between terranes from north of the Kun Lun Fault into the Qang Tang of central Tibet, and a contrast among their intracrustal images and compositions. In the north, P and S velocities are consistent with a dominantly felsic composition and show that only the upper crust thickened. South of the Kun Lun Fault a thicker crust made of two layers could result from the superposition of the originally thin crust of the Bayan Har terrane on the lower part of the crust of the domain to the north, which upper crust it shoved and thickened. Different modes of crustal thickening, either by thickening of individual layers or superpositions and imbrication among them appear to work jointly to raise the topography.     

Metamorphic petrology and zircon geochronology of high-grade rocks from the central Mozambique Belt of Tanzania: crustal recycling of Archean and Palaeoproterozoic material during the Pan-African orogeny

New data on the metamorphic petrology and zircon geochronology of high‐grade rocks in the central Mozambique Belt (MB) of Tanzania show that this part of the orogen consists of Archean and Palaeoproterozoic material that was structurally reworked during the Pan‐African event. The metamorphic rocks are characterized by a clockwise P–T path, followed by strong decompression, and the time of peak granulite facies metamorphism is similar to other granulite terranes in Tanzania. The predominant rock types are mafic to intermediate granulites, migmatites, granitoid orthogneisses and kyanite/sillimanite‐bearing metapelites. The meta‐granitoid rocks are of calc‐alkaline composition, range in age from late Archean to Neoproterozoic, and their protoliths were probably derived from magmatic arcs during collisional processes. Mafic to intermediate granulites consist of the mineral assemblage garnet–clinopyroxene–plagioclase–quartz–biotite–amphibole ± K‐feldspar ± orthopyroxene ± oxides. Metapelites are composed of garnet‐biotite‐plagioclase ± K‐feldspar ± kyanite/sillimanite ± oxides. Estimated values for peak granulite facies metamorphism are 12–13 kbar and 750–800 °C. Pressures of 5–8 kbar and temperatures of 550–700 °C characterize subsequent retrogression to amphibolite facies conditions. Evidence for a clockwise P–T path is provided by late growth of sillimanite after kyanite in metapelites. Zircon ages indicate that most of the central part of the MB in Tanzania consists of reworked ancient crust as shown by Archean (c. 2970–2500 Ma) and Palaeoproterozoic (c. 2124–1837 Ma) protolith ages. Metamorphic zircon from metapelites and granitoid orthogneisses yielded ages of c. 640 Ma which are considered to date peak regional granulite facies metamorphism during the Pan‐African orogenic event. However, the available zircon ages for the entire MB in East Africa and Madagascar also document that peak metamorphic conditions were reached at different times in different places. Large parts of the MB in central Tanzania consist of Archean and Palaeoproterozoic material that was reworked during the Pan‐African event and that may have been part of the Tanzania Craton and Usagaran domain farther to the west.     

Intercomparison of Stratospheric Chemistry Models under Polar Vortex Conditions

Several stratospheric chemistry modules from box, 2-D or 3-D models, have been intercompared. The intercomparison was focused on the ozone loss and associated reactive species under the conditions found in the cold, wintertime Arctic and Antarctic vortices. Comparisons of both gas phase and heterogeneous chemistry modules show excellent agreement between the models under constrained conditions for photolysis and the microphysics of polar stratospheric clouds. While the mean integral ozone loss ranges from 4–80% for different 30–50 days long air parcel trajectories, the mean scatter of model results around these values is only about ±1.5%. In a case study, where the models employed their standard photolysis and microphysical schemes, the variation around the mean percentage ozone loss increases to about ±7%. This increased scatter of model results is mainly due to the different treatment of the PSC microphysics and heterogeneous chemistry in the models, whereby the most unrealistic assumptions about PSC processes consequently lead to the least representative ozone chemistry. Furthermore, for this case study the model results for the ozone mixing ratios at different altitudes were compared with a measured ozone profile to investigate the extent to which models reproduce the stratospheric ozone losses. It was found that mainly in the height range of strong ozone depletion all models underestimate the ozone loss by about a factor of two. This finding corroborates earlier studies and implies a general deficiency in our understanding of the stratospheric ozone loss chemistry rather than a specific problem related to a particular model simulation.     

Finite volume method on unstructured grids in 3D with applications to the simulation of gravity waves

Summary ?A code for the simulation of atmospheric flows in 3D is presented. The underlying mathematical model is fully compressible, it takes gravity into account but Coriolis forces, turbulence and viscosity are neglected. The general numerical code consists of a finite volume discretization on unstructured hexahedral grids in 3D. The code is presently being investigated on applications to the calculation of atmospheric gravity waves on a mesh which has a structured type and is locally refined near the orography. We develop two schemes, the main difference between them lies in the different discretizations for the mass fluxes. We show that both schemes resolve typical structures of gravity waves in potential flow, linear hydrostatic motion and nonlinear non-hydrostatic regime. We compare advantages and disadvantages of the developed schemes. Received April 20, 2001; revised September 10, 2001     

Evolution of the northeastern Kaapvaal craton and the central Limpopo belt,South Africa,based on single-zircon ages (extended abstract)

Geologie en Mijnbouw -     

Hypolimnetic density currents traced by sulphur hexafluoride (SF6)

The artificial tracer sulphur hexafluoride (SF₆) has been used to study the density-driven deep water exchange between two sill-separated basins of Lake Lucerne, Gersauersee and Urnersee. The sources of the density gradients between the two basins are (1) salinity differences between the major inlets due to the different geology of their drainage areas, and (2) temperature differences due to spatial variation of wind forcing. Wind speeds are generally larger in Urnersee, especially in spring during the so-called Föhn events, when winds blow from the south. In contrast, Gersauersee is protected form these winds. In spring 1989, a total of 630 g of SF₆ was released at 80 to 120 m depth in the small Treib Basin located between Urnersee and Gersauersee. During about 100 days the distribution of SF₆ in the lake was determined by gaschromatography. Two models are used to quantify the exchange flow, (1) a one-box mass balance model for SF₆ in the deep part of Treib Basin, and (2) a one-dimensional diffusion/advection model describing the temporal and vertical temperature variation in Urnersee. According to the first model, the flow into the deep hypolimnion of Urnersee, decreases from 21·10⁶ m³·d^?1 at the end of March to about 8·10⁶ m³·d^?1 in late April. The second model yields similar flow rates. The decrease of the flow rate during spring, confirmed by both approaches, is consistent (1) with the decreasing strength of the density gradient above the sill during spring and early summer, and (2) with hydrographic information collected in Lake Lucerne during other years.     

Ambiguous Moment Tensors and Radiation Patterns in Anisotropic Media with Applications to the Modeling of Earthquake Mechanisms in W-Bohemia

Anisotropic material properties are usually neglected during inversions for source parameters of earthquakes. In general anisotropic media, however, moment tensors for pure-shear sources can exhibit significant non-double-couple components. Such effects may be erroneously interpreted as an indication for volumetric changes at the source. Here we investigate effects of anisotropy on seismic moment tensors and radiation patterns for pure-shear and tensile-type sources. Anisotropy can significantly influence the interpretation of the source mechanisms. For example, the orientation of the slip within the fault plane may affect the total seismic moment. Also, moment tensors due to pure-shear and tensile faulting can have similar characteristics depending on the orientation of the elastic tensor. Furthermore, the tensile nature of an earthquake can be obscured by near-source anisotropic properties. As an application, we consider effects of inhomogeneous anisotropic properties on the seismic moment tensor and the radiation patterns of a selected type of micro-earthquakes observed in W-Bohemia. The combined effects of near-source and along-path anisotropy cause characteristic amplitude distortions of the P, S1 and S2 waves. However, the modeling suggests that neither homogeneous nor inhomogeneous anisotropic properties alone can explain the observed large non-double-couple components.The results also indicate that a correct analysis of the source mechanism, in principle, is achievable by application of anisotropic moment tensor inversion.