Spatio – Temporal Dynamics of Soil Organic Carbon Stock in Greek Croplands: A Long-Term Assessment

1. Introduction

The intricate interplay of SOC within terrestrial ecosystems is a cornerstone in comprehending soil health, nutrient cycling, and the broader dynamics of the global carbon cycle. The complexity of SOC's role in influencing soil fertility, water retention, and overall ecosystem health underscores its critical significance in the pursuit of sustainable agriculture and effective climate change mitigation strategies [1,2,3,4,5]. In this work, the spatio-temporal dynamics of SOC stock in Greek croplands are examined, where diverse agroecosystems and varying environmental conditions beckon to explore the nuanced relationships between SOC and key soil properties [6,7].

Extending our gaze beyond local landscapes, SOC's significance transcends boundaries, resonating with global implications for climate change mitigation. SOC emerges as a crucial ally in the ongoing battle against climate change by sequestering atmospheric carbon dioxide and actively participating in greenhouse gas mitigation efforts on a global scale [2,8,9,10]. The increasing global emphasis on carbon sequestration reinforces the imperative for a nuanced understanding of SOC dynamics within specific regional contexts.

The canvas of Greek croplands unfolds as a mosaic of agroecological diversity, woven by the threads of climatic variability, topographical nuances, and intricate land management practices [11,12]. From the expansive plains of Thessaly to the terraced slopes of Crete, the multifaceted nature of Greek croplands prompts a comprehensive investigation. This exploration transcends mere scientific inquiry; it is an endeavor to decipher the influences shaping SOC dynamics within the unique context of Greek agriculture.

The relationship between anthropogenic activities and SOC stock unfolds as a tightrope walk in the socio-economic fabric of the Mediterranean, especially in Greece. While croplands contribute substantially to the socio-economic tapestry, human interventions and land management practices can wield a transformative influence on SOC stock [13,14]. This delicate balance between agricultural productivity and soil health becomes the focal point, necessitating a meticulous examination of the dynamic interplay within SOC.

The classification bestowed by Soil Order emerges as a key determinant, weaving a complex tapestry of SOC variability. Drawing from documented evidence, Histosols in wetland areas boast elevated SOC content, emphasizing their organic-rich nature [15,16,17,18]. In contrast, the Regosol order may present a different narrative, with lower SOC stock attributed to shallower profiles and limited organic inputs [19,20]. The influence of Soil Order becomes a narrative thread guiding us through the diverse SOC landscapes within Greek croplands.

The concentration of organic carbon in surface horizons is a well-documented phenomenon; however, the depths of the soil profile emerge as substantial carbon reservoirs with far-reaching implications for long-term sequestration potential [21,22,23]. This vertical perspective becomes integral for holistic assessments of carbon stock.

Coarse materials, such as gravel and rocks, are pivotal components introducing heterogeneity in the soil narrative, influencing water movement, drainage, and organic matter distribution, thus creating microenvironments that significantly impact microbial activity and decomposition rates [22,24,25]. Expanding our understanding of the intricate role played by these coarse materials adds layers of complexity to SOC distribution in Greek croplands.

Navigating the microscopic landscape of SOC dynamics, Soil Texture takes center stage as a defining factor in shaping carbon stock. The relative proportions of sand, silt, and clay particles influence SOC stock by orchestrating water retention, aeration, and microbial activity [26,27]. Clayey soils, characterized by higher Cation Exchange Capacity (CEC), emerge as crucial custodians of organic matter retention, contributing significantly to higher SOC stock compared to their sandy counterparts. This microscopic world of soil particles, often overlooked, unfolds as a vital driver in the SOC narrative within Greek croplands, emphasizing the need for a holistic understanding of soil texture dynamics for effective land management.

The interplay between Drainage and Slope introduces a narrative of erosion risks and SOC reallocation. Steeper slopes, more prone to erosion, become areas where topsoil and organic carbon-rich horizons are more vulnerable for removal. Understanding the relationship between slope, drainage, and SOC stock becomes pivotal, offering a compass to predict erosion risks and implement effective soil conservation measures [2,5,20]. The topographical intricacies of Greek croplands infuse heightened relevance into these factors.

Within the spectrum of acidity and alkalinity, Soil pH plays an important role in microbial activity, nutrient availability, and organic matter stability. The variation in microbial decomposition rates with pH levels adds layers of complexity to the SOC dynamics [28,29].

In the synthesis of existing knowledge and the unveiling of novel insights, this manuscript embarks on a scientific endeavor to enrich our understanding of the spatio-temporal dynamics of SOC stock in Greek croplands. The comprehensive consideration of an array of soil properties, as outlined above, becomes the lens through which we seek to enhance our comprehension of the intricate relationships governing SOC sequestration and turnover. This synthesis offers not just a compilation of data, but a narrative woven from diverse soil properties, providing valuable insights for sustainable land management practices and climate change mitigation strategies in the unique context of Greek croplands.

The objectives of this study can be summarized as follows: 1. Examine Soil Properties: Investigate a range of soil properties, including pH, soil texture, cation exchange capacity (CEC), drainage, slope, and the presence of coarse materials, to understand their influence on soil organic carbon (SOC) dynamics, 2. Understand SOC Loss Patterns: Analyze how SOC loss patterns vary across different soil orders, soil depths, and other categories such as pH and soil texture, 3. Contextualize Within Greek Agriculture: Focus on the unique context of Greek agriculture to provide specific insights into the factors influencing SOC dynamics in this region, 4. Balance Agricultural Productivity and Soil Health: Highlight the critical need to examine the dynamic interplay of soil properties to achieve a balance between maintaining agricultural productivity and ensuring soil health, and, 5. Provide Nuanced Perspectives: Offer detailed perspectives on how various soil properties interact and contribute to SOC dynamics, helping to inform better soil management practices in Greek agriculture.

2. Materials and Methods

2.1. Data Sources

This study focuses on the different landscapes of Greece, offering a rich tapestry of agroecological diversity shaped by climatic variations and intricate land management practices. The Greek croplands serve as the canvas for investigating SOC dynamics. To examine the spatio-temporal intricacies of SOC stock, we rely on comprehensive datasets from authoritative sources.

The SOC stock map [30] from the European Soil Data Centre (ESDAC) by the Joint Research Centre (JRC, European Commission) [31,32,33] serves as a foundational dataset, providing detailed insights into the distribution and variability of SOC across Greek croplands (Figure 1a). This map, derived from advanced modeling and monitoring techniques, offers a robust foundation for understanding the SOC landscape. Additionally, we integrate the SOC stock map provided by the Food and Agriculture Organization of the United Nations (FAO) [34,35], which was prepared by our team at the Institute of Soil and Water Resources of ELGO DIMITRA (as the National Focal Point of Global Soil Partnership - FAO) (Figure 1b). The FAO's SOC stock map contributes valuable perspectives, enhancing the comprehensiveness of our study. These datasets with a spatial resolution of 1 Km or 30 Arc-Second, meticulously compiled by international organizations, enable us to explore the nuanced relationships between SOC and key soil properties in the unique context of Greek agriculture. The Greek Soil Map with a 1:30.000 scale, is given in Figure 2 [36]. This study endeavors to contribute novel insights into the dynamics of SOC stock, fostering a deeper understanding of sustainable land management practices and facilitating effective climate change mitigation strategies in the Greek cropland terrain. The two SOC maps refer to the topsoil (0-30 depth).

2.2. Methodology

This research employs a comprehensive methodology to investigate the spatio-temporal dynamics of SOC stock in Greek croplands. The study design encompasses a multi-faceted approach to capture the diverse agroecological nuances and varying environmental conditions across the Greek landscape. The temporal scope of the study spans 11 years (2010-2021), enabling the analysis of SOC variations over time. The integration of authoritative datasets, including the SOC maps from the European Soil Data Centre (ESDAC) by the Joint Research Centre (JRC) and the Food and Agriculture Organization (FAO), forms the cornerstone of our methodology.

2.2.1. Data Collection

