Long-term soil moisture content estimation using satellite and climate data in agricultural area of Mongolia

The purpose of this study is to estimate long-term SMC and find its relation with soil moisture (SM) of climate station in different depths and NDVI for the growing season. The study area is located in agricultural regions in the North of Mongolia. The Pearson’s correlation methodology was used in this study. We used MODIS and SPOT satellite data and 14 years data for precipitation, temperature and SMC of 38 climate stations. The estimated SMC from this methodology were compared with SM from climate data and NDVI. The estimated SMC was compared with SM of climate stations at a 10-cm depth (r² = 0.58) and at a 50-cm depth (r² = 0.38), respectively. From the analysis, it can be seen that the previous month’s SMC affects vegetation growth of the following month, especially from May to August. The methodology can be an advantageous indicator for taking further environmental analysis in the region.     

Crop area mapping in West Africa using landscape stratification of MODIS time series and comparison with existing global land products

In Africa, food security early warning systems use satellite-derived data concerning crop conditions and agricultural production. Such systems can be improved if they are provided with a more reliable estimation of the cultivated area at national scale. This paper evaluates the potential of using time series from the MODerate resolution Imaging Spectroradiometer MOD13Q1 (16-day composite of normalized difference vegetation index at 250 m resolution) to extract cultivated areas in the fragmented rural landscapes of Mali. To this end, we first stratified Southern Mali into 13 rural landscapes based on the spatio-temporal variability of NDVI and textural indices, using an object-oriented classification scheme.The accuracy of the resulting map (MODIS_crop) and how it compares with existing coarse-resolution global land products (GLC2000 Africa, GLOBCOVER, MODIS V05 and ECOCLIMAP-II), was then assessed against six crop/non-crop maps derived from SPOT 2.5 m resolution images used as references. For crop areal coverage, the MODIS_crop cultivated map was successful in assessing the overall cultivated area at five out of the six validation sites (less than 6% of the absolute difference), while in terms of crop spatial distribution, the producer accuracy was between 33.1% and 80.8%. This accuracy was linearly correlated with the mean patch size index calculated on the SPOT crop maps (r² = 0.8). Using the Pareto boundary as an accuracy assessment method at the study sites, we showed that (i) 20-40% of the classification crop error was due to the spatial resolution of the MODIS sensor (250 m), and that (ii) compared to MODIS V05, which otherwise performed better than the other existing products, MODIS_crop generally minimized omission-commission errors. A spatial validation of the different products was carried out using SPOT image classifications as reference. In the corresponding error matrices, the fraction of correctly classified pixels for our product was 70%, compared to 58% for MODIS V05, while it ranged between 40% and 51% for the GLC2000, the ECOCLIMAP-II and the GLOBCOVER.     

Estimating fractional vegetation cover and the vegetation index of bare soil and highly dense vegetation with a physically based method

Normalized difference vegetation index (NDVI) of highly dense vegetation (NDVI_v) and bare soil (NDVI_s), identified as the key parameters for Fractional Vegetation Cover (FVC) estimation, are usually obtained with empirical statistical methods However, it is often difficult to obtain reasonable values of NDVI_v and NDVI_s at a coarse resolution (e.g., 1 km), or in arid, semiarid, and evergreen areas. The uncertainty of estimated NDVI_s and NDVI_v can cause substantial errors in FVC estimations when a simple linear mixture model is used. To address this problem, this paper proposes a physically based method. The leaf area index (LAI) and directional NDVI are introduced in a gap fraction model and a linear mixture model for FVC estimation to calculate NDVI_v and NDVI_s. The model incorporates the Moderate Resolution Imaging Spectroradiometer (MODIS) Bidirectional Reflectance Distribution Function (BRDF) model parameters product (MCD43B1) and LAI product, which are convenient to acquire. Two types of evaluation experiments are designed 1) with data simulated by a canopy radiative transfer model and 2) with satellite observations. The root-mean-square deviation (RMSD) for simulated data is less than 0.117, depending on the type of noise added on the data. In the real data experiment, the RMSD for cropland is 0.127, for grassland is 0.075, and for forest is 0.107. The experimental areas respectively lack fully vegetated and non-vegetated pixels at 1 km resolution. Consequently, a relatively large uncertainty is found while using the statistical methods and the RMSD ranges from 0.110 to 0.363 based on the real data. The proposed method is convenient to produce NDVI_v and NDVI_s maps for FVC estimation on regional and global scales.     

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.     

