Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Evaluation of GPM-IMERG Precipitation Product at Multiple Spatial and Sub-Daily Temporal Scales over Mainland China
Next Article in Special Issue
UAV-Hyperspectral Imaging to Estimate Species Distribution in Salt Marshes: A Case Study in the Cadiz Bay (SW Spain)
Previous Article in Journal
Mapping the Distribution and Dynamics of Coniferous Forests in Large Areas from 1985 to 2020 Combining Deep Learning and Google Earth Engine
Previous Article in Special Issue
Potential of Satellite Spectral Resolution Vegetation Indices for Estimation of Canopy Chlorophyll Content of Field Crops: Mitigating Effects of Leaf Angle Distribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cropland Productivity Evaluation: A 100 m Resolution Country Assessment Combining Earth Observation and Direct Measurements

1
Centre for Agricultural Research, Institute for Soil Sciences, Hermann Otto 15., H-1022 Budapest, Hungary
2
Institute of Advanced Studies, Chernel u. 14., H-9730 Kőszeg, Hungary
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(5), 1236; https://doi.org/10.3390/rs15051236
Submission received: 6 December 2022 / Revised: 21 February 2023 / Accepted: 22 February 2023 / Published: 23 February 2023
(This article belongs to the Special Issue Crops and Vegetation Monitoring with Remote/Proximal Sensing)

Abstract

:
A methodology is presented for the quantitative assessment of soil biomass productivity at 100 m spatial resolution on a national scale. The traditional land evaluation approach—where crop yield is the dependent variable—was followed using measured yield and net primary productivity data derived from satellite images, together with digital soil and climate maps. In addition to characterizing of soil biomass productivity based on measured data, the weight of soil properties on productivity was also quantified to provide measured soil health and soil quality indicators as an information base for designing sustainable land management practices. To produce these results, we used only the Random Forest method for our calculations. The study considers high-input agriculture, which is predominant in the country. Biomass productivity indices for the main crops (wheat, maize and sunflowers) and general productivity indices were calculated for the whole agricultural area of Hungary. Results can be implemented in cadastral systems, in applied in agricultural and rural development programs. The assessment can be repeated for monitoring purposes to support general monitoring objectives as well as for reporting in relation to the United Nations Sustainable Development Goals. However, on the basis of the results, we also propose a method for periodically updating the assessment, which can also be used for monitoring biomass productivity in the context of climate change, land degradation and the development of cultivation technology.

1. Introduction

A key natural resource that ensures food security, ecological security and sustainable development is cultivable land. Recently, the importance of soil has been increasingly put into focus as the general public also become more aware of it as a non-renewable resource that can be lost quickly if improperly used or managed with very little chance of regeneration. Despite the critical importance of soil productivity, not only as indicator, but also in sustaining life on Earth, knowledge of the spatial and temporal variability of soil from regional to global scales is limited or fragmented. For the creation of effective agricultural and food policies at the regional levels, accurate soil productivity predictions are essential. The limited information on soil productivity hinders national (Farmers’ Soil Conservation Programme, National Rural Development Programme) and international (EU Soil Mission) programs to monitoring its changes and build future scenarios on it.
The Sustainable Development Goals (SDGs) of the United Nations’ Agenda 2030 framework include targets that recommend direct consideration of land and soil resources [1,2,3], which were adopted by all United Nations member states in 2015. Soil resources are linked to the SDGs through several soil functions [2], of which the biomass productivity function is at the core of SDGs 2.3 and 2.4., which explicitly target the sustainable increases in agricultural productivity. Furthermore, biomass productivity is proposed as an indicator of land degradation [4], which is linked to SGD 15.3 [5].
Biomass productivity is conditioned by inherent soil properties, climatic and management factors, thus variable in both space and time [6]. Spatial variability of soil productivity is traditionally assessed within the broad framework of land evaluation [7]. However, land evaluation should also include socio-economic components [8], which are not necessary for soil productivity evaluation. Nevertheless, soil is an integral part of the land with a distinct spatial location and therefore biophysical characteristics of the studied sites, such as climate and relief conditions, need to be taken into account when assessing its productivity [9].
The aim of classical quantitative land evaluation is to establish productivity indices based on actual yields in order to reflect production potentials for taxation and planning purposes [10,11,12,13,14,15,16,17,18]. A similar quantitative approach can be applied to reveal soil biomass productivity, its drivers and changes for monitoring purposes.
Dynamic and simulation models [7,16,19,20,21,22] can provide an alternative to classical productivity evaluation, but their validation still requires measured biomass or yield data. Advantages of the classical data-driven assessment, i.e., where yield is the dependent variable and biophysical factors are independent inputs, are high reliability, explicit spatial validity and easy interpretation. Process-based modeling and statistical modeling are also two frequently employed techniques for forecasting crop yield responses to climate variability. Process-based crop models are effective for predicting crop yields because they simulate physiological processes of crop growth and development in response to environmental factors and management techniques, especially at the field scale [23]. Traditional regression techniques have some drawbacks that can be addressed by statistical modeling techniques based on machine-learning algorithms. Machine-learning techniques have been used increasingly in recent years as niche-based classification modeling tools [24,25,26]. For our analysis we selected the Random Forest (RF) technique [27,28], which uses the Classification and Regression Trees method as the basis for growing multiple classification trees. The study considers high-input agriculture, which is predominant in the country and uses time series information (measured crop yield statistics and satellite-derived biomass productivity indicators).
A scientific-based biomass productivity assessment should be based on a numerical assessment of production potential based on statistical studies. Previous national land evaluation techniques were estimation procedures, which inevitably introduced classification errors. Since the only objective measure of land quality is yield over time, our method is designed with yield as the dependent variable and environmental factors (soil, climate, topography) that affect yield as the independent variables. The method must be designed in such a way that the parameterization process can be repeated as the amount of available data increases, so that the land classification system can be easily revised and refined at any time on the basis of changes in production conditions.
Based on the above considerations, we performed a detailed study with country coverage with the following aims: (i) to identify main soil and climatic determinants of biomass productivity, (ii) to quantify the weights of soil and climatic factors of productivity for the main crop types (wheat, maize, sunflowers), (iii) to produce crop-specific and general productivity maps for all agricultural land of the country, and (iv) to propose a methodology for integrated monitoring of biomass productivity.

2. Materials and Methods

Soil biomass productivity evaluation must be based on biomass data and the assessment of the environmental and management factors influencing it. This requires biomass data, and geographical and management data, including soil, topography, climate and fertilizer data. Country-wide implementation of the agricultural biomass productivity model can only be based on information that is available for the full agricultural area of the country. To ensure the best possible spatial detail to develop and implement a new productivity model, data of dependent variables (measured yield and remotely sensed biomass indicators) as well as independent variables on soil properties were collected at parcel scale and implemented at soil property maps of 100 m resolution.