The primary data sources for this study are the SOC maps from ESDAC (JRC) and FAO. These maps provide a spatially explicit representation of SOC content, allowing for a detailed exploration of SOC variability within Greek croplands. The ESDAC (JRC) map, derived from advanced modeling techniques and remote sensing technologies, offers a high-resolution depiction of SOC distribution. Concurrently, the FAO SOC map, a globally recognized resource, contributes additional insights, ensuring the robustness of our analysis.

More specifically, the soil data for the model came from the European Soil Database at the European Soil Data Centre (ESDAC), focusing on the top-soil layer (0–30 cm) with properties like texture, bulk density, pH, drainage class, and rock content at a 1 km × 1 km grid resolution. CENTURY was chosen for the pan-European SOC assessment due to its integration of crop growth routines, successful testing in Europe, ability to simulate management practices, and reduced computational time. Model results were validated against inventories from EIONET and 20,000 soil samples from the 2009 LUCAS survey, aiming to create a harmonized top-soil data set for the EU.

The FAO SOC map methodology utilizes the SCORPAN model framework for digital soil mapping, which predicts SOC stocks based on soil-forming factors. Key variables for estimating SOC stocks include environmental parameters: climate data, thematic maps, digital terrain data, geomorphometry, and soil data. Modelling techniques, such as random forests, support vector machines, and regression-kriging are employed to estimate SOC stocks. The performance of the Greek National SOC map is evaluated using both internal and external accuracy metrics.

The issues raised by comparing datasets from different years with different methodology can be mitigating by applying analysis of uncertainties in each SOC stock map and evaluation of the quality and accuracy of both datasets. Moreover, checking for any biases, errors, or inconsistencies using statistical methods is required to quantify the reliability of the data. All the above issues were addressed in both 2010 and 2021 SOC stock maps [30,34]. Therefore, any inconsistence produced by comparison of these two datasets was effectively minimized because of their robustness and proper methodological execution.

2.2.2. Spatiο - Temporal Analysis

Spatial analysis forms a crucial component of our methodology, leveraging Geographic Information System (GIS) tools to analyze and interpret SOC patterns. We delineate the study area into distinct regions, considering factors such as agroecological zones, soil orders, and topographical characteristics. By overlaying the SOC maps with additional spatial datasets, such as land use and slope, we aim to discern relationships between SOC stock and key environmental variables.

To explore temporal trends in SOC stock, we conduct a temporal analysis using multiple datasets representing different time points. This temporal dimension allows us to identify patterns, anomalies, and potential drivers of SOC variations over time. We examine the influence of climatic factors, land management practices, and other temporal variables on SOC dynamics within the Greek cropland terrain.

2.2.3. Statistical Analysis

Normality and equality of variances were tested with Q-Q plots and Levene’s test. Kruskal Wallis test (non-parametric equivalent to One Way Analysis of Variance) was used to determine whether significant differences existed between groups of variable ΔSOC stock (SOC 2021 - SOC 2010) for each factor. If there were significant differences, Dunn’s post hoc test was performed (with Holm correction for multiple statistical tests). In each case, the Kruskal Wallis test assesses whether there are significant differences between the groups defined by the respective factor. The p-value is a measure of the evidence against a null hypothesis, with smaller p-values suggesting stronger evidence against the null hypothesis. A p-value less than 0.05 (95% confidence interval) indicates statistically significant differences between the groups for each factor. Boxplots for variable SOC stock 2010 and 2021 present the median, minimum, maximum, 1st and 3rd quantile and outliers. The median is reported instead of the mean in all cases. JASP (0.18.3.0) [37] was used for statistical analysis of the results.

2.2.4. Synthesis and Interpretation

This involves synthesizing the diverse datasets and drawing meaningful interpretations. The synthesis aims to examine the complex interplay of SOC dynamics within Greek croplands, considering the spatial and temporal dimensions, as well as the intricate relationships with key soil properties. The interpretation of results contributes to a comprehensive understanding of the factors influencing SOC sequestration and turnover in the context of sustainable land management practices and climate change mitigation strategies in Greek agriculture.

