Deriving Crop Phenology Metrics and Their Trends Using Times Series NOAA-AVHRR NDVI Data

In this study, an attempt has been made to derive the spatial patterns of temporal trends in phenology metrics and productivity of crops grown, at disaggregated level in Indo-Gangetic Plains of India (IGP), which are helpful in understanding the impact of climatic, ecological and socio-economic drivers. The NOAA-AVHRR NDVI PAL dataset from 1981 to 2001 was stacked as per the crop year and subjected to Savitzky-Golay filtering. For crop pixels, maximum and minimum values of normalized difference vegetation index (NDVI), their time of occurrence and total duration of kharif (June-October) and rabi (November–April) crop seasons were derived for each crop year and later subjected to pixel-wise regression with time to derive the rate and direction of change. The maximum NDVI value showed increasing trends across IGP during both kharif and rabi seasons indicating a general increase in productivity of crops. The trends in time of occurrence of peak NDVI during kharif dominated with rice showed that the maximum vegetative growth stage was happening early with time during study period across most of Punjab, North Haryana, Parts of Central and East Uttar Pradesh and some parts of Bihar and West Bengal. Only central parts of Haryana showed a delay in occurrence of maximum vegetative stage with time. During rabi, no significant trends in occurrence of peak NDVI were observed in most of Punjab and Haryana except in South Punjab and North Haryana where early occurrence of peak NDVI with time was observed. Most parts of Central and Eastern Uttar Pradesh, North Bihar and West Bengal showed a delay in occurrence of peak NDVI with time. In general, the rice dominating system was showing an increase in duration with time in Punjab, Haryana, Western Uttar Pradesh, Central Uttar Pradesh and South Bihar whereas in some parts of North Bihar and West Bengal a decrease in the duration with time was also observed. During rabi season, except Punjab, the wheat dominating system was showing a decreasing trend in crop duration with time.     

Detecting winter wheat phenology with SPOT-VEGETATION data in the North China Plain

Monitoring phenological change in agricultural land improves our understanding of the adaptation of crops to a warmer climate. Winter wheat–maize and winter wheat–cotton double-cropping are practised in most agricultural areas in the North China Plain. A curve-fitting method is presented to derive winter wheat phenology from SPOT-VEGETATION S10 normalized difference vegetation index (NDVI) data products. The method uses a double-Gaussian model to extract two phenological metrics, the start of season (SOS) and the time of maximum NDVI (MAXT). The results are compared with phenological records at local agrometeorological stations. The SOS and MAXT have close agreement with in situ observations of the jointing date and milk-in-kernel date respectively. The phenological metrics detected show spatial variations that are consistent with known phenological characteristics. This study indicates that time-series analysis with satellite data could be an effective tool for monitoring the phenology of crops and its spatial distribution in a large agricultural region.     

Derivation of crop phenological parameters using multi-date SPOT-VGT-NDVI data: A case study for Punjab

Crop phenological parameters, such as the start and end time of the crop growth, the total length of the growing season, time of peak vegetation and rate of greening and senescence are important for planning crop management and crop diversification/intensification. Multi-temporal remote sensing data provides opportunity to characterize the crop phenology at regional level. This study was conducted during the kharif season of the year 2001–02 for Punjab. The ten-day Normalised Difference Vegetation Index (NDVI) composite products, with 1 km spatial resolution, available from the Vegetation sensor onboard SPOT4 were used for the study. Twenty-one temporal datasets from May 1, 2001 to November 21, 2001 were used. Logical modelling approach was followed to compute the minimum and maximum NDVI, the amplitude of NDVI, the threshold NDVI during sowing and harvest, the crop duration, integrated NDVI and skewness of profile. The analysis showed that before July beginning, in the whole of Punjab, sowing/planting was over. It was found that the crop emergence in the eastern part of the state started earlier than the western region. The maximum NDVI, which represented peak vegetative stage, was above 0.7 and occurred mostly during August. The duration of crops ranged between 90–140 days, with majority between 110–120 days. Total integrated NDVI in Punjab was generally above 60. Using principal component analysis and divergence analysis seven best metrics were selected for crop discrimination.     

Extracting Regional Pattern of Wheat Sowing Dates Using Multispectral and High Temporal Observations from Indian Geostationary Satellite

