Using geostatistical Bayesian Updating to integrate airborne radiometrics and soil geochemistry to improve mapping for mineral exploration

Mineral exploration programmes around the world use data from remote sensing, geophysics and direct sampling. On a regional scale, the combination of airborne geophysics and ground-based geochemical sampling can aid geological mapping and economic minerals exploration. The fact that airborne geophysical and traditional soil-sampling data are generated at different spatial resolutions means that they are not immediately comparable due to their different sampling density. Several geostatistical techniques, including indicator cokriging and collocated cokriging, can be used to integrate different types of data into a geostatistical model. With increasing numbers of variables the inference of the cross-covariance model required for cokriging can be demanding in terms of effort and computational time. In this paper a Gaussian-based Bayesian updating approach is applied to integrate airborne radiometric data and ground-sampled geochemical soil data to maximise information generated from the soil survey, to enable more accurate geological interpretation for the exploration and development of natural resources. The Bayesian updating technique decomposes the collocated estimate into a production of two models: prior and likelihood models. The prior model is built from primary information and the likelihood model is built from secondary information. The prior model is then updated with the likelihood model to build the final model. The approach allows multiple secondary variables to be simultaneously integrated into the mapping of the primary variable. The Bayesian updating approach is demonstrated using a case study from Northern Ireland where the history of mineral prospecting for precious and base metals dates from the 18th century. Vein-hosted, strata-bound and volcanogenic occurrences of mineralisation are found. The geostatistical technique was used to improve the resolution of soil geochemistry, collected one sample per 2 km2, by integrating more closely measured airborne geophysical data from the GSNI Tellus Survey, measured over a footprint of 65 x 200 m. The directly measured geochemistry data were considered as primary data in the Bayesian approach and the airborne radiometric data were used as secondary data. The approach produced more detailed updated maps and in particular maximized information on mapped estimates of zinc, copper and lead. Greater delineation of an elongated northwest/southeast trending zone in the updated maps strengthened the potential to investigate stratabound base metal deposits.

A reconnaissance study of the metamorphic petrology and volcanic geochemistry of a portion of the Island Lake greenstone belt, Manitoba /

The Island Lake greenstone belt is one of the major Archean supracrustal exposures in the northwestern part of the Superior Province of the Canadian Shield. This belt is subdivided into two units: 1) a lower sequence characterised by pillowed to massive, locally pyroclastic, basalt to andesite with a thin central zone of felsic derivatives, all of which are interbedded with and overlain by thick sequences of turbidite facies rock; 2) the upper unit which consists of thick stratified conglomerate overlain by thickly bedded arkose and feldspathic greywacke. Reconnaissance sampling traverses were completed across both the strike of the belt and along its margins with adjacent granitoids. Most of the belt is within the greenschist metamorphic f acies with amphibolite facies occurring in certain areas near t he margins. A post-tectonic, low pressure thermal event may be responsible for the development of a unit of cordierite schi s t which stretches southeastwards from the east end of Cochrane Bay. Volcanism is cyclical in nature changing from tholeiitic to calc-alkaline. There is a general progression in the character of the lavas from mafic t o felsic with stratigraphic height. Chemica l d a ta sugges t that h i gh level fractionation of a mantle- derived ' dry' magma i s t he s ource of the thole i iti c lavas. Contamination of this magma with 'we t' sia l and subsequent fractionation may be r esponsi b l e for the calcalkaline phases .Observations of stratigraphic relationships (in particular the contact between the supracrustals and the granitoids) coupled with the metamorphic and chemical studies, allow the construction of a preliminary model for the evolution of the Island Lake greenstone belt. The following sequential development is suggested: 1) a platform stage characterised by the subaqueous effusion of mafic to intermediate lavas of alternating tholeiitic and calc-alkaline affinities; 2) an edifice stage marked by the eruption of felsic calc-alkaline rocks; 3) an erosional stage characterised by the deposit~on of thick sequences of turbidite facies rocks; 4) the impingement of granitic masses into the margins of the greenstone belt, which was probably related to a downward warping of the supracrustal pilei 5) the erosion of sialic massifs surrounding and within the greenstone belt and of early supracrustal piles, to give the clastic upper unit.

