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 234U/238U 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 234U/238U ratios in vadose-zone waters resulting from lengthened non-flushing time as prescribed by the model provide an interpretative basis for using 234U/238U 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.