The methodology outlined above integrates advanced remote sensing, GIS, statistical analyses, and field data to provide a holistic investigation into the intricate dynamics of SOC stock in Greek croplands. This multi-dimensional approach ensures the robustness and reliability of our findings, facilitating valuable contributions to the scientific understanding of soil health and carbon cycling in agroecosystems.

3. Results and Discussion

3.1. Soil Orders

The Figure 3 presents the changes in topsoil SOC stock from 2010 to 2021 across various Soil Orders within Greek croplands. For 2010, Gleysols had the highest SOC stock (81.48 ±14.76 t C/ha st. dev, N=115) while Histosols had the lowest (51.97 ± 16.34 t C/ha st. dev, N=14). For 2021, SOC stock was highest in Histosols (73.62 ± 24.23 t C/ha st. dev, N=14) and lowest in Leptosols (43.63 ± 17.54 t C/ha st. dev, N=845). Histosols and Gleysols are the rarest soil orders in Greece (representing only the 1‰ and 7‰ respectively) and their size was not enough to accurately capture changes in these soil orders. The distribution of soil orders in Greek croplands are presenting in Figure 4.

In Figure 5, the values represent changes in SOC stock (ΔSOC), measured in metric tons of carbon per hectare (t C/ha), during the specified time within the corresponding soil types. Positive values indicate an increase in SOC stock, while negative values suggest a decrease. Post hoc tests revealed that 30 Soil Order pairs were significantly different, and 6 pairs were not significantly different. The following pairs were not significantly different: Cambisol-Gleysol, Cambisol-Luvisol, Fluvisol-Gleysol, Gleysol-Luvisol, Gleysol-Regosol and Histosol-Vertisol.

Calcisols (ΔSOC: -28.78 t C/ha, N=1076) and Leptosols (ΔSOC: -36.55 t C/ha, N=840) exhibited the most substantial decrease in SOC stock, indicating a negative change over the study period. This reduction may be attributed to factors such as changes in land management practices, erosion, or climatic and environmental conditions not favoring carbon sequestration.

Cambisols (ΔSOC: -15.69 t C/ha, N=3438) and Luvisols (ΔSOC: -16.21 t C/ha, N=902) also show a negative change in SOC stock, potentially influenced by factors like agricultural practices, soil structure, or changes in land use. Fluvisols (ΔSOC: -5.23 t C/ha, N=2293), Gleysols (ΔSOC: -7.59 t C/ha, N=109) and Regosols (ΔSOC: -11.71 t C/ha, N=389) indicate a negative change in SOC stock though to a lesser extent compared to Calcisols and Cambisols. The moderate decrease might be influenced by factors such as water dynamics, or land management practices. Drainage of wetland areas associated with Gleysols may be responsible for soil organic carbon loss.

Finally, Vertisols (ΔSOC: 3.65 t C/ha, N=305) and Histosols (ΔSOC: 19.13 t C/ha, N=14) show an increase in SOC stock, indicating gain of organic carbon. Factors such as water holding capacity and changes in land use may contribute to the observed increase in SOC stock.

Overall, at country level, ΔSOC stock was -15.02 (± 22.66 st. dev) t C/ha, N=9366. [39] found SOC stock change ranging between -5 and 1.5 t C/ha for countries in the Mediterranean during 2009-2018. Soil type affects soil organic carbon, Umbrisols and Podzols had higher organic carbon content than Leptosols and Fluvisols [38]. Soil type, pedogenic information and SOC depth distribution should be included in SOC inventory studies [40].

3.2. Soil Depth