Palynology and foraminiferal geochemistry of the Lower Pleistocene Olduvai Subchron (ca. 1.8 Ma) in DSDP Hole 603C, western North Atlantic

Marine palynology and benthic and planktonic foraminiferal geochemistry are combined to reveal long- and short-term (Milankovitch-scale) paleoceanographic changes across the upper half of the Olduvai Subchron (ca. 1.86--1.77 Ma, lower Pleistocene) in DSDP Hole 603C from the lower New Jersey continental rise. Planktonic foraminiferal Mg/Ca ratios reveal annual sea-surface temperatures between 14.5° and 25°C, whereas modern values vary between 16° and 20°e. Despite evidence of downslope transport in much of the studied interval, dinoflagellate cyst and acritarch assemblages appear to reflect fluctuating temperate to subtropical water masses. These assemblages comprise both neritic and oceanic species, and are marked by a transition upsection from warm conditions, dominated by Lingulodinium machaerophorum, Polysphaeridium zoharyi and Cymatiosphaera? invaginata, to cooler conditions dominated by Filisphaera filifera. Combining dinoflagellate cyst proxies with planktonic foraminiferal geochemistry allows downslope transport events to be recognized during glacial episodes, and events dominated by intensified bottom-water circulation during interglacial episodes. Sixtytwo in-situ dinoflagellate cyst and acritarch taxa were recorded including several not previously described.

Mesoproterozoic rapakivi granites of the Rondonia Tin Province, southwestern border of the Amazonian craton, Brazil - I. Reconnaissance U-Pb geochronology and regional implications

Rapakivi granites and associated mafic and ultramafic rocks in the Rondonia Tin Province, southwestern Amazonian craton, Brazil were emplaced during six discrete episodes of magmatism between ca 1600 and 970 Ma. The seven rapakivi granite suites emplaced at this time were the Serra da Providencia Intrusive Suite (U-Pb ages between 1606 and 1532 Ma); Santo Antonio Intrusive Suite(U-Pb age 1406 Ma), Teotonio Intrusive Suite (U-Pb age 1387 Ma); Alto Candeias Intrusive Suite (U-Pb ages between 1346 and 1338 Ma); Sao Lourenco-Caripunas Intrusive Suite (U-Pb ages between 1314 and 1309 Ma); Santa Clara Intrusive Suite (U-Pb ages between 1082 and 1074 Ma); and Younger Granites of Rondonia (U-Pb ages between 998 and 974 Ma). The Serra da Providencia Intrusive Suite intruded the Paleoproterozoic (1.80 to 1.70 Ga) Rio Negro-Juruena crust whereas the other suites were emplaced into the 1.50 to 1.30 Ga Rondonia-San Ignacio crust. Their intrusion was contemporaneous with orogenic activity in other parts of the southwestern Amazonian craton, except for the oldest, Serra da Providencia Intrusive Suite. Orogenic events coeval with emplacement of the Serra da Providencia Intrusive Suite are not clearly recognized in the region. The Santo Antonio, Teotonio, Alto Candeias and Sao Lourenco-Caripunas Intrusive Suites are interpreted to represent extensional anorogenic magmatism associated with the terminal stages of the Rondonian-San Ignacio orogeny. At least the Sao Lourenco-Caripunas rapakivi granites and coeval intra-continental rift sedimentary rocks may, in contrast, represent the products of extensional tectonics and rifting preceding the Sunsas/Aguapei orogeny (1.25 to 1.0 Ga). The two youngest rapakivi suites, the Santa Clara Intrusive Suite and Younger Granites of Rondonia, seemingly represent inboard magmatism in the Rondonian-San Ignacio Province during a younger episode of reworking in the Rio Negro-Juruena Province during the waning stages of the collisional 1.1 to 1.0 Ga Sunsas/Aguapei orogeny. The six intra-plate rapakivi granite episodes in the southwestern part of the Amazonian craton form three broad periods of anorogenic magmatism that have age-correlative events composed of similar rocks and geologic environments in eastern Laurentia and Baltica, although the exact timing of magmatism appears slightly different. Recognition of lithologic and chronological correlations between various cratons provide important constraints to models explaining the interplay between rapakivi granite magmatism and deep crustal evolution of an early Mesoproterozoic supercontinent. They are, furthermore, important to plate tectonic models for the assembly, dispersal and reassembly of Amazonia, Laurentia and Baltica in the Mesoproterozoic and Neoproterozoic. (C) 1999 Elsevier B.V. B.V. All rights reserved.

