Soil organic carbon stocks in European croplands and grasslands: How much have we lost in the past decade?

The EU Soil Strategy 2030 aims to increase soil organic carbon (SOC) in agricultural land to enhance soil health and support biodiversity as well as to offset greenhouse gas emissions through soil carbon sequestration. Therefore, the quantification of current SOC stocks and the spatial identification of the main drivers of SOC changes is paramount in the preparation of agricultural policies aimed at enhancing the resilience of agricultural systems in the EU. In this context, changes of SOC stocks (Δ SOCs) for the EU + UK between 2009 and 2018 were estimated by fitting a quantile generalized additive model (qGAM) on data obtained from the revisited points of the Land Use/Land Cover Area Frame Survey (LUCAS) performed in 2009, 2015 and 2018. The analysis of the partial effects derived from the fitted qGAM model shows that land use and land use change observed in the 2009, 2015 and 2018 LUCAS campaigns (i.e. continuous grassland [GGG] or cropland [CCC], conversion grassland to cropland (GGC or GCC) and vice versa [CGG or CCG]) was one of the main drivers of SOC changes. The CCC was the factor that contributed to the lowest negative change on Δ SOC with an estimated partial effect of -0.04 ± 0.01 g C kg-1 year-1 , while the GGG the highest positive change with an estimated partial effect of 0.49 ± 0.02 g C kg-1 year-1 . This confirms the C sequestration potential of converting cropland to grassland. However, it is important to consider that local soil and environmental conditions may either diminish or enhance the grassland's positive effect on soil C storage. In the EU + UK, the estimated current (2018) topsoil (0-20 cm) SOC stock in agricultural land below 1000 m a.s.l was 9.3 Gt, with a Δ SOC of -0.75% in the period 2009-2018. The highest estimated SOC losses were concentrated in central-northern countries, while marginal losses were observed in the southeast.

Including land management in a European carbon model with lateral transfer to the oceans.

The use of cover crops (CCs) is a promising cropland management practice with multiple benefits, notably in reducing soil erosion and increasing soil organic carbon (SOC) storage. However, the current ability to represent these factors in land surface models remains limited to small scales or simplified and lumped approaches due to the lack of a sediment-carbon erosion displacement scheme. This precludes a thorough understanding of the consequences of introducing a CC into agricultural systems. In this work, this problem was addressed in two steps with the spatially distributed CE-DYNAM model. First, the historical effect of soil erosion, transport, and deposition on the soil carbon budget at a continental scale in Europe was characterized since the early industrial era, using reconstructed climate and land use forcings. Then, the impact of two distinct policy-oriented scenarios for the introduction of CCs were evaluated, covering the European cropping systems where surface erosion rates or nitrate susceptibility are critical. The evaluation focused on the increase in SOC storage and the export of particulate organic carbon (POC) to the oceans, compiling a continental-scale carbon budget. The results indicated that Europe exported 1.95 TgC/year of POC to the oceans in the last decade, and that CCs can contribute to reducing this amount while increasing SOC storage. Compared to the simulation without CCs, the additional rate of SOC storage induced by CCs peaked after 10 years of their adoption, followed by a decrease, and the cumulative POC export reduction stabilized after around 13 years. The findings indicate that the impacts of CCs on SOC and reduced POC export are persistent regardless of their spatial allocation adopted in the scenarios. Together, the results highlight the importance of taking the temporal aspect of CC adoption into account and indicate that CCs alone are not sufficient to meet the targets of the 4‰ initiative. Despite some known model limitations, which include the lack of feedback of erosion on the net primary productivity and the representation of carbon fluxes with an emulator, the current work constitutes the first approach to successfully couple a distributed routing scheme of eroded carbon to a land carbon model emulator at a reasonably high resolution and continental scale. SHORT ABSTRACT: A spatially distributed model coupling erosion, transport, and deposition to the carbon cycle was developed. Then, it was used to simulate the impact of cover crops on both erosion and carbon, to show that cover crops can simultaneously increase organic carbon storage and reduce particulate organic carbon export to the oceans. The results seemed persistent regardless of the spatial distribution of cover crops.

Land use and climate change impacts on global soil erosion by water (2015-2070).

