Archaeological age constraints from extrusion ages of obsidian: Examples from the Middle Awash,Ethiopia

Extrusion ages of archaeological obsidian, especially as determined by the ⁴⁰Ar/³⁹Ar method, can provide reliable maximum ages for tool manufacture. In at least one case in the Middle Awash of Ethiopia, freshly extruded obsidian was used for tool making, resulting in useful maximum ages for site occupation. Hydration resulting in mobility of K and/or Ar in glass, and recoil artifacts produced by neutron irradiation, fatally affect most glass shards from volcanic ashes. The much lower surface area to volume ratio of most archaeological obsidian, however, indicates that the affected areas can be manually removed prior to analysis and the recoil and hydration problems can be easily overcome. A more important issue in dating obsidian is that of possible mass-dependent kinetic isotope fractionation during or subsequent to quenching of volcanic glasses. This is evidenced in some cases by sub-atmospheric initial ⁴⁰Ar/³⁶Ar ratios, and more generally in sub-atmospheric ³⁸Ar/³⁶Ar. Resulting bias can be avoided through the use of isochron ages, which do not entail the assumption of an initial value of ⁴⁰Ar/³⁶Ar as is required for plateau ages. Since step heating of glasses often yields limited variability in ⁴⁰Ar:³⁹Ar:³⁶Ar (and therefore little spread on isochrons), another approach is to use an average value for initial ⁴⁰Ar/³⁶Ar, with concomitantly larger uncertainty than is associated with atmospheric ⁴⁰Ar/³⁶Ar, when calculating a plateau age. The ³⁸Ar/³⁶Ar of an un-irradiated subset of our samples validates the inference of kinetic fractionation, and potentially provides a basis for determining initial ⁴⁰Ar/³⁶Ar in samples that fail to yield isochrons, but only in samples lacking magmatic excess ⁴⁰Ar. These approaches allow us to reliably apply the ⁴⁰Ar/³⁹Ar method to volcanic glasses, which has resulted in maximum ages for archaeological sites that are not amenable to traditional geochronological methods. ⁴⁰Ar/³⁹Ar geochronology can also provide information on the geological provenance of the raw material used for tool making, especially when combined with geochemical data.     

Statistical Methods for Jointly Estimating the Decay Constant of 40K and the Age of a Dating Standard

To simultaneously evaluate the decay constant of ⁴⁰K () and the age of a standard (t _std) using isotopic data from geologic materials, we applied a series of statistical methods. The problem of estimating the most probable intercept of many nonlinear curves in and t _std space is formulated by an errors-in-variables nonlinear regression model. Then a maximum likelihood method is applied to the model for a point estimate, which is equivalent to the nonlinear least square method when measurement error distributions are Gaussian. Uncertainties and confidence regions of the estimates can be approximated using three methods: the asymptotic normal approximation, the parametric bootstrap method and Bonferroni confidence regions. Five pairs of published data for samples with ages from 2 ka to 4.5 Ga were used to estimate and the age of Fish Canyon sanidine (t _FCs). The statistical procedure yields most probable estimates of (5.4755 ± 0.0170 × 10^–10 (1)/year) and t _FCs (28.269 ± 0.0661 (1) Ma) which are in between previously published values. These results indicate the power of our approach to provide improved constraints on these parameters, although the preliminary nature of some of the input data require further review before the values can be adopted.     

Decay constants in geochronology

Geologic time is fundamental to the Earth Sciences, and progress in many disciplines depends critically on our ability to measure time with increasing accuracy and precision. Isotopic geochronology makes use of the decay of radioactive nuclides as a help to quantify the histories of rock, minerals, and other materials. Both accuracy and precision of radioisotopic ages are, at present, limited by those of radioactive decay constants. Modem mass spectrometers can measure isotope ratios with a precision of 10-4 or better. On the other hand, the uncertainties associated with direct half-life determinations are, in most cases, still at the percent level. The present short note briefly summarizes progress and problems that have been encountered during the Working Group‘s activity.     

Re–Os geochronology of a Mesoproterozoic sedimentary succession,Taoudeni basin,Mauritania: Implications for basin-wide correlations and Re–Os organic-rich sediments systematics

