A new method integrating high-precision U-Pb geochronology with zircon trace element analysis (U-Pb TIMS-TEA)

Increased precision in isotope-dilution thermal ionization mass spectrometry (ID-TIMS) U-Pb geochronology has revealed age complexities in zircon populations that require new tools for understanding how the growth of zircon is related to geologic processes. U and Pb are routinely separated from other elements in dated minerals by ion exchange separation prior to TIMS isotope measurement. We develop a method in which trace elements in the exact same volume of zircon are redissolved and analyzed using solution nebulization inductively coupled plasma sector-field mass spectrometry with matrix-matched external liquid calibration. Using <0.5 ml solution, resulting concentrations are between <1 ppt for elements such as Ti, Nb and Ta and tens of ppb for Zr. By analyzing a series of standard solutions, zircons and procedural blanks, we show that accurate measurements are performed on Zr, Hf, Y, Sc, and the HREE while low-concentration elements can be measured accurately to <5 ppt. We performed combined U-Pb ID-TIMS geochronology with trace element analysis (here called U-Pb TIMS-TEA) on zircons from eight volcanic rocks comprising several volcanic systems and one metamorphic sample. Similar to previous in situ trace element analyses, zircon geochemistry is distinct between different samples and records petrogenetic processes such as fractional crystallization, assimilation and/or magma mixing. Unique from in situ analysis, U-Pb TIMS-TEA can trace geochemical evolution in accessory minerals with adequate age precision to resolve magmatic processes in rocks at least 200 million years old. This provides a means to identify auto-, ante- and xenocrystic zircon and lead to more robust age interpretations in ID-TIMS U-Pb geochronology. One suite of Cretaceous andesitic zircons shows correlations in geochemistry and absolute time that record evolution of a magmatic system over ∼250 ka prior to eruption. Future work will combine U-Pb TIMS-TEA with solution isotopic analysis of Nd, Sr and Hf and will be applied to a host of datable minerals such as monazite, sphene, apatite, rutile, xenotime, and baddeleyite. These combined tools will provide access to an improved understanding of a wide range of igneous and metamorphic processes as a function of time.     

Determining accurate temperature-time paths from U-Pb thermochronology: An example from the Kaapvaal craton, southern Africa

Thermochronology has revolutionized our understanding of the establishment and evolution of lithospheric thermal structure. However, many potential benefits provided by the application of diffusion theory to thermochronology have yet to be fully exploited. This study uses apatite (T_c = 450-550 °C) and titanite (T_c = 550-650 °C) U-Pb ID-TIMS thermochronology at the single- to sub-grain scale to separate the variable effects of volume diffusion of Pb from metamorphic (over)growth above and below the T_c of a mineral. Data are presented from two ca. 3227 Ma tonalite samples from north and south of the Barberton Greenstone Belt (BGB), southern Africa. Two distinct populations of apatite from a sample north of the BGB record fast cooling followed by metamorphic growth ∼10 Myr later. Both apatite and titanite dates from south of the BGB show a strong correlation with the grain size and record 100 Myr of post-emplacement cooling. Complex core-rim zoning observed in cathodoluminescence images of apatite is interpreted to reflect metamorphic overgrowth above the T_c. The age and topology of grain size versus date curves from titanite and apatite are used in combination with a finite-difference numerical model to show that slow, non-linear, cooling and not thermal resetting is responsible for the observed distribution. The thermal histories from either side of the BGB are very different and provide unique insight into the BGB’s tectonic evolution: a ∼70 Myr period of apparent stability after ca. 3.2 Ga terrane assembly was followed by fast exhumation south of the BGB that led to lower-crustal melting and intrusion of granitic batholiths ca. 3.14-3.10 Ga.     

Arsenic attenuation in tailings at a former Cu–W–As mine,SW Finland