Soil erosion is a major global soil degradation threat to land, freshwater, and oceans. Wind and water are the major drivers, with water erosion over land being the focus of this work; excluding gullying and river bank erosion. Improving knowledge of the probable future rates of soil erosion, accelerated by human activity, is important both for policy makers engaged in land use decision-making and for earth-system modelers seeking to reduce uncertainty on global predictions. Here we predict future rates of erosion by modeling change in potential global soil erosion by water using three alternative (2.6, 4.5, and 8.5) Shared Socioeconomic Pathway and Representative Concentration Pathway (SSP-RCP) scenarios. Global predictions rely on a high spatial resolution Revised Universal Soil Loss Equation (RUSLE)-based semiempirical modeling approach (GloSEM). The baseline model (2015) predicts global potential soil erosion rates of [Formula: see text] Pg yr-1, with current conservation agriculture (CA) practices estimated to reduce this by ∼5%. Our future scenarios suggest that socioeconomic developments impacting land use will either decrease (SSP1-RCP2.6-10%) or increase (SSP2-RCP4.5 +2%, SSP5-RCP8.5 +10%) water erosion by 2070. Climate projections, for all global dynamics scenarios, indicate a trend, moving toward a more vigorous hydrological cycle, which could increase global water erosion (+30 to +66%). Accepting some degrees of uncertainty, our findings provide insights into how possible future socioeconomic development will affect soil erosion by water using a globally consistent approach. This preliminary evidence seeks to inform efforts such as those of the United Nations to assess global soil erosion and inform decision makers developing national strategies for soil conservation.

Can N2 O emissions offset the benefits from soil organic carbon storage?

To respect the Paris agreement targeting a limitation of global warming below 2°C by 2100, and possibly below 1.5°C, drastic reductions of greenhouse gas emissions are mandatory but not sufficient. Large-scale deployment of other climate mitigation strategies is also necessary. Among these, increasing soil organic carbon (SOC) stocks is an important lever because carbon in soils can be stored for long periods and land management options to achieve this already exist and have been widely tested. However, agricultural soils are also an important source of nitrous oxide (N2 O), a powerful greenhouse gas, and increasing SOC may influence N2 O emissions, likely causing an increase in many cases, thus tending to offset the climate change benefit from increased SOC storage. Here we review the main agricultural management options for increasing SOC stocks. We evaluate the amount of SOC that can be stored as well as resulting changes in N2 O emissions to better estimate the climate benefits of these management options. Based on quantitative data obtained from published meta-analyses and from our current level of understanding, we conclude that the climate mitigation induced by increased SOC storage is generally overestimated if associated N2 O emissions are not considered but, with the exception of reduced tillage, is never fully offset. Some options (e.g. biochar or non-pyrogenic C amendment application) may even decrease N2 O emissions.

Integrated management for sustainable cropping systems: Looking beyond the greenhouse balance at the field scale.

Cover crops (CC) promote the accumulation of soil organic carbon (SOC), which provides multiple benefits to agro-ecosystems. However, additional nitrogen (N) inputs into the soil could offset the CO2 mitigation potential due to increasing N2 O emissions. Integrated management approaches use organic and synthetic fertilizers to maximize yields while minimizing impacts by crop sequencing adapted to local conditions. The goal of this work was to test whether integrated management, centered on CC adoption, has the potential to maximize SOC stocks without increasing the soil greenhouse gas (GHG) net flux and other agro-environmental impacts such as nitrate leaching. To this purpose, we ran the DayCent bio-geochemistry model on 8,554 soil sampling locations across the European Union. We found that soil N2 O emissions could be limited with simple crop sequencing rules, such as switching from leguminous to grass CC when the GHG flux was positive (source). Additional reductions of synthetic fertilizers applications are possible through better accounting for N available in green manures and from mineralization of soil reservoirs while maintaining cash crop yields. Therefore, our results suggest that a CC integrated management approach can maximize the agro-environmental performance of cropping systems while reducing environmental trade-offs.

Mapping LUCAS topsoil chemical properties at European scale using Gaussian process regression.

This paper presents the second part of the mapping of topsoil properties based on the Land Use and Cover Area frame Survey (LUCAS). The first part described the physical properties (Ballabio et al., 2016) while this second part includes the following chemical properties: pH, Cation Exchange Capacity (CEC), calcium carbonates (CaCO3), C:N ratio, nitrogen (N), phosphorus (P) and potassium (K). The LUCAS survey collected harmonised data on changes in land cover and the state of land use for the European Union (EU). Among the 270,000 land use and cover observations selected for field visit, approximately 20,000 soil samples were collected in 24 EU Member States in 2009 together with more than 2000 samples from Bulgaria and Romania in 2012. The chemical properties maps for the European Union were produced using Gaussian process regression (GPR) models. GPR was selected for its capacity to assess model uncertainty and the possibility of adding prior knowledge in the form of covariance functions to the model. The derived maps will establish baselines that will help monitor soil quality and provide guidance to agro-environmental research and policy developments in the European Union.

