Glacier shrinkage, a notable consequence of climate change, is expected to intensify, particularly in high-elevation areas. While plant diversity and soil microbial communities have been studied, research on soil organic matter (SOM) and soil protein function dynamics in glacier forefields is limited. This limited understanding, especially regarding the link between microbial protein functions and biogeochemical functions, hampers our knowledge of soil-ecosystem processes along chronosequences. This study aims to elucidate the mechanistic relationships among soil bacterial protein functions, SOM decomposition, and environmental factors such as plant density and soil pH to advance understanding of the processes driving ecosystem succession in glacier forefields over time. Proteomic analysis showed that as ecosystems matured, the dominant protein functions transition from primarily managing cellular and physiological processes (biological controllers) to orchestrating broader ecological processes (ecosystem regulators) and increasingly include proteins involved in the degradation and utilization of OM. This shift was driven by plant density and pH, leading to increased ecosystem complexity and stability. Our confirmatory path analysis findings indicate that plant density is the main driver of soil process evolution, with plant colonization directly affecting pH, which in turn influenced nutrient metabolizing protein abundance, and SOM decomposition rate. Nutrient availability was primarily influenced by plant density, nutrient metabolizing proteins, and SOM decomposition, with SOM decomposition increasing with site age. These results underscore the critical role of plant colonization and pH in guiding soil ecosystem trajectories, revealing complex mechanisms and emphasizing the need for ongoing research to understand long-term ecosystem resilience and carbon sequestration.