Soil organic carbon mapping using remote sensing techniques and multivariate regression model

Soil organic carbon (SOC) is an important aspect of soil quality and plays an imperative role in soil productivity in the agriculture ecosystems. The present study was applied to estimate the SOC stock using space-borne satellite data (Landsat 4–5 Thematic Mapper [TM]) and ground verification in the Medinipur Block, Paschim Medinipur District and West Bengal in India. In total, 50 soil samples were collected randomly from the region according to field surveys using a hand-held Global Positioning System (GPS) unit to estimate the surface SOC concentrations in the laboratory. Bare soil index (BSI) and normalized difference vegetation ndex (NDVI) were explored from TM data. The satellite data-derived indices were used to estimate spatial distribution of SOC using multivariate regression model. The regression analysis was performed to determine the relationship between SOC and spectral indices (NDVI and BSI) and compared the observed SOC (field measure) to predict SOC (estimated from satellite images). Goodness fit test was performed to determine the significance of the relationship between observed and predicted SOC at p ≤ 0.05 level. The results of regression analysis between observed SOC and NDVI values showed significant relationship (R² = 0.54; p < 0.0075). A significant statistical relationship (r = ?0.72) was also observed between SOC and BSI. Finally, our model showed nearly 71% of the variance of SOC distribution could be explained by SOC and NDVI values. The information from this study has advanced our understanding of the ongoing ecological development that affects SOC dissemination and might be valuable for effective soil management.     

Fusing MODIS with Landsat 8 data to downscale weekly normalized difference vegetation index estimates for central Great Basin rangelands,USA

Data fused from distinct but complementary satellite sensors mitigate tradeoffs that researchers make when selecting between spatial and temporal resolutions of remotely sensed data. We integrated data from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor aboard the Terra satellite and the Operational Land Imager sensor aboard the Landsat 8 satellite into four regression-tree models and applied those data to a mapping application. This application produced downscaled maps that utilize the 30-m spatial resolution of Landsat in conjunction with daily acquisitions of MODIS normalized difference vegetation index (NDVI) that are composited and temporally smoothed. We produced four weekly, atmospherically corrected, and nearly cloud-free, downscaled 30-m synthetic MODIS NDVI predictions (maps) built from these models. Model results were strong with R² values ranging from 0.74 to 0.85. The correlation coefficients (r ≥ 0.89) were strong for all predictions when compared to corresponding original MODIS NDVI data. Downscaled products incorporated into independently developed sagebrush ecosystem models yielded mixed results. The visual quality of the downscaled 30-m synthetic MODIS NDVI predictions were remarkable when compared to the original 250-m MODIS NDVI. These 30-m maps improve knowledge of dynamic rangeland seasonal processes in the central Great Basin, United States, and provide land managers improved resource maps.     

Generation of fraction images from AVHRR data using linear mixing model

The coarse resolution satellite data have been widely used for regional and global studies as they provide high temporal frequency. The information contained in the coarse resolution pixels are mostly mixture of several components. The extraction of information contained in a pixel find its role in Geosphere-Biosphere context. The present study address the utility of constrained least square model applied to coarse spatial resolution data from NOAA-AVHRR for generating fraction images of vegetation, soil and water/shade. The red and near-infrared channels have been used to run the constrained least square model to generate fraction images. The derived fraction images are related to normalised difference vegetation index (NDVI) for model validation. The results suggest that vegetation fraction components are strongly correlated with NDVI values (r²=0.98). The soil fractions (r²=?0.84) and water/shade fractions (4²=?0.97) are negatively correlated with NDVI. The relationship between the fraction images and NDVI show the potential of the model in deriving sub-pixel component information using coarse resolution data.     

