Numerical modeling of marine radiocarbon reservoir ages during the last deglaciation

A critical problem in radiocarbon dating is the spatial and temporal variability of marine reservoir ages (MRAs). We assessed the MRA evolution during the last deglaciation by numerical modeling, applying a self-consistent iteration scheme in which an existing radiocarbon chronology (derived by Hughen et al., Quat. Sci. Rev., 25, pp. 3216-3227, 2006) was readjusted by transient, 3-D simulations of marine and atmospheric Delta14C. To estimate the uncertainties regarding the ocean ventilation during the last deglaciation, we considered various ocean overturning scenarios which are based on different climatic background states (PD: modern climate, GS: LGM climate conditions). Minimum and maximum MRAs are included in file 'MRAminmax_21-14kaBP.nc'. Three further files include MRAs according to equilibrium simulations of the preindustrial ocean (file 'C14age_preindustrial.nc'; this is an update of our results published in 2005) and of the glacial ocean (files 'C14age_spinupLGM_GS.nc' and 'C14age_spinupLGM_PD.nc').

(Table 1) Apparent ages from 40Ar/39Ar heating steps for DSDP Hole 72-516F basalts

Submarine basalts are difficult to date accurately by the potassium-argon method. Dalrymple and Moore (1968) and Dymond (1970), for example, showed that, when the conventional K-Ar method is used, pillow lavas may contain excess 40Ar. Use of the 40Ar/39Ar step-heating method has not overcome the problem, as had been hoped, and has produced some conflicting results. Ozima and Saito (1973) concluded that the excess 40Ar is retained only in high temperature sites, but Seidemann (1978) found that it could be released at all temperatures. Furthermore, addition of potassium, from seawater, to the rock after it has solidified can result in low ages (Seidemann, 1977), the opposite effect to that of excess 40Ar. Thus, apparent ages may be either greater or less than the age of extrusion. Because of this discouraging record, the present study was approached pragmatically, to investigate whether self-consistent results can be obtained by the 40Ar/39Ar step-heating method.

Physical properties of ODP Sites 123-765 and 123-766

During Ocean Drilling Program Leg 123, two sites were drilled in the deep Indian Ocean. Physical properties were measured in soft Quaternary and Lower Cretaceous sediments to relatively fresh, glass-bearing pillow lavas and massive basalts. Porosities ranged from 89% near the seafloor to 1.6% for the dense basalts. This self-consistent set of measurements permitted some descriptive models of physical properties to be more rigorously tested than before. Predictive relationships between porosity and compressional-wave velocity have generally been based upon the Wyllie time average equation. However, this equation does not adequately describe the actual relationship between these two parameters, and many have attempted to improve it. In most cases, models were derived by testing them against a set of data representing a relatively narrow range of porosity values. Similarly, the use of the Wyllie equation has often been justified by a pseudolinear fit to the data over a narrow range of porosity values. The limitations of the Wyllie relationship have been re-emphasized here. A semi-empirical acoustic impedance equation is developed that provides a more accurate porosity-velocity transform, using realistic material parameters, than has hitherto been possible. A closer correlation can be achieved with this semi-empirical relationship than with more theoretically based equations. In addition, a satisfactory empirical equation can be used to describe the relationship between thermal conductivity and porosity. If enough is known about core sample lithologies to provide estimates of the matrix and pore water parameters, then these predictive equations enable one to describe completely the behavior of a saturated rock core in terms of compressional-wave velocity, thermal conductivity, porosity, and bulk density.