Decoupled crystallization and eruption histories of the rhyolite magmatic system at Tarawera volcano revealed by zircon ages and growth rates

We obtained U–Th disequilibrium age data on zircons from each of the four rhyolite eruptions that built Tarawera volcano in the last 22 ka within the Okataina Volcanic Center (OVC), caldera, New Zealand. Secondary ion mass spectrometry analyses on unpolished euhedral crystal faces that lack resorption features show that crystal growth variously terminated from near-eruption age to ~100 ka prior to eruption. Age-depth profiling of crystals reveals long periods of continuous (~34 ka) and discontinuous growth (~90 ka). Growth hiatuses of up to ~40 ka duration occur, but do not all relate to obvious resorption surfaces. Age differences up to similar magnitude are found on opposing faces of some crystals suggesting episodes of partial exposure to melts. These features are best explained by periodic, complete, or partial, sub-solidus storage and/or inclusion in larger crystal phases, followed by rapid liberation prior to eruption. This is supported by high abundances of U and Th (~500 − >2,000 ppm) in some zircons consistent with periods of high crystallinity (>70%) in the magmatic system, based on crystal/melt partitioning. Contemporaneous but contrasting rim-ward trends of these elements within crystals, even in the same lava hand sample, require synchronous growth in separate melt bodies and little connectivity within the system, but also significant crystal transport and mixing prior to eruption. Many crystals record continuity of growth through the preceding ~60 ka OVC caldera-collapse and subsequent eruptions from Tarawera. This demonstrates a decoupling between eruption triggers, such as shallow crustal extension and mafic intrusion, and the crystallization state of the OVC silicic magmatic system. The data highlights the need to distinguish between the time for accumulation of eruptible magma and the long-term magma residence time based on the age of crystals with high closure temperatures, when assessing the potential for catastrophic eruptions.     

Resolution of GPR bowtie antennas: An experimental approach

Since target reflections directly depend on the emitted pulse characteristics, a key factor for carrying out a successful GPR survey is to know as much as possible about the transmission features of the antennas used. This information is very important in order to choose the right antennas and set the appropriate configuration parameters for a specific survey. With this in mind this paper deals with the development of a set of laboratory experiments on the resolution capabilities of three bowtie antennas at frequencies of 500, 800 and 1000 MHz. Results from these measurements give a first estimation of the resolution of the antennas under test, showing the advantage of performing experiments rather than relying only on theoretical assumptions. The results are also expressed in terms of the central wavelength for each antenna and compared with some theoretical estimations proposed in the specialized bibliography.     

Seismic refraction traveltime inversion for static corrections in a glaciated shield rock environment: a case study

A traveltime inversion technique is applied to model the upper ∼40 m of the subsurface of a glaciated shield rock area in order to calculate static corrections for a multi-azimuth multi-depth walk-away vertical seismic profile and a surface seismic reflection profile. First break information from a seismic refraction survey is used in conjunction with a ray-tracing program and an iterative damped least-squares inversion algorithm to create a two-dimensional model of the subsurface. The layout of the seismic survey required crooked seismic lines and substantial gaps in the source and receiver coverage to be accounted for. Additionally, there is substantial topographical variation and a complex geology consisting of glaciofluvial sediment and glacial till overlying a crystalline bedrock. The resolution and reliability of the models is measured through a parameter perturbation technique, normalized χ² values, root means square traveltime residuals and comparison to known geology.     

Holocene and Pleistocene fringing reef growth and the role of accommodation space and exposure to waves and currents (Bora Bora,Society Islands,French Polynesia)

