Origins of cold-wet-oxidizing to hot-dry-reducing rhyolite magma cycles and distribution in the Taupo Volcanic Zone, New Zealand

Slab-derived aqueous fluid components (Ba, Cl) correlate well with oxygen fugacity, and other well-defined characteristics of silicic magmas in the Taupo Volcanic Zone (TVZ) between a cold-wet-oxidizing magma type [R1: amphibole ± biotite; high Sr, low Zr and FeO*/MgO, depleted middle rare-earth elements (MREE)] and a hot-dry-reducing magma type (R2: orthopyroxene ± clinopyroxene; low Sr, high Zr, and FeO*/MgO, less depleted MREE). Oxygen fugacity was obtained from analysis of Fe–Ti oxides and ranges between −0.04 and +2.1 log units (ΔQFM, where QFM = quartz + fayalite + magnetite buffer) and is positively correlated with the bulk-rock Ba/La ratio, indicating that slab-derived fluid is the oxidizing agent in the rhyolites. Chlorine contents in hornblende also correlate with the bulk-rock Ba/La ratio. Hence, high-fluid flux typically correlates with the R1 and low-fluid flux with R2 rhyolite magma types. A geochemical evolution and distribution can be tracked in time and space throughout the central region of the TVZ from 550 ka to present and has revealed two distinct magmatic cycles that vary in length. The first cycle included widespread R1 type magmatism across the central TVZ beginning ca. 550 ka and was directly associated with previously unreported dome-building and ignimbrite-forming volcanism, and led to a voluminous (>3,000 km³) ignimbrite ‘flare-up’ between ca. 340 and 240 ka. The second cycle began roughly 180 ka, erupting ca. 800 km³ of magma, and continues to the present. The duration, rate, and composition of magma production within these cycles appears to be governed by the flux of fluid released from the subducting slab, while the distribution of magmas may be governed more by extension along the central rift axis. Shorter cycles have also been identified and are unrelated to subduction processes, but occur following large, caldera-forming events.     

Timing and thermochemical constraints on multi-element mineralisation at the Nori/RA Cu–Mo–U prospect, Great Bear magmatic zone, Northwest Territories, Canada

The timing of Cu–Mo–U mineralisation at the Nori/RA prospect in the Paleoproterozoic Great Bear magmatic zone has been investigated using Re–Os molybdenite and ⁴⁰Ar–³⁹Ar biotite geochronology. The Re–Os molybdenite ages presented are the first robust sulphide mineralisation ages derived from the Great Bear magmatic zone. Cu–Mo–U mineralisation is hosted in early to syn-deformational hydrothermal veins consisting of quartz and K-feldspar or more commonly tourmaline-biotite-quartz-K-feldspar, with associated wall-rock alteration assemblages being predominantly biotite. Sulphide and oxide minerals consist of chalcopyrite, molybdenite and uraninite with lesser pyrite and magnetite. Elevated light rare earth elements and tungsten concentrations associated with the Cu–Mo–U mineralisation have also been reported at the prospect by previous workers. Molybdenite and uraninite occur intimately in dravitic tourmaline growth zones and at grain margins, attesting to their syngenetic nature (with respect to hydrothermal veining). Two molybdenite separates yield Re–Os model ages of 1,874.4 ± 8.7 (2σ) and 1,872.4 ± 8.8 Ma (2σ) with a weighted average model age of 1,873.4 ± 6.1 Ma (2σ). Laser step heating of biotite from the marginal alteration of the wall-rock adjacent to the veins yields a ⁴⁰Ar–³⁹Ar maximum cooling age of 1,875 ± 8 Ma (MSWD = 3.8; 2σ), indistinguishable from the Re–Os molybdenite model age and a previously dated ‘syn-tectonic’ aplitic dyke in the region. Dravitic tourmaline hosts abundant primary liquid–vapour–solid-bearing fluid inclusions. Analytical results indicate liquid–vapour homogenisation at >260°C constraining the minimum temperature of mineralisation. The solids, which are possibly trapped, did not homogenise with the liquid–vapour by 400°C. Salinities in the inclusions are variable. Raman spectra identify that at least some of the solids are calcite and anhydrite. Raman spectra also confirm the vapour phases contain some CO₂; whereas clathrates or CH₄ was not observed or detected. Quartz grains only host secondary fluid inclusions, which fluoresce under ultraviolet light, indicating trapped hydrocarbons. We speculate that these resulted from Phanerozoic fluid circulation through the Proterozoic basement. The collective interpretation of the age, hydrothermal character and associated metals, high temperature and variable salinity suggests that the Nori/RA Cu–Mo–U mineralisation can be linked with the earliest stages of plutonism in the Great Bear magmatic zone. From a regional perspective, the mineralisation may pre-date the extensive multi-element mineralisation now recognised as part of the iron oxide copper–gold (IOCG) spectrum of deposits. As IOCG provinces generally contain a variety of mineralisation styles, we interpret this as the earliest phase of the extensive mineralising system.     