Lateral carbon transfer from erosion in noncroplands matters.

This study combines two unprecedentedly high resolution (250 × 250 m) maps of soil erosion (inter-rill and rill processes) and soil organic carbon to calculate a global estimate of erosion-induced organic carbon (C) displacement. The results indicate a gross C displacement by soil erosion of 2.5-0.3+0.5 Pg C/year. The greatest share of displaced C (64%) comes from seminatural lands and forests. This suggests that lateral C transfer from erosion in noncroplands may play a more important role than previously assumed.

Quantifying the erosion effect on current carbon budget of European agricultural soils at high spatial resolution.

The idea of offsetting anthropogenic CO2 emissions by increasing global soil organic carbon (SOC), as recently proposed by French authorities ahead of COP21 in the 'four per mil' initiative, is notable. However, a high uncertainty still exits on land C balance components. In particular, the role of erosion in the global C cycle is not totally disentangled, leading to disagreement whether this process induces lands to be a source or sink of CO2. To investigate this issue, we coupled soil erosion into a biogeochemistry model, running at 1 km(2) resolution across the agricultural soils of the European Union (EU). Based on data-driven assumptions, the simulation took into account also soil deposition within grid cells and the potential C export to riverine systems, in a way to be conservative in a mass balance. We estimated that 143 of 187 Mha have C erosion rates <0.05 Mg C ha(-1) yr(-1), although some hot-spot areas showed eroded SOC >0.45 Mg C ha(-1) yr(-1). In comparison with a baseline without erosion, the model suggested an erosion-induced sink of atmospheric C consistent with previous empirical-based studies. Integrating all C fluxes for the EU agricultural soils, we estimated a net C loss or gain of -2.28 and +0.79 Tg yr(-1) of CO2 eq, respectively, depending on the value for the short-term enhancement of soil C mineralization due to soil disruption and displacement/transport with erosion. We concluded that erosion fluxes were in the same order of current carbon gains from improved management. Even if erosion could potentially induce a sink for atmospheric CO2, strong agricultural policies are needed to prevent or reduce soil erosion, in order to maintain soil health and productivity.

A new baseline of organic carbon stock in European agricultural soils using a modelling approach.

Proposed European policy in the agricultural sector will place higher emphasis on soil organic carbon (SOC), both as an indicator of soil quality and as a means to offset CO2 emissions through soil carbon (C) sequestration. Despite detailed national SOC data sets in several European Union (EU) Member States, a consistent C stock estimation at EU scale remains problematic. Data are often not directly comparable, different methods have been used to obtain values (e.g. sampling, laboratory analysis) and access may be restricted. Therefore, any evolution of EU policies on C accounting and sequestration may be constrained by a lack of an accurate SOC estimation and the availability of tools to carry out scenario analysis, especially for agricultural soils. In this context, a comprehensive model platform was established at a pan-European scale (EU + Serbia, Bosnia and Herzegovina, Croatia, Montenegro, Albania, Former Yugoslav Republic of Macedonia and Norway) using the agro-ecosystem SOC model CENTURY. Almost 164 000 combinations of soil-climate-land use were computed, including the main arable crops, orchards and pasture. The model was implemented with the main management practices (e.g. irrigation, mineral and organic fertilization, tillage) derived from official statistics. The model results were tested against inventories from the European Environment and Observation Network (EIONET) and approximately 20 000 soil samples from the 2009 LUCAS survey, a monitoring project aiming at producing the first coherent, comprehensive and harmonized top-soil data set of the EU based on harmonized sampling and analytical methods. The CENTURY model estimation of the current 0-30 cm SOC stock of agricultural soils was 17.63 Gt; the model uncertainty estimation was below 36% in half of the NUTS2 regions considered. The model predicted an overall increase of this pool according to different climate-emission scenarios up to 2100, with C loss in the south and east of the area (involving 30% of the whole simulated agricultural land) compensated by a gain in central and northern regions. Generally, higher soil respiration was offset by higher C input as a consequence of increased CO2 atmospheric concentration and favourable crop growing conditions, especially in northern Europe. Considering the importance of SOC in future EU policies, this platform of simulation appears to be a very promising tool to orient future policymaking decisions.