Monitoring of Agricultural crops using remote sensing data is an emerging tool in recent years. Spatial determination of sowing date is an important input of any crop model. Geostationary satellite has the capability to provide data at high temporal interval to monitor vegetation throughout the entire growth period. A study was conducted to estimate the sowing date of wheat crop in major wheat growing states viz. Punjab, Haryana, Uttar Pradesh (UP), Madhya Pradesh (MP), Rajasthan and Bihar. Data acquired by Charged Couple Detector (CCD) onboard Indian geostationary satellite INSAT 3A have continental (Asia) coverage at 1 km?×?1 km spatial resolution in optical spectral bands with high temporal frequency. Daily operational Normalized Difference Vegetation Index (NDVI) product from INSAT 3A CCD available through Meteorological and Oceanographic Satellite Data Archival Centre (MOSDAC) was used to estimate sowing date of wheat crop in selected six states. Daily NDVI data acquired from September 1, 2010 to December 31, 2010 were used in this study. A composite of 7 days was prepared for further analysis of temporal profile of NDVI. Spatial wheat crop map derived from AWiFS (56 m) were re-sampled at INSAT 3A CCD parent resolution and applied over each 7 day composite. The characteristic temporal profiles of 7 day NDVI composite was used to determine sowing date. NDVI profile showed decreasing trend during maturity of kharif crop, minimum value after harvest and increasing trend after emergence of wheat crop. A mathematical model was made to capture the persistent positive slope of NDVI profile after an inflection point. The change in behavior of NDVI profile was detected on the basis of change in NDVI threshold of 0.3 and sowing date was estimated for wheat crop in six states. Seven days has been deducted after it reached to threshold value with persistent positive slope to get sowing date. The clear distinction between early sowing and late sowing regions was observed in study area. Variation of sowing date was observed ranging from November 1 to December 20. The estimated sowing date was validated with the reported sowing date for the known wheat crop regions. The RMSD of 3.2 (n?=?45) has been observed for wheat sowing date. This methodology can also be applied over different crops with the availability of crop maps.     

On the prognostics of growth profile based yield modelling

In order to improve the prognostics of yield forecasts two approaches have been explored using NDVI-based growth profiles for wheat crop of 1987-88 and 1990-91 seasons for some districts of Punjab and Haryana. Correlation of yield with variables based on profile area segments and with product of profile segment area and time to peak occurrence of growth cycle have been investigated. While the correlations are low and inconsistent for area variables, the îndex time product moment (IIPM) variable shows consistent and significant correlations and advances the date of forecast by 45-50 days over other approaches.     

Differences between cropland and rangeland MODIS phenology (start-of-season) in Mali

Start-of-season data are more and more used to qualify the land surface phenology trends in relation with climate variability and, more rarely, with human land management. In this paper, we compared the phenology of rangeland vs cropped land in the Sahel belt of Africa, using the only currently available global phenology product (MODIS MCD12Q2 – Land Cover Dynamics Yearly), and an enhanced crop mask of Mali. The differences in terms of start-of-season (SOS) are spatially (north south gradient), and temporally (10 years, 2001–2009) analyzed in bioclimatic terms. Our results show that globally the MODIS MCD12Q2 SOS dates of croplands and rangelands differ, and that these differences depend on the bioclimatic zone. In Sahelian and Guinean regions, cropland vegetation begins to grow earlier than rangeland vegetation (8-day and 4-day advance, respectively). Between, in the Sudanian and Sudano-Sahelian parts of Mali, rangeland vegetation greens about one week earlier than croplands. These results are discussed in the context of the land surface heterogeneity at MODIS scale, and in the context of the natural vegetation ecology. These results could help interpreting phenological trends in climate change analysis.     

Wheat Monitoring by Using Kernel Based Possibilistic <Emphasis Type="Italic">c</Emphasis>-Means Classifier: Bi-sensor Temporal Multi-spectral Data