Evaluating and comparing Sentinel 2A and Landsat-8 Operational Land Imager (OLI) spectral indices for estimating fire severity in a Mediterranean pine ecosystem of Greece

The main purpose of this study is to explore the relationship between three field-based fire severity indices (Composite Burn Index-CBI, Geometrically structure CBI, weighted CBI) and spectral indices derived from Sentinel 2A and Landsat-8 OLI imagery on a recent large fire in Thasos, Greece. We employed remotely sensed indices previously used from the remote sensing fire community (Normalized Difference Vegetation Index (NDVI), Normalized Burn Ratio (NBR), differenced NDVI, differenced NBR, relative differenced NBR, Relativized Burn Ratio) and seven Sentinel 2A-specific indices considering the availability of spectral information recorded in the red-edge spectral region. The statistical correlation indicated a slightly stronger relationship between the differenced NBR and the GeoCBI for both Sentinel 2A (r = 0.872) and Landsat-8 OLI (r = 0.845) imagery. Predictive local thresholds of dNBR values showed slightly higher classification accuracy for Sentinel 2A (73.33%) than Landsat-8 OLI (71.11%), suggesting the adequacy of Sentinel 2A for forest fire severity assessment and mapping in Mediterranean pine ecosystems. The evaluation of the classification thresholds calculated in this study over other fires with similar pre-fire conditions could contribute in the operational mapping and reconstruction of the historical patterns of fire severity over the Eastern Mediterranean region.     

Developing water quality retrieval models with in situ hyperspectral data in Poyang Lake,China

Applying remote sensing techniques to develop the retrieval models and further to obtain the spatiotemporal information of water quality parameters is necessary for understanding, managing, and protecting lake ecosystems. This study aimed to calibrate and validate the retrieval models for estimating the concentrations of chlorophyll a (C_CHL), suspended particulate matter (C_SPM), and dissolved organic carbon (C_DOC) with the in situ hyperspectral measurements in Poyang Lake, China in 2010 and 2011. The model calibration and validation results indicated that: (1) for C_CHL retrieval, significantly strong and moderate correlations existed between the measured and estimated values (with the correlation coefficient r = 0.92 and r = 0.76) using the exponential model and the three-band model, respectively, with biased estimation observed for the exponential model; (2) for retrieving C_SPM, there was a strong correlation between the measured and estimated values (r = 0.95) using the exponential model; and (3) no significant correlation between measured and estimated C_DOC values was found with our developed models. More work is needed to allow the water quality of Poyang Lake to be accurately and steadily estimated, especially for C_CHL and C_DOC.     

Analysis of Seasonal Variability of Vegetation and Methane Concentration over India using SPOT-VEGETATION and ENVISAT-SCIAMACHY Data

This paper highlights the spatial and temporal variability of atmospheric columnar methane (CH₄) concentration over India and its correlation with the terrestrial vegetation dynamics. SCanning IMaging Absorption spectrometer for Atmospheric CHartographY (SCIAMACHY) on board ENVIronmental SATellite (ENVISAT) data product (0.5° × 0.5°) was used to analyze the atmospheric CH₄ concentration. Satellite Pour l'Observation de la Terre (SPOT)-VEGETATION sensor’s Normalized Difference Vegetation Index (NDVI) product, aggregated at 0.5° × 0.5° grid level for the same period (2004 and 2005), was used to correlate the with CH₄ concentration. Analysis showed mean monthly CH₄ concentration during the Kharif season varied from 1,704 parts per billion volume (ppbv) to 1,780 ppbv with the lowest value in May and the highest value in September. Correspondingly, mean NDVI varied from 0.28 (May) to 0.53 (September). Analysis of correlation between CH₄ concentration and NDVI values over India showed positive correlation (r = 0.76; n = 6) in Kharif season. Further analysis using land cover information showed characteristic low correlation in natural vegetation region and high correlation in agricultural area. Grids, particularly falling in the Indo-Gangetic Plains showed positive correlation. This could be attributed to the rice crop which is grown as a predominant crop during this period. The CH₄ concentration pattern matched well with growth pattern of rice with the highest concentration coinciding with the peak growth period of crop in the September. Characteristically low correlation was observed (r = 0.1; n = 6) in deserts of Rajasthan and forested Himalayan ecosystem. Thus, the paper emphasizes the synergistic use of different satellite based data in understanding the variability of atmospheric CH₄ concentration in relation to vegetation.     

