A 1951–80 global land precipitation climatology for the evaluation of general circulation models

Previous evaluations of model precipitation fields have suffered from two weaknesses; they have used only mean observed climatologies which have prevented an explicit evaluation of interannual variability, and they have generally failed to quantify the significance of differences between model and observed fields. To rectify these weaknesses, a global precipitation climatology is required which is designed with model evaluation in mind. This paper describes such a climatology representative of the period 1951–80. The climatology is based on historical gauge-precipitation measurements from over 2500 land-based station time series representing over 28% of the Earth's surface. It is necessarily biased towards terrestrial areas. The climatology (CRU5180) is derived from month-by-month gridbox precipitation estimates at 5° resolution. Although other global precipitation climatologies exist, this is the first one to have used a consistent reference period for each station, and to include the details of interannual variability. Fields of mean seasonal and annual precipitation and mean temporal variability are presented, and the variability of global-mean precipitation over 1951–80 assessed. The resulting mean monthly global precipitation fields are compared briefly with two other observed climatologies used for model evaluation, those prepared by Jaeger and Legates and Willmott. The global and hemispheric means, mean seasonal cycles, and spatial patterns of the three cimatologies are compared. Although based on a smaller set of stations than Legates and Willmott, the CRU5180 precipitation estimates agree closely with their uncorrected climatology.     

Petrogenesis and tectonic setting of bimodal volcanism in the Sakoli Mobile Belt, Central Indian shield

The Sakoli Mobile Belt comprises bimodal volcanic rocks that include metabasalt, rhyolite, tuffs, and epiclastic rocks with metapelites, quartzite, arkose, conglomerate, and banded iron formation (BIF). Mafic volcanic rocks are tholeiitic to quartz‐tholeiitic with normative quartz and hypersthene. SiO₂ shows a large compositional gap between the basic and acidic volcanics, depicting their bimodal nature. Both the volcanics have distinct geochemical trends but display some similarity in terms of enriched light rare earth element–large ion lithophile element characteristics with positive anomalies for U, Pb, and Th and distinct negative anomalies for Nb, P, and Ti. These characteristics are typical of continental rift volcanism. Both the volcanic rocks show strong negative Sr and Eu anomalies indicating fractionation of plagioclases and K‐feldspars, respectively. The high Fe/Mg ratios for the basic rocks indicate their evolved nature. Whole rock Sm–Nd isochrons for the acidic volcanic rocks indicate an age of crystallization for these volcanic rocks at about 1675 ± 180 Ma (initial ¹⁴³Nd/¹⁴⁴Nd = 0.51017 ± 0.00017, mean square weighted deviate [MSWD] = 1.6). The εNd_t (t = 2000 Ma) varies between ?0.19 and +2.22 for the basic volcanic rock and between ?2.85 and ?4.29 for the acidic volcanic rocks. Depleted mantle model ages vary from 2000 to 2275 Ma for the basic and from 2426 to 2777 Ma for the acidic volcanic rocks, respectively. These model ages indicate that protoliths for the acidic volcanic rocks probably had a much longer crustal residence time. Predominantly basaltic magma erupted during the deposition of the Dhabetekri Formation and part of it pooled at crustal or shallower subcrustal levels that probably triggered partial melting to generate the acidic magma. The influence of basic magma on the genesis of acidic magma is indicated by the higher Ni and Cr abundance at the observed silica levels of the acidic magma. A subsequent pulse of basic magma, which became crustally contaminated, erupted as minor component along with the dominantly acidic volcanics during the deposition of the Bhiwapur Formation.     

Creating a Population of 12‐Digit Headwater Basins within the Albemarle‐Pamlico Estuary System

Ecological research within the US Environmental Protection Agency's Office of Research and Development has recently changed its focus to quantifying and mapping ecosystem services provided to humans. Our local research group has been charged to develop a regional assessment of several ecosystem services in the Albemarle‐Pamlico Estuary System (APES). Time, data, and funding constraints precluded explicit modeling of the entire APES, so in Phase 1 of our research plan we chose to model ecosystem services in a random sample of headwater catchments. After observing numerous inconsistencies between the National Hydrography Dataset‐Plus (NHDPlus) stream network and the Virginia/North Carolina Watershed Boundary Dataset 12‐digit HUC coverage, we began by creating modified 12‐digit hydrologic units (HUCs) by aggregating smaller catchments delineated within the NHDPlus. In defining our population of interest (headwater 12‐digit HUCs with perennial, natural, wadeable pour points), we generally excluded HUCs that had multiple pour points, no pour points, or whose pour points were intermittent streams, artificial segments, ditches/canals, or lentic systems (lakes and reservoirs). After taking these actions, 318 HUCs remained and a stratified random sample (Omnerik ecoregions as strata) of 50 HUCs was chosen from this population.     