The existence of mixed pixels in the satellite images has always been an area of concern. Adding to the challenge is an occurrence of non-linearity between the classes, which is generally overlooked. The study makes an attempt to solve the two frequently occurring problems by kernel based fuzzy approach. This research work deals with Possibilistic c-Means (PCM) classifier with local, global, spectral angle and hyper tangent kernels for wheat crop (Triticum aestivum) identification in Haridwar, Uttarakhand, India. The multi-temporal vegetation index data of Formosat-2 have been used which covers the whole phenology of wheat crop. The additional sensor Landsat-8 OLI imagery has been filled the crucial gap of Formosat-2 temporal datasets. Nine types of proposed kernels based PCM classifier have been applied on three temporal datasets (four, five and six date combinations) to classify two classes early sown and late sown wheat crop. These test results have been concluded that at optimized weighted constant KMOD and polynomial kernel was found effective to separate wheat crop. The five and six date combination were sufficient to discriminate early sown and late sown wheat crop.     

Relation of wheat yield with parameters derived from a spectral growth profile

An attempt has been made to generate crop growth profiles using multi-date NOAA AVHRR data of wheat-growing season of 1987–88 for the districts of Punjab and Haryana states of India. A profile model proposed by Badhwar was fitted to the multi-date Normalised Difference Vegetation Index (NDVI) values obtained from geographically referenced samples in each district. A novel approach of deriving a set of physiologically meaningful profile parameters has been outlined and the relation of these parameters with district wheat yields has been studied in order to examine the potential of growth profiles for crop-yield modelling. The parameter ‘area under the profile’ is found to be the best estimator of yield. However, with such a parameter time available for prediction gets reduced. Combination of different profile parameters shows improvement in correlation but lacks the consistency for individual state data.     

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.     

A comparative analysis of the NDVIg and NDVI3g in monitoring vegetation phenology changes in the Northern Hemisphere

Phenology is a sensitive and critical feature of vegetation and is a good indicator for climate change studies. The global inventory modelling and mapping studies (GIMMS) normalized difference vegetation index (NDVI) has been the most widely used data source for monitoring of the vegetation dynamics over large geographical areas in the past two decades. With the release of the third version of the NDVI (GIMMS NDVI3g) recently, it is important to compare the NDVI3g data with those of the previous version (NDVIg) to link existing studies with future applications of the NDVI3g in monitoring vegetation phenology. In this study, the three most popular satellite start of vegetation growing season (SOS) extraction methods were used, and the differences between SOSg and SOS3g arising from the methods were explored. The amplitude and the peak values of the NDVI3g are higher than those of the NDVIg curve, which indicated that the SOS derived from the NDVIg (SOSg) was significantly later than that derived from the NDVI3g (SOS3g) based on all the methods, for the whole northern hemisphere. In addition, SOSg and SOS3g both showed an advancing trend during 1982–2006, but that trend was more significant with SOSg than with SOS3g in the results from all three methods. In summary, the difference between SOSg and SOS3g (in the multi-year mean SOS, SOS change slope and the turning point in the time series) varied among the methods and was partly related to latitude. For the multi-year mean SOS, the difference increased with latitude intervals in the low latitudes (0–30°N) and decreased in the mid- and high-latitude intervals. The GIMMS NDVI3g data-sets seemed more sensitive than the GIMMS NDVIg in detecting information about the ground, and the SOS3g data were better correlated both with the in situ observations and the SOS derived from the Moderate Resolution Imaging Spectroradiometer NDVI. For the northern hemisphere, previous satellite measures (SOS derived from GIMMS NDVIg) may have overestimated the advancing trend of the SOS by an average of 0.032 d yr^–1.     

Spectral wheat growth profile in Punjab using IRS WiFS data

A functional form of crop spectral profile suggested by Badhwar was applied to district-wise wheat Normalised Difference Vegetation Index (NDVI) values relatively normalised by Pseudo-Invariant Feature (urban and built-up) NDVI values, derived from Wide Field Sensor (WiFS) onboard Indian Remote Sensing Satellites (IRS) for 17 dates during 1999–2000 rabi season. The goodness of overall profile fitting and the three basic parameters i.e., crop emergence date (T_o), and crop specific parameters (a and P) was found to be statistically significant. While a corresponds to profile progressive growth rate, β corresponds to profile decay rate. A comparison with earlier studies in Punjab using NOAA-AVHRR indicated improvement in relation between peak NDVI and wheat yield. The estimated time of spectral emergence and profile-derived peak NDVI follow the observed behaviour of shortened crop pre-anthesis period with delayed sowing.     

