Combining in situ isotopic, trace element and textural analyses of quartz from four magmatic-hydrothermal ore deposits

This study couples in situ ¹⁶O, ¹⁷O and ¹⁸O isotope and in situ trace element analyses to investigate and characterize the geochemical and textural complexity of magmatic-hydrothermal quartz crystals. Euhedral quartz crystals contemporaneous with mineralization were obtained from four magmatic-hydrothermal ore deposits: El Indio Au–Ag–Cu deposit; Summitville Au–Ag–Cu deposit; North Parkes Cu–Au deposit and Kingsgate quartz-Mo–Bi–W deposit. The internal features of the crystals were imaged using cathodoluminescence and qualitative electron microprobe maps. Quantitative isotopic data were collected in situ using 157 nm laser ablation inductively coupled plasma mass spectrometry (for 40 trace elements in quartz) and sensitive high-resolution ion microprobe (for 3 isotopes in quartz). Imaging revealed fine oscillatory zoning, sector zoning, complex “macromosaic” textures and hidden xenocrystic cores. In situ oxygen isotope analyses revealed a δ¹⁸O range of up to 12.4 ± 0.3 ‰ in a single crystal—the largest isotopic range ever ascribed to oscillatory zonation in quartz. Some of these crystals contain a heavier δ¹⁸O signature than expected by existing models. While sector-zoned crystals exhibited strong trace element variations between faces, no evidence for anisotropic isotope fractionation was found. We found: (1) isotopic heterogeneity in hydrothermal quartz crystals is common and precludes provenance analysis (e.g., δD–δ¹⁸O) using bulk analytical techniques, (2) the trace element signature of quartz is not an effective pathfinder toward noble metal mineralization and (3) in three of the four samples, both textural and isotopic data indicate non-equilibrium deposition of quartz.     

Precessional forcing of tropical vegetation carbon storage

Since the Mid Pleistocene Revolution, which occurred about one million years ago, global temperatures have fluctuated with a quasi‐periodicity of ca. 100 ka. The pattern of past change in the extent of woodlands, and therefore by inference vegetation carbon storage, has been demonstrated to have a strong positive link with this global temperature change at high and mid latitudes. However, understanding of climate systems and ecosystem function indicates that the pattern of woodland change at low latitudes may follow a fundamentally different pattern. We present output from the intermediate complexity model GENIE‐1, comprising a single transient simulation over the last 800 ka and a 174‐member ensemble of 130 ka transient simulations over the last glacial cycle. These simulations suggest that while vegetation carbon storage in mid–high northern latitudes robustly follows the characteristic ca. 100 ka cycle, this signal is not a robust feature of tropical vegetation, which is subject to stronger direct forcing by the precessional (21 ka) orbital cycle (albeit with a highly uncertain response). We conclude that the correlation of palaeoenvironmental records from low latitudes with global temperature change must be done with caution. Copyright © 2011 John Wiley & Sons, Ltd.     

High‐Precision,High‐Accuracy Oxygen Isotope Measurements of Zircon Reference Materials with the SHRIMP‐SI

Analytical protocols for SHRIMP‐SI oxygen isotope analysis (δ¹⁸O) of a suite of zircon reference materials (RMs) are presented. Data reduction involved a robust estimate of uncertainties associated with the individual spot as well as for groups where the spot data are combined. The repeatability of δ¹⁸O measurements is dependent on both the analytical conditions and the choice of the primary reference material. Under optimised conditions, repeatability was often better than 0.4‰ (2s) allowing sample uncertainties to be obtained to better than 0.2‰ (at 95% confidence limit). Single spot uncertainty combined the within‐spot precision with the scatter associated with repeated measurements of the primary zircon reference material during a measurement session. The uncertainty for individual spots measured under optimised conditions was between 0.3 and 0.4‰ (at 95% confidence). The analytical protocols described were used to assess a variety of zircon RMs that have been used for geochronology and for which laser fluorination oxygen isotope data are available (Temora 2, FC1, R33, QGNG and Ple?ovice), as well as zircons that have been used as RMs for trace element or other types of determination (Mud Tank, Monastery, 91500, AS57, AS3, KIM‐5, OG1, SL13, CZ3 and several other Sri Lankan zircons). Repeated analyses over nine sessions and seven different mounts show agreement within analytical uncertainty for Temora 2, FC1, R33, QGNG, Ple?ovice and 91500, when normalised to Mud Tank. For existing ion microprobe mounts with these materials, an appropriate δ¹⁸O can be determined. However, care should be taken when using zircons from the Duluth Complex (i.e., FC1, AS57 and AS3) as reference materials as our data indicated an excess scatter on δ¹⁸O values associated with low‐U zircon grains.     