Yield Estimation Model and Water Productivity of Wheat Crop (Triticum aestivum) in an Irrigation Command Using Remote Sensing and GIS

Crop yield estimation has an important role on economy development and its accuracy and speed influence yield price and helps in deciding the excess or deficit production conditions. The water productivity evaluates the irrigation command through water use efficiency (WUE). Remote sensing (RS) and geographical information system (GIS) techniques were used for crop yield and water productivity estimation of wheat crop (Triticum aestivum) grown in Tarafeni South Main Canal (TSMC) irrigation command of West Bengal State in India. One IRS P6 image and four wide field sensor (WiFS) images for different months of winter season were used to determine the Normalized Difference Vegetation Index (NDVI) and Soil Adjusted Vegetation Index (SAVI) for area under wheat crop. The temporally and spatially distributed spectral growth profile and AREASUM of NDVI (A_NDVI) and SAVI (A_SAVI) with time after sowing of wheat crop were developed and correlated with actual crop yield of wheat (Y_act). The developed relationships between A_SAVI and Y_act resulted high correlation in comparison to that of A_NDVI. Using the developed model the RS based wheat yield (Y_RS) predicted from A_SAVI varied on entire TSMC irrigation command from 22.67 to 33.13 q ha⁻¹ respectively, which gave an average yield of 26.50 q ha⁻¹. The RS generated yield based water use efficiency (WUE_YRS) for water supplied from canal of TSMC irrigation command was found to be 6.69 kg ha⁻¹ mm⁻¹.     

Quantification of LST and CO2 levels using Landsat-8 thermal bands on urban environment

Land surface temperature (LST) is an important aspect in global to regional change studies, for control of climate change and balancing of high temperature. Urbanization is one of the influencing factors increasing land surface and atmospheric temperature, by the emission of greenhouse gases (e.g. CO₂, NO and methane). In the present study, LST was derived from Landsat-8 of multitemporal data sets to analyse the spatial structure of the urban thermal environment in relation to the urban surface characteristics and land use–land cover (LULC). LST is influenced by the greenhouse gases i.e. CO₂ plays an important role in increasing the earth’s surface temperature. In order to provide the evidence of influence of CO₂ on LST, the relationship between LST, air temperature and CO₂ was analysed. Landsat-8 satellite has two thermal bands, 10 and 11. These bands were used to accurately to calculate the temperature over the study area. Results showed that the strength of correlation between ground monitoring data and satellite data was high. Based on correlation values of each month April (R² = 0.994), May (R² = 0.297) and June (R² = 0.934), observed results show that band 10 was significantly correlating with air temperature. Relationship between LST and CO₂ levels were obtained from linear regression analysis. band 11 was correlating significantly with CO₂ values in each of the months April (R² = 0.217), May (R² = 0.914) and June, (R² = 0.934), because band 11 is closer to the 15-micron band of CO₂. From the results, it was observed that band 10 can be used for calculating air temperature and band 11 can be used for estimation of greenhouse gases.     

Malaria incidence in Nairobi,Kenya and dekadal trends in NDVI and climatic variables

The primary objective of this research was to determine if the remotely-sensed metric, Normalised Difference Vegetation Index (NDVI) and ground-collected dekadal climatological variables were useful predictors of future malaria outbreaks in an epidemic-prone area of Nairobi, Kenya. Data collected consisted of 36 dekadal (10-day) periods for the variables rainfall, temperature and NDVI along with yearly documented malaria admissions in 2003 for Nairobi, Kenya. Linear regression models were built for malaria cases reported in Nairobi, Kenya, as the dependent variable and various time-based groupings of temperature, rainfall and NDVI data from the dekads in both the current and the previous month as the independent variables. Data from 2003 show that malaria incidence in any given month is best predicted (R²  = 0.881, p < 0.001) by the average NDVI for the 30 days including the final two dekads of the previous month and first dekad of the current month, and by the average rainfall for the 30 days including the three dekads of rainfall data from the prior month. Forecasting an outbreak in an epidemic zone would allow public health entities to plan for and disseminate resources to the general public such as antimalarials and insecticide impregnated bed nets. In addition, vector control measures could be implemented to slow the rate of transmission in the impacted population.     