Computerized Database of Salt Affected Soils for Agro‐climatic Regions in the Indo–Gangetic Plain of India Using GIS

Abstract

Indo_Gangetic Plain (IGP) of India that stretched from the foothills of Himalayas near the Punjab State to the Gangetic delta in West Bengal State was known for highly fertile soil and favorable climatic condition for highest production of rice‐wheat. Appearance of soil salinity in large areas of IGP caused a major concern due to loss of productivity. The salt affected soils maps of India (NRSA 1997) showed vast areas of salt affected soils distributed along the Gangetic Plain covering the States of Haryana, Punjab, Uttar Pradesh, Bihar and West Bengal. In the analogue form, these maps contain voluminous data were difficult to handle without messing the whole dataset. An attempt was made to prepare a digitized database of salt affected soils to facilitate easy access, retrieval and map calculations required for reclamation and management of salt affected soil. The salt affected soils maps on 1:250, 000 scale were digitized for the States of Punjab, Haryana, Uttar Pradesh, Bihar and West Bengal using ILWIS. GIS. The Survey of India topomap was used for geo‐referencing and basemap preparation overlaying thematic layers for administrative and political boundaries, infrastructure, irrigation and drainage and settlements. The attribute data on physiography and the soil characteristics were stored in an attribute table and linked with the digitized polygons to prepare a relational database. Combining geo‐referenced (State) maps of Haryana, Punjab, Uttar Pradesh, Bihar and West Bengal using GIS, a composite map for Indo‐Gangetic plain was prepared. Four Agroclimatic regions (ACRs) and seventeen Agroclimatic zones (ACZs) were identified in the Indo‐Gangetic Plain (The Planning Commission of India) for planning and development of natural resources at regional level. The boundaries of ACZs and ACRs were delineated from the primary (master) database of IGP using ILWIS.GIS. The distribution of SAS polygons at regional and zonal level was delineated superimposing digitized boundaries of ACRs and ACZs over the master database of IGP. The state‐wise, region‐wise and zone‐wise extent of SAS was calculated. Soils were essentially saline at Lower‐ and Middle Gangetic Plain regions but highly variable and complex saline‐sodic in the Upper‐ and Trans‐Gangetic Plain regions. The area statistics showed that maximum SAS area occurred in ACR V (Upper Gangetic Plain) in Uttar Pradesh (UP) followed by ACR IV (Middle Gangetic Plain) in UP and Bihar, ACR III (Lower Gangetic Plain) in West Bengal and ACR VI (Trans‐Gangetic Plain) of Haryana and Punjab. Such database in digital format provides geo‐referenced, easy to access and retrievable, relational database comprising of thematic and attribute information of salt affected soils at state, regional and zonal level to facilitate overlay and map calculation of related data such as water quality, climatic, landform etc, useful for planning and decision making in reclamation and management of salt affected soils in IGP and other similar regions.     

A Spatial Database of Cropping System and its Characteristics to Aid Climate Change Impact Assessment Studies

Cropping system study is not only useful to understand the overall sustainability of agricultural system, but also it helps in generating many important parameters which are useful in climate change impact assessment. Considering its importance, Space Applications Centre, took up a project for mapping and characterizing major cropping systems of Indo-Gangetic Plains of India. The study area included the five states of Indo-Gangetic Plains (IGP) of India, i.e. Punjab, Haryana, Uttar Pradesh, Bihar and West Bengal. There were two aspects of the study. The first aspect included state and district level cropping system mapping using multi-date remote sensing (IRS-AWiFS and Radarsat ScanSAR) data. The second part was to characterize the cropping system using moderate spatial resolution multi-date remote sensing data (SPOT VGT NDVI) and ground survey. The remote sensing data was used to compute three cropping system performance indices (Multiple Cropping Index, Area Diversity Index and Cultivated Land Utilization Index). Ground survey was conducted using questionnaires filled up by 1,000 farmers selected from 103 villages based on the cropping systems map. Apart from ground survey, soil and water sampling and quality analysis were carried out to understand the effect of different cropping systems and their management practices. The results showed that, rice-wheat was the major cropping system of the IGP, followed by Rice-Fallow-Fallow and Maize-Wheat. Other major cropping systems of IGP included Sugarcane based, Pearl millet-Wheat, Rice-Fallow-Rice, Cotton-Wheat. The ground survey could identify 77 cropping systems, out of which 38 are rice-based systems. Out of these 77 cropping systems, there were 5 single crop systems, occupying 6.5% coverage (of all cropping system area), 56 double crop systems with 72.7% coverage, and 16 triple crop systems with 20.8% coverage. The cropping system performance analysis showed that the crop diversity was found to be highest in Haryana, while the cropping intensity was highest in Punjab state.     

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.     