Predicted changes in fire weather suggest increases in lightning fire initiation and future area burned in the mixedwood boreal forest

Forecasting future fire activity as a function of climate change is a step towards understanding the future state of the western mixedwood boreal ecosystem. We developed five annual weather indices based on the Daily Severity Rating (DSR) of the Canadian Forest Fire Weather Index System and estimated their relationship with annual, empirical counts of lightning fire initiation for 588 landscapes in the mixedwood boreal forest in central-eastern Alberta, Canada from data collected between 1983 and 2001 using zero-inflated negative binomial regression models. Two indices contributed to a parsimonious model of initiation; these were Seasonal Severity Rating (SSR), and DSR-sequence count. We used parameter estimates from this model to predict lightning fire initiation under weather conditions predicted in 1 × CO₂ (1975–1985), 2 × CO₂ (2040–2049) and 3 × CO₂ (2080–2089) conditions simulated by the Canadian Regional Climate Model (CRCM). We combined predicted initiation rates for these conditions with existing empirical estimates of the number of fire initiations that grow to be large fires (fire escapes) and the fire size distribution for the region, to predict the annual area burned by lightning-caused fires in each of the three climate conditions. We illustrated a 1.5-fold and 1.8-fold increase of lightning fire initiation by 2040–2049 and 2080–2089 relative to 1975–1985 conditions due to changes in fire weather predicted by the CRCM; these increases were calculated independent of changes in lightning activity. Our simulations suggested that weather-mediated increases in initiation frequency could correspond to a substantial increase in future area burned with 1.9-fold and 2.6-fold increases in area burned in 2040–2049 and 2080–2089 relative to 1975–1985 conditions, respectively. We did not include any biotic effects in these estimates, though future patterns of initiation and fire growth will be regulated not only by weather, but also by vegetation and fire management.     

Effects of pipe outlet blocking on hydrological functioning in a degraded blanket peatland

Peatland restoration practitioners are keen to understand the role of drainage via natural soil pipes, especially where erosion has released large quantities of fluvial carbon in stream waters. However, little is known about pipe-to-stream connectivity and whether blocking methods used to impede flow in open ditch networks and gullies also work on pipe networks. Two streams in a heavily degraded blanket bog (southern Pennines, UK) were used to assess whether impeding drainage from pipe networks alters the streamflow responses to storm events, and how such intervention affects the hydrological functioning of the pipe network and the surrounding peat. Pipeflow was impeded in half of the pipe outlets in one stream, either by inserting a plug-like structure in the pipe-end or by the insertion of a vertical screen at the pipe outlet perpendicular to the direction of the predicted pipe course. Statistical response variable η² showed the overall effects of pipe outlet blocking on stream responses were small with η² = 0.022 for total storm runoff, η² = 0.097 for peak discharge, η² = 0.014 for peak lag, and η² = 0.207 for response index. Both trialled blocking methods either led to new pipe outlets appearing or seepage occurring around blocks within 90 days of blocking. Discharge from four individual pipe outlets was monitored for 17 months before blocking and contributed 11.3% of streamflow. Pipe outlets on streambanks with headward retreat produced significantly larger peak flows and storm contributions to streamflow compared to pipe outlets that issued onto straight streambank sections. We found a distinctive distance-decay effect of the water table around pipe outlets, with deeper water tables around pipe outlets that issued onto straight streambanks sections. We suggest that impeding pipeflow at pipe outlets would exacerbate pipe development in the gully edge zone, and propose that future pipe blocking efforts in peatlands prioritize increasing the residence time of pipe water by forming surface storage higher up the pipe network.     

Patterns and drivers of peat topographic changes determined from Structure-from-Motion photogrammetry at field plot and laboratory scales

Little is known about the spatial and temporal variability of peat erosion nor some of its topographic and weather-related drivers. We present field and laboratory observations of peat erosion using Structure-from-Motion (SfM) photogrammetry. Over a 12 month period, 11 repeated SfM surveys were conducted on four geomorphological sites of 18–28 m² (peat hagg, gully wall, riparian area and gully head) in a blanket peatland in northern England. A net topographic change of –14 to +30 mm yr^–1 for the four sites was observed during the whole monitoring period. Cold conditions in the winter of 2016 resulted in highly variable volume change (net surface topographic rise first and lowering afterwards) via freeze–thaw processes. Long periods of dry conditions in the summer of 2017 led to desiccation and drying and cracking of the peat surface and a corresponding surface lowering. Topographic changes were mainly observed over short-term intervals when intense rainfall, flow wash, needle-ice production or surface desiccation was observed. In the laboratory, we applied rainfall simulations on peat blocks and compared the peat losses quantified by traditional sediment flux measurements with SfM derived topographic data. The magnitude of topographic change determined by SfM (mean value: 0.7 mm, SD: 4.3 mm) was very different to the areal average determined by the sediment yield from the blocks (mean value: –0.1 mm, SD: 0.1 mm). Topographic controls on spatial patterns of topographic change were illustrated from both field and laboratory surveys. Roughness was positively correlated to positive topographic change and was negatively correlated to negative topographic change at field plot scale and laboratory macroscale. Overall, the importance of event-scale change and the direct relationship between surface roughness and the rate of topographic change are important characteristics which we suggest are generalizable to other environments. © 2018 John Wiley & Sons, Ltd.     

A comparison of stream water temperature regimes from open and afforested moorland,Yorkshire Dales,northern England

Despite the known importance of water temperature for river ecosystems, the thermal regime of streams and rivers can be heavily modified by afforestation. Although the nature of the heat budget affecting streams in forested catchments shows high variability in space and time, most of the studies of stream temperature response to afforestation have lacked replication among streams. This study examined the impacts of coniferous forest plantations on stream water temperature at six sites (three forested and three open moorland) in the Yorkshire Dales, northern England. Our aim was to test the hypothesis that afforestation would alter the thermal regime of streams, leading to reduced year‐round thermal variability, and cooler summer/warmer winter water temperatures, relative to streams flowing across open moorland. Data collected from April 2007 to March 2009 showed similar thermal dynamics among all six streams over the study period, although variability in forested streams was markedly lower as expected. Mean and maximum daily water temperatures were significantly higher in open moorland streams for much of the year but while some forested streams were warmer than individual moorland streams during winter months (November to February), there was considerable overlap in water temperature between moorland and forest streams. Most stream temperature records showed evidence of low/no winter flow and freezing. These results contrast with many previous studies that have reported warmer temperatures in forested versus open moorland streams during winter, a finding that most likely reflects site‐specific hydrological, geomorphological and climatological influences on water temperature in addition to afforestation. This study demonstrates the need for replication of hydrological monitoring when examining the effects of basin‐scale management practices and provides further evidence for changes in stream thermal regime following afforestation, a practice that is likely to increase in future due to growing demands for increased forest cover in the UK uplands. Copyright © 2010 John Wiley & Sons, Ltd.     

The dynamics of natural pipe hydrological behaviour in blanket peat

Natural soil pipes are found in peatlands, but little is known about their hydrological role. This paper presents the most complete set of pipe discharge data to date from a deep blanket peatland in Northern England. In a 17.4‐ha catchment, we identified 24 perennially flowing and 60 ephemerally flowing pipe outlets. Eight pipe outlets along with the catchment outlet were continuously gauged over an 18‐month period. The pipes in the catchment were estimated to produce around 13.7% of annual streamflow, with individual pipes often producing large peak flows (maximum peak of 3.8 l s^?1). Almost all pipes, whether ephemerally or perennially flowing, shallow or deep (outlets > 1 m below the peat surface), showed increased discharge within a mean of 3 h after rainfall commencement and were dominated by stormflow, indicating good connectivity between the peatland surface and the pipes. However, almost all pipes had a longer period between the hydrograph peak and the return to base flow compared with the stream (mean of 23.9 h for pipes, 19.7 h for stream). As a result, the proportion of streamflow produced by the pipes at any given time increased at low flows and formed the most important component of stream discharge for the lowest 10% of flows. Thus, a small number of perennially flowing pipes became more important to the stream system under low‐flow conditions and probably received water via matrix flow during periods between storms. Given the importance of pipes to streamflow in blanket peatlands, further research is required into their wider role in influencing stream water chemistry, water temperature and fluvial carbon fluxes, as well as their role in altering local hydrochemical cycling within the peat mass itself. Enhanced piping within peatlands caused by environmental change may lead to changes in the streamflow regime with larger low flows and more prolonged drainage of the peat. Copyright © 2012 John Wiley & Sons, Ltd.     

Stochastic Structural Modeling

A consistent stochastic model for faults and horizons is described. The faults are represented as a parametric invertible deformation operator. The faults may truncate each other. The horizons are modeled as correlated Gaussian fields and are represented in a grid. Petrophysical variables may be modeled in a reservoir before faulting in order to describe the juxtaposition effect of the faulting. It is possible to condition the realization on petrophysics, horizons, and fault plane observations in wells in addition to seismic data. The transmissibility in the fault plane may also be included in the model. Four different methods to integrate the fault and horizon models in a common model is described. The method is illustrated on an example from a real petroleum field with 18 interpreted faults that are handled stochastically.