2.1. Study Area

Hungary is located in Central Europe and the Carpathian Basin, which is a part of the Pannonian biogeographic region (45°43′ to 48°35′N and 16°06′ to 22°53′E). The country is 93,033 km2 and has an elevation range between 77 and 1014 m above sea level, and agricultural lands are typically located between 77 and 350 m altitude. Agriculture is the dominant land use, with non-irrigated arable land (Figure 1) accounting for 61% of the country’s total area [29]. Winter wheat (Triticum aestivum), maize (Zea mays) and sunflowers (Helianthus annuus), which have been sown on up to 80% of Hungary’s arable land in recent decades, were selected for the productivity assessment.

2.2. Databases

2.2.1. National Plot-Level Field Soil, Fertilization and Yield Databases (AIIR Field Database)

The AIIR database [31] contains crop type, yield, fertilization and soil information for each cultivated parcel, summing up to 80,000 cultivated parcels of Hungary for 5 years (1985–1989). The data were provided by the Central Plant and Soil Conservation Service (Budapest) for the purpose of land evaluation research. The sampling for the soil tests was carried out in such a way that the parcels were divided into 12 ha sections and then, along the diagonals of the selected sections, soil samples were taken from at least 20 locations using the so-called parallel sampling method. The subsamples were taken homogenized, so that an average sample was taken from the subplots of each agricultural field. For areas with a slope greater than 12%, average samples were taken separately for each (upper, middle, lower) section of the slope, taking into account erosion and different soil nutrient supply. The database was digitized in 2000 and in 2014 was upgraded to a modern geo-spatial database (point data with coordinates). We have selected the points that still fall on arable land at the time of our study. The database includes the following three major types of data:
-
Basic data of the parcels (location, size, land user);
-
Soil taxonomical and laboratory analysis data (soil type and subtype, pH, texture, organic matter, nitrogen, phosphorus and potassium content);
-
Agricultural management data (crop type, yield, date of sowing, fertilization and harvest, fertilizer doses);
-
Crop type and yield data.
Distribution of data by soil types is presented in Table 1.

2.2.2. Remote Sensing Derived Biomass Productivity Indicators

Long term (2003–2018) time series remote sensing data were used to derive mean gross primary productivity (GPP) values as proposed by Jin and Eklundh (2014) [34]. The MODIS dataset (MOD17) [35] was used at a nominal 500 m spatial resolution to produce GPP datasets for the whole country. It is important to note that crop yields and GPP represent different aspects of productivity. However, in managed cropland there is a strong correlation between the two [36]. We used the normalized productivity (value range 1–100) as the target variable, and all of our results were normalized between 1 and 100, making it easier to integrate into our model.

2.2.3. Time Series Meteorological Data

The Central-European FORESEE meteorological database [37], which covers the whole area of the country with a 0.1 × 0.1 degree grid, was used to derive mean temperature and total precipitation at monthly scales (between 1951 and 2013). Mean temperature and precipitation values were linked to the spatial units (100 m pixels) of the assessment. The downscaling was performed by the bilinear resampling method.

2.2.4. Topographic Data

The Shuttle Radar Topography Mission [38,39] provides a dataset of 30 m resolution grid cells as the basis for the digital elevation model (DEM). SRTM mapped Earth’s topography between 56 degrees south and 60 degrees north of the equator. SRTM has a vertical accuracy of 5.3 m (RMSE) in Hungary [40]. The SRTM-derived DEM was used to include a topographic component to the land evaluation model.

2.2.5. Land Use Data

The Hungarian coverage of the CORINE [30] land cover database for the year 2018 was used to delineate croplands in the country. The 1:100,000 scale datasets have a minimum mapping unit of 25 ha for patches and a minimum width of 100 m for linear elements. A total of 44 land cover and land use categories are included in the dataset, 28 of which are appropriate for Hungary [30]. All assessments and the map visualization of the results were based on the cropland areas (see Figure 1).

2.2.6. Map Series of Soil Types and Soil Properties

The unified national soil type and soil property maps of Pásztor et al. (2020, 2018, 2017, 2015) [41,42,43,44] provided the soil information base for the assessment. A total of 41 soil types, belonging to 9 main soil type groups of the Hungarian Soil Classification System [45], are covered by the dataset. Soil chemical and physical data include pH, calcium carbonate content, organic matter content and texture. The map series are all produced at a 100 m resolution and can be viewed on the dosoremi.hu website. The 100 m resolution of the soil maps was considered to be sufficiently detailed for parcel-scale productivity evaluation, and therefore this spatial resolution defined the resolution of the assessment. There is a slight difference in the semantic component of the soil type maps and the soil type information in the AIIR dataset (Table 1). There are soil types in the national soil map with areas covering <1% of the country that are not available in the AIIR dataset, or which are available only with a very limited sample size. These were not sufficient for statistical tests. This minor inconsistency required an expert-based modification of the final evaluation system.

2.3. Data Preparation

A quality and consistency check of the AIIR dataset was carried out in the first phase of the data preparation to filter out typos and false records. Inconsistent records (outliers), such as soil samples with high carbonate content and low pH, were excluded from the dataset. We then selected those records from the AIIR dataset that corresponded to agricultural parcels of intensive (i.e., high fertilizer use) cultivation. The selection was made based on the amount of fertilizers applied, and records containing at least 125 kg × ha−1 of nitrogen and 30 kg × ha−1 of active phosphorus input were kept. In this way, the analysis of the current assessment focused on data from intensively cultivated fields.
Winter wheat (Triticum aestivum), maize (Zea mays) and sunflowers (Helianthus annuus), which have been sown on up to 80% of the croplands in Hungary [46] in recent decades, were selected for the productivity evaluation. In order to establish a common basis for the analysis, the yield data of these three main crops from each parcel of the dataset were normalized to a scale of 1 to 100. For the same reason, the GPP values were also normalized to a scale of 1 to 100. Normalization was applied to all wheat, maize and sunflower yield data in the five years covered by the AIIR database and to all cropland pixels in the GPP dataset.
The AIIR database with normalized yield data and the normalized GPP dataset were integrated with the climate geodatabase into a single geodatabase using geographical coordinates as unique identifiers. The result was a georeferenced dataset created to include all soil, climate, management and yield data. Productivity analysis was carried out using information of the georeferenced pixels, including their geographical coordinates.
The GPP data, originally produced at 500 m resolution, were downscaled to 100 m resolution and normalized to values between 1 and 100. The downscaling was performed by the nearest neighbor resampling method The SRTM data, which were originally produced at 30 m resolution were generalized to 100 m resolution using the bilinear interpolation technique.
All datasets were converted to the Uniform National Projection System (EOV) to create a coherent geodatabase.

2.4. Assessment and Implementation Methods