Potential carbon sequestration of European arable soils estimated by modelling a comprehensive set of management practices.

Bottom-up estimates from long-term field experiments and modelling are the most commonly used approaches to estimate the carbon (C) sequestration potential of the agricultural sector. However, when data are required at European level, important margins of uncertainty still exist due to the representativeness of local data at large scale or different assumptions and information utilized for running models. In this context, a pan-European (EU + Serbia, Bosnia and Herzegovina, Montenegro, Albania, Former Yugoslav Republic of Macedonia and Norway) simulation platform with high spatial resolution and harmonized data sets was developed to provide consistent scenarios in support of possible carbon sequestration policies. Using the CENTURY agroecosystem model, six alternative management practices (AMP) scenarios were assessed as alternatives to the business as usual situation (BAU). These consisted of the conversion of arable land to grassland (and vice versa), straw incorporation, reduced tillage, straw incorporation combined with reduced tillage, ley cropping system and cover crops. The conversion into grassland showed the highest soil organic carbon (SOC) sequestration rates, ranging between 0.4 and 0.8 t C ha(-1)  yr(-1) , while the opposite extreme scenario (100% of grassland conversion into arable) gave cumulated losses of up to 2 Gt of C by 2100. Among the other practices, ley cropping systems and cover crops gave better performances than straw incorporation and reduced tillage. The allocation of 12 to 28% of the European arable land to different AMP combinations resulted in a potential SOC sequestration of 101-336 Mt CO2 eq. by 2020 and 549-2141 Mt CO2 eq. by 2100. Modelled carbon sequestration rates compared with values from an ad hoc meta-analysis confirmed the robustness of these estimates.

European croplands under climate change: Carbon input changes required to increase projected soil organic carbon stocks.

Increasing soil organic carbon (SOC) stocks in agricultural systems is a pivotal strategy for promoting soil health and mitigating climate change. Global initiatives have set ambitious targets, aspiring to achieve an annual SOC stock increase of 4 ‰. In the European Union, the recently approved Nature Restoration Law aims to increase SOC stock trends in the top 30 cm of cropland mineral soils. However, current monitoring and reporting practices in some countries rely on simplistic SOC models with default parameters, which may not provide reliable predictions. In this paper, we study the feasibility of a 4 ‰ target in European croplands (i.e., an aspirational target proposed by The international "4 per 1000" Initiative), through estimations of required C input changes. To ensure robust predictions, we propose a novel calibration approach that links model parameters to pedo-climatic variables via statistical relationships from 16 long-term experiments. The effectiveness of the method is evaluated for three SOC models across 4281 sites from the European LUCAS soil survey. Our findings demonstrate that the statistical calibration of the multi-model ensemble improves the accuracy of 2015 and 2018 SOC stock predictions, compared to default parameterization. This improvement was however mainly due to the substantial enhancement of one of the models. According to the weighted multi-model mean, median C input changes to reach a 4 ‰ target for Northern, Central, and Southern Europe stand at 1.85, 1.20, and 0.13 Mg C ha-1 yr-1 under RCP 2.6, and 2.21, 1.26, and -0.10 Mg C ha-1 yr-1 under RCP 6.0, respectively. To achieve the aspirational 4 ‰ target, estimated C input change requirements exceed the predicted changes in net primary productivity under RCP 2.6 and RCP 6.0. This emphasizes the importance of strategic land-use and land-management interventions to enhance SOC stocks.

Quantifying the impact of an abrupt reduction in mineral nitrogen fertilization on crop yield in the European Union.

