Real-time retrieval of Leaf Area Index from MODIS time series data

Real-time/near real-time inversion of land surface biogeophysical variables from satellite observations is required to monitor rapid land surface changes, and provide the necessary input for numerical weather forecasting models and decision support systems. This paper develops a new inversion method for the real-time estimation of the Leaf Area Index (LAI) of land surfaces from MODIS time series reflectance data (MOD09A1). It consists of a series of procedures, including time series data smoothing, data quality control and real-time estimation of LAI. After the historical LAI time series is smoothed by a multi-step Savitzky-Golay filter to determine the upper LAI envelope, a Seasonal Auto-Regressive Integrated Moving Average (SARIMA) model is used to derive the LAI climatology. Based on the climatology from the SARIMA model to evolve LAI in time, a dynamic model is then constructed and used to provide the short-range forecast of LAI. Predictions from this model are used with Ensemble Kalman Filter (EnKF) techniques to recursively update biophysical variables as new observations arrive. The validation results produced using MODIS surface reflectance data and field-measured LAI data at eight BELMANIP sites show that the real-time inversion method is able to efficiently produce a relatively smooth LAI product. In addition, the accuracy is significantly improved over the MODIS LAI product.     

Impact of a satellite-derived leaf area index monthly climatology in a global numerical weather prediction model

The leaf area index (LAI), defined as the one-sided green leaf area per unit ground area, is used in many numerical weather prediction (NWP) models as an indicator of the vegetation development state, which is of paramount importance to characterize land evaporation, photosynthesis, and carbon-uptake processes. LAI is often simply represented by lookup tables, dependent on the vegetation type and seasons. However, global LAI datasets derived from remote sensing observations have more recently become available. These products are based on sensors such as the Advanced Very High Resolution Radiometer (AVHRR) or the Moderate Resolution Imaging Spectroradiometer (MODIS), onboard polar orbiting satellites that can cover the entire globe within typically 3 days and with a spatial resolution of the order of 1 km.

We examine the meteorological impact of satellite-derived LAI products on near-surface air temperature and humidity, which comes both from the stomatal transpiration of leaves and from the intercepted water on the surface of leaves, re-evaporating into the atmosphere.

Two distinct monthly LAI climatology datasets derived respectively from AVHRR and MODIS sensors are tested. A set of forecasts and data assimilation experiments with the integrated forecasting system of the European Centre for Medium-range Weather Forecasts is performed with the monthly LAI climatology datasets as opposed to a vegetation-dependent constant LAI. The monthly LAI is shown to improve the forecasts of near-surface (screen-level) air temperature and relative humidity through its effect on evapotranspiration, with the largest impact obtained over needleleaf forests, crops, and grassland. At longer time-scales, the introduction of the monthly LAI is shown to have a positive impact on the model climate particularly during the boreal spring, where the LAI climatology has a large seasonal cycle.     

Assessment of the MODIS LAI product for Australian ecosystems

The leaf area index (LAI) product from the Moderate Resolution Imaging Spectroradiometer (MODIS) is important for monitoring and modelling global change and terrestrial dynamics at many scales. The algorithm relies on spectral reflectances and a six biome land cover classification. Evaluation of the specific behaviour and performance of the product for regions of the globe such as Australia are needed to assist with product refinement and validation. We made an assessment of Collection 4 of the MODIS LAI product using four approaches: (a) assessment against a continental scale Structural Classification of Australian Vegetation (SCAV); (b) assessment against a continental scale land use classification (LUC); (c) assessment against historical field-based measurement of LAI collected prior to the Terra Mission; and (d) direct comparison of MODIS LAI with coincident field measurements of LAI, mostly from hemispherical photography. The MODIS LAI product produced a wide variety of geographically and structurally specific temporal response profiles between different classes and even for sub-groups within classes of the SCAV. Historical and concurrent field measurements indicated that MODIS LAI was giving reasonable estimates for LAI for most cover types and land use types, but that major overestimation of LAI occurs in some eastern Australian open forests and woodlands. The six biome structural land cover classification showed some significant deviations in class allocation compared to the SCAV particularly where grasslands are allocated to shrubland, savanna woodlands are allocated to shrubland, savanna and broadleaf forest, and open forests are allocated to savanna and broadleaf forest. The land cover and LAI products could benefit from some additional examination of Australian data addressing the structural representation of Eucalypt canopies in the “space of canopy realisation” for savanna and broadleaf forest classes.     

Correcting MM estimates for “fat” data sets

Regression MM estimates require the estimation of the error scale, and the determination of a constant that controls the efficiency. These two steps are based on the asymptotic results that are derived assuming that the number of predictors p remains fixed while the number of observations n tends to infinity, which means assuming that the ratio p/n is “small”. However, many high-dimensional data sets have a “large” value of p/n (say, ≥0.2). It is shown that the standard asymptotic results do not hold if p/n is large; namely that (a) the estimated scale underestimates the true error scale, and (b) that even if the scale is correctly estimated, the actual efficiency can be much lower than the nominal one. To overcome these drawbacks simple corrections for the scale and for the efficiency controlling constant are proposed, and it is demonstrated that these corrections improve on the estimate’s performance under both normal and contaminated data.     

Reprocessing the MODIS Leaf Area Index products for land surface and climate modelling

Land surface and climate modelling requires continuous and consistent Leaf Area Index (LAI). High spatiotemporal resolution and long-time record data are more in demand nowadays and will continue to be in the future. MODIS LAI products meet these requirements to some degree. However, due to the presence of cloud and seasonal snow cover, the instrument problems and the uncertainties of retrieval algorithm, the current MODIS LAI products are spatially and temporally discontinuous and inconsistent, which limits their application in land surface and climate modelling. To improve the MODIS LAI products on a global scale, we considered the characteristics of the MODIS LAI data and made the best use of quality control (QC) information, and developed an integrated two-step method to derive the improved MODIS LAI products effectively and efficiently on a global scale. First, we used the modified temporal spatial filter (mTSF) method taking advantage of background values and QC information at each pixel to do a simple data assimilation for relatively low quality data. Then we applied the post processing-TIMESAT (A software package to analyze time-series of satellite sensor data) Savitzky-Golay (SG) filter to get the final result. We implemented the method to 10 years of the MODIS Collection 5 LAI data. In comparison with the LAI reference maps and the MODIS LAI data, our results showed that the improved MODIS LAI data are closer to the LAI reference maps in magnitude and also more continuous and consistent in both time-series and spatial domains. In addition, simple statistics were used to evaluate the differences between the MODIS LAI and the improved MODIS LAI.     

Flow pattern around an appended tanker hull form in simple maneuvering conditions

A detailed computational study is presented of the flow pattern around the Esso Osaka with rudder in simple maneuvering conditions: “static rudder” and “pure drift”. The objectives are: (1) apply RANS for maneuvering simulation; (2) perform verification and validation on field quantities; (3) characterize flow pattern; and (4) correlate behavior of the integral quantities with the flow field. The general-purpose code CFDSHIP-IOWA is used. The free surface is neglected and the two-equation k-ω turbulence model is used. The levels of verification of the velocity components for the “straight-ahead”, “static rudder” and “pure drift” conditions show ranges from 5.5% to 28.3% of free stream, U₀, for the axial velocity U and 2.5-29.1%U₀ for the cross flow (V, W). Qualitative validation against limited experimental data shows encouraging results with respect to trends and levels. The flow pattern is characterized by fore and aft body bilge and side vortices, which are similar for “straight-ahead” and “static rudder” conditions, except in close vicinity of the rudder. The “pure drift” condition shows strong asymmetry on windward vs. leeward sides and a more complex vortex system with additional bilge vortices. Similarities and differences with data for other tanker, container, and surface combatant hulls and relation between flow pattern and forces and moments are discussed. Future work focuses on influence of propeller.     

