Retrieval of leaf water content spanning the visible to thermal infrared spectra

The objective of this study was to investigate the entire spectra (from visible to the thermal infrared; 0.390–14.0 μm) to retrieve leaf water content in a consistent manner. Narrow-band spectral indices (calculated from all possible two band combinations) and a partial least square regression (PLSR) were used to assess the strength of each spectral region. The coefficient of determination (R²) and root mean square error (RMSE) were used to report the prediction accuracy of spectral indices and PLSR models. In the visible-near infrared and shortwave infrared (VNIR–SWIR), the most accurate spectral index yielded R² of 0.89 and RMSE of 7.60%, whereas in the mid infrared (MIR) the highest R² was 0.93 and RMSE of 5.97%. Leaf water content was poorly predicted using two-band indices developed from the thermal infrared (R² = 0.33). The most accurate PLSR model resulted from MIR reflectance spectra (R² = 0.96, RMSE = 4.74% and RMSE cross validation RMSE_CV = 6.17%) followed by VNIR–SWIR reflectance spectra (R² = 0.91, RMSE = 6.90% and RMSE_CV = 7.32%). Using thermal infrared (TIR) spectra, the PLSR model yielded a moderate retrieval accuracy (R² = 0.67, RMSE = 13.27% and RMSE_CV = 16.39%). This study demonstrated that the mid infrared (MIR) and shortwave infrared (SWIR) domains were the most sensitive spectral region for the retrieval of leaf water content.     

Estimating winter wheat biomass by assimilating leaf area index derived from fusion of Landsat-8 and MODIS data

A sufficient number of satellite acquisitions in a growing season are essential for deriving agronomic indicators, such as green leaf area index (GLAI), to be assimilated into crop models for crop productivity estimation. However, for most high resolution orbital optical satellites, it is often difficult to obtain images frequently due to their long revisit cycles and unfavorable weather conditions. Data fusion algorithms, such as the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) and the Enhanced STARFM (ESTARFM), have been developed to generate synthetic data with high spatial and temporal resolution to address this issue. In this study, we evaluated the approach of assimilating GLAI into the Simple Algorithm for Yield Estimation model (SAFY) for winter wheat biomass estimation. GLAI was estimated using the two-band Enhanced Vegetation Index (EVI2) derived from data acquired by the Operational Land Imager (OLI) onboard the Landsat-8 and a fusion dataset generated by blending the Moderate-Resolution Imaging Spectroradiometer (MODIS) data and the OLI data using the STARFM and ESTARFM models. The fusion dataset had the temporal resolution of the MODIS data and the spatial resolution of the OLI data. Key parameters of the SAFY model were optimised through assimilation of the estimated GLAI into the crop model using the Shuffled Complex Evolution-University of Arizona (SCE-UA) algorithm. A good agreement was achieved between the estimated and field measured biomass by assimilating the GLAI derived from the OLI data (GLAI_L) alone (R² = 0.77 and RMSE = 231 g m⁻²). Assimilation of GLAI derived from the fusion dataset (GLAI_F) resulted in a R² of 0.71 and RMSE of 193 g m⁻² while assimilating the combination of GLAI_L and GLAI_F led to further improvements (R² = 0.76 and RMSE = 176 g m⁻²). Our results demonstrated the potential of using the fusion algorithms to improve crop growth monitoring and crop productivity estimation when the number of high resolution remote sensing data acquisitions is limited.     

Improving the prediction of arsenic contents in agricultural soils by combining the reflectance spectroscopy of soils and rice plants

Visible and near-infrared reflectance spectroscopy provides a beneficial tool for investigating soil heavy metal contamination. This study aimed to investigate mechanisms of soil arsenic prediction using laboratory based soil and leaf spectra, compare the prediction of arsenic content using soil spectra with that using rice plant spectra, and determine whether the combination of both could improve the prediction of soil arsenic content. A total of 100 samples were collected and the reflectance spectra of soils and rice plants were measured using a FieldSpec3 portable spectroradiometer (350–2500 nm). After eliminating spectral outliers, the reflectance spectra were divided into calibration (n = 62) and validation (n = 32) data sets using the Kennard-Stone algorithm. Genetic algorithm (GA) was used to select useful spectral variables for soil arsenic prediction. Thereafter, the GA-selected spectral variables of the soil and leaf spectra were individually and jointly employed to calibrate the partial least squares regression (PLSR) models using the calibration data set. The regression models were validated and compared using independent validation data set. Furthermore, the correlation coefficients of soil arsenic against soil organic matter, leaf arsenic and leaf chlorophyll were calculated, and the important wavelengths for PLSR modeling were extracted. Results showed that arsenic prediction using the leaf spectra (coefficient of determination in validation, R_v² = 0.54; root mean square error in validation, RMSE_v = 12.99 mg kg⁻¹; and residual prediction deviation in validation, RPD_v = 1.35) was slightly better than using the soil spectra (R_v² = 0.42, RMSE_v = 13.35 mg kg⁻¹, and RPD_v = 1.31). However, results also showed that the combinational use of soil and leaf spectra resulted in higher arsenic prediction (R_v² = 0.63, RMSE_v = 11.94 mg kg⁻¹, RPD_v = 1.47) compared with either soil or leaf spectra alone. Soil spectral bands near 480, 600, 670, 810, 1980, 2050 and 2290 nm, leaf spectral bands near 700, 890 and 900 nm in PLSR models were important wavelengths for soil arsenic prediction. Moreover, soil arsenic showed significantly positive correlations with soil organic matter (r = 0.62, p < 0.01) and leaf arsenic (r = 0.77, p < 0.01), and a significantly negative correlation with leaf chlorophyll (r = −0.67, p < 0.01). The results showed that the prediction of arsenic contents using soil and leaf spectra may be based on their relationships with soil organic matter and leaf chlorophyll contents, respectively. Although RPD of 1.47 was below the recommended RPD of >2 for soil analysis, arsenic prediction in agricultural soils can be improved by combining the leaf and soil spectra.     

Plant phenolics and absorption features in vegetation reflectance spectra near 1.66 μm

Past laboratory and field studies have quantified phenolic substances in vegetative matter from reflectance measurements for understanding plant response to herbivores and insect predation. Past remote sensing studies on phenolics have evaluated crop quality and vegetation patterns caused by bedrock geology and associated variations in soil geochemistry. We examined spectra of pure phenolic compounds, common plant biochemical constituents, dry leaves, fresh leaves, and plant canopies for direct evidence of absorption features attributable to plant phenolics. Using spectral feature analysis with continuum removal, we observed that a narrow feature at 1.66 μm is persistent in spectra of manzanita, sumac, red maple, sugar maple, tea, and other species. This feature was consistent with absorption caused by aromatic CH bonds in the chemical structure of phenolic compounds and non-hydroxylated aromatics. Because of overlapping absorption by water, the feature was weaker in fresh leaf and canopy spectra compared to dry leaf measurements. Simple linear regressions of feature depth and feature area with polyphenol concentration in tea resulted in high correlations and low errors (% phenol by dry weight) at the dry leaf (r² = 0.95, RMSE = 1.0%, n = 56), fresh leaf (r² = 0.79, RMSE = 2.1%, n = 56), and canopy (r² = 0.78, RMSE = 1.0%, n = 13) levels of measurement. Spectra of leaves, needles, and canopies of big sagebrush and evergreens exhibited a weak absorption feature centered near 1.63 μm, short ward of the phenolic compounds, possibly consistent with terpenes. This study demonstrates that subtle variation in vegetation spectra in the shortwave infrared can directly indicate biochemical constituents and be used to quantify them. Phenolics are of lesser abundance compared to the major plant constituents but, nonetheless, have important plant functions and ecological significance. Additional research is needed to advance our understanding of the spectral influences of plant phenolics and terpenes relative to dominant leaf biochemistry (water, chlorophyll, protein/nitrogen, cellulose, and lignin).     

Assessment of leaf carotenoids content with a new carotenoid index: Development and validation on experimental and model data

Leaf carotenoids content (LCar) is an important indicator of plant physiological status. Accurate estimation of LCar provides valuable insight into early detection of stress in vegetation. With spectroscopy techniques, a semi-empirical approach based on spectral indices was extensively used for carotenoids content estimation. However, established spectral indices for carotenoids that generally rely on limited measured data, might lack predictive accuracy for carotenoids estimation in various species and at different growth stages. In this study, we propose a new carotenoid index (CARI) for LCar assessment based on a large synthetic dataset simulated from the leaf radiative transfer model PROSPECT-5, and evaluate its capability with both simulated data from PROSPECT-5 and 4SAIL and extensive experimental datasets: the ANGERS dataset and experimental data acquired in field experiments in China in 2004. Results show that CARI was the index most linearly correlated with carotenoids content at the leaf level using a synthetic dataset (R² = 0.943, RMSE = 1.196 μg/cm²), compared with published spectral indices. Cross-validation results with CARI using ANGERS data achieved quite an accurate estimation (R² = 0.545, RMSE = 3.413 μg/cm²), though the RBRI performed as the best index (R² = 0.727, RMSE = 2.640 μg/cm²). CARI also showed good accuracy (R² = 0.639, RMSE = 1.520 μg/cm²) for LCar assessment with leaf level field survey data, though PRI performed better (R² = 0.710, RMSE = 1.369 μg/cm²). Whereas RBRI, PRI and other assessed spectral indices showed a good performance for a given dataset, overall their estimation accuracy was not consistent across all datasets used in this study. Conversely CARI was more robust showing good results in all datasets. Further assessment of LCar with simulated and measured canopy reflectance data indicated that CARI might not be very sensitive to LCar changes at low leaf area index (LAI) value, and in these conditions soil moisture influenced the LCar retrieval accuracy.     

Comparison of MODIS-based models for retrieving suspended particulate matter concentrations in Poyang Lake,China

Suspended particulate matter (SPM) is a key parameter describing water quality, and developing the retrieval model of SPM concentration (C_SPM) is fundamental for obtaining the spatiotemporal information of C_SPM and further for understanding, managing and protecting aquatic ecosystems. This study aimed to compare moderate resolution imaging spectroradiometer (MODIS)-based C_SPM retrieval models in order to find the optimal model for improving the C_SPM estimation in Poyang Lake. The C_SPM measurements on 27 September 2007 and their coincident MODIS Terra image were used to calibrate retrieval models with the least-squares technique. The C_SPM measurements on 31 August 2012 and the MODIS Terra image on 30 August 2012 were applied to validate the calibrated models, and the correlation coefficient (r) between the measured and estimated C_SPM values, the root mean square error (RMSE) and relative root mean square error (RRMSE) of estimation as well as the model bias evaluation result were compared to determine the optimal model for estimating the C_SPM values of Poyang Lake from MODIS images. Model calibration showed that, after two samples were removed, the exponential models of blue, green and red band, the linear model of infrared band, the cubic model of red band as well as the exponential model of red minus infrared band explained about 92%, 88%, 90%, 89%, 90% and 76% of the variation of C_SPM, respectively; while model validation indicated that, after removing two samples, the exponential models of blue and green band got biased C_SPM estimations, the agreement between the measured and estimated C_SPM values was not very high (r = <0.8) for the models with single red and infrared band, and the exponential model of red minus infrared band got the best result among all calibrated models (r = 0.87, RMSE = 22.1 mg/l, RRMSE = 52.8%). We concluded that the exponential model of red minus infrared band obtained stable C_SPM estimation and was the optimal model for C_SPM estimation in this study, and more independent datasets should be obtained to further validate our finding for improving the C_SPM estimation in Poyang Lake.     

A radiative transfer model-based method for the estimation of grassland aboveground biomass

This paper presents a novel method to derive grassland aboveground biomass (AGB) based on the PROSAILH (PROSPECT + SAILH) radiative transfer model (RTM). Two variables, leaf area index (LAI, m²m⁻², defined as a one-side leaf area per unit of horizontal ground area) and dry matter content (DMC, gcm⁻², defined as the dry matter per leaf area), were retrieved using PROSAILH and reflectance data from Landsat 8 OLI product. The result of LAI × DMC was regarded as the estimated grassland AGB according to their definitions. The well-known ill-posed inversion problem when inverting PROSAILH was alleviated using ecological criteria to constrain the simulation scenario and therefore the number of simulated spectra. A case study of the presented method was applied to a plateau grassland in China to estimate its AGB. The results were compared to those obtained using an exponential regression, a partial least squares regression (PLSR) and an artificial neural networks (ANN). The RTM-based method offered higher accuracy (R² = 0.64 and RMSE = 42.67 gm⁻²) than the exponential regression (R² = 0.48 and RMSE = 41.65 gm⁻²) and the ANN (R² = 0.43 and RMSE = 46.26 gm⁻²). However, the proposed method offered similar performance than PLSR as presented better determination coefficient than PLSR (R² = 0.55) but higher RMSE (RMSE = 37.79 gm⁻²). Although it is still necessary to test these methodologies in other areas, the RTM-based method offers greater robustness and reproducibility to estimate grassland AGB at large scale without the need to collect field measurements and therefore is considered the most promising methodology.     

Retrieving aboveground biomass of wetland Phragmites australis (common reed) using a combination of airborne discrete-return LiDAR and hyperspectral data

Wetland biomass is essential for monitoring the stability and productivity of wetland ecosystems. Conventional field methods to measure or estimate wetland biomass are accurate and reliable, but expensive, time consuming and labor intensive. This research explored the potential for estimating wetland reed biomass using a combination of airborne discrete-return Light Detection and Ranging (LiDAR) and hyperspectral data. To derive the optimal predictor variables of reed biomass, a range of LiDAR and hyperspectral metrics at different spatial scales were regressed against the field-observed biomasses. The results showed that the LiDAR-derived H_p99 (99th percentile of the LiDAR height) and hyperspectral-calculated modified soil-adjusted vegetation index (MSAVI) were the best metrics for estimating reed biomass using the single regression model. Although the LiDAR data yielded a higher estimation accuracy compared to the hyperspectral data, the combination of LiDAR and hyperspectral data produced a more accurate prediction model for reed biomass (R² = 0.648, RMSE = 167.546 g/m², RMSE_r = 20.71%) than LiDAR data alone. Thus, combining LiDAR data with hyperspectral data has a great potential for improving the accuracy of aboveground biomass estimation.     

Uncertainty analysis of remote sensing of colored dissolved organic matter: Evaluations and comparisons for three rivers in North America

The uncertainties involved in remote sensing inversion of CDOM (Colored Dissolved Organic Matter) were analyzed in estuarine and coastal regions of three North American rivers: Mississippi, Hudson, and Neponset. Water optical and biogeochemical properties, including CDOM absorption and above-surface spectra, were collected in very high resolution. CDOM’s concentrations (a_g(440), absorption coefficient at 440 nm) were inverted from EO-1 Hyperion images, using a quasi-analytical algorithm for CDOM (QAA-CDOM). Uncertainties are classified to five levels, in which the underwater measurement uncertainty (level 1), image preprocessing uncertainty (level 4) and inverse model uncertainty (level 5) were evaluated. Results indicate that at level 1, in situ CDOM measurement is significant with 0.1 in the unit of QSU and 0.01 in the unit of a_g(440) (m⁻¹). At level 4, surface wave is a potential uncertainty source for high-resolution images in estuarine and coastal regions. The remote sensing reflectance of wavy water is about 10 times of the truth. At level 5, the overall uncertainty of QAA-CDOM inversion is 0.006 m⁻¹, with accuracy R² = 0.77, k = 1.1 and RMSE_log = 0.33 m⁻¹. The correlations between uncertainties and other water properties indicate that the large uncertainty in some rivers, such as the Neponset and Atchafalaya, might be caused by high-concentration chlorophyll or sediments. The relationships among the three level uncertainties show that the level 1 uncertainty generally does not propagate into level 4 and 5, but the large uncertainty at level 4 usually introduce large uncertainty at level 5.     

Determination of the crop row orientations from Formosat-2 multi-temporal and panchromatic images

This paper presents a technique developed for the retrieval of the orientation of crop rows, over anthropic lands dedicated to agriculture in order to further improve estimate of crop production and soil erosion management. Five crop types are considered: wheat, barley, rapeseed, sunflower, corn and hemp. The study is part of the multi-sensor crop-monitoring experiment, conducted in 2010 throughout the agricultural season (MCM’10) over an area located in southwestern France, near Toulouse. The proposed methodology is based on the use of satellite images acquired by Formosat-2, at high spatial resolution in panchromatic and multispectral modes (with spatial resolution of 2 and 8 m, respectively). Orientations are derived and evaluated for each image and for each plot, using directional spatial filters (45° and 135°) and mathematical morphology algorithms. “Single-date” and “multi-temporal” approaches are considered. The single-date analyses confirm the good performances of the proposed method, but emphasize the limitation of the approach for estimating the crop row orientation over the whole landscape with only one date. The multi-date analyses allow (1) determining the most suitable agricultural period for the detection of the row orientations, and (2) extending the estimation to the entire footprint of the study area. For the winter crops (wheat, barley and rapeseed), best results are obtained with images acquired just after harvest, when surfaces are covered by stubbles or during the period of deep tillage (0.27 > R² > 0.99 and 7.15° > RMSE > 43.02°). For the summer crops (sunflower, corn and hemp), results are strongly crop and date dependents (0 > R² > 0.96, 10.22° > RMSE > 80°), with a well-marked impact of flowering, irrigation equipment and/or maximum crop development. Last, the extent of the method to the whole studied zone allows mapping 90% of the crop row orientations (more than 45,000 ha) with an error inferior to 40°, associated to a confidence index ranging from 1 to 5 for each agricultural plot.     

Application of 3D triangulations of airborne laser scanning data to estimate boreal forest leaf area index

We propose 3D triangulations of airborne Laser Scanning (ALS) point clouds as a new approach to derive 3D canopy structures and to estimate forest canopy effective LAI (LAI_e). Computational geometry and topological connectivity were employed to filter the triangulations to yield a quasi-optimal relationship with the field measured LAI_e. The optimal filtering parameters were predicted based on ALS height metrics, emulating the production of maps of LAI_e and canopy volume for large areas. The LAI_e from triangulations was validated with field measured LAI_e and compared with a reference LAI_e calculated from ALS data using logarithmic model based on Beer’s law. Canopy transmittance was estimated using All Echo Cover Index (ACI), and the mean projection of unit foliage area (β) was obtained using no-intercept regression with field measured LAI_e. We investigated the influence species and season on the triangulated LAI_e and demonstrated the relationship between triangulated LAI_e and canopy volume. Our data is from 115 forest plots located at the southern boreal forest area in Finland and for each plot three different ALS datasets were available to apply the triangulations. The triangulation approach was found applicable for both leaf-on and leaf-off datasets after initial calibration. Results showed the Root Mean Square Errors (RMSEs) between LAI_e from triangulations and field measured values agreed the most using the highest pulse density data (RMSE = 0.63, the coefficient of determination (R²) = 0.53). Yet, the LAI_e calculated using ACI-index agreed better with the field measured LAI_e (RMSE = 0.53 and R² = 0.70). The best models to predict the optimal alpha value contained the ACI-index, which indicates that within-crown transmittance is accounted by the triangulation approach. The cover indices may be recommended for retrieving LAI_e only, but for applications which require more sophisticated information on canopy shape and volume, such as radiative transfer models, the triangulation approach may be preferred.     

Determining iron content in Mediterranean soils in partly vegetated areas,using spectral reflectance and imaging spectroscopy

The possibility of quantifying iron content in the topsoil of the slopes of the El Hacho Mountain complex in Southern Spain using imaging spectroscopy is investigated. Laboratory, field and airborne spectrometer (ROSIS) data are acquired, in combination with soil samples, which are analysed for dithionite extractable iron (Fe_d) content. Analysis of the properties of two iron related absorption features present in laboratory spectra demonstrates good relations, especially between the standard deviation (S.D.) of the values in an absorption feature and the Fe_d content (R² = 0.67) as well as the ratio based Redness Index (R² = 0.51). Such derived relations are less strong for the ROSIS data (R² for S.D. = 0.26 and R² for Redness Index = 0.22). The spatial distribution of iron in vegetated areas shows a strong sensitivity of these relations with the presence of vegetation. A combination of both methods shows that the overestimation of the Fe_d content with the one method is (partly) compensated by the underestimation with the other method.     

Upscaling plot-scale soil respiration in winter wheat and summer maize rotation croplands in Julu County,North China

Soil respiration (R_s) data from 45 plots were used to estimate the spatial patterns of R_s during the peak growing seasons of winter wheat and summer maize in Julu County, North China, by combining satellite remote sensing data, field-measured data, and a support vector regression (SVR) model. The observed R_s values were well reproduced by the model at the plot scale, with a root-mean-square error (RMSE) of 0.31 μmol CO₂ m⁻² s⁻¹ and a coefficient of determination (R²) of 0.73. No significant difference was detected between the prediction accuracy of the SVR model for winter wheat and summer maize. With forcing from satellite remote sensing data and gridded soil property data, we used the SVR model to predict the spatial distributions of R_s during the peak growing seasons of winter wheat and summer maize rotation croplands in Julu County. The SVR model captured the spatial variations of R_s at the county scale. The satellite-derived enhanced vegetation index was found to be the most important input used to predict R_s. Removal of this variable caused an RMSE increase from 0.31 μmol CO₂ m⁻² s⁻¹ to 0.42 μmol CO₂ m⁻² s⁻¹. Soil properties such as soil organic carbon (SOC) content and soil bulk density (SBD) were the second most important factors. Their removal led to an RMSE increase from 0.31 μmol CO₂ m⁻² s⁻¹ to 0.37 μmol CO₂ m⁻² s⁻¹. The SVR model performed better than multiple regression in predicting spatial variations of R_s in winter wheat and summer maize rotation croplands, as shown by the comparison of the R² and RMSE values of the two algorithms. The spatial patterns of R_s are better captured using the SVR model than performing multiple regression, particularly for the relatively high and relatively low R_s values at the center and northeast study areas. Therefore, SVR shows promise for predicting spatial variations of R_s values on the basis of remotely sensed data and gridded soil property data at the county scale.     

Comparison of UAV and WorldView-2 imagery for mapping leaf area index of mangrove forest

Unmanned Aerial Vehicle (UAV) remote sensing has opened the door to new sources of data to effectively characterize vegetation metrics at very high spatial resolution and at flexible revisit frequencies. Successful estimation of the leaf area index (LAI) in precision agriculture with a UAV image has been reported in several studies. However, in most forests, the challenges associated with the interference from a complex background and a variety of vegetation species have hindered research using UAV images. To the best of our knowledge, very few studies have mapped the forest LAI with a UAV image. In addition, the drawbacks and advantages of estimating the forest LAI with UAV and satellite images at high spatial resolution remain a knowledge gap in existing literature. Therefore, this paper aims to map LAI in a mangrove forest with a complex background and a variety of vegetation species using a UAV image and compare it with a WorldView-2 image (WV2).In this study, three representative NDVIs, average NDVI (AvNDVI), vegetated specific NDVI (VsNDVI), and scaled NDVI (ScNDVI), were acquired with UAV and WV2 to predict the plot level (10 × 10 m) LAI. The results showed that AvNDVI achieved the highest accuracy for WV2 (R² = 0.778, RMSE = 0.424), whereas ScNDVI obtained the optimal accuracy for UAV (R² = 0.817, RMSE = 0.423). In addition, an overall comparison results of the WV2 and UAV derived LAIs indicated that UAV obtained a better accuracy than WV2 in the plots that were covered with homogeneous mangrove species or in the low LAI plots, which was because UAV can effectively eliminate the influence from the background and the vegetation species owing to its high spatial resolution. However, WV2 obtained a slightly higher accuracy than UAV in the plots covered with a variety of mangrove species, which was because the UAV sensor provides a negative spectral response function(SRF) than WV2 in terms of the mangrove LAI estimation.     

A practical algorithm for estimating surface soil moisture using combined optical and thermal infrared data

Surface soil moisture (SSM) is a critical variable for understanding the energy and water exchange between the land and atmosphere. A multi-linear model was recently developed to determine SSM using ellipse variables, namely, the center horizontal coordinate (x₀), center vertical coordinate (y₀), semi-major axis (a) and rotation angle (θ), derived from the elliptical relationship between diurnal cycles of land surface temperature (LST) and net surface shortwave radiation (NSSR). However, the multi-linear model has a major disadvantage. The model coefficients are calculated based on simulated data produced by a land surface model simulation that requires sufficient meteorological measurements. This study aims to determine the model coefficients directly using limited meteorological parameters rather than via the complicated simulation process, decreasing the dependence of the model coefficients on meteorological measurements. With the simulated data, a practical algorithm was developed to estimate SSM based on combined optical and thermal infrared data. The results suggest that the proposed approach can be used to determine the coefficients associated with all ellipse variables based on historical meteorological records, whereas the constant term varies daily and can only be determined using the daily maximum solar radiation in a prediction model. Simulated results from three FLUXNET sites over 30 cloud-free days revealed an average root mean square error (RMSE) of 0.042 m³/m³ when historical meteorological records were used to synchronously determine the model coefficients. In addition, estimated SSM values exhibited generally moderate accuracies (coefficient of determination R² = 0.395, RMSE = 0.061 m³/m³) compared to SSM measurements at the Yucheng Comprehensive Experimental Station.     

Canopy foliar nitrogen retrieved from airborne hyperspectral imagery by correcting for canopy structure effects

A statistical relationship between canopy mass-based foliar nitrogen concentration (%N) and canopy bidirectional reflectance factor (BRF) has been repeatedly demonstrated. However, the interaction between leaf properties and canopy structure confounds the estimation of foliar nitrogen. The canopy scattering coefficient (the ratio of BRF and the directional area scattering factor, DASF) has recently been suggested for estimating %N as it suppresses the canopy structural effects on BRF. However, estimation of %N using the scattering coefficient has not yet been investigated for longer spectral wavelengths (>855 nm). We retrieved the canopy scattering coefficient for wavelengths between 400 and 2500 nm from airborne hyperspectral imagery, and then applied a continuous wavelet analysis (CWA) to the scattering coefficient in order to estimate %N. Predictions of %N were also made using partial least squares regression (PLSR). We found that %N can be accurately retrieved using CWA (R² = 0.65, RMSE = 0.33) when four wavelet features are combined, with CWA yielding a more accurate estimation than PLSR (R² = 0.47, RMSE = 0.41). We also found that the wavelet features most sensitive to %N variation in the visible region relate to chlorophyll absorption, while wavelet features in the shortwave infrared regions relate to protein and dry matter absorption. Our results confirm that %N can be retrieved using the scattering coefficient after correcting for canopy structural effect. With the aid of high-fidelity airborne or upcoming space-borne hyperspectral imagery, large-scale foliar nitrogen maps can be generated to improve the modeling of ecosystem processes as well as ecosystem-climate feedbacks.     

A new calibration of the effective scattering albedo and soil roughness parameters in the SMOS SM retrieval algorithm

This study focuses on the calibration of the effective vegetation scattering albedo (ω) and surface soil roughness parameters (H_R, and N_Rp, p = H,V) in the Soil Moisture (SM) retrieval from L-band passive microwave observations using the L-band Microwave Emission of the Biosphere (L-MEB) model. In the current Soil Moisture and Ocean Salinity (SMOS) Level 2 (L2), v620, and Level 3 (L3), v300, SM retrieval algorithms, low vegetated areas are parameterized by ω = 0 and H_R = 0.1, whereas values of ω = 0.06 − 0.08 and H_R = 0.3 are used for forests. Several parameterizations of the vegetation and soil roughness parameters (ω, H_R and N_Rp, p = H,V) were tested in this study, treating SMOS SM retrievals as homogeneous over each pixel instead of retrieving SM over a representative fraction of the pixel, as implemented in the operational SMOS L2 and L3 algorithms. Globally-constant values of ω = 0.10, H_R = 0.4 and N_Rp = −1 (p = H,V) were found to yield SM retrievals that compared best with in situ SM data measured at many sites worldwide from the International Soil Moisture Network (ISMN). The calibration was repeated for collections of in situ sites classified in different land cover categories based on the International Geosphere-Biosphere Programme (IGBP) scheme. Depending on the IGBP land cover class, values of ω and H_R varied, respectively, in the range 0.08–0.12 and 0.1–0.5. A validation exercise based on in situ measurements confirmed that using either a global or an IGBP-based calibration, there was an improvement in the accuracy of the SM retrievals compared to the SMOS L3 SM product considering all statistical metrics (R = 0.61, bias = −0.019 m³ m⁻³, ubRMSE = 0.062 m³ m⁻³ for the IGBP-based calibration; against R = 0.54, bias = −0.034 m³ m⁻³ and ubRMSE = 0.070 m³ m⁻³ for the SMOS L3 SM product). This result is a key step in the calibration of the roughness and vegetation parameters in the operational SMOS retrieval algorithm. The approach presented here is the core of a new forthcoming SMOS optimized SM product.     

Deriving leaf mass per area (LMA) from foliar reflectance across a variety of plant species using continuous wavelet analysis

Leaf mass per area (LMA), the ratio of leaf dry mass to leaf area, is a trait of central importance to the understanding of plant light capture and carbon gain. It can be estimated from leaf reflectance spectroscopy in the infrared region, by making use of information about the absorption features of dry matter. This study reports on the application of continuous wavelet analysis (CWA) to the estimation of LMA across a wide range of plant species. We compiled a large database of leaf reflectance spectra acquired within the framework of three independent measurement campaigns (ANGERS, LOPEX and PANAMA) and generated a simulated database using the PROSPECT leaf optical properties model. CWA was applied to the measured and simulated databases to extract wavelet features that correlate with LMA. These features were assessed in terms of predictive capability and robustness while transferring predictive models from the simulated database to the measured database. The assessment was also conducted with two existing spectral indices, namely the Normalized Dry Matter Index (NDMI) and the Normalized Difference index for LMA (ND_LMA).Five common wavelet features were determined from the two databases, which showed significant correlations with LMA (R²: 0.51–0.82, p < 0.0001). The best robustness (R² = 0.74, RMSE = 18.97 g/m² and Bias = 0.12 g/m²) was obtained using a combination of two low-scale features (1639 nm, scale 4) and (2133 nm, scale 5), the first being predominantly important. The transferability of the wavelet-based predictive model to the whole measured database was either better than or comparable to those based on spectral indices. Additionally, only the wavelet-based model showed consistent predictive capabilities among the three measured data sets. In comparison, the models based on spectral indices were sensitive to site-specific data sets. Integrating the ND_LMA spectral index and the two robust wavelet features improved the LMA prediction. One of the bands used by this spectral index, 1368 nm, was located in a strong atmospheric water absorption region and replacing it with the next available band (1340 nm) led to lower predictive accuracies. However, the two wavelet features were not affected by data quality in the atmospheric absorption regions and therefore showed potential for canopy-level investigations. The wavelet approach provides a different perspective into spectral responses to LMA variation than the traditional spectral indices and holds greater promise for implementation with airborne or spaceborne imaging spectroscopy data for mapping canopy foliar dry biomass.     

Scaling estimates of vegetation structure in Amazonian tropical forests using multi-angle MODIS observations

Detailed knowledge of vegetation structure is required for accurate modelling of terrestrial ecosystems, but direct measurements of the three dimensional distribution of canopy elements, for instance from LiDAR, are not widely available. We investigate the potential for modelling vegetation roughness, a key parameter for climatological models, from directional scattering of visible and near-infrared (NIR) reflectance acquired from NASA’s Moderate Resolution Imaging Spectroradiometer (MODIS). We compare our estimates across different tropical forest types to independent measures obtained from: (1) airborne laser scanning (ALS), (2) spaceborne Geoscience Laser Altimeter System (GLAS)/ICESat, and (3) the spaceborne SeaWinds/QSCAT. Our results showed linear correlation between MODIS-derived anisotropy to ALS-derived entropy (r² = 0.54, RMSE = 0.11), even in high biomass regions. Significant relationships were also obtained between MODIS-derived anisotropy and GLAS-derived entropy (0.52 ≤ r² ≤ 0.61; p < 0.05), with similar slopes and offsets found throughout the season, and RMSE between 0.26 and 0.30 (units of entropy). The relationships between the MODIS-derived anisotropy and backscattering measurements (σ⁰) from SeaWinds/QuikSCAT presented an r² of 0.59 and a RMSE of 0.11. We conclude that multi-angular MODIS observations are suitable to extrapolate measures of canopy entropy across different forest types, providing additional estimates of vegetation structure in the Amazon.