Next Article in Journal
Water Saving and Yield of Potatoes under Partial Root-Zone Drying Drip Irrigation Technique: Field and Modelling Study Using SALTMED Model in Saudi Arabia
Next Article in Special Issue
Cover Crop Contributions to Improve the Soil Nitrogen and Carbon Sequestration in Almond Orchards (SW Spain)
Previous Article in Journal
Do We Need New Crops for Arid Regions? A Review of Fruit Species Domestication in Israel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geostatistical Analysis of Soil C/N Deficiency and Its Effect on Agricultural Land Management of Major Crops in Eastern Croatia

1
Faculty of Agrobiotechnical Sciences Osijek, Josip Juraj Strossmayer University of Osijek, Vladimira Preloga 1, 31000 Osijek, Croatia
2
Department of Biology, Josip Juraj Strossmayer University of Osijek, Cara Hadrijana 8/A, 31000 Osijek, Croatia
3
Faculty of Geodesy, University of Zagreb, Kačićeva 26, 10000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
Agronomy 2020, 10(12), 1996; https://doi.org/10.3390/agronomy10121996
Submission received: 29 November 2020 / Revised: 15 December 2020 / Accepted: 17 December 2020 / Published: 18 December 2020

Abstract

:
Interpolation accuracy of the soil carbon-to-nitrogen ratio (C/N) has a direct impact on soil maps for agricultural land management planning, affecting the financial gains of farmers and environmental contamination. This study aimed to evaluate the effect of different ordinary kriging parameters on the interpolation accuracy and spatial variability of soil C/N. These values were necessary to determine the effect of potential soil C/N deficiency on five major crops in the study area during 2017–2019, containing maize, wheat, sunflower, rapeseed and soybean parcels. The Gaussian model with logarithmic transformation produced the higher outer accuracy, representing the prediction accuracy. Average values of final soil C/N rasters resulted as 11.35 at 0–10 cm soil depth and 12.60 at 20–30 cm soil depth, representing values 11.0% higher at a deeper soil layer. The moderate soil N deficiency was observed for all major crops, having soil C/N values mainly higher than 12:1, compared to the optimal value range of 9:1 to 12:1. These observations suggest a need for the adjustment of agricultural land management plans through the regionalization of agricultural production to prevent further soil degradation and ensure optimal crop development conditions.

1. Introduction

The cultivation of agricultural soils inevitably accelerates the processes of humus decomposition, which results in the reduction of organic matter content [1]. The rate of organic matter content reduction is dependent on the land management system, and every agrotechnical operation should consider the amount of soil organic matter [2]. Soil organic matter represents a major source of most nutrients and consists of significant quantities of carbon (C) in combination with smaller quantities of nitrogen (N) [3]. The soil carbon-to-nitrogen ratio (C/N) has a direct impact on the agricultural residues’ degradation, as well as on the circulation of N in soils [4]. Climate changes with the increase of annual air temperatures and lower precipitation indirectly stimulate soil organic carbon mineralization, obstructing the input of soil organic carbon and gradually decreasing its quantity [5]. Soil C/N is one of the primary indicators of these effects, reflecting the severity of soil degradation as the lower capacity of soil for transfer of important nutrients to vegetation and microorganism development [6]. The intensive agricultural production caused by the growing necessity for higher crop quantity demands invariably affects dynamics of carbon and nitrogen, endangering the sustainability of present cropping systems [7]. Recent research by Bogunović et al. [8] supported the claim that the intensive agricultural production negatively influences these dynamics, specifically reducing soil organic matter and soil organic carbon stocks. To ensure the sustainability of agricultural production in such areas, fertilization, weed control and soil tillage methods should be included in agricultural land management plans, as they directly or indirectly affect long-term soil degradation [9,10]. Furthermore, the regulation of soil C/N is necessary for optimal agricultural land management and the regionalization of crop production, having a major impact on crop development and growth [11]. This is especially relevant for organic farming, as it mainly depends on the circulation of soil organic matter through crop residues within an agricultural parcel [12]. Crop yield produced by organic farming is generally lower comparing to the conventional approach, which reinforces the need for soil C/N management for sustainable organic farming [13]. As per Eurostat data [14], the average percentage of total agricultural area utilized for organic farming was 7.50% in 2019, which steadily increased compared to 5.64% in 2012. There is still a large heterogeneity of organic farming policies within the European Union (EU), as these percentages during 2018 vary from 2.43% for Romania and 24.08% for Austria. Croatia in among the EU countries with the highest trend of organic farming expansion, growing from 2.40% in 2012 to 6.94% in 2018 of total agricultural area utilized for organic farming and is expected to surpass the EU average in forthcoming years.
Reliable soil maps are fundamental decision-making bases for soil C/N management [15]. A Geographic Information System (GIS) combines the tools for collecting and processing the field samples, their processing by geostatistical methods, geospatial analysis and the creation of thematic maps [16]. These data are commonly implemented in various suitability analyses as a part of GIS-based multicriteria analysis studies for agricultural and environmental purposes [17]. In soil analysis, it is possible to collect only a limited number of samples, which in such form do not give a full representation of soil properties. The spatial variability of these soil values should also be considered, which indicates a necessity for spatial interpolation methods [18]. Geostatistical methods enable the interpolation of values between sample points used for soil map creation [19]. Soil maps are vital for precision farming, and their quality depends on the quality and density of sample points and interpolation parameters [20]. Kriging was proven to be superior in most cases to deterministic interpolation techniques for interpolation of soil chemical properties [21]. Ordinary kriging (OK) is the most commonly used kriging method and was successfully applied for the interpolation of various soil properties [22,23]. Normality and stationarity of input data are two requirements for the interpolation using OK [24]. These properties, as well as other properties like the number of samples and value ranges, determine the optimal interpolation parameters and interpolation accuracy [25]. Therefore, an evaluation of interpolation parameters is necessary for accurate creation of soil maps [26].
The hypothesis of the study is that the properties of input samples and interpolation parameters have a major impact on interpolation accuracy and the predicted spatial variability of soil C/N using OK. The specific objectives of the study were to: (1) evaluate normality, stationarity and other relevant properties of input samples, (2) perform spatial interpolation using OK with different interpolation parameters and assess their accuracy, and (3) evaluate the areas of soil C/N deficiency and determine its effect on major crops in the study area.

2. Materials and Methods

2.1. Study Area and Data Acquisition

The study area was the County of Osijek-Baranja, situated in eastern Croatia (Figure 1). Agricultural land covers 2082 km2 (50.1%) of the total 4155 km2 county area, according to the Paying Agency for Agriculture, Fisheries and Rural Development data for the year 2019. As reported by the Croatian Bureau of Statistics in 2016, 14% of Croatia’s agricultural area is being utilized in the County of Osijek-Baranja, this being the traditionally first-ranked county in agricultural production. According to the same source, more than 20% of cereals, 46% of sugar beet and 17% of orchards from the total country’s agricultural area were cultivated in County. The recent crop suitability studies in the study area [25,27] showed that the current agricultural land management plans could be upgraded to allow considerably better regionalization of agricultural production for individual crops. The update of obsolete agricultural land management according to climate changes is crucial for a stable crop yield and the potential export of high-quality crops [28]. The traditionally intensive conventional tillage in eastern Croatia, with ploughing as the main tillage procedure, mainly affects soil organic matter [29]. As per the European soil database v2.0 (https://esdac.jrc.ec.europa.eu/content/european-soil-database-v20-vector-and-attribute-data), orthic luvisol and calcaric fluvisol are the most common soil types, covering 59.3% of the study area combined. A moderately warm and rainy climate type class ‘Cfwbx’ is present in Osijek-Baranja County, classified according to a Köppen scale [30]. Historical climate data (1899–2018) for the weather station Osijek near the centroid of the county showed an average annual air temperature of 11.1 °C and the average total annual precipitation of 692.9 mm [31]. The average summer air temperature calculated from June to September during the same period was 19.8 °C.
Four sample data sets were obtained from the Croatian Agency for the Environment and Nature soil data WFS service, collected in the field during the year 2016. Pre-treatment of soil samples was conducted according to ISO 11464:2006, while the elemental analysis using dry combustion for soil C and N was performed according to ISO 10694:1995 and ISO 13878:1998. A total of 31 samples were collected in the agricultural land on 0–10 cm soil depth for C (C 0–10 cm) and N (N 0–10 cm) and 64 samples on 20–30 cm soil depth for C (C 20–30 cm) and N (N 20–30 cm). These samples were collected as the supplementary data to the national soil sampling in 2013, causing an uneven sampling distribution at the 0–10 cm and 20–30 cm soil layers. Soil sampling at the 10–20 cm soil layer was omitted from the process, which is the reason why it is not contained in the input dataset. Soil sample values labeled as ‘C’ designate soil organic carbon and ‘N’ designate total soil nitrogen. These values were successfully implemented for the evaluation of soil C/N properties in [32,33]. The measurement unit for both soil C and N was g 100 g−1. These sampling depths were selected as soil C/N generally remains consistent with the increase of soil depth and is mainly influenced by the crops, while about of 80% of the soil organic carbon is located at the 0–40 cm soil depth [34].
Agricultural parcel data was received by the Croatian Paying Agency for Agriculture, Fisheries and Rural Development from 2017, 2018 and 2019 for five major crops in Osijek-Baranja County. Maize, wheat, sunflower, rapeseed and soybean consistently amounted to about two thirds of total county cultivated crop area, having a major importance on the county’s economy. Previous studies proved that soil C/N values indicate the yield of maize and wheat [32], sunflower [35], rapeseed [36] and soybean [37]. Annual cultivated area and the part of total county agricultural area per crop are displayed in Table 1. The agricultural parcels used were visualized using three example subsets, as shown in Figure 2. These example subsets cover three representative areas regarding the land-cover and land-use systems: a western part consisting of fragmented agricultural parcels surrounded by forests, a central part characterized by the intensive agricultural production, and an eastern part in Baranja region with large parcels of a major agricultural company in the proximity of wetland. The total cultivated agricultural area in the study area was 209,874 ha in 2017, 209,702 ha in 2018 and 208,163 ha in 2019.

2.2. Geostatistical Analysis

The geostatistical analysis was conducted in three steps: evaluation of normality and stationarity of input samples, evaluation of properties and accuracy of the interpolation results, and interpolation of soil C/N with the assessment of spatial variability.
Evaluation of input samples consisted of tested for data normality and stationarity, as required assumptions for kriging [24]. The normality test was used for the assessment of sample values distribution, whereas stationarity was used for the assessment of spatial variability distribution. Descriptive statistics consisting of the mean, minimum, maximum, coefficient of variation (CV), skewness and kurtosis were used for the evaluation of data normality and stationarity. The distribution of data was additionally examined using the Shapiro–Wilk normality test [38]. A borderline p-value of 0.05 was used for the evaluation of the null hypothesis of normal data distribution with 95% confidence. The Shapiro–Wilk test was performed using Python version 3.7.3 (Scotts Valley, CA, USA) with library SciPy version 1.3.0. Stationarity (Austin, TX, USA) was evaluated by Thiessen polygons, indicating differences in local variations of the sample data. Thiessen polygons were also used to indicate locations of sparse sample data, having notably larger area polygons than others.
Interpolation was performed using the OK method, which is based on the calculation of spatial autocorrelation using the variogram calculated according to the following formula (Equation (1)) [19]:
γ ( h ) =   1 2 m ( h ) i = 1 n ( h ) [ z ( x i )     z ( x i + h ) ] 2 ,
where γ ( h ) is the semivariance, z ( x i ) and z ( x i + h ) are the input values at locations x i and x i + h , n(h) is the number of paired comparisons at distance h. Evaluated mathematical models for the fitting to empirical variogram were linear, square root, Gaussian and spherical [39]. The formula for the prediction of values based on OK method was (Equation (2)) [19]:
z ( x ) =   i = 1 n λ i z ( x i ) ,
where z   ( x ) is predicted value at location x and λ i is the weight calculated using the variogram. The interpolation using logarithmic transformation was evaluated in the research, being recommended in the absence of normal distribution of input data [40]. Each of the four mathematical models was evaluated with and without logarithmic transformation for the quantification of the impact of a normal distribution on spatial variability. A total of 9 lags were used in the interpolation for 0–10 samples, while 12 lags were applied for the interpolation using 20–30 cm samples. Geostatistical calculations were conducted in SAGA GIS version 7.3.0 (Hamburg, Germany). The spatial resolution of all interpolation results was set to 250 m.
Nugget (n), sill (s) and range (r) values were calculated for each iteration using different OK interpolation parameters for the determination of interpolation characteristics. The degree of spatial dependence was calculated for an indication of the spatial autocorrelation strength. Soil C/N rasters were calculated by dividing interpolated C and N rasters of the corresponding interpolation parameters. By combining the four evaluated input samples, mathematical models and logarithmic transformation, a total of 16 interpolation variants using OK were calculated. Descriptive statistics consisting of mean, minimum, maximum and CV were calculated for the final soil C/N values to quantify value ranges and spatial variability produced by various parameters.

2.3. Accuracy Assessment of Interpolated Results

Coefficients of determination were used for the accuracy assessment, representing inner ( R I 2 ) and outer accuracy ( R O 2 ). This procedure was successfully applied in a study by Jurišić et al. [26]. Inner accuracy showed the correlation between input sample values and interpolation results on an identical location. This value quantified the preservation of the ground-truth values in the interpolation results. Outer accuracy represented the prediction accuracy of the particular interpolation parameters at an unknown location. Outer accuracy was evaluated using a cross-validation method with a 2-fold procedure. The calculation was performed by the interpolation with missing random two samples per iteration and repeating the same procedure until all of the input samples were evaluated.

2.4. Soil Carbon-to-Nitrogen Ratio (C/N) Deficiency Evaluation for Five Major Crops during 2017–2019

Soil C/N interpolation results which produced the highest outer accuracy were determined as representative and selected for further analysis. The evaluation of soil C/N dynamics for five major crops was performed for 2017 as soil sampling was conducted the previous year. This was achieved by the calculation of mean and standard deviation values for each of five major crops, covering all agricultural parcels within the study area of a particular crop. The analysis was performed individually for 0–10 cm and 20–30 cm soil depths, as the 0–10 cm soil layer is particularly susceptible to soil degradation and reflects the areas of the most intensive agricultural production [41]. Based on the same soil C/N interpolation results, the projections for 2018 and 2019 were created using the same approach to determine the suitability of current crop management plans according to the present state for soil degradation. The soil C/N ratios in ranges of 10:1 to 12:1 are regarded as the standard values for the agriculturally important cultivated soils, allowing an optimal balance of carbon oxidation process and nitrogen immobilization [42]. The soil C/N of intensely cultivated soils is generally expected to produce lower soil C/N compared to the uncultivated soils, based on the same source. The threshold values of soil C/N for the balance of these processes differed from 9.7 and 10.8, depending on the soil texture, C and N soil contents [43]. Marchetti et al. [44] accordingly proposed the soil C/N classification in three classes: optimal conditions (9 ≤ C/N ≤ 12), dominant oxidation (C/N < 9) and nitrogen deficiency (C/N > 12). The severity of soil degradation according to five major crops was evaluated according to the soil C/N classification in these classes.

3. Results

Higher input mean values for both C and N were observed in the top soil layer (Table 2). All samples possessed a moderate variability, with an average value of 0.53. High skewness values indicated a non-normal value distribution for all samples. Kurtosis values showed a platykurtic distribution, except for N 20–30 cm sample. Shapiro–Wilk test results suggested the absence of normal distribution of values for all input samples, as all of the p-values were below the border value. C 0–10 cm sample had the closest values distribution to normal, having the lowest skewness value and the highest p-value of observed samples. Both N samples produced very low W- and p-values.
Thiessen polygons suggested that all input samples possessed a moderate stationarity. The mild transition of neighboring colors indicated the lower local variation between the two polygons. Minor anomalies in the properties of Thiessen polygons were observed in the central area of the County, having a high local variation (Figure 3). The highest amounts of C and N are generally grouped on the eugley soils in the northeastern part of the study area. Areas with mainly cultivated agricultural land in the central part had the lowest values. Polygons for samples in 20–30 cm soil depth have an approximately even area through the County, indicating even uniform spatial distribution of samples. For samples in 0–10 cm soil depth, polygons have varying areas, the largest being in the central part.
Evaluation of properties and accuracy of the interpolation results produced by various interpolation parameters was performed according to statistical values for 0–10 cm (Table 3) and 20–30 cm samples (Table 4). Interpolation with no transformation produced the higher spatial dependence coefficient compared to interpolation with logarithmic transformation in 12 out of 16 interpolation variants. All four cases with a higher spatial dependence using logarithmic transformation occurred in highly skewed N samples. An overall average spatial dependence using logarithmic transformation was considerably higher for N samples (0.888), compared to C samples (0.572). A similar tendency occurred for both inner and outer accuracies since N samples produced a higher average R I 2 of 11.5% and average R O 2 of 9.2% compared to C samples with the application of logarithmic transformation. The difference of average spatial dependence using logarithmic transformation between N and C samples of 0.316 is considerably higher than in the case of using no transformation (0.178). Logarithmic transformation produced lower R I 2 values compared to no transformation, with average respective values of 78.4% and 87.4% for C samples, alongside 89.9% and 93.1% for N samples. However, the application of logarithmic transformation allowed a more accurate prediction of variables, having higher R O 2 than no transformation of 8.5% for C and 22.9% for N samples. A spherical model produced a higher average R I 2 compared to other mathematical models, both observing the chemical properties (86.9% for C and 96.5% for N) and soil sampling depths (97.4% for 0–10 cm and 86.0% for 20–30 cm). Similarly to the spherical model for inner accuracy, the Gaussian model resulted in the highest average R O 2 for the chemical properties (56.7% for C and 61.9% for N) and soil sampling depths (58.5% for 0–10 cm and 60.1% for 20–30 cm). The effect of sparser sampling density was present in the outer interpolation accuracy values, with interpolation results of 20–30 cm samples having 3.2% higher average R O 2 compared to results of the 0–10 cm samples. The variograms of mathematical models that produced the highest inner and outer accuracy per input sample are shown in Figure 4.
Figure 5 shows the soil C/N rasters from the 16 variants using different interpolation parameters. An average soil C/N interpolated for 0–10 cm soil depth resulted in 11.35, while the average value for 20–30 cm soil depth was 12.60 (Table 5). The interpolation results with a logarithmic transformation had larger average value ranges, but a smaller average CV compared to no transformation. Gaussian model produced the highest value ranges and CV, representing high local variability in the interpolation results. The linear and spherical models performed similarly in terms of spatial variability, smoothing minor variabilities in soil C/N and having a CV close to 0.25. The square root model retained the lowest amount of local spatial soil C/N variability with the mean CV of 0.20. Final soil C/N results at 0–10 cm soil depth produced more variability than 20–30 cm soil depth, with a 0.05 higher average CV.
The combined scatterplots represent the correlation between interpolation results calculated using the same interpolation method. Scatterplots with Pearson’s correlation coefficient are presented in the bottom-left part of the figures with no transformation, while the results of interpolation with logarithmic transformation are located in the upper-right part. Histograms of particular interpolation results are located in the middle of the figure diagonally. Combined scatterplots for 0–10 cm and 20–30 cm soil depths are shown in Figure 6 and Figure 7, respectively. Linear and square root mathematical models produced the highest average correlation of all four interpolation sets, followed by the combination of linear and spherical mathematical models. Gaussian model produced the distinctive results, having an average correlation coefficient to other interpolation results as low as 0.173 for 0–10 cm samples with the application of logarithmic transformation.
The soil C/N balance per major crops was evaluated according to the values displayed in Figure 8. Values from 2017 indicate a close distribution of soil C/N at 0–10 cm soil layer for maize, wheat, sunflower and rapeseed, with the mean C/N values of 12.71, 13.25, 12.42 and 12.56 to 1, respectively. Soil C/N values for soybean largely deviated from these observations, having a mean of 14.40 to 1, alongside the largest heterogeneity with the standard deviation of 7.26. These values are generally above the limit of 12:1, indicating a moderate soil N deficiency. Projections for soil C/N at 0–10 cm soil layer during 2018 and 2019 remained consistent with the 2017 values, with a slight gradual decrease of mean soil C/N for sunflower and rapeseed. Maize, wheat and rapeseed mean soil C/N values at the 20–30 cm soil layer from 2017 retained the same trend as 0–10 cm, having similar values with slightly lower standard deviation. The largest difference resulted for soybean, having a 2.62 lower mean soil C/N value at the higher soil depth. These values retained the same lower projected trend during 2018 and 2019, having the projected mean soil C/N of 11.45 to 1, which is a decrease of 0.33 compared to 2017 values. By contrast, sunflower resulted in the highest mean soil C/N value of 13.12 to 1 at the 20–30 cm soil layer during 2017 with similar projections for 2018 and 2019.

4. Discussion

The interpolation accuracy and spatial variability of the predicted soil C/N were conditioned by the properties of input samples as noted in [45], confirming the hypothesis of the study. The higher spatial dependence was reflected in higher inner and outer accuracy values, based on the average spatial dependence of C and N input samples. While the application of no transformation allowed higher spatial dependence for the majority of input samples, in the case of the very low normality and stationarity, logarithmic transformation produced better results. The necessity for the objective quantification of data stationarity appeared during the evaluation of the properties of input samples, as the interpretation of Thiessen polygons heavily relies on the human subjective assessment. While logarithmic transformation had a worse effect on the preservation of ground-truth data in the interpolation results compared to no transformation, it produced considerably higher prediction accuracy, which increased with the skewness of input data. Consequentially, more different transformation types offer a potential of improving the interpolation accuracy of OK [46]. The applied procedure of the evaluation of the properties of input samples supports automatization to reduce the processing time, as every step besides stationarity evaluation was quantitatively and objectively expressed.
Kriging has been successfully applied in soil-related studies with a similar number of input samples but restricted to a county- or national-scale mapping. Schloeder et al. [47] used OK with 46 soil samples, having low stationarity and irregular sample spacing. Ding and Yu [48] successfully applied kriging on 68 samples, in combination with remote-sensing data for estimation of soil salinity. Testing of deterministic interpolation methods, such as inverse distance weighted, might be beneficial for highly skewed input samples with low sampling density [26]. The application of no transformation and a spherical mathematical model produced the highest average inner interpolation accuracy. This was reflected in the largest retention of ground-truth values from input samples in the interpolation results. Outer accuracy received much more attention in recent studies than inner accuracy [22,23,38,40], indicating an accuracy for the prediction of values at an unknown location. Logarithmic transformation and the Gaussian model were the best-performing parameters considering outer accuracy. The Gaussian model also resulted in the highest accuracy in a similar study by Fashi et al. [49], where spherical, exponential, linear and Gaussian models were tested. Besides the highest outer accuracy, the Gaussian model allowed the highest flexibility in fitting to the empirical variogram, having consistently high fitting correlation values in a wide range. These observations were made by examining a generally non-normal data with moderate stationarity, so additional evaluation using different input samples is necessary.
Average soil C/N values from the interpolation results in this study differ with the observations that soil C/N of most soils decreases with soil depth [50,51]. The lowest soil C/N values on both soil depths were observed in the central part of the study area, which is characterized by intensive agricultural production. The highest values resulted mainly in the northeastern part of the study area for both soil depths, with more heterogeneity in the case of 20–30 cm soil depth. Soil sampling density and the application of logarithmic transformation caused by varying input samples distribution had the highest effect on the spatial variability of soil C/N in the interpolation results, as noted in [27]. Evaluated mathematical models produced contrasting interpolation results, primarily in terms of the local variabilities of soil C/N. The Gaussian model, which had the highest average outer accuracy, produced the highest mean soil C/N for both soil depths. It also resulted in the highest value range and CV, which translated in the highlighting of spatial variabilities. The square root model produced the smoothest interpolation results with the smallest value range, being suitable for the mapping on a very small scale. The spherical model highlighted local spatial variabilities for most of the evaluated mathematical models, retaining the ground-truth data near to those of input samples. The selection of a mathematical model for fitting to empirical variogram had the highest impact on spatial variabilities of soil C/N, reinforcing the need for the evaluation of various alternatives for the mapping of soil features [52].
The analysis of soil C/N values at 0–10 cm and 20–30 cm soil depths determined that the current soil C/N in the study area is generally above the 12:1 value, indicating a widespread N deficiency. The major crops are commonly located at the intensively cultivated parts of the study area, which showed the greater N deficiency compared to the less cultivated area. The area of the highest soil C/N deficiency in the central part of the study area especially coincides with the area of intensive agricultural production, having an increased risk of soil degradation. Uneven distribution of soil C/N at the different soil depths could produce unexpected long-term crop development and growth, most notably for soybean and sunflower. A main observation remains for the soil C/N values for soybean at the 20–30 cm soil layer, which showed a decrease in the projections for 2018 and 2019 and the presence of vertically stratified soil C/N in the 0–10 cm and 20–30 cm soil layers. Given the non-compliance of crop rotation systems based on the Croatian Paying Agency for Agriculture, Fisheries and Rural Development data were observed at the multiple locations in the study area during 2017–-2019. This is a possible cause of the lower soil C/N for soybean at 20-30 cm compared to the 0–10 cm soil layer, as well as for the large soil C/N value range for soybean, as high N soil content might occur if a legume crop was cultivated prior to soybean. These inconsistencies suggest the need for the adjustment of agricultural land management plans at the county level. This approach of regionalization of agricultural production could ensure sustainable organic farming and maximum possible yield, with possible financial gains for farmers and lower environmental contamination. One of the possible ways of accomplishing this is by the distribution of agricultural subsidies for individual crops according to their land suitability, which in itself ensures quality climate and soil conditions for crop development [53]. GIS-based multicriteria analyses with the analytic hierarchy process (AHP) method represent a reliable and flexible solution for cropland suitability calculation globally and could assist in resolving this problem [25,54]. In addition, the constant development of remote-sensing technology offers a reliable data source for the modeling of crop suitability through vegetation indices both at macro scales using satellite imagery [55] and micro scales using unmanned aerial vehicles [56].

5. Conclusions

The hypothesis of the study that the selection of interpolation parameters and the properties of input samples has a major impact on the interpolation accuracy and spatial variability was confirmed. Therefore, an evaluation process of interpolation parameters and accuracy assessment is recommended in similar studies. Non-normal values distribution and moderate stationarity were observed during the evaluation of all four input samples. Logarithmic transformation negatively affected spatial dependence values compared to no transformation, which results in lower average inner accuracy. However, it enabled a more accurate prediction of values at an unknown location, represented by outer accuracy. The spherical model produced the highest ability to retain ground-truth data in the interpolation results, while the Gaussian model enabled the highest outer accuracy. The Gaussian model allowed the highest flexibility of tested mathematical models for fitting to empirical variograms, having a wide range with a high fitting coefficient of determination. Interpolated soil C/N values resulted in a higher average value at 20–30 cm soil depth, compared to 0–10 cm soil depths. The higher variability of interpolated soil C/N was observed at a 0–10 cm soil depth. This observation was contrary to some recent studies that soil C/N lowers in a deeper layer of the soil. The Gaussian model produced the highest spatial variability of predicted soil C/N, also having the largest value interval. The spherical model enabled the highest local spatial variabilities in the interpolation results, retaining the input sample values in the proximity of samples. The analysis of soil C/N distribution for five major crops in the study area (maize, wheat, sunflower, rapeseed and soybean) indicated a moderate soil N deficiency. The soil C/N values generally resulted as higher than 12:1, which differs from the optimal soil C/N range from 9:1 to 12:1 and generally indicates a moderate soil degradation risk. These observations suggest a need for the modification of present agricultural land management plans through regionalization of agricultural production to prevent further soil degradation and to ensure the best possible conditions for crop development.

Author Contributions

Conceptualization, M.J., D.R. and M.G.; methodology, D.R.; software, D.R.; validation, M.G.; formal analysis, M.J., D.R., S.K. and I.P.; investigation, M.J., D.R. and I.P.; resources, M.J., S.K. and I.P.; data curation, D.R.; writing—original draft preparation, D.R. and M.G.; writing—review and editing, M.J., D.R., S.K., I.P. and M.G.; visualization, D.R.; supervision, M.J. and M.G.; project administration, M.J.; funding acquisition, M.J., D.R., I.P. and M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This work was supported by the Faculty of Agrobiotechnical Sciences Osijek as a part of the scientific project “AgroGIT—technical and technological crop production systems, GIS, and land conservation”. This work was supported by the University of Zagreb as a part of the scientific project: “Advanced photogrammetry and remote sensing methods for environmental change monitoring” (Grant No. RS4ENVIRO).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Scotti, R.; Bonanomi, G.; Scelza, R.; Zoina, A.; Rao, M.A. Organic amendments as sustainable tool to recovery fertility in intensive agricultural systems. J. Soil Sci. Plant Nut. 2015, 15, 333–352. [Google Scholar] [CrossRef] [Green Version]
  2. Peter, P.C. Soil Organic Matter/Carbon Dynamics in Contrasting Tillage and Land Management Systems: A Case for Smallholder Farmers with Degraded and Marginal Soils. Commun. Soil Sci. Plant Anal. 2017, 48, 2013–2031. [Google Scholar] [CrossRef]
  3. Puttaso, P.; Namanusart, W.; Thumanu, K.; Kamolmanit, B.; Brauman, A.; Lawongsa, P. Assessing the Effect of Rubber (Hevea brasiliensis (Willd. ex A. Juss.) Muell. Arg.) Leaf Chemical Composition on Some Soil Properties of Differently Aged Rubber Tree Plantations. Agronomy 2020, 10, 1871. [Google Scholar] [CrossRef]
  4. Truong, T.H.H.; Marschner, P. Presence of wheat straw in soil influences nutrient availability and leaching in soil mulched with high or low C/N organic materials. Arch. Agron. Soil Sci. 2020, 1–12. [Google Scholar] [CrossRef]
  5. Zhang, Z.C.; Hou, G.; Liu, M.; Wei, T.X.; Sun, J. Degradation induces changes in the soil C: N: P stoichiometry of alpine steppe on the Tibetan Plateau. J. Mt. Sci. 2019, 16, 2348–2360. [Google Scholar] [CrossRef]
  6. Dlamini, P.; Chivenge, P.; Manson, A.; Chaplot, V. Land degradation impact on soil organic carbon and nitrogen stocks of sub-tropical humid grasslands in South Africa. Geoderma 2014, 235, 372–381. [Google Scholar] [CrossRef]
  7. Aher, G.; Cihacek, L.J.; Cooper, K. An evaluation of C and N on fresh and aged crop residue from mixed long-term no-till cropping systems. J. Plant. Nutr. 2017, 40, 177–186. [Google Scholar] [CrossRef]
  8. Bogunovic, I.; Viduka, A.; Magdic, I.; Telak, L.J.; Francos, M.; Pereira, P. Agricultural and Forest Land-Use Impact on Soil Properties in Zagreb Periurban Area (Croatia). Agronomy 2020, 10, 1331. [Google Scholar] [CrossRef]
  9. Meseldžija, M.; Rajković, M.; Dudić, M.; Vranešević, M.; Bezdan, A.; Jurišić, A.; Ljevnaić-Mašić, B. Economic Feasibility of Chemical Weed Control in Soybean Production in Serbia. Agronomy 2020, 10, 291. [Google Scholar] [CrossRef] [Green Version]
  10. Dekemati, I.; Simon, B.; Bogunovic, I.; Kisic, I.; Kassai, K.; Kende, Z.; Birkás, M. Long Term Effects of Ploughing and Conservation Tillage Methods on Earthworm Abundance and Crumb Ratio. Agronomy 2020, 10, 1552. [Google Scholar] [CrossRef]
  11. McAllister, C.H.; Beatty, P.H.; Good, A.G. Engineering nitrogen use efficient crop plants: The current status. Plant. Biotechnol. J. 2012, 10, 1011–1025. [Google Scholar] [CrossRef] [PubMed]
  12. Leifeld, J. How sustainable is organic farming? Agric. Ecosyst. Environ. 2012, 150, 121–122. [Google Scholar] [CrossRef]
  13. Boone, L.; Roldán-Ruiz, I.; Muylle, H.; Dewulf, J. Environmental sustainability of conventional and organic farming: Accounting for ecosystem services in life cycle assessment. Sci. Total Environ. 2019, 695, 133841. [Google Scholar] [CrossRef] [PubMed]
  14. Eurostat, Organic Crop Area by Agricultural Production Methods and Crops (from 2012 Onwards). Available online: https://appsso.eurostat.ec.europa.eu/nui/submitViewTableAction.do (accessed on 28 November 2020).
  15. Nabiollahi, K.; Golmohamadi, F.; Taghizadeh-Mehrjardi, R.; Kerry, R.; Davari, M. Assessing the effects of slope gradient and land use change on soil quality degradation through digital mapping of soil quality indices and soil loss rate. Geoderma 2018, 318, 16–28. [Google Scholar] [CrossRef]
  16. Fu, W.; Zhao, K.; Tunney, H.; Zhang, C. Using GIS and geostatistics to optimize soil phosphorus and magnesium sampling in temperate grassland. Soil. Sci. 2013, 178, 240–247. [Google Scholar] [CrossRef]
  17. Gašparović, I.; Gašparović, M. Determining Optimal Solar Power Plant Locations Based on Remote Sensing and GIS Methods: A Case Study from Croatia. Remote Sens. 2019, 11, 1481. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, Z.; Yu, D.; Wang, X.; Pan, Y.; Zhang, G.; Shi, X. Influence of the Selection of Interpolation Method on Revealing Soil Organic Carbon Variability in the Red Soil Region, China. Sustainability 2018, 10, 2290. [Google Scholar] [CrossRef] [Green Version]
  19. Oliver, M.A.; Webster, R. Basic Steps in Geostatistics: The Variogram and Kriging; SpringerBriefs in Agriculture; Springer International Publishing: Cham, Switzerland, 2015. [Google Scholar]
  20. Oliver, M.A.; Bishop, T.F.A.; Marchant, B.P. Precision Agriculture for Sustainability and Environmental Protection; Routledge: Abingdon, UK, 2013. [Google Scholar]
  21. Benslama, A.; Khanchoul, K.; Benbrahim, F.; Boubehziz, S.; Chikhi, F.; Navarro-Pedreño, J. Monitoring the Variations of Soil Salinity in a Palm Grove in Southern Algeria. Sustainability 2020, 12, 6117. [Google Scholar] [CrossRef]
  22. Dai, F.; Zhou, Q.; Lv, Z.; Wang, X.; Liu, G. Spatial prediction of soil organic matter content integrating artificial neural network and ordinary kriging in Tibetan Plateau. Ecol. Indic. 2014, 45, 184–194. [Google Scholar] [CrossRef]
  23. Chabala, L.M.; Mulolwa, A.; Lungu, O. Application of ordinary kriging in mapping soil organic carbon in Zambia. Pedosphere 2017, 27, 338–343. [Google Scholar] [CrossRef]
  24. Song, J.J.; Kwon, S.; Lee, G. Incorporation of parameter uncertainty into spatial interpolation using Bayesian trans-Gaussian kriging. Adv. Atmos. Sci. 2015, 32, 413–423. [Google Scholar] [CrossRef]
  25. Radočaj, D.; Jurišić, M.; Gašparović, M.; Plaščak, I. Optimal Soybean (Glycine max L.) Land Suitability Using GIS-Based Multicriteria Analysis and Sentinel-2 Multitemporal Images. Remote Sens. 2020, 12, 1463. [Google Scholar] [CrossRef]
  26. Jurišić, M.; Plaščak, I.; Antonić, O.; Radočaj, D. Suitability Calculation for Red Spicy Pepper Cultivation (Capsicum annum L.) Using Hybrid GIS-Based Multicriteria Analysis. Agronomy 2020, 10, 3. [Google Scholar] [CrossRef] [Green Version]
  27. Radočaj, D.; Jurišić, M.; Zebec, V.; Plaščak, I. Delineation of Soil Texture Suitability Zones for Soybean Cultivation: A Case Study in Continental Croatia. Agronomy 2020, 10, 823. [Google Scholar] [CrossRef]
  28. Zhao, J.; Zhan, X.; Jiang, Y.; Xu, J. Variations in climatic suitability and planting regionalization for potato in northern China under climate change. PLoS ONE 2018, 13, e0203538. [Google Scholar] [CrossRef] [PubMed]
  29. Bogunovic, I.; Trevisani, S.; Seput, M.; Juzbasic, D.; Durdevic, B. Short-range and regional spatial variability of soil chemical properties in an agro-ecosystem in eastern Croatia. Catena 2017, 154, 50–62. [Google Scholar] [CrossRef]
  30. Croatian Bureau of Statistics, Statistical Yearbook of the Republic of Croatia. 2018. Available online: https://www.dzs.hr/Hrv_Eng/ljetopis/2018/sljh2018.pdf (accessed on 28 November 2020).
  31. Croatian Meteorological and Hydrological Service, Monthly Values and Extremes, Osijek. Available online: https://meteo.hr/klima_e.php?section=klima_podaci&param=k1&Grad=osijek (accessed on 8 December 2020).
  32. Liu, N.; Li, Y.; Cong, P.; Wang, J.; Guo, W.; Pang, H.; Zhang, L. Depth of straw incorporation significantly alters crop yield, soil organic carbon and total nitrogen in the North China Plain. Soil Tillage Res. 2021, 205, 104772. [Google Scholar] [CrossRef]
  33. Metwally, M.S.; Shaddad, S.M.; Liu, M.; Yao, R.-J.; Abdo, A.I.; Li, P.; Jiao, J.; Chen, X. Soil Properties Spatial Variability and Delineation of Site-Specific Management Zones Based on Soil Fertility Using Fuzzy Clustering in a Hilly Field in Jianyang, Sichuan, China. Sustainability 2019, 11, 7084. [Google Scholar] [CrossRef] [Green Version]
  34. Gan, Y.T.; Liang, B.C.; Liu, L.P.; Wang, X.Y.; McDonald, C.L. C: N ratios and carbon distribution profile across rooting zones in oilseed and pulse crops. Crop. Pasture Sci. 2011, 62, 496–503. [Google Scholar] [CrossRef]
  35. Rosner, K.; Bodner, G.; Hage-Ahmed, K.; Steinkellner, S. Long-term Soil Tillage and Cover Cropping Affected Arbuscular Mycorrhizal Fungi, Nutrient Concentrations, and Yield in Sunflower. Agron. J. 2018, 110, 2664–2672. [Google Scholar] [CrossRef]
  36. Liu, Z.; Chen, X.; Jing, Y.; Li, Q.; Zhang, J.; Huang, Q. Effects of biochar amendment on rapeseed and sweet potato yields and water stable aggregate in upland red soil. Catena 2014, 123, 45–51. [Google Scholar] [CrossRef]
  37. Ruffo, M.L.; Bullock, D.G.; Bollero, G.A. Soybean yield as affected by biomass and nitrogen uptake of cereal rye in winter cover crop rotations. Agron. J. 2004, 96, 800–805. [Google Scholar] [CrossRef]
  38. Mendes, M.P.; Ribeiro, L. Nitrate probability mapping in the northern aquifer alluvial system of the river Tagus (Portugal) using Disjunctive Kriging. Sci. Total Environ. 2010, 408, 1021–1034. [Google Scholar] [CrossRef] [PubMed]
  39. Isaaks, E.H.; Srivastava, R.M. An Introduction to Applied Geostatistics; Oxford University Press: New York, NY, USA, 1989. [Google Scholar]
  40. Panday, D.; Maharjan, B.; Chalise, D.; Shrestha, R.K.; Twanabasu, B. Digital soil mapping in the Bara district of Nepal using kriging tool in ArcGIS. PLoS ONE 2018, 13, e0206350. [Google Scholar] [CrossRef]
  41. Zhang, Z.; Liu, M.; Sun, J.; Wei, T. Degradation leads to dramatic decrease in topsoil but not subsoil root biomass in an alpine meadow on the Tibetan Plateau, China. J. Arid Land 2020, 1–13. [Google Scholar] [CrossRef]
  42. United States Department of Agriculture, Soil Survey Laboratory Information Manual. Available online: https://www.nrcs.usda.gov/Internet/FSE_DOCUMENTS/nrcs142p2_052226.pdf (accessed on 8 December 2020).
  43. Thomsen, I.K.; Petersen, B.M.; Bruun, S.; Jensen, L.S.; Christensen, B.T. Estimating soil C loss potentials from the C to N ratio. Soil Biol. Biochem. 2008, 40, 849–852. [Google Scholar] [CrossRef]
  44. Marchetti, A.; Piccini, C.; Francaviglia, R.; Mabit, L. Spatial distribution of soil organic matter using geostatistics: A key indicator to assess soil degradation status in central Italy. Pedosphere 2012, 22, 230–242. [Google Scholar] [CrossRef]
  45. Sun, X.L.; Yang, Q.; Wang, H.L.; Wu, Y.J. Can regression determination, nugget-to-sill ratio and sampling spacing determine relative performance of regression kriging over ordinary kriging? Catena 2019, 181, 104092. [Google Scholar] [CrossRef]
  46. Varouchakis, E.A.; Hristopulos, D.T.; Karatzas, G.P. Improving kriging of groundwater level data using nonlinear normalizing transformations—A field application. Hydrol. Sci. J. 2012, 57, 1404–1419. [Google Scholar] [CrossRef] [Green Version]
  47. Schloeder, C.A.; Zimmerman, N.E.; Jacobs, M.J. Comparison of methods for interpolating soil properties using limited data. Soil Sci. Soc. Am. J. 2001, 65, 470–479. [Google Scholar] [CrossRef]
  48. Ding, J.; Yu, D. Monitoring and evaluating spatial variability of soil salinity in dry and wet seasons in the Werigan–Kuqa Oasis, China, using remote sensing and electromagnetic induction instruments. Geoderma 2014, 235, 316–322. [Google Scholar] [CrossRef]
  49. Fashi, F.H.; Sharifi, F.; Kamali, K. Modelling infiltration and geostatistical analysis of spatial variability of sorptivity and transmissivity in a flood spreading area. Span. J. Agric. Res. 2014, 12, 277–288. [Google Scholar] [CrossRef] [Green Version]
  50. Jenkinson, D.S.; Poulton, P.R.; Bryant, C. The turnover of organic carbon in subsoils. Part 1. Natural and bomb radiocarbon in soil profiles from the Rothamsted long-term field experiments. Eur. J. Soil Sci. 2008, 59, 391–399. [Google Scholar] [CrossRef]
  51. Rumpel, C.; Kögel-Knabner, I. Deep soil organic matter—A key but poorly understood component of terrestrial C cycle. Plant Soil 2011, 338, 143–158. [Google Scholar] [CrossRef]
  52. Tesfahunegn, G.B.; Tamene, L.; Vlek, P.L. Catchment-scale spatial variability of soil properties and implications on site-specific soil management in northern Ethiopia. Soil Tillage Res. 2011, 117, 124–139. [Google Scholar] [CrossRef]
  53. Yang, L.; Wang, L.; Huang, J.; Mansaray, L.R.; Mijiti, R. Monitoring policy-driven crop area adjustments in northeast China using Landsat-8 imagery. Int. J. Appl. Earth Obs. Geoinf. 2019, 82, 101892. [Google Scholar] [CrossRef]
  54. Singha, C.; Swain, K.C.; Swain, S.K. Best Crop Rotation Selection with GIS-AHP Technique Using Soil Nutrient Variability. Agriculture 2020, 10, 213. [Google Scholar] [CrossRef]
  55. Radočaj, D.; Obhođaš, J.; Jurišić, M.; Gašparović, M. Global Open Data Remote Sensing Satellite Missions for Land Monitoring and Conservation: A Review. Land 2020, 9, 402. [Google Scholar] [CrossRef]
  56. Gašparović, M.; Zrinjski, M.; Barković, Đ.; Radočaj, D. An automatic method for weed mapping in oat fields based on UAV imagery. Comput. Electron. Agric. 2020, 173, 105385. [Google Scholar] [CrossRef]
Figure 1. The location of Osijek-Baranja County and soil samples at 0–10 and 20–30 cm soil layers.
Figure 1. The location of Osijek-Baranja County and soil samples at 0–10 and 20–30 cm soil layers.
Agronomy 10 01996 g001
Figure 2. Annual area covered by five major crops in the Osijek-Baranja County during 2017–2019 and the three enlarged example subsets.
Figure 2. Annual area covered by five major crops in the Osijek-Baranja County during 2017–2019 and the three enlarged example subsets.
Agronomy 10 01996 g002
Figure 3. Thiessen polygons for the stationarity evaluation of input samples.
Figure 3. Thiessen polygons for the stationarity evaluation of input samples.
Agronomy 10 01996 g003
Figure 4. Variograms of mathematical models producing highest inner and outer accuracy.
Figure 4. Variograms of mathematical models producing highest inner and outer accuracy.
Agronomy 10 01996 g004
Figure 5. Soil carbon-to-nitrogen ratio (C/N) maps calculated using various interpolation parameters.
Figure 5. Soil carbon-to-nitrogen ratio (C/N) maps calculated using various interpolation parameters.
Agronomy 10 01996 g005
Figure 6. Combined scatterplot of interpolation results in 0–10 cm soil depth.
Figure 6. Combined scatterplot of interpolation results in 0–10 cm soil depth.
Agronomy 10 01996 g006
Figure 7. Combined scatterplot of interpolation results in 20–30 cm soil depth.
Figure 7. Combined scatterplot of interpolation results in 20–30 cm soil depth.
Agronomy 10 01996 g007
Figure 8. Mean and standard deviation values of soil C/N for five major crops for 2017 and projected values for 2018 and 2019.
Figure 8. Mean and standard deviation values of soil C/N for five major crops for 2017 and projected values for 2018 and 2019.
Agronomy 10 01996 g008
Table 1. Annual crop coverage statistics of five major crops in the Osijek-Baranja County.
Table 1. Annual crop coverage statistics of five major crops in the Osijek-Baranja County.
Crop201720182019
Area Cultivated (ha)Part of Total Agricultural Area (%)Area Cultivated (ha)Part of Total Agricultural Area (%)Area Cultivated (ha)Part of Total Agricultural Area (%)
Maize42,36820.1940,23219.1956,54927.17
Wheat32,94115.7040,48619.3142,04820.20
Sunflower22,34710.6523,81811.3620,1799.69
Rapeseed20,98110.0021,03910.0314,1376.79
Soybean16,6427.9315,1177.2113,1996.34
Total135,27964.46140,69267.09146,11270.19
Table 2. Descriptive statistics and Shapiro–Wilk normality test for input samples.
Table 2. Descriptive statistics and Shapiro–Wilk normality test for input samples.
ValuesC 0–10 cmN 0–10 cmC 20–30 cmN 20–30 cm
Mean (g 100 g−1)2.580.241.980.16
Minimum (g 100 g−1)0.970.090.740.04
Maximum (g 100 g−1)5.880.665.420.48
CV0.480.520.590.53
Skewness0.951.501.341.87
Kurtosis0.531.821.263.61
Shapiro–WilkW0.810.470.850.42
p0.029>0.001>0.001>0.001
Table 3. Statistical properties of interpolation results for 0–10 cm soil depth.
Table 3. Statistical properties of interpolation results for 0–10 cm soil depth.
Soil ElementModelTransformation MethodnsSpatial Dependencer (m) R I 2   ( % ) R O 2   ( % )
Clinearnone0.3592.3670.8482534895.753.5
logarithmic0.0640.3450.81493.658.8
square rootnone0.2982.1590.8623747196.042.8
logarithmic0.0250.3420.92796.855.1
Gaussiannone0.1362.2700.9402534897.158.7
logarithmic0.0410.3300.87695.059.5
sphericalnone0.0981.7980.9452476899.254.3
logarithmic0.0220.2800.92198.559.1
Nlinearnone0.0370.1880.8034298195.842.9
logarithmic0.1470.6720.78180.565.5
square rootnone0.0350.1350.7415508997.438.0
logarithmic0.1580.5140.69381.163.0
Gaussiannone0.0090.1900.9534298193.849.5
logarithmic0.0710.6490.89188.166.3
sphericalnone0.0030.1000.9703604999.750.9
logarithmic0.0870.4170.79192.263.3
Table 4. Statistical properties of interpolation results for 20–30 cm soil depth.
Table 4. Statistical properties of interpolation results for 20–30 cm soil depth.
Soil ElementModelTransformation MethodnsSpatial Dependencer (m) R I 2   ( % ) R O 2   ( % )
Clinearnone0.7051.5360.54121,19472.057.5
logarithmic0.2200.2970.25954.556.2
square rootnone0.3191.4580.78120,01789.057.0
logarithmic0.1960.2860.31564.558.7
Gaussiannone0.8061.4610.44823,54968.544.4
logarithmic0.2430.2850.14755.567.3
sphericalnone0.6231.4710.57622,39081.443.1
logarithmic0.2180.3200.31968.564.7
Nlinearnone0.0040.0860.95325,90495.640.7
logarithmic0.0140.4890.97193.074.6
square rootnone0.0020.0770.97441,21195.538.6
logarithmic0.0030.4450.99389.079.3
Gaussiannone0.0010.0830.98831,79168.961.5
logarithmic0.0010.4440.99899.380.1
sphericalnone0.0010.0590.98328,84997.947.8
logarithmic0.0060.3560.98196.278.7
Table 5. Descriptive statistics of the final soil C/N interpolation results.
Table 5. Descriptive statistics of the final soil C/N interpolation results.
Soil DepthModelTransformation MethodMeanMinMaxCV
0–10 cmlinearnone10.804.6920.260.31
logarithmic11.283.7620.850.26
square rootnone10.814.9319.280.25
logarithmic11.285.8219.860.20
Gaussiannone11.453.9231.700.48
logarithmic12.773.8837.350.46
sphericalnone11.254.1623.910.32
logarithmic11.174.2421.770.22
20–30 cmlinearnone12.054.4722.200.24
logarithmic11.704.4426.760.25
square rootnone12.045.5322.570.14
logarithmic11.544.6825.740.21
Gaussiannone12.883.9731.940.39
logarithmic11.863.1531.020.36
sphericalnone12.104.1726.150.23
logarithmic11.623.6528.600.25
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jurišić, M.; Radočaj, D.; Krčmar, S.; Plaščak, I.; Gašparović, M. Geostatistical Analysis of Soil C/N Deficiency and Its Effect on Agricultural Land Management of Major Crops in Eastern Croatia. Agronomy 2020, 10, 1996. https://doi.org/10.3390/agronomy10121996

AMA Style

Jurišić M, Radočaj D, Krčmar S, Plaščak I, Gašparović M. Geostatistical Analysis of Soil C/N Deficiency and Its Effect on Agricultural Land Management of Major Crops in Eastern Croatia. Agronomy. 2020; 10(12):1996. https://doi.org/10.3390/agronomy10121996

Chicago/Turabian Style

Jurišić, Mladen, Dorijan Radočaj, Stjepan Krčmar, Ivan Plaščak, and Mateo Gašparović. 2020. "Geostatistical Analysis of Soil C/N Deficiency and Its Effect on Agricultural Land Management of Major Crops in Eastern Croatia" Agronomy 10, no. 12: 1996. https://doi.org/10.3390/agronomy10121996

APA Style

Jurišić, M., Radočaj, D., Krčmar, S., Plaščak, I., & Gašparović, M. (2020). Geostatistical Analysis of Soil C/N Deficiency and Its Effect on Agricultural Land Management of Major Crops in Eastern Croatia. Agronomy, 10(12), 1996. https://doi.org/10.3390/agronomy10121996

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop