An archaean granulite-grade tonalite-trondhjemite-granite suite from scourie,NW Scotland: Geochemistry and origin

The Lewisian complex of the Scourie-Badcall area is composed predominantly of banded tonalitic gneiss which intrudes layered gabbro-ultramafic complexes. Intrusive into both gabbro and tonalitic gneiss are homogeneous acid sheets which are trondhjemitic to granitic in composition. All rocks were subjected to granulite facies metamorphism. Smooth continuous trends on chemical variation diagrams suggest that the evolution of these rocks was dominated by fractional crystallisation. A scheme is proposed whereby a tonalitic melt was parental to trondhjemite and granite. Variation within tonalites was a function of the fractional crystallisation of hornblende and plagioclase, and trondhjemite was derived from tonalite by the fractional crystallisation of hornblende and/or plagioclase. Granite and granodiorite represent residual liquids which evolved along the quartz-feldspar cotectic surface; they were derived by the fractional crystallisation of plagioclase from a trondhjemite liquid. Some trondhjemitic sheets are quartz-plagioclase residues from which a granitic melt was removed. The associated gabbros and ultramafic rocks are not directly related to the proposed fractional crystallisation scheme and are not crystal residues removed from the tonalitic melt. Tonalites were probably derived from a basaltic source by partial melting or fractional crystallisation with either hornblende and/or garnet as residual phases.     

The breakdown of potassium feldspar at high water pressures

The equilibrium position of the reaction between sanidine and water to form “sanidine hydrate” has been determined by reversal experiments on well characterised synthetic starting materials in a piston cylinder apparatus. The reaction was found to lie between four reversed brackets of 2.35 and 2.50 GPa at 450 °C, 2.40 and 2.59 GPa at 550 °C, 2.67 and 2.74 GPa at 650 °C, and 2.70 and 2.72 GPa at 680 °C. Infrared spectroscopy showed that the dominant water species in sanidine hydrate was structural H₂O. The minimum quantity of this structural H₂O, measured by thermogravimetric analysis, varied between 4.42 and 5.85 wt% over the pressure range of 2.7 to 3.2 GPa and the temperature range of 450 to 680 °C. Systematic variation in water content with pressure and temperature was not clearly established. The maximum value was below 6.07 wt%, the equivalent of 1 molecule of H₂O per formula unit. The water could be removed entirely by heating at atmospheric pressure to produce a metastable, anhydrous, hexagonal KAlSi₃O₈ phase (“hexasanidine”) implying that the structural H₂O content of sanidine hydrate can vary. The unit cell parameters for sanidine hydrate, measured by powder X-ray diffraction, were a = 0.53366 (±0.00022) nm and c = 0.77141 (±0.00052) nm, and those for hexasanidine were a = 0.52893 (±0.00016) nm and c = 0.78185 (±0.00036) nm. The behaviour and properties of sanidine hydrate appear to be analogous to those of the hydrate phase cymrite in the equivalent barium system. The occurrence of sanidine hydrate in the Earth would be limited to high pressure but very low temperature conditions and hence it could be a potential reservoir for water in cold subduction zones. However, sanidine hydrate would probably be constrained to granitic rock compositions at these pressures and temperatures. Received: 6 May 1997 / Accepted: 2 October 1997     

Integrated provenance analysis of a convergent retroarc foreland system: U–Pb ages,heavy minerals,Nd isotopes,and sandstone compositions of the Middle Magdalena Valley basin,northern Andes,Colombia