LEAD DISTRIBUTION AND ISOTOPE SIGNATURE IN BOTTOM SEDIMENTS OF GUAMA RIVER AND GUAJARA BAY (BELEM - PARA)

LEAD DISTRIBUTION AND ISOTOPE SIGNATURE IN BOTTOM SEDIMENTS OF GUAMA RIVER AND GUAJARA BAY (BELEM - PARA). Lead analyses in bottom sediments from the hydrographic system of Belem (Para) indicated low contents of this metal for the sediments from the Guama river, with no significant anthropogenic contribution. A concentration of 18.1 +/- 1.5 mg kg(-1) and Pb-206/Pb-207 isotopic signature of 1.196 +/- 0.002 are assigned for Pb from natural sources. On the other hand, the significant increase of Pb contents in the sediments from the Guajara bay, together with the decrease of Pb-206/Pb-207 ratios (1.172 < Pb-206/Pb-207 < 1.188) point to an anthropogenic lead contribution, originated by the industrial and urban activities of the city of Belem.

Provenance signature of a forearc basin modified by spreading ridge subduction: Detrital zircon geochronology and detrital modes from the Paleogene Arkose Ridge Formation, southern Alaska

Upper Paleocene–Eocene boulder conglomerate, cross-stratified sandstone, and laminated carbonaceous mudstone of the Arkose Ridge Formation exposed in the southern Talkeetna Mountains record fluvial-lacustrine deposition proximal to the volcanic arc in a forearc basin modified by Paleogene spreading ridge subduction beneath southern Alaska. U-Pb ages of detrital zircon grains and modal analyses were obtained from stratigraphic sections spanning the 2,000 m thick Arkose Ridge Formation in order to constrain the lithology, age, and location of sediment sources that provided detritus. Detrital modes from 24 conglomerate beds and 54 sandstone thin sections aredominated by plutonic and volcanic clasts and plagioclase feldspar with minor quartz, schist, hornblende, argillite, and metabasalt. Westernmost sandstone and conglomerate strata contain <5% volcanic clasts whereas easternmost sandstone and conglomerate strata contain 40 to >80% volcanic clasts. Temporally, eastern sandstones andconglomerates exhibit an upsection increase in volcanic detritus from <40 to >80% volcanic clasts. U-Pb ages from >1400 detrital zircons in 15 sandstone samples reveal three main populations: late Paleocene–Eocene (60-48 Ma; 16% of all grains), Late Cretaceous–early Paleocene (85–60 Ma; 62%) and Jurassic–Early Cretaceous (200–100 Ma; 12%). A plot of U/Th vs U-Pb ages shows that >97% of zircons are <200 Ma and>99% of zircons have <10 U/Th ratios, consistent with mainly igneous source terranes. Strata show increased enrichment in late Paleocene–Eocene detrital zircons from <2% in the west to >25% in the east. In eastern sections, this younger age population increases temporally from 0% in the lower 50 m of the section to >40% in samples collected >740 m above the base. Integration of the compositional and detrital geochronologic data suggests: (1) Detritus was eroded mainly from igneous sources exposed directly north of the Arkose Ridge Formation strata, mainly Jurassic–Paleocene plutons and Paleocene–Eocenevolcanic centers. Subordinate metamorphic detritus was eroded from western Mesozoic low-grade metamorphic sources. Subordinate sedimentary detritus was eroded from eastern Mesozoic sedimentary sources. (2) Eastern deposystems received higher proportions of juvenile volcanic detritus through time, consistent with construction of adjacent slab-window volcanic centers during Arkose Ridge Formation deposition. (3)Western deposystems transported detritus from Jurassic–Paleocene arc plutons that flank the northwestern basin margin. (4) Metasedimentary strata of the Chugach accretionaryprism, exposed 20-50 km south of the Arkose Ridge Formation, did not contribute abundant detritus. Conventional provenance models predict reduced input of volcanic detritus to forearc basins during exhumation of the volcanic edifice and increasing exposure ofsubvolcanic plutons (Dickinson, 1995; Ingersoll and Eastmond, 2007). In the forearc strata of these conventional models, sandstone modal analyses record progressive increases upsection in quartz and feldspar concomitant with decreases in lithic grains, mainly volcanic lithics. Additionally, as the arc massif denudes through time, theyoungest detrital U-Pb zircon age populations become significantly older than the age of forearc deposition as the arc migrates inboard or ceases magmatism. Westernmost strata of the Arkose Ridge Formation are consistent with this conventional model. However, easternmost strata of the Arkose Ridge Formation contain sandstone modes that record an upsection increase in lithic grains accompanied by a decrease in quartz and feldspar, and detrital zircon age populations that closely match the age of deposition. This deviation from the conventional model is due to the proximity of the easternmost strata to adjacent juvenile volcanic rocks emplaced by slab-window volcanic processes. Provenance data from the Arkose Ridge Formation show that forearc basins modified by spreading ridge subduction may record upsection increases in non-arc, syndepositional volcanic detritusdue to contemporaneous accumulation of thick volcanic sequences at slab-window volcanic centers. This change may occur locally at the same time that other regions of the forearc continue to receive increasing amounts of plutonic detritus as the remnant arc denudes, resulting in complex lateral variations in forearc basin petrofacies and chronofacies.