Investigating the vegetation and agricultural responses to El Nino/Southern Oscillation using AVHRR data

This study explores the possible linkages of El Nino/Southern Oscillation (ENSO) with vegetation and rainfall patterns, vegetation activity and food grain yields, in arid and semi-arid regions of western India. A sequence of 20-year (1981–2000) monthly maximum Normalized Difference Vegetation Index (NDVI) data from the Advanced Very High Resolution Radiometer (AVHRR) and monthly rainfall from 160 stations were examined to study the seasonal patterns and their relation to ENSO activity. In addition, a direct (ENSO-crop yield) linkage and an intermediate (ENSO-NDVI) linkage of agricultural responses to ENSO were also investigated. The results indicate below-normal seasonal NDVI and rainfall associated with El Nino (warm) events, except during 1997, while positive anomalies occur during La Nina (cold) events. Sea surface temperature (SST) anomalies from NINO 3 region (5°N–5°S; 150°W–90°W), as an indicator of ENSO were significantly correlated with NDVI anomalies, rainfall anomalies and yield anomalies but the Southern Oscillation Index (SOI) was significantly related to NDVI anomalies only. NDVI anomaly patterns correspond to rainfall variability including that associated with ENSO activity. The observed strong intermediate linkage between yield anomalies and NDVI anomaly signal (r = 0.609) indicates that NDVI is an ideal index for understanding and analysing agricultural response to ENSO climate teleconnections.     

Winter cover crops in Dutch maize fields: Variability in quality and its drivers assessed from multi-temporal Sentinel-2 imagery

Planting a cover crop between the main cropping seasons is an agricultural management measure with multiple potential benefits for sustainable food production. In the maize production system of the Netherlands, an effective establishment of a winter cover crop is important for reducing nitrogen leaching to groundwater. Cover crop establishment after maize cultivation is obliged by law for sandy soils and consequently implemented on nearly all maize fields, but the winter-time vegetative ground cover varies significantly between fields. The objectives of this study are to assess the variability in winter vegetative cover and evaluate to what extent this variability can be explained by the timing of cover crop establishment and weather conditions in two growing seasons (2017–2018). We used Sentinel-2 satellite imagery to construct NDVI time series for fields known to be cultivated with maize within the province of Overijssel. We fitted piecewise logistic functions to the time series in order to estimate cover crop sowing date and retrieve the fitted NDVI value for 1 December (NDVI_Dec). We used NDVI_Dec to represent the quality of cover crop establishment at the start of the winter season. The Sentinel-2 estimated sowing dates compared reasonably with ground reference data for eight fields (RMSE = 6.6 days). The two analysed years differed considerably, with 2018 being much drier and warmer during summer. This drought resulted in an earlier estimated cover crop sowing date (on average 19 days) and an NDVI_Dec value that was 0.2 higher than in 2017. Combining both years and all fields, we found that Sentinel-2 retrieved sowing dates could explain 55% of the NDVI_Dec variability. This corresponded to a positive relationship (R² = 0.50) between NDVI_Dec and the cumulative growing degree days (GDD) between sowing date and 1 December until reaching 400 GDD. Based on cumulative GDD derived from two weather stations within Overijssel, we found that on average for the past three decades a sowing date of 19 September (± 7 days) allowed to attain these 400 GDD; this provides support for the current legislation that states that from 2019 onwards a cover crop should be sown before 1 October. To meet this deadline, while simultaneously ascertaining a harvest-ready main crop, in practice implies that undersowing of the cover crop during spring will gain importance. Our results show that Sentinel-2 NDVI time series can assess the effectiveness and timing of cover crop growth for small agricultural fields, and as such has potential to inform regulatory frameworks as well as farmers with actionable information that may help to reduce nitrogen leaching.     