Mapping crop phenology using NDVI time-series derived from HJ-1 A/B data

With the availability of high frequent satellite data, crop phenology could be accurately mapped using time-series remote sensing data. Vegetation index time-series data derived from AVHRR, MODIS, and SPOT-VEGETATION images usually have coarse spatial resolution. Mapping crop phenology parameters using higher spatial resolution images (e.g., Landsat TM-like) is unprecedented. Recently launched HJ-1 A/B CCD sensors boarded on China Environment Satellite provided a feasible and ideal data source for the construction of high spatio-temporal resolution vegetation index time-series. This paper presented a comprehensive method to construct NDVI time-series dataset derived from HJ-1 A/B CCD and demonstrated its application in cropland areas. The procedures of time-series data construction included image preprocessing, signal filtering, and interpolation for daily NDVI images then the NDVI time-series could present a smooth and complete phenological cycle. To demonstrate its application, TIMESAT program was employed to extract phenology parameters of crop lands located in Guanzhong Plain, China. The small-scale test showed that the crop season start/end derived from HJ-1 A/B NDVI time-series was comparable with local agro-metrological observation. The methodology for reconstructing time-series remote sensing data had been proved feasible, though forgoing researches will improve this a lot in mapping crop phenology. Last but not least, further studies should be focused on field-data collection, smoothing method and phenology definitions using time-series remote sensing data.     

Evaluations and comparisons of rule-based and machine-learning-based methods to retrieve satellite-based vegetation phenology using MODIS and USA National Phenology Network data

Vegetation phenology is a sensitive indicator that reflects the vegetation–atmosphere interactions and vegetation processes under global atmospheric changes. Fast-developing remote sensing technologies that monitor the land surface at high spatial and temporal resolutions have been widely used in vegetation phenology retrieval and analysis at a large scale. While researchers have developed many phenology retrieving methods based on remote sensing data, the relationships and differences among the phenology retrieving methods are unclear, and there is a lack of evaluation and comparison with the field phenology recoding data. In this study, we evaluated and compared eight phenology retrieving methods using Moderate Resolution Imaging Spectroradiometer (MODIS) and the USA National Phenology Network data from across North America. The studied phenology retrieving methods included six commonly used rule-based methods (i.e., amplitude threshold, the first-order derivative, the second-order derivative, the third-order derivative, the relative change curvature, and the curvature change rate) and two newly developed machine learning methods (i.e., neural network and random forest). At the large scale, the start of the season (SOS) values, derived by all methods, had similar spatial distributions; however, the retrieved values had large uncertainties in each pixel, and the end of the season (EOS) inverted values were largely different among methods. At the site scale, the SOS and EOS values extracted by the rule-based methods all had significant positive correlations with the field phenology observations. Among the rule-based methods, the amplitude threshold method performed the best. The machine learning methods outperformed the rule-based methods in terms of retrieving the SOS when assessed using the field observations. Our study highlighted that there were large differences among the methods in retrieving the vegetation phenology from satellite data and that researchers must be cautious in selecting an appropriate method for analyzing the satellite-retrieved phenology. Our results also demonstrated the importance of field phenology observations and the usefulness of the machine learning methods in understanding the satellite-based land surface phenology. These findings provide a valuable reference for the future development of global and regional phenology products.     

Crop polarisation behavior derived from Multi-Temporal adeos-polder data