Bias associated with detrital zircon geochronology and thermochronometry

Detrital studies that utilize zircon U–Pb geochronology and fission-track (FT) thermochronometry are subject to a range of potential sources of bias that should be properly evaluated and minimized. Some of them are common to any single-grain mineral analysis (e.g., variable bedrock mineral fertility, hydraulic sorting during transport, selective grain loss during sample processing), whereas others are intrinsic to zircon, and are related to radiation damage and age discordance. In this article, we quantify the impact of intrinsic bias on detrital studies thanks to the analysis of modern detritus shed from the European Alps, and illustrate the general implications on geological interpretations. We show that detrital zircon U–Pb age distributions based on statistically robust datasets are highly reproducible and representative of the parent bedrock ages in the catchment. Arbitrary or selective removal of discordant grain ages can be minimized by using the Kolmogorov–Smirnov test to identify an appropriate cutoff level. Loss of metamict (α-damaged) zircon has a minor impact on data representativeness, and is mainly controlled by regional metamorphism rather than by mechanical abrasion during river transport. Zircon FT grain-age distributions were found to have poor reproducibility, although age spectra are consistent with bedrock data. However, unlike the U–Pb datasets, U-rich zircon grains (> 1000 ppm) are systematically missed, and undatable grains may exceed 50%. We identify two major sources of distribution bias specific to zircon FT datasets: (i) sediment sources dominated by U-rich zircon grains are markedly underrepresented in the detrital record, because such grains often have uncountable high densities of fission tracks (“U concentration bias”); (ii) sediment sources that shed zircon grains with high levels of α-damage are underrepresented, because these grains are lost during chemical etching for FT revelation (“etching bias”). In the case of multimethod dating on the same grains (e.g., FT and U–Pb double dating), bias affecting detrital zircon FT dating propagates to the entire dataset. These effects may not impact on exhumation-rate studies that utilize the youngest grain ages (i.e., lag-time approach). However, they represent a limiting factor for conventional provenance studies, and generally preclude application of zircon FT dating to sediment budget calculations.

