Assessment of satellite ocean color products of MERIS,MODIS and SeaWiFS along the East China Coast (in the Yellow Sea and East China Sea)

The validation of satellite ocean-color products is an important task of ocean-color missions. The uncertainties of these products are poorly quantified in the Yellow Sea (YS) and East China Sea (ECS), which are well known for their optical complexity and turbidity in terms of both oceanic and atmospheric optical properties. The objective of this paper is to evaluate the primary ocean-color products from three major ocean-color satellites, namely the Moderate Resolution Imaging Spectroradiometer (MODIS), Medium Resolution Imaging Spectrometer (MERIS), and Sea-viewing Wide Field-of-view Sensor (SeaWiFS). Through match-up analysis with in situ data, it is found that satellite retrievals of the spectral remote sensing reflectance R_rs(λ) at the blue-green and green bands from MERIS, MODIS and SeaWiFS have the lowest uncertainties with a median of the absolute percentage of difference (APD_m) of 15–27% and root-mean-square-error (RMS) of 0.0021–0.0039 sr⁻¹, whereas the R_rs(λ) uncertainty at 412 nm is the highest (APD_m 47–62%, RMS 0.0027–0.0041 sr⁻¹). The uncertainties of the aerosol optical thickness (AOT) τ_a, diffuse attenuation coefficient for downward irradiance at 490 nm K_d(490), concentrations of suspended particulate sediment concentration (SPM) and Chlorophyll a (Chl-a) were also quantified. It is demonstrated that with appropriate in-water algorithms specifically developed for turbid waters rather than the standard ones adopted in the operational satellite data processing chain, the uncertainties of satellite-derived properties of K_d(490), SPM, and Chl-a may decrease significantly to the level of 20–30%, which is true for the majority of the study area. This validation activity advocates for (1) the improvement of the atmosphere correction algorithms with the regional aerosol optical model, (2) switching to regional in-water algorithms over turbid coastal waters, and (3) continuous support of the dedicated in situ data collection effort for the validation task.     

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.     

A visible band index for remote sensing leaf chlorophyll content at the canopy scale

Leaf chlorophyll content is an important variable for agricultural remote sensing because of its close relationship to leaf nitrogen content. The triangular greenness index (TGI) was developed based on the area of a triangle surrounding the spectral features of chlorophyll with points at (670 nm, R₆₇₀), (550 nm, R₅₅₀), and (480 nm, R₄₈₀), where R_λ is the spectral reflectance at wavelengths of 670, 550 and 480, respectively. The equation is TGI = −0.5[(670 − 480)(R₆₇₀ − R₅₅₀) − (670 − 550)(R₆₇₀ − R₄₈₀)]. In 1999, investigators funded by NASA's Earth Observations Commercialization and Applications Program collaborated on a nitrogen fertilization experiment with irrigated maize in Nebraska. Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data and Landsat 5 Thematic Mapper (TM) data were acquired along with leaf chlorophyll meter and other data on three dates in July during late vegetative growth and early reproductive growth. TGI was consistently correlated with plot-averaged chlorophyll-meter values at the spectral resolutions of AVIRIS, Landsat TM, and digital cameras. Simulations using the Scattering by Arbitrarily Inclined Leaves (SAIL) canopy model indicate an interaction among TGI, leaf area index (LAI) and soil type at low crop LAI, whereas at high LAI and canopy closure, TGI was only affected by leaf chlorophyll content. Therefore, TGI may be the best spectral index to detect crop nitrogen requirements with low-cost digital cameras mounted on low-altitude airborne platforms.     

Modelling high arctic percent vegetation cover using field digital images and high resolution satellite data

In this study, digital images collected at a study site in the Canadian High Arctic were processed and classified to examine the spatial-temporal patterns of percent vegetation cover (PVC). To obtain the PVC of different plant functional groups (i.e., forbs, graminoids/sedges and mosses), field near infrared-green-blue (NGB) digital images were classified using an object-based image analysis (OBIA) approach. The PVC analyses comparing different vegetation types confirmed: (i) the polar semi-desert exhibited the lowest PVC with a large proportion of bare soil/rock cover; (ii) the mesic tundra cover consisted of approximately 60% mosses; and (iii) the wet sedge consisted almost exclusively of graminoids and sedges. As expected, the PVC and green normalized difference vegetation index (GNDVI; (R_NIR − R_Green)/(R_NIR + R_Green)), derived from field NGB digital images, increased during the summer growing season for each vegetation type: i.e., ∼5% (0.01) for polar semi-desert; ∼10% (0.04) for mesic tundra; and ∼12% (0.03) for wet sedge respectively. PVC derived from field images was found to be strongly correlated with WorldView-2 derived normalized difference spectral indices (NDSI; (R_x − R_y)/(R_x + R_y)), where R_x is the reflectance of the red edge (724.1 nm) or near infrared (832.9 nm and 949.3 nm) bands; R_y is the reflectance of the yellow (607.7 nm) or red (658.8 nm) bands with R²’s ranging from 0.74 to 0.81. NDSIs that incorporated the yellow band (607.7 nm) performed slightly better than the NDSIs without, indicating that this band may be more useful for investigating Arctic vegetation that often includes large proportions of senescent vegetation throughout the growing season.     

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.     

Retrieval of tea polyphenol at leaf level using spectral transformation and multi-variate statistical approach

In the present study, field based hyperspectral data was used to estimate the tea (Camellia sinensis L.) polyphenol at Deha Tea garden of Assam state, India. Leaf reflectance spectra were first filtered for noise and then transformed into normalized and first derivative reflectance for further analysis. Stepwise discriminant analysis was carried out to select sensitive bands for a range of polyphenol concentration by minimizing the effects of other factors such as age of the bushes and management practices. The wavelengths at 358, 369, 484, 845, 916, 1387, 1420, 1435, 1621 and 2294 nm were identified as sensitive to tea polyphenol, among which 2294 nm was found to be the most recurring band. The noise removed selected bands, their transformed derivatives and principal components were regressed with the tea polyphenol using univariate and multi-variate analysis. In univariate analysis the correlation was very poor with RMSE more than 3.0. A significant improvement in R² values were observed when multivariate analyses like stepwise multiple linear regression (SMLR) and partial least square regression (PLSR) was carried out. The PLSR of first derivative reflectance was most accurate (R² = 0.81 and RMSE = 1.39 mg g⁻¹) among all the uni- and multivariate analysis for predicting the polyphenol of fresh tea leaves.     

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.     

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.     

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.     

Remote quantification of phycocyanin in potable water sources through an adaptive model

Cyanobacterial blooms in water supply sources in both central Indiana USA (CIN) and South Australia (SA) are a cause of great concerns for toxin production and water quality deterioration. Remote sensing provides an effective approach for quick assessment of cyanobacteria through quantification of phycocyanin (PC) concentration. In total, 363 samples spanning a large variation of optically active constituents (OACs) in CIN and SA waters were collected during 24 field surveys. Concurrently, remote sensing reflectance spectra (R_rs) were measured. A partial least squares–artificial neural network (PLS–ANN) model, artificial neural network (ANN) and three-band model (TBM) were developed or tuned by relating the R_rs with PC concentration. Our results indicate that the PLS–ANN model outperformed the ANN and TBM with both the original spectra and simulated ESA/Sentinel-3/Ocean and Land Color Instrument (OLCI) and EO-1/Hyperion spectra. The PLS–ANN model resulted in a high coefficient of determination (R²) for CIN dataset (R² = 0.92, R: 0.3–220.7 μg/L) and SA (R² = 0.98, R: 0.2–13.2 μg/L). In comparison, the TBM model yielded an R² = 0.77 and 0.94 for the CIN and SA datasets, respectively; while the ANN obtained an intermediate modeling accuracy (CIN: R² = 0.86; SA: R² = 0.95). Applying the simulated OLCI and Hyperion aggregated datasets, the PLS–ANN model still achieved good performance (OLCI: R² = 0.84; Hyperion: R² = 0.90); the TBM also presented acceptable performance for PC estimations (OLCI: R² = 0.65, Hyperion: R² = 0.70). Based on the results, the PLS–ANN is an effective modeling approach for the quantification of PC in productive water supplies based on its effectiveness in solving the non-linearity of PC with other OACs. Furthermore, our investigation indicates that the ratio of inorganic suspended matter (ISM) to PC concentration has close relationship to modeling relative errors (CIN: R² = 0.81; SA: R² = 0.92), indicating that ISM concentration exert significant impact on PC estimation accuracy.     

Hyperspectral-based predictive modelling of grapevine water status in the Portuguese Douro wine region

In this study, hyperspectral reflectance (HySR) data derived from a handheld spectroradiometer were used to assess the water status of three grapevine cultivars in two sub-regions of Douro wine region during two consecutive years. A large set of potential predictors derived from the HySR data were considered for modelling/predicting the predawn leaf water potential (Ψ_pd) through different statistical and machine learning techniques. Three HySR vegetation indices were selected as final predictors for the computation of the models and the in-season time trend was removed from data by using a time predictor. The vegetation indices selected were the Normalized Reflectance Index for the wavelengths 554 nm and 561 nm (NRI_554;561), the water index (WI) for the wavelengths 900 nm and 970 nm, and the D1 index which is associated with the rate of reflectance increase in the wavelengths of 706 nm and 730 nm. These vegetation indices covered the green, red edge and the near infrared domains of the electromagnetic spectrum. A large set of state-of-the-art analysis and statistical and machine-learning modelling techniques were tested. Predictive modelling techniques based on generalized boosted model (GBM), bagged multivariate adaptive regression splines (B-MARS), generalized additive model (GAM), and Bayesian regularized neural networks (BRNN) showed the best performance for predicting Ψ_pd, with an average determination coefficient (R²) ranging between 0.78 and 0.80 and RMSE varying between 0.11 and 0.12 MPa. When cultivar Touriga Nacional was used for training the models and the cultivars Touriga Franca and Tinta Barroca for testing (independent validation), the models performance was good, particularly for GBM (R² = 0.85; RMSE = 0.09 MPa). Additionally, the comparison of Ψ_pd observed and predicted showed an equitable dispersion of data from the various cultivars. The results achieved show a good potential of these predictive models based on vegetation indices to support irrigation scheduling in vineyard.     

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.     

Vegetation chlorophyll estimates in the Amazon from multi-angle MODIS observations and canopy reflectance model

As a preparatory study for future hyperspectral missions that can measure canopy chemistry, we introduce a novel approach to investigate whether multi-angle Moderate Resolution Imaging Spectroradiometer (MODIS) data can be used to generate a preliminary database with long-term estimates of chlorophyll. MODIS monthly chlorophyll estimates between 2000 and 2015, derived from a fully coupled canopy reflectance model (ProSAIL), were inspected for consistency with eddy covariance fluxes, tower-based hyperspectral images and chlorophyll measurements. MODIS chlorophyll estimates from the inverse model showed strong seasonal variations across two flux-tower sites in central and eastern Amazon. Marked increases in chlorophyll concentrations were observed during the early dry season. Remotely sensed chlorophyll concentrations were correlated to field measurements (r² = 0.73 and r² = 0.98) but the data deviated from the 1:1 line with root mean square errors (RMSE) ranging from 0.355 μg cm⁻² (Tapajós tower) to 0.470 μg cm⁻² (Manaus tower). The chlorophyll estimates were consistent with flux tower measurements of photosynthetically active radiation (PAR) and net ecosystem productivity (NEP). We also applied ProSAIL to mono-angle hyperspectral observations from a camera installed on a tower to scale modeled chlorophyll pigments to MODIS observations (r² = 0.73). Chlorophyll pigment concentrations (Chl_A+B) were correlated to changes in the amount of young and mature leaf area per month (0.59 ≤  r² ≤  0.64). Increases in MODIS observed Chl_A+B were preceded by increased PAR during the dry season (0.61 ≤ r² ≤  0.62) and followed by changes in net carbon uptake. We conclude that, at these two sites, changes in LAI, coupled with changes in leaf chlorophyll, are comparable with seasonality of plant productivity. Our results allowed the preliminary development of a 15-year time series of chlorophyll estimates over the Amazon to support canopy chemistry studies using future hyperspectral sensors.     

Estimating dissolved organic carbon concentration in turbid coastal waters using optical remote sensing observations

Dissolved Organic Carbon (DOC) is an important component in the global carbon cycle. It also plays an important role in influencing the coastal ocean biogeochemical (BGC) cycles and light environment. Studies focussing on DOC dynamics in coastal waters are data constrained due to the high costs associated with in situ water sampling campaigns. Satellite optical remote sensing has the potential to provide continuous, cost-effective DOC estimates. In this study we used a bio-optics dataset collected in turbid coastal waters of Moreton Bay (MB), Australia, during 2011 to develop a remote sensing algorithm to estimate DOC. This dataset includes data from flood and non-flood conditions. In MB, DOC concentration varied over a wide range (20–520 μM C) and had a good correlation (R² = 0.78) with absorption due to coloured dissolved organic matter (CDOM) and remote sensing reflectance. Using this data set we developed an empirical algorithm to derive DOC concentrations from the ratio of R_rs(412)/R_rs(488) and tested it with independent datasets. In this study, we demonstrate the ability to estimate DOC using remotely sensed optical observations in turbid coastal waters.     

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.     

Evaluating suitability of MODIS-Terra images for reproducing historic sediment concentrations in water bodies: Lake Tana,Ethiopia

Government and NGO funded conservation programs are being implemented in developing countries with the potential benefit of reduced sediment inflow into fresh water lakes. However, these claims are difficult to verify due to limited historical sediment concentration data in lakes and rivers. Remote sensing can potentially aid in monitoring sediment concentration. With almost daily availability over the past ten years and consistent atmospheric correction applied to the images, Moderate Resolution Imaging Spectroradiometer (MODIS) 250 meter images are potential resources capable of monitoring future concentrations and reconstructing historical sediment concentration records. In this paper, site-specific relationships are developed between reflectance in near-infrared (NIR) images and three factors: total suspended solids (TSS), turbidity and Secchi depth for Lake Tana near the mouth of the Gumara River. The first two sampling campaigns on November 27, 2010 and May 13, 2011 are used in calibration. Reflectance in the NIR varies linearly with turbidity (R² = 0.89) and TSS (R² = 0.95). Secchi depth fit best to an exponential relation with R² of 0.74. The relationships are validated using a third sample set collected on November 7, 2011 with RMSE of 11 Nephelometric Turbidity Units (NTU) for Turbidity, 16.5 mg l⁻¹ for TSS and 0.12 meters for Secchi depth. The MAE was 10% for TSS, 14% for turbidity and 0.1% for Secchi depth. Using the relationship for TSS, a 10-year time series of sediment concentration in Lake Tana near the Gumara River was plotted. It was found that after the severe drought of 2002 and 2003 the concentration in the lake increased significantly. The results showed that MODIS images are potential cost effective tools to monitor suspended sediment concentration and obtain a past history of concentration for evaluating the effect of best management practices.     

Evaluation of a three-band model for estimating chlorophyll-a concentration in tidal reaches of the Pearl River Estuary, China

Accurate assessment of phytoplankton chlorophyll-a (Chla) concentration in turbid waters by means of remote sensing was challenging due to the optical complexity of turbid waters. Recently, a conceptual model containing reflectance in three spectral bands in the red and near-infrared range of the spectrum was suggested for retrieving Chla concentrations in turbid productive waters. The objective of this paper was to evaluate the performance of this three-band model to estimate Chla concentration in the Pearl River Estuary (PRE), China. Reflectance spectra of surface water and water samples were collected concurrently. The samples contained variable Chla (4.80-92.60 mg/m³) and total suspended solids (0.4-55.2 mg/L dry wt). Colored dissolved organic matter (CDOM) absorption at 400 nm was 0.40-1.41 m⁻¹; turbidity ranged from 4 to 25 NTU (Nephelometric Turbidity Units). The three-band model was spectrally calibrated by iterative and least-square linear regression methods to select the optimal spectral bands for the most accurate Chla estimation. Strong linear relationships (R²=0.81, RMSE=1.4 mg/m³, N=32) were established between measured Chla and the levels obtained from the calibrated three-band model [R⁻¹(684)-R⁻¹(690)]×R(718), where R(λ) was the reflectance at wavelength λ. The calibrated three-band model was independently validated (R²=0.9521, RMSE=6.44 mg/m³, N=16) and applied to retrieve Chla concentrations from the calibrated EO-1 Hyperion reflectance data in the PRE on December 21, 2006. The EO-1 Hyperion-derived Chla concentrations were further validated using synchronous in situ data collected on the same day (R²=0.64, RMSE=2 mg/m³, N=9). The spatial tendency of Chla distribution mapping by Hyperion showed gradually increased concentrations of Chla farther from the river mouths (although decreasing from east to west), which were disturbed by the combination of river outlets and tidal current in Lingding Bay of the PRE. This observation conformed to previous observations and studies, and could reasonably be explained by geographical changes. Also, results indicated that the slope of the three-band regression line decreased as the Chla concentration increased, resulting in the first sensitive band of the three-band model to move towards short wavelengths. These findings validated the rationale behind the conceptual model and demonstrated the robustness of this algorithm for Chla retrieval from in situ data and the Hyperion satellite sensor in turbid estuarine waters of the PRE, China.     

Understanding the optical responses of leaf nitrogen in Mediterranean Holm oak (Quercus ilex) using field spectroscopy

The direct estimation of nitrogen (N) in fresh vegetation is challenging due to its weak influence on leaf reflectance and the overlaps with absorption features of other compounds. Different empirical models relate in this work leaf nitrogen concentration ([N]_Leaf) on Holm oak to leaf reflectance as well as derived spectral indices such as normalized difference indices (NDIs), the three bands indices (TBIs) and indices previously used to predict leaf N and chlorophyll. The models were calibrated and assessed their accuracy, robustness and the strength of relationship when other biochemicals were considered. Red edge was the spectral region most strongly correlated with [N]_Leaf, whereas most of the published spectral indexes did not provide accurate estimations. NDIs and TBIs based models could achieve robust and acceptable accuracies (TBI_{1310,1720,730}: R² = 0.76, [0.64,0.86]; RMSE (%) = 9.36, [7.04,12.83]). These models sometimes included indices with bands close to absorption features of N bonds or nitrogenous compounds, but also of other biochemicals. Models were independently and inter-annually validated using the bootstrap method, which allowed discarding those models non-robust across different years. Partial correlation analysis revealed that spectral estimators did not strongly respond to [N]_Leaf but to other leaf variables such as chlorophyll and water, even if bands close to absorption features of N bonds or compounds were present in the models.     

Predicting grain protein content of winter wheat using remote sensing data based on nitrogen status and water stress

Advanced site-specific knowledge of grain protein content of winter wheat from remote sensing data would provide opportunities to manage grain harvest differently, and to maximize output by adjusting input in fields. In this study, remote sensing data were utilized to predict grain protein content. Firstly, the leaf nitrogen content at winter wheat anthesis stage was proved to be significantly correlated with grain protein content (R² = 0.36), and spectral indices significantly correlated to leaf nitrogen content at anthesis stage were potential indicators for grain protein content. The vegetation index, VI_green, derived from the canopy spectral reflectance at green and red bands, was significantly correlated to the leaf nitrogen content at anthesis stage, and also highly significantly correlated to the final grain protein content (R² = 0.46). Secondly, the external conditions, such as irrigation, fertilization and temperature, had important influence on grain quality. Water stress at grain filling stage can increase grain protein content, and leaf water content is closely related to irrigation levels, therefore, the spectral indices correlated to leaf water content can be potential indicators for grain protein content. The spectral reflectance of TM channel 5 derived from canopy spectra or image data at grain filling stage was all significantly correlated to grain protein content (R² = 0.31 and 0.37, respectively). Finally, not only this study proved the feasibility of using remote sensing data to predict grain protein content, but it also provided a tentative prediction of the grain protein content in Beijing area using the reflectance image of TM channel 5.