This paper reports a study on multi-temporal polarized response of wheat crop from spaceborne ADEOS-POLDER sensor over a homogeneous wheat region of Punjab, India. Both the polarized as well as total reflectance of wheat were observed at different scattering angles for two spectral bands i.e. 670 nm and 865 nm during crop growth from November to April in rabi 1996-97 season. Results show that sun-target-viewing geometry plays an important role in polarization property. The top of atmosphere (TOA) polarized reflectance is found to decrease exponentially with increasing scattering angle. Polarized reflectance of crop was found to be an order of magnitude smaller in comparison to the total reflectance. An attempt was also made to model the observed polarized behavior over an agricultural area using a theoretical simplified crop reflectance model and accounting for atmospheric molecular (Rayleigh) contribution in the single scattering approximation. It was found that there was a decrease in the polarized reflectance at the grain filling (heading) stage of wheat crop. This is in accordance with ground- based observations and can be due to the reduction in the specular component of the reflected light during post-heading stage of the crop.     

青藏高原纳木错湖冰物候变化遥感监测与模拟

湖冰物候是反映区域气候变化的直观指标.由于青藏高原湖冰物候的地面观测不足,遥感与模拟成为动态监测湖冰物候变化并揭示其变化机理的重要途径.本文以纳木错为例,通过不同遥感方法获取了纳木错2000年-2015年湖冰物候的动态变化.在此基础上,将遥感与物理基础清晰的湖泊过程模型相结合,重建了纳木错1963年-2018年的湖冰物...     

Characterizing spring phenology of temperate broadleaf forests using Landsat and Sentinel-2 time series

Vegetation phenology has a great impact on land-atmosphere interactions like carbon cycling, albedo, and water and energy exchanges. To understand and predict these critical land-atmosphere feedbacks, it is crucial to measure and quantify phenological responses to climate variability, and ultimately climate change. Coarse-resolution sensors such as MODIS and AVHRR have been useful to study vegetation phenology from regional to global scales. These sensors are, however, not capable of discerning phenological variation at moderate spatial scales. By offering increased observation density and higher spatial resolution, the combination of Landsat and Sentinel-2 time series might provide the opportunity to overcome this limitation.In this study, we analyzed the potential of combined Sentinel-2 and Landsat time series for estimating start of season (SOS) of broadleaf forests across Germany for the year 2018. We tested two common statistical modeling approaches (logistic and generalized additive models using thin plate splines) and the two most commonly used vegetation indices, the Normalized Difference Vegetation Index (NDVI) and the Enhanced Vegetation Index (EVI).We found strong agreement between SOS estimates from logistic and spline models (r_EVI = 0.86; r_NDVI = 0.65), whereas agreement was higher for EVI than for NDVI (RMSD_EVI = 3.07, RMSD_NDVI = 5.26 days). The choice of vegetation index thus had a higher impact on the results than the fitting method. The EVI-based SOS also showed higher correlation with ground observations compared to NDVI (r_EVI = 0.51, r_NDVI = 0.42). Data density played an important role in estimating land surface phenology. Models combining Sentinel-2A/B, with an average cloud-free observation frequency of 12 days, were largely consistent with the combined Landsat and Sentinel-2 models, suggesting that Sentinel-2A/B may be sufficient to capture SOS for most areas in Germany in 2018. However, in non-overlapping swath areas and mountain areas, observation frequency was significantly lower, underlining the need to combine Landsat and Sentinel-2 for consistent SOS estimates over large areas. Our study demonstrates that estimating SOS of temperate broadleaf forests at medium spatial resolution has become feasible with combined Landsat and Sentinel-2 time series.     

利用NOAA NDVI数据集监测冬小麦生育期的研究

探索了利用NDVI研究作物生育期的方法，对黄淮海冬麦区的返青期、抽穗期、成熟期进行了估测，并利用地面实际观测资料进行了验证。结果表明，NDVI数据对大范围农作物生育期监测是可行的。冬小麦遥感反青期由南到北依次推迟，符合春季绿波由南到北推移规律。对冬小麦遥感生育期年际变化分析表明，黄淮海平原返青期变化相对较大，而抽穗期和成熟期变化较小。根据历年月平均温度与返青期分析，冬小麦返青日期与2月份平均温度密切相关。对于局部地区，利用5d合成1km分辨率数据，且按农业生态分区分别制定生育期判别标准，估测效果将更好。