A simple assessment approach for winter wheat loss risk impacted by water stress

In this study, an empirical assessment approach for the risk of crop loss due to water stress was developed and used to evaluate the risk of winter wheat loss in China, the United States, Germany, France and the United Kingdom. We combined statistical and remote sensing data on crop yields with climate data and cropland distribution to model the effect of water stress from 1982 to 2011. The average value of winter wheat loss due to water stress for the three European countries was about ?931 kg/ha, which was higher than that in China (?570 kg/ha) and the United States (?367 kg/ha). Our study has important implications for the operational assessment of crop loss risk at a country or regional scale. Future studies should focus on using higher spatial resolution remote sensing data, combining actual evapotranspiration to estimate water stress, improving the method for downscaling of statistical crop yield data and establishing more sophisticated zoning methods.     

Winter wheat yield forecasting in Ukraine based on Earth observation,meteorological data and biophysical models

Ukraine is one of the most developed agriculture countries and one of the biggest crop producers in the world. Timely and accurate crop yield forecasts for Ukraine at regional level become a key element in providing support to policy makers in food security. In this paper, feasibility and relative efficiency of using moderate resolution satellite data to winter wheat forecasting in Ukraine at oblast level is assessed. Oblast is a sub-national administrative unit that corresponds to the NUTS2 level of the Nomenclature of Territorial Units for Statistics (NUTS) of the European Union. NDVI values were derived from the MODIS sensor at the 250 m spatial resolution. For each oblast NDVI values were averaged for a cropland map (Rainfed croplands class) derived from the ESA GlobCover map, and were used as predictors in the regression models. Using a leave-one-out cross-validation procedure, the best time for making reliable yield forecasts in terms of root mean square error was identified. For most oblasts, NDVI values taken in April–May provided the minimum RMSE value when comparing to the official statistics, thus enabling forecasts 2–3 months prior to harvest. The NDVI-based approach was compared to the following approaches: empirical model based on meteorological observations (with forecasts in April–May that provide minimum RMSE value) and WOFOST crop growth simulation model implemented in the CGMS system (with forecasts in June that provide minimum RMSE value). All three approaches were run to produce winter wheat yield forecasts for independent datasets for 2010 and 2011, i.e. on data that were not used within model calibration process. The most accurate predictions for 2010 were achieved using the CGMS system with the RMSE value of 0.3 t ha⁻¹ in June and 0.4 t ha⁻¹ in April, while performance of three approaches for 2011 was almost the same (0.5–0.6 t ha⁻¹ in April). Both NDVI-based approach and CGMS system overestimated winter wheat yield comparing to official statistics in 2010, and underestimated it in 2011. Therefore, we can conclude that performance of empirical NDVI-based regression model was similar to meteorological and CGMS models when producing winter wheat yield forecasts at oblast level in Ukraine 2–3 months prior to harvest, while providing minimum requirements to input datasets.     

Regional yield estimation for winter wheat with MODIS-NDVI data in Shandong, China

The significance of crop yield estimation is well known in agricultural management and policy development at regional and national levels. The primary objective of this study was to test the suitability of the method, depending on predicted crop production, to estimate crop yield with a MODIS-NDVI-based model on a regional scale. In this paper, MODIS-NDVI data, with a 250 m resolution, was used to estimate the winter wheat (Triticum aestivum L.) yield in one of the main winter-wheat-growing regions. Our study region is located in Jining, Shandong Province. In order to improve the quality of remote sensing data and the accuracy of yield prediction, especially to eliminate the cloud-contaminated data and abnormal data in the MODIS-NDVI series, the Savitzky–Golay filter was applied to smooth the 10-day NDVI data. The spatial accumulation of NDVI at the county level was used to test its relationship with winter wheat production in the study area. A linear regressive relationship between the spatial accumulation of NDVI and the production of winter wheat was established using a stepwise regression method. The average yield was derived from predicted production divided by the growing acreage of winter wheat on a county level. Finally, the results were validated by the ground survey data, and the errors were compared with the errors of agro-climate models. The results showed that the relative errors of the predicted yield using MODIS-NDVI are between −4.62% and 5.40% and that whole RMSE was 214.16 kg ha⁻¹ lower than the RMSE (233.35 kg ha⁻¹) of agro-climate models in this study region. A good predicted yield data of winter wheat could be got about 40 days ahead of harvest time, i.e. at the booting-heading stage of winter wheat. The method suggested in this paper was good for predicting regional winter wheat production and yield estimation.     