Neotethys closure history of Anatolia: insights from 40Ar–39Ar geochronology and P–T estimation in high pressure metasedimentary rocks

The multiple high-pressure (HP), low-temperature (LT) metamorphic units of Western and Central Anatolia offer a great opportunity to investigate the subduction- and continental accretion-related evolution of the eastern limb of the long-lived Aegean subduction system. Recent reports of the HP–LT index mineral Fe-Mg-carpholite in three metasedimentary units of the Gondwana-derived Anatolide–Tauride continental block (namely the Afyon Zone, the Ören Unit and the southern Menderes Massif) suggest a more complicated scenario than the single-continental accretion model generally put forward in previous studies. This study presents the first isotopic dates (white mica 40Ar–39Ar geochronology), and where possible are combined with P–T estimates (chlorite thermometry, phengite barometry, multi-equilibrium thermobarometry), on carpholite-bearing rocks from these three HP–LT metasedimentary units. It is shown that, in the Afyon Zone, carpholite-bearing assemblages were retrogressed through greenschist-facies conditions at c. 67–62 Ma. Early retrograde stages in the Ören Unit are dated to 63–59 Ma. In the Kurudere–Nebiler Unit (HP Mesozoic cover of the southern Menderes Massif), HP retrograde stages are dated to c. 45 Ma, and post-collisional cooling to c. 26 Ma. These new results support that the Ören Unit represents the westernmost continuation of the Afyon Zone, whereas the Kurudere–Nebiler Unit correlates with the Cycladic Blueschist Unit of the Aegean Domain. In Western Anatolia, three successive HP–LT metamorphic belts thus formed: the northernmost Tavşanlı Zone (c. 88–82 Ma), the Ören–Afyon Zone (between 70 and 65 Ma), and the Kurudere–Nebiler Unit (c. 52–45 Ma). The southward younging trend of the HP–LT metamorphism from the upper and internal to the deeper and more external structural units, as in the Aegean Domain, points to the persistence of subduction in Western Anatolia between 93–90 and c. 35 Ma. After the accretion of the Menderes–Tauride terrane, in Eocene times, subduction stopped, leading to continental collision and associated Barrovian-type metamorphism. Because, by contrast, the Aegean subduction did remain active due to slab roll-back and trench migration, the eastern limb (below Southwestern Anatolia) of the Hellenic slab was dramatically curved and consequently teared. It therefore is suggested that the possibility for subduction to continue after the accretion of buoyant (e.g. continental) terranes probably depends much on palaeogeography.

Coupling of LMS with a fs-laser ablation ion source: elemental and isotope composition measurements

Mass spectrometric analysis of elemental and isotopic compositions of several NIST standards is performed by a miniature laser ablation/ionisation reflectron-type time-of-flight mass spectrometer (LMS) using a fs-laser ablation ion source (775 nm, 190 fs, 1 kHz). The results of the mass spectrometric studies indicate that in a defined range of laser irradiance (fluence) and for a certain number of accumulations of single laser shot spectra, the measurements of isotope abundances can be conducted with a measurement accuracy at the per mill level and at the per cent level for isotope concentrations higher and lower than 100 ppm, respectively. Also the elemental analysis can be performed with a good accuracy. The LMS instrument combined with a fs-laser ablation ion source exhibits similar detection efficiency for both metallic and non-metallic elements. Relative sensitivity coefficients were determined and found to be close to one, which is of considerable importance for the development of standard-less instruments. Negligible thermal effects, sample damage and excellent characteristics of the fs-laser beam are thought to be the main reason for substantial improvement of the instrumental performance compared to other laser ablation mass spectrometers.