Holocene fringing reef development around Bora Bora is controlled by variations in accommodation space (as a function of sea‐level and antecedent topography) and exposure to waves and currents. Subsidence ranged from 0 to 0·11 m kyr^?1, and did not create significant accommodation space. A windward fringing reef started to grow 8·7 kyr bp , retrograded towards the coast over a Pleistocene fringing reef until ca 6·0 kyr bp , and then prograded towards the lagoon after sea‐level had reached its present level. The retrograding portion of the reef is dominated by corals, calcareous algae and microbialite frameworks; the prograding portion is largely detrital. The reef is up to 13·5 m thick and accreted vertically with an average rate of 3·12 m kyr^?1. Lateral growth amounts to 13·3 m kyr^?1. Reef corals are dominated by an inner Pocillopora assemblage and an outer Acropora assemblage. Both assemblages comprise thick crusts of coralline algae. Palaeobathymetry suggests deposition in 0 to 10 m depth. An underlying Pleistocene fringing reef formed during the sea‐level highstand of Marine Isotope Stage 5e, and is also characterized by the occurrence of corals, coralline algal crusts and microbialites. A previously investigated, leeward fringing reef started to form contemporaneously (8·78 kyr bp ), but is thicker (up to 20 m) and solely prograded throughout the Holocene. A shallow Pocillopora assemblage and a deeper water Montipora assemblage were identified, but detrital facies dominate. At the Holocene reef base, only basalt was recovered. The Holocene windward–leeward differences are a consequence of less accommodation space on the eastern island side that eventually led to a more complex reef architecture. As a result of higher rates of exposure and flushing, the reef framework on the windward island side is more abundant and experienced stronger cementation. In the Pleistocene, the environmental conditions on the leeward island side were presumably unfavourable for fringing reef growth.     

Triassic synthems of southern South America (southwestern Gondwana) and the Western Caucasus (the northern Neotethys), and global tracing of their boundaries

Global tracing of the key surfaces of Triassic deposits may contribute significantly to the understanding of the common patterns in their accumulation. We attempt to define synthems – disconformity-bounded sedimentary complexes – in the Triassic successions of southern South America (southwestern Gondwana, Brazil and Argentina) and the Western Caucasus (the northern Neotethys, Russia), and then to trace their boundaries in the adjacent regions and globally. In southern South America, a number of synthems have been recognized – the Cuyo Basin: the Río Mendoza–Cerro de las Cabras Synthem (Olenekian–Ladinian) and the Potrerillos–Cacheuta–Río Blanco Synthem (Carnian–Rhaetian); the Ischigualasto Basin: the Ischichuca-Los Rastros Synthem (Anisian–Ladinian) and the Ischigualasto–Los Colorados Synthem (Carnian–Rhaetian); the Chaco–Paraná Basin: the Sanga do Cabral Synthem (Induan), the Santa Maria 1 Synthem (Ladinian), the Santa Maria 2 Synthem (Carnian), and the Caturrita Synthem (Norian); western Argentina: the Talampaya Synthem (Lower Triassic) and the Tarjados Synthem (Olenekian?). In the Western Caucasus, three common synthems have been distinguished: WC-1 (Induan–Anisian), WC-2 (uppermost Anisian–Carnian), and WC-3 (Norian–lower Rhaetian). The lower boundary of WC-1 corresponds to a hiatus whose duration seems to be shorter than that previously postulated. The synthem boundaries that are common to southwestern Gondwana and the Western Caucasus lie close to the base and top of the Triassic. The Lower Triassic, Ladinian, and Upper Triassic disconformities are traced within the studied basins of southern South America, and the first two are also established in South Africa. The Upper Triassic disconformity is only traced within the entire Caucasus, whereas all synthem boundaries established in the Western Caucasus are traced partly within Europe. In general, the synthem boundaries recognized in southern South America and the Western Caucasus are correlated to the global Triassic sequence boundaries and sea-level falls. Although regional peculiarities are superimposed on the appearance of global events in the Triassic synthem architecture, the successful global tracing suggests that planetary-scale mechanisms of synthem formation existed and that they were active in regions dominated by both marine and non-marine sedimentation.     

Magma evolution in the Purico ignimbrite complex, northern Chile: evidence for zoning of a dacitic magma by injection of rhyolitic melts following mafic recharge

