Next Article in Journal
Tea Tree Pest Detection Algorithm Based on Improved Yolov7-Tiny
Next Article in Special Issue
Spatial Prediction and Mapping of Soil Water Content by TPE-GBDT Model in Chinese Coastal Delta Farmland with Sentinel-2 Remote Sensing Data
Previous Article in Journal
Combinatorial Olfactory Signaling in Short-Distance Determines Host Plant Recognition in Locust
Previous Article in Special Issue
A New Method for Estimating Irrigation Water Use via Soil Moisture
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing the Within-Field Heterogeneity Using Rapid-Eye NDVI Time Series Data

1
Crop Science Group, Institute of Crop Science and Resource Conservation (INRES), University of Bonn, Katzenburgweg 5, 53115 Bonn, Germany
2
Research Center Jülich, Wilhelm-Johnen-Straße, 52428 Jülich, Germany
3
Department of Agricultural Sciences, University of Helsinki, Koetilantie 5, FI-00014 Helsinki, Finland
*
Author to whom correspondence should be addressed.
Agriculture 2023, 13(5), 1029; https://doi.org/10.3390/agriculture13051029
Submission received: 18 February 2023 / Revised: 4 May 2023 / Accepted: 6 May 2023 / Published: 9 May 2023
(This article belongs to the Special Issue Remote Sensing Technologies in Agricultural Crop and Soil Monitoring)

Abstract

:
(1) Background: The relation between the sub-field heterogeneity of soil properties and high-resolution satellite time series data might help to explain spatiotemporal patterns of crop growth, but detailed field studies are seldom. (2) Methods: Normalized Difference Vegetation Index (NDVI) data derived from satellite time series images were used to identify changes in the spatial distribution of winter triticale (×Triticosecale), winter rye (Secale cereale) and winter barley (Hordeum vulgare) growth (2015 to 2020) for a field in north-eastern Germany. NDVI patterns (quartiles) that remained persistent over time were identified and it was tested if spatially heterogeneous soil characteristics such as water holding capacity and altitude could explain them. (3) Results: A statistically significant relationship between elevation and soil classes with NDVI values was found in most cases. The lowest NDVI quartiles, considered as representing the poorest growth conditions, were generally found in the depressions with the lowest water holding capacity. These areas showed temporally stable spatial patterns, especially during the pre-harvest period. Over the six-year period, up to 80% of the grid cells with the lowest NDVI values were spatially consistent over time. Differences in the climatic water balance were rather low but could contribute to explaining spatial patterns, such as the lower clustering of values in the wettest year. (4) Conclusions: High-resolution satellite NDVI time series are a valuable information source for precise land management in order to optimize crop management with respect to yield and ecosystem services.

1. Introduction

Heterogeneous growing conditions within a field plot can benefit from modern precision agriculture in terms of spatially optimized management [1] and illustrate the high importance of acknowledging the interaction between the soil conditions, topography and cereal production [2,3,4,5,6,7]. For an efficient and non-destructive assessment of the heterogeneity of soil conditions of arable land and plant growth characteristics, the use of remote sensing data is becoming increasingly established [8,9]. Plant growth across landscapes or within a single field can show temporally consistent spatial patterns, such as low- and high-yielding areas. Depending on the spatial resolution, vegetation indices (VI) allow for assessing such patterns and their change over time. Several studies have successfully used vegetation indices based on remote sensing data to monitor and predict plant growth at the field scale [4,5,10,11,12]. However, the direct link between within-field changes in soil properties, topography and plant growth and their change over time has seldom been studied in more detail. It therefore often remains unclear if spatial growth differences strongly vary over a longer time period and are governed by external factors such as meteorological conditions, or are mainly influenced by intrinsic factors, such as soil parameters, which would lead to temporally consistent spatial growth patterns.
This study aims at providing an example and a workflow for using high-resolution satellite time series (Rapid Eye, Planet Scope) to assess spatially consistent vegetation patterns and their relation to soil characteristics at the field scale. The test site is an arable field of about six ha located in the state of Brandenburg, Germany. Using satellite data from the years 2015 to 2020 the following research hypotheses were tested: (1) the spatial distribution of Normalized Difference Vegetation Index (NDVI) mainly depends on soil heterogeneity and is consistent over time, (2) annual weather conditions can slightly modify the spatial pattern of NDVI. Years with dry conditions lead to more pronounced spatial patterns.

2. Materials and Methods

The study area is located in Großmutz, Brandenburg, Germany, between continental (low precipitation, low humidity) and maritime (high precipitation, high humidity) climates [13] and is classified as a moderately dry lowland climate. The long-term averages between 1991 and 2020 (WHO climate normal) of the closest climate station of the German Weather Service (DWD) in Neuruppin-Alt Rupin (Station ID: 96) are 623 mm of annual precipitation and 9.6 °C mean air temperature [14,15]. The field has an area of six ha and is characterized by spatial heterogeneity in terms of soil properties, topography and crop growth [16]. The altitude decreases from the southeastern corner (58 m a.s.l.) to the northern corner (51 m a.s.l.) (Figure 1). A soil reconnaissance map by [17] classifies the investigated area as an anhydromorphic substrate type “ASL w” (sandy ground moraine with undulated relief). The dominant soil types can be classified as Arenic Cambisol or Haplic Luvisol [18]. The agricultural field is conventionally managed and not irrigated. The cereal crops cultivated during the observation period are listed in Table 1. Summer oat (harvested in 2018) was excluded from this study to facilitate interannual comparisons between winter cereals.
The meteorological data used in this study were based on station data provided by the DWD (Germany’s national meteorological service (Deutscher Wetterdienst, DWD)) [18]. The climatic water balance (CWB in mm) was calculated for each day by subtracting the potential evaporation from the precipitation. Daily water balances were summed up to monthly values.
Soil samples were taken up to one m depth with a Pürckhauer soil auger along nine transect lines as described in [16]. The transect orientation followed the slope from the southwestern side to the northern corner with a spatial distance of 50 m from each other. The outer transects were 25 m apart from the edge of the field. Within the transects, samples were taken in 15 m intervals and each coordinate of the 87 sampling points was recorded with a Trimble Geo 7× Handheld GNSS System (Sunnyvale, CA, USA). Following the Munsell soil color charts [19] and the FAO (Food and Agriculture Organization of the United Nations) guidelines for soil description [20], the physical properties of the soil profiles were estimated in moist conditions by manually wetting soil probes. In general, all soil profiles down to 1 m depth were composed of a layer with sandy soil material (predominantly Loamy Sand or Sand) with variable thickness (35 to 100 cm) and a subsequent loam layer (Sandy loam to Loam) with variable thickness (0–65 cm). Based on the spatial variability of the vertical position (depth) of the loamy layer, three soil classes were derived as an indicator for differences in the available water capacity (AWC) of the soils [16]. Soil class 1 (S1) was characterized by a shallow upper boundary of the loamy subsoil layer (<60 cm depth), class 2 (S2), a deeper upper boundary of the loamy layer at 60–100 cm depth and class 3 (S3) was sandy throughout 100 cm soil depth (no loamy layer). Assuming that sandy soils have a low water holding capacity [21] that is improved by the presence of an underlying loamy layer, soil class 1 is assumed to show the highest and class 3 the lowest AWC. An interpolated soil map was created by kriging the soil sample points. To take the different altitude levels into account, which describe the position of the drilling points along the slopes (catena) at our site, a freely available digital terrain model (DTM) provided by the Amt für Landesvermessung und Geobasisinformation Brandenburg (LGB, https://geobasis-bb.de/lgb/de/geodaten/3d-produkte/gelaendemodell/, last visited: 17 October 2022) was used to further distinguish between three altitude classes: A1: <53.5 m, A2: 53.5–55.37 m, A3: >55.37 m a.s.l. Both classifications are shown in Figure 1.
The main type of satellite imagery used in this paper originates from the RapidEye satellite constellation. RapidEye was owned by Planet (Planet Labs, Inc., San Francisco, CA, USA) and the data were provided in cooperation with ESA’s Third-Party Missions Programme. The RapidEye satellite program consisted of 5 satellites collecting data in five spectral bands: red: 630–685 nm, green: 520–590 nm, blue: 440–510 nm, red edge: 690–730 nm and near-infrared (NIR): 760–850 nm. The resolution was 5 × 5 m per grid cell on orthorectified imagery resulting from 6.5 m GSD (Ground Sampling Distance) at the nadir [22]. Only images captured within the main growing season after winter were selected, which was between April and June of each year. Cloud coverage was visually checked. Since RapidEye’s operation was discontinued in spring 2020, PlanetScope imagery was applied as the satellite data source for the 2020 data. PlanetScope was chosen because the images were delivered by the same provider (Planet Labs, Inc., San Francisco, CA, USA), had a similar high temporal and spatial resolution compared to RapidEye and were available free of charge, including an atmospheric correction. The first 12 satellites were launched on 22 June 2016; therefore, PlanetScope was not able to cover the whole investigation period either. The images used were taken with either sensor type PS2 or PS2.SD [22]. The sensor PS2 consists of a four-band frame imager (red: 590–670 nm, green: 500–590 nm, blue: 455–515 nm, NIR: 780–860 nm) with a split-frame for the visible RGB (red/green/blue) and the NIR. PS2.SD consisted of a four-band frame imager (red: 650–682 nm, green: 547–585 nm, blue: 464–517 nm, NIR: 846–888 nm) with butcher-block filter providing separate stripes for blue, green, red and NIR. The GSD of 3.9 m leads to a higher resolution (3 × 3 m). Both RapidEye and PlanetScope provided atmospherically corrected imagery.
To calculate the Normalized Difference Vegetation Index (NDVI) [23,24], the Red (R) and Near Infrared (NIR) bands were used:
NDVI = NIR R / NIR + R
Subsequently, data were limited to the scene with the highest NDVI values in each month (maximum value approach) assuming that these images were least affected by cloud and haze. Image selection was undertaken with a preference for scenes captured in the middle of the month and around the same day each month to allow a constant observation frequency and to increase the comparability across years. However, this was not always feasible due to the limited availability of cloud-free satellite images. The 14 dates used are listed in Table 2.
Due to the described inherent different spatial and spectral resolutions between RapidEye and PlanetScope imagery, and to reduce sources of bias in the NDVI data, PlanetScope data were only used for statistical data analysis when NDVI values were normalized and divided into quartiles. By using normalized NDVI values and a quantile-based approach, we evaluated general trends and patterns and accounted for the uncertainty in NDVI values due to interannual variations in crop phenological development and crop-specific optical properties.
For all statistical analyses and graphs, the statistical software R (version 3.6.3, Windows 10 [25]) was used. Masks were applied to the satellite scenes to select corresponding grid cells for each altitude class (Figure 1) and to subsequently extract information on soil classes (Figure 1) within the altitude classes (package: “raster” [26]). NDVI values were not normally distributed, which was tested with the Shapiro–Wilk test. Therefore, the non-parametric Kruskal–Wallis test was used to test for significant differences in the NDVI values between classes for each scene. Here, an analysis of variance was used to test whether independent samples belong to the same population [27]. As a post-hoc test, the non-parametric Wilcoxon signed-rank test was applied whenever significant differences were found to specify which groups significantly differed by testing two paired samples for their central tendencies [28]. To avoid type 1 errors in the Wilcoxon test, the p-value was adjusted with the Bonferroni method [29].
To compare NDVI values across years, the NDVI values were normalized by dividing all the NDVI values of each grid cell with the mean of all grid cells. Subsequently, the relative NDVI values of each scene could be classified into four different classes. To do so, the values were split into quartiles, and each quartile was assigned its own value (1 to 4). This processing resulted in images showing the classified relative NDVI values, ranging from low (class Q1) to high (class Q4) NDVI values.
First, the non-parametric Mann–Kendall trend test [30] was applied to test for the presence of a monotonic trend in the time series over time (R-package: “Kendall” [31]). Tests were applied for each month separately. This test is commonly used in environmental research [32,33]. Subsequently, a correlation matrix was created by using the “corrplot” package [34]. In the first step, 10% of the pixels per scene were subsampled and pairwise correlations between those pixels were calculated across the years, using the non-parametric Spearman rank correlation using the global “cor” function implemented in R. To quantify the changes in the spatial distribution and size of NDVI clusters between years, cross-tabulations (crosstabs) for all possible combinations of months across all years were computed.

3. Results

3.1. Climate

During the vegetation period, the total precipitation was between 99 and 117 mm. Only in 2017 (220 mm) was the precipitation sum comparatively high and significantly higher than the long-term average value (149 mm). However, this mainly resulted from higher precipitation during the late growing season in June and July, thus, after the observed period used for our study (Figure 2). In general, the climatic water balances (CWB) were rather positive in autumn (October/November) and winter until February/March, followed by negative CWBs during spring and summer (Figure 2). The observation years did not include the main Central European drought year, 2018 [35]. However, spring 2020 was drier compared with the other years in March (CWB = −21 mm). Further, the summed CWB for the months of April to June revealed the driest conditions in 2019 (CWB = −232 mm), followed by 2020 (−205 mm), 2016 (−196 mm), 2015 (−183 mm) and rather wet conditions in 2017 (−57 mm).

3.2. Variation of the NDVI in Altitude and Soil Classes

The total number of grid cells (n) of a single RapidEye scene was 2181. The total number of grid cells per category (soil, altitude) and the resulting area are shown in Table 3. Across all years, the lowest and highest NDVI values were recorded for altitude class A1 in June 2019 (0.47, SE: 0.002) and for A3 in May 2015 (0.88, SE: 0.002), respectively. With respect to the different soil classes, the distribution of the highest (S2, April 2017: 0.87, SE < 0.001) and lowest NDVI values (S3, June 2019: 0.49, SE: 0.002) was slightly less pronounced. In addition, as shown in Table 4, A1 overlapped most frequently with S3 (554 of 751 same grid cells), A2 with S1 (461 of 806same grid cells) and A3 with S1 (393 of 624 same grid cells).
Results from the Kruskal–Wallis test indicated statistically significant differences (p < 0.05) of NDVI values between soil classes and altitude classes. To test for the source of these differences and to identify non-significant class pairs, the Wilcoxon signed-rank test was performed for all scenes. Out of 66 possible pairs of combinations, only 12 were not highly significant (p ≥ 0.05, Table 5). Four of the six non-significant soil class combinations were found between S1 and S2 and four of six of the altitude classes between A2 and A3.

3.3. Spatial Distribution of NDVI over Time

Based on the quartile classification of NDVI data, results showed that the northernmost area of the field had consistently lower NDVI values (most pixels belong to the first quartile Q1 in each scene) (Figure 3 and Figure 4) as compared with the rest of the field. In April and May 2017, the spatial clustering of NDVI values was much lower compared to all other years (Figure 3). A figure with all NDVI quartile maps can be found in the supplementary material (Supplementary Material Figure S1).
The Mann–Kendall trend test ruled out the presence of temporal trends in the NDVI distribution. Comparing monthly data across years using correlation analyses, it was noticeable that, overall, data from April had the lowest correlation values (mean: 0.26, SE: 0.039, n: 10), indicating the highest interannual differences in spatial patterns during early spring. The highest interannual agreement was found for June recordings (mean: 0.64, SE: 0.025, n: 10). The cross-tabulation matrix (Table 6) showed that the category with the lowest NDVI values (Q1) had the highest spatial consistency across years (60–70%), followed by the class with the highest NDVI value (Q4: 42–49%). Consistency in the assignment to the different NDVI quartiles was lowest in April (max. 45%) and highest in June (up to 77%).

4. Discussion

4.1. Spatial Distribution of NDVI

Since the differences in average NDVI values between soil classes as well as between altitude classes were statistically significant in 60 out of 66 possible pairwise comparisons, each soil type and altitude class could be assigned to an average high or low NDVI value for each scene. The fact that both soil type and elevation affected crop growth is consistent with the literature [4,12]. Due to the large spatial overlap of classes, final statements about the ranking of the influence of elevation or soil on NDVI values cannot be made. A study that has been carried out in the same field evaluated the influence of altitude and soil on plant growth as well as harvest yield in 2020 and showed that both factors have a major impact on plant growth depending on the crop growth stage [16]. A combination of soil type S3 (no loamy layer in the subsoil, thus, lower available water capacity) and elevation class A1 (downslope) appears to represent the worst conditions for plant growth and grain yield at harvest [16]. This combination is true for most of the grid cells in the experimental field. These findings were confirmed by our study in terms of low NDVI values in the late growing season, which has been associated with lower grain yields [36,37,38,39]. Results on altitude effects contrast with those from studies where low-lying field sections benefit from their topographic position in terms of water availability [39]. This highlights the need for knowledge of vertical soil profile properties for understanding topography-induced effects on soil moisture and NDVI patterns [40].
A lower influence of the different soil and elevation classes on the NDVI values of rye and thus a spatially more heterogeneous distribution of rye NDVI values in 2017 compared to the other years might be attributed to both, (1) a higher summer water balance compared with the other years and (2) the different response of winter cereals to environmental growing conditions [41]. Winter rye, grown in 2017, is well known as a less demanding crop in terms of environmental conditions as compared with wheat; thus, NDVI patterns are likely less clustered.
In summary, results on the temporally persistent spatial distribution of NDVI classes suggest the presence of generally high- and low-yielding areas. In consequence, a spatially targeted crop management, e.g., with respect to planting density or fertilizer and pesticide application, will most likely allow for more resource-efficient land management at our site.
In general, the spatial scale of Rapid-Eye satellite data was suitable for obtaining information on this area. Using UAV data, future studies might test for an optimum spatial scale to assess the within-field heterogeneity for guiding specific management strategies [39] and help in reducing uncertainty in interannual comparisons caused by species-specific phenological development. This study used data from different winter cereals grown within the same fields within a crop rotation. Plants might have differed in terms of their physiological development and greenness in each year [42]. While this was accounted for by using normalized data and a quantile-based approach, uncertainty resulting from species-specific plant and canopy characteristics cannot be completely ruled out.

4.2. Influence of Weather Conditions on Spatial Patterns of NDVI at Field Scale

Due to the limited number of observation years with rather low differences in the climatic water balance (Figure 2), conclusive statements on the influence of different weather conditions (in particular, drought stress) on the spatial distribution of NDVI values have to be considered with care. Even the driest year (2019, April to June: CWB = −232 mm) did not lead to pronounced differences in the spatial distribution of NDVI values compared to the other years. The year with the highest CWB (2017, −57 mm) showed a more heterogeneous distribution (lowest clustering of NDVI values) (thereby rather confirming). However, in 2017, the field was cultivated with winter rye, which is known for being a robust and low-demanding crop. In consequence, the lower spatial clustering of NDVI values in 2017 could be either due to more humid conditions or due to the higher adaptability of winter rye to changes in soil conditions.
Under Mediterranean climate conditions and generally in arid or semi-arid climates, spatial dynamics in NDVI values have been related to differences in the soil water capacity and to drought indices [43,44,45]. Contrastingly, at our site, under mild temperate conditions with optimum fertilizer supply and using varieties optimized for local conditions, the influence of weather conditions might not be assessed at the sub-field scale for annual cereal crops. Satellite NDVI values are affected by both biomass and greenness. These factors might have differed between crops and varieties each year, thereby adding uncertainty to the results. As an example, in spring, when NDVI values are not yet saturated, nitrogen limitations have been found to be critical for canopy greenness, whereas the importance of soil water availability increases in later summer months. Such observations illustrate the complexity and nonlinearity of environmental effects on NDVI patterns and highlight the importance of spatially detailed studies at the sub-field scale to understand environmental effects on vegetation dynamics as observed from space [39].

5. Conclusions

For a field in North-Eastern Germany, spatial patterns of NDVI values from winter cereals are affected by the spatial distribution of soil available water capacity, as derived from the depth to a loamy layer, and topographic position. The quartile distribution of normalized NDVI values illustrated temporally stable spatial patterns across years, especially in the areas with poor growing conditions, which can be linked to high- and low-yielding zones. Our study highlights the strong influence of the spatial variability of soil and terrain characteristics, especially for those field areas that have poor growth conditions (consistently low NDVI values). Differences in the climatic water balance can contribute to explaining interannual differences in spatial patterns, i.e., clustering of values (with lower clustering in wetter years); however, species-specific differences in the sensitivity to water stress have to be considered. Satellite data provided by RapidEye or PlanetScope with a high spatial resolution are suitable to detect within-field differences in agricultural growing conditions for supporting optimized variable rate management.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agriculture13051029/s1, Figure S1: NDVI values divided into quartiles (1–4). 1 (light yellow) corresponds to quartile 1 (lowest NDVI values) and 4 (dark green) to quartile 4 (highest NDVI values). Matching the year in the headline: (A): April, (B): May, (C): June, (D): April, (E): May, (F) June. Quartile classification of NDVI data showed that the northernmost area of the field had consistently lower NDVI values (Q1).

Author Contributions

Conceptualization, J.M., A.T., H.A. and T.G.; methodology, J.M. and A.T.; software, J.M.; resources, A.T.; data curation, J.M.; writing—original draft preparation, J.M.; writing—review and editing, A.T., H.A. and T.G.; visualization, J.M.; supervision, T.G.; project administration, T.G.; funding acquisition, H.A. and T.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the German Federal Ministry of Education and Research (BMBF) through the Digital Agriculture Knowledge and Information System (DAKIS) Project, grant number 031B0729E, and by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy-EXC 2070-390732324 (PhenoRob). The author held a study grant from the Friedrich Ebert Foundation.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

First: we would like to thank Gunther Krauss for his assistance with questions regarding R and data analysis. Furthermore, J.M. would like to thank Simon Haberstroh for his support during the preparation of the master’s thesis on the topic of the paper. Without both of you, this article would never have been written.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Yuzugullu, O.; Lorenz, F.; Fröhlich, P.; Liebisch, F. Understanding Fields by Remote Sensing: Soil Zoning and Property Mapping. Remote Sens. 2020, 12, 1116. [Google Scholar] [CrossRef]
  2. Kravchenko, A.N.; Thelen, K.D.; Bullock, D.G.; Miller, N.R. Relationship among Crop Grain Yield, Topography, and Soil Electrical Conductivity Studied with Cross-Correlograms. Agron. J. 2003, 95, 1132–1139. [Google Scholar] [CrossRef]
  3. Huang, X.; Wang, L.; Yang, L.; Kravchenko, A.N. Management Effects on Relationships of Crop Yields with Topography Represented by Wetness Index and Precipitation. Agron. J. 2008, 100, 1463–1471. [Google Scholar] [CrossRef]
  4. Thelemann, R.; Johnson, G.; Sheaffer, C.; Banerjee, S.; Cai, H.; Wyse, D. The Effect of Landscape Position on Biomass Crop Yield. Agron. J. 2010, 102, 513–522. [Google Scholar] [CrossRef]
  5. Kumhálová, J.; Kumhála, F.; Kroulík, M.; Matějková, Š. The Impact of Topography on Soil Properties and Yield and the Effects of Weather Conditions. Precis. Agric. 2011, 12, 813–830. [Google Scholar] [CrossRef]
  6. Xue, W.; Huang, L.; Yu, F.-H. Spatial Heterogeneity in Soil Particle Size: Does It Affect the Yield of Plant Communities with Different Species Richness? J. Plant Ecol. 2016, 9, 608–615. [Google Scholar] [CrossRef]
  7. Beuschel, R.; Piepho, H.-P.; Joergensen, R.G.; Wachendorf, C. Similar Spatial Patterns of Soil Quality Indicators in Three Poplar-Based Silvo-Arable Alley Cropping Systems in Germany. Biol. Fertil. Soils 2019, 55, 1–14. [Google Scholar] [CrossRef]
  8. Sahoo, R.; Ray, S.S.; Manjunath, K. Hyperspectral Remote Sensing of Agriculture. Curr. Sci. 2015, 108, 848–859. [Google Scholar]
  9. Kumhálová, J.; Novák, P.; Madaras, M. Monitoring Oats and Winter Wheat Within-Field Spatial Variability by Satellite Images. Sci. Agric. Bohem. 2018, 49, 127–135. [Google Scholar] [CrossRef]
  10. Jiang, P.; Thelen, K.D. Effect of Soil and Topographic Properties on Crop Yield in a North-Central Corn–Soybean Cropping System. Agron. J. 2004, 96, 252–258. [Google Scholar] [CrossRef]
  11. Stadler, A.; Rudolph, S.; Kupisch, M.; Langensiepen, M.; van der Kruk, J.; Ewert, F. Quantifying the Effects of Soil Variability on Crop Growth Using Apparent Soil Electrical Conductivity Measurements. Eur. J. Agron. 2015, 64, 8–20. [Google Scholar] [CrossRef]
  12. He, D.; Wang, E. On the Relation between Soil Water Holding Capacity and Dryland Crop Productivity. Geoderma 2019, 353, 11–24. [Google Scholar] [CrossRef]
  13. Deutscher Wetterdienst (DWD). Klimareport Brandenburg. Fakten Bis Zur Gegenwart—Erwartungen Für Die Zukunft 2020; Deutscher Wetterdienst (DWD): Offenbach am Main, Germany, 2019. [Google Scholar]
  14. Deutscher Wetterdienst (DWD). Mean Precipitation Neuruppin 1991–2020; Deutscher Wetterdienst (DWD): Offenbach am Main, Germany. Available online: https://opendata.dwd.de/climate_environment/CDC/observations_germany and https://opendata.dwd.de/climate_environment/CDC/grids_germany/ (accessed on 8 February 2023).
  15. Deutscher Wetterdienst (DWD). Mean Temperature Neuruppin 1991–2020; Deutscher Wetterdienst (DWD): Offenbach am Main, Germany. Available online: https://opendata.dwd.de/climate_environment/CDC/observations_germany and https://opendata.dwd.de/climate_environment/CDC/grids_germany/ (accessed on 8 February 2023).
  16. Habib-ur-Rahman, M.; Raza, A.; Ahrends, H.E.; Hüging, H.; Gaiser, T. Impact of In-Field Soil Heterogeneity on Biomass and Yield of Winter Triticale in an Intensively Cropped Hummocky Landscape under Temperate Climate Conditions. Precis. Agric. 2022, 23, 912–938. [Google Scholar] [CrossRef]
  17. Kopp, D. Die Böden Des Nordostdeutschen Tieflands Und Ihre Zusammenwirkung Mit Relief, Klima Und Vegetation; BGR, Bundesanstalt für Geowissenschaften und Rohstoffe: Berlin, Germany, 2003. [Google Scholar]
  18. IUSS Working Group WRB. World Reference Base for Soil Resources 2006; World Soil Resources Reports; FAO: Rome, Italy, 2006. [Google Scholar]
  19. Munsell Color (Firm). Munsell Soil Color Charts: With Genuine Munsell Color Chips; 2009 year revised; Munsell Color: Grand Rapids, MI, USA, 2010. [Google Scholar]
  20. Baxter, S. Guidelines for Soil Description. Rome: Food and Agriculture Organization of the United Nations. (2006), pp. 108, US$40.00. ISBN 92-5-1055-21-1. Exp. Agric. 2007, 43, 263–264. [Google Scholar] [CrossRef]
  21. Huang, J.; Hartemink, A.E. Soil and Environmental Issues in Sandy Soils. Earth-Sci. Rev. 2020, 208, 103295. [Google Scholar] [CrossRef]
  22. Planet Labs PBC. Planet Labs Planet Imagery Product Specifications; Planet Labs PBC: San Francisco, CA, USA, 2022. [Google Scholar]
  23. Rouse, J.W.; Haas, R.H.; Schell, J.A.; Deering, D.W. Monitoring Vegetation Systems in the Great Plains with ERTS. NASA Spec. Publ. 1974, 351, 309. [Google Scholar]
  24. Tucker, C.J. Red and Photographic Infrared Linear Combinations for Monitoring Vegetation. Remote Sens. Environ. 1979, 8, 127–150. [Google Scholar] [CrossRef]
  25. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021. [Google Scholar]
  26. Hijmans, R.J. Raster: Geographic Data Analysis and Modeling. 2022. Available online: http://CRAN.R-project.org/package=raster (accessed on 3 May 2023).
  27. Kruskal, W.H.; Wallis, W.A. Use of Ranks in One-Criterion Variance Analysis. J. Am. Stat. Assoc. 1952, 47, 583–621. [Google Scholar] [CrossRef]
  28. Wilcoxon, F. Individual Comparisons by Ranking Methods. Biom. Bull. 1945, 1, 80–83. [Google Scholar] [CrossRef]
  29. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B (Methodol.) 1995, 57, 289–300. [Google Scholar] [CrossRef]
  30. Mann, H.B. Nonparametric Tests Against Trend. Econometrica 1945, 13, 245–259. [Google Scholar] [CrossRef]
  31. McLeod, A.I. Kendall: Kendall Rank Correlation and Mann-Kendall Trend Test. Am. J. Clim. Chang. 2018, 7, 1. [Google Scholar]
  32. Kumar, S.; Machiwal, D.; Dayal, D. Spatial Modelling of Rainfall Trends Using Satellite Datasets and Geographic Information System. Hydrol. Sci. J. 2017, 62, 1636–1653. [Google Scholar] [CrossRef]
  33. Hipel, K.W.; McLeod, A.I. Time Series Modelling of Water Resources and Environmental Systems. In Developments in Water Science, 1st ed.; Elsevier: Amsterdam, The Netherlands, 1974; Volume 45, Available online: https://www.elsevier.com/books/time-series-modelling-of-water-resources-and-environmental-systems/hipel/978-0-444-89270-6 (accessed on 2 April 2023).
  34. Rakovec, O.; Samaniego, L.; Hari, V.; Markonis, Y.; Moravec, V.; Thober, S.; Hanel, M.; Kumar, R. The 2018–2020 Multi-Year Drought Sets a New Benchmark in Europe. Earth’s Future 2022, 10, e2021EF002394. [Google Scholar] [CrossRef]
  35. Boken, V.K.; Shaykewich, C.F. Improving an Operational Wheat Yield Model Using Phenological Phase-Based Normalized Difference Vegetation Index. Int. J. Remote Sens. 2002, 23, 4155–4168. [Google Scholar] [CrossRef]
  36. Wall, L.; Larocque, D.; Léger, P. The Early Explanatory Power of NDVI in Crop Yield Modelling. Int. J. Remote Sens. 2008, 29, 2211–2225. [Google Scholar] [CrossRef]
  37. Mkhabela, M.S.; Bullock, P.; Raj, S.; Wang, S.; Yang, Y. Crop Yield Forecasting on the Canadian Prairies Using MODIS NDVI Data. Agric. For. Meteorol. 2011, 151, 385–393. [Google Scholar] [CrossRef]
  38. Stoy, P.C.; Khan, A.M.; Wipf, A.; Silverman, N.; Powell, S.L. The Spatial Variability of NDVI within a Wheat Field: Information Content and Implications for Yield and Grain Protein Monitoring. PLoS ONE 2022, 17, e0265243. [Google Scholar] [CrossRef]
  39. Fry, J.E.; Guber, A.K. Temporal Stability of Field-Scale Patterns in Soil Water Content across Topographically Diverse Agricultural Landscapes. J. Hydrol. 2020, 580, 124260. [Google Scholar] [CrossRef]
  40. McGoverin, C.M.; Snyders, F.; Muller, N.; Botes, W.; Fox, G.; Manley, M. A Review of Triticale Uses and the Effect of Growth Environment on Grain Quality. J. Sci. Food Agric. 2011, 91, 1155–1165. [Google Scholar] [CrossRef]
  41. Prabhakara, K.; Hively, W.D.; McCarty, G.W. Evaluating the Relationship between Biomass, Percent Groundcover and Remote Sensing Indices across Six Winter Cover Crop Fields in Maryland, United States. Int. J. Appl. Earth Obs. Geoinf. 2015, 39, 88–102. [Google Scholar] [CrossRef]
  42. Diacono, M.; Castrignanò, A.; Troccoli, A.; De Benedetto, D.; Basso, B.; Rubino, P. Spatial and Temporal Variability of Wheat Grain Yield and Quality in a Mediterranean Environment: A Multivariate Geostatistical Approach. Field Crops Res. 2012, 131, 49–62. [Google Scholar] [CrossRef]
  43. Vicente-Serrano, S.M.; Gouveia, C.; Camarero, J.J.; Beguería, S.; Trigo, R.; López-Moreno, J.I.; Azorín-Molina, C.; Pasho, E.; Lorenzo-Lacruz, J.; Revuelto, J.; et al. Response of Vegetation to Drought Time-Scales across Global Land Biomes. Proc. Natl. Acad. Sci. USA 2013, 110, 52–57. [Google Scholar] [CrossRef] [PubMed]
  44. Peled, E.; Dutra, E.; Viterbo, P.; Angert, A. Technical Note: Comparing and Ranking Soil Drought Indices Performance over Europe, through Remote-Sensing of Vegetation. Hydrol. Earth Syst. Sci. 2010, 14, 271–277. [Google Scholar] [CrossRef]
  45. Piedallu, C.; Chéret, V.; Denux, J.P.; Perez, V.; Azcona, J.S.; Seynave, I.; Gégout, J.C. Soil and Climate Differently Impact NDVI Patterns According to the Season and the Stand Type. Sci. Total Environ. 2019, 651, 2874–2885. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Spatial distribution of the three soil and altitude classes. (A): Altitude classes. A1: <53.5 m a.s.l., A2: 53.5–55.37 m a.s.l., A3: >55.37 m a.s.l. (B): Soil classes. S1: Upper boundary of loamy layer at <60 cm depth, S2: Upper boundary of loamy layer at 60–100 cm depth, S3: No loamy layer. EPSG: 32633.
Figure 1. Spatial distribution of the three soil and altitude classes. (A): Altitude classes. A1: <53.5 m a.s.l., A2: 53.5–55.37 m a.s.l., A3: >55.37 m a.s.l. (B): Soil classes. S1: Upper boundary of loamy layer at <60 cm depth, S2: Upper boundary of loamy layer at 60–100 cm depth, S3: No loamy layer. EPSG: 32633.
Agriculture 13 01029 g001
Figure 2. Air temperature (T) [°C], precipitation (P) [mm] and climatic water balance (CWB) [mm] of the experimental site. (A): 2015, (B): 2016, (C): 2017, (D): 2019, (E): 2020.
Figure 2. Air temperature (T) [°C], precipitation (P) [mm] and climatic water balance (CWB) [mm] of the experimental site. (A): 2015, (B): 2016, (C): 2017, (D): 2019, (E): 2020.
Agriculture 13 01029 g002
Figure 3. NDVI values divided into quartiles (1–4). 1 (light yellow) corresponds to quartile 1 (lowest NDVI values) and 4 (dark green) to quartile 4 (highest NDVI values). Triticale: (A): April 2016, (B): May 2016, (C): June 2016. Rye: (D): April 2017, (E): May 2017, (F): June 2017. Quartile classification of NDVI data showed that the northernmost area of the field had consistently lower NDVI values (Q1). In April and May 2017, the NDVI values seemed to be more heterogeneously distributed compared to all other years.
Figure 3. NDVI values divided into quartiles (1–4). 1 (light yellow) corresponds to quartile 1 (lowest NDVI values) and 4 (dark green) to quartile 4 (highest NDVI values). Triticale: (A): April 2016, (B): May 2016, (C): June 2016. Rye: (D): April 2017, (E): May 2017, (F): June 2017. Quartile classification of NDVI data showed that the northernmost area of the field had consistently lower NDVI values (Q1). In April and May 2017, the NDVI values seemed to be more heterogeneously distributed compared to all other years.
Agriculture 13 01029 g003
Figure 4. NDVI values divided into quartiles (1–4) of all June acquisitions. 1 (light yellow) corresponds to quartile 1 (lowest NDVI values) and 4 (dark green) to quartile 4 (highest NDVI values). (A): June 2015, (B): June 2016, (C): June 2017, (D): June 2019, (E): June 2020. Quartile classification of NDVI data showed that the northernmost area of the field had consistently lower NDVI values in June (Q1).
Figure 4. NDVI values divided into quartiles (1–4) of all June acquisitions. 1 (light yellow) corresponds to quartile 1 (lowest NDVI values) and 4 (dark green) to quartile 4 (highest NDVI values). (A): June 2015, (B): June 2016, (C): June 2017, (D): June 2019, (E): June 2020. Quartile classification of NDVI data showed that the northernmost area of the field had consistently lower NDVI values in June (Q1).
Agriculture 13 01029 g004
Table 1. Cultivated crops on the experimental site in chronological order. The growing season here refers to the time span between the sowing date in the first year and the harvest date in the subsequent year.
Table 1. Cultivated crops on the experimental site in chronological order. The growing season here refers to the time span between the sowing date in the first year and the harvest date in the subsequent year.
Growing SeasonSpecieCultivar
2014/15Winter triticale (×Triticosecale)Talentro
2015/16Winter triticale (×Triticosecale)Lombardo
2016/17Winter rye (Secale cereale)KWS Bono
2017/18Summer oat (Avena sativa)Poseidon
2018/19Winter barley (Hordeum vulgare)Anja
2019/20Winter triticale (×Triticosecale)Lombardo
Table 2. Overview of the used satellite images with the satellite identifier and the sensor type. No usable images were available for May 2019.
Table 2. Overview of the used satellite images with the satellite identifier and the sensor type. No usable images were available for May 2019.
DateSatellite Constellation
15 April 2015RapidEye
24 May 2015RapidEye
12 June 2015RapidEye
23 April 2016RapidEye
11 May 2016RapidEye
08 June 2016RapidEye
22 April 2017RapidEye
19 May 2017RapidEye
11 June 2017RapidEye
23 April 2019RapidEye
05 June 2019RapidEye
21 April 2020PlanetScope
21 May 2020PlanetScope
15 June 2020PlanetScope
Table 3. Size [m2] and grid cell number of each soil (S1–S3) and altitude (A1–A3) class with the definition of each class.
Table 3. Size [m2] and grid cell number of each soil (S1–S3) and altitude (A1–A3) class with the definition of each class.
ClassDefinitionNumbers of Grid Cell [n]Area [m2]
S1Loamy layer at <60 cm depth95123.775
S2Loamy layer at 60–100 cm depth3518.775
S3No loamy layer87921.975
A1<53.5 m a.s.l.75218.800
A253.5–55.37 m a.s.l.80420.100
A3>55.37 m a.s.l.62515.625
Table 4. Cross-table to show the overlapping grid cells of the altitude (A1–A3) and soil (S1–S3) classes. In the columns: S1–S3, in the rows: A1–A3. A1: lowest altitude, A3: highest altitude. S1: highest water holding capacity, S3: lowest water holding capacity.
Table 4. Cross-table to show the overlapping grid cells of the altitude (A1–A3) and soil (S1–S3) classes. In the columns: S1–S3, in the rows: A1–A3. A1: lowest altitude, A3: highest altitude. S1: highest water holding capacity, S3: lowest water holding capacity.
ClassS1S2S3
A197100554
A2461171174
A339380151
Table 5. All combinations of not highly significant (p-value ≥ 0.001) differences within the soil and altitude classes calculated with the Wilcoxon signed-rank test.
Table 5. All combinations of not highly significant (p-value ≥ 0.001) differences within the soil and altitude classes calculated with the Wilcoxon signed-rank test.
YearMonthCombinationp-Value
2015AprilS2–S3≥0.05
A2–A31
MayS1–S2≥0.05
A2–A3≥0.05
2016JuneS1–S21
2017AprilS1–S31
A1–A3≥0.05
MayA2–A31
JuneS1–S21
A2–A31
2019AprilS1–S2≥0.05
A1–A2≥0.05
Table 6. Cross table considering all months (April–June). (A): All years, (B): Only years with triticale planted (2015/2016/2020). Each crosstab shows the percentage of grid cells per defined class (1–4) in the columns for one, and in the rows for another, scenery. One hundred percent correspond to all grid cells per row (i.e., quartiles). With this approach, it was possible to quantify the number of spatially equally distributed grid cell classes between two sceneries as well as the change in classes. The results of the cross table showed that the category with the lowest NDVI values (Q1) remained most frequently the same over the years, followed by the class with the highest NDVI value. Consistency in the NDVI category was lowest in April (max. 45%) and highest in June (up to 77%).
Table 6. Cross table considering all months (April–June). (A): All years, (B): Only years with triticale planted (2015/2016/2020). Each crosstab shows the percentage of grid cells per defined class (1–4) in the columns for one, and in the rows for another, scenery. One hundred percent correspond to all grid cells per row (i.e., quartiles). With this approach, it was possible to quantify the number of spatially equally distributed grid cell classes between two sceneries as well as the change in classes. The results of the cross table showed that the category with the lowest NDVI values (Q1) remained most frequently the same over the years, followed by the class with the highest NDVI value. Consistency in the NDVI category was lowest in April (max. 45%) and highest in June (up to 77%).
(A)
Quartile1 [%]2 [%]3 [%]4 [%]
1 to 60.119.311.78.8
2 to1933.227.120.7
3 to11.927.132.828.2
4 to9.120.327.842.7
(B)
Quartile1 [%]2 [%]3 [%]4 [%]
1 to70.218.46.84.6
2 to16.639.726.717
3 to7.727.336.628.4
4 to5.714.729.150.4
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

Mohr, J.; Tewes, A.; Ahrends, H.; Gaiser, T. Assessing the Within-Field Heterogeneity Using Rapid-Eye NDVI Time Series Data. Agriculture 2023, 13, 1029. https://doi.org/10.3390/agriculture13051029

AMA Style

Mohr J, Tewes A, Ahrends H, Gaiser T. Assessing the Within-Field Heterogeneity Using Rapid-Eye NDVI Time Series Data. Agriculture. 2023; 13(5):1029. https://doi.org/10.3390/agriculture13051029

Chicago/Turabian Style

Mohr, Jasper, Andreas Tewes, Hella Ahrends, and Thomas Gaiser. 2023. "Assessing the Within-Field Heterogeneity Using Rapid-Eye NDVI Time Series Data" Agriculture 13, no. 5: 1029. https://doi.org/10.3390/agriculture13051029

APA Style

Mohr, J., Tewes, A., Ahrends, H., & Gaiser, T. (2023). Assessing the Within-Field Heterogeneity Using Rapid-Eye NDVI Time Series Data. Agriculture, 13(5), 1029. https://doi.org/10.3390/agriculture13051029

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