Estimation of LAI and above-ground biomass in deciduous forests: Western Ghats of Karnataka, India

This study demonstrates the potentials of IRS P6 LISS-IV high-resolution multispectral sensor (IGFOV  6 m)-based estimation of biomass in the deciduous forests in the Western Ghats of Karnataka, India. Regression equations describing the relationship between IRS P6 LISS-IV data-based vegetation index (NDVI) and field measured leaf area index (ELAI) and estimated above-ground biomass (EAGB) were derived. Remote sensing (RS) data-based leaf area index (PLAI) image is generated using regression equation based on NDVI and ELAI (r² = 0.68, p ≤ 0.05). RS-based above-ground biomass (PAGB) image was generated based on regression equation developed between PLAI and EAGB (r² = 0.63, p ≤ 0.05). The mean value of estimated above-ground biomass and RS-based above-ground biomass in the study area are 280(±72.5) and 297.6(±55.2) Mg ha⁻¹, respectively. The regression models generated in the study between NDVI and LAI; LAI and biomass can also help in generating spatial biomass map using RS data alone. LISS-IV-based estimation of biophysical parameters can also be used for the validation of various coarse resolution satellite products derived from the ground-based measurements alone.     

Estimating vegetation parameter for soil erosion assessment in an alpine catchment by means of QuickBird imagery

Soil erosion rates in alpine regions are related to high spatial variability complicating assessment of risk and damages. A crucial parameter triggering soil erosion that can be derived from satellite imagery is fractional vegetation cover (FVC). The objective of this study is to assess the applicability of normalized differenced vegetation index (NDVI), linear spectral unmixing (LSU) and mixture tuned matched filtering (MTMF) in estimating abundance of vegetation cover in alpine terrain. To account for the small scale heterogeneity of the alpine landscape we used high resolved multispectral QuickBird imagery (pixel resolution = 2.4 m) of a site in the Urseren Valley, Central Swiss Alps (67 km²). A supervised land-cover classification was applied (total accuracy 93.3%) prior to the analysis in order to stratify the image. The regression between ground truth FVC assessment and NDVI as well as MTMF-derived vegetation abundance was significant (r² = 0.64, r² = 0.71, respectively). Best results were achieved for LSU (r² = 0.85). For both spectral unmixing approaches failed to estimate bare soil abundance (r² = 0.39 for LSU, r² = 0.28 for MTMF) due to the high spectral variability of bare soil at the study site and the low spectral resolution of the QuickBird imagery. The LSU-derived FVC map successfully identified erosion features (e.g. landslides) and areas prone to soil erosion. FVC represents an important but often neglected parameter for soil erosion risk assessment in alpine grasslands.     

Generation of high spatial and temporal resolution NDVI and its application in crop biomass estimation

Abstract

While data like HJ-1 CCD images have advantageous spatial characteristics for describing crop properties, the temporal resolution of the data is rather low, which can be easily made worse by cloud contamination. In contrast, although Moderate Resolution Imaging Spectroradiometer (MODIS) can only achieve a spatial resolution of 250 m in its normalised difference vegetation index (NDVI) product, it has a high temporal resolution, covering the Earth up to multiple times per day. To combine the high spatial resolution and high temporal resolution of different data sources, a new method (Spatial and Temporal Adaptive Vegetation index Fusion Model [STAVFM]) for blending NDVI of different spatial and temporal resolutions to produce high spatial–temporal resolution NDVI datasets was developed based on Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). STAVFM defines a time window according to the temporal variation of crops, takes crop phenophase into consideration and improves the temporal weighting algorithm. The result showed that the new method can combine the temporal information of MODIS NDVI and spatial difference information of HJ-1 CCD NDVI to generate an NDVI dataset with both high spatial and high temporal resolution. An application of the generated NDVI dataset in crop biomass estimation was provided. An average absolute error of 17.2% was achieved. The estimated winter wheat biomass correlated well with observed biomass (R ² of 0.876). We conclude that the new dataset will improve the application of crop biomass estimation by describing the crop biomass accumulation in detail. There is potential to apply the approach in many other studies, including crop production estimation, crop growth monitoring and agricultural ecosystem carbon cycle research, which will contribute to the implementation of Digital Earth by describing land surface processes in detail.     