Methods for the visualization of digital elevation models for landform mapping

Digital elevation models (DEMs) are increasingly used for landform mapping, particularly with the growing availability of national and global datasets. In this paper we describe a variety of techniques that can visualize a DEM. We then compare five techniques to ascertain which performs the most complete and unbiased visualization. We assess the visualization techniques by comparing landforms mapped from them against a detailed morphological map (derived from mapping of multi‐azimuth relief‐shaded DEMs cross‐checked with stereo aerial photographs). Results show that no single visualization method provides complete and unbiased mapping. The relief‐shaded visualizations are particularly prone to azimuth biasing, although they can highlight subtle landforms. We recommend curvature visualization for initial mapping as this provides a non‐illuminated (and therefore unbiased) image. Initial mapping can then be supplemented with data from relief‐shaded visualizations. Copyright © 2005 John Wiley & Sons, Ltd.     

Quantifying the rate and depth dependence of bioturbation based on optically‐stimulated luminescence (OSL) dates and meteoric 10Be

Both the rate and the vertical distribution of soil disturbance modify soil properties such as porosity, particle size, chemical composition and age structure; all of which play an important role in a soil's biogeochemical functioning. Whereas rates of mixing have been previously quantified, the nature of bioturbation's depth dependence remains poorly constrained. Here we constrain, for the first time, the relationship between mixing rate and depth in a bioturbated soil in northeast Queensland, Australia using a novel method combining OSL (optically‐stimulated luminescence) ages and meteoric beryllium‐10 (¹⁰Be) inventories. We find that the best fit mixing rate decreases non‐linearly with increasing soil depth in this soil and the characteristic length scale of 0.28 m over which the mixing coefficient decays is comparable to reported rooting depth coefficients. In addition we show that estimates of surface mixing rates from OSL data are highly dependent on erosion rate and that erosion rate must be constrained if accurate mixing rates are to be quantified. We calculate surface diffusion‐like mixing coefficients of 1.8 × 10^?4 and 2.1 × 10^?4 m² yr^?1 for the studied soil for two different estimates of soil erosion. © 2014 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Geospatial Information Integration for Authoritative and Crowd Sourced Road Vector Data

This article describes results from a research project undertaken to explore the technical issues associated with integrating unstructured crowd sourced data with authoritative national mapping data. The ultimate objective is to develop methodologies to ensure the feature enrichment of authoritative data, using crowd sourced data. Users increasingly find that they wish to use data from both kinds of geographic data sources. Different techniques and methodologies can be developed to solve this problem. In our previous research, a position map matching algorithm was developed for integrating authoritative and crowd sourced road vector data, and showed promising results ( Anand et al. 2010 ). However, especially when integrating different forms of data at the feature level, these techniques are often time consuming and are more computationally intensive than other techniques available. To tackle these problems, this project aims at developing a methodology for automated conflict resolution, linking and merging of geographical information from disparate authoritative and crowd‐sourced data sources. This article describes research undertaken by the authors on the design, implementation, and evaluation of algorithms and procedures for producing a coherent ontology from disparate geospatial data sources. To integrate road vector data from disparate sources, the method presented in this article first converts input data sets to ontologies, and then merges these ontologies into a new ontology. This new ontology is then checked and modified to ensure that it is consistent. The developed methodology can deal with topological and geometry inconsistency and provide more flexibility for geospatial information merging.     

U–Pb and Th–Pb dating of apatite by LA-ICPMS