On the relationship of NDVI with leaf area index in a deciduous forest site

Numerous studies have reported on the relationship between the normalized difference vegetation index (NDVI) and leaf area index (LAI), but the seasonal and annual variability of this relationship has been less explored. This paper reports a study of the NDVI-LAI relationship through the years from 1996 to 2001 at a deciduous forest site. Six years of LAI patterns from the forest were estimated using a radiative transfer model with input of above and below canopy measurements of global radiation, while NDVI data sets were retrieved from composite NDVI time series of various remote sensing sources, namely NOAA Advanced Very High Resolution Radiometer (AVHRR; 1996, 1997, 1998 and 2000), SPOT VEGETATION (1998-2001), and Terra MODIS (2001). Composite NDVI was first used to remove the residual noise based on an adjusted Fourier transform and to obtain the NDVI time-series for each day during each year.The results suggest that the NDVI-LAI relationship can vary both seasonally and inter-annually in tune with the variations in phenological development of the trees and in response to temporal variations of environmental conditions. Strong linear relationships are obtained during the leaf production and leaf senescence periods for all years, but the relationship is poor during periods of maximum LAI, apparently due to the saturation of NDVI at high values of LAI. The NDVI-LAI relationship was found to be poor (R² varied from 0.39 to 0.46 for different sources of NDVI) when all the data were pooled across the years, apparently due to different leaf area development patterns in the different years. The relationship is also affected by background NDVI, but this could be minimized by applying relative NDVI.Comparisons between AVHRR and VEGETATION NDVI revealed that these two had good linear relationships (R²=0.74 for 1998 and 0.63 for 2000). However, VEGETATION NDVI data series had some unreasonably high values during beginning and end of each year period, which must be discarded before adjusted Fourier transform processing. MODIS NDVI had values greater than 0.62 through the entire year in 2001, however, MODIS NDVI still showed an “M-shaped” pattern as observed for VEGETATION NDVI in 2001. MODIS enhanced vegetation index (EVI) was the only index that exhibited a poor linear relationship with LAI during the leaf senescence period in year 2001. The results suggest that a relationship established between the LAI and NDVI in a particular year may not be applicable in other years, so attention must be paid to the temporal scale when applying NDVI-LAI relationships.     

Retrieving soil temperature profile by assimilating MODIS LST products with ensemble Kalman filter

Proper estimation of initial state variables and model parameters are vital importance for determining the accuracy of numerical model prediction. In this work, we develop a one-dimensional land data assimilation scheme based on ensemble Kalman filter and Common Land Model version 3.0 (CoLM). This scheme is used to improve the estimation of soil temperature profile. The leaf area index (LAI) is also updated dynamically by MODIS LAI production and the MODIS land surface temperature (LST) products are assimilated into CoLM. The scheme was tested and validated by observations from four automatic weather stations (BTS, DRS, MGS, and DGS) in Mongolian Reference Site of CEOP during the period of October 1, 2002 to September 30, 2003. Results indicate that data assimilation improves the estimation of soil temperature profile about 1 K. In comparison with simulation, the assimilation results of soil heat fluxes also have much improvement about 13 W m^− 2 at BTS and DGS and 2 W m^− 2 at DRS and MGS, respectively. In addition, assimilation of MODIS land products into land surface model is a practical and effective way to improve the estimation of land surface variables and fluxes.     

Solving the minimum sum-of-squares clustering problem by hyperbolic smoothing and partition into boundary and gravitational regions