Figure 6 presents the changes in topsoil (0-30 cm) SOC stock from 2010 to 2021 across soils with various Soil Depths within Greek croplands. In Figure 7 ΔSOC stock is presented across soils with different depths. In the post hoc tests, significant differences were observed in 12 pairs, and not significant differences in 3 pairs. The following pairs were not significantly different: (0-15, 15-30), (0-15, 30-60) and (100-150, >150). There is a substantial decrease in topsoil SOC stock in soils with 0-15 cm depth (ΔSOC: -28.32 t C/ha, N=97), 15-30 cm depth (ΔSOC: -33.88 t C/ha, N=603) and 30-60 cm depth (ΔSOC: -25.46 t C/ha, N=1354). Such a decrease can be attributed to various factors, including changes in land management practices, vegetation cover, or climatic conditions that do not favor carbon sequestration in the upper soil layers. Deeper soils lost carbon to a lesser extent, 60-100 cm depth (ΔSOC: -14.91 t C/ha, N=1986), 100-150 cm depth (ΔSOC: -10.95 t C/ha, N=1961) and > 150 cm depth (ΔSOC: -10.41 t C/ha, N=3365). These findings suggest that soil depth affects topsoil (0-30 cm) organic carbon. Shallow soils lost more topsoil carbon compared to topsoil carbon lost in soils with deep profiles. The topsoil in shallow soils could be more vulnerable to carbon loss due to differences in root activities [41] and decomposition dynamics [42].

3.3. Soil Coarse Materials

Figure 8 presents the SOC stock from 2010 to 2021 across different percentages of Coarse Materials within Greek croplands. In Figure 9 ΔSOC stock is presented across soils with different coarse material percentages. Post hoc tests revealed that all coarse material pairs (ΔSOC) were significantly different. Category 0 represents areas with no coarse materials, ΔSOC stock in these areas from 2010 to 2021 was -22.56 t C/ha, N=746. In areas where the content of coarse materials is less than 20%, ΔSOC SOC stock was -26.15 t C/ha, N=2200. This suggests that soil with a moderate presence of coarse materials experienced a more substantial decrease in SOC over the specified time. For areas with coarse material content ranging from 20% to 60%, ΔSOC was -9.8 t C/ha, N=6420. This category represents soils with a significant proportion of coarse materials, and the change in SOC stock was comparatively lower than in areas with less coarse material. Generally, soils with no coarse materials and soils with less than 20% coarse materials have experienced higher decreases in SOC compared to soils with higher coarse material content (20% to 60%). This information is valuable as it indicates that the texture of the soil, particularly the presence of coarse materials, plays a role in influencing SOC dynamics over the specified time frame.

3.4. Soil Texture

Figure 10 presents the SOC stock from 2010 to 2021 across different Soil Texture classes, Figure 11 presents ΔSOC (SOC₂₀₂₁-SOC₂₀₁₀). Five out of six pairwise comparisons were significantly different. Differences between moderately coarse and fine soils were not significantly different. The negative ΔSOC value of -13.07 t C/ha, N=2741 suggests a decrease in SOC content within soils with fine texture from 2010 to 2021. The relatively high ΔSOC value of -18.82 t C/ha, N=4505 indicates a substantial decrease in SOC stock within soils of medium texture over the specified timeframe.

Moderately coarse-textured soils, containing a higher proportion of sand, tend to have faster drainage and leaching rates. The negative ΔSOC value of -10.42 t C/ha, N=1082 suggests a decrease in SOC content, possibly influenced by factors such as land management practices and vegetation cover. The negative ΔSOC value of -7.51 t C/ha, N=1026 indicates a modest decrease in SOC stock within soils of coarse texture, reflecting the complex interplay of factors influencing carbon dynamics. A negative trend in SOC stock changes across all soil texture categories from 2010 to 2021. The magnitude of change varies, with soils of medium texture exhibiting the highest decrease, followed by moderately fine texture, fine texture, and moderately coarse texture. Coarse-textured soils show a more modest decrease.

Soils with medium texture lost more carbon compared to finer and coarser soils. On the contrary, fine soils and medium texture soils were expected to retain more carbon compared to coarser soils. A positive correlation has been found between SOC, clay and silt, while a negative correlation for SOC and sand [43]. Soils with fine texture retain more water than coarser soils [44] influencing microbial activity. Soils with moderately fine and fine textures typically have a higher surface area due to smaller particle sizes (clay and silt), this can enhance organic matter retention [45].

