Neural network estimation of chlorophyll a from MERIS full resolution data for the coastal waters of Galician rias (NW Spain)

In typical Case 2 waters, accurate remote sensing retrieval of chlorophyll a (chla) is still a challenging task. In this study, focusing on the Galician rias (ΝW Spain), algorithms based on neural network (NN) techniques were developed for the retrieval of chla concentration in optically complex waters, using Medium Resolution Imaging Spectrometer (MERIS) data. There is considerable interest in the accurate estimation of chla for the Galician rias, because of the economic and social importance of the extensive culture of mussels, and the high frequency of harmful algal events. Fifteen MERIS full resolution (FR) cloud-free images paired with in situ chla data (for 2002-2004 and 2006-2008) were used for the development and validation of the NN. The scope of NN was established from the clusters obtained using fuzzy c-mean (FCM) clustering techniques applied to the satellite-derived data. Three different NNs were developed: one including the whole data set, and two others using only points belonging to one of the clusters. The input data for these latter two NNs was chosen depending on the quality level, defined on the basis of quality flags given to each data set. The fitting results were fairly good and proved the capability of the tool to predict chla concentrations in the study area. The best prediction was given for the NN trained with high-quality data using the most abundant cluster data set. The performance parameters in the validation set of this NN were R² = 0.86, mean percentage error (MPE) = − 0.14, root mean square error (RMSE) = 0.75 mg m^− 3, and relative RMSE = 66%. The NN developed in this study detected accurately the peaks of chla, in both training and validation sets. The performance of the Case-2-Regional (C2R) algorithm, routinely used for MERIS data, was also tested and compared with our best performing NN and the sea-truthing data. Results showed that this NN outperformed the C2R, giving much higher R² and lower RMSE values.This study showed that the combination of in situ data and NN technology improved the retrieval of chla in Case 2 waters, and could be used to obtain more accurate chla maps. A local-based algorithm for the chla retrieval from an ocean colour sensor with the characteristics of MERIS would be a great support in the quantitative monitoring and study of harmful algal events in the coastal waters of the Rias Baixas. The limitations and possible improvements of the developed chla algorithms are also discussed.     

Towards a practical remote-sensing model of suspended sediment concentrations in turbid waters using MERIS measurements

A new empirical index, termed the normalized suspended sediment index (NSSI), is proposed to predict total suspended sediment (TSS) concentrations in inland turbid waters using Medium Resolution Imaging Spectrometer (MERIS) full-resolution (FR) 300 m data. The algorithm is based on the normalized difference between two MERIS spectral bands, 560 and 760 nm. NSSI shows its potential in application to our study region – Poyang Lake – the largest freshwater lake in China. An exponential function (R² = 0.90, p < 0.01) accurately explained the variance in the in situ data and showed better performance for the TSS range 10–524 mg l^?1. The algorithm was then validated with TSS estimates using an atmospheric-corrected MERIS FR image. The validation showed that the NSSI algorithm was a more robust TSS algorithm than the band-ratio algorithms. Findings of this research imply that NSSI can be successfully used on MERIS images to obtain TSS in Poyang Lake. This work provided a practical remote-sensing approach to estimate TSS in the optically and hydrologically complex Poyang Lake and the method can be easily extended to other similar waters.     

Validation of MERIS LAI and FAPAR products for winter wheat-sown test fields in North-East Bulgaria

Progress in deriving land surface biophysical parameters in a spatially explicit manner using remotely sensed data has greatly enhanced our ability to model ecosystem processes and monitor crop development. A multitude of satellite sensors and algorithms have been used to generate ready-to-use maps of various biophysical parameters. Validation of these products for different vegetation types is needed to assess their reliability and consistency. While most of the current satellite biophysical products have spatial resolution of one kilometre, a recent effort utilizing data from the Medium Resolution Imaging Spectrometer (MERIS) provided leaf area index (LAI), Fraction of Absorbed Photosynthetically Active Radiation (FAPAR), and other canopy parameters in a resolution as fine as 300 m over the European continent. This resolution would be more appropriate for application at the regional scale, particularly for crop monitoring. This higher-resolution MERIS product has been evaluated in a limited number of studies to date. This article aims to validate LAI and FAPAR from the MERIS 10-day composite BioPar BP-10 product over winter wheat fields in northeast Bulgaria. The ground measurements of LAI and FAPAR were up-scaled and 30 m resolution reference raster layers were created using empirical relationships with Landsat TM (RMSE = 0.06 and RMSE = 0.22 for FAPAR and LAI, respectively). MERIS FAPAR and LAI were found to have significant correlation with FAPAR and LAI from the reference raster layers (R² = 0.84 and R² = 0.78, respectively). When MERIS Green LAI was calculated (incorporating the fraction of vegetation and brown vegetation cover from the BioPar BP-10 product), better correspondence with LAI values from the reference raster layer was achieved, with RMSE and bias reduced by 30–35%. The results from this study confirm the findings of previous validations showing that MERIS Green LAI tends to overestimate LAI values lower than 1. As a conclusion of the study, the BioPar BP-10 product was found to provide reliable estimates of FAPAR and acceptably accurate estimates of LAI for winter wheat crops in North-East Bulgaria.     

Remote sensing of cyanobacteria-dominant algal blooms and water quality parameters in Zeekoevlei, a small hypertrophic lake, using MERIS

Eutrophication and cyanobacterial algal blooms present an increasing threat to the health of freshwater ecosystems and to humans who use these resources for drinking and recreation. Remote sensing is being used increasingly as a tool for monitoring these phenomena in inland and near-coastal waters. This study uses the Medium Resolution Imaging Spectrometer (MERIS) to view Zeekoevlei, a small hypertrophic freshwater lake situated on the Cape Flats in Cape Town, South Africa, dominated by Microcystis cyanobacteria. The lake's small size, highly turbid water, and covariant water constituents present a challenging case for both algorithm development and atmospheric correction. The objectives of the study are to assess the optical properties of the lake, to evaluate various atmospheric correction procedures, and to compare the performance of empirical and semi-analytical algorithms in hypertrophic water. In situ water quality parameter and radiometric measurements were made simultaneous to MERIS overpasses. Upwelling radiance measurements at depth 0.66 m were corrected for instrument self-shading and processed to water-leaving reflectance using downwelling irradiance measurements and estimates of the vertical attenuation coefficient for upward radiance, K_u, generated from a simple bio-optical model estimating the total absorption, a(λ), and backscattering coefficients, b_b(λ). The normalised water-leaving reflectance was used for assessing the accuracy of image-based Dark Object Subtraction and 6S Radiative Transfer Code atmospheric correction procedures applied to MERIS. Empirical algorithms for estimating chlorophyll a (Chl a), Total Suspended Solids (TSS), Secchi Disk depth (z_SD) and absorption by CDOM (a_CDOM) were derived from simultaneously collected in situ and MERIS measurements. The empirical algorithms gave high correlation coefficient values, although they have a limited ability to separate between signals from covariant water constituents. The MERIS Neural Network algorithms utilised in the standard Level 2 Case 2 waters product and Eutrophic Lakes processor were also used to derive water constituent concentrations. However, these failed to produce reasonable comparisons with in situ measurements owing to the failure of atmospheric correction and divergence between the optical properties and ranges used to train the algorithms and those of Zeekoevlei. Maps produced using the empirical algorithms effectively show the spatial and temporal variability of the water quality parameters during April 2008. On the basis of the results it is argued that MERIS is the current optimal sensor for frequent change detection applications in inland waters. This study also demonstrates the considerable potential value for simple TOA algorithms for hypertrophic systems. It is recommended that regional algorithm development be prioritized in southern Africa and that remote sensing be integrated into future operational water quality monitoring systems.     

Calibration and validation of a generic multisensor algorithm for mapping of total suspended matter in turbid waters

Mapping of total suspended matter concentration (TSM) can be achieved from space-based optical sensors and has growing applications related to sediment transport. A TSM algorithm is developed here for turbid waters, suitable for any ocean colour sensor including MERIS, MODIS and SeaWiFS. Theory shows that use of a single band provides a robust and TSM-sensitive algorithm provided the band is chosen appropriately. Hyperspectral calibration is made using seaborne TSM and reflectance spectra collected in the southern North Sea. Two versions of the algorithm are considered: one which gives directly TSM from reflectance, the other uses the reflectance model of Park and Ruddick (2005) to take account of bidirectional effects.Applying a non-linear regression analysis to the calibration data set gave relative errors in TSM estimation less than 30% in the spectral range 670-750 nm. Validation of this algorithm for MODIS and MERIS retrieved reflectances with concurrent in situ measurements gave the lowest relative errors in TSM estimates, less than 40%, for MODIS bands 667 nm and 678 nm and for MERIS bands 665 nm and 681 nm. Consistency of the approach in a multisensor context (SeaWiFS, MERIS, and MODIS) is demonstrated both for single point time series and for individual images.     

Coupled retrieval of aerosol optical thickness, columnar water vapor and surface reflectance maps from ENVISAT/MERIS data over land

An algorithm for the derivation of atmospheric parameters and surface reflectance data from MEdium Resolution Imaging Specrometer Instrument (MERIS) on board ENVIronmental SATellite (ENVISAT) images has been developed. Geo-rectified aerosol optical thickness (AOT), columnar water vapor (CWV) and spectral surface reflectance maps are generated from MERIS Level-1b data over land. The algorithm has been implemented so that AOT, CWV and reflectance products are provided on an operational manner, making no use of ancillary parameters apart from those attached to MERIS products. For this reason, it has been named Self-Contained Atmospheric Parameters Estimation from MERIS data (SCAPE-M). The fundamental basis of the algorithm and applicable error figures are presented in the first part of this paper. In particular, errors of ± 0.03, ± 4% and ± 8% have been estimated for AOT, CWV and surface reflectance retrievals, respectively, by means of a sensitivity analysis based on a synthetic data set simulated under a usual MERIS scene configuration over land targets. The assumption of a fixed aerosol model, the coarse spatial resolution of the AOT product and the neglection of surface reflectance directional effects were also identified as limitations of SCAPE-M. Validation results are detailed in the second part of the paper. Comparison of SCAPE-M AOT retrievals with data from AErosol RObotic NETwork (AERONET) stations showed an average Root Mean Square Error (RMSE) of 0.05, and an average correlation coefficient R² of about 0.7-0.8. R² values grew up to more than 0.9 in the case of CWV after comparison with the same stations. A good correlation is also found with the MERIS Level-2 ESA CWV product. Retrieved surface reflectance maps have been successfully compared with reflectance data derived from the Compact High Resolution Imaging Spectrometer (CHRIS) on board the PRoject for On-Board Autonomy (PROBA) in the first place. Reflectance retrievals have also been compared with reflectance data derived from MERIS images by the Bremen AErosol Retrieval (BAER) method. A good correlation in the red and near-infrared bands was found, although a considerably higher proportion of pixels was successfully processed by SCAPE-M.     

Lidar-based mapping of leaf area index and its use for validating GLOBCARBON satellite LAI product in a temperate forest of the southern USA

Lidar provides enhanced abilities to remotely map leaf area index (LAI) with improved accuracies. We aim to further explore the capability of discrete-return lidar for estimating LAI over a pine-dominated forest in East Texas, with a secondary goal to compare the lidar-derived LAI map and the GLOBCARBON moderate-resolution satellite LAI product. Specific problems we addressed include (1) evaluating the effects of analysts and algorithms on in-situ LAI estimates from hemispherical photographs (hemiphoto), (2) examining the effectiveness of various lidar metrics, including laser penetration, canopy height and foliage density metrics, to predict LAI, (3) assessing the utility of integrating Quickbird multispectral imagery with lidar for improving the LAI estimate accuracy, and (4) developing a scheme to co-register the lidar and satellite LAI maps and evaluating the consistency between them. Results show that the use of different analysts or algorithms in analyzing hemiphotos caused an average uncertainty of 0.35 in in-situ LAI, and that several laser penetration metrics in logarithm models were more effective than other lidar metrics, with the best one explaining 84% of the variation in the in-situ LAI (RMSE = 0.29 LAI). The selection of plot size and height threshold in calculating laser penetration metrics greatly affected the effectiveness of these metrics. The combined use of NDVI and lidar metrics did not significantly improve estimation over the use of lidar alone. We also found that mis-registration could induce a large artificial discrepancy into the pixelwise comparison between the coarse-resolution satellite and fine-resolution lidar-derived LAI maps. By compensating for a systematic sub-pixel shift error, the correlation between two maps increased from 0.08 to 0.85 for pines (n = 24 pixels). However, the absolute differences between the two LAI maps still remained large due to the inaccuracy in accounting for clumping effects. Overall, our findings imply that lidar offers a superior tool for mapping LAI at local to regional scales as compared to optical remote sensing, accuracies of lidar-estimate LAI are affected not only by the choice of models but also by the absolute accuracy of in-situ reference LAI used for model calibration, and lidar-derived LAI maps can serve as reliable references for validating moderate-resolution satellite LAI products over large areas.     

Comparison of different vegetation indices for the remote assessment of green leaf area index of crops

Many algorithms have been developed for the remote estimation of biophysical characteristics of vegetation, in terms of combinations of spectral bands, derivatives of reflectance spectra, neural networks, inversion of radiative transfer models, and several multi-spectral statistical approaches. However, the most widespread type of algorithm used is the mathematical combination of visible and near-infrared reflectance bands, in the form of spectral vegetation indices. Applications of such vegetation indices have ranged from leaves to the entire globe, but in many instances, their applicability is specific to species, vegetation types or local conditions. The general objective of this study is to evaluate different vegetation indices for the remote estimation of the green leaf area index (Green LAI) of two crop types (maize and soybean) with contrasting canopy architectures and leaf structures. Among the indices tested, the chlorophyll Indices (the CI_Green, the CI_Red-edge and the MERIS Terrestrial Chlorophyll Index, MTCI) exhibited strong and significant linear relationships with Green LAI, and thus were sensitive across the entire range of Green LAI evaluated (i.e., 0.0 to more than 6.0 m²/m²). However, the CI_Red-edge was the only index insensitive to crop type and produced the most accurate estimations of Green LAI in both crops (RMSE = 0.577 m²/m²). These results were obtained using data acquired with close range sensors (i.e., field spectroradiometers mounted 6 m above the canopy) and an aircraft-mounted hyperspectral imaging spectroradiometer (AISA). As the CI_Red-edge also exhibited low sensitivity to soil background effects, it constitutes a simple, yet robust tool for the remote and synoptic estimation of Green LAI. Algorithms based on this index may not require re-parameterization when applied to crops with different canopy architectures and leaf structures, but further studies are required for assessing its applicability in other vegetation types (e.g., forests, grasslands).     

Evaluation of the MODIS LAI product using independent lidar-derived LAI: A case study in mixed conifer forest

This study presents an alternative assessment of the MODIS LAI product for a 58,000 ha evergreen needleleaf forest located in the western Rocky Mountain range in northern Idaho by using lidar data to model (R² = 0.86, RMSE = 0.76) and map LAI at higher resolution across a large number of MODIS pixels in their entirety. Moderate resolution (30 m) lidar-based LAI estimates were aggregated to the resolution of the 1-km MODIS LAI product and compared to temporally-coincident MODIS retrievals. Differences in the MODIS and lidar-derived values of LAI were grouped and analyzed by several different factors, including MODIS retrieval algorithm, sun/sensor geometry, and sub-pixel heterogeneity in both vegetation and terrain characteristics. Of particular interest is the disparity in the results when MODIS LAI was analyzed according to algorithm retrieval class. We observed relatively good agreement between lidar-derived and MODIS LAI values for pixels retrieved with the main RT algorithm without saturation for LAI LAI ≤ 4. Moreover, for the entire range of LAI values, considerable overestimation of LAI (relative to lidar-derived LAI) occurred when either the main RT with saturation or back-up algorithm retrievals were used to populate the composite product regardless of sub-pixel vegetation structural complexity or sun/sensor geometry. These results are significant because algorithm retrievals based on the main radiative transfer algorithm with or without saturation are characterized as suitable for validation and subsequent ecosystem modeling, yet the magnitude of difference appears to be specific to retrieval quality class and vegetation structural characteristics.     

Uncertainty estimates for the FAPAR operational products derived from MERIS — Impact of top-of-atmosphere radiance uncertainties and validation with field data

This paper discusses the accuracy of the operational Medium Resolution Imaging Spectrometer (MERIS) Level 2 land product which corresponds to the Fraction of Absorbed Photosynthetically Active Radiation (FAPAR). The FAPAR value is estimated from daily MERIS spectral measurements acquired at the top-of-atmosphere, using a physically based approach. The products are operationally available at the reduced spatial resolution, i.e. 1.2 km, and can be computed at the full spatial resolution, i.e. at 300 m, from the top-of-atmosphere MERIS data by using the same algorithm. The quality assessment of the MERIS FAPAR products capitalizes on the availability of five years of data acquired globally. The actual validation exercise is performed in two steps including, first, an analysis of the accuracy of the FAPAR algorithm itself with respect to the spectral measurements uncertainties and, second, with a direct comparison of the FAPAR time series against ground-based estimations as well as similar FAPAR products derived from other optical sensor data. The results indicate that the impact of top-of-atmosphere radiance uncertainties on the operational MERIS FAPAR products accuracy is expected to be at about 5-10% and the agreement with the ground-based estimates over different canopy types is achieved within ± 0.1.     

Mapping tree and shrub leaf area indices in an ombrotrophic peatland through multiple endmember spectral unmixing

Leaf area index (LAI) is an important parameter used by most process-oriented ecosystem models. LAI of forest ecosystems has routinely been mapped using spectral vegetation indices (SVI) derived from remote sensing imagery. The application of SVI-based approaches to map LAI in peatlands presents a challenge, mainly due to peatlands characteristic multi-layer canopy comprising shrubs and open, discontinuous tree canopies underlain by a continuous ground cover of different moss species, which reduces the greenness contrast between the canopy and the background.Our goal is to develop a methodology to map tree and shrub LAI in peatlands and similar ecosystems based on multiple endmember spectral mixture analysis (MESMA). This new mapping method is validated using LAI field measurements from a precipitation-fed (ombrotrophic) peatland near Ottawa, Ontario, Canada. We demonstrate first that three commonly applied SVI are not suitable for tree and shrub LAI mapping in ombrotrophic peatlands. Secondly, we demonstrate for a three-endmember model the limitations of traditional linear spectral mixture analysis (SMA) due to the unique and widely varying spectral characteristics of Sphagnum mosses, which are significantly different from vascular plants. Next, by using a geometric-optical radiative transfer model, we determine the nature of the equation describing the empirical relationship between shadow fraction and tree LAI using nonlinear ordinary least square (OLS) regression. We then apply this equation to describe the empirical relationships between shadow and shrub fractions obtained from mixture decomposition with SMA and MESMA, respectively, and tree and shrub LAI, respectively. Less accurate fractions obtained from SMA result in weaker relationships between shadow fraction and tree LAI (R² = 0.61) and shrub fraction and shrub LAI (R² = 0.49) compared to the same relationships based on fractions obtained from MESMA with R² = 0.75 and R² = 0.68, respectively. Cross-validation of tree LAI (R² = 0.74; RMSE = 0.48) and shrub LAI (R² = 0.68; RMSE = 0.42) maps using fractions from MESMA shows the suitability of this approach for mapping tree and shrub LAI in ombrotrophic peatlands. The ability to account for a spectrally varying, unique Sphagnum moss ground cover during mixture decomposition and a two layer canopy is particularly important.     

Utility of an image-based canopy reflectance modeling tool for remote estimation of LAI and leaf chlorophyll content at the field scale

This paper presents a physically-based approach for estimating critical variables describing land surface vegetation canopies, relying on remotely sensed data that can be acquired from operational satellite sensors. The REGularized canopy reFLECtance (REGFLEC) modeling tool couples leaf optics (PROSPECT), canopy reflectance (ACRM), and atmospheric radiative transfer (6SV1) model components, facilitating the direct use of at-sensor radiances in green, red and near-infrared wavelengths for the inverse retrieval of leaf chlorophyll content (C_ab) and total one-sided leaf area per unit ground area (LAI). The inversion of the canopy reflectance model is constrained by assuming limited variability of leaf structure, vegetation clumping, and leaf inclination angle within a given crop field and by exploiting the added radiometric information content of pixels belonging to the same field. A look-up-table with a suite of pre-computed spectral reflectance relationships, each a function of canopy characteristics, soil background effects and external conditions, is accessed for fast pixel-wise biophysical parameter retrievals. Using 1 m resolution aircraft and 10 m resolution SPOT-5 imagery, REGFLEC effectuated robust biophysical parameter retrievals for a corn field characterized by a wide range in leaf chlorophyll levels and intermixed green and senescent leaf material. Validation against in-situ observations yielded relative root-mean-square deviations (RMSD) on the order of 10% for the 1 m resolution LAI (RMSD = 0.25) and C_ab (RMSD = 4.4 μg cm^− 2) estimates, due in part to an efficient correction for background influences. LAI and C_ab retrieval accuracies at the SPOT 10 m resolution were characterized by relative RMSDs of 13% (0.3) and 17% (7.1 μg cm^− 2), respectively, and the overall intra-field pattern in LAI and C_ab was well established at this resolution. The developed method has utility in agricultural fields characterized by widely varying distributions of model variables and holds promise as a valuable operational tool for precision crop management. Work is currently in progress to extend REGFLEC to regional scales.     

Remote sensing of water quality in an Australian tropical freshwater impoundment using matrix inversion and MERIS images

The purpose of this study was to investigate how semi-analytical inversion techniques developed for the remote sensing of water quality parameters (chlorophyll a, tripton and coloured dissolved organic matter (CDOM)) in inland waters could be adapted or improved for application to Australian tropical and sub-tropical water bodies. The Matrix Inversion Method (MIM) with a semi-analytic model of the anisotropy of the in-water light field was applied to MERIS images of Burdekin Falls Dam, Australia, a tropical freshwater impoundment. Specific attention was required to improve the atmospheric correction of the MERIS data. The performance of the conventional three band exact solution of the MIM was compared to that of over-determined solutions that used constant and differential weighting for each sensor band.The results of the application of the MIM algorithm showed that the best weighting scheme had a mean chlorophyll a retrieval difference of 1.0 μgl^− 1, the three band direct matrix inversion scheme had a mean difference of 4.2 μgl^− 1 and the constant weight scheme had a mean difference of 5.5 μgl^− 1. For tripton, the best performed weighting scheme had a mean difference of 1.2 mgl^− 1, the three band scheme had a mean difference of 3.4 mgl^− 1 and the constant weight scheme had a mean difference of 1.8 mgl^− 1. For the CDOM retrieval, the mean difference was found to be 0.12 m^− 1 for the best performed weighting scheme, 0.25 m^− 1 for the three band scheme and 0.52 m^− 1 for the constant weight scheme. It was found that significant improvements in the accuracy and precision of retrieved water quality parameter values can be obtained by using differentially weighted, over-determined systems of equations, rather than exact solutions. These more reliable estimates of water quality parameters will allow water resource managers to improve their monitoring regimes.     

Chlorophyll retrieval with MERIS Case-2-Regional in perialpine lakes

Semi-analytical remote sensing applications for eutrophic waters are not applicable to oligo- and mesotrophic lakes in the perialpine area, since they are insensitive to chlorophyll concentration variations between 1 and 10 mg/m³. The neural network based Case-2-Regional algorithm for MERIS was developed to fill this gap, along with the ICOL adjacency effect correction algorithm. The algorithms are applied to a collection of 239 satellite images from 2003-2008, and the results are compared to experimental and official water quality data collected in six perialpine lakes in the same period. It is shown that remote sensing estimates can provide an adequate supplementary data source to in situ data series of the top 5 m water layer, provided that a sufficient number of matchups for a site specific maximum temporal offset are available.     

Intercomparison of ocean colour band-ratio algorithms for chlorophyll concentration in the Subtropical Front east of New Zealand

This study intercompared the performance of eight band-ratio chlorophyll-a algorithms which together can be used to process measurements from the ocean colour satellite sensors CZCS, OCTS, SeaWiFS, MODIS, MERIS, and GLI. The study area included Subtropical, Subtropical Front and Subantarctic waters east of New Zealand, and Case 1 waters of the New Zealand northeast continental shelf. Over 170 co-incident measurements of spectral normalised water-leaving radiance and near-surface concentration of chlorophyll-a were made on nine research voyages between 1998 and 2000. The studentised bootstrap method was used to identify statistically significant bias in algorithm products relative to in situ measurements. The band-ratio algorithms used by CZCS, OCTS and SeaWiFS missions systematically underestimated chlorophyll-a concentration in the offshore regions by between 21% and 45%, but showed no systematic bias in the continental shelf waters. The band-ratio algorithms applicable to the MODIS and MERIS sensors had no clear bias with respect to in situ measurements in offshore waters, but had a positive bias of 20% over the continental shelf. The proposed GLI band-ratio algorithm led to estimates that were negatively biased with respect to in situ measurement offshore (− 30%), and positively biased over the continental shelf (20%). The results were consistent with unusually high values of absorption in the blue part of the spectrum (443-490 nm) compared to the green part (∼ 550 nm) by phytoplankton pigments in the offshore waters, and high chlorophyll-specific absorption over the continental shelf.     

Remote estimation of chl-a concentration in turbid productive waters — Return to a simple two-band NIR-red model?

Today the water quality of many inland and coastal waters is compromised by cultural eutrophication in consequence of increased human agricultural and industrial activities. Remote sensing is widely applied to monitor the trophic state of these waters. This study investigates the performance of near infrared-red models for the remote estimation of chlorophyll-a concentrations in turbid productive waters and evaluates several near infrared-red models developed within the last 34 years. Three models were calibrated for a dataset with chlorophyll-a concentrations from 0 to 100 mg m⁻³ and validated for independent and statistically different datasets with chlorophyll-a concentrations from 0 to 100 mg m⁻³ and 0 to 25 mg m⁻³ for the spectral bands of the MEdium Resolution Imaging Spectrometer (MERIS) and MODerate resolution Imaging Spectroradiometer (MODIS). The MERIS two-band model estimated chlorophyll-a concentrations slightly more accurately than the more complex models, with mean absolute errors of 2.3 mg m⁻³ for chlorophyll-a concentrations from 0 to 100 mg m⁻³ and 1.2 mg m⁻³ for chlorophyll-a concentrations from 0 to 25 mg m⁻³. Comparable results from several near infrared-red models with different levels of complexity, calibrated for inland and coastal waters around the world, indicate a high potential for the development of a simple universally applicable near infrared-red algorithm.     

Evaluation of compositing period and AVHRR and MERIS combination for improvement of spring phenology detection in deciduous forests

Accurate and precise detection of phenology events is needed to assess trends in seasonal vegetation development indicative of climate or other environmental change processes. In this research, detection accuracy of start of season (SOS) phenology for deciduous forest across Eastern Canada was assessed using satellite time series and in situ PlantWatch observations. Several aspects were evaluated regarding performance of phenology information extraction: 1) effect of compositing period, 2) individual performance of the Advanced Very High Resolution Radiometer (AVHRR) and the Medium Resolution Imaging Spectrometer (MERIS) sensors, and 3) performance for these sensors combined. The AVHRR and MERIS sensors were used as they are overlapping operational missions with planned future continuity. Three approaches to utilizing the multi-sensor data were tested: 1) inter-calibrating NDVI data between sensors and using the multi-sensor data stream to detect SOS, 2) combining independently derived SOS estimates from AVHRR and MERIS based on a weighted average, and 3) combining approaches 1 and 2. Comparison with in situ observations of leaf out and first bloom showed that combining independent SOS estimates from AVHRR and MERIS was better than using the inter-calibrated multi-sensor data. Combining SOS estimates from both sensors reduced error by 1-2 days compared to the individual sensor results. Composite periods from 7 to 11 days produced the best results for leaf out with a mean absolute error (MAE) of 5 days. Results for first bloom were not as good as those for leaf out, producing a MAE of 6.5 days. For first bloom, compositing periods greater than 11 days did not increase error at the same rate as seen for leaf out. However, the larger MAE observed for first bloom may have masked this effect.     

Multitemporal-patch ensemble inversion of coupled surface-atmosphere radiative transfer models for land surface characterization

In the last few years, encouraging results using radiative transfer model inversion techniques were obtained for land biophysical variables retrieval. However, the inversion of radiative transfer models is a severely ill-posed problem that may lead to significant uncertainties in the biophysical variables estimates. Improvement of performances of the inversion process requires more information to be exploited including better radiative transfer models, exploitation of proper prior information on the distribution of canopy and atmosphere variables, knowledge of uncertainties in satellite measurements, as well as possible spatial and temporal constraints. In this study we focus on the use of coupled atmosphere-surface radiative transfer models (SMAC + SAIL + PROSPECT) to estimate some key biophysical variables from top of atmosphere canopy reflectance data. The inversion is achieved over an ensemble of pixels belonging to a spatial window where aerosol properties are supposed to be constant, and over a temporal window of few days where vegetation state is assumed not to vary. The ensemble inversion scheme accounting for the spatial and temporal constraints is described. Top of atmosphere reflectance observations are simulated for 13 bands within the visible and near infrared domains. The coupled model is inverted with a variational method implementation dedicated to solve very large inverse problems. It is based on the use of the adjoint model and a Quasi-Newton optimisation technique with BFGS update.The multitemporal-patch inversion approach exploiting the spatial and temporal constraints is compared to the classical instantaneous-local inversion applied on single pixel and date. The ‘ensemble’ approach shows significant performance improvements when retrieving aerosol optical thickness τ₅₅₀ and some canopy characteristics (LAI, LAI × C_ab and ALA). Conclusions are drawn on the interest of such approaches, and perspectives are given, with due attention to their applicability within operational algorithms.     

Estimation of instantaneous net surface longwave radiation from MODIS cloud-free data

This paper develops a statistical regression method to estimate the instantaneous Downwelling Surface Longwave Radiation (DSLR) for cloud-free skies using only the satellite-based radiances measured at the Top Of the Atmosphere (TOA), and subsequently combines the DSLR with the MODIS land surface temperature/emissivity products (MOD11_L2) to estimate the instantaneous Net Surface Longwave Radiation (NSLR). The proposed method relates the DSLR directly to the TOA radiances in the MODIS Thermal InfraRed (TIR) channels provided that the terrain altitude and the satellite Viewing Zenith Angle (VZA) are known. The simulation analysis shows that the instantaneous DSLR could be estimated by the proposed method with the Root Mean Square Error (RMSE) of 12.4 W/m² for VZA = 0 and terrain altitude z = 0 km. Similar results are obtained for the other VZAs and altitudes. Considering the MODIS instrumental errors of 0.25 K for the TOA brightness temperatures in channels 28, 33 and 34, and of 0.05 K for channels 29 and 31, and of 0.35 K for channel 36, the overall retrieval accuracy in terms of the RMSE is decreased to 13.1 W/m² for the instantaneous DSLR. Moreover, a comparison of MODIS derived DSLR and NSLR are done with the field measurements made at six sites of the Surface Radiation Budget Network (SURFRAD) in the United States for days with cloud-free conditions at the moment of MODIS overpass in 2006. The results show that the bias, RMSE and the square of the correlation coefficient (R²) between the MODIS derived DSLR with the proposed method and the field measured DSLR are 20.3 W/m², 30.1 W/m² and 0.91 respectively, and bias = 11.7 W/m², RMSE = 26.1 W/m² and R² = 0.94 for NSLR. In addition, the scheme proposed by Bisht et al. [Bisht, G., Venturini, V., Islam, S., & Jiang, L. (2005). Estimation of the net radiation using MODIS (Moderate Resolution Imaging Spectroradiometer) data for clear-sky days. Remote Sensing of Environment, 97, 52-67], which requires the MODIS atmospheric profile product (MOD07) and also the MODIS land surface temperature/emissivity products (MOD11_L2) as inputs, is used to estimate the instantaneous DSLR and NSLR for comparison with the field measurements as well as the MODIS derived DSLR and NSLR using our proposed method. The results of the comparisons show that, at least for our cases, our proposed method for estimating DSLR from the MODIS radiances at the TOA and the resultant NSLR gives results comparable to those estimated with Bisht et al.'s scheme [Bisht, G., Venturini, V., Islam, S., & Jiang, L. (2005). Estimation of the net radiation using MODIS (Moderate Resolution Imaging Spectroradiometer) data for clear-sky days. Remote Sensing of Environment, 97, 52-67].