This article considers the minimum sum-of-squares clustering (MSSC) problem. The mathematical modeling of this problem leads to a min-sum-min formulation which, in addition to its intrinsic bi-level nature, has the significant characteristic of being strongly nondifferentiable. To overcome these difficulties, the proposed resolution method, called hyperbolic smoothing, adopts a smoothing strategy using a special C^∞ differentiable class function. The final solution is obtained by solving a sequence of low dimension differentiable unconstrained optimization subproblems which gradually approach the original problem. This paper introduces the method of partition of the set of observations into two nonoverlapping groups: “data in frontier” and “data in gravitational regions”. The resulting combination of the two methodologies for the MSSC problem has interesting properties, which drastically simplify the computational tasks.     

植被MODIS-LAI的温度降水响应

本研究首先利用贵州省2001全年基于MODIS遥感影像反演的8天一景的LAI平均资料,采用高光谱遥感数据处理分析方法对研究区地物进行分类,得到研究区植被覆盖类型结果;然后基于此分类结果,分别选取研究区19个气象站点附近各种类型植被,结合2001年贵州省温度和降水气象数据,建立了各种植被LAI与温度及降水的相关关系,用以揭示植被的季节变化对气候变化的敏感性特征响应,并从地球物理学、生态学、气候学等多个角度进一步阐明了出现相关特征的原因和物理机制。结果表明:研究区植被覆盖与温度随季节变化有显著的相关,而与降水的相关则由于植被对水分的滞后性反映则不甚显著;其中常绿林和落叶林由于较为稳固的生态环境和较少受到一些人为因素干扰而与温度和降水具有更高的相关性,而农田和草地对温度和降水的响应则更为复杂和敏感。     

Leaf area index estimation with MODIS reflectance time series and model inversion during full rotations of Eucalyptus plantations

The leaf area index (LAI) of fast-growing Eucalyptus plantations is highly dynamic both seasonally and inter-annually, and is spatially variable depending on pedo-climatic conditions. LAI is very important in determining the carbon and water balance of a stand, but is difficult to measure during a complete stand rotation and at large scales. Remote-sensing methods allowing the retrieval of LAI time series with accuracy and precision are therefore necessary. Here, we tested two methods for LAI estimation from MODIS 250m resolution red and near-infrared (NIR) reflectance time series. The first method involved the inversion of a coupled model of leaf reflectance and transmittance (PROSPECT4), soil reflectance (SOILSPECT) and canopy radiative transfer (4SAIL2). Model parameters other than the LAI were either fixed to measured constant values, or allowed to vary seasonally and/or with stand age according to trends observed in field measurements. The LAI was assumed to vary throughout the rotation following a series of alternately increasing and decreasing sigmoid curves. The parameters of each sigmoid curve that allowed the best fit of simulated canopy reflectance to MODIS red and NIR reflectance data were obtained by minimization techniques. The second method was based on a linear relationship between the LAI and values of the GEneralized Soil Adjusted Vegetation Index (GESAVI), which was calibrated using destructive LAI measurements made at two seasons, on Eucalyptus stands of different ages and productivity levels. The ability of each approach to reproduce field-measured LAI values was assessed, and uncertainty on results and parameter sensitivities were examined. Both methods offered a good fit between measured and estimated LAI (R² = 0.80 and R² = 0.62 for model inversion and GESAVI-based methods, respectively), but the GESAVI-based method overestimated the LAI at young ages.     

Monitoring boreal forest leaf area index across a Siberian burn chronosequence: a MODIS validation study