MODIS植被指数的美国玉米单产遥感估测

针对中国开展的国外农作物产量遥感估测大多依靠中低分辨率耕地信息、省级(州级)或国家级作物产量统计数据的现状,本文以美国玉米为例,探讨利用多年中高分辨率作物分布信息、时序遥感植被指数和县级作物产量统计数据开展国外重点地区作物单产遥感估测技术研究,以期进一步提高中国对国外农作物产量监测精度和精细化水平。首先,利用美国农业部国家农业统计局(NASS/USDA)生产的作物分布数据(CDL)获得多个年份玉米空间分布图,并对相应年份250 m分辨率16天合成的MODIS-NDVI时序数据进行掩膜处理,统计获得每年各县域内玉米主要生育期NDVI均值;其次,以各州为估产区,以多年县级玉米统计单产和县域内玉米主要生育期NDVI均值为基础,建立各州玉米主要生育期NDVI与玉米单产间关系模型;然后,通过主要生育期玉米单产和玉米植被指数间拟合程度,筛选确定各州玉米最佳估产期和最佳估产模型。最终,利用最佳估产模型实现美国各州玉米单产估测和全国玉米单产推算。其中,建模数据覆盖时间为2007年—2010年,验证数据为2011年。结果表明,应用最佳估产模型的2011年美国各州玉米单产估测相对误差在-4.16%—4.92%,均方根误差在148.75—820.93 kg/ha,各州估测结果计算获得全国玉米单产的相对误差仅为2.12%,均方根误差为285.57 kg/ha。可见,本研究的作物单产遥感估测技术方法具有一定可行性,可准确估测全球重点地区作物单产信息。     

A study on trade-offs between spatial resolution and temporal sampling density for wheat yield estimation using both thermal and calendar time

Within-season forecasting of crop yields is of great economic, geo-strategic and humanitarian interest. Satellite Earth Observation now constitutes a valuable and innovative way to provide spatio-temporal information to assist such yield forecasts. This study explores different configurations of remote sensing time series to estimate of winter wheat yield using either spatially finer but temporally sparser time series (5daily at 100 m spatial resolution) or spatially coarser but denser (300 m and 1 km at daily frequency) time series. Furthermore, we hypothesised that better yield estimations could be made using thermal time, which is closer to the crop physiological development. Time series of NDVI from the PROBA-V instrument, which has delivered images at a spatial resolution of 100 m, 300 m and 1 km since 2013, were extracted for 39 fields for field and 56 fields for regional level analysis across Northern France during the growing season 2014-2015. An asymmetric double sigmoid model was fitted on the NDVI series of the central pixel of the field. The fitted model was subsequently integrated either over thermal time or over calendar time, using different baseline NDVI thresholds to mark the start and end of the cropping season. These integrated values were used as a predictor for yield using a simple linear regression and yield observations at field level. The dependency of this relationship on the spatial pixel purity was analysed for the 100 m, 300 m and 1 km spatial resolution. At field level, depending on the spatial resolution and the NDVI threshold, the adjusted R² ranged from 0.20 to 0.74; jackknifed – leave-one-field-out cross validation – RMSE ranged from 0.6 to 1.07 t/ha and MAE ranged between 0.46 and 0.90 t/ha for thermal time analysis. The best results for yield estimation (adjusted R² = 0.74, RMSE =0.6 t/ha and MAE =0.46 t/ha) were obtained from the integration over thermal time of 100 m pixel resolution using a baseline NDVI threshold of 0.2 and without any selection based on pixel purity. The field scale yield estimation was aggregated to the regional scale using 56 fields. At the regional level, there was a difference of 0.0012 t/ha between thermal and calendar time for average yield estimations. The standard error of mean results showed that the error was larger for a higher spatial resolution with no pixel purity and smaller when purity increased. These results suggest that, for winter wheat, a finer spatial resolution rather than a higher revisit frequency and an increasing pixel purity enable more accurate yield estimations when integrated over thermal time at the field scale and at the regional scale only if higher pixel purity levels are considered. This method can be extended to larger regions, other crops, and other regions in the world, although site and crop-specific adjustments will have to include other threshold temperatures to reflect the boundaries of phenological activity. In general, however, this methodological approach should be applicable to yield estimation at the parcel and regional scales across the world.     