Nearly half a century after mine closure, release of As from the Ylöjärvi Cu–W–As mine tailings in groundwater and surface water run-off was observed. Investigations by scanning electron microscopy (SEM), electron microprobe analysis (EMPA), synchrotron-based micro-X-ray diffraction (μ-XRD), micro-X-ray absorption near edge structure (μ-XANES) and micro-extended X-ray absorption fine structure (μ-EXAFS) spectroscopy, and a sequential extraction procedure were performed to assess As attenuation mechanisms in the vadose zone of this tailings deposit. Results of SEM, EMPA, and sequential extractions indicated that the precipitation of As bearing Fe(III) (oxy)hydroxides (up to 18.4 wt.% As₂O₅) and Fe(III) arsenates were important secondary controls on As mobility. The μ-XRD, μ-XANES and μ-EXAFS analyses suggested that these phases correspond to poorly crystalline and disordered As-bearing precipitates, including arsenical ferrihydrite, scorodite, kaňkite, and hydrous ferric arsenate (HFA). The pH within 200 cm of the tailings surface averaged 5.7, conditions which favor the precipitation of ferrihydrite. Poorly crystalline Fe(III) arsenates are potentially unstable over time, and their transformation to ferrihydrite, which contributes to As uptake, has potential to increase the As adsorption capacity of the tailings. Arsenic mobility in tailings pore water at the Ylöjärvi mine will depend on continued arsenopyrite oxidation, dissolution or transformation of secondary Fe(III) arsenates, and the As adsorption capacity of Fe(III) (oxy)hydroxides within this tailings deposit.     

Curve-Fitting Schemes for Fragmentation Data

Fragmentation measurements in the form of sieve passing and mass fraction data were used to test the capability of three different distributions to fit the observed data over a wide range in fragment size and mass. These distributions were based on Rosin-Rammler, lognormal and simple sigmoidal (S-shaped) functions, having 2 input parameters for the single-component versions and 5 input parameters for the two-component versions. Provided convergence was achieved in the non-linear curve-fitting technique, the two-component versions always provided superior fits to the observed data. However, these versions were very sensitive to variations in the values chosen for the input parameters. In this particular regard, the two-component sigmoidal function was the most robust. The present results also show that the two-component lognormal function provided the best fit to the fragmentation data in a general sense, and the two-component Rosin-Rammler function provided the worst fit. However, there was not a significant difference between any of the three methods.     

Triggered-lightning properties inferred from measured currents and very close electric fields

Triggered-lightning properties, including dart-leader charge density, return-stroke propagation speed, dart-leader electric potential, dart-leader propagation speed, and dart-leader current, inferred from return-stroke current and very close electric field measurements, are presented. Although most of the estimates are based on relatively crude models, they are all generally in good agreement with independent measurements and/or theoretical considerations found in the literature. The results are likely to be applicable to subsequent strokes in natural lightning.     

Sedimentary processes and facies of the waterlaid Anvil Spring Canyon alluvial fan, Death Valley, California

The Anvil Spring Canyon fan of the Panamint Range piedmont in central Death Valley was built entirely by water-flow processes, as revealed by an analysis of widespread 2- to 12-m-high stratigraphic cuts spanning the 9·7 km radial length of this 2·5–5·0° sloping fan. Two facies deposited from fan sheetfloods dominate the fan from apex to toe. The main one (60–95% of cuts) consists of sandy, granular, fine to medium pebble gravel that regularly and sharply alternates with cobbly coarse to very coarse pebble gravel in planar couplets 5–25 cm thick oriented parallel to the fan surface. The other facies (0–25% of cuts) comprises 10- to 60-cm-thick, wedge-planar and wedge-trough beds of pebbly sand and sandy pebble gravel in backsets sloping 3–28°. Both facies are interpreted as resulting from rare, sediment-charged flash floods from the catchment, and were deposited by supercritical standing waves of expanding sheetfloods on the fan. Standing waves were repeatedly initiated, enlarged, migrated, and then terminated either by gradually rejoining the flood or by more violent breakage and washout. The frequent autocyclic growth and destruction of standing waves during an individual sheetflood resulted in the deposition of multiple coarse and fine couplet and backset sequences 50–250 cm thick across the active depositional lobe of the fan. Erosional intensity during washout of the standing wave determined whether early-phase backset-bed deposits or washout-phase sheetflood couplet deposits were selectively preserved in a given cycle. Two minor facies are also found in the Anvil fan. Pebble–cobble gravel lags (0–20% of cuts) are present above erosional scours into the sheetflood couplet and backset deposits. They consist of coarse gravel concentrated through fine-fraction winnowing of the host sheetflood facies by sediment-deficient water flows. This reworking occurred during recessional flood stage or from non-catastrophic discharge during the long intervals between major flash floods. This facies is common at the surface, giving rise to a ‘braided-stream’ appearance. However, it is stratigraphically limited, present as thin, continuous to discontinuous beds or lenses that bound 50- to 250-cm-thick sheetflood sequences. The other minor facies of the Anvil fan consists of clast-supported and imbricated, thickly stratified, pebbly, cobbly, boulder gravel present in narrow, radially aligned ribbons nested within sheetflood deposits. This facies is interpreted as representing deposition in the incised channel of the fan, a subenvironment characterized by greater flow competence resulting from maintained depth from channel-wall confinement, and by more frequent water flows and winnowing events caused by its direct connection with the catchment feeder channel.     

Cause of dominance by sheetflood vs. debris-flow processes on two adjoining alluvial fans, Death Valley, California

Two large, adjoining alluvial fans of the Panamint Range piedmont, Death Valley, California, are composed of different facies assemblages deposited by contrasting sedimentary processes. The Anvil Spring fan was built solely by water-flow processes (incised-channel floods and sheetfloods), whereas the neighbouring Warm Spring fan has been constructed principally by debris flows. The boundary between these fans delineates a sharp provincial piedmont contact between sheetflood-dominated fans to the south and debris-flow-dominated fans to the north. Factors such as climate, catchment area, fan area, catchment relief, aspect, vegetation types and density, and neotectonic setting are essentially identical for these two fans. The key difference between them is that their catchments are underlain by dissimilar bedrock types, which weather to yield distinctive sediment suites. Weathering of the granite and andesite of the Anvil fan catchment produces significant volumes of medium to very coarse sand, granules, pebbles, cobbles and boulders, but minimal silt and clay. In contrast, the shale, quartzite and dolomite that dominate bedrock in the Warm Spring catchment weather to yield a wide suite of sedimentary particles spanning from clay to boulders. The abundance of mud, and the unsorted character of the yielded sediment, cause precipitation-induced slope failures in the Warm Spring catchment to transform readily into debris flows. This propensity is due to the low permeability of the colluvial sediment, which causes added water to become trapped quickly and pore pressure to rise rapidly, promoting transformations to debris flows. In contrast, the limited volume of sediment finer than medium sand yielded from the Anvil fan catchment causes the colluvium to have high permeability. This factor prevents the transformation of wet colluvium to a debris flow during hydrologically triggered slope failures, instead maintaining sediment transport as entrained bed load or suspended load in a water flow.     

Modelling study of three-dimensional circulation and particle movement over the Sable Gully of Nova Scotia

The Sable Gully is a broad deep underwater canyon located to the east of Sable Island on the edge of the Scotian Shelf. Being the home of many marine species including the endangered Northern Bottlenose Whale, the Gully was designated as a marine protected area (MPA) in 2004. Better understanding of physical environmental conditions over this MPA is needed for sustainable ecosystem management. In this study, a multi-nested ocean circulation model and a particle tracking model are used to examine the three-dimensional (3D) circulation and movement of particles carried passively by the flow over the Sable Gully. The 3D circulation model is driven by tides, wind, and surface heat/freshwater fluxes. The model performance is assessed by comparing the results with the previous numerical tidal results and current meter observations made in the Gully. The simulated tidal circulation over the Gully and adjacent waters is relatively strong on shallow banks and relatively weak on the continental slope. Below the depth of the Gully rim ( ～ 200 m), the tidal currents are constrained by the thalweg of the Gully and amplified toward the Gully head. The simulated subtidal circulation in the Gully has a complex spatial structure and significant seasonal variability. The simulated time-dependent 3D flow fields are then used in a particle tracking model to study the particle movements, downstream and upstream areas, and residence time of the Gully. Based on the movements of particles released at the depth of the Gully rim and tracked forward in time, the e-folding residence time is estimated to be about 7 and 13 days in February and August 2006, respectively. The Gully flanks are identified as high retention areas with the typical residence time of 10 and 20 days in February and August 2006, respectively. Tracking particles with and without tides reveals that tidal circulation reduces the value of residence time in the Gully, particularly along the Gully flanks.