The exceptionally well-preserved sedimentary rocks of the Taoudeni basin, NW Africa represent one of the world's most widespread (> 1 M km²) Proterozoic successions. Hitherto, the sedimentary rocks were considered to be Mid Tonian based on Rb–Sr illite and glauconite geochronology of the Atar Group. However, new Re–Os organic-rich sediment (ORS) geochronology from two drill cores indicates that the Proterozoic Atar Group is ～ 200 Ma older (1107 ± 12 Ma, 1109 ± 22 Ma and 1105 ± 37 Ma). The Re–Os geochronology suggests that the Rb–Sr geochronology records the age of diagenetic events possibly associated with the Pan African collision.The new Re–Os geochronology data provide absolute age constraints for recent carbon isotope chemostratigraphy which suggests that the Atar Group is Mesoproterozoic and not Neoproterozoic.The new Re–Os ORS geochronology supports previous studies that suggest that rapid hydrocarbon generation (flash pyrolysis) from contact metamorphism of a dolerite sill does not significantly disturb the Re–Os ORS systematics. Modelled contact conditions suggest that the Re–Os ORS systematics remain undisturbed at ～ 650 °C at the sill/shale contact and ≥ 280 °C 20 m from the sill/shale contact.Moreover, the Re–Os geochronology indicates that the West African craton has a depositional history that predates 1100 Ma and that ORS can be correlated on a basin-wide scale. In addition, the Re–Os depositional ages for the ORS of the Taoudeni basin are comparable to those of ORS from the São Francisco craton, suggesting that these cratons are correlatable. This postulate is further supported by identical Os_i values for the Atar Group and the Vazante Group of the São Francisco craton.     

Silicic magmas from the continental Cameroon Volcanic Line (Oku, Bambouto and Ngaoundere): 40Ar-39Ar dates, petrology, Sr-Nd-O isotopes and their petrogenetic significance

The intraplate Cameroon Volcanic Line (CVL) straddles the African-South Atlantic continent-ocean boundary and is composed mainly of alkaline basic volcanic rocks. Voluminous silicic volcanics characterize the continental sector of the CVL. We present here new geochemical, isotopic (Sr-Nd-O) and ⁴⁰Ar/³⁹Ar geochronological data on the main silicic volcanic centres of the Western (Mt. Oku, Sabga and Mt. Bambouto) and Eastern (Ngaoundere plateau) Cameroon Highlands. The silicic volcanism of Mt. Oku, Sabga and Mt. Bambouto occurred between 25 and 15 Ma and is represented by voluminous quartz-normative trachytes and minor rhyolitic ignimbrites. At Mt. Bambouto central volcano about 700 m of silicic volcanics erupted in less than 2.7 million years. These silicic volcanics are associated with slightly to moderately alkaline basalts and minor basanites. In general, onset of the silicic volcanism migrated from NE (Oku: 25 Ma) to SW (Sabga: 23 Ma; Bambouto: 18 Ma; and Mt. Manengouba: 12 Ma). The silicic volcanism of the Ngaoundere plateau (eastern branch of the CVL) is instead dominated by nepheline-normative trachytes which are associated with strongly alkaline basalts and basanitic rocks. These Ne-trachytes are younger (11-9 Ma) than the Q-trachytes of the Western Highlands. The least differentiated silicic volcanics are isotopically similar (⁸⁷Sr/⁸⁶Sr < 0.70380; ¹⁴³Nd/¹⁴⁴Nd > 0.51278) to the associated alkaline basalts suggesting differentiation processes without appreciable interaction with crustal materials. Such interactions may, however, have played some role in the genesis of the most evolved silicic volcanics which have ⁸⁷Sr/⁸⁶Sr as high as 0.705–0.714. Fractional crystallization is the preferred mechanism for genesis of the silicic melts of both Western and Eastern Highlands, as shown by modeling major and trace element variations. The genesis of the least evolved Q-trachytes from the Western Highlands, starting from slightly to moderately alkaline basalts, is compatible with fractionation of dominantly plagioclase, clinopyroxene and magnetite. Crystal fractionation may have occurred at low pressure and at QFM buffer f _O2conditions. Parental magmas of the Ngaoundere Ne-trachytes are likely instead to have been strongly alkaline basalts which evolved through crystal fractionation at higher P (6-2 kbar) and f _O2 (QFM + 2). The migration (25 to 12 Ma) of the silicic volcanism from NE to SW in the continental sector of the CVL is reminiscent of that (31-5 Ma) of the onset of the basic volcanism in the oceanic sector (Principe to Pagalu islands) of the CVL. These ages, and that (11-9 Ma) of the silicic volcanism of the Ngaoundere plateau, indicate that the Cameroon Volcanic Line as a whole may not be easily interpreted as the surface expression of hot-spot magmatism. Received: 24 February 1998 / Accepted: 22 September 1998     