Sediment provenance analysis remains a powerful method for testing hypotheses on the temporal and spatial evolution of uplifted source regions, but issues such as recycling, nonunique sources, and pre- and post-depositional modifications may complicate interpretation of results from individual provenance techniques. Convergent retroarc systems commonly contain sediment sources that are sufficiently diverse (continental magmatic arc, fold–thrust belt, and stable craton) to enable explicit provenance assessments. In this paper, we combine detrital zircon U–Pb geochronology, heavy mineral identification, Nd isotopic analyses, conventional sandstone petrography, and paleocurrent measurements to reconstruct the clastic provenance history of a long-lived sedimentary basin now exposed in an intermontane zone of the northern Andean hinterland of Colombia. The Middle Magdalena Valley basin, situated between the Central Cordillera and Eastern Cordillera, contains a 5–10 km-thick succession of Upper Cretaceous to Quaternary fill. The integrated techniques show a pronounced change in provenance during the Paleocene transition from the lower to upper Lisama Formation. We interpret this as a shift from an eastern cratonic source to a western Andean source composed of magmatic-arc rocks uplifted during initial shortening of the Central Cordillera. The appearance of detrital chloritoid and a shift to more negative ε_Nd(t=0) values in middle Eocene strata of the middle La Paz Formation are attributed to shortening-related exhumation of a continental basement block (La Cira–Infantas paleohigh), now buried, along the axis of the Magdalena Valley. The diverse provenance proxies also show distinct changes during middle to late Eocene deposition of the Esmeraldas Formation that likely reflect initial rock uplift and exhumation of the fold–thrust belt defining the Eastern Cordillera. Upsection, detrital zircon U–Pb ages and heavy mineral assemblages for Oligocene and younger clastic deposits indicate that the Mesozoic sedimentary cover of the Eastern Cordillera was recycled during continued Cenozoic shortening. Our multidisciplinary provenance study refines the tectonic history of the Colombian Andes and demonstrates that uncertainties related to sediment recycling, nonunique sources, source heterogeneity, and climate in interpreting provenance data can be minimized via an integrated approach.     

Seismic interferometry and ambient noise tomography in the British Isles

Traditional methods of imaging the Earth's subsurface using seismic waves require an identifiable, impulsive source of seismic energy, for example an earthquake or explosive source. Naturally occurring, ambient seismic waves form an ever-present source of energy that is conventionally regarded as unusable since it is not impulsive. As such it is generally removed from seismic data and subsequent analysis. A new method known as seismic interferometry can be used to extract useful information about the Earth's subsurface from the ambient noise wavefield. Consequently, seismic interferometry is an important new tool for exploring areas which are otherwise seismically quiescent, such as the British Isles in which there are relatively few strong earthquakes. One of the possible applications of seismic interferometry is ambient noise tomography (ANT). ANT is a way of using interferometry to image subsurface seismic velocity variations using seismic (surface) waves extracted from the background ambient vibrations of the Earth. To date, ANT has been used successfully to image the Earth's crust and upper-mantle on regional and continental scales in many locations and has the power to resolve major geological features such as sedimentary basins and igneous and metamorphic cores. Here we provide a review of seismic interferometry and ANT, and show that the seismic interferometry method works well within the British Isles. We illustrate the usefulness of the method in seismically quiescent areas by presenting the first surface wave group velocity maps of the Scottish Highlands using only ambient seismic noise. These maps show low velocity anomalies in sedimentary basins such as the Moray Firth, and high velocity anomalies in igneous and metamorphic centres such as the Lewisian complex. They also suggest that the Moho shallows from south to north across Scotland which agrees with previous geophysical studies in the region.     

Initial Characterisation of SK01 Sanidine: A Candidate 40Ar/39Ar Reference Material

The precision and accuracy of ⁴⁰Ar/³⁹Ar dates are ultimately linked to co-irradiated reference materials of known age. Here we provide new data from the SK01 sanidine, which was analysed in three different laboratories to evaluate it as a ⁴⁰Ar/³⁹Ar reference material. Aliquots of 5 mg, incrementally heated in two laboratories, yielded indistinguishable results with a weighted mean age of 27.58 ± 0.06 Ma (95% confidence level). Single-crystal step heating and single-crystal total fusion analyses of SK01 sanidine were undertaken in the third laboratory to further test the intracrystalline homogeneity. For the seven step-heating analyses, six crystals yielded nearly concordant age spectra with ⁴⁰Ar/³⁹Ar ages ranging from 26.853 ± 0.094 Ma to 26.963 ± 0.067 Ma, whereas one crystal gave an older age of 27.774 ± 0.071 Ma with a slightly discordant age spectrum. Twenty-three single-crystal total fusion analyses yielded ⁴⁰Ar/³⁹Ar ages ranging from 27.070 ± 0.108 Ma to 27.736 ± 0.062 Ma with a dispersion of ~ 3.8%. The older ages from single-crystal total fusion dates are interpreted to reflect an inherited or excess argon component in some crystals. This is an initial characterisation of the SK01 sanidine, and additional work needs to be conducted to further evaluate the age dispersion so that it can be utilised as a ⁴⁰Ar/³⁹Ar reference material.     