The 1.3 Ma Purico complex is part of an extensive Neogene-Pleistocene ignimbrite province in the central Andes. Like most other silicic complexes in the province, Purico is dominated by monotonous intermediate ash-flow sheets and has volumetrically minor lava domes. The Purico ignimbrites (total volume 80-100 km³) are divided into a Lower Purico Ignimbrite (LPI) with two extensive flow units, LPI I and LPI II; and a smaller Upper Purico Ignimbrite (UPI) unit. Crystal-rich dacite is the dominant lithology in all the Purico ignimbrites and in the lava domes. It is essentially the only lithology present in the first LPI flow unit (LPI I) and in the Upper Purico Ignimbrite, but the LPI II flow unit is unusual for its compositional diversity. It constitutes a stratigraphic sequence with a basal fall-out deposit containing rhyolitic pumice (68-74 wt% SiO₂) overlain by ignimbrite with dominant crystal-rich dacitic pumice (64-66 wt% SiO₂). Rare andesitic and banded pumice (60-61 wt% SiO₂) are also present in the uppermost part of the flow unit. The different compositional groups of pumice in LPI II flow unit (rhyolite, andesite, dacite) have initial Nd and Sr isotopic compositions that are indistinguishable from each other and from the dominant dacitic pumice ()Nd=-6.7 to -7.2 and ⁸⁷Sr/⁸⁶Sr=0.7085-0.7090). However, two lines of evidence show that the andesite, dacite and rhyolite pumices do not represent a simple fractionation series. First, melt inclusions trapped in sequential growth zones of zoned plagioclase grains in the rhyolite record fractionation trends in the melt that diverge from those shown by dacite samples. Second, mineral equilibrium geothermometry reveals that dacites from all ignimbrite flow units and from the domes had relatively uniform and moderate pre-eruptive temperatures (780-800 °C), whereas the rhyolites and andesites yield consistently higher temperatures (850-950 °C). Hornblende geobarometry and pressure constraints from H₂O and CO₂ contents in melt inclusions indicate upper crustal (4-8 km) magma storage conditions. The petrologic evidence from the LPI II system thus indicates an anomalously zoned magma chamber with a rhyolitic cap that was hotter than, and chemically unrelated to, the underlying dacite. We suggest that the hotter rhyolite and andesite magmas are both related to an episode of replenishment in the dacitic Purico magma chamber. Rapid and effective crystal fractionation of the fresh andesite produced a hot rhyolitic melt whose low density and viscosity permitted ascent through the chamber without significant thermal and chemical equilibration with the resident dacite. Isotopic and compositional variations in the Purico system are typical of those seen throughout the Neogene ignimbrite complexes of the Central Andes. These characteristics were generated at moderate crustal depths (<30 km) by crustal melting, mixing and homogenization involving mantle-derived basalts. For the Purico system, assimilation of at least 30% mantle-derived material is required.     

The dark side of zircon: textural,age, oxygen isotopic and trace element evidence of fluid saturation in the subvolcanic reservoir of the Island Park-Mount Jackson Rhyolite,Yellowstone (USA)

The Island Park-Mount Jackson series in the Yellowstone volcanic field, Wyoming (USA), is a suite of rhyolitic domes and lavas that erupted between the caldera-forming eruptions of the Mesa Fall Tuff (1.3 Ma) and the Lava Creek Tuff (0.6 Ma). Combined zircon U/Pb geochronology, Raman spectroscopy, oxygen isotopic and trace elemental compositions document storage conditions of these magmas between consecutive supereruptions. Based on comparison with co-erupted melt compositions and textural criteria, four zircon compositional groups are identified that record different stages along a continuous magmatic evolution from trace element-poor rhyolite at high temperatures to extremely fractionated rhyolite where zircon trace elements are highly enriched (e.g., >?1000 ppm U). These latter zircon domains are dark in cathodoluminescence images and show broadened Raman peaks relative to near-endmember zircon, indicating that substitution of non-stoichiometric trace elements into zircon leads to distortion of the crystal lattice. Some of these zircon domains contain inclusions of U-Th-REE-phases, likely originating from coupled dissolution–reprecipitation of metastable trace element-rich zircon in the presence of a fluid phase. Rhyolite-MELTS simulations indicate that at the conditions required to produce the observed enrichment in trace elements, a fluid phase is likely present. These findings illustrate that zircons can be assembled from a variety of co-existing magmatic environments in the same magma reservoir, including near-solidus volatile-rich melts close to the magmatic–hydrothermal transition.     