Integration of maximum crop response with machine learning regression model to timely estimate crop yield

Timely and reliable estimation of regional crop yield is a vital component of food security assessment, especially in developing regions. The traditional crop forecasting methods need ample time and labor to collect and process field data to release official yield reports. Satellite remote sensing data is considered a cost-effective and accurate way of predicting crop yield at pixel-level. In this study, maximum Enhanced Vegetation Index (EVI) during the crop-growing season was integrated with Machine Learning Regression (MLR) models to estimate wheat and rice yields in Pakistan's Punjab province. Five MLR models were compared using a fivefold cross-validation method for their predictive accuracy. The study results revealed that the regression model based on the Gaussian process outperformed over other models. The best performing model attained coefficient of determination (R2), Root Mean Square Error (RMSE, t/ ha), and Mean Absolute Error (MAE, t/ha) of 0.75, 0.281, and 0.236 for wheat; 0.68, 0.112, and 0.091 for rice, respectively. The proposed method made it feasible to predict wheat and rice 6– 8 weeks before the harvest. The early prediction of crop yield and its spatial distribution in the region can help formulate efficient agricultural policies for sustainable social, environmental, and economic progress.     

An analysis of cropland mask choice and ancillary data for annual corn yield forecasting using MODIS data

The Midwestern United States is one of the world’s most important corn-producing regions. Monitoring and forecasting of corn yields in this intensive agricultural region are important activities to support food security, commodity markets, bioenergy industries, and formation of national policies. This study aims to develop forecasting models that have the capability to provide mid-season prediction of county-level corn yields for the entire Midwestern United States. We used multi-temporal MODIS NDVI (normalized difference vegetation index) 16-day composite data as the primary input, with digital elevation model (DEM) and parameter-elevation relationships on independent slopes model (PRISM) climate data as additional inputs. The DEM and PRISM data, along with three types of cropland masks were tested and compared to evaluate their impacts on model predictive accuracy. Our results suggested that the use of general cropland masks (e.g., summer crop or cultivated crops) generated similar results compared with use of an annual corn-specific mask. Leave-one-year-out cross-validation resulted in an average R² of 0.75 and RMSE value of 1.10 t/ha. Using a DEM as an additional model input slightly improved performance, while inclusion of PRISM climate data appeared not to be important for our regional corn-yield model. Furthermore, our model has potential for real-time/early prediction. Our corn yield esitmates are available as early as late July, which is an improvement upon previous corn-yield prediction models. In addition to annual corn yield forecasting, we examined model uncertainties through spatial and temporal analysis of the model's predictive error distribution. The magnitude of predictive error (by county) appears to be associated with the spatial patterns of corn fields in the study area.     

小麦单产的遥感—气象综合模式研究

本文以北京顺义县为例,以气象因子与垂直植被指数(PVI)作为参数,用灰色模型G(0,2)和逐段订正模型即阶乘模型,建立冬小表遥感信息-气象因子综合模型。计算结果表明,改进后的综合模型其平均精度比单纯的遥感信息模型提高近7％,个别年份达到10％以上。     

MODIS与TM冬小麦种植面积遥感测量一致性研究——小区域实验研究

针对大范围冬小麦种植面积遥感测量业务化运行中所存在的不同空间分辨率遥感影像的获取能力、空间分辨率与测量精度之间相互制约的现实问题,选择典型实验区,采取全覆盖的多时相低分辨率MODIS数据与中高分辨率TM样区数据相结合的方法,以支持向量机（SVM）为主要技术手段,通过选取不同比例的TM样本量,对MODIS进行混合像元分解,并对MODIS冬小麦测量结果与TM测量结果进行一致性分析,进而提出了一套可用于大范围冬小麦种植面积遥感测量业务化推广的识别和精度检验方法。研究结果表明：当TM样本量大于5%时,相对于TM识别结果而言,MODIS像元精度可以稳定在82%以上;当TM样本量大于40%时,区域精度可以稳定在97%。因此,从理论上讲,在实际业务化运行中,只要能够获得监测区40%的中高分辨率的影像,采用本研究提出的多尺度冬小麦种植面积测量方法,基本上可以满足业务化运行的精度要求。     

利用双变量同化与PCA-Copula法的冬小麦单产估测

为了进一步提高冬小麦产量估测的精度,基于集合卡尔曼滤波算法和粒子滤波（particle filter, PF）算法,对CERES–Wheat模型模拟的冬小麦主要生育期条件植被温度指数（vegetation temperature condition index,VTCI）、叶面积指数（leaf area index, LAI）和中分辨率成像光谱仪（moderate-resolution imaging spectroradiometer, MODIS）数据反演的VTCI、LAI进行同化,利用主成分分析与Copula函数结合的方法构建单变量和双变量的综合长势监测指标,建立冬小麦单产估测模型,并通过对比分析选择最优模型,对2017—2020年关中平原的冬小麦单产进行估测。结果表明,单点尺度的同化VTCI、同化LAI均能综合反映MODIS观测值和模型模拟值的变化特征,且PF算法具有更好的同化效果;区域尺度下利用PF算法得到的同化VTCI和LAI所构建的双变量估产模型精度最高,与未同化VTCI和LAI构建的估产模型精度相比,研究区各县（区）的冬小麦估测单产与实际单产的均方根误差降低了56.25 kg/hm2,平均相对误差降低了1.51%,表明该模型能有效提高产量估测的精度,应用该模型进行大范围的冬小麦产量估测具有较好的适用性。     

Wheat crop production estimation using satellite data

This paper reports acreage, yield and production forecasting of wheat crop using remote sensing and agrometeorological data for the 1998–99 rabi season. Wheat crop identification and discrimination using Indian Remote Sensing (IRS) ID LISS III satellite data was carried out by supervised maximum likelihood classification. Three types of wheat crop viz. wheat-1 (high vigour-normal sown), wheat-2 (moderate vigour-late sown) and wheat-3 (low vigour-very late sown) have been identified and discriminated from each other. Before final classification of satellite data spectral separability between classes were evaluated. For yield prediction of wheat crop spectral vegetation indices (RVI and NDVI), agrometeorological parameters (ET_max and TD) and historical crop yield (actual yield) trend analysis based linear and multiple linear regression models were developed. The estimated wheat crop area was 75928.0 ha. for the year 1998–99, which sowed ?2.59% underestimation with land record commissioners estimates. The yield prediction through vegetation index based and vegetation index with agrometeorological indices based models were 1753 kg/ha and 1754 kg/ha, respectively and have shown relative deviation of 0.17% and 0.22%, the production estimates from above models when compared with observed production show relative deviation of ?2.4% and ?2.3% underestimations, respectively.     

Estimating aboveground and organ biomass of plant canopies across the entire season of rice growth with terrestrial laser scanning

Non-destructive and accurate estimation of crop biomass is crucial for the quantitative diagnosis of growth status and timely prediction of grain yield. As an active remote sensing technique, terrestrial laser scanning (TLS) has become increasingly available in crop monitoring for its advantages in recording structural properties. Some researchers have attempted to use TLS data in the estimation of crop aboveground biomass, but only for part of the growing season. Previous studies rarely investigated the estimation of biomass for individual organs, such as the panicles in rice canopies, which led to the poor understanding of TLS technology in monitoring biomass partitioning among organs. The objective of this study was to investigate the potential of TLS in estimating the biomass for individual organs and aboveground biomass of rice and to examine the feasibility of developing universal models for the entire growing season. The field plots experiments were conducted in 2017 and 2018 and involved different nitrogen (N) rates, planting techniques and rice varieties. Three regression approaches, stepwise multiple linear regression (SMLR), random forest regression (RF) and linear mixed-effects (LME) modeling, were evaluated in estimating biomass with extensive TLS and biomass data collected at multiple phenological stages of rice growth across the entire season. The models were calibrated with the 2017 dataset and validated independently with the 2018 dataset.The results demonstrated that growth stage in LME modeling was selected as the most significant random effect on rice growth among the three candidates, which were rice variety, growth stage and planting technique. The LME models grouped by growth stage exhibited higher validation accuracies for all biomass variables over the entire season to varying degrees than SMLR models and RF models. The most pronounced improvement with a LME model was obtained for panicle biomass, with an increase of 0.74 in R² (LME: R² = 0.90, SMLR: R² = 0.16) and a decrease of 1.15 t/ha in RMSE (LME: RMSE =0.79 t/ha, SMLR: RMSE =2.94 t/ha). Compared to SMLR and RF, LME modeling yielded similar estimation accuracies of aboveground biomass for pre-heading stages, but significantly higher accuracies for post-heading stages (LME: R² = 0.63, RMSE =2.27 t/ha; SMLR: R² = 0.42, RMSE =2.42 t/ha; RF: R² = 0.57, RMSE =2.80 t/ha). These findings implied that SMLR was only suitable for the estimation of biomass at pre-heading stages and LME modeling performed remarkably well across all growth stages, especially for post-heading. The results suggest coupling TLS with LME modeling is a promising approach to monitoring rice biomass at post-heading stages at high accuracy and to overcoming the saturation of canopy reflectance signals encountered in optical remote sensing. It also has great potential in the monitoring of other crops in cloud-cover conditions and the instantaneous prediction of grain yield any time before harvest.     

