Sorption regulates the dispersion of pesticides from cropped areas to surrounding water bodies as well as their persistence. Assessing the risk of water contamination and evaluating the efficiency of mitigation measures, requires fine-resolution sorption data and a good knowledge of its drivers. This study aimed to assess the potential of a new approach combining chemometric and soil metabolomics to estimate the adsorption and desorption coefficients of a range of pesticides. It also aims to identify and characterise key components of soil organic matter (SOM) driving the sorption of these pesticides. We constituted a dataset of 43 soils from Tunisia, France and Guadeloupe (West Indies), covering extensive ranges of texture, organic carbon and pH. We performed untargeted soil metabolomics by liquid chromatography coupled with high-resolution mass spectrometry (UPLC-HRMS). We measured the adsorption and desorption coefficients of three pesticides namely glyphosate, 2,4-D and difenoconazole for these soils. We developed Partial Least Square Regression (PLSR) models for the prediction of the sorption coefficients from the RT-m/z matrix and conducted further ANOVA analyses to identify, annotate and characterise the most significant constituents of SOM in the PLSR models. The curated metabolomics matrix yielded 1213 metabolic markers. The prediction performance of the PLSR models was generally high for the adsorption coefficients Kdads (0.3 < R2 < 0.8) and for the desorption coefficients Kfdes (0.6 < R2 < 0.8) but low for ndes (0.03 < R2 < 0.3). The most significant features in the predictive models were annotated with a confidence level of 2 or 3. The molecular descriptors of these putative compounds suggest that the pool of SOM compounds driving glyphosate sorption is reduced compared to 2,4-D and difenoconazole, and these compounds are generally more polar. This approach can provide estimates of the adsorption and desorption coefficients of pesticides, including polar pesticide, for contrasted pedoclimates.