Surface complexation modeling of arsenic mobilization from goethite: Interpretation of an in-situ experiment

Sorption has been recognized as a predominant process for the mobility and transport of arsenic (As) in groundwater. However, the model-based description of the chemical and physical mechanisms controlling As interaction with mineral surfaces under natural hydrochemical conditions remains a formidable challenge. In this study we assess and compare the ability of existing surface complexation models (SCMs), the diffuse double layer model (DDL) and the charge distribution multisite complexation model (CD-MUSIC), to simulate As desorption from goethite in groundwater. We consider the outcomes of an in-situ experiment recently performed in an arsenic-contaminated aquifer of Northern China where As-loaded goethite-coated sand was deployed in 7 monitoring wells. Determination and measurements of As-surface species were carried out over a time period of 80 days in all monitoring wells. Simultaneous calibration of the models with the measurements in the 7 wells was performed to obtain single sets of surface complexation parameters for the DDL and the CD-MUSIC models, respectively. Although the general dynamic pattern of the As release at the site is captured by both models, the approach based on the CD-MUSIC agrees best with the field experimental data without modifications of the surface complexation database compiled from previous studies on goethite. Conversely, calibration of the DDL affinity constants led to substantial improvement in the agreement between model simulations and the considered field dataset.