Apatite is a common U- and Th-bearing accessory mineral in igneous and metamorphic rocks, and a minor but widespread detrital component in clastic sedimentary rocks. U–Pb and Th–Pb dating of apatite has potential application in sedimentary provenance studies, as it likely represents first cycle detritus compared to the polycyclic behavior of zircon. However, low U, Th and radiogenic Pb concentrations, elevated common Pb and the lack of a U–Th–Pb apatite standard remain significant challenges in dating apatite by LA-ICPMS, and consequently in developing the chronometer as a provenance tool.This study has determined U–Pb and Th–Pb ages for seven well known apatite occurrences (Durango, Emerald Lake, Kovdor, Mineville, Mud Tank, Otter Lake and Slyudyanka) by LA-ICPMS. Analytical procedures involved rastering a 10 μm spot over a 40 × 40 μm square to a depth of 10 μm using a Geolas 193 nm ArF excimer laser coupled to a Thermo ElementXR single-collector ICPMS. These raster conditions minimized laser-induced inter-element fractionation, which was corrected for using the back-calculated intercept of the time-resolved signal. A Tl–U–Bi–Np tracer solution was aspirated with the sample into the plasma to correct for instrument mass bias. External standards (Ple?ovice and 91500 zircon, NIST SRM 610 and 612 silicate glasses and STDP5 phosphate glass) along with Kovdor apatite were analyzed to monitor U–Pb, Th–Pb, U–Th and Pb–Pb ratiosCommon Pb correction employed the ²⁰⁷Pb method, and also a ²⁰⁸Pb correction method for samples with low Th/U. The ²⁰⁷Pb and ²⁰⁸Pb corrections employed either the initial Pb isotopic composition or the Stacey and Kramers model and propagated conservative uncertainties in the initial Pb isotopic composition. Common Pb correction using the Stacey and Kramers (1975) model employed an initial Pb isotopic composition calculated from either the estimated U–Pb age of the sample or an iterative approach. The age difference between these two methods is typically less than 2%, suggesting that the iterative approach works well for samples where there are no constraints on the initial Pb composition, such as a detrital sample. No ²⁰⁴Pb correction was undertaken because of low ²⁰⁴Pb counts on single collector instruments and ²⁰⁴Pb interference by ²⁰⁴Hg in the argon gas supply.Age calculations employed between 11 and 33 analyses per sample and used a weighted average of the common Pb-corrected ages, a Tera–Wasserburg Concordia intercept age and a Tera–Wasserburg Concordia intercept age anchored through common Pb. The samples in general yield ages consistent (at the 2σ level) with independent estimates of the U–Pb apatite age, which demonstrates the suitability of the analytical protocol employed. Weighted mean age uncertainties are as low as 1–2% for U- and/or Th-rich Palaeozoic–Neoproterozoic samples; the uncertainty on the youngest sample, the Cenozoic (31.44 Ma) Durango apatite, ranges from 3.7–7.6% according to the common Pb correction method employed. The accurate and relatively precise common Pb-corrected ages demonstrate the U–Pb and Th–Pb apatite chronometers are suitable as sedimentary provenance tools. The Kovdor carbonatite apatite is recommended as a potential U–Pb and Th–Pb apatite standard as it yields precise and reproducible ²⁰⁷Pb-corrected, ²³²Th–²⁰⁸Pb, and common Pb-anchored Tera–Wasserburg Concordia intercept ages.     

Wide-angle seismic velocities in heterogeneous crust

Seismic velocities measured by wide-angle surveys are commonly used to constrain material composition in the deep crust. Therefore, it is important to understand how these velocities are affected by the presence of multiscale heterogeneities. The effects may be characterised by the scale of the heterogeneity relative to the dominant seismic wavelength (λ); what is clear is that heterogeneities of all scales and strengths bias wide-angle velocities to some degree. Waveform modelling was used to investigate the apparent wide-angle P -wave velocities of different heterogeneous lower crusts. A constant composition (50 per cent felsic and 50 per cent ultramafic) was formed into a variety of 1- and 2-D heterogeneous arrangements and the resulting wide-angle seismic velocity was estimated. Elastic, 1-D models produced the largest velocity shift relative to the true average velocity of the medium (which is the velocity of an isotropic mixture of the two components). Thick (width > λ) horizontal layers, as a result of Fermat's Principle, provided the largest increase in velocity; thin (width ≪λ) vertical layers produced the largest decrease in velocity. Acoustic 2-D algorithms were shown to be inadequate for modelling the kinematics of waves in bodies with multiscale heterogeneities. Elastic, 2-D modelling found velocity shifts (both positive and negative) that were of a smaller magnitude than those produced by 1-D models. The key to the magnitude of the velocity shift appears to be the connectivity of the fast (and/or slow) components. Thus, the models with the highest apparent levels of connectivity between the fast phases, the 1-D layers, produced the highest-magnitude velocity shifts. To understand the relationship between measured seismic velocities and petrology in the deep crust it is clear that high-resolution structural information (which describes such connectivity) must be included in any modelling.     

Muscovite dehydration melting: Reaction mechanisms,microstructures, and implications for anatexis