Contemporary crop production in Europe relies on nitrogen (N) fertilization. Fertilizer prices soared in 2021-2022, and remained at historical high levels in 2023. These high prices invoked an immediate concern on the possible consequences for Europe's food production. In this study, we use a biogeochemical model framework to estimate the impact of reducing mineral N fertilization on crop yields in the European Union (EU). First, crop yields simulated with the biogeochemical DayCent model are evaluated against subnational yield data averaged for 2015-2018 reported by Eurostat and National Statistical Institutes in the EU for soft wheat, barley, grain maize and rapeseed. Then, we simulate three different scenarios where mineral N fertilization across the EU is abruptly reduced by respectively 5, 15 and 25 %, and compare yields to the projected baseline for contemporary conditions (2019-2022). The model evaluation gives r2 values ranging from 0.28 (rapeseed) to 0.61 (soft wheat) and root mean square errors (RMSE) ranging from 0.6 (rapeseed) to 1.95 t ha-1 (maize). The model shows a reduction in yield per crop at the EU level up to 2.1, 6.4 and 11.2 % with the 5, 15 and 25 % reduction scenario, respectively. Different crops show different percentage reduction in yield following a reduction in mineral N fertilization, showing a legacy effect over the years and depending on the availability of organic fertilizer. The strongest relative yield reduction occurs for soft wheat for all three scenarios. Even with 25 % drop in mineral N fertilization, maize yield in the Netherlands, Belgium and Denmark is not significantly reduced, because of the high N surplus and large share of organic fertilization in these countries. This process-based modelling study provides spatially explicit, high resolution information on the response of crop yields to N fertilizer input reductions, helping policy-makers in decision-making on food security and environmentally-friendly food systems.

A unifying modelling of multiple land degradation pathways in Europe.

Land degradation is a complex socio-environmental threat, which generally occurs as multiple concurrent pathways that remain largely unexplored in Europe. Here we present an unprecedented analysis of land multi-degradation in 40 continental countries, using twelve dataset-based processes that were modelled as land degradation convergence and combination pathways in Europe's agricultural (and arable) environments. Using a Land Multi-degradation Index, we find that up to 27%, 35% and 22% of continental agricultural (~2 million km2) and arable (~1.1 million km2) lands are currently threatened by one, two, and three drivers of degradation, while 10-11% of pan-European agricultural/arable landscapes are cumulatively affected by four and at least five concurrent processes. We also explore the complex pattern of spatially interacting processes, emphasizing the major combinations of land degradation pathways across continental and national boundaries. Our results will enable policymakers to develop knowledge-based strategies for land degradation mitigation and other critical European sustainable development goals.

From regional to parcel scale: A high-resolution map of cover crops across Europe combining satellite data with statistical surveys.

The reformed Common Agricultural Policy of 2023-2027 aims to promote a more sustainable and fair agricultural system in the European Union. Among the proposed measures, the incentivized adoption of cover crops to cover the soil during winter provides numerous benefits such as improved soil structure and reduced nutrient leaching and erosion. Despite this recognized importance, the availability of spatial data on cover crops is scarce. The increasing availability of field parcel declarations in the European Union has not yet filled this data gap due to its insufficient information content, limited public availability and a lack of standardization at continental scale. At present, the best information available is regionally aggregated survey data, which although indicative, hinders the development of spatially accurate studies. In this work, we propose a statistical model relating Sentinel-1 data to the existence of cover crops at the 100-m spatial resolution over the entirety of the European Union and United Kingdom and estimate its parameters using the spatially aggregated survey data. To validate the method in a spatially-explicit way, predictions were compared against farmers' registered declarations in France, where the adoption of cover crops is widespread. The results indicate a good agreement between predictions and parcel-level data. When interpreted as a binary classifier, the model yielded an Area Under the Curve (AUC) of 0.74 for the whole country. When the country was divided into five regions for the evaluation of regional biases, the AUC values were 0.77, 0.75, 0.74, 0.70, and 0.65 for the North, Center, West, East, and South regions respectively. Despite limitations such as the lack of data for validation outside France, and the non-standardized nomenclature for cover crops among Member States, this work constitutes the first effort to obtain a relevant cover crop map at a European scale for researchers and practitioners.

Patterns in soil microbial diversity across Europe.

Factors driving microbial community composition and diversity are well established but the relationship with microbial functioning is poorly understood, especially at large scales. We analysed microbial biodiversity metrics and distribution of potential functional groups along a gradient of increasing land-use perturbation, detecting over 79,000 bacterial and 25,000 fungal OTUs in 715 sites across 24 European countries. We found the lowest bacterial and fungal diversity in less-disturbed environments (woodlands) compared to grasslands and highly-disturbed environments (croplands). Highly-disturbed environments contain significantly more bacterial chemoheterotrophs, harbour a higher proportion of fungal plant pathogens and saprotrophs, and have less beneficial fungal plant symbionts compared to woodlands and extensively-managed grasslands. Spatial patterns of microbial communities and predicted functions are best explained when interactions among the major determinants (vegetation cover, climate, soil properties) are considered. We propose guidelines for environmental policy actions and argue that taxonomical and functional diversity should be considered simultaneously for monitoring purposes.

Improving the phosphorus budget of European agricultural soils.

Despite phosphorus (P) being crucial for plant nutrition and thus food security, excessive P fertilization harms soil and aquatic ecosystems. Accordingly, the European Green Deal and derived strategies aim to reduce P losses and fertilizer consumption in agricultural soils. The objective of this study is to calculate a soil P budget, allowing the quantification of the P surpluses/deficits in the European Union (EU) and the UK, considering the major inputs (inorganic fertilizers, manure, atmospheric deposition, and chemical weathering) and outputs (crop production, plant residues removal, losses by erosion) for the period 2011-2019. The Land Use/Cover Area frame Survey (LUCAS) topsoil data include measured values for almost 22,000 samples for both available and total P. With advanced machine learning models, we developed maps for both attributes at 500 m resolution. We estimated the available P for crops at a mean value of 83 kg ha-1 with a clear distinction between North and South. The ratio of available P to the total P is about 1:17. The inorganic fertilizers and manure contribute almost equally as P inputs (mean 16 ± 2 kg P ha-1 yr-1 at 90 % confidence level) to agricultural soils, with high regional variations depending on farming practices, livestock density, and cropping systems. The P outputs came mainly from the exportation by the harvest of crop products and residues (97.5 %) and, secondly, by erosion. Using a sediment distribution model, we quantified the P fluxes to river basins and sea outlets. In the EU and UK, we estimated an average surplus of 0.8 kg P ha-1 yr-1 with high variability between countries with some regional variations. The P annual budget at regional scale showed ample possibility to improve P management by both reducing inputs in regions with high surplus (and P soil available) and rebalancing fertilization in those at risk of soil fertility depletion.

Soil erosion is unlikely to drive a future carbon sink in Europe.

Understanding of the processes governing soil organic carbon turnover is confounded by the fact that C feedbacks driven by soil erosion have not yet been fully explored at large scale. However, in a changing climate, variation in rainfall erosivity (and hence soil erosion) may change the amount of C displacement, hence inducing feedbacks onto the land C cycle. Using a consistent biogeochemistry-erosion model framework to quantify the impact of future climate on the C cycle, we show that C input increases were offset by higher heterotrophic respiration under climate change. Taking into account all the additional feedbacks and C fluxes due to displacement by erosion, we estimated a net source of 0.92 to 10.1 Tg C year-1 from agricultural soils in the European Union to the atmosphere over the period 2016-2100. These ranges represented a weaker and stronger C source compared to a simulation without erosion (1.8 Tg C year-1), respectively, and were dependent on the erosion-driven C loss parameterization, which is still very uncertain. However, when setting a baseline with current erosion rates, the accelerated erosion scenario resulted in 35% more eroded C, but its feedback on the C cycle was marginal. Our results challenge the idea that higher erosion driven by climate will lead to a C sink in the near future.

Copper distribution in European topsoils: An assessment based on LUCAS soil survey.

Copper (Cu) distribution in soil is influenced by climatic, geological and pedological factors. Apart from geological sources and industrial pollution, other anthropogenic sources, related to the agricultural activity, may increase copper levels in soils, especially in permanent crops such as olive groves and vineyards. This study uses 21,682 soil samples from the LUCAS topsoil survey to investigate copper distribution in the soils of 25 European Union (EU) Member States. Generalized Linear Models (GLM) were used to investigate the factors driving copper distribution in EU soils. Regression analysis shows the importance of topsoil properties, land cover and climate in estimating Cu concentration. Meanwhile, a copper regression model confirms our hypothesis that different agricultural management practices have a relevant influence on Cu concentration. Besides the traditional use of copper as a fungicide for treatments in several permanent crops, the combined effect of soil properties such as high pH, soil organic carbon and clay, with humid and wet climatic conditions favours copper accumulation in soils of vineyards and tree crops. Compared to the overall average Cu concentration of 16.85 mg kg-1, vineyards have the highest mean soil Cu concentration (49.26 mg kg-1) of all land use categories, followed by olive groves and orchards. Gaussian Process Regression (GPR) combined with kriging were used to map copper concentration in topsoils and to evidence the presence of outliers. GPR proved to be performant in predicting Cu concentration, especially in combination with kriging, accounting for 66% of Cu deviance. The derived maps are novel as they include information about the importance of topsoil properties in the copper mapping process, thus improving its accuracy. Both models highlight the influence of land management practices in copper concentration and the strong correlation between topsoil copper and vineyards.