Modeling non-steady state radioisotope transport in the vadose zone - A case study using uranium isotopes at Peña Blanca, Mexico

Current models using U- and Th-series disequilibria to study radioisotope transport in groundwater systems mostly consider a steady-state situation. These models have limited applicability to the vadose zone (UZ) where the concentration and migratory behavior of radioisotopes in fluid are often transitory. We present here, as a first attempt of its kind, a model simulating the non-steady state, intermittent fluid transport in vadose layers. It provides quantitative constraints on in-situ migration of dissolved and colloidal radioisotopes in terms of retardation factor and rock-water interaction (or water transit) time. For uranium, the simulation predicts that intermittent flushing in the UZ leads to a linear relationship between reciprocal U concentration and ²³⁴U/²³⁸U ratio in percolating waters, with the intercept and slope bearing information on the rates of dissolution and α-recoil of U isotopes, respectively. The general validity of the model appears to be borne out by the measurement of uranium isotopes in UZ waters collected at various times over a period during 1995-2006 from a site in the Peña Blanca mining district, Mexico, where the Nopal I uranium deposit is located. Enhanced ²³⁴U/²³⁸U ratios in vadose-zone waters resulting from lengthened non-flushing time as prescribed by the model provide an interpretative basis for using ²³⁴U/²³⁸U in cave calcites to reconstruct the regional changes in hydrology and climate. We also provide a theoretical account of the model's potential applications using radium isotopes.