Assimilation of preexisting Pleistocene intrusions at Long Valley by periodic magma recharge accelerates rhyolite generation: rethinking the remelting model

Rhyolite flows and tuffs from the Long Valley area of California, which were erupted over a two-million-year time period, exhibit systematic trends in Nd, Hf, and Pb isotopes, trace element composition, erupted volume, and inferred magma residence time that provide evidence for a new model for the production of large volumes of silica-rich magma. Key constraints come from geochronology of zircon crystal populations combined with a refined eruption chronology from Ar–Ar geochronology; together these data give better estimates of magma residence time that can be evaluated in the context of changing magma compositions. Here, we report Hf, Nd, and Sr isotopes, major and trace element compositions, ⁴⁰Ar/³⁹Ar ages, and U–Pb zircon ages that combined with existing data suggest that the chronology and geochemistry of Long Valley rhyolites can be explained by a dynamic interaction of crustal and mantle-derived magma. The large volume Bishop Tuff represents the culmination of a period of increased mantle-derived magma input to the Long Valley volcanic system; the effect of this input continued into earliest postcaldera time. As the postcaldera evolution of the system continued, new and less primitive crustal-derived magmas dominated the system. A mixture of varying amounts of more mafic mantle-derived and felsic crustal-derived magmas with recently crystallized granitic plutonic materials offers the best explanation for the observed chronology, secular shifts in Hf and Nd isotopes, and the apparently low zircon crystallization and saturation temperatures as compared to Fe–Ti oxide eruption temperatures. This scenario in which transient crustal magma bodies remained molten for varying time periods, fed eruptions before solidification, and were then remelted by fresh recharge provides a realistic conceptual framework that can explain the isotopic and geochemical evidence. General relationships between crustal residence times and magma sources are that: (1) precaldera rhyolites had long crustal magma residence times and high crustal affinity, (2) the caldera-related Bishop Tuff and early postcaldera rhyolites have lower crustal affinity and short magma residence times, and (3) later postcaldera rhyolites again have stronger crustal signatures and longer magma residence times.     

A lattice Boltzmann model for noble gas diffusion in solids: The importance of domain shape and diffusive anisotropy and implications for thermochronometry

Thermochronometry based on radiogenic noble gases is critically dependent upon accurate knowledge of the kinetics of diffusion. With few exceptions, complex natural crystals are represented by ideal geometries such as infinite sheets, infinite cylinders, or spheres, and diffusivity is assumed to be isotropic. However, the physical boundaries of crystals generally do not conform to ideal geometries and diffusion within some crystals is known to be anisotropic. Our failure to incorporate such complexities into diffusive models leads to inaccuracies in both thermal histories and diffusion parameters calculated from fractional release data. To address these shortcomings we developed a code based on the lattice Boltzmann (LB) method to model diffusion from complex 3D geometries having isotropic, temperature-independent anisotropic, and temperature-dependent anisotropic diffusivity. In this paper we outline the theoretical basis for the LB code and highlight several advantages of this model relative to more traditional finite difference approaches. The LB code, along with existing analytical solutions for diffusion from simple geometries, is used to investigate the affect of intrinsic crystallographic features (e.g., crystal topology and diffusion anisotropy) on calculated diffusion parameters and a novel method for approximating thermal histories from crystals with complex topologies and diffusive anisotropy is presented.