Landscapes containing differing amounts of ecological disturbance provide an excellent opportunity to validate and better understand the emerging Moderate Resolution Imaging Spectrometer (MODIS) vegetation products. Four sites, including 1‐year post‐fire coniferous, 13‐year post‐fire deciduous, 24‐year post‐fire deciduous, and >100 year old post‐fire coniferous forests, were selected to serve as a post‐fire chronosequence in the central Siberian region of Krasnoyarsk (57.3°N, 91.6°E) with which to study the MODIS leaf area index (LAI) and vegetation index (VI) products. The collection 4 MODIS LAI product correctly represented the summer site phenologies, but significantly underestimated the LAI value of the >100 year old coniferous forest during the November to April time period. Landsat 7‐derived enhanced vegetation index (EVI) performed better than normalized difference vegetation index (NDVI) to separate the deciduous and conifer forests, and both indices contained significant correlation with field‐derived LAI values at coniferous forest sites (r ² = 0.61 and r ² = 0.69, respectively). The reduced simple ratio (RSR) markedly improved LAI prediction from satellite measurements (r ² = 0.89) relative to NDVI and EVI. LAI estimates derived from ETM+ images were scaled up to evaluate the 1 km resolution MODIS LAI product; from this analysis MODIS LAI overestimated values in the low LAI deciduous forests (where LAI<5) and underestimated values in the high LAI conifer forests (where LAI>6). Our results indicate that further research on the MODIS LAI product is warranted to better understand and improve remote LAI quantification in disturbed forest landscapes over the course of the year.     

Combining vegetation index and model inversion methods for the extraction of key vegetation biophysical parameters using Terra and Aqua MODIS reflectance data

Accurate estimates of vegetation biophysical variables are valuable as input to models describing the exchange of carbon dioxide and energy between the land surface and the atmosphere and important for a wide range of applications related to vegetation monitoring, weather prediction, and climate change. The present study explores the benefits of combining vegetation index and physically based approaches for the spatial and temporal mapping of green leaf area index (LAI), total chlorophyll content (TC_ab), and total vegetation water content (VWC). A numerical optimization method was employed for the inversion of a canopy reflectance model using Terra and Aqua MODIS multi-spectral, multi-temporal, and multi-angle reflectance observations to aid the determination of vegetation-specific physiological and structural canopy parameters. Land cover and site-specific inversion modeling was applied to a restricted number of pixels to build multiple species- and environmentally dependent formulations relating the three biophysical properties of interest to a number of selected simpler spectral vegetation indices (VI). While inversions generally are computationally slow, the coupling with the simple and computationally efficient VI approach makes the combined retrieval scheme for LAI, TC_ab, and VWC suitable for large-scale mapping operations. In order to facilitate application of the canopy reflectance model to heterogeneous forested areas, a simple correction scheme was elaborated, which was found to improve forest LAI predictions significantly and also provided more realistic values of leaf chlorophyll contents.The inversion scheme was designed to enable biophysical parameter retrievals for land cover classes characterized by contrasting canopy architectures, leaf inclination angles, and leaf biochemical constituents without utilizing calibration measurements. Preliminary LAI validation results for the Island of Zealand, Denmark (57°N, 12°E) provided confidence in the approach with root mean square (RMS) deviations between estimates and in-situ measurements of 0.62, 0.46, and 0.63 for barley, wheat, and deciduous forest sites, respectively. Despite the independence on site-specific in-situ measurements, the RMS deviations of the automated approach are in the same range as those established in other studies employing field-based empirical calibration.Being completely automated and image-based and independent on extensive and impractical surface measurements, the retrieval scheme has potential for operational use and can quite easily be implemented for other regions. More validation studies are needed to evaluate the usefulness and limitations of the approach for other environments and species compositions.     

Variational retrieval of leaf area index from MODIS time series data: examples from the Heihe river basin,north-west China