Assessment of undiscovered nickel sulphide resources,Kalgoorlie Terrane,Western Australia: Part 1. Deposit and endowment density models

The use of mineral deposit density regression models to estimate the number of undiscovered deposits is gaining acceptance in mineral resources assessments. The deposit density regression models currently in use are based on well-established power law relationships between deposit density (deposits/km²) and the areal extent of the host rocks in well explored regions (control areas) worldwide. Although these generalized or global deposit density models can generate guideline estimates that are useful at the terrane scale, locally-derived terrane-based deposit density regression models may potentially yield more relevant estimates at the terrane scale. Using 12 selected komatiite-defined control areas in the Kalgoorlie Terrane, Western Australia, we found that the size (km²) of the control areas had power law relationships with (i) nickel sulphide deposit density, and (ii) nickel endowment density (nickel metal/km²). Regression analyses showed that both power law relationships are statistically significant at the 5% level. This suggests that nickel sulphide deposit and endowment density models could be used to estimate the number of undiscovered nickel sulphide deposits and amount of nickel metal endowment in less explored komatiites in the Kalgoorlie Terrane. This study shows that global geological relationships can be viably downscaled onto local geological terranes thereby supporting the hypothesis that the processes of mineral deposit formation and preservation are scale-independent and self-similar.     

Multi-criteria, multi-objective and uncertainty analysis for agro-energy spatial modelling

The worldwide increase in the use of biomass as a Renewable Energy Source raises the issue of introducing crops dedicated to energy production into rural landscapes. The purpose of this paper is to set-up a GIS based multi-criteria approach to assess a range of possibilities for perennial energy crops conversion. The presented method was implemented at the regional level in the Yorkshire and the Humber Region in Northern UK. The first phase of the study aims to set-up a land capability model for the specific purpose of assessing the potential of different typologies of perennial energy crops, on the basis of specific pedo-climatic and topographic factors. The model output illustrates a range of potentials for energy crop conversion that can be explored in the given landscape. In the second phase a uncertainty analysis of the land capability model was performed through a simulation approach in order to interpret the influence of assumptions and uncertainty on input data and model parameters. The last phase of the study allows allocating the energy crop conversion area according to specific environmental constraints, nature protection targets, food production priorities and land capability values. The land allocation output gives a rather restrictive energy crop penetration scenario, where more than half of the conversion area is allocated to cropping systems with low land degradation potential. This scenario represents a preliminary regional analysis of the energy crop potential in terms of theoretically available conversion areas. The final results also show that the areas with highest environmental risks correspond to the areas with both the lowest suitability for energy crop cultivation and the highest model uncertainty.     

Rapid topographic change measured by high-resolution satellite radar at Soufriere Hills Volcano,Montserrat, 2008–2010