Identification and mapping of winter wheat by integrating temporal change information and Kullback–Leibler divergence

Crop acreage and its spatial distribution are a base for agriculture related works. Current research combining medium and low spatial resolution images focuses on data fusion and unmixing methods. The purpose of the former is to generate synthetic fine spatial resolution data instead of directly solving the problem. In the latter, high-resolution data is only used to provide endmembers and the result is usually an abundance map rather than the true spatial distribution data. To solve this problem, this paper designs a conceptual model which divides the study area into different types of pixels at a MODIS 250 m scale. Only three types of pixels contain winter wheat, i.e., pure winter wheat pixels (P_A), the mixed pixels comprising winter wheat and other vegetation (M_A) and the mixed pixels comprising winter wheat and other crops (M_B). Different strategies are used in processing them. (1) Within the pure cultivated land pixels, the Kullback–Leibler (KL) divergence is employed to analyze the similarity between unknown pixels and the pure winter wheat samples on the temporal change characteristics of NDVI. Further P_A is identified. (2) For M_A, a proposed reverse unmixing method is firstly used to extract the temporal change information of cultivated land components, after which winter wheat is identified from the cultivated land components as previously described. (3) For M_B which only appears on the border of P_A, a mask is created by expanding the P_A and temporal difference is utilized to identify winter wheat under the mask. Finally, these three results are integrated at a TM scale with the aid of 25 m resolution land use data. We applied the proposed solution and obtained a good result in the main agricultural area of the Yiluo River Basin. The identified winter wheat planting acreage is 161,050.00 hm². The result is validated based on the five-hundred random validation points. Overall accuracy is 94.80% and Kappa coefficient is 0.85. This demonstrates that the temporal information reflecting crop growth is also an important indicator, and the KL divergence makes it more convenient in identifying winter wheat. This research provided a new perspective for the combination of low and medium spatial resolution remote sensing images. The proposed solution can also be effectively applied in other places and countries for the crop which has a clear temporal change characteristic that is different from others.     

Predicting agricultural drought in eastern Rajasthan of India using NDVI and standardized precipitation index

In the present study, prediction of agricultural drought has been addressed through prediction of agricultural yield using a model based on NDVI-SPI. It has been observed that the meteorological drought index SPI with different timescale is correlated with NDVI at different lag. Also NDVI of current fortnight is correlated with NDVI of previous lags. Based on the correlation coefficients, the Multiple Regression Model was developed to predict NDVI. The NDVI of current fortnight was found highly correlated with SPI of previous fortnight in semi-arid and transitional zones. The correlation between NDVI and crop yield was observed highest in first fortnight of August. The RMSE of predicted yield in drought year was found to be about 17.07 kg/ha which was about 6.02 per cent of average yield. In normal year, it was 24 kg/Ha denoting about 2.1 per cent of average yield.     

统计数据总量约束下全局优化阈值的冬小麦分布制图

