Purpose Evaluation of soil functionality in Life Cycle Assessment (LCA) has progressively gained importance, although only a small cluster of studies deliver detailed guidelines on how to calculate quantified indicators. In addition, there is a lack of bibliography assessing impacts on the pedosphere due to spatially differentiated land use changes (LUC). In this study, an automated geospatial simulation of LUC based on crop rotation probabilities in Luxembourg was implemented in the programming environment R. Furthermore, this method based on coupling LCA and geographic information system (GIS) was used to calculate changes in soil functionality by implementing both the soil organic carbon (SOC) method and the Land Use Indicator Value Calculation Tool (LANCA®). The developed R script was then applied to a case study dealing with maize production for bioenergy purposes in Luxembourg.