Using machine learning to explore oxyanion adsorption ability of goethite with different specific surface area

In this study, we developed prediction models for the adsorption of divalent and trivalent oxyanions on goethite based on machine learning algorithms. After verifying the reliability of the models, the importance of goethite specific surface area (SSA) and the average oxyanion adsorption capacities of goethite with different SSAs were calculated by shapley additive explanations (SHAP) importance analysis and partial dependence (PD) analysis. Despite there were differences in the feature importance of divalent and trivalent oxyanions, the contribution of goethite's SSA to the adsorption amount ranked the fourth based on SHAP importance, indicating SSA played the important role in oxyanion adsorption. Meanwhile, the PD values of SSA and the optimized complexation constants from surface complexation modeling (SCM) both indicated a non-monotonic relationship between the goethite with different SSA and its oxyanions binding capacity. When the total site concentration and crystal face composition were used as the machine learning model input features, the SHAP importance values of crystal faces and the PD decomposition results indicated that the (001) face showed the crucial influence on oxyanions adsorption amount. These findings demonstrated the important role of crystal face composition in goethite's adsorption ability, and provided a theoretical explanation for the variations of oxyanions adsorption amount on different SSA goethite.