Uncertainty and sensitivity in a bank stability model: implications for estimating phosphorus loading

Eutrophication of aquatic ecosystems is one of the most pressing water quality concerns in the United States and around the world. Bank erosion has been largely overlooked as a source of nutrient loading, despite field studies demonstrating that this source can account for the majority of the total phosphorus load in a watershed. Substantial effort has been made to develop mechanistic models to predict bank erosion and instability in stream systems; however, these models do not account for inherent natural variability in input values. To quantify the impacts of this omission, uncertainty and sensitivity analyses were performed on the Bank Stability and Toe Erosion Model (BSTEM), a mechanistic model developed by the US Department of Agriculture – Agricultural Research Service (USDA‐ARS) that simulates both mass wasting and fluvial erosion of streambanks. Generally, bank height, soil cohesion, and plant species were found to be most influential in determining stability of clay (cohesive) banks. In addition to these three inputs, groundwater elevation, stream stage, and bank angle were also identified as important in sand (non‐cohesive) banks. Slope and bank height are the dominant variables in fluvial erosion modeling, while erodibility and critical shear stress had low sensitivity indices; however, these indices do not reflect the importance of critical shear stress in determining the timing of erosion events. These results identify important variables that should be the focus of data collection efforts while also indicating which less influential variables may be set to assumed values. In addition, a probabilistic Monte‐Carlo modeling approach was applied to data from a watershed‐scale sediment and phosphorus loading study on the Missisquoi River, Vermont to quantify uncertainty associated with these published results. While our estimates aligned well with previous deterministic modeling results, the uncertainty associated with these predictions suggests that they should be considered order of magnitude estimates only. Copyright © 2016 John Wiley & Sons, Ltd.     

Near‐surface hydraulic conductivity of northern hemisphere glaciers

The hydrology of near‐surface glacier ice remains a neglected aspect of glacier hydrology despite its role in modulating meltwater delivery to downstream environments. To elucidate the hydrological characteristics of this near‐surface glacial weathering crust, we describe the design and operation of a capacitance‐based piezometer that enables rapid, economical deployment across multiple sites and provides an accurate, high‐resolution record of near‐surface water‐level fluctuations. Piezometers were employed at 10 northern hemisphere glaciers, and through the application of standard bail–recharge techniques, we derive hydraulic conductivity (K) values from 0.003 to 3.519 m day^?1, with a mean of 0.185 ± 0.019 m day^?1. These results are comparable to those obtained in other discrete studies of glacier near‐surface ice, and for firn, and indicate that the weathering crust represents a hydrologically inefficient aquifer. Hydraulic conductivity correlated positively with water table height but negatively with altitude and cumulative short‐wave radiation since the last synoptic period of either negative air temperatures or turbulent energy flux dominance. The large range of K observed suggests complex interactions between meteorological influences and differences arising from variability in ice structure and crystallography. Our data demonstrate a greater complexity of near‐surface ice hydrology than hitherto appreciated and support the notion that the weathering crust can regulate the supraglacial discharge response to melt production. The conductivities reported here, coupled with typical supraglacial channel spacing, suggest that meltwater can be retained within the weathering crust for at least several days. Not only does this have implications for the accuracy of predictive meltwater run‐off models, but we also argue for biogeochemical processes and transfers that are strongly conditioned by water residence time and the efficacy of the cascade of sediments, impurities, microbes, and nutrients to downstream ecosystems. Because continued atmospheric warming will incur rising snowline elevations and glacier thinning, the supraglacial hydrological system may assume greater importance in many mountainous regions, and consequently, detailing weathering crust hydraulics represents a research priority because the flow path it represents remains poorly constrained.     

