Two-dimensional numerical modeling of Sr transport in an unsaturated Chinese loess under artificial sprinkling

⁹⁰Sr is a fission byproduct of uranium and plutonium, and it presents a major health problem in the environment. A field test on the transport of various radionuclides including ⁹⁰Sr in an unsaturated Chinese loess was conducted under artificial rain conditions from July 1997 to August 2000. The vertical concentration distribution of ⁹⁰Sr displayed an unusual profile of double concentration peaks, which were separated by a thin (0.7 cm) source layer. In order to interpret the double-peak concentration profile, the transport of ³H and ⁹⁰Sr in the unsaturated Chinese loess under artificial sprinkling conditions was simulated using WATERM, a numerical code for simulating flow field, and NESOR, also a numerical code but for simulating nuclide migration. The models were able to adequately simulate the double-peak concentration profile. The observation suggested that the fine arenaceous quartz layer, though 0.7 cm thick, formed a capillary barrier together with the local loess, which prevented water from penetrating. A significant discrepancy was observed between the model-fitted distribution coefficient (K_d) of ⁹⁰Sr and that determined from independent laboratory experiments, which can be attributed to a number of factors such as the capillary barrier effect, solution-to-solid ratio and soil water content. Therefore, when the model is used for predictive purposes where K_d is used as an input parameter, K_d must be determined under well controlled conditions by taking into account these factors as well as the heterogeneity in the field.