Leaf area index (LAI) products retrieved from observations acquired on one occasion have obvious discontinuity in the time series owing to cloud coverage and other factors, and the accuracy may not meet the needs of many applications. Effectively utilizing data assimilation techniques to retrieve biophysical parameters from the time series of remote-sensing data has attracted special interest. The data assimilation technique is based on a reasonable consideration of dynamic change rules of biophysical parameters and time series observational quantities, thereby improving the quality of retrieved profiles. In this article, a variational assimilation procedure for retrieving LAI from the time series of remote-sensing data is developed. The procedure is based on the formulation of an objective function. A dynamic model is constructed based on the climatology from multi-year Moderate Resolution Imaging Spectroradiometer (MODIS) LAI data to evolve LAI in time, and a radiative transfer model is coupled with the dynamic model to simulate a time series of surface reflectances. A shuffled complex evolution method (developed at the University of Arizona; SCE-UA) optimization algorithm is then used to minimize the objective function and estimate the dynamic model states and the parameters of the coupled model from the MODIS reflectance data with a higher quality in a given time window. The variational assimilation method is applied to the MODIS surface reflectance data for the whole of 2008 at the Heihe river basin to produce regional LAI mapping results. The ground LAI data measured in situ are used to develop algorithms to estimate LAI from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) surface reflectance, and ASTER LAI maps are produced for each ASTER scene using the algorithms developed. Then the ASTER LAI maps are aggregated to compare with the new LAI products. It is found that the variational assimilation method is able to produce temporal continuous LAI products and that accuracy has been improved over the MODIS LAI standard product.     

Penalized spline estimation for functional coefficient regression models

The functional coefficient regression models assume that the regression coefficients vary with some “threshold” variable, providing appreciable flexibility in capturing the underlying dynamics in data and avoiding the so-called “curse of dimensionality” in multivariate nonparametric estimation. We first investigate the estimation, inference, and forecasting for the functional coefficient regression models with dependent observations via penalized splines. The P-spline approach, as a direct ridge regression shrinkage type global smoothing method, is computationally efficient and stable. With established fixed-knot asymptotics, inference is readily available. Exact inference can be obtained for fixed smoothing parameter λ, which is most appealing for finite samples. Our penalized spline approach gives an explicit model expression, which also enables multi-step-ahead forecasting via simulations. Furthermore, we examine different methods of choosing the important smoothing parameter λ: modified multi-fold cross-validation (MCV), generalized cross-validation (GCV), and an extension of empirical bias bandwidth selection (EBBS) to P-splines. In addition, we implement smoothing parameter selection using mixed model framework through restricted maximum likelihood (REML) for P-spline functional coefficient regression models with independent observations. The P-spline approach also easily allows different smoothness for different functional coefficients, which is enabled by assigning different penalty λ accordingly. We demonstrate the proposed approach by both simulation examples and a real data application.     

On intra-annual EVI variability in the dry season of tropical forest: A case study with MODIS and hyperspectral data

Considerable controversy is associated with dry season increases in the Enhanced Vegetation Index (EVI), observed using the Moderate Resolution Imaging Spectroradiometer (MODIS), compared with field-based estimates of decreasing plant productivity. Here, we investigate potential causes of intra-annual variability by comparing EVI from mature forest with field-measured Leaf Area Index (LAI) to validate space-based observations. EVI was calculated from 19 nadir and off-nadir Hyperion images in the 2005 dry season, and inspected for consistency with MODIS observations from 2004 to 2009. The objective was to evaluate the possible influence of the view-illumination geometry and of canopy foliage and leaf flush on the EVI. Spectral mixture models were used to evaluate the relationship between EVI and the shade fraction, a measure that varies with pixel brightness. MODIS LAI values were compared with LAI estimated using hemispherical photographs taken in two field campaigns in the dry season. To keep LAI and leaf flush conditions as constant variables and vary solar illumination, we used airborne Hyperspectral Mapper (Hymap) data acquired over mature forest from another region on the same day but with two distinct solar zenith angles (SZA) (29° and 53°). Results showed that intra-annual variability in MODIS and nadir Hyperion EVI in the dry season of tropical forest were driven by solar illumination effects rather than changes in LAI. The reflectance of the MODIS and Hyperion blue, red and near infrared (NIR) bands was higher at the end of the dry season because of the predominance of sunlit canopy components for the sensors due to decreasing SZA from June (44°) to September (26°). Because EVI was highly correlated with the reflectance of the NIR band used to generate it (r of + 0.98 for MODIS and + 0.88 for Hyperion), this vegetation index followed the general NIR pattern, increasing with smaller SZA towards the end of the dry season. Hyperion EVI was inversely correlated with the shade fraction (r = − 0.93). Changes in canopy foliage detected from MODIS LAI data were not consistent with LAI estimates from hemispherical photographs. Although further research is necessary to measure the impact of leaf flush on intra-annual EVI variability in the Querência region, analysis of Hymap data with fixed LAI and leaf flush conditions confirmed the influence of the illumination effects on the EVI.     