大范围、长时间和高精度农作物空间分布基础农业科学数据的准确获取对资源、环境、生态、气候变化和国家粮食安全等问题研究具有重要现实意义和科学意义。本文针对传统阈值法农作物识别过程中阈值设置存在灵巧性差和自动化程度低等弱点,以中国粮食主产区黄淮海平原内河北省衡水市景县为典型实验区,首次将全局优化算法应用于阈值模型中阈值优化选取,开展了利用全局优化算法改进基于阈值检测的农作物分布制图方法创新研究。以冬小麦为研究对象,国产高分一号(GF-1)为主要遥感数据源,在作物面积统计数据为总量控制参考标准和全局参数优化的复合型混合演化算法SCE-UA (Shuffled Complex Evolution-University of Arizona)支持下,提出利用时序NDVI数据开展阈值模型阈值参数自动优化的冬小麦空间分布制图方法。最终,获得实验区冬小麦阈值模型最优参数,并利用优化后的阈值参数对冬小麦空间分布进行提取。通过地面验证表明,利用本研究所提方法获取的冬小麦识别结果分类精度均达到较高水平。其中冬小麦识别结果总量精度达到了99.99%,证明本研究所提阈值模型参数优化方法冬小麦提取分类结果总量控制效果良好;同时,与传统的阈值法、最大似然和支持向量机等分类方法相比,本研究所提阈值模型参数优化法区域冬小麦作物分类总体精度和Kappa系数分别都有所提高,其中,总体精度分别提高4.55%、2.43%和0.15%,Kappa系数分别提高0.12、0.06和0.01,这体现出SCE-UA全局优化算法对提高阈值模型冬小麦空间分布识别精度具有一定优势。以上研究结果证明了利用本研究所提基于作物面积统计数据总量控制以及SCE-UA全局优化算法支持下阈值模型参数优化作物分布制图方法的有效性和可行性,可获得高精度冬小麦作物空间分布制图结果,这对提高中国冬小麦空间分布制图精度和自动化水平具有一定意义,也可为农作物面积农业统计数据降尺度恢复重建和大范围区域作物空间分布制图研究提供一定技术参考。     

Assessing gaps in irrigated agricultural productivity through satellite earth observations—A case study of the Fergana Valley,Central Asia

Improving crop area and/or crop yields in agricultural regions is one of the foremost scientific challenges for the next decades. This is especially true in irrigated areas because sustainable intensification of irrigated crop production is virtually the sole means to enhance food supply and contribute to meeting food demands of a growing population. Yet, irrigated crop production worldwide is suffering from soil degradation and salinity, reduced soil fertility, and water scarcity rendering the performance of irrigation schemes often below potential. On the other hand, the scope for improving irrigated agricultural productivity remains obscure also due to the lack of spatial data on agricultural production (e.g. crop acreage and yield). To fill this gap, satellite earth observations and a replicable methodology were used to estimate crop yields at the field level for the period 2010/2014 in the Fergana Valley, Central Asia, to understand the response of agricultural productivity to factors related to the irrigation and drainage infrastructure and environment. The results showed that cropping pattern, i.e. the presence or absence of multi-annual crop rotations, and spatial diversity of crops had the most persistent effects on crop yields across observation years suggesting the need for introducing sustainable cropping systems. On the other hand, areas with a lower crop diversity or abundance of crop rotation tended to have lower crop yields, with differences of partly more than one t/ha yield. It is argued that factors related to the infrastructure, for example, the distance of farms to the next settlement or the density of roads, had a persistent effect on crop yield dynamics over time. The improvement potential of cotton and wheat yields were estimated at 5%, compared to crop yields of farms in the direct vicinity of settlements or roads. In this study it is highlighted how remotely sensed estimates of crop production in combination with geospatial technologies provide a unique perspective that, when combined with field surveys, can support planners to identify management priorities for improving regional production and/or reducing environmental impacts.     

Regional-scale monitoring of cropland intensity and productivity with multi-source satellite image time series

