Next Article in Journal
The Influence of Shadow Effects on the Spectral Characteristics of Glacial Meltwater
Next Article in Special Issue
Impacts of Heat and Drought on Gross Primary Productivity in China
Previous Article in Journal
Instance Segmentation for Large, Multi-Channel Remote Sensing Imagery Using Mask-RCNN and a Mosaicking Approach
Previous Article in Special Issue
Estimating Wheat Grain Yield Using Sentinel-2 Imagery and Exploring Topographic Features and Rainfall Effects on Wheat Performance in Navarre, Spain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seasonal Cropland Trends and Their Nexus with Agrometeorological Parameters in the Indus River Plain

Department of Geography, Hong Kong Baptist University, Hong Kong, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(1), 41; https://doi.org/10.3390/rs13010041
Submission received: 24 November 2020 / Revised: 18 December 2020 / Accepted: 21 December 2020 / Published: 24 December 2020
(This article belongs to the Special Issue Remote Sensing of Plant-Climate Interactions)

Abstract

:
The fine-scale insights of existing cropland trends and their nexus with agrometeorological parameters are of paramount importance in assessing future food security risks and analyzing adaptation options under climate change. This study has analyzed the seasonal cropland trends in the Indus River Plain (IRP), using multi-year remote sensing data. A combination of Sen’s slope estimator and Mann–Kendall test was used to quantify the existing cropland trends. A correlation analysis between enhanced vegetation index (EVI) and 9 agrometeorological parameters, derived from reanalysis and remote sensing data, was conducted to study the region’s cropland-climate nexus. The seasonal trend analysis revealed that more than 50% of cropland in IRP improved significantly from the year 2003 to 2018. The lower reaches of the IRP had the highest fraction of cropland, showing a significant decreasing trend during the study period. The nexus analysis showed a strong correlation of EVI with the evaporative stress index (ESI) during the water-stressed crop season. Simultaneously, it exhibited substantial nexus of EVI with actual evapotranspiration (AET) during high soil moisture crop season. Temperature and solar radiation had a negative linkage with EVI response. In contrast, a positive correlation of rainfall with EVI trends was spatially limited to the IRP’s upstream areas. The relative humidity had a spatially broad positive correlation with EVI compare to other direct climatic parameters. The study concluded that positive and sustainable growth in IRP croplands could be achieved through effective agriculture policies to address spatiotemporal AET anomalies.

Graphical Abstract

1. Introduction

The cropland is a complex ecosystem that is subject to many biotic (diseases, pathogens, weeds, etc.) and abiotic (drought, salinity, waterlogging, etc.) stresses, as well as human-induced land-use changes and environmental pollution. Moreover, climate change is directly (by temperature, precipitation, and carbon dioxide changes) or indirectly (by resource degradation and shifts in crop suitability) affecting cropland’s productivity and is becoming an inevitable threat to its sustainability. The negative impacts of climate change on crop productivity can be eased by adopting efficient strategies and making the cropping system more resilient towards extreme events [1]. However, past studies have reported variable impacts of climate parameters on crop productivity subjected to land use patterns and geographic location [2,3,4]. Therefore, local-scale existing trends in cropland and its nexus with climate parameters are of extreme importance to layout efficient adaptation policies to limit climate change adversities on the cropping system. Regional cropland surveillance has been revolutionized since the first Earth observation satellite launch more than five decades ago. The space-borne remote sensing allowed the collection of accurate data to explore the different Earth’s ecosystems on a large scale. The free access to the remote sensing data with moderate spatial resolutions provided the developing countries a way to integrate local-scale ground information into their regional policymaking for better sustainability of the resources.
In recent years, several large-scale studies have utilized remote sensing data to analyze the region’s spatiotemporal vegetation dynamics and its nexus with different climatic parameters. For example, Zhang et al. [5] studied the vegetation dynamics and its driving forces in the upstream region of China’s three main rivers from 1982 to 2012. The net primary productivity (NPP) and its trends were estimated and analyzed using the normalized difference vegetation index (NDVI). The study concluded that NPP in a large part of the area increased significantly, and solar radiation was a key climatic parameter contributing to this change. Similarly, Zhang et al. [6] observed a significant greening trend in global vegetation from 2001 to 2015 and attributed it to enhanced carbon sinkage in the land. Whereas Pan et al. [7], using multi-source remote sensing data, reported an increasing trend of vegetation browning hidden under overall vegetation greening in the northern hemisphere due to drought events. Two studies conducted in the Yangtze River basin concluded that vegetation trends are predominately influenced by human activities [8,9].
The cropland trend analysis may include various components of anthropogenic activities based on the region’s characteristics to better understand the existing patterns and their implication in the future. For example, Wang et al. [10] recently studied the trends of croplands in China’s loess plateau by analyzing the spatiotemporal patterns of field observed evapotranspiration and found a significant upward trend. The factorial analysis of variance between ET and leaf area index (LAI) pointed out that the enhanced cropping intensity was the main driver for the observed cropland trend. Many studies have also analyzed the effect of urban expansion on croplands around the big cities [11,12,13,14]. Most land-use change studies that analyzed the urban expansion on agricultural land have adopted an image classification strategy using multi-temporal images during the study period.
Apart from land-use change studies, the methods used for trend analysis of croplands are primarily based on different parametric and nonparametric statistical tests. The Mann–Kendall trend test is a nonparametric test that has been used in various regional and global studies to detect monotonic trends in vegetative areas [8,15]. The magnitude of trends can be estimated using simple linear regression slope or more robust Theil-Sen’s slope estimator [16,17]. The Pearson correlation coefficient (r) with a significance level can be used to illustrate the relationship between variables and to identify the driving factors of the observed pattern [18,19]. However, the parameters required by the correlation analysis are often limited by their availability from the sources of statistics.
This study aims (1) to identify the existing seasonal trends in cropland of Indus River Plain (IRP) using vegetation index (VI) data from 2003 to 2018; (2) to correlate the temporal response of VI with various agrometeorological parameters to establish existing cropland-climate nexus in the region; (3) to analyze the dominant driving factor and underline causes. The results provide support to agricultural policymaking for improving cropland productivity and regional food security and better understand the potential impact of future climate change.

2. Materials

2.1. Study Area

The Indus river basin lies between latitude 24–37°N and longitude 66–83°E. It encompasses the territories of four countries, including China, Afghanistan, Pakistan, and India. The Indus River Plain (IRP) starts at the foothills of the upper Indus basin. It is a vast expanse of fertile soil covering about 619,000 km2 from Himalayan Piedmont in the north to the Arabian Sea in the south. The IRP can be further divided into upper and lower IRP primarily based on river physiography. Five primary tributaries of the Indus river traverse through the upper IRP and converge into one mainstream before entering to lower IRP. Thal and Thar are two deserts that exist in the IRP (Figure 1).
About 90% of the vegetative area in IRP consists of croplands and forest cover in the north. The fraction of cropland between Pakistan, India, and Afghanistan is 63%, 36%, and 1%, respectively. The cropland in four provinces of Pakistan (Punjab, Sindh, Balochistan, and Khyber Pakhtunkhwa) and three states of India (Punjab, Haryana, and Rajasthan) is the focus of this study. The study area has two major cropping seasons. The Kharif season starts in April-May and lasts till October-November. The Rabi season starts in November-December and ends in April-May of the next year. The study period of this study is from November 2002 to October 2018. The geographical location, administrative units, and significant geographical features of the Indus river basin are shown in Figure 1.

2.2. Remote Sensing Data

The Moderate Resolution Imaging Spectroradiometer (MODIS) is a sensor aboard the Terra and Aqua satellites. There are multiple standard data products of MODIS that can be used for various global to regional studies. Two data products from MODIS were used in this research. One is the enhanced vegetation index (EVI). The MODIS-EVI data from Terra (MOD13A1 (V6)) and Aqua (MYD13A1 (V6)) at 500 m spatial resolution was used to examine cropland trends and to evaluate cropland-climate nexus. By temporally interlacing the 16-day MOD13A1 and MYD13A1, an EVI data set with the 8-day temporal resolution was composed.
The second MODIS-based data product is the actual evapotranspiration (AET). The data product MOD16A2 (V6) with 500 m spatial and 8-day temporal resolution was used to extract AET data. The 8-day MYD16A2 (V6) was also used to fill the data gaps in the primary data product MOD16A2. The MODIS-AET product is based on the Penman-Monteith equation that employs daily meteorological data and vegetation biophysical properties like LAI (Leaf Area Index) and albedo derived from remote sensing [20]. The spatial extent and description of tiles of the MODIS products are shown in Figure 1. The preprocessing has applied to the remote sensing data sets including format conversion, mosaic, subsetting, and reprojection.

2.3. Agrometeorological Data

Due to the sparse spread of meteorological stations in the region and limited data availability, climate parameters from a newly available reanalysis dataset (AgERA5) were used for cropland-climate nexus. The AgERA5 provides gridded climate parameters at a spatial resolution of 10 km from 1979 to 2018. The AgERA5 dataset is based on fifth-generation climate reanalysis data (ERA5) produced by independent intergovernmental organization ECMWF (European Centre for Medium-Range Weather Forecasts). The detailed characteristics and general configuration of ERA5 global reanalysis can be found in Hersbach et al. [21]. Several recent studies have indicated better performance of ERA5-derived climate parameters when compared with other climate reanalysis datasets, satellite-derived, and ground-observed climate parameters [22,23,24,25].
Five climatic parameters, including total daily precipitation (PCP), mean daily temperature (TMP), average daily relative humidity (RH), average daily wind speed (WND), and mean daily solar radiation (SLR) were derived from AgERA5. The Hargreaves method [26] was used to estimate historic potential evapotranspiration (PET) using TMP and SLR of AgERA5. Additionally, two indices were also calculated to indicate water deficiency across the study area. The Aridity Index (AI) was calculated using the ratio of PCP and PET to highlight the natural deficiency of water at a given location, as defined by Barrow [27]. The Evaporative Stress Index (ESI) was derived as ESI = AET/PET and it is structurally similar to the crop coefficient (Kc) defined by Allen et al. [28].

3. Methods

3.1. Overall Methodology

The stepwise procedure of this study is shown in Figure 2. The analytical procedure has two main components: (a) identifying existing croplands trends in the study area, and (b) analyzing the linkage between cropland response to agrometeorological parameters to understand cropland-climate nexus. For the first part, the seasonal stacks of EVI for the Rabi and Kharif seasons were prepared. The Mann–Kendall (MK) test was used to detect EVI’s monotonic trend in both seasons, whereas Sen’s slope estimator was used to quantify the seasonal trend. The results of the MK-test and Sen’s slope were combined to identify different classes of EVI trends.
For the second part, a Pearson correlation analysis between EVI and 9 agrometeorological parameters, including PCP, TMP, SLR, WND, RH, AET, PET, ESI, and AI, was conducted at a monthly level for both seasons. The alternate hypothesis of correlation greater than zero at the significance level of 0.05 was used in estimating the p-value. The spatial analysis of the correlation level between EVI and the agrometeorological parameters was conducted for each season to identify the key parameters responsible for the observed seasonal EVI trend. The spatial resolution of remote sensing data, EVI and AET, was down-scaled to 10 km by the pixel aggregation method, using average value, to match the climatic data spatial resolution.

3.2. Trend Analysis

The Mann–Kendall (MK) is a nonparametric test used to indicate upward or downward monotonic trends shown by a given data set during a period of time. The null hypothesis in the MK test assumes no monotonic trend [29,30]. At a specified significance level (α), the Z-score value rejects the null hypothesis and accepts the alternative hypothesis of the monotonic trend presence. In general, α = 0.05 indicates a statistical significance. The MK test calculates the statistic S as:
S = i = 1 n 1 j = i + 1 n sgn ( x j x i )
where sgn ( x j x i ) is a sign function as
sgn ( x j x i ) = { 1 , x i < x j 0 , x i = x j 1 , x i > x j }
The variance is computed as
V A R ( S ) = 1 18 [ n ( 1 ) ( 2 n + 5 ) p = 1 g t p ( t p 1 ) ( 2 t p + 5 ) ]
where n is the extent of the time series data, x i and x j are the data values in the time series i and j ( j > i ) , respectively, g is the number of tied groups and t p is a tie group number that has equal values in p th group. If there are no tie groups, tie groups’ summation process can be ignored [31]. The Z statistic of the MK test is defined as
Z = { S 1 V A R ( S ) i f S > 0 0 i f S = 0 S + 1 V A R ( S ) i f S < 0 }
The value of Z ranges from ( , + ) . Positive and negative values of Z value indicate increasing and decreasing trends, respectively. At significance level α = 0.05 value of | Z | = ±1.960. For | Z | > 1.960, the null hypothesis is rejected, and a significant trend exists in the time series data.
The Theil–Sen estimator is another robust nonparametric technique used to quantify the slope of the observed trend. Sen’s slope is a robust parameter because the outliers in the data least affect the slope value. Sen’s slope estimator makes a pair-wise combination of total n ( n 1 ) / 2 data pairs and calculates slopes for each of the pairs [32]. The median of all the slopes reflects the change in data per unit time. Mathematically Sen’s slope is expressed as:
β = m e d i a n ( x j x i j i ) , 2002 i < j 2018
where x i and x j are the data values at times i and j , respectively. The sign of β represents the data trend, while its value indicates the change per unit time. In a given period, if β > 0 , the trend of given data is increasing, while β < 0 , a decreasing trend is shown. The combination of Sen’s slope values ( β ) and MK test, Z -score values at the confidence interval of 95% were used to classify the cropland trends in the region. The cropland area with β > 0 and Z > 1.96 (p < 0.05) and area β < 0 and Z < −1.96 (p < 0.05) was classified as the significantly improved and significantly decreased trend, respectively. The area with, β > 0 or β < 0 with Z values between −1.96 to 1.96 (p ≥ 0.05) was termed as having insignificantly improved or decreased cropland trends.

4. Results

4.1. Spatial Variation of Seasonal Cropland Trends

4.1.1. Rabi Season

Based on the classification scheme above, in the study area, 58% of the total cropland during the Rabi season has significantly improved over the past 16 years. In comparison, only 2% of the cropland suffered from significant degradation. A large portion (40%) of the cropland was under insignificant trends during Rabi season. The spatial variation of Z-score, value of Theil-Sen slope, and cropland trend classification during Rabi season is shown in Figure 3. From Figure 3, it is clear that a large portion of cropland in the central region of Punjab (Pak) have shown insignificant decreasing trends. Various clusters of croplands in Sindh province have also demonstrated an insignificant decreasing trend. The croplands at the outskirts of IRP have improved significantly, especially in Thal and Thar regions’ drylands. In Punjab’s arid regions (Pak) and KP, the croplands have shown mostly positive trends with some large portion showing statistically insignificant improvement. More than 70% of the cropland in Indian Punjab and Haryana has shown a significant positive trend. Whereas, 25% of the cropland in Rajasthan showed an insignificant positive trend. The fraction of the cropland area in each category of trend class within a province is shown in Table 1.

4.1.2. Kharif Season

The cropland of IRP during the Kharif season has shown spatially different trends in various parts of regions as compared to the Rabi season. Nearly 50% of the total cropland has shown a significantly increasing trend, and 2% of cropland shows a negative (p < 0.05) trend. About 47% of the croplands are under insignificant trends during this season. The spatial variation of z-score, value of Theil-Sen slope, and cropland trend classes during the Kharif season is shown in Figure 4. The analysis within each province revealed that Balochistan and Sindh have the lowest fraction of the cropland showing a significantly increasing trend. Moreover, both provinces have the highest fraction of croplands with a significantly decreasing trend compared to other provinces/states. The fraction of the cropland area in each category of Kharif season’s trend class at the province/state level is presented in Table 2.

4.2. Nexus of Cropland and Agrometeorological Parameters

4.2.1. Rabi Season

The pair-wise correlation matrix between EVI and 9 agrometeorological parameters is shown in Figure 5. The red color in Figure 5 indicates that Pearson’s correlation (r) value is positive and significant (p ≤ 0.05). The positive and significant nexus of EVI with AET and ESI was among the highest (r = 0.78 and r = 0.70, respectively). From direct climate parameters, EVI correlated with RH and PCP (r = 0.45 and r = 0.18, respectively, at p ≤ 0.05). The nexus of EVI with other agrometeorological parameters and their mutual correlation can be seen in Figure 5. To further analyze the EVI and climate interactions, pixel-level correlations were divided into five categories using r and p. The description of each category is presented in Table 3. Each category’s spatial analysis was done at the basin and sub-basin (state/province) level to identify the primary driver of seasonal cropland trend.
The spatial analysis of the Rabi season’s correlation at basin level revealed that 63% of the IRP croplands have a strong correlation between EVI and ESI. The strong nexus of EVI with AET is at 53% of the total IRP’s cropland. No cropland area demonstrated strong nexus of EVI with PCP and RH. The fractions of cropland area at the basin and sub-basin level with defined nexus classes of EVI and agrometeorological parameters are shown in Figure 6. The nexus class analysis at the province/state level found that the highest cropland fraction of Punjab (Ind, Pak), Haryana, Sindh, and Rajasthan had strong nexus between EVI and ESI compared to AET. In comparison, Balochistan and KP’s cropland had the highest fraction of non-significant nexus of EVI with all agrometeorological parameters. Further details of cropland fractions at defined nexus levels between EVI and various climatic parameters within each province/state can be seen in Figure 6. The pixel-level spatial variation of Rabi season’s cropland-climate nexus is shown in Figure 7. Figure 7 depicted that low altitude croplands of Balochistan and KP had better nexus with ESI and AET. However, these two provinces’ high-altitude croplands were better correlated with TMP, PET, and SLR. The correlation between EVI and TMP was not shown in Figure 7 due to the nexus between TMP and PET.

4.2.2. Kharif Season

For the Kharif season, the pair-wise correlation matrix is shown in Figure 8. Like Rabi season, the EVI is strongly correlated with AET and ESI (r = 0.75 and r = 0.71, respectively, at p ≤ 0.05). Another positive and significant correlation of Kharif season’s EVI is with RH (r = 0.53). The correlation of AI and PCP with EVI is r = 0.42 and r = 0.41, respectively (Figure 8). The spatial analysis of nexus categories at basin level discovered that area fraction under strong nexus class was highest between EVI and RH (37%). The area fraction of total cropland under strong nexus class for AET and ESI was 36% and 28%, respectively. However, the cumulative fraction of strong and relatively strong nexus for Kharif season’ cropland was between EVI and AET (Figure 9).
The fraction of cropland-climate nexus class at basin and sub-basin level for the Kharif season is shown in Figure 9. The spatial analysis at the province/state level revealed that the strong nexus class of EVI with RH was limited to Punjab (Pak and Ind) and Haryana regions. For the rest of the sub-basins, the fraction of strong and relatively strong EVI nexus existed between AET and ESI. Overall, AI compared to PCP had better nexus with Kharif season’s EVI at basin and sub-basin level. The positive nexus of IRP’s seasonal EVI was negligible with agrometeorological parameters such as TMP, SLR, WND, and PET. The spatial variation of the Kharif season’s cropland-climate nexus is shown in Figure 10. The spatial nexus of TMP and PET with EVI was similar; therefore, only EVI interactions with PET are shown in Figure 10.

5. Discussion

The seasonal trend analysis based on EVI indicated a significantly positive trend in more than 50% of the IRP croplands since 2003. This finding of cropland greening is in line with previous regional and global studies. For example, Sarmah et al. [33] analyzed seasonal vegetation trends in South Asia from 2001 to 2013 using VI data from three remote sensing sources. The study reported a significant greening trend in croplands of IRP along with other cropland dominated semi-arid regions. Similarly, Gao et al. [34] analyzed global agricultural trends using remote sensing-derived LAI datasets. The study concluded that agricultural greening trends in developing regions were higher than developed regions, mainly due to cropland expansion and cropping season extension. In the current study, the positive trend of EVI at the peripheries of IRP indicates cropland expansion. Whereas, from the visual analysis of Figure 3 and Figure 4, a negative trend in croplands around the major cities, primarily due to urban sprawl, can be noticed. A large fraction of IRP’s cropland based on time-series MODIS-EVI demonstrated insignificant (increasing and decreasing) trends during both crop seasons. These insignificant trends indicate irregular signals of cropland’s temporal phenology. The prime reason for the signal inconsistency is coarser pixel data of MODIS compared to single farm size and cropland diversity [35].
The correlation analysis of 9 agrometeorological parameters revealed that AET and ESI had throughout strong nexus with EVI in both cropping seasons. The cropland trend and nexus results imply that AET and ESI in IRP’s cropland have also improved significantly during the study period. The AET improvement in regional cropland fall in line with past global-scale research of Javadian et al. [36]. The study identified a significant increase of AET in regional croplands compared to no significant increase in natural vegetation worldwide. Similarly, Anderson et al. [37] studied the correlation of cropland productivity with satellite-driven LAI, PCP, and ESI anomalies. The study concluded that during the flash drought events, ESI was better able to capture soil moisture deficiency and its impact on grain yield. The IRP is a water stress region, especially during the Rabi season when it receives limited rainfall and has minimum flow in rivers. The spatial analysis of Rabi season’s cropland-climate nexus categories at the sub-basin level revealed that ESI strongly correlates with EVI in major crop-producing regions. This means that during the water-limited cropping season, the ESI, which reflects the relative stress in AET, is a better indicator of cropland’s phenology trends.
Among the direct climate parameters, RH had the highest positive nexus with EVI temporal response. The IRP’s areas with significantly positive EVI trends suggest that those regions’ RH level has also increased. This rise in RH is linked to intensive irrigation practices in the Indus river basin, as reported by Krakauer et al. [38]. The spatiotemporal nexus of EVI had a relatively strong linkage with climate dryness and wetness reflected as AI compared to direct PCP rates. The findings of this research demonstrated that AI in the north-eastern side of IRP has improved during the study period, especially during the Kharif season. Kimura [39] studied the satellite-driven global AI trends from 2001 to 2013 and reported a transition towards wetness in various parts of large river basins globally, including irrigated croplands of the Indus basin.
The agrometeorological parameters that include TMP, SLR, WND, and PET had a significantly negative nexus with EVI seasonal trend in the IRP. However, their mutual correlation was significantly positive (Figure 5 and Figure 8). A future increase in any of these agrometeorological parameters resulting from global warming would harm the regional cropland’s EVI response. For example, a rise in PET would increase the crop irrigation requirement to maintain optimal AET rates and ensure stable cropland productivity [40]. The more exploitation of water resources would further degrade the basin’s ecosystem services and aggravate the existing water disputes at local and regional scales. In addition, future rise in TMP would lead to a rapid depletion of snow cover in the upper Indus basin and would limit the primary source of surface irrigation in the IRP [41]. Further investigations of regional snowfall trends are also required to study the linkage of snow melting and sustainability of IRP’s cropland trends, keeping in view the strong relationship between EVI, AET, and ESI established in this study.
To sum up, the EVI trends in IRP is primarily affected by seasonal AET rates and relative stress in AET. Among other farm management practices, irrigation is exceptionally crucial to improve AET rates. The IRP is a heavily irrigated cropland with spatiotemporally skewed rainfall patterns. Regional crop water demands are met by conjunctive irrigation of surface and groundwater [42]. The surface water highly depends on natural conditions such as precipitation and snow melting. Similarly, the improved infrastructure for better water resource management is linked to the region’s economic conditions. The positive EVI trends in the Thal and Thar’s drylands indicate an increase in soil moisture, which is more related to the accessibility to a water supply through canal networks and adaptation of rainwater conservation techniques than an increase in direct rainfall.

6. Conclusions

This study has demonstrated the effective use of satellite remote sensing and fine resolution climate reanalysis data to evaluate the regional cropland trends and identify its agrometeorological drivers. This study concluded that most parts of IRP’s cropland have significantly improved during the study period based on MODIS-EVI trends. The nexus analysis revealed that ESI had strong nexus with EVI during water limiting crop season, whereas for water excess crop season, AET was strongly correlated with EVI response. This research reiterates the importance of evapotranspiration and its stresses in controlling the cropland trends, especially in intensively irrigated regions with an arid and semi-arid climate. Knowledge of the existing cropland-climate nexus can improve the interpretation of future climate change impacts on regional cropland. This approach can be applied to other data-scarce regions to plan robust measures against changing climate and ensure sustainable productivity of regional cropland.

Author Contributions

Conceptualization, Q.Z.; Data curation, A.I.; Formal analysis, A.I.; Funding acquisition, Q.Z.; Methodology, A.I.; Supervision, Q.Z.; Writing—original draft, A.I.; Writing—review & editing, Q.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the Natural Science Foundation of China (NSFC) General Research Grant (41971386) Hong Kong Research Grant Council (RGC) General Research Fund (12301820). The work is a part of Ph.D. research funded by Hong Kong Ph.D. Fellowship Scheme (HKPFS).

Institutional Review Board Statement

Not applicable for this study.

Informed Consent Statement

Not applicable for this study.

Data Availability Statement

The input data used in this research can be accessed freely from online sources. Remote sensing data used in this study can be downloaded via various web portals of NASA, e.g., https://earthdata.nasa.gov/. The climate reanalysis is available at climate data store website of European Centre for Medium-Range Weather Forecasts (https://cds.climate.copernicus.eu/). All the derived data is included in the manuscript. The trend and correlation analysis was done in MATLAB interface and scripts can be shared upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Raza, A.; Razzaq, A.; Mehmood, S.S.; Zou, X.; Zhang, X.; Lv, Y.; Xu, J. Impact of climate change on crops adaptation and strategies to tackle its outcome: A review. Plants 2019, 8, 34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ray, D.K.; Gerber, J.S.; Macdonald, G.K.; West, P.C. Climate variation explains a third of global crop yield variability. Nat. Commun. 2015, 6, 5989. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Lobell, D.B.; Schlenker, W.; Costa-Roberts, J. Climate trends and global crop production since 1980. Science 2011, 333, 616–620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Kukal, M.S.; Irmak, S. Climate-Driven Crop Yield and Yield Variability and Climate Change Impacts on the U.S. Great Plains Agricultural Production. Sci. Rep. 2018, 8, 1–18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Zhang, Y.; Zhang, C.; Wang, Z.; Chen, Y.; Gang, C.; An, R.; Li, J. Vegetation dynamics and its driving forces from climate change and human activities in the Three-River Source Region, China from 1982 to 2012. Sci. Total Environ. 2016, 563, 210–220. [Google Scholar] [CrossRef]
  6. Zhang, Y.; Song, C.; Band, L.E.; Sun, G.; Li, J. Reanalysis of global terrestrial vegetation trends from MODIS products: Browning or greening? Remote Sens. Environ. 2017, 191, 145–155. [Google Scholar] [CrossRef] [Green Version]
  7. Pan, N.; Feng, X.; Fu, B.; Wang, S.; Ji, F.; Pan, S. Increasing global vegetation browning hidden in overall vegetation greening: Insights from time-varying trends. Remote Sens. Environ. 2018, 214, 59–72. [Google Scholar] [CrossRef]
  8. Qu, S.; Wang, L.; Lin, A.; Zhu, H.; Yuan, M. What drives the vegetation restoration in Yangtze River basin, China: Climate change or anthropogenic factors? Ecol. Indic. 2018, 90, 438–450. [Google Scholar] [CrossRef]
  9. Qu, S.; Wang, L.; Lin, A.; Yu, D.; Yuan, M.; Li, C. Distinguishing the impacts of climate change and anthropogenic factors on vegetation dynamics in the Yangtze River Basin, China. Ecol. Indic. 2020, 108, 105724. [Google Scholar] [CrossRef]
  10. Wang, F.; Liang, W.; Fu, B.; Jin, Z.; Yan, J.; Zhang, W.; Fu, S.; Yan, N. Changes of cropland evapotranspiration and its driving factors on the loess plateau of China. Sci. Total Environ. 2020, 728, 138582. [Google Scholar] [CrossRef]
  11. Tang, J.; Di, L. Past and future trajectories of farmland loss due to rapid urbanization using Landsat imagery and the Markov-CA model: A case study of Delhi, India. Remote Sens. 2019, 11, 180. [Google Scholar] [CrossRef] [Green Version]
  12. Gumma, M.K.; Mohammad, I.; Nedumaran, S.; Whitbread, A.; Lagerkvist, C.J. Urban sprawl and adverse impacts on agricultural land: A case study on Hyderabad, India. Remote Sens. 2017, 9, 1136. [Google Scholar] [CrossRef] [Green Version]
  13. Xu, Y.; Yu, L.; Zhao, F.R.; Cai, X.; Zhao, J.; Lu, H.; Gong, P. Tracking annual cropland changes from 1984 to 2016 using time-series Landsat images with a change-detection and post-classification approach: Experiments from three sites in Africa. Remote Sens. Environ. 2018, 218, 13–31. [Google Scholar] [CrossRef]
  14. Liu, F.; Zhang, Z.; Zhao, X.; Wang, X.; Zuo, L.; Wen, Q.; Yi, L.; Xu, J.; Hu, S.; Liu, B. Chinese cropland losses due to urban expansion in the past four decades. Sci. Total Environ. 2019, 650, 847–857. [Google Scholar] [CrossRef]
  15. Fassnacht, F.E.; Schiller, C.; Kattenborn, T.; Zhao, X.; Qu, J. A Landsat-based vegetation trend product of the Tibetan Plateau for the time-period 1990–2018. Sci. Data 2019, 6, 1–11. [Google Scholar] [CrossRef]
  16. Liu, Y.; Li, L.; Chen, X.; Zhang, R.; Yang, J. Temporal-spatial variations and influencing factors of vegetation cover in Xinjiang from 1982 to 2013 based on GIMMS-NDVI3g. Glob. Planet. Chang. 2018, 169, 145–155. [Google Scholar] [CrossRef]
  17. Lamchin, M.; Lee, W.K.; Jeon, S.W.; Wang, S.W.; Lim, C.H.; Song, C.; Sung, M. Long-term trend and correlation between vegetation greenness and climate variables in Asia based on satellite data. Sci. Total Environ. 2018, 618, 1089–1095. [Google Scholar] [CrossRef]
  18. Zhao, Q.; Ma, X.; Liang, L.; Yao, W. Spatial-temporal variation characteristics of multiple meteorological variables and vegetation over the loess plateau region. Appl. Sci. 2020, 10, 1000. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, Y.; Wang, X.; Li, C.; Cai, Y.; Yang, Z.; Yi, Y. NDVI dynamics under changing meteorological factors in a shallow lake in future metropolitan, semiarid area in North China. Sci. Rep. 2018, 8, 1–13. [Google Scholar] [CrossRef] [Green Version]
  20. Running, S.; Mu, Q.; Zhao, M.; Moreno, A. MODIS Global Terrestrial Evapotranspiration (ET) Product (NASA MOD16A2/A3) NASA Earth Observing System MODIS Land Algorithm; NASA: Washington, DC, USA, 2017.
  21. Hersbach, H.; Bell, B.; Berrisford, P.; Hirahara, S.; Horányi, A.; Muñoz-Sabater, J.; Nicolas, J.; Peubey, C.; Radu, R.; Schepers, D.; et al. The ERA5 global reanalysis. Q. J. R. Meteorol. Soc. 2020, 146, 1999–2049. [Google Scholar] [CrossRef]
  22. Urraca, R.; Huld, T.; Gracia-Amillo, A.; Martinez-de-Pison, F.J.; Kaspar, F.; Sanz-Garcia, A. Evaluation of global horizontal irradiance estimates from ERA5 and COSMO-REA6 reanalyses using ground and satellite-based data. Sol. Energy 2018, 164, 339–354. [Google Scholar] [CrossRef]
  23. Nogueira, M. Inter-comparison of ERA-5, ERA-interim and GPCP rainfall over the last 40 years: Process-based analysis of systematic and random differences. J. Hydrol. 2020, 583, 124632. [Google Scholar] [CrossRef]
  24. Tarek, M.; Brissette, F.P.; Arsenault, R. Evaluation of the ERA5 reanalysis as a potential reference dataset for hydrological modelling over North America. Hydrol. Earth Syst. Sci. 2020, 24, 2527–2544. [Google Scholar] [CrossRef]
  25. Belmonte Rivas, M.; Stoffelen, A. Characterizing ERA-Interim and ERA5 surface wind biases using ASCAT. Ocean Sci. 2019, 15, 831–852. [Google Scholar] [CrossRef] [Green Version]
  26. Hargreaves, G.H. Moisture availability and crop production. Trans. ASAE 1975, 18, 0980–0984. [Google Scholar] [CrossRef]
  27. Barrow, C.J. World atlas of desertification (United nations environment programme), edited by N. Middleton and D. S. G. Thomas. Edward Arnold, London. Land Degrad. Dev. 1992, 3, 249. [Google Scholar] [CrossRef]
  28. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration—Guidelines for Computing Crop Water Requirements—FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998. [Google Scholar]
  29. Kendall, M.G. Rank Correlation Methods. In Biometrika; Oxford University Press: Oxford, UK, 1957; p. 298. [Google Scholar]
  30. Mann, H.B. Nonparametric Tests Against Trend. Econometrica 1945, 13, 245–259. [Google Scholar] [CrossRef]
  31. Kisi, O.; Ay, M. Comparison of Mann-Kendall and innovative trend method for water quality parameters of the Kizilirmak River, Turkey. J. Hydrol. 2014, 513, 362–375. [Google Scholar] [CrossRef]
  32. Sen, P.K. Estimates of the Regression Coefficient Based on Kendall’s Tau. J. Am. Stat. Assoc. 1968, 63, 1379–1389. [Google Scholar] [CrossRef]
  33. Sarmah, S.; Jia, G.; Zhang, A.; Singha, M. Assessing seasonal trends and variability of vegetation growth from NDVI3G, MODIS NDVI and EVI over South Asia. Remote Sens. Lett. 2018, 9, 1195–1204. [Google Scholar] [CrossRef]
  34. Gao, X.; Liang, S.; He, B. Detected global agricultural greening from satellite data. Agric. For. Meteorol. 2019, 276, 107652. [Google Scholar] [CrossRef]
  35. Nabil, M.; Zhang, M.; Bofana, J.; Wu, B.; Stein, A.; Dong, T.; Zeng, H.; Shang, J. Assessing factors impacting the spatial discrepancy of remote sensing based cropland products: A case study in Africa. Int. J. Appl. Earth Obs. Geoinf. 2020, 85, 102010. [Google Scholar] [CrossRef]
  36. Javadian, M.; Behrangi, A.; Smith, W.K.; Fisher, J.B. Global trends in evapotranspiration dominated by increases across large cropland regions. Remote Sens. 2020, 12, 1221. [Google Scholar] [CrossRef] [Green Version]
  37. Anderson, M.C.; Zolin, C.A.; Sentelhas, P.C.; Hain, C.R.; Semmens, K.; Tugrul Yilmaz, M.; Gao, F.; Otkin, J.A.; Tetrault, R. The Evaporative Stress Index as an indicator of agricultural drought in Brazil: An assessment based on crop yield impacts. Remote Sens. Environ. 2016, 174, 82–99. [Google Scholar] [CrossRef]
  38. Krakauer, N.Y.; Cook, B.I.; Puma, M.J. Effect of irrigation on humid heat extremes. Environ. Res. Lett. 2020, 15, 094010. [Google Scholar] [CrossRef]
  39. Kimura, R. Global detection of aridification or increasing wetness in arid regions from 2001 to 2013. Nat. Hazards 2020, 103, 2261–2276. [Google Scholar] [CrossRef]
  40. Dimri, A.P.; Kumar, D.; Chopra, S.; Choudhary, A. Indus River Basin: Future climate and water budget. Int. J. Climatol. 2019, 39, 395–406. [Google Scholar] [CrossRef] [Green Version]
  41. Adnan, M.; Nabi, G.; Saleem Poomee, M.; Ashraf, A. Snowmelt runoff prediction under changing climate in the Himalayan cryosphere: A case of Gilgit River Basin. Geosci. Front. 2017, 8, 941–949. [Google Scholar] [CrossRef] [Green Version]
  42. Usman, M.; Abbas, A.; Saqib, Z.A. Conjunctive use of water and its management for enhanced productivity of major crops across tertiary canal irrigation system of indus basin in Pakistan. Pak. J. Agric. Sci. 2016, 53, 257–264. [Google Scholar] [CrossRef]
Figure 1. The geographical location of the Indus river basin, its plains, and administrative boundaries, MODIS tiles coverage, and major geographical features.
Figure 1. The geographical location of the Indus river basin, its plains, and administrative boundaries, MODIS tiles coverage, and major geographical features.
Remotesensing 13 00041 g001
Figure 2. The stepwise procedure involved in completing this research.
Figure 2. The stepwise procedure involved in completing this research.
Remotesensing 13 00041 g002
Figure 3. The spatial variation of z-score, the value of Theil-Sen slope, and cropland trend classification during Rabi season.
Figure 3. The spatial variation of z-score, the value of Theil-Sen slope, and cropland trend classification during Rabi season.
Remotesensing 13 00041 g003
Figure 4. The spatial variation of z-score, value of Theil-Sen slope, and cropland trend classes during the Kharif season.
Figure 4. The spatial variation of z-score, value of Theil-Sen slope, and cropland trend classes during the Kharif season.
Remotesensing 13 00041 g004
Figure 5. The pair-wise correlation matrix between enhanced vegetation index (EVI) and nine agrometeorological parameters during Rabi season.
Figure 5. The pair-wise correlation matrix between enhanced vegetation index (EVI) and nine agrometeorological parameters during Rabi season.
Remotesensing 13 00041 g005
Figure 6. The fractions of cropland area at basin and sub-basin level with defined nexus classes. The class fraction is the Rabi season’s correlation between EVI and 9 agrometeorological parameters as indicated on the x-axis.
Figure 6. The fractions of cropland area at basin and sub-basin level with defined nexus classes. The class fraction is the Rabi season’s correlation between EVI and 9 agrometeorological parameters as indicated on the x-axis.
Remotesensing 13 00041 g006
Figure 7. The pixel-level spatial variation of Rabi season’s cropland-climate nexus.
Figure 7. The pixel-level spatial variation of Rabi season’s cropland-climate nexus.
Remotesensing 13 00041 g007
Figure 8. The pair-wise correlation matrix between EVI and nine agrometeorological parameters during Kharif season.
Figure 8. The pair-wise correlation matrix between EVI and nine agrometeorological parameters during Kharif season.
Remotesensing 13 00041 g008
Figure 9. The fraction of cropland-climate nexus class at basin and sub-basin level for Kharif season. The class fraction is the seasonal correlation between EVI and 9 agrometeorological parameters as indicated on the x-axis.
Figure 9. The fraction of cropland-climate nexus class at basin and sub-basin level for Kharif season. The class fraction is the seasonal correlation between EVI and 9 agrometeorological parameters as indicated on the x-axis.
Remotesensing 13 00041 g009
Figure 10. The pixel-level spatial variation of Kharif season’s cropland-climate nexus.
Figure 10. The pixel-level spatial variation of Kharif season’s cropland-climate nexus.
Remotesensing 13 00041 g010
Table 1. The fraction of provincial cropland area in each trend category during Rabi season.
Table 1. The fraction of provincial cropland area in each trend category during Rabi season.
Sub-Basin/Trend ClassIncrease (p < 0.05)Decrease (p < 0.05)Increase (p > 0.05)Decrease (p > 0.05)
Balochistan56.750.798.3634.10
KP48.241.7212.4537.59
Punjab (Pak)55.542.0916.7925.58
Sindh39.876.6723.0130.44
Haryana76.820.906.8015.48
Punjab (Ind)71.030.9813.1514.85
Rajasthan67.281.076.4325.21
Table 2. The fraction of provincial cropland area in each trend category during Kharif season.
Table 2. The fraction of provincial cropland area in each trend category during Kharif season.
Sub-Basin/Trend ClassIncrease (p < 0.05)Decrease (p < 0.05)Increase (p > 0.05)Decrease (p > 0.05)
Balochistan26.955.3344.2923.44
KP47.351.0442.309.31
Punjab (Pak)52.461.8834.6710.99
Sindh27.046.2636.2030.50
Haryana53.540.9829.6515.84
Punjab (Ind)68.151.1115.4615.28
Rajasthan61.890.1435.962.01
Table 3. The description of cropland-climate nexus categories used in the spatial analysis.
Table 3. The description of cropland-climate nexus categories used in the spatial analysis.
No.Nexus LevelR-Valuep-Value
1Strong>0.8 to 1≤0.05
2Relatively Strong>0.6 to 0.8≤0.05
3Moderate>0.4 to 0.6≤0.05
4Weak>0 to 0.4≤0.05
5No Nexus<0>0.05
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, Q.; Ismaeel, A. Seasonal Cropland Trends and Their Nexus with Agrometeorological Parameters in the Indus River Plain. Remote Sens. 2021, 13, 41. https://doi.org/10.3390/rs13010041

AMA Style

Zhou Q, Ismaeel A. Seasonal Cropland Trends and Their Nexus with Agrometeorological Parameters in the Indus River Plain. Remote Sensing. 2021; 13(1):41. https://doi.org/10.3390/rs13010041

Chicago/Turabian Style

Zhou, Qiming, and Ali Ismaeel. 2021. "Seasonal Cropland Trends and Their Nexus with Agrometeorological Parameters in the Indus River Plain" Remote Sensing 13, no. 1: 41. https://doi.org/10.3390/rs13010041

APA Style

Zhou, Q., & Ismaeel, A. (2021). Seasonal Cropland Trends and Their Nexus with Agrometeorological Parameters in the Indus River Plain. Remote Sensing, 13(1), 41. https://doi.org/10.3390/rs13010041

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