High-resolution satellite radar observations of erupting volcanoes can yield valuable information on rapidly changing deposits and geomorphology. Using the TerraSAR-X (TSX) radar with a spatial resolution of about 2 m and a repeat interval of 11 days, we show how a variety of techniques were used to record some of the eruptive history of the Soufriere Hills Volcano, Montserrat between July 2008 and February 2010. After a 15-month pause in lava dome growth, a vulcanian explosion occurred on 28 July 2008 from a vent that was hidden by dense cloud. We were able to show the civil authorities using TSX difference images of surface roughness change that this explosion had not disrupted the dome sufficiently to warrant continuation of a previous, precautionary evacuation. Change difference images also proved to be valuable in mapping new pyroclastic flow deposits: the valley-occupying block-and-ash component tended to increase backscatter and the marginal surge deposits to reduce it, with the pattern reversing after the event due to erosion and deposition. By comparing east- and west-looking images acquired 12 h apart, the deposition of some individual pyroclastic flows can be inferred from change differences. Some of the narrow upper sections of valleys draining the volcano received many tens of metres of rockfall and pyroclastic flow deposits over periods of a few weeks. By measuring the changing radar shadows cast by these valleys in TSX images the changing depth of infill by deposits could be estimated. In addition to using the amplitude data from the radar images we also used their phase information within the InSAR technique to calculate the topography during a period of no surface activity. This enabled areas of transient topography, crucial for directing future flows, to be captured.     

Inherent sea ice predictability in the rapidly changing Arctic environment of the Community Climate System Model, version 3

Seasonal predictions of Arctic sea ice have typically been based on statistical regression models or on results from ensemble ice model forecasts driven by historical atmospheric forcing. However, in the rapidly changing Arctic environment, the predictability characteristics of summer ice cover could undergo important transformations. Here global coupled climate model simulations are used to assess the inherent predictability of Arctic sea ice conditions on seasonal to interannual timescales within the Community Climate System Model, version 3. The role of preconditioning of the ice cover versus intrinsic variations in determining sea ice conditions is examined using ensemble experiments initialized in January with identical ice?Cocean?Cterrestrial conditions. Assessing the divergence among the ensemble members reveals that sea ice area exhibits potential predictability during the first summer and for winter conditions after a year. The ice area exhibits little potential predictability during the spring transition season. Comparing experiments initialized with different mean ice conditions indicates that ice area in a thicker sea ice regime generally exhibits higher potential predictability for a longer period of time. In a thinner sea ice regime, winter ice conditions provide little ice area predictive capability after approximately 1?year. In all regimes, ice thickness has high potential predictability for at least 2?years.     

Assessing water quality in Marine Protected Areas from Southern California, USA

Despite the regulatory mandate to maintain “natural water quality”, there are ?271 storm drain discharges that potentially threaten the 14 designated marine water quality protected areas in Southern California called Areas of Special Biological Significance (ASBS). After sampling 35 site-events, the geomean concentrations of total suspended solids, nutrients, total and dissolved trace metals, and polycyclic aromatic hydrocarbons in the ocean following storm events were similar between reference drainages and ASBS discharge sites. Concentrations of chlorinated hydrocarbons were nondetectable and no post-storm sample exhibited significant toxicity to the endemic purple sea urchin (Strongylocentrotus purpuratus) near ASBS discharge sites. A reference-based threshold was developed and, despite the similarities in average concentrations, there were some individual ASBS discharge sites that were greater than reference background. Cumulatively across all ASBS, the constituents that were most frequently greater than the reference-based threshold were nutrients and general constituents, followed by dissolved and total trace metals.     

Currents in a small channel on a sandy tidal flat

Channels affect drainage and bed stresses on tidal flats. Here, near-bottom currents observed on a sandy tidal flat are compared with those observed 35 m away inside a shallow (≈0.3 m deep) channel. For water depths between 0.5 and 2.5 m (when both current meters are submerged), current speeds 0.13 m above the bed on the flat are about 30% greater than those observed 0.13 m above the bed in the channel, and are approximately equal to those observed 0.58 m above the channel bed (0.26 m above the flat elevation). Flow directions on the flat are similar to those in the channel. For flows directed across the channel axis, the ratio of speeds increases from about 1.3 to about 2.2 with increasing water depth. The corresponding ratio of the vertical velocity variances (a proxy for turbulence) decreases from about 1.5 to about 0.2, suggesting that the turbulence near the bed of the channel is greater than that near the bed of the flat for water depths greater than about 1.0 m. Drag coefficients estimated with the vertical velocity variance are approximately 70% larger in the channel than over the visually smoother flat, consistent with prior studies suggesting that channels may increase tidal-flat roughness. For flows directed along the channel axis (in the cross-flat direction), the ratio of speeds (1.2) is similar to the ratio predicted by a cross-flat momentum (along-channel) balance.     

Post-glacial regional climate variability along the East Antarctic coastal margin—Evidence from shallow marine and coastal terrestrial records

We review the post-glacial climate variability along the East Antarctic coastline using terrestrial and shallow marine geological records and compare these reconstructions with data from elsewhere. Nearly all East Antarctic records show a near-synchronous Early Holocene climate optimum (11.5–9 ka BP), coinciding with the deglaciation of currently ice-free regions and the optimum recorded in Antarctic ice and marine sediment cores. Shallow marine and coastal terrestrial climate anomalies appear to be out of phase after the Early Holocene warm period, and show complex regional patterns, but an overall trend of cooling in the terrestrial records. A Mid to Late Holocene warm period is present in many East Antarctic lake and shallow coastal marine records. Although there are some differences in the regional timing of this warm period, it typically occurs somewhere between 4.7 and 1 ka BP, which overlaps with a similar optimum found in Antarctic Peninsula terrestrial records. The differences in the timing of these sometimes abrupt warm events in different records and regions points to a number of mechanisms that we have yet to identify. Nearly all records show a neoglacial cooling from 2 ka BP onwards. There is no evidence along the East Antarctic coastline for an equivalent to the Northern Hemisphere Medieval Warm Period and there is only weak circumstantial evidence in a few places for a cool event crudely equivalent in time to the Northern Hemisphere's Little Ice Age. There is a need for well-dated, high resolution climate records in coastal East Antarctica and particularly in Terre Adélie, Dronning Maud Land and Enderby Land to fully understand the regional climate anomalies, the disparity between marine and terrestrial records, and to determine the significance of the heterogeneous temperature trends being measured in the Antarctic today.     

Short-term variation in thermal stratification complicates estimation of lake metabolism

Previous studies have used sondes to measure diel changes in dissolved oxygen and thereby estimate gross primary production (GPP), ecosystem respiration (R), and net ecosystem production (NEP). Most of these studies estimate rates for the surface layer and require knowing the depth of the mixed layer (Z_mix), which is usually determined from discrete daily or weekly temperature profiles. However, Z_mix is dynamic, as the thermal structure of lakes may change at scales of minutes rather than days or weeks. We studied two thermally stratified lakes that exhibited intermittent microstratification in the mixed layer. We combined sonde-based estimates of metabolism with high-frequency measurements of stratification using thermistor chains to determine how the short-term dynamics of stratification affect metabolic rates. We calculated estimates of metabolism using time series of Z_mix measured at seasonal, weekly, daily, and 5-min intervals. Areal rates of GPP and R were up to 24 and 29% less, respectively, using the 5-min measurements of Z_mix rather than weekly Z_mix, while NEP was not significantly different. These reduced areal rates are mostly the consequence of the reduction in the depth of the mixed layer. Microstratification occurred frequently in both lakes and affected volumetric rates in one lake where R was significantly lower, NEP was significantly higher, and GPP was marginally lower compared to days without microstratification. Hence, microstratification not only affects the depth of the mixed layer, but also alters the processes that influence photosynthesis and respiration. Future studies should consider microstratification and possibly employ multiple sondes with thermistor chains that enable integrating metabolic rates to a specific depth, rather than assuming a stable upper mixed layer as the basis for calculations.