Systematic atomistic studies of surface reactivity for amorphous materials have not been possible in the past because of the complexity of these materials and the lack of the computer power necessary to draw representative statistics. With the emergence and popularization of machine learning (ML) approaches in materials science, systematic (and accurate) studies of the surface chemistry of disordered materials are now coming within reach. In this paper, we show how the reactivity of amorphous carbon (a-C) surfaces can be systematically quantified and understood by a combination of ML interatomic potentials, ML clustering techniques, and density functional theory calculations. This methodology allows us to process large amounts of atomic data to classify carbon atomic motifs on the basis of their geometry and quantify their reactivity toward hydrogen- and oxygen-containing functionalities. For instance, we identify subdivisions of sp and sp2 motifs with markedly different reactivities. We therefore draw a comprehensive, both qualitative and quantitative, picture of the surface chemistry of a-C and its reactivity toward -H, -O, -OH, and -COOH. While this paper focuses on a-C surfaces, the presented methodology opens up a new systematic and general way to study the surface chemistry of amorphous and disordered materials.