Dehydration melting of muscovite in metasedimentary sequences is the initially dominant mechanism of granitic melt generation in orogenic hinterlands. In dry (vapour-absent) crust, muscovite reacts with quartz to produce K-feldspar, sillimanite, and monzogranitic melt. When water vapour is present in excess, sillimanite and melt are the primary products of muscovite breakdown, and any K-feldspar produced is due to melt crystallization. Here we document the reaction mechanisms that control nucleation and growth of K-feldspar, sillimanite, and silicate melt in the metamorphic core of the Himalaya, and outline the microstructural criteria used to distinguish peritectic K-feldspar from K-feldspar grains formed during melt crystallization. We have characterized four stages of microstructural evolution in selected psammitic and pelitic samples from the Langtang and Everest regions: (a) K-feldspar nucleates epitaxially on plagioclase while intergrowths of fibrolitic sillimanite and the remaining hydrous melt components replace muscovite. (b) In quartzofeldspathic domains, K-feldspar replaces plagioclase by K⁺–Na⁺ cation exchange, while melt and intergrowths of sillimanite+quartz form in the aluminous domains. (c) At 7–8 vol.% melt generation, the system evolves from a closed to open system and all phases coarsen by up to two orders of magnitude, resulting in large K-feldspar porphyroblasts. (d) Preferential crystallization of residual melt on K-feldspar porphyroblasts and coarsened quartz forms an augen gneiss texture with a monzogranitic-tonalitic matrix that contains intergrowths of sillimanite+tourmaline+muscovite+apatite. Initial poikiloblasts of peritectic K-feldspar trap fine-grained inclusions of quartz and biotite by replacement growth of matrix plagioclase. During subsequent coarsening, peritectic K-feldspar grains overgrow and trap fabric-aligned biotite, resulting in a core to rim coarsening of inclusion size. These microstructural criteria enable a mass balance of peritectic K-feldspar and sillimanite to constrain the amount of free H₂O present during muscovite dehydration. The resulting modal proportion of K-feldspar in the Himalayan metamorphic core requires vapour-absent conditions during muscovite dehydration melting and leucogranite formation, indicating that the generation of large volumes of granitic melts in orogenic belts is not necessarily contingent on an external source of fluids.     

Lithogeochemical localisation of disseminated gold in the White River area, Yukon, Canada

Gold mineralisation in the White River area, 80 km south of the highly productive Klondike alluvial goldfield, is hosted in amphibolite facies gneisses in the same Permian metamorphic pile as the basement for the Klondike goldfield. Hydrothermal fluid which introduced the gold was controlled by fracture systems associated with middle Cretaceous to early Tertiary extensional faults. Gold deposition occurred where highly fractured and chemically reactive rocks allowed intense water–rock interaction and hydrothermal alteration, with only minor development of quartz veins. Felsic gneisses were sericitised with recrystallisation of hematite and minor arsenic mobility, and extensively pyritised zones contain gold and minor arsenic (ca 10 ppm). Graphitic quartzites (up to 5 wt.% carbon) caused chemical reduction of mineralising fluids, with associated recrystallisation of metamorphic minerals (graphite, pyrrhotite, pyrite, chalcopyrite) in host rocks and veins, and introduction of arsenic (up to 1 wt.%) to form arsenopyrite in veins and disseminated through host rock. Veins have little or no hydrothermal quartz, and up to 19 wt.% carbon as graphite. Late-stage oxidation of arsenopyrite in some graphitic veins has formed pharmacosiderite. Gold is closely associated with disseminated and vein sulphides in these two rock types, with grades of up to 3 ppm on the metre scale. Other rock types in the White River basement rocks, including biotite gneiss, hornblende gneiss, pyroxenite, and serpentinite, have not developed through-going fracture systems because of their individual mineralogical and rheological characteristics, and hence have been little hydrothermally altered themselves, have little hydrothermal gold, and have restricted flow of fluids through the rock mass. Some small post-metamorphic quartz veins (metre scale) have been intensely fractured and contain abundant gold on fractures (up to 40 ppm), but these are volumetrically minor. The style of gold mineralisation in the White River area is younger than, and distinctly different from, that of the Klondike area. Some of the mineralised zones in the White River area resemble, mineralogically and geochemically, nearby coeval igneous-hosted gold deposits, but this resemblance is superficial only. The White River mineralisation is an entirely new style of Yukon gold deposit, in which host rocks control the mineralogy and geochemistry of disseminated gold, without quartz veins.