Retrieval of high-resolution sea surface temperature data for Sendai Bay, Japan, using the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER)

We studied sea surface temperature (SST) retrieval algorithms for Sendai Bay, using output from the thermal-infrared channels of the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) on board Terra. While the highest resolutions of other satellite SST products are about 1 km, the ASTER thermal-infrared channels provide 90-m spatial resolution. To develop the ASTER algorithm, we employed statistical methods in which SSTs retrieved from the thermal-infrared measurements were tuned against the Moderate Resolution Imaging Spectroradiometer (MODIS) SST product with a 1-km spatial resolution. Terra also carries a MODIS sensor, which observed the same area as the ASTER sensor at the same time. The MODIS SST was validated around Sendai Bay, revealing a bias of −0.15 °C and root mean-square difference (RMSD) of 0.67 °C against in situ SSTs. Taking into account the spatial-resolution difference between ASTER and MODIS, match-up was generated only if the variability of ASTER brightness temperatures (T₁₃) was small in a pixel of MODIS SST (MP). The T₁₃ within one MP was about 121 pixels. The standard deviation (σ₁₃) of T₁₃ was calculated for each cloud-free MP, and the threshold of σ₁₃ for choosing match-up MPs was decided by analyzing the σ₁₃ histogram of one ASTER image. The 15 synchronous pairs of ASTER/MODIS images are separated into two groups of 8 pairs called set (A) and 7 pairs called set (B). Using the common procedure, the match-ups are generated for set (A) and set (B). The former is used for developing the ASTER Multi-Channel SST (MCSST) algorithm, and the latter for validation of the developed ASTER SST. Analysis of the whole 15 pairs indicated that ASTER SST does not depend on the satellite zenith angle. We concluded that, using Akaike's information criterion with set (A) match-ups, the multiple regression formula with all five thermal-infrared channels was adequate for the ASTER SST retrieval. Validation of ASTER SST using match-up set (B) indicated a bias of 0.101 °C and RMSD of 0.455 °C.     

Modeling MODIS LAI time series using three statistical methods

Leaf Area Index (LAI) is one of the most important variables characterizing land surface vegetation and dynamics. Many satellite data, such as the Moderate Resolution Imaging Spectroradiometer (MODIS), have been used to generate LAI products. It is important to characterize their spatial and temporal variations by developing mathematical models from these products. In this study, we aim to model MODIS LAI time series and further predict its future values by decomposing the LAI time series of each pixel into several components: trend, intra-annual variations, seasonal cycle, and stochastic stationary or irregular parts. Three such models that can characterize the non-stationary time series data and predict the future values are explored, including Dynamic Harmonics Regression (DHR), STL (Seasonal-Trend Decomposition Procedure based on Loess), and Seasonal ARIMA (AutoRegressive Intergrated Moving Average) (SARIMA). The preliminary results using six years (2001-2006) of the MODIS LAI product indicate that all these methods are effective to model LAI time series and predict 2007 LAI values reasonably well. The SARIMA model gives the best prediction, DHR produces the smoothest curve, and STL is more sensitive to noise in the data. These methods work best for land cover types with pronounced seasonal variations.     

Mapping leaf chlorophyll and leaf area index using inverse and forward canopy reflectance modeling and SPOT reflectance data

Reflectance data in the green, red and near-infrared wavelength region were acquired by the SPOT high resolution visible and geometric imaging instruments for an agricultural area in Denmark (56°N, 9°E) for the purpose of estimating leaf chlorophyll content (C_ab) and green leaf area index (LAI). SPOT reflectance observations were atmospherically corrected using aerosol data from MODIS and profiles of air temperature, humidity and ozone from the Atmospheric Infrared Sounder (AIRS), and used as input for the inversion of a canopy reflectance model. Computationally efficient inversion schemes were developed for the retrieval of soil and land cover-specific parameters which were used to build multiple species and site dependent formulations relating the two biophysical properties of interest to vegetation indices or single spectral band reflectances. Subsequently, the family of model generated relationships, each a function of soil background and canopy characteristics, was employed for a fast pixel-wise mapping of C_ab and LAI.The biophysical parameter retrieval scheme is completely automated and image-based and solves for the soil background reflectance signal, leaf mesophyll structure, specific dry matter content, Markov clumping characteristics, C_ab and LAI without utilizing calibration measurements.Despite the high vulnerability of near-infrared reflectances (ρ_nir) to variations in background properties, an efficient correction for background influences and a strong sensitivity of ρ_nir to LAI, caused LAI-ρ_nir relationships to be very useful and preferable over LAI-NDVI relationships for LAI prediction when LAI > 2. Reflectances in the green waveband (ρ_green) were chosen for producing maps of C_ab.The application of LAI-NDVI, LAI-ρ_nir and C_ab-ρ_green relationships provided reliable quantitative estimates of C_ab and LAI for agricultural crops characterized by contrasting architectures and leaf biochemical constituents with overall root mean square deviations between estimates and in-situ measurements of 0.74 for LAI and 5.0 μg cm^− 2 for C_ab.The results of this study illustrate the non-uniqueness of spectral reflectance relationships and the potential of physically-based inverse and forward canopy reflectance modeling techniques for a reasonably fast and accurate retrieval of key biophysical parameters at regional scales.     

Estimating the leaf area index in Indian tropical forests using Landsat-8 OLI data

Leaf area index (LAI) is a key vegetation biophysical parameter and is extensively used in modelling of phenology, primary production, light interception, evapotranspiration, carbon, and nitrogen dynamics. In the present study, we attempt to spatially characterize LAI for natural forests of Western Ghats India, using ground based and Landsat-8 Operational Land Imager (OLI) sensor satellite data. For this, 41 ground-based LAI measurements were carried out across a gradient of tropical forest types, viz. dry, moist, and evergreen forests using LAI-2200 plant canopy analyser, during the month of March 2015. Initially, measured LAI values were regressed with 15 spectral variables, including nine spectral vegetation indices (SVIs) and six Landsat-8 surface reflectance (ρ) variables using univariate correlation analysis. Results showed that the red (ρ_red), near-infrared (ρ_NIR), shortwave infrared (ρ_SWIR1, ρ_SWIR2) reflectance bands (R² > 0.6), and all SVIs (R² > 0.7) except simple ratio (SR) have the highest and second highest coefficient of determination with ground-measured LAI. In the second step, to select significant (high R², low root mean square error (RMSE), and p-level < 0.05) SVIs to determine the best representative model, stepwise multiple linear regression (SMLR) was implemented. The results indicate that the SMLR model predicted LAI with better coefficient of determination (R² = 0.83, RMSE = 0.78) using normalized difference vegetation index, enhanced vegetation index, and soil-adjusted vegetation index variables compared to the univariate approach. The predicted SMLR model was used to estimate a spatial map of LAI. It is desirable to evaluate the stability and potentiality of regional LAI models in natural forest ecosystems against the operationally accepted Moderate Resolution Imaging Spectroradiometer (MODIS) global LAI product. To do this, the Landsat-8 pixel-based LAI map was resampled to 1 km resolution and compared with the MODIS derived LAI map. Results suggested that Landsat-8 OLI-based VIs provide significant LAI maps at moderate resolution (30 m) as well as coarse resolution (1 km) for regional climate models.