Up until now, only a small number of studies have been dedicated to the binding processes of As(III) with organic matter (OM) via ionic Fe(III) bridges; none was interested in Fe (II). Complexation isotherms were carried out with As(III), Fe(II) or Fe(III) and Leonardite humic acid (HA). Although PHREEQC/Model VI, implemented with OM thiol groups, reproduced the experimental datasets with Fe(III), the poor fit between the experimental and modeled Fe(II) data suggested another binding mechanism for As(III) to OM. PHREEQC/Model VI was modified to take various possible As(III)-Fe(II)-OM ternary complex conformations into account. The complexation of As(III) as a mononuclear bidentate complex to a bidentate Fe (II)-HA complex was evidenced. However, the model needed to be improved since the distribution of the bidentate sites appeared to be unrealistic with regards to the published XAS data. In the presence of Fe (III), As(III) was bound to thiol groups which are more competitive with regards to the low density of formed Fe(III)-HA complexes. Based on the new data and previously published results, we propose a general scheme describing the various As(III)-Fe-MO complexes that are able to form in Fe and OM-rich waters. (C) 2016 Elsevier Inc. All rights reserved.