Micro-scale S isotope studies of the Kharaelakh intrusion, Noril’sk region, Siberia: Constraints on the genesis of coexisting anhydrite and sulfide minerals

The coexistence of magmatic anhydrite and sulfide minerals in non-arc-related mafic magmas has only rarely been documented. Likewise the S isotope fractionation between sulfate and sulfide in mafic rocks has infrequently been measured. In the Kharaelakh intrusion associated with the world-famous Noril’sk ore district in Siberia coexisting magmatic anhydrite and sulfide minerals have been identified. Sulfur isotope compositions of the anhydrite-sulfide assemblages have been measured via both ion microprobe and conventional analyses to help elucidate the origin of the anhydrite-sulfide pairs. Magmatic anhydrite and chalcopyrite are characterized by δ³⁴S values between 18.8‰ and 22.8‰, and 9.3‰ and 13.2‰, respectfully. Coexisting anhydrite and chalcopyrite show Δ values that fall between 8.5‰ and 11.9‰. Anhydrite in the Kharaelakh intrusion is most readily explained by the assimilation of sulfate from country rocks; partial reduction to sulfide led to mixing between sulfate-derived sulfide and sulfide of mantle origin. The variable anhydrite and sulfide δ³⁴S values are a function of differing degrees of sulfate reduction, variable mixing of sulfate-derived and mantle sulfide, incomplete isotopic homogenization of the magma, and a lack of uniform attainment of isotopic equilibrium during subsolidus cooling. The δ³⁴S values of sulfide minerals have changed much less with cooling than have anhydrite values due in large part to the high sulfide/sulfate ratio. Variations in both sulfide and anhydrite δ³⁴S values indicate that isotopically distinct domains existed on a centimeter scale. Late stage hydrothermal anhydrite and pyrite also occur associated with Ca-rich hydrous alteration assemblages (e.g., thomsonite, prehnite, pectolite, epidote, xonotlite). δ³⁴S values of secondary hydrothermal anhydrite and pyrite determined by conventional analyses are in the same range as those of the magmatic minerals. Anhydrite-pyrite Δ values are in the 9.1-10.1‰ range, and are smaller than anticipated for the low temperatures indicated by the silicate alteration assemblages. The small Δ values are suggestive of either sulfate-sulfide isotopic disequilibrium or closure of the system to further exchange between ∼550 and 600 °C. Our results confirm the importance of the assimilation of externally derived sulfur in the generation of the elevated δ³⁴S values in the Kharaelakh intrusion, but highlight the sulfur isotopic variability that may occur in magmatic systems. In addition, our results confirm the need for more precise experimental determination of sulfate-sulfide sulfur isotope fractionation factors in high-T systems.     

Nekton use of salt marsh creeks in the upper Tejo estuary

The use of the Tejo estuary, Portugal, salt marsh creeks by nekton was examined based on sampling surveys with a fyke net from September 1998 until August 2001. From the 20 taxa (14 fish species, 5 decapod crustacean species, and 1 cephalopod species) identified in the studied creeks, 16 were regularly caught throughout the sampling period. The shrimpPalaemonetes varians was the most numerically abundant species in the creeks, while the biomass was dominated by the mulletLiza ramada. The nekton assemblage was mainly represented by marine-estuarine opportunist species, comprising 85% of the total. A high seasonality was detected on the species abundance patterns: the most abundant species (P. varians, Crangon crangon, L. ramada, Pomatoschistus microps, Syngnathus sp., andAnguilla anguilla) occurred throughout the sampling period,Sardina pilchardus, Dicentrarchus, labrax, andAtherina boyeri were particularly abundant in spring and summer, andEngraulis encrasicholus, Liza aurata, Gambusia holbrooki, Palaemon longirostris, andPalaemon serratus were most abundant in autumn and winter.L. ramada occurred in the tidal creeks in high numbers during neap tides, while the majority of the remaining taxa were most abundant during spring tides, suggesting a differential pattern of habitat use occording to species.