A GPP assimilation model for the southeastern Tibetan Plateau based on CO2 eddy covariance flux tower and remote sensing data

The gross primary production (GPP) at individual CO₂ eddy covariance flux tower sites (GPP_Tower) in Dali (DL), Wenjiang (WJ) and Linzhi (LZ) around the southeastern Tibetan Plateau were determined by the net ecosystem exchange of CO₂ (NEE) and ecosystem respiration (R_e). The satellite remote sensing-VPM model estimates of GPP values (GPP_MODIS) used the satellite-derived 8-day surface reflectance product (MOD09A1), including satellite-derived enhanced vegetation index (EVI) and land surface water index (LSWI). In this paper, we assembled a subset of flux tower data at these three sites to calibrate and test satellite-VPM model estimated GPP_MODIS, and introduced the satellite data and site-level environmental factors to develop four new assimilation models. The new assimilation models’ estimates of GPP values were compared with GPP_MODIS and GPP_Tower, and the final optimum model among the four assimilation models was determined and used to calibrate GPP_MODIS. The results showed that GPP_MODIS had similar temporal variations to the GPP_Tower, but GPP_MODlS were commonly higher in absolute magnitude than GPP_Tower with relative error (RE) about 58.85%. While, the assimilation models’ estimates of GPP values (GPP_MODEL) were much more closer to GPP_Tower with RE approximately 6.98%, indicating that the capacity of the simulation in the new assimilation model was greatly improved, the R² and root mean square error (RMSE) of the new assimilation model were 0.57–4.90% higher and 0.74–2.47 g C m⁻² s⁻¹ lower than those of the GPP_MODIS, respectively. The assimilation model was used to predicted GPP dynamics around the Tibetan Plateau and showed a reliable result compared with other researches. This study demonstrated the potential of the new assimilation model for estimating GPP around the Tibetan Plateau and the performances of site-level biophysical parameters in related to satellite-VPM model GPP.     

Assessment and validation of evapotranspiration using SEBAL algorithm and Lysimeter data of IARI agricultural farm,India

Evapotranspiration (ET) is a vital process in land surface atmosphere research. In this study, Surface Energy Balance Algorithm for Land (SEBAL) for the assessment of ET (for 23 December 2010, 8 January 2011, 24 January 2011, 9 February 2011, 25 February 2011, 29 March 2011 and 14 April 2011) from LANDSAT7-ETM+ and validation with Lysimeter data set is illustrated. It is based on the evaporative fraction concept, and it has been applied to LANDSAT7-ETM + (30 m resolution) data acquired over the Indian Agricultural Research Institute’s agricultural farm land. The ET from SEBAL was compared with Lysimeter ET using four statistical tests (root-mean-square error (RMSE), relative root-mean-square error (R-RMSE), mean absolute error (MAE), and normalized root-mean square error (NRMSE)), and each test showed a good correlation between the predicted and observed ET values. Results from this study revealed that the RMSE of crop-growing period was 0.51 mm d^?1 for ET_SEBAL, i.e. ET_SEBAL having good accuracy with respect to observed ET_Lysimeter. Results were also validated using R-RMSE test, which also proved that ET_SEBAL data are having good accuracy with respect to observed ET_Lysimeter as R-RMSE of crop-growing period is 0.19 mm d^?1. MAE (0.19), NRMSE (0.21) and r² (0.91) tests indicated that model prediction is significant, and model can be effectively used for the estimation of ET from SEBAL as input of remote sensing data sets. Finally, the SEBAL has been useful for remote agricultural land where ground-based data (Lysimeter data) are not available for daily ET (ET_24 h) estimation. The temporal study of the ET_24 h values analysed has revealed that the highest ET_24 h values are owing to the higher development (high greenness) of crops, whereas the lower values are related to the lower development (low greenness) or null crop.