Evaluating the capabilities of Sentinel-2 for quantitative estimation of biophysical variables in vegetation

The red edge position (REP) in the vegetation spectral reflectance is a surrogate measure of vegetation chlorophyll content, and hence can be used to monitor the health and function of vegetation. The Multi-Spectral Instrument (MSI) aboard the future ESA Sentinel-2 (S-2) satellite will provide the opportunity for estimation of the REP at much higher spatial resolution (20 m) than has been previously possible with spaceborne sensors such as Medium Resolution Imaging Spectrometer (MERIS) aboard ENVISAT. This study aims to evaluate the potential of S-2 MSI sensor for estimation of canopy chlorophyll content, leaf area index (LAI) and leaf chlorophyll concentration (LCC) using data from multiple field campaigns. Included in the assessed field campaigns are results from SEN3Exp in Barrax, Spain composed of 35 elementary sampling units (ESUs) of LCC and LAI which have been assessed for correlation with simulated MSI data using a CASI airborne imaging spectrometer. Analysis also presents results from SicilyS2EVAL, a campaign consisting of 25 ESUs in Sicily, Italy supported by a simultaneous Specim Aisa-Eagle data acquisition. In addition, these results were compared to outputs from the PROSAIL model for similar values of biophysical variables in the ESUs. The paper in turn assessed the scope of S-2 for retrieval of biophysical variables using these combined datasets through investigating the performance of the relevant Vegetation Indices (VIs) as well as presenting the novel Inverted Red-Edge Chlorophyll Index (IRECI) and Sentinel-2 Red-Edge Position (S2REP). Results indicated significant relationships between both canopy chlorophyll content and LAI for simulated MSI data using IRECI or the Normalised Difference Vegetation Index (NDVI) while S2REP and the MERIS Terrestrial Chlorophyll Index (MTCI) were found to have the strongest correlation for retrieval of LCC.     

Retrieval of leaf area index in alpine wetlands using a two-layer canopy reflectance model

In this paper, we focused on the retrieval of the LAI in an alpine wetland located in western part of China in late August and early July 2011. A two-layer canopy reflectance model (ACRM) was used to establish the relationships between the LAI and the reflectance of near-infrared (NIR) and red (RED) wavebands. The reflectance data were derived from Landsat TM L1T product and the Terra and Aqua MODIS 16-day and 8-day composite reflectance products (MOD/MYD09) at 250 m resolution. Due to the lack of the information about some major input parameters for ACRM, which are sensitive to model outputs in the reflectance of NIR and RED wavebands, the inverse problem was ill-posed. To overcome this problem, a method of increasing the sensitivity of the LAI while reducing the influence of other model free parameters based on the study of free parameters’ sensitivity to the ACRM outputs and the region’s features was studied. The area of interest was divided into two parts using the approximately statistic normalized difference vegetation index (NDVI) value around 0.5. One part was sparse vegetation (0.1 < NDVI < 0.5), which is more sensitive to soil background effects and less sensitive to the canopy biophysical and biochemical variables. The other part was dense vegetation (0.5 ≤ NDVI < 1.0), which is less sensitive to soil background effects and more sensitive to plant canopies and leaf parameters. Then, the relationships of ρ_nir–LAI and ρ_red–LAI were established using a look-up table algorithm for the two parts. Furthermore, a regularization technique for fast pixel-wise retrieval was introduced to reduce the elements of LUT sets while maintaining a relatively high accuracy. The results were very promising compared to the field measured LAI values that the correlation (R²) of the measured LAI values and retrieved LAI values reached 0.95, and the root-mean-square deviation (RMSD) was 0.33 for late August, 2011, while the R² reached 0.82 and RMSD was 0.25 for early July 2011.     

Retrieval of wheat leaf area index from AWiFS multispectral data using canopy radiative transfer simulation

Accurate representation of leaf area index (LAI) from high resolution satellite observations is obligatory for various modelling exercises and predicting the precise farm productivity. Present study compared the two retrieval approach based on canopy radiative transfer (CRT) method and empirical method using four vegetation indices (VI) (e.g. NDVI, NDWI, RVI and GNDVI) to estimate the wheat LAI. Reflectance observations available at very high (56 m) spatial resolution from Advanced Wide-Field Sensor (AWiFS) sensor onboard Indian Remote Sensing (IRS) P6, Resourcesat-1 satellite was used in this study. This study was performed over two different wheat growing regions, situated in different agro-climatic settings/environments: Trans-Gangetic Plain Region (TGPR) and Central Plateau and Hill Region (CPHR). Forward simulation of canopy reflectances in four AWiFS bands viz. green (0.52–0.59 μm), red (0.62–0.68 μm), NIR (0.77–0.86 μm) and SWIR (1.55–1.70 μm) were carried out to generate the look up table (LUT) using CRT model PROSAIL from all combinations of canopy intrinsic variables. An inversion technique based on minimization of cost function was used to retrieve LAI from LUT and observed AWiFS surface reflectances. Two consecutive wheat growing seasons (November 2005–March 2006 and November 2006–March 2007) datasets were used in this study. The empirical models were developed from first season data and second growing season data used for validation. Among all the models, LAI-NDVI empirical model showed the least RMSE (root mean square error) of 0.54 and 0.51 in both agro-climatic regions respectively. The comparison of PROSAIL retrieved LAI with in situ measurements of 2006–2007 over the two agro-climatic regions produced substantially less RMSE of 0.34 and 0.41 having more R² of 0.91 and 0.95 for TGPR and CPHR respectively in comparison to empirical models. Moreover, CRT retrieved LAI had less value of errors in all the LAI classes contrary to empirical estimates. The PROSAIL based retrieval has potential for operational implementation to determine the regional crop LAI and can be extendible to other regions after rigorous validation exercise.     

Estimation of leaf area index using PROSAIL based LUT inversion,MLRA-GPR and empirical models: Case study of tropical deciduous forest plantation,North India

Forests play a vital role in biological cycles and environmental regulation. To understand the key processes of forest canopies (e.g., photosynthesis, respiration and transpiration), reliable and accurate information on spatial variability of Leaf Area Index (LAI), and its seasonal dynamics is essential. In the present study, we assessed the performance of biophysical parameter (LAI) retrieval methods viz. Look-Up Table (LUT)-inversion, MLRA-GPR (Machine Learning Regression Algorithm- Gaussian Processes Regression) and empirical models, for estimating the LAI of tropical deciduous plantation using ARTMO (Automated Radiative Transfer Models Operator) tool and Sentinel-2 satellite images. The study was conducted in Central Tarai Forest Division, Haldwani, located in the Uttarakhand state, India. A total of 49 ESUs (Elementary Sampling Unit) of 30 m × 30 m size were established based on variability in composition and age of plantation stands. In-situ LAI was recorded using plant canopy imager during the leaf growing, peak and senescence seasons. The PROSAIL model was calibrated with site-specific biophysical and biochemical parameters before used to the predicted LAI. The plantation LAI was also predicted by an empirical approach using optimally chosen Sentinel-2 vegetation indices. In addition, Sentinel-2 and MODIS LAI products were evaluated with respect to LAI measurements. MLRA-GPR offered best results for predicting LAI of leaf growing (R² = 0.9, RMSE = 0.14), peak (R² = 0.87, RMSE = 0.21) and senescence (R² = 0.86, RMSE = 0.31) seasons while LUT inverted model outperformed VI’s based parametric regression model. Vegetation indices (VIs) derived from 740 nm, 783 nm and 2190 nm band combinations of Sentinel-2 offered the best prediction of LAI.     

Retrieval of boreal forest LAI using a forest reflectance model and empirical regressions

Spectral invariants provide a novel approach for characterizing canopy structure in forest reflectance models and for mapping biophysical variables using satellite images. We applied a photon recollision probability (p) based forest reflectance model (PARAS) to retrieve leaf area index (LAI) from fine resolution SPOT HRVIR and Landsat ETM+ satellite data. First, PARAS was parameterized using an extensive database of LAI-2000 measurements from five conifer-dominated boreal forest sites in Finland, and mixtures of field-measured forest understory spectra. The selected vegetation indices (e.g. reduced simple ratio, RSR), neural networks and kNN method were used to retrieve effective LAI (L_e) based on reflectance model simulations. For comparison, we established empirical vegetation index-LAI regression models for our study sites. The empirical RSR–L_e regression performed best when applied to an independent test site in southern Finland [RMSE 0.57 (24.2%)]. However, the difference to the best reflectance model based retrievals produced by neural networks was only marginal [RMSE 0.59 (25.1%)]. According to this study, the PARAS model provides a simple and flexible modelling tool for calibrating algorithms for LAI retrieval in conifer-dominated boreal forests. The advantage of PARAS is that it directly uses field measurements to parameterize canopy structure (LAI-2000, hemispherical photographs) and optical properties of foliage and understory.     

Estimation of vegetation LAI from hyperspectral reflectance data: Effects of soil type and plant architecture

The retrieval of canopy biophysical variables is known to be affected by confounding factors such as plant type and background reflectance. The effects of soil type and plant architecture on the retrieval of vegetation leaf area index (LAI) from hyperspectral data were assessed in this study. In situ measurements of LAI were related to reflectances in the red and near-infrared and also to five widely used spectral vegetation indices (VIs). The study confirmed that the spectral contrast between leaves and soil background determines the strength of the LAI–reflectance relationship. It was shown that within a given vegetation species, the optimum spectral regions for LAI estimation were similar across the investigated VIs, indicating that the various VIs are basically summarizing the same spectral information for a given vegetation species. Cross-validated results revealed that, narrow-band PVI was less influenced by soil background effects (0.15 ≤ RMSE_cv ≤ 0.56). The results suggest that, when using remote sensing VIs for LAI estimation, not only is the choice of VI of importance but also prior knowledge of plant architecture and soil background. Hence, some kind of landscape stratification is required before using hyperspectral imagery for large-scale mapping of vegetation biophysical variables.     

Principal component inversion technique for the retrieval of leaf area index

A simulation study has been carried out to investigate the Principal Component Inversion (PCI) technique for the retrieval of leaf area index (LAI). The PROSAIL model has been used for the forward analysis, i.e., estimation of reflectance for various combinations of LAI, soil reflectance, leaf angle distribution (ø₁), chlorophyll a+b concentration (C_ab), etc. Independent test on sample with LAI range 0.1-7.0 indicated that the retrieved LAI from PCI has higher accuracy (RMSE=0.137) than the classical NDVI-LAI empirical approach (RMSE=1.139). The study needs to be extended to cover retrieval from different types of soil and simultaneous retrieval of different biophysical parameters viz., LAI, C_ab, ø₁ to test the wider applicability of the PCI technique.     

Comparative analysis of different retrieval methods for mapping grassland leaf area index using airborne imaging spectroscopy

Fine scale maps of vegetation biophysical variables are useful status indicators for monitoring and managing national parks and endangered habitats. Here, we assess in a comparative way four different retrieval methods for estimating leaf area index (LAI) in grassland: two radiative transfer model (RTM) inversion methods (one based on look-up-tables (LUT) and one based on predictive equations) and two statistical modelling methods (one partly, the other entirely based on in situ data). For prediction, spectral data were used that had been acquired over Majella National Park in Italy by the airborne hyperspectral HyMap instrument. To assess the performance of the four investigated models, the normalized root mean squared error (nRMSE) and coefficient of determination (R²) between estimates and in situ LAI measurements are reported (n = 41). Using a jackknife approach, we also quantified the accuracy and robustness of empirical models as a function of the size of the available calibration data set. The results of the study demonstrate that the LUT-based RTM inversion yields higher accuracies for LAI estimation (R² = 0.91, nRMSE = 0.18) as compared to RTM inversions based on predictive equations (R² = 0.79, nRMSE = 0.38). The two statistical methods yield accuracies similar to the LUT method. However, as expected, the accuracy and robustness of the statistical models decrease when the size of the calibration database is reduced to fewer samples. The results of this study are of interest for the remote sensing community developing improved inversion schemes for spaceborne hyperspectral sensors applicable to different vegetation types. The examples provided in this paper may also serve as illustrations for the drawbacks and advantages of physical and empirical models.     

Soil moisture retrieval from airborne L-band passive microwave using high resolution multispectral data

For the soil moisture retrieval from passive microwave sensors, such as ESA’s Soil Moisture and Ocean Salinity (SMOS) and the NASA Soil Moisture Active and Passive (SMAP) mission, a good knowledge about the vegetation characteristics is indispensable. Vegetation cover is a principal factor in the attenuation, scattering and absorption of the microwave emissions from the soil; and has a direct impact on the brightness temperature by way of its canopy emissions. Here, brightness temperatures were measured at three altitudes across the TERENO (Terrestrial Environmental Observatories) Rur catchment site in Germany to achieve a range of spatial resolutions using the airborne Polarimetric L-band Multibeam Radiometer 2 (PLMR2). The L-band Microwave Emission of the Biosphere (L-MEB) model which simulates microwave emissions from the soil–vegetation layer at L-band was used to retrieve surface soil moisture for all resolutions. A Monte Carlo approach was developed to simultaneously estimate soil moisture and the vegetation parameter b’ describing the relationship between the optical thickness τ and the Leaf Area Index (LAI). LAI was retrieved from multispectral RapidEye imagery and the plant specific vegetation parameter b′ was estimated from the lowest flight altitude data for crop, grass, coniferous forest, and deciduous forest. Mean values of b’ were found to be 0.18, 0.07, 0.26 and 0.23, respectively. By assigning the estimated b′ to higher flight altitude data sets, a high accuracy soil moisture retrieval was achieved with a Root Mean Square Difference (RMSD) of 0.035 m³ m⁻³ when compared to ground-based measurements.     

Evaluation of SPOT imagery for the estimation of grassland biomass

In many regions, a decrease in grasslands and change in their management, which are associated with agricultural intensification, have been observed in the last half-century. Such changes in agricultural practices have caused negative environmental effects that include water pollution, soil degradation and biodiversity loss. Moreover, climate-driven changes in grassland productivity could have serious consequences for the profitability of agriculture. The aim of this study was to assess the ability of remotely sensed data with high spatial resolution to estimate grassland biomass in agricultural areas. A vegetation index, namely the Normalized Difference Vegetation Index (NDVI), and two biophysical variables, the Leaf Area Index (LAI) and the fraction of Vegetation Cover (fCOVER) were computed using five SPOT images acquired during the growing season. In parallel, ground-based information on grassland growth was collected to calculate biomass values. The analysis of the relationship between the variables derived from the remotely sensed data and the biomass observed in the field shows that LAI outperforms NDVI and fCOVER to estimate biomass (R² values of 0.68 against 0.30 and 0.50, respectively). The squared Pearson correlation coefficient between observed and estimated biomass using LAI derived from SPOT images reached 0.73. Biomass maps generated from remotely sensed data were then used to estimate grass reserves at the farm scale in the perspective of operational monitoring and forecasting.     

The Effect of Radiometric Resolution on the Retrieval of Leaf Area Index from Agricultural Crops

Motivated by the increasing availability of remote sensing data of high radiometric resolution, a study was conducted to determine whether high-resolution data (10 bits or more) yielded more accurate vegetation leaf area index (LAI) information than low-radiometric-resolution multispectral data (7-bit or less). The study evaluated the performance of simulated 12-bit LISS-III sensor data (derived from EO-1 hyperspectral Hyperion data) with original 7-bit LISS-III sensor data for the estimation of LAI of major agricultural crops (e.g., cotton, sugar cane, and rice). There was no significant improvement in the correlation coefficient encountered when using the high-radiometric-resolution (12-bit) LISS-III data versus the low-radiometric-resolution (7-bit) LISS-III data for the retrieval of LAI. The retrieval of LAI of agricultural crops met with moderate success, with overall correlation coefficients of around 0.55. These results suggest that satellite data of very high radiometric resolution may not be a required for remote measurement of LAI. Spectral bandwidth, band placement, and the method of retrieving biophysical parameters may be more important.     

冠层反射光谱对植被理化参数的全局敏感性分析

植被理化参数与许多有关植物物质能量交换的生态过程密切相关,定量分析植被反射光谱对理化参数的敏感性是遥感反演理化参数含量的前提。本文采用EFAST(Extended Fourier Amplitude Sensitivity Test)全局敏感性分析方法,利用PROSAIL辐射传输模型分析了冠层疏密程度对叶片生化组分含量、冠层结构以及土壤背景等多种参数敏感性的影响,并对植被理化参数反演所需先验知识的精度问题进行了初步探讨。研究表明:(1)对于较为稠密的冠层,可见光波段的冠层反射率主要受叶绿素含量的影响,近红外和中红外波段的冠层反射率主要受干物质量和含水量的影响;(2)对于稀疏的冠层,LAI是影响400—2500 nm波段范围内冠层反射率的最重要参数,土壤湿度次之,叶片生化参数对冠层反射率的敏感性较低;(3)在已知稀疏冠层LAI的情况下进一步确定土壤的干湿状态,可显著提高冠层反射率对叶绿素含量的敏感度,有助于稀疏冠层叶绿素含量的反演。     

A Spatio-Temporal Enhancement Method for medium resolution LAI (STEM-LAI)

Satellite remote sensing has been used successfully to map leaf area index (LAI) across landscapes, but advances are still needed to exploit multi-scale data streams for producing LAI at both high spatial and temporal resolution. A multi-scale Spatio-Temporal Enhancement Method for medium resolution LAI (STEM-LAI) has been developed to generate 4-day time-series of Landsat-scale LAI from existing medium resolution LAI products. STEM-LAI has been designed to meet the demands of applications requiring frequent and spatially explicit information, such as effectively resolving rapidly evolving vegetation dynamics at sub-field (30 m) scales. In this study, STEM-LAI is applied to Moderate Resolution Imaging Spectroradiometer (MODIS) based LAI data and utilizes a reference-based regression tree approach for producing MODIS-consistent, but Landsat-based, LAI. The Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) is used to interpolate the downscaled LAI between Landsat acquisition dates, providing a high spatial and temporal resolution improvement over existing LAI products. STARFM predicts high resolution LAI by blending MODIS and Landsat based information from a common acquisition date, with MODIS data from a prediction date. To demonstrate its capacity to reproduce fine-scale spatial features observed in actual Landsat LAI, the STEM-LAI approach is tested over an agricultural region in Nebraska. The implementation of a 250 m resolution LAI product, derived from MODIS 1 km data and using a scale consistent approach based on the Normalized Difference Vegetation Index (NDVI), is found to significantly improve accuracies of spatial pattern prediction, with the coefficient of efficiency (E) ranging from 0.77–0.94 compared to 0.01–0.85 when using 1 km LAI inputs alone. Comparisons against an 11-year record of in-situ measured LAI over maize and soybean highlight the utility of STEM-LAI in reproducing observed LAI dynamics (both characterized by r² = 0.86) over a range of plant development stages. Overall, STEM-LAI represents an effective downscaling and temporal enhancement mechanism that predicts in-situ measured LAI better than estimates derived through linear interpolation between Landsat acquisitions. This is particularly true when the in-situ measurement date is greater than 10 days from the nearest Landsat acquisition, with prediction errors reduced by up to 50%. With a streamlined and completely automated processing interface, STEM-LAI represents a flexible tool for LAI disaggregation in space and time that is adaptable to different land cover types, landscape heterogeneities, and cloud cover conditions.     

Applying different inversion techniques to retrieve stand variables of summer barley with PROSPECT + SAIL

This study describes the retrieval of state variables (LAI, canopy chlorophyll, water and dry matter contents) for summer barley from airborne HyMap data by means of a canopy reflectance model (PROSPECT + SAIL). Three different inversion techniques were applied to explore the impact of the employed method on estimation accuracies: numerical optimization (downhill simplex method), a look-up table (LUT) and an artificial neural network (ANN) approach. By numerical optimization (Num Opt), reliable estimates were obtained for LAI and canopy chlorophyll contents (LAI × C_ab) with r² of 0.85 and 0.94 and RDP values of 1.81 and 2.65, respectively. Accuracies dropped for canopy water (LAI × C_w) and dry matter contents (LAI × C_m). Nevertheless, the range of leaf water contents (C_w) was very narrow in the studied plant material. Prediction accuracies generally decreased in the order Num Opt > LUT > ANN. This decrease in accuracy mainly resulted from an increase in offset in the obtained values, as the retrievals from the different approaches were highly correlated. The same decreasing order in accuracy was found for the difference between the measured spectra and those reconstructed from the retrieved variable values. The parallel application of the different inversion techniques to one collective data set was helpful to identify modelling uncertainties, as shortcomings of the retrieval algorithms themselves could be separated from uncertainties in model structure and parameterisation schemes.     

Retrieval of remotely sensed LAI using Landsat ETM+ data and ground measurements of solar radiation and vegetation structure: Implication of leaf inclination angle

A time series of leaf area index (LAI) of a managed birch forest in Germany (near Dresden) has been developed based on 16-day normalized difference vegetation index (NDVI) data from the Landsat ETM+ sensor at 30 m resolution. The Landsat ETM+ LAI was retrieved using a modified physical radiative transfer (RTM) model which establishes a relationship between LAI, fractional vegetation cover (fC), and given patterns of surface reflectance, view-illumination conditions and optical properties of vegetation. In situ measurements of photosynthetically active radiation (PAR) and vegetation structure parameters using hemispherical photography (HSP) served for calibration of model parameters, while data from litter collection at the study site provided the ground-based estimates of LAI for validation of modelling results. Influence of view-illumination conditions on optical properties of canopy was simulated by a view angle geometry model incorporating the solar zenith angle and the sensor viewing angle. Effects of intra-annual and inter-annual variability of structural properties of the canopy on the light extinction coefficient were simulated by implementing variability of the leaf inclination angle (LIA), which was confirmed in the study site. The results revealed good compatibility of the produced Landsat ETM+ LAI data set with the litter-estimated LAI. The results also showed high sensitivity of the LAI retrieval algorithm to variability of structural properties of the canopy: the implementation of LIA dynamics into the LAI retrieval algorithm significantly improved the model accuracy.     

An Algorithm to Produce Temporally and Spatially Continuous MODIS-LAI Time Series

Ecological and climate models require high-quality consistent biophysical parameters as inputs and validation sources. NASA's moderate resolution imaging spectroradiometer (MODIS) biophysical products provide such data and have been used to improve our understanding of climate and ecosystem changes. However, the MODIS time series contains occasional lower quality data, gaps from persistent clouds, cloud contamination, and other gaps. Many modeling efforts, such as those used in the North American Carbon Program, that use MODIS data as inputs require gap-free data. This letter presents the algorithm used within the MODIS production facility to produce temporally smoothed and spatially continuous biophysical data for such modeling applications. We demonstrate the algorithm with an example from the MODIS-leaf-area-index (LAI) product. Results show that the smoothed LAI agrees with high-quality MODIS LAI very well. Higher R-squares and better linear relationships have been observed when high-quality retrieval in each individual tile reaches 40% or more. These smoothed products show similar data quality to MODIS high-quality data and, therefore, can be substituted for low-quality retrievals or data gaps.     

Modelling the spectral response of the desert tree Prosopis tamarugo to water stress

In this paper, we carried out a laboratory experiment to study changes in canopy reflectance of Tamarugo plants under controlled water stress. Tamarugo (Prosopis tamarugo Phil.) is an endemic and endangered tree species adapted to the hyper-arid conditions of the Atacama Desert, Northern Chile. Observed variation in reflectance during the day (due to leaf movements) as well as changes over the experimental period (due to water stress) were successfully modelled by using the Soil-Leaf-Canopy (SLC) radiative transfer model. Empirical canopy reflectance changes were mostly explained by the parameters leaf area index (LAI), leaf inclination distribution function (LIDF) and equivalent water thickness (EWT) as shown by the SLC simulations. Diurnal leaf movements observed in Tamarugo plants (as adaptation to decrease direct solar irradiation at the hottest time of the day) had an important effect on canopy reflectance and were explained by the LIDF parameter. The results suggest that remote sensing based assessment of this desert tree should consider LAI and canopy water content (CWC) as water stress indicators. Consequently, we tested fifteen different vegetation indices and spectral absorption features proposed in literature for detecting changes of LAI and CWC, considering the effect of LIDF variations. A sensitivity analysis was carried out using SLC simulations with a broad range of LAI, LIDF and EWT values. The Water Index was the most sensitive remote sensing feature for estimating CWC for values less than 0.036 g/cm², while the area under the curve for the spectral range 910–1070 nm was most sensitive for values higher than 0.036 g/cm². The red-edge chlorophyll index (CI_red-edge) performed the best for estimating LAI. Diurnal leaf movements had an effect on all remote sensing features tested, particularly on those for detecting changes in CWC.     

Inversion of PROSAIL Model for Retrieval of Plant Biophysical Parameters

The current development of satellite technology particularly in the sensors like POLDER and MISR, has emphasized more on directional reflectance measurements (i.e. spectral reflectance of the target measured from different view zenith and azimuth angles) of the earth surface features mainly the vegetation for retrieval of biophysical parameters at regional scale using radiative transfer models. This approach being physical process based and uses directional reflectance measurement has been found to better and more reliable compared to the conventional statistical approach used till date and takes care of anisotropic nature (i.e. reflectance from the target is different if measured from different view angles) of the target. Keeping this in view a field experiment was conducted in mustard crop to evaluate the radiative transfer model for biophysical parameter retrieval through its inversion with the objectives set as (i) to relate canopy biophysical parameters and geometry to its bidirectional reflectance, (ii) to evaluate a canopy reflectance model to best represent the radiative transfer within the canopy for its inversion and (iii) to retrieve crop biophysical parameters through inversion of the model. Two varieties of the mustard crop (Brassica juncea L) were grown with two nitrogen treatments. The bidirectional reflectance data obtained at 5 nm interval for a range of 400–1100 nm were integrated to IRS LISS–II sensor’s four band values using Newton Cotes Integration technique. Biophysical parameters like leaf area index, leaf chlorophyll content, leaf length, plant height and average leaf inclination angle, biomass etc were estimated synchronizing with the bi-directional reflectance measurements. Radiative transfer model PROSAIL model was validated and its inversion was done to retrieve LAI and ALA. Look Up Table (LUT) of Bidirectional reflectance distribution function (BRDF) was prepared simulating through PROSAIL model varying only LAI (0.2 interval from 1.2 to 5.4 ) and ALA (5° interval from 40° to 55°) parameters and inversion was done using a merit function and numerical optimization technique given by Press et al. (1986). The derived LAI and ALA values from inversion were well matched with observed one with RMSE 0.521 and 5.57, respectively.     

Relationship of bidirectional reflectance of wheat with biophysical parameters and its radiative transfer modeling using PROSAIL

The algorithms for deriving vegetation biophysical parameters rely on the understanding of bi-directional interaction of radiation and its subsequent linkages with canopy radiative transfer models and their inversion. In this study, an attempt has been made to define the geometry of sensor and source position to best relate plant biophysical parameters with bidirectional reflectance of wheat varieties varying in canopy architecture and to validate the performance of PROSAIL (PROSPECT+SAIL) canopy radiative transfer model. A field experiment was conducted with two wheat cultivars varying in canopy geometry and phenology. The bidirectional measurements between 400nm–1100nm at 5nm interval were recorded every week at six view azimuth and four view zenith positions using spectro-radiometer. Canopy biophysical parameters were recorded synchronous to bi-directional reflectance measurements. The broadband reflectances were used to compute the NDVIs which were subsequently related to leaf area index and biomass. Results showed that the bidirectional reflectance increased with increase in view zenith from 20⁰ to 60⁰ irrespective of the sensor azimuth. For a given view zenith, the reflectance was observed to be maximum at 150⁰ azimuth where the difference between the sun and sensor azimuth was least. For sun azimuth of 160⁰ and zenith of 52⁰, the view geometry defined by 150⁰ azimuth and 50⁰ zenith corresponded to hotspot position. The measured bidirectional NDVI had significant logarithmic relationship with LAI and linear relationship with biomass for both the varieties of wheat and maximum correlation of NDVI with LAI and with biomass was obtained at the hotspot position. The PROSAIL validation results showed that the model simulated well the overall shape of spectra for all combination of view zenith and azimuth position for both wheat varieties with overall RMSE less than 5 per cent. The hotspot and dark spot positions were also well simulated and hence model performance may be suitable for deriving wheat biophysical parameters using satellite derived reflectances.