3.5. Soil Cation Exchange Capacity

Figure 12 presents the SOC stock from 2010 to 2021 across different ranges of Cation Exchange Capacity (CEC). CEC is a crucial parameter for soils, representing the soil's ability to retain and exchange cations, including nutrients like calcium, magnesium, and potassium. The values in the Figure 13 represent changes in SOC stock (ΔSOC) for different CEC ranges. Post hoc tests revealed that CEC pairs (3-8, >16) and (8-16, >16) were significantly different while the pair (3-8, 8-16) was not significantly different.

The first category represents soils with a CEC ranging from 3 to 8. Change in SOC stock for these soils over the specified period was ΔSOC: -2.17 t C/ha, N=170. Soils with a CEC between 8 and 16 fall into the second category. Change in SOC stock for this range is higher, ΔSOC: -4.8 t C/ha, N=1863. Next category includes soils with a CEC exceeding 16. For these soils, change in SOC stock was the highest among the presented categories, ΔSOC: -9.92 t C/ha, N=3961. CEC is affected by many factors such as clay, pH, and organic carbon [46]. Understanding the relationship between CEC and SOC dynamics has implications for sustainable land management practices. Soils with higher CEC may be less resilient to organic carbon loss.

3.6. Soil Drainage

Figure 14 presents the SOC stock from 2010 to 2021 across different drainage classes. In contrast to other soil variables, most post hoc tests (14 out of 21) did not show significant differences between ΔSOC in different drainage classes. The following pairs (7 out of 21) showed significant differences: (good drainage - moderate drainage), (good drainage - very good drainage), (good drainage – very poor drainage), (moderate drainage – very good drainage), (poor drainage – very good drainage), (perm. water table <50 – very good drainage), (very good drainage – very poor drainage).

Excessively drained soils typically have a rapid drainage rate, which can influence SOC dynamics. The substantial decrease in SOC stock (ΔSOC: -26.42 t C/ha, N=4329) in Figure 15 suggests a reduction in SOC stock in excessively drained soils, as expected. Well-drained soils exhibit good drainage characteristics, facilitating aeration and microbial activity. The negative ΔSOC value (-8.42 t C/ha, N=2430) indicates a decrease in SOC stock within well-drained areas from 2010 to 2021. Moderately drained soils strike a balance between drainage and water retention. The moderate decrease in SOC stock (ΔSOC: -3.79 t C/ha, N=1131) suggests a negative trend in carbon sequestration within these areas. The interplay between drainage and water availability may contribute to the observed changes in SOC over the specified timeframe. Poorly drained soils often experience waterlogging, influencing microbial activity and slowing organic matter decomposition [47]. The negative ΔSOC value (-3.53 t C/ha, N=205) indicates a decrease in SOC stock within poorly drained areas. Changes in drainage conditions, land use, or other environmental factors may contribute to the observed reduction in SOC levels. Very poorly drained soils typically suffer from prolonged water saturation, impacting SOC dynamics. The negative ΔSOC value (-3.4 t C/ha, N=1105) suggests a decrease in SOC stock within very poorly drained areas. These changes reflect the complex interplay of drainage conditions, land management practices, and environmental factors influencing the carbon dynamics in different soil types. The results underscore the importance of considering soil drainage characteristics when assessing SOC changes and implementing sustainable land management strategies.

3.7. Soil Slope

Figure 16 provides a snapshot of SOC stock and Figure 17 of SOC stock changes (ΔSOC) across different slope categories. The SOC stock values denote the amount of organic carbon present in the soil within each slope category. The highest SOC stock in 2010 was observed in 35-50% slope (93.17± 7 t C/ha st. dev, N=11) and lowest in 0-2% slope (56.47 ± 15.48 t C/ha st. dev, N=2808). In 2021, highest SOC stock was observed in 0-2% slope (55.82 ± 13.1 t C/ha st. dev, N=2824) and lowest in 25-35% (45.54 ± 15.35 t C/ha st. dev, N=35). In the post hoc tests, 13 ΔSOC stock pairs were significantly different and 8 pairs were not significantly different. The following slope pairs were not significantly different: (12-18%, 18-25%), (12-18%, 25-35%), (12-18%, 35-50%), (18-25%, 25-35%), (18-25%, 35-50%), (25-35%, 35-50%), (25-35%, 6-12%) and (35-50%, 6-12%).

