农业生态环境-灌溉Agro-ecosystem & Environment—Irrigation
Monitoring agricultural drought using remote sensing data is crucial for precision irrigation in modern agriculture. Utilizing unmanned aerial vehicle (UAV) remote sensing, we explored the applicability of an empirical crop water stress index (CWSI) based on canopy temperature and three-dimensional drought indices (TDDI) constructed from surface temperature (Ts), air temperature (Ta) and five vegetation indices (VIs) for monitoring the moisture status of dryland crops. Three machine learning algorithms (random forest regression (RFR), support vector regression, and partial least squares regression) were used to compare the performance of the drought indices for vegetation moisture content (VMC) estimation in sorghum and maize. The main results of the study were as follows: (1) Comparative analysis of the drought indices revealed that Ts-Ta-normalized difference vegetation index (TDDIn) and Ts-Ta-enhanced vegetation index (TDDIe) were more strongly correlated with VMC compared with the other indices. The indices exhibited varying sensitivities to VMC under different irrigation regimes; the strongest correlation observed was for the TDDIe index with maize under the fully irrigated treatment (r=−0.93). (2) Regarding spatial and temporal characteristics, the TDDIn, TDDIe and CWSI indices showed minimal differences. Over the experimental period, with coefficients of variation were 0.25, 0.18 and 0.24, respectively. All three indices were capable of effectively characterizing the moisture distribution in dryland maize and sorghum crops, but the TDDI indices more accurately monitored the spatial distribution of crop moisture after a rainfall or irrigation event. (3) For prediction of the moisture content of single crops, RFR models based on TDDIn and TDDIe estimated VMC most accurately (R2>0.7), and the TDDIn-based model predicted VMC with the highest accuracy when considering multiple-crop samples, with R2 and RMSE of 0.62 and 14.26%, respectively. Thus, TDDI proved more effective than the CWSI in estimating crop water content.
Rapid and accurate acquisition of soil organic matter (SOM) information in cultivated land is important for sustainable agricultural development and carbon balance management. This study proposed a novel approach to predict SOM with high accuracy using multiyear synthetic remote sensing variables on a monthly scale. We obtained 12 monthly synthetic Sentinel-2 images covering the study area from 2016 to 2021 through the Google Earth Engine (GEE) platform, and reflectance bands and vegetation indices were extracted from these composite images. Then the random forest (RF), support vector machine (SVM) and gradient boosting regression tree (GBRT) models were tested to investigate the difference in SOM prediction accuracy under different combinations of monthly synthetic variables. Results showed that firstly, all monthly synthetic spectral bands of Sentinel-2 showed a significant correlation with SOM (P<0.05) for the months of January, March, April, October, and November. Secondly, in terms of single-monthly composite variables, the prediction accuracy was relatively poor, with the highest R2 value of 0.36 being observed in January. When monthly synthetic environmental variables were grouped in accordance with the four quarters of the year, the first quarter and the fourth quarter showed good performance, and any combination of three quarters was similar in estimation accuracy. The overall best performance was observed when all monthly synthetic variables were incorporated into the models. Thirdly, among the three models compared, the RF model was consistently more accurate than the SVM and GBRT models, achieving an R2 value of 0.56. Except for band 12 in December, the importance of the remaining bands did not exhibit significant differences. This research offers a new attempt to map SOM with high accuracy and fine spatial resolution based on monthly synthetic Sentinel-2 images.
Improved simulation of winter wheat yield in North China Plain by using PRYM-Wheat integrated dry matter distribution coefficient
The accurate simulation of regional-scale winter wheat yield is important for national food security and the balance of grain supply and demand in China. Presently, most remote sensing process models use the “biomass×harvest index (HI)” method to simulate regional-scale winter wheat yield. However, spatiotemporal differences in HI contribute to inaccuracies in yield simulation at the regional scale. Time-series dry matter partition coefficients (Fr) can dynamically reflect the dry matter partition of winter wheat. In this study, Fr equations were fitted for each organ of winter wheat using site-scale data. These equations were then coupled into a process-based and remote sensing-driven crop yield model for wheat (PRYM-Wheat) to improve the regional simulation of winter wheat yield over the North China Plain (NCP). The improved PRYM-Wheat model integrated with the fitted Fr equations (PRYM-Wheat-Fr) was validated using data obtained from provincial yearbooks. A 3-year (2000–2002) averaged validation showed that PRYM-Wheat-Fr had a higher coefficient of determination (R²=0.55) and lower root mean square error (RMSE=0.94 t ha–1) than PRYM-Wheat with a stable HI (abbreviated as PRYM-Wheat-HI), which had R² and RMSE values of 0.30 and 1.62 t ha–1, respectively. The PRYM-Wheat-Fr model also performed better than PRYM-Wheat-HI for simulating yield in verification years (2013–2015). In conclusion, the PRYM-Wheat-Fr model exhibited a better accuracy than the original PRYM-Wheat model, making it a useful tool for the simulation of regional winter wheat yield.
Combining field data and modeling to better understand maize growth response to phosphorus (P) fertilizer application and soil P dynamics in calcareous soils
A hybrid CNN-LSTM model for diagnosing rice nutrient levels at the rice panicle initiation stage
Nitrogen (N) and potassium (K) are two key mineral nutrient elements involved in rice growth. Accurate diagnosis of N and K status is very important for the rational application of fertilizers at a specific rice growth stage. Therefore, we propose a hybrid model for diagnosing rice nutrient levels at the early panicle initiation stage (EPIS), which combines a convolutional neural network (CNN) with an attention mechanism and a long short-term memory network (LSTM). The model was validated on a large set of sequential images collected by an unmanned aerial vehicle (UAV) from rice canopies at different growth stages during a two-year experiment. Compared with VGG16, AlexNet, GoogleNet, DenseNet, and inceptionV3, ResNet101 combined with LSTM obtained the highest average accuracy of 83.81% on the dataset of Huanghuazhan (HHZ, an indica cultivar). When tested on the datasets of HHZ and Xiushui 134 (XS134, a japonica rice variety) in 2021, the ResNet101-LSTM model enhanced with the squeeze-and-excitation (SE) block achieved the highest accuracies of 85.38 and 88.38%, respectively. Through the cross-dataset method, the average accuracies on the HHZ and XS134 datasets tested in 2022 were 81.25 and 82.50%, respectively, showing a good generalization. Our proposed model works with the dynamic information of different rice growth stages and can efficiently diagnose different rice nutrient status levels at EPIS, which are helpful for making practical decisions regarding rational fertilization treatments at the panicle initiation stage.
Quantifying the agreement and accuracy characteristics of four satellite-based LULC products for cropland classification in China
Various land use and land cover (LULC) products have been produced over the past decade with the development of remote sensing technology. Despite the differences in LULC classification schemes, there is a lack of research on assessing the accuracy of their application to croplands in a unified framework. Thus, this study evaluated the spatial and area accuracies of cropland classification for four commonly used global LULC products (i.e., MCD12Q1 V6, GlobCover2009, FROM-GLC and GlobeLand30) based on the harmonised FAO criterion, and quantified the relationships between four factors (i.e., slope, elevation, field size and crop system) and cropland classification agreement. The validation results indicated that MCD12Q1 and GlobeLand30 performed well in cropland classification regarding spatial consistency, with overall accuracies of 94.90 and 93.52%, respectively. The FROM-GLC showed the worst performance, with an overall accuracy of 83.17%. Overlaying the cropland generated by the four global LULC products, we found the proportions of complete agreement and disagreement were 15.51 and 44.72% for the cropland classification, respectively. High consistency was mainly observed in the Northeast China Plain, the Huang-Huai-Hai Plain and the northern part of the Middle-lower Yangtze Plain, China. In contrast, low consistency was detected primarily on the eastern edge of the northern and semiarid region, the Yunnan-Guizhou Plateau and southern China. Field size was the most important factor for mapping cropland. For area accuracy, compared with China Statistical Yearbook data at the provincial scale, the accuracies of different products in descending order were: GlobeLand30, FROM-GLC, MCD12Q1, and GlobCover2009. The cropland classification schemes mainly caused large area deviations among the four products, and they also resulted in the different ranks of spatial accuracy and area accuracy among the four products. Our results can provide valuable suggestions for selecting cropland products at the national or provincial scale and help cropland mapping and reconstruction, which is essential for food security and crop management, so they can also contribute to achieving the Sustainable Development Goals issued by the United Nations.
Studying the significant impacts on vegetation of drought due to global warming is crucial in order to understand its dynamics and interrelationships with temperature, rainfall, and normalized difference vegetation index (NDVI). These factors are linked to excesses drought frequency and severity on the regional scale, and their effect on vegetation remains an important topic for climate change study. East Asia is very sensitive and susceptible to climate change. In this study, we examined the effect of drought on the seasonal variations of vegetation in relation to climate variability and determined which growing seasons are most vulnerable to drought risk; and then explored the spatio-temporal evolution of the trend in drought changes in East Asia from 1982 to 2019. The data were studied using a series of several drought indexes, and the data were then classified using a heat map, box and whisker plot analysis, and principal component analysis. The various drought indexes from January to August improved rapidly, except for vegetation health index (VHI) and temperature condition index (TCI). While these indices were constant in September, they increased again in October, but in December, they showed a descending trend. The seasonal and monthly analysis of the drought indexes and the heat map confirmed that the East Asian region suffered from extreme droughts in 1984, 1993, 2007, and 2012 among the study years. The distribution of the trend in drought changes indicated that more severe drought occurred in the northwestern region than in the southeastern area of East Asia. The drought tendency slope was used to describe the changes in drought events during 1982–2019 in the study region. The correlations among monthly precipitation anomaly percentage (NAP), NDVI, TCI, vegetation condition index (VCI), temperature vegetation drought index (TVDI), and VHI indicated considerably positive correlations, while considerably negative correlations were found among the three pairs of NDVI and VHI, TVDI and VHI, and NDVI and TCI. This ecological and climatic mechanism provides a good basis for the assessment of vegetation and drought-change variations within the East Asian region. This study is a step forward in monitoring the seasonal variation of vegetation and variations in drought dynamics within the East Asian region, which will serve and contribute to the better management of vegetation, disaster risk, and drought in the East Asian region.
Accurate estimation of regional winter wheat yields is essential for understanding the food production status and ensuring national food security. However, using the existing remote sensing-based crop yield models to accurately reproduce the inter-annual and spatial variations in winter wheat yields remains challenging due to the limited ability to acquire irrigation information in water-limited regions. Thus, we proposed a new approach to approximating irrigations of winter wheat over the North China Plain (NCP), where irrigation occurs extensively during the winter wheat growing season. This approach used irrigation pattern parameters (IPPs) to define the irrigation frequency and timing. Then, they were incorporated into a newly-developed process-based and remote sensing-driven crop yield model for winter wheat (PRYM–Wheat), to improve the regional estimates of winter wheat over the NCP. The IPPs were determined using statistical yield data of reference years (2010–2015) over the NCP. Our findings showed that PRYM–Wheat with the optimal IPPs could improve the regional estimate of winter wheat yield, with an increase and decrease in the correlation coefficient (R) and root mean square error (RMSE) of 0.15 (about 37%) and 0.90 t ha–1 (about 41%), respectively. The data in validation years (2001–2009 and 2016–2019) were used to validate PRYM–Wheat. In addition, our findings also showed R (RMSE) of 0.80 (0.62 t ha–1) on a site level, 0.61 (0.91 t ha–1) for Hebei Province on a county level, 0.73 (0.97 t ha–1) for Henan Province on a county level, and 0.55 (0.75 t ha–1) for Shandong Province on a city level. Overall, PRYM–Wheat can offer a stable and robust approach to estimating regional winter wheat yield across multiple years, providing a scientific basis for ensuring regional food security.
Remote sensing has been increasingly used for precision nitrogen management to assess the plant nitrogen status in a spatial and real-time manner. The nitrogen nutrition index (NNI) can quantitatively describe the nitrogen status of crops. Nevertheless, the NNI diagnosis for cotton with unmanned aerial vehicle (UAV) multispectral images has not been evaluated yet. This study aimed to evaluate the performance of three machine learning models, i.e., support vector machine (SVM), back propagation neural network (BPNN), and extreme gradient boosting (XGB) for predicting canopy nitrogen weight and NNI of cotton over the whole growing season from UAV images. The results indicated that the models performed better when the top 15 vegetation indices were used as input variables based on their correlation ranking with nitrogen weight and NNI. The XGB model performed the best among the three models in predicting nitrogen weight. The prediction accuracy of nitrogen weight at the upper half-leaf level (R2=0.89, RMSE=0.68 g m–2, RE=14.62% for calibration and R2=0.83, RMSE=1.08 g m–2, RE=19.71% for validation) was much better than that at the all-leaf level (R2=0.73, RMSE=2.20 g m–2, RE=26.70% for calibration and R2=0.70, RMSE=2.48 g m–2, RE=31.49% for validation) and at the plant level (R2=0.66, RMSE=4.46 g m–2, RE=30.96% for calibration and R2=0.63, RMSE=3.69 g m–2, RE=24.81% for validation). Similarly, the XGB model (R2=0.65, RMSE=0.09, RE=8.59% for calibration and R2=0.63, RMSE=0.09, RE=8.87% for validation) also outperformed the SVM model (R2=0.62, RMSE=0.10, RE=7.92% for calibration and R2=0.60, RMSE=0.09, RE=8.03% for validation) and BPNN model (R2=0.64, RMSE=0.09, RE=9.24% for calibration and R2=0.62, RMSE=0.09, RE=8.38% for validation) in predicting NNI. The NNI predictive map generated from the optimal XGB model can intuitively diagnose the spatial distribution and dynamics of nitrogen nutrition in cotton fields, which can help farmers implement precise cotton nitrogen management in a timely and accurate manner
Rapid and large area acquisition of nitrogen (N) deficiency status is important for achieving the optimal fertilization of rice. Most existing studies, however, focus on the use of unmanned aerial vehicle (UAV) remote sensing to diagnose N nutrition in rice, while there are fewer studies on the quantitative description of the degree of N deficiency in rice, and the effects of the critical N concentration on the spectral changes in rice have rarely been explored. Therefore, based on the canopy spectral data obtained by remotely-sensed UAV hyperspectral images, the N content in rice was obtained through field sampling. The construction method of the rice curve for the northeastern critical N concentration was studied, and on this basis, N deficiency was determined. Taking the spectrum of the critical N concentration state as the standard spectrum, the spectral reflectivity data were transformed by the ratios and differences, and the feature extraction of the spectral data was carried out by the successive projections algorithm (SPA). Finally, by taking the characteristic band as the input variable and N deficiency as the output variable, a set of multivariate linear regression (MLR), long short-term memory (LSTM) inversion models based on extreme learning machine (ELM), and the non-dominated sorting genetic algorithm III extreme learning machine (NSGA-III-ELM) were constructed. The results showed two key aspects of this system: 1) The correlation between the N deficiency data and original spectrum was poor, but the correlation between the N deficiency data and N deficiency could be improved by a difference change and ratio transformation; 2) The inversion results based on the ratio spectrum and NSGA-III-ELM algorithm were the best, as the R2 values of the training set and validation set were 0.852 and 0.810, and the root mean square error (RMSE) values were 0.291 and 0.308, respectively. From the perspective of the spectral data, the inversion accuracy of the ratio spectrum was better than the accuracy of the original spectrum or difference spectrum. At the algorithm level, the model inversion results based on LSTM algorithms showed a serious overfitting phenomenon and poor inversion effect. The inversion accuracy based on the NSGA-III-ELM algorithm was better than the accuracy of the MLR algorithm or the ELM algorithm. Therefore, the inversion model based on the ratio spectrum and NSGA-III-ELM algorithm could effectively invert the N deficiency in rice and provide critical technical support for accurate topdressing based on the N status in the rice.