Public health impacts of climate change in Washington State: projected mortality risks due to heat events and air pollution

Illness and mortality related to heat and worsening air quality are core public health concerns associated with climate change projections. We examined the historical relationship between age- and cause-specific mortality rates from 1980 through 2006 and heat events at the 99th percentile of humidex values in the historic period from January 1, 1970 to December 31, 2006 in the greater Seattle area (King, Pierce and Snohomish counties), Spokane County, the Tri-Cities (Benton and Franklin counties) and Yakima County; the relative risks of mortality during heat events were applied to population and climate projections for Washington State to calculate number of deaths above the baseline (1980–2006) expected during projected heat events in 2025, 2045 and 2085. Three different warming scenarios were used in the analysis. Relative risks for the greater Seattle area showed a significant dose-response relationship between heat event duration and daily mortality rates for non-traumatic deaths for persons ages 45 and above, typically peaking at four days of exposure to humidex values above the 99th percentile. The largest number of projected excess deaths in all years and scenarios for the Seattle region was found for age 65 and above. Under the middle warming scenario, this age group is expected to have 96, 148 and 266 excess deaths from all non-traumatic causes in 2025, 2045 and 2085, respectively. We also examined projected excess deaths due to ground-level ozone concentrations at mid century (2045–2054) in King and Spokane counties. Current (1997–2006) ozone measurements and mid-twenty-first century ozone projections were coupled with dose-response data from the scientific literature to produce estimates overall and cardiopulmonary mortality. Daily maximum 8-h ozone concentrations are forecasted to be 16–28% higher in the mid twenty-first century compared to the recent decade of 1997–2006. By mid-century in King County the non-traumatic mortality rate related to ozone was projected to increase from baseline (0.026 per 100,000; 95% confidence interval 0.013–0.038) to 0.033 (95% CI 0.017–0.049). For the same health outcome in Spokane County, the baseline period rate of 0.058 (95% CI 0.030–0.085) was estimated increase to 0.068 (95% CI 0.035–0.100) by mid-century. The cardiopulmonary death rate per 100,000 due to ozone was estimated to increase from 0.011 (95% CI 0.005–0.017) to 0.015 (0.007–0.022) in King County, and from 0.027 (95% CI 0.013–0.042) to 0.032 (95% CI 0.015–0.049) in Spokane County. Public health interventions aimed at protecting Washington’s population from excessive heat and increased ozone concentrations will become increasingly important for preventing deaths, especially among older adults. Furthermore, heat and air quality related illnesses that do not result in death, but are serious nevertheless, may be reduced by the same measures.     

Adjusting to Policy Expectations in Climate Change Modeling

This paper surveys and interprets the attitudes of scientists to the use of flux adjustments in climate projections with coupled Atmosphere Ocean General Circulation Models. The survey is based largely on the responses of 19 climate modellers to several questions and a discussion document circulated in 1995. We interpret the responses in terms of the following factors: the implicit assumptions which scientists hold about how the environmental policy process deals with scientific uncertainty over human-related global warming; the different scientific styles that exist in climate research; and the influence of organisations, institutions, and policy upon research agendas. We find evidence that scientists' perceptions of the policy process do play a role in shaping their scientific practices. In particular, many of our respondents expressed a preference for keeping discussion of the issue of flux adjustments within the climate modeling community, apparently fearing that climate contrarians would exploit the issue in the public domain. While this may be true, we point to the risk that such an approach may backfire. We also identify assumptions and cultural commitments lying at a deeper level which play at least as important a role as perceptions of the policy process in shaping scientific practices. This leads us to identify two groups of scientists, pragmatists and purists, who have different implicit standards for model adequacy, and correspondingly are or are not willing to use flux adjustments.