The ΔSOC values increase as the slope becomes steeper, suggesting a potential negative correlation between slope steepness and SOC stock. Steeper slopes are often more susceptible to erosion, which can lead to the removal of carbon in the topsoil and, consequently, a decrease in SOC stock [48]. SOC content has been found to decrease with increasing slope gradient [49]. However, the slope effect on SOC is not always pronounced and significant, it also depends on land use [50]. Changes in land use, agricultural practices, or afforestation efforts may influence SOC dynamics. For instance, implementing conservation practices on steeper slopes could mitigate SOC stock reduction.

3.8. Soil pH

Figure 18 reflects the SOC stock from 2010 till 2021 and Figure 19 the changes in SOC stock (ΔSOC) associated with varying pH levels in Greek croplands. In the post hoc tests, 8 ΔSOC stock pairs were significantly different and 2 pairs were not significantly different. The following pH pairs were not significantly different: (4.5-5.5, 5.6-6.9) and (8-8.5, >8.5).

In the pH range between 4.5-5.5, a small increase in SOC stock was observed, ΔSOC: + 2.35 t C/ha, N=129. This indicates a relatively stable trend in SOC content within this pH category. Croplands with a pH between 5.6-6.9 experienced a small decrease in SOC stock, ΔSOC: -1 t C/ha, N=677. Moderately alkaline (pH 7-7.9) and alkaline (pH 8-8.5) soils exhibited a substantial decrease in SOC stock, ΔSOC: -7.39 t C/ha, N=2284 and ΔSOC: -10.95 t C/ha, N=2777. Highly alkaline soils, with a pH greater than 8.5, had the largest SOC stock decrease -17,37 t C/ha, N=131. The trends suggest that acidic soils retained more soil carbon compared to alkaline soils. That could be explained by the effect of pH on microbial activity (and decomposition) or changes land use intensity. [51] observed a SOC stock increase at pH 4.2-6.5 and decrease at pH 6.5-9.2. Positive correlations between SOC and pH have also been found [52], meaning that the relationship between SOC and pH is context dependent.

5. Conclusions

The comprehensive investigation into the spatio-temporal dynamics of SOC stock in Greek croplands has yielded valuable insights, shaping our understanding of the intricate relationships governing SOC sequestration and turnover. This study, guided by a multidimensional exploration of various soil properties, ranging from pH to soil texture, has provided nuanced perspectives on the factors influencing SOC dynamics within the unique context of Greek agriculture. SOC loss patterns differ across soil order, soil depth, coarse materials, soil texture, CEC, drainage, slope, and pH categories. The study highlights the critical need for a particular examination of this dynamic interplay to strike a balance between agricultural productivity and soil health. The classification bestowed by soil order emerges as a key determinant, weaving a complex tapestry of SOC variability. Various microscopic factors, from soil texture to pH, play pivotal roles in SOC dynamics. In addition, landscape factors such as slope play a very important role in soil organic loss rate. These insights emphasize the significance of understanding the regional intricacies and tailoring soil management strategies accordingly.

The relationship between anthropogenic activities and SOC stock reveals the socio-economic fabric of the Mediterranean, particularly in Greece. Although croplands significantly contribute to the regional socio-economic tapestry, human interventions and land management practices wield transformative influences on SOC stock, necessitating an examination of this dynamic interplay.

Lastly, this scientific endeavor not only synthesizes existing knowledge but unveils novel insights into the spatio-temporal dynamics of SOC stock in Greek croplands. The narrative woven from diverse soil properties provides a robust foundation for sustainable land management practices and climate change mitigation strategies tailored to the unique context of Greek agriculture. The findings underscore the importance of localized approaches in addressing global challenges, ensuring that strategies are both effective and adaptable to the specificities of regional ecosystems.