2.4.1. Model Development

Soil biomass productivity assessment is the process of establishing relationships between soil properties and yields. Data-mining methods are tools for revealing hidden relationships in datasets structured by input variables. In soil assessment, data mining can help to identify the most important factors in yield formation and establish the weights of these factors. For our analysis we chose the Random Forest technique [27,28], which uses the Classification and Regression Trees method as a basis for growing multiple classification trees. For this operation, the database is divided into a series of training and test datasets to establish and validate relationships, respectively. Each training dataset (80% of the dataset) is a randomly selected subset that is used to develop a tree model using randomly selected predictors. The remaining data (10%) after the random selection of the subset (test data, 10% of dataset) are used to validate the developed model [47]. We used the createDataPartition function from the caret package to select data randomly. The generalized error of the forest depends on two parameters: how accurate each individual classifier is and how independent the different classifiers are from each other (i.e., the strength of each tree in the forest and the correlation between them). The Random Forest analysis was performed with the ranger R package [48]. The long-term means normalized productivity index (MNPI), taking into account both the measured AIIR data and the GPP data, was computed by taking the average of the two normalized datasets. The Random Forest operation was performed with the MNPI as the dependent variable and the environmental (soil, climate) variables as explanatory variables (Figure 2). First, the assessment was carried out separately for winter wheat, maize and sunflowers in order to evaluate crop-specific productivity of Hungarian croplands using the MNPI data of these crops. As a result, crop-specific productivity indices were produced for the three main crops. As our overall interest was to establish the MNPI for each Hungarian parcel at 100 m resolution, three parallel models were developed for the three major crops (wheat, maize, sunflowers) based on the crop-specific entries of the normalized yield data, and a fourth, a general productivity model, was developed based on the MNPI. As a result, both crop-specific (weighted means, wheat 40%, maize 40% and sunflowers 20%) and general productivity indices were assigned to climate and soil property combinations. Due to the limited information for some minor soil types (i.e., occupying area < 0.5% of agricultural lands), statistical testing could not be successfully performed for these soils. To assess the productivity evaluation of these soils, two evaluation approaches were applied and their results were combined. Firstly, an expert-based judgement was carried out. Productivity indices were established considering those of closely related soils in the Hungarian soil taxonomy using information from previous land evaluation systems [49], related literature [50,51,52,53,54] and expert knowledge. Secondly, a statistical test based solely on the GPP data was carried out to evaluate the effect of soil properties and climate, although without statistically significant results, but for orientation purposes. The relative importance of the explanatory variables was calculated. We analyzed the importance of all variables using the imp function of bclust package in R [55]. Relative importance was calculated by dividing the importance score of each variable by the largest importance score of the variable, and then multiplying by 100. Harmonizing the results of the two approaches ensured the consistency across the system, even for parcels with soils that make up a small proportion of the country’s croplands. The theoretical range of the final productivity indices was set between 1 and 100, corresponding to the normalized yield values of the test dataset and following the indexing approach of traditional soil productivity evaluation of Hungary [51]. Model validation was performed using normalized yield data as independent variables of the test dataset. The test dataset included a randomly selected 10% of the data and a predict function of the ranger package was used. We calculated the correlation coefficient to show the relationship between the observed and the predicted values, the mean absolute error (MAE) to show the distance of the predicted values from the observed values [56], and the mean absolute percentage error (MAPE) to show the percentage of error between observed and predicted values [57].

2.4.2. Spatial Implementation

Soil, meteorological and digital terrain maps were used for the spatial implementation of the soil biomass productivity model, i.e., to produce soil productivity maps. The developed model provides productivity indices on a scale of 1 to 100 for several combinations of climate and soil properties in the country. Basic input layers for the spatial implementation include detailed (100 m) soil type and soil property maps and climate data. Slope correction coefficients (see Appendix A Table A1) from the previous official Hungarian land evaluation model [49] were applied to produce the final productivity indices. The coefficients reflect the effect of slope angle and slope direction on productivity. The SRTM digital topographic data were used to implement the correction coefficients and to produce the final maps. Presentation of the results covers all cropland areas of the country at a 100 m resolution.

3. Results

3.1. Model Development and Estimation Efficiency

The general, country-wide productivity model using biophysical explanatory variables explains up to 40% of the biomass productivity in the country (R2 = 0.402). This model fit can be considered adequate for a country scale assessment, especially for a country with a wide variety of soil types from salt-affected soils to Arenosols, Luvisols and chernozems. The efficiency of the crop-specific models is best for wheat, followed by maize and sunflowers, in the order of the available sample size, respectively (Figure 3 and Table 2). Results were statistically significant at the 0.01 level.
The combination of measured and satellite-driven data for the general productivity model development gave almost the same model fit as the crop-specific one for wheat, which was based on a large sample size of measured yields. The descriptive power of sunflower productivity estimation was not as strong (Figure 3D). The MAPE results are as follows: all cropland 19.28%, wheat 18.07%, maize 19.17% and sunflowers 29.81%. The most accurate prediction based on the MAPE and MAE results was for wheat followed by the maize and sunflower predictions.

3.2. Baseline Biomass Productivity Indices and Map for Croplands of Hungary

By implementing the biomass productivity model on the national soil, climate and topographic geodatabase, a new soil biomass productivity map was produced (Figure 4). The map shows the general productivity of croplands. In the same process, crop-specific productivity maps were also produced. While the crop-specific productivity indices and maps can be used for planning land use and cropping, the general productivity map provides an overview of the spatial pattern of biomass potential of agricultural parcels in the country.
The map confirms the empirical knowledge that the most fertile areas are on chernozem soils in the east and on various loamy soils in the west of the country. Sandy soils, whether in the western, the central or the eastern part of the country, perform rather poorly. This phenomenon is typical of a country where water supply is the main climatic factor limiting crop production.
The mean productivity index for all the croplands is 64.7, with a standard deviation of 13.4, reflecting the dominancy of medium-to-good land within the agricultural areas of the country in terms of the spatial extent (Figure 4A).
The mean productivity index after slope coefficient correction for all the croplands of the country is 58.9, with a standard deviation of 18.5, reflecting the dominancy of medium-to-good land within the agricultural areas of the country (Figure 5).

3.3. Soil and Climatic Determinants of Biomass Productivity in Hungary

In general, humus content (which also reflects organic matter content), pH, CaCO3 content, soil type and soil texture were the most important soil-based input parameters for predicting wheat productivity (Figure 6A). However, the geographical location was found to be an even more important explanatory variable. While this information suggests the importance of climate, the measured climatic variables ranked lower in the importance list. The mean temperatures of January, February, December and June (in the order of importance) are the most important thermal parameters for productivity. Regarding precipitation, the amounts in November, August, June and December are the most important determinants.
In the case of maize, on the other hand, the measured climatic variables were found to be of high importance, together with humus, pH, CaCO3 and soil type, while the location was not considered to be relevant. These differences suggest the appropriateness of the crop-specific evaluation approach. There are also differences in the mean temperature and precipitation. Figure 6B shows that the precipitation (July, November, May, August and October) is a more important factor than temperature (most important in January and February) in the case of maize.
Based on Figure 6C, the variables of the sunflower prediction show a completely different pattern. The most important variables are the location and mean temperature in February, followed by humus content, mean temperature in January, March and November, pH, precipitation in November and CaCO3. Soil type and texture are the least important variables. Climatic variables (mainly amount of precipitation) have a more significant effect on the sunflower yield amount than the soil type and texture.

4. Discussion

Biomass productivity is a dynamic property that changes over time, partly due to changing climatic conditions (within and between years) and changing soil properties (pH, organic matter, soil nutrients content, etc.), but also due to new crop varieties and advances in crop management having an important influence. The effect of the changes in the biophysical factors may be synergistic or in the opposite direction. Nevertheless, it is possible to estimate the weight of the factors in productivity on a reasonable time scale. Twenty to thirty years seem to be an adequate time scale for estimating soil biomass productivity and for identifying the weights of different factors in it. A moving timeframe with intervals of 3 to 6 years can be proposed for the updating of the biomass productivity indices. If the system is to be used for monitoring purposes, soil biomass productivity will need to be compared on the basis of different time periods, e.g., on moving time windows or trends and supplemented by the monitoring of soil properties, subject to degradation. The moving time window for biomass productivity monitoring can be harmonized with the periodic assessments in soil monitoring, i.e., 3–6 years.
Our validation results (predicted data vs. measured data in the test set) showed that there is a significant difference in the prediction accuracy between the different crops. The sunflower model has a lower performance in calculating biomass productivity, which may be due to the fewer number of data available for model training. Furthermore, sunflowers are a cash crop grown on very diverse soils and not so much linked to bioclimatic factors and soil parameters [58,59], while the R2 value (0.41) for the biomass productivity map can be regarded as adequate for a national estimate. The MAPE value indicates that only 18.06% of the model is inaccurate, that is, above the accuracy of results published in other studies [60,61]. The MAE indicator values indicate an average deviation of 6.78, which is not outstanding on a scale of 1 to 100. Cheng et al. (2022) [25] found similar R2 values in the case of maize and wheat based on MODIS GPP values, with stronger correlation in the case of maize. However, other studies in the case of wheat presented lower values [24,61,62]. Although the performance of the sunflower productivity model is rather low, its inclusion in the assessment provides a more comprehensive overview of plant-specific productivities and their differences, including major factors and the varying weights of the factors in plant-specific productivities. Furthermore, the inclusion of an additional plant-specific model extends the potential of the applied method to provide a general soil productivity assessment considering multiple crops, which is often needed in land use planning.
A new soil biomass productivity map was created by applying the biomass productivity model using national soil, climate and topography geodatabase. Crop-specific productivity maps, which shall be the ultimate source of multicriteria land use planning [63], were also produced using the same technique. The spatial distribution of biomass potential is shown on the general productivity map, which can be used to plan land use and crop production. Further to that, weights of individual soil and climate parameters of crop-specific productivity indices were also derived. For our final model results (soil biomass productivity) we applied a correction that takes into account the topography. Slope angle and orientation both matter for solar radiation to be taken into account. Our solution for incorporating terrain indices from an earlier national land evaluation model [51] was tailored for Hungarian conditions (average slope under 2.3%), instead of applying more complex methods [64,65] used in other pedoclimatic conditions. Random Forest is considered to be an appropriate method to predict crop-specific biomass productivity, as proven by Jeon et al. [23] and as also highlighted in our country assessment. Results of RF-based models can be applied to plan agricultural land use in order to increase the yield and make it sustainable, without environmental side effects. One of the most important and interesting results in our perspective is the quantification of the relative importance of explanatory variables, which best reflects the different edaphic and climatic needs of the observed crop species. For wheat, soil characteristics are the most important factors, while temperature and precipitation are less important [66]. In case of maize, soil parameters are still important but temperature and precipitation have more importance than in the case of wheat, highlighting that, even in a relatively small country like Hungary, climate tolerance of plants is a differentiating factor. This observation becomes more evident when studying sunflowers, where the importance of mean temperature and precipitation outweighs those of soil type and soil textures as earlier presented by Kern et al. (2018) in case studies from Hungary. Nevertheless climatic variables, such as precipitation in October and November and temperature in January and February are also important for winter wheat [66]. Our results also show the importance of summer rainfall totals (May, June, July) for maize, while for sunflowers the most important parameters are spring and autumn temperatures. We have to emphasize that it is often difficult to compare our results with those of other researchers, because the bioclimatic variables of the study area differ. For example, the work of Vannoppen and Gobin (2018) from northern Belgium, investigating the importance of climatic variables in winter wheat yield estimation, found similar parameters to be important, but in a different order. While in Hungary, the mean temperature in January and the amount of precipitation in November are the most important, in Belgium, winter precipitation is the most important [67]. The model fit can be further improved by adding information on management factors such as nutrient levels and fertilizer inputs [52,68].
Soil plays an important role in increasing crop production. The soil science community is trying to define the appropriate indicators. The presented analysis on the importance of variables in calculating productivity also provides a good basis for SDG indicators, as the related target of SDG is to improve land and soil quality progressively. Addressing soil health and soil quality are the main criteria for achieving sustainable agriculture. Climate change largely affects the minimum and maximum temperatures and the amount of precipitation per month [69,70,71,72]. Our results suggest that these variables are also important for winter wheat, maize and sunflowers, and that changes in these variables could change soil productivity in the future.
We established a baseline prediction model for biomass productivity applicable for Hungarian croplands using Earth observation data and yield statistics, identified the importance of soil and climatic determinants of biomass productivity, and proposed a methodology for integrated monitoring of biomass productivity.

5. Conclusions

Our present assessment shows the long-term productivity of soils in Hungary. Long-term productivity in this context means the mean productivity of the last three decades. A period of 20 to 30 years was found to be an adequate time scale for estimating the productivity of soil biomass and for identifying the weights of different factors in it, and also as prospective baseline and threshold values of soil health and soil quality indicators, which can be used in land degradation and soil improvement assessment. A new generalized biomass productivity map was created on a 100 m resolution, which can be implemented in the cadastral system and in multipurpose land use planning programs. The general map of productivity was produced from crop-specific productivity maps by applying biomass productivity models on the country-scale soil, climate and topography geodatabase. Soil properties and characteristics play the most important roles in wheat biomass productivity, while maize has a more significant relationship with precipitation. In the case of sunflowers, soil type and texture are less important factors. The spatial pattern of biomass potential is shown on the general productivity map at 100 m resolution. This map can be used to plan land use in general and agricultural production in croplands. Climate change largely affects the minimum and maximum temperatures, their variability and the amount of precipitation and its temporal distribution, which all have considerable impact on soil biomass productivity. The most important climatic variables for crops deserve particular attention in the next decade, particularly in developing adaptation strategies. We believe that our soil–climate-based land productivity models will help in developing new methods for such adaptation. However, in order to measure changes in biomass production potential, further assessment is required, including trend analysis and the analysis of the effects of changing combinations of soil properties. Nevertheless, the proposed methodology, in addition to possible applications in cadastral systems and in land use planning and agricultural development programs, is also applicable to the integrated monitoring of biomass productivity, which is in line with the goals related to the UN SDGs.

Author Contributions

Conceptualization, G.T., N.C. and B.S.; methodology, G.T., B.S., N.C. and G.S.; resources, L.P. and A.L.; writing—original draft preparation, B.S., G.T. and N.C.; writing—review and editing, N.C., B.S., T.H., A.L., J.M., L.P., G.S., K.T. and G.T.; visualization, N.C.; supervision, B.S. and G.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to private property.

Acknowledgments

We thank the anonymous referees for their valuable recommendations and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Coefficients to modify the computed productivity based on slope relief and orientation [49].
Table A1. Coefficients to modify the computed productivity based on slope relief and orientation [49].
Slope (%)South,
South-West
West,
South-East
East,
North-West
North-EastNorth
111110.98
21110.980.96
3110.980.960.94
410.980.960.940.92
50.980.960.940.920.9
60.960.940.920.90.88
70.940.920.90.880.86
80.920.90.880.860.84
90.90.880.860.840.82
100.880.860.840.820.8
110.860.840.820.80.78
120.840.820.80.780.76
130.820.80.780.760.74
140.80.780.760.740.72
150.780.760.740.720.7
160.760.740.720.70.68
170.740.720.70.680.66
180.720.70.680.660.64
190.70.680.660.640.62
200.680.660.640.620.6
210.660.640.620.60.58
220.640.620.60.580.56
230.620.60.580.560.54
240.60.580.560.540.52
250.580.560.540.520.5
250.560.540.520.50.48

References

  1. United Nations’ Agenda Sustainable Development Goals (SDGs). Available online: https://www.undp.org/sustainable-development-goals?utm_source=EN&utm_medium=GSR&utm_content=US_UNDP_PaidSearch_Brand_English&utm_campaign=CENTRAL&c_src=CENTRAL&c_src2=GSR&gclid=Cj0KCQjw2MWVBhCQARIsAIjbwoPU19Uvs4z3V0arAu_3QfBuDppDDaLFi5wfsb9husx4Hdj7FZ62lk (accessed on 4 December 2022).
  2. Keesstra, S.D.; Bouma, J.; Wallinga, J.; Tittonell, P.; Smith, P.; Cerdà, A.; Montanarella, L.; Quinton, J.N.; Pachepsky, Y.; Van Der Putten, W.H.; et al. The significance of soils and soil science towards realization of the United Nations sustainable development goals. SOIL 2016, 2, 111–128. [Google Scholar] [CrossRef]
  3. Bouma, J.; Montanarella, L.; Evanylo, G. The challenge for the soil science community to contribute to the implementation of the UN Sustainable Development Goals. Soil Use Manag. 2019, 35, 538–546. [Google Scholar] [CrossRef]
  4. Sims, N.C.; Newnham, G.J.; England, J.R.; Guerscham, J.; Cox, S.J.D.; Roxburgh, S.H.; Viscara Rossel, R.A.; Fritz, S.; Wheeler, I. Good Practice Guidance; SDG Indicator 15.3.1, Proportion of Land That Is Degraded Over Total Land Area Version 2.0; United Nations Conventions to Combat Desertification: Bonn, Germany, 2021. [Google Scholar]
  5. UNCCD SDG Indicator 15.3. Available online: https://knowledge.unccd.int/ldn/ldn-monitoring/sdg-indicator-1531 (accessed on 4 December 2022).
  6. McBride, R.A.; Bober, M.L. Quantified evaluation of agricultural soil capability at the local scale: A GIS-assisted case study from Ontario, Canada. Soil Use Manag. 1993, 9, 58–65. [Google Scholar] [CrossRef]
  7. Rossiter, D.G. A theoretical framework for land evaluation. Geoderma 1996, 72, 165–190. [Google Scholar] [CrossRef]
  8. FAO. A Framework for Land Evaluation; Food and Agriculture Organization: Rome, Italy, 1976; p. 32. [Google Scholar]
  9. Kumar, R.; Kalbende, A.R.; Landey, R.J. Soil Evaluation for Agricultural Land Use—II. Productivity Potential Appraisal. J. Indian Soc. Soil Sci. 1984, 32, 467–472. [Google Scholar]
  10. Schroers, J.O. Zur Entwicklung der Landnutzung auf Grenzstandorten in Abhängikeit Agrarmarktpolitischer, Agrarstrukturpolitischer und Produktions-Technologischer Rahmenbedingungen; Justus-Liebig-Universität: Gießen, Germany, 2006. [Google Scholar]
  11. Esch, E.; Mccann, K.; Kamm, C.; Arce, B.; Carroll, O.; Dolezal, A.; Mazzorato, A.; Noble, D.; Fraser, E.; Fryxell, J.; et al. Rising farm costs, marginal land cropping, and ecosystem service markets. Preprint 2021. [Google Scholar] [CrossRef]
  12. Esch, E.; MacDougall, A.S.; Esch, E.; MacDougall, A.S. More at the Margin: Leveraging ECOSYSTEM Services on Marginal Lands to Improve Agricultural Sustainability and Slow Trends of Farming Costs Outpacing Yield Gains; American Geophysical Union: Washington, DC, USA, 2018. [Google Scholar]
  13. Gopalakrishnan, G.; Negri, M.C.; Snyder, S.W. A Novel Framework to Classify Marginal Land for Sustainable Biomass Feedstock Production. J. Environ. Q. 2011, 40, 1593–1600. [Google Scholar] [CrossRef]
  14. Riquier, J.; Bramao, D.L.; Cornet, J.P. A new system of soil appraisal in terms of actual and potential productivity. FAO Soil Resour. Dev. Conserv. 1970, 38, 31–33. [Google Scholar]
  15. Sys, C. Land Evaluation. I-II-III; State University of Ghent: Ghent, Belgium, 1985. [Google Scholar]
  16. Van Lanen, H.A.J.; Broeke, M.J.D.H.; Bouma, J.; de Groot, W.J.M. A mixed qualitative/quantitative physical land evaluation methodology. Geoderma 1992, 55, 37–54. [Google Scholar] [CrossRef]
  17. Tóth, G. Evaluation of cropland productivity in Hungary with the D-e-Meter land evaluation system. Agrokémia Talajt 2011, 60, 161–174. [Google Scholar]
  18. Tóth, G.; Gardi, C.; Bódis, K.; Ivits, É.; Aksoy, E.; Jones, A.; Jeffrey, S.; Petursdottir, T.; Montanarella, L. Continental-scale assessment of provisioning soil functions in Europe. Ecol. Process. 2013, 2, 32. [Google Scholar] [CrossRef]
  19. Burrough, P.A. Fuzzy mathematical methods for soil survey and land evaluation. J. Soil Sci. 1989, 40, 477–492. [Google Scholar] [CrossRef]
  20. Godev, G.; Klestov, V. Statistical evaluation of soil fertility at given plant environment system. In Proceedings of the UNDP/FAO 472 Meeting of Panel of Experts on Land Productivity Evaluation, Sofia, Bulgaria, 27 September 1971. [Google Scholar]
  21. Trashliev, H.; Godev, G.; Krastanov, S.; Klevstov, A.; Kabakchiev, I.; Hershkovich, E.; Dilkov, D. Assessment of ecological conditions for wheat and maize in Bulgaria by means of multivariate regression analysis. In Proceedings of the UNDP/FAO 472 Meeting of Panel of Experts on Land Productivity Evaluation, Sofia, Bulgaria, 27 September 1971. [Google Scholar]
  22. Bonfante, A.; Terribile, F.; Bouma, J. Refining physical aspects of soil quality and soil health when exploring the effects of soil degradation and climate change on biomass production: An Italian case study. SOIL 2019, 5, 1–14. [Google Scholar] [CrossRef]
  23. Jeong, J.H.; Resop, J.P.; Mueller, N.D.; Fleisher, D.H.; Yun, K.; Butler, E.E.; Timlin, D.J.; Shim, K.M.; Gerber, J.S.; Reddy, V.R.; et al. Random Forests for Global and Regional Crop Yield Predictions. PLoS ONE 2016, 11, e0156571. [Google Scholar] [CrossRef]
  24. Roell, Y.E.; Beucher, A.; Møller, P.G.; Greve, M.B.; Greve, M.H. Comparing a Random Forest Based Prediction of Winter Wheat Yield to Historical Yield Potential. Agronomy 2020, 10, 395. [Google Scholar] [CrossRef]
  25. Cheng, M.; Jiao, X.; Shi, L.; Penuelas, J.; Kumar, L.; Nie, C.; Wu, T.; Liu, K.; Wu, W.; Jin, X. High-resolution crop yield and water productivity dataset generated using random forest and remote sensing. Sci. Data 2022, 9, 641. [Google Scholar] [CrossRef] [PubMed]
  26. Guo, Y.; Xia, H.; Zhao, X.; Qiao, L.; Du, Q.; Qin, Y. Early-season mapping of winter wheat and garlic in Huaihe basin using Sentinel-1/2 and Landsat-7/8 imagery. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2022, 1–10. [Google Scholar] [CrossRef]
  27. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  28. Cutler, A.; Cutler, D.R.; Stevens, J.R. Random Forests. In Ensemble Machine Learning; Springer: Boston, MA, USA, 2012; pp. 157–175. [Google Scholar] [CrossRef]
  29. Farkas, J.Z.; Lennert, J. Modelling and predicting of the land use change in Hungary. In Climate Change—Society—Economy: Long-Term Processes and Trends in Hungary; Czirfusz, M., Hoyk, E., Suvák, A., Eds.; Publikon: Pécs, Hungary, 2015; pp. 193–221. ISBN 978-615-5457-62-3. [Google Scholar]
  30. EEA. Corine Land Cover (CLC) 2018; EEA: Copenhagen, Denmark, 2018. [Google Scholar]
  31. Kocsis, M.; Tóth, G.; Berényi-Üveges, J.; Makó, A. Presentation of soil data from the National Pedological and Crop Production Database (NPCPD) and investigations on spatial representativeness. Agrokémia Talajt. 2014, 63, 223–248. [Google Scholar] [CrossRef]
  32. Stefanovits, P.; Michéli, E. Talajgenetika, Talajosztályozás II; Gödöllői Agrártudományi Egyetem: Gödöllő, Hungary, 1989. [Google Scholar]
  33. FAO. World Reference Base for Soil Resources. In International Soil Classification System for Naming Soils and Creating Legends for Soil Maps; FAO: Rome, Italy, 2015. [Google Scholar]
  34. Jin, H.; Eklundh, L. A physically based vegetation index for improved monitoring of plant phenology. Remote Sens. Environ. 2014, 152, 512–525. [Google Scholar] [CrossRef]
  35. NASA Modis Dataset. Available online: https://modis.gsfc.nasa.gov/ (accessed on 4 December 2022).
  36. He, M.; Kimball, J.S.; Maneta, M.P.; Maxwell, B.D.; Moreno, A.; Beguería, S.; Wu, X. Regional Crop Gross Primary Productivity and Yield Estimation Using Fused Landsat-MODIS Data. Remote Sens. 2018, 10, 372. [Google Scholar] [CrossRef]
  37. Dobor, L.; Barcza, Z.; Hlásny, T.; Havasi; Horváth, F.; Ittzés, P.; Bartholy, J. Bridging the gap between climate models and impact studies: The FORESEE Database. Geosci. Data J. 2015, 2, 1–11. [Google Scholar] [CrossRef]
  38. SRTM Shuttle Radar Topography Mission (SRTM) Global. Available online: https://portal.opentopography.org/datasetMetadata?otCollectionID=OT.042013.4326.1 (accessed on 4 December 2022).
  39. Rabus, B.; Eineder, M.; Roth, A.; Bamler, R. The shuttle radar topography mission—A new class of digital elevation models acquired by spaceborne radar. ISPRS J. Photogramm. Remote Sens. 2003, 57, 241–262. [Google Scholar] [CrossRef]
  40. Józsa, E.; Fábián, S.Á.; Kovács, M. An evaluation of EU-DEM in comparison with ASTER GDEM, SRTM and contour-based DEMs over the Eastern Mecsek Mountains. Hungarian Geogr. Bull. 2014, 63, 401–423. [Google Scholar] [CrossRef]
  41. Pásztor, L.; Laborczi, A.; Bakacsi, Z.; Szabó, J.; Illés, G. Compilation of a national soil-type map for Hungary by sequential classification methods. Geoderma 2018, 311, 93–108. [Google Scholar] [CrossRef]
  42. Pásztor, L.; Laborczi, A.; Takács, K.; Szatmári, G.; Dobos, E.; Illés, G.; Bakacsi, Z.; Szabó, J. Compilation of novel and renewed, goal oriented digital soil maps using geostatistical and data mining tools. Hungarian Geogr. Bull. 2015, 64, 49–64. [Google Scholar] [CrossRef]
  43. Pásztor, L.; Laborczi, A.; Takács, K.; Szatmári, G.; Fodor, N.; Illés, G.; Farkas-Iványi, K.; Bakacsi, Z.; Szabó, J. Compilation of Functional Soil Maps for the Support of Spatial Planning and Land Management in Hungary. In Soil Mapping and Process Modeling for Sustainable Land Use Management; Pereira, P., Brevik, E.C., Munoz-Rojas, M., Miller, B.A., Eds.; Elsevier: Amsterdam, The Netherlands, 2017; pp. 293–317. [Google Scholar]
  44. Pásztor, L.; Laborczi, A.; Takács, K.; Illés, G.; Szabó, J.; Szatmári, G. Progress in the elaboration of GSM conform DSM products and their functional utilization in Hungary. Geoderma Reg. 2020, 21, e00269. [Google Scholar] [CrossRef]
  45. Baranyai, F.; Fekete, A.; Kovács, I. A Magyarországi Tápanyag-Vizsgálatok Eredményei; Mezőgazdasági Kiadó: Budapest, Hungary, 1987; pp. 8–79. [Google Scholar]
  46. KSH Központi Statisztikai Hivatal (Central Statistical Office). Available online: https://www.ksh.hu (accessed on 4 December 2022).
  47. Farmanov, N.; Amankulova, K.; Szatmari, J.; Sharifi, A.; Abbasi-Moghadam, D.; Mirhossein-Nejad, M.; Mucsi, L. Crop Type Classification by DESIS Hyperspectral Imagery and Machine Learning Algorithms. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2023, 16, 1576–1588. [Google Scholar] [CrossRef]
  48. Wright, M.N.; Ziegler, A. Ranger: A Fast Implementation of Random Forests for High Dimensional Data in C++ and R. J. Stat. Softw. 2017, 77, 1–17. [Google Scholar] [CrossRef]
  49. Mezőgazdasági és Élelmezésügyi Minisztérium (MÉM) Táblázatok a Földértékelés Végrehajtásához; MÉM: Budapest, Hungary, 1982.
  50. Fekete, Z. Direktívák a Gyakorlati Földértékeléshez; Mezőgazdasági Kiadó: Budapest, Hungary, 1965. [Google Scholar]
  51. Fórizsné, J.; Máté, F.; Stefanovits, P. Talajbonitáció—Földértékelés. MTA Agrártudományok Osztályának Közleményei 1972, 30, 359–378. [Google Scholar]
  52. Hermann, T.; Kismányoky, T.; Tóth, G. A humuszellátottság hatása a kukorica (Zea mays L.) termésére csernozjom és barna erdőtalajú termőhelyeken, különböző évjáratokban. Növénytermelés 2014, 63, 1–18. [Google Scholar]
  53. Máté, F. Megjegyzések a talajok termékenységük szerinti osztályozásához. Agrokémia Talajt. 1960, 9, 419–426. [Google Scholar]
  54. Tóth, G.; Máté, F. Jellegzetes dunántúli talajok főbb növényenkénti relatív termékenysége. Agrokémia Talajt. 1999, 48, 172. [Google Scholar]
  55. Nia, V.P.; Davison, A.C. High-Dimensional Bayesian Clustering with Variable Selection: The R Package bclust. J. Stat. Softw. 2012, 47, 1–22. [Google Scholar] [CrossRef]
  56. Willmott, C.; Matsuura, K. Advantages of the mean absolute error (MAE) over the root mean square error (RMSE) in assessing average model performance. Clim. Res. 2005, 30, 79–82. [Google Scholar] [CrossRef]
  57. De Myttenaere, A.; Golden, B.; Grand, B.L.; Rossi, F. Mean Absolute Percentage Error for regression models. Neurocomputing 2016, 192, 38–48. [Google Scholar] [CrossRef]
  58. Seddaiu, G.; Iocola, I.; Farina, R.; Orsini, R.; Iezzi, G.; Roggero, P.P. Long term effects of tillage practices and N fertilization in rainfed Mediterranean cropping systems: Durum wheat, sunflower and maize grain yield. Eur. J. Agron. 2016, 77, 166–178. [Google Scholar] [CrossRef]
  59. Sadras, V.O.; Calviño, P.A. Quantification of Grain Yield Response to Soil Depth in Soybean, Maize, Sunflower, and Wheat. Agron. J. 2001, 93, 577–583. [Google Scholar] [CrossRef]
  60. Feng, P.; Wang, B.; Liu, D.L.; Waters, C.; Xiao, D.; Shi, L.; Yu, Q. Dynamic wheat yield forecasts are improved by a hybrid approach using a biophysical model and machine learning technique. Agric. For. Meteorol. 2020, 285–286, 107922. [Google Scholar] [CrossRef]
  61. Kheir, A.M.S.; Ammar, K.A.; Amer, A.; Ali, M.G.M.; Ding, Z.; Elnashar, A. Machine learning-based cloud computing improved wheat yield simulation in arid regions. Comput. Electron. Agric. 2022, 203, 107457. [Google Scholar] [CrossRef]
  62. Reeves, M.C.; Zhao, M.; Running, S.W. Usefulness and limits of MODIS GPP for estimating wheat yield. Int. J. Remote Sens. 2007, 26, 1403–1421. [Google Scholar] [CrossRef]
  63. Mishra, U.; Torn, M.S.; Fingerman, K. Miscanthus biomass productivity within US croplands and its potential impact on soil organic carbon. GCB Bioenergy 2013, 5, 391–399. [Google Scholar] [CrossRef]
  64. Xie, X.; Li, A.; Tian, J.; Wu, C.; Jin, H. A fine spatial resolution estimation scheme for large-scale gross primary productivity (GPP) in mountain ecosystems by integrating an eco-hydrological model with the combination of linear and non-linear downscaling processes. J. Hydrol. 2023, 616, 128833. [Google Scholar] [CrossRef]
  65. Sabetraftar, K.; Mackey, B.; Croke, B. Sensitivity of modelled gross primary productivity to topographic effects on surface radiation: A case study in the Cotter River Catchment, Australia. Ecol. Modell. 2011, 222, 795–803. [Google Scholar] [CrossRef]
  66. Kern, A.; Barcza, Z.; Marjanović, H.; Árendás, T.; Fodor, N.; Bónis, P.; Bognár, P.; Lichtenberger, J. Statistical modelling of crop yield in Central Europe using climate data and remote sensing vegetation indices. Agric. For. Meteorol. 2018, 260–261, 300–320. [Google Scholar] [CrossRef]
  67. Vannoppen, A.; Gobin, A. Estimating Farm Wheat Yields from NDVI and Meteorological Data. Agron. 2021, 11, 946. [Google Scholar] [CrossRef]
  68. Debreczeni, B.; Németh, T.; Tóth, G. Nutrient factor of land quality. In Land Evaluation and Land Use Information; Gaál, Z., Máté, F., Tóth, G., Eds.; Veszprémi Egyetem: Keszthely, Hungary, 2003; pp. 39–48. ISBN 963949525. [Google Scholar]
  69. Farooq, M.; Bramley, H.; Palta, J.A.; Siddique, K.H.M. Heat Stress in Wheat during Reproductive and Grain-Filling Phases. Crit. Rev. Plant Sci. 2011, 30, 491–507. [Google Scholar] [CrossRef]
  70. Hatfield, J.L.; Boote, K.J.; Kimball, B.A.; Ziska, L.H.; Izaurralde, R.C.; Ort, D.; Thomson, A.M.; Wolfe, D. Climate Impacts on Agriculture: Implications for Crop Production. Agron. J. 2011, 103, 351–370. [Google Scholar] [CrossRef]
  71. Lobell, D.B.; Schlenker, W.; Costa-Roberts, J. Climate trends and global crop production since 1980. Science 2011, 333, 616–620. [Google Scholar] [CrossRef]
  72. Lobell, D.B.; Cahill, K.N.; Field, C.B. Historical effects of temperature and precipitation on California crop yields. Clim. Change 2007, 81, 187–203. [Google Scholar] [CrossRef]
Figure 1. Arable land areas of Hungary (study area), based on Corine Land Cover 2018 dataset [30].
Figure 1. Arable land areas of Hungary (study area), based on Corine Land Cover 2018 dataset [30].
Remotesensing 15 01236 g001
Figure 2. Flowchart of land evaluation modeling process.
Figure 2. Flowchart of land evaluation modeling process.
Remotesensing 15 01236 g002
Figure 3. Scatter plot of observed vs estimated biomass productivity of total cropland area (A), wheat (B), maize (C) and sunflowers (D). Results were significant at the 0.01 level.
Figure 3. Scatter plot of observed vs estimated biomass productivity of total cropland area (A), wheat (B), maize (C) and sunflowers (D). Results were significant at the 0.01 level.
Remotesensing 15 01236 g003
Figure 4. Croplands’ land evaluation values range between 0 and 100 in the case of general biomass productivity of arable lands without slope correction coefficients (A) and separately wheat (B), maize (C) and sunflowers (D).
Figure 4. Croplands’ land evaluation values range between 0 and 100 in the case of general biomass productivity of arable lands without slope correction coefficients (A) and separately wheat (B), maize (C) and sunflowers (D).
Remotesensing 15 01236 g004
Figure 5. Croplands’ land evaluation values range between 0 and 100 in the case of general biomass productivity of arable lands after slope correction.
Figure 5. Croplands’ land evaluation values range between 0 and 100 in the case of general biomass productivity of arable lands after slope correction.
Remotesensing 15 01236 g005
Figure 6. Overall importance of explanatory variables in predicting wheat (A), maize (B) and sunflower (C) biomass productivity.
Figure 6. Overall importance of explanatory variables in predicting wheat (A), maize (B) and sunflower (C) biomass productivity.
Remotesensing 15 01236 g006
Table 1. Main features of the AIIR dataset, based on Hungarian [32] and World Reference Base for Soil Resources [33] classification.
Table 1. Main features of the AIIR dataset, based on Hungarian [32] and World Reference Base for Soil Resources [33] classification.
Soil Taxonomical Unit of Major Agricultural SoilsNo. of Parcels CoveredArea (ha)Area (%)
Hungarian classificationWRB 2014
Lessivated brown forest soil (non-podzolic)Haplic Luvisol11,062385,04810.06
Raman-type brown forest soilHaplic Cambisol6567270,2397.06
Rust-brown sandy forest soilArenic Cambisol2988114,8723
Typical calcareous chernozemHaplic Chernozems3792228,2405.96
Great Plains calcareous chernozemHaplic Chernozems2042120,1233.14
Carbonated meadow chernozemGleyic Chernozems5540330,2008.63
Non-carbonated meadow chernozemLuvic Chernozems2021108,1492.83
Carbonated meadow soilCalcic Vertisols3952184,8534.83
Non-carbonated meadow soilHaplic Vertisols3460151,3943.96
Carbonated alluvial meadow soilGleyic Fluvisols3129142,5353.73
Non-carbonated alluvial meadow soilDystric Fluvisols4658179,1014.68
Carbonated humic alluvial soilCalcaric Fluvisols121051,7201.35
Non-carbonated humic alluvial soilDystric Fluvisols158450,7891.33
Carbonated humic sandy soilCalcaric Cambisols3714138,0443.61
Non-carbonated humic sandy soilDystric Cambisols245875,6561.98
major soils in total58,1772,530,96366.2
other soils28,5171,295,46733.8
86,6953,826,430100
Table 2. Test validation results of all cropland, wheat, maize and sunflowers. R2: correlation coefficient, R: Pearson correlation, MAPE: mean absolute percentage error, MAE: mean absolute error, N: number of pairs.
Table 2. Test validation results of all cropland, wheat, maize and sunflowers. R2: correlation coefficient, R: Pearson correlation, MAPE: mean absolute percentage error, MAE: mean absolute error, N: number of pairs.
R2RMAPE (%)MAEN
All cropland0.40.6319.287.334381
Wheat0.410.6418.066.782631
Maize0.350.5919.177.931646
Sunflower0.270.5229.8111.7104
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Csikós, N.; Szabó, B.; Hermann, T.; Laborczi, A.; Matus, J.; Pásztor, L.; Szatmári, G.; Takács, K.; Tóth, G. Cropland Productivity Evaluation: A 100 m Resolution Country Assessment Combining Earth Observation and Direct Measurements. Remote Sens. 2023, 15, 1236. https://doi.org/10.3390/rs15051236

AMA Style

Csikós N, Szabó B, Hermann T, Laborczi A, Matus J, Pásztor L, Szatmári G, Takács K, Tóth G. Cropland Productivity Evaluation: A 100 m Resolution Country Assessment Combining Earth Observation and Direct Measurements. Remote Sensing. 2023; 15(5):1236. https://doi.org/10.3390/rs15051236

Chicago/Turabian Style

Csikós, Nándor, Brigitta Szabó, Tamás Hermann, Annamária Laborczi, Judit Matus, László Pásztor, Gábor Szatmári, Katalin Takács, and Gergely Tóth. 2023. "Cropland Productivity Evaluation: A 100 m Resolution Country Assessment Combining Earth Observation and Direct Measurements" Remote Sensing 15, no. 5: 1236. https://doi.org/10.3390/rs15051236

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