In the context of growing populations and limited resources, the sustainable intensification of agricultural production is of great importance to achieve food security. As the need to support management at a range of spatial scales grows, decision-support tools appear increasingly important to enable the timely and regular assessment of agricultural production over large areas and identify priorities for improving crop production in low-productivity regions. Understanding productivity patterns requires the timely provision of gapless, spatial information about agricultural productivity. In this study, dense 30-m time series covering the 2004–2014 period were generated from Landsat and MODerate-resolution Imaging Spectroradiometer (MODIS) satellite images over the irrigated cropped area of the Fergana Valley, Central Asia. A light-use efficiency model was combined with machine learning classifiers to assess the crop yield at the field level. The classification accuracy of land cover maps reached 91% on average. Crop yield and acreage estimates were in good agreement (R² = 0.812 and 0.871, respectively) with reported yields and acreages at the district level. Several indicators of cropland intensity and productivity were derived on a per-field basis and used to highlight homogeneous regions in terms of productivity by means of clustering. Results underlined that regions with lower water-use efficiency were not only located further away from irrigation canals and intake points, but also had limited access to markets and roads. The results underline that yield could be increased by roughly 1.0 and 1.4 t/ha for cotton and wheat, respectively, if the access to water would be optimized in some of the regions. The minimum calibration requirement of the method and the fusion of multi-sensor data are keys to cope with the constraints of operational crop monitoring and guarantee a sustained and timely delivery of the agricultural indicators to the user community. The results of this study can form the baseline to support regional land- and water-resource management.     

PhenoCrop: An integrated satellite-based framework to estimate physiological growth stages of corn and soybeans

The accurate and timely estimates of crop physiological growth stages are essential for efficient crop management and precise modeling of agricultural systems. Satellite remote sensing has been widely used to retrieve vegetation phenology metrics at local to global scales. However, most of these phenology metrics (e.g., green-up) are different from crop growth stages (e.g., emergence) used in crop management and modeling. As such, an integrated framework referred to as PhenoCrop was developed to: 1) establish a connection between remote sensing-derived phenology metrics and key crop growth stages based on Wang and Engle plant phenology model and 2) use fused MODIS-Landsat 30 m 8-day reflectance data generated using Kalman Filter-based data fusion technique to produce onset dates of key growth stages of corn (Zea mays L.) and soybeans (Glycine max L.) at 30 m spatial resolution. In this paper, we described the PhenoCrop framework, and tested its performance for the State of Nebraska for 2012–2016 by comparison to observations of estimated key growth stages at four experimental sites, and state-level statistical data from Crop Progress Reports (CPRs) published by the United States Department of Agriculture’s (USDA) National Agricultural Statistical Services (NASS). In addition, to evaluate the suitability of using coarse or high spatial resolution satellite imagery, fused MODIS-Landsat-based estimates were compared with those produced using EOS MODIS 250 m (MOD9Q1) reflectance data.The PhenoCrop estimates captured the typical spatial trends of gradual delay in the progression of the growing season from southeast to northwest Nebraska. Also inter-annual differences due to factors such as weather fluctuations and change in management strategies (e.g., early season in 2012) were evident in the estimates. Validation results revealed that average root mean square error (RMSE) of the state-level estimates of corn and soybean growth stages ranged from 1.10 to 4.20 days and from 3.81 to 7.89 days, respectively, while pixel level estimates had a RMSE ranging from 3.72 to 8.51 days for corn and 4.76–9.51 days for soybean growth stages. Although MODIS 250 m based estimates showed similar general spatial patterns observed in the fused MODIS-Landsat based estimates, the accuracy and ability to capture field scale variations was improved with fused MODIS-Landsat data. Overall, results showed the ability of PhenoCrop framework to provide reliable estimates of crop growth stages that can be highly useful in crop modeling and crop management during the growing season.     

Relationship of wheat grain yield with spectral indices

An experiment was conducted during 1996–97 and 1997–98 to study spectral indices and their relationships with grain yield of wheat. Variations of ratio vegetation index (RVI), normalized differences vegetation index (NDVI). difference vegetation index (DVI), transformed vegetation index (TVI), perpendicular vegetation index (PVI) and greenness vegetation index (GVI) have been studied at anthesis stage under different moisture and nitrogen levels. Spectral indices were correlated with crop parameters and it was found that GVI was the best index for yield estimation (r = 0.91 ).     

面向农作物监测的遥感信息处理技术研究

以区域性主要农作物种类识别、长势分析与产量估算及农业种植结构现状监测为主要研究对象,开展适于农业管理部门业务化运行的卫星遥感信息处理的关键技术研究。从分析主要作物类型识别的遥感物理依据入手,提出了卫星遥感数据处理及专题信息提取的基本技术框架、主要农作物类型及种植面积信息的提取方法以及主要粮食作物长势分析和产量估算模型,并对结果进行了简要的精度分析。