Soils provide multiple ecosystem services (ES) including food production, biomass-energy and contribution to climate change regulation through carbon storage. However, implementing changes in land use and management can create trade-offs between services1. Therefore, new land planning approaches that consider ES relationships are needed to optimise overall ES delivery. Such an approach should be based on an understanding of the variability of soil bio-physico-chemical properties, both spatially within a given region and vertically across a soil profile, since these properties affect ES delivery2. Our study estimates the spatial repartition of ES delivery for various combinations of soil types and land occupation choices within a 300 km2 region dominated by agricultural and forested lands (Meuse / Haute Marne, France). Based on information from a 1/50 000 soil map and the description of 84 soil profiles, we separate the study area into homogeneous zones characterised by specific soil, land occupation and topographical properties. In each zone, soil functions and ES are rated on a scale of 0-3 using the DESTISOL decision support tool3. Maps of soil ecosystem services are then compared for different land cover scenarios. Firstly, we compare the evolution in soil ES delivery between 1850 and 2018. Secondly, we define realistic scenarios of land planning over the next 50 years, based on climatic, socioeconomic and public decision factors; finally, we propose scenarios that optimise each components of bundle of ES delivery at the regional scale. 1. Turkelboom, F. et al. When we cannot have it all: Ecosystem services trade-offs in the context of spatial planning. Ecosyst. Serv. 29, 566–578 (2018). 2. Martín, M. Á. et al. Soil structure and function in a changing world: Characterization and scaling. Geoderma 287, 1–3 (2017). 3. RPF 2008. 3. Blanchart, A. et al. Towards an operational methodology to optimize ecosystem services provided by urban soils. Landsc. Urban Plan. 176, 1–9 (2018).