Integrated analysis of PALSAR/Radarsat-1 InSAR and ENVISAT altimeter data for mapping of absolute water level changes in Louisiana wetlands

Interferometric Synthetic Aperture Radar (InSAR) has been used to detect relative water level changes in wetlands. We developed an innovative method to integrate InSAR and satellite radar altimetry for measuring absolute or geocentric water level changes and applied the methodology to remote areas of swamp forest in coastal Louisiana. Coherence analysis of InSAR pairs suggested that the HH polarization is preferred for this type of observation, and polarimetric analysis can help to identify double-bounce backscattering areas in the wetland. ENVISAT radar altimeter-measured 18-Hz (along-track sampling of 417 m) water level data processed with regional stackfile method have been used to provide vertical references for water bodies separated by levees. The high-resolution (~ 40 m) relative water changes measured from ALOS PALSAR L-band and Radarsat-1 C-band InSAR are then integrated with ENVISAT radar altimetry to obtain absolute water level. The resulting water level time series were validated with in situ gauge observations within the swamp forest. We anticipate that this new technique will allow retrospective reconstruction and concurrent monitoring of water conditions and flow dynamics in wetlands, especially those lacking gauge networks.     

Preliminary results of ENVISAT RA-2-derived water levels validation over the Amazon basin

Since the launch of the ENVISAT satellite in 2002, the Radar Altimetry Mission provides systematic observations of the Earth topography. Among the different goals of the ENVISAT Mission, one directly concerns land hydrology: the monitoring of the water levels of lakes, wetlands, and rivers. The ENVISAT Geophysical Data Records products contain, over different type of surfaces, altimeter ranges derived from four specialized algorithms or retrackers. However, none of the retrackers are intended to the processing of the radar echoes over continental waters. A validation study is necessary to assess the performances of the different ENVISAT-derived water levels to monitor inland waters. We have selected four test-zones over the Amazon basin to achieve this validation study. We compare first the performances of these retracking algorithms to deliver reliable water levels for land hydrology. Comparisons with in-situ gauge stations showed that Ice-1 algorithm, based on the Offset Centre of Gravity technique, provides the more accurate water stages. Second, we examine the potentiality to combine water levels derived from different sensors (Topex/Poseidon, ERS-1 and -2, GFO).     

Impact of footprint diameter and off-nadir pointing on the precision of canopy height estimates from spaceborne lidar

A spaceborne lidar mission could serve multiple scientific purposes including remote sensing of ecosystem structure, carbon storage, terrestrial topography and ice sheet monitoring. The measurement requirements of these different goals will require compromises in sensor design. Footprint diameters that would be larger than optimal for vegetation studies have been proposed. Some spaceborne lidar mission designs include the possibility that a lidar sensor would share a platform with another sensor, which might require off-nadir pointing at angles of up to 16°. To resolve multiple mission goals and sensor requirements, detailed knowledge of the sensitivity of sensor performance to these aspects of mission design is required.This research used a radiative transfer model to investigate the sensitivity of forest height estimates to footprint diameter, off-nadir pointing and their interaction over a range of forest canopy properties. An individual-based forest model was used to simulate stands of mixed conifer forest in the Tahoe National Forest (Northern California, USA) and stands of deciduous forests in the Bartlett Experimental Forest (New Hampshire, USA). Waveforms were simulated for stands generated by a forest succession model using footprint diameters of 20 m to 70 m. Off-nadir angles of 0 to 16° were considered for a 25 m diameter footprint diameter.Footprint diameters in the range of 25 m to 30 m were optimal for estimates of maximum forest height (R² of 0.95 and RMSE of 3 m). As expected, the contribution of vegetation height to the vertical extent of the waveform decreased with larger footprints, while the contribution of terrain slope increased. Precision of estimates decreased with an increasing off-nadir pointing angle, but off-nadir pointing had less impact on height estimates in deciduous forests than in coniferous forests. When pointing off-nadir, the decrease in precision was dependent on local incidence angle (the angle between the off-nadir beam and a line normal to the terrain surface) which is dependent on the off-nadir pointing angle, terrain slope, and the difference between the laser pointing azimuth and terrain aspect; the effect was larger when the sensor was aligned with the terrain azimuth but when aspect and azimuth are opposed, there was virtually no effect on R² or RMSE. A second effect of off-nadir pointing is that the laser beam will intersect individual crowns and the canopy as a whole from a different angle which had a distinct effect on the precision of lidar estimates of height, decreasing R² and increasing RMSE, although the effect was most pronounced for coniferous crowns.     

From satellite altimetry to ocean topography A survey of data processing techniques

Abstract

Satellite altimetry was introduced as a brief experiment on Skylab in 1973. Since this experiment, operational radar altimeters have been flown on GEOS-3 and Seasat and currently the Geosat altimetry misssion is in progress. For the early 1990s additional altimeters are planned for launch on the European ERS-1 and ERS-2 satellites, the U.S.-French Topex/Poseidon satellite, the U.S.NROSS satellite and the Japanese MOS-2 satellite. This paper presents a brief survey of altimetry missions and of some important aspects of altimeter data processing. In addition, the topic of orbit computation accuracy is addressed in somewhat more detail and some examples of geophysical products obtained from Seasat altimeter measurements are discussed. These examples do not necessarily represent products of the highestachievablequality, but they serve to illustrate the possibilities and limitations of present-day processing techniques.     

Long-term ERS/ENVISAT deformation time-series generation at full spatial resolution via the extended SBAS technique

We extend the small baseline subset (SBAS) differential synthetic aperture radar (SAR) interferometry (DInSAR) approach to allow the generation of deformation time-series by processing, at the full spatial resolution scale, long sequences of European Remote Sensing (ERS-1 and ERS-2) and Environmental Satellite (ENVISAT) SAR data acquired with the same illumination geometry. In particular, we avoid the generation of ERS/ENVISAT cross-interferograms, which are severely affected by noise phenomena due to the carrier frequency separation of the two SAR systems, and we focus on single-platform interferograms only (i.e. ERS/ERS and ENVISAT/ENVISAT interferograms) that are properly combined by applying the singular value decomposition (SVD)-based SBAS approach. Moreover, we exploit the Doppler centroid variations of the post-2000 acquisitions of the ERS-2 sensor and the carrier frequency difference between the ERS-1/2 and the ENVISAT systems, in order to maximize the number of investigated SAR pixels and to improve their geocoding. The presented results, achieved on two data sets relevant to the Napoli Bay area and to the Murge region, both located in southern Italy, confirm the effectiveness of the extended SBAS technique and demonstrate the relevance of deformation analysis carried out at the scale of single buildings or human-made structures with more than 15 years of ERS and ENVISAT acquisitions.     

Use of ENVISAT/ASAR wide-swath data for timely rice fields mapping in the Mekong River Delta

Because of the importance of rice for the global food security and because of the role of inundated paddy fields in greenhouse gases emissions, monitoring the rice production world-wide has become a challenging issue for the coming years. Local rice mapping methods have been developed previously in many studies by using the temporal change of the backscatter from C-band synthetic aperture radar (SAR) co-polarized data. The studies indicated in particular the need of a high observation frequency. In the past, the operational use of these methods has been limited by the small coverage and the poor acquisition frequency of the available data (ERS-1/2, Radarsat-1). In this paper, the method is adapted for the first time to map rice at large scale, by using wide-swath images of the Advanced SAR (ASAR) instrument onboard ENVISAT. To increase the observation frequency, data from different satellite tracks are combined. The detection of rice fields is achieved by exploiting the high backscatter increase at the beginning of the growing cycle, which allows the production of rice maps early in the season (in the first 50 days). The method is tested in the Mekong delta in Vietnam. The mapping results are compared to existing rice maps in the An Giang province, with a good agreement (higher than 81%). The rice planted areas are retrieved from the maps and successfully validated with the official statistics available at each province (R² = 0.92). These results show that the method is useful for large scale early mapping of rice areas, using current and future C band wide-swath SAR data.     

Water level dynamics of Amazon wetlands at the watershed scale by satellite altimetry

In this study we used satellite altimetry to characterize the time and space variations in water stored in or circulating through rivers, floodplains, wetlands and lakes in the major sub-basins of the Amazon basin. Using a specific methodology to rigorously select original three-dimensional (3D) data from an Environmental Satellite (ENVISAT) mission, water level time series were calculated at the crossing path of the satellite tracks with the water bodies. We took advantage of the continuous sampling of the water level along the satellite track segments that cross the watershed to analyse both spatial and temporal relationships between: (i) the river and its floodplain and (ii) different basins. This work evidences in particular the existence of water leaking between the Negro and Solimões basins at the high water stage. It highlights that the phenomenon of a secondary flood peak occurring in the water level series in the Solimões basin at rising water, known as repiquete, is caused by the rain equatorial regime of the northern upstream tributaries of the Solimões River, but is disconnected from the same phenomenon occurring within the Rio Negro basin.     

An assessment of chlorophyll-a algorithms available for SeaWiFS in coastal and open areas of the Bay of Bengal and Arabian Sea

Three ocean colour algorithms, OC4v6, Carder and OC5 were tested for retrieving Chlorophyll-a (Chla) in coastal areas of the Bay of Bengal and open ocean areas of the Arabian Sea. Firstly, the algorithms were run using ~ 80 in situ Remote Sensing Reflectance, (R_rs(λ)) data collected from coastal areas during eight cruises from January 2000 to March 2002 and the output was compared to in situ Chla. Secondly, the algorithms were run with ~ 20 SeaWiFS R_rs(λ) and the results were compared with coincident in situ Chla. In both cases, OC5 exhibited the lowest log₁₀-RMS, bias, had a slope close to 1 and this algorithm appears to be the most accurate for both coastal and open ocean areas. Thirdly the error in the algorithms was regressed against Total Suspended Material (TSM) and Coloured Dissolved Organic Material (CDOM) data to assess the co-variance with these parameters. The OC5 error did not co-vary with TSM and CDOM. OC4v6 tended to over-estimate Chla > 2 mg m⁻³ and the error in OC4v6 co-varied with TSM. OC4v6 was more accurate than the Carder algorithm, which over-estimated Chla at concentrations > 1 mg m⁻³ and under-estimated Chla at values < 0.5 mg m⁻³. The error in Carder Chla also co-varied with TSM. The algorithms were inter-compared using > 5500 SeaWiFS R_rs(λ) data from coastal to offshore transects in the Northern Bay of Bengal. There was good agreement between OC4v6 and OC5 in open ocean waters and in coastal areas up to 2 mg m⁻³. There was a strong divergence between Carder and OC5 in open ocean and coastal waters. OC4v6 and Carder tended to over-estimate Chla in coastal areas by a factor of 2 to 3 when TSM > 25 g m⁻³. We strongly recommend the use of OC5 for coastal and open ocean waters of the Bay of Bengal and Arabian Sea. A Chla time series was generated using OC5 from 2000 to 2003, which showed that concentrations at the mouths of the Ganges reach a maxima (~ 5 mg m⁻³) in October and November and were 0.08 mg m⁻³ further offshore increasing to 0.2 mg m⁻³ during December. Similarly in early spring from February to March, Chla was 0.08 to 0.2 mg m⁻³ on the east coast of the Bay.     

Large area forest stem volume mapping in the boreal zone using synergy of ERS-1/2 tandem coherence and MODIS vegetation continuous fields

ERS-1/2 tandem coherence was reported to have high potential for the mapping of boreal forest stem volume (e.g. Santoro et al., 2002, 2007a; Wagner et al., 2003; Askne & Santoro, 2005). Large-scale application of the data for forest stem volume mapping, however, is hindered by the variability of coherence with meteorological and environmental acquisition conditions. The traditional way of stem volume retrieval is based on the training of models, relating coherence to stem volume, with the aid of forest inventory data which is generally available for a few small test sites but not for large areas. In this paper a new approach is presented that allows model training using the MODIS Vegetation Continuous Fields canopy cover product (Hansen et al., 2003) without further need for ground data. A comparison of the new approach with the traditional regression-based and ground-data dependent model training is presented in this paper for a multi-seasonal ERS-1/2 tandem dataset covering several well known Central Siberian forest sites. As a test scenario for large-area application, the approach was applied to a multi-seasonal ERS-1/2 tandem dataset of 223 ERS-1 and ERS-2 image pairs covering Northeast China (~ 1.5 million km²) to map four stem volume classes (0-20, 20-50, 50-80, and > 80 m³/ha).     

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.     

Assessment of the impacts of surface topography, off-nadir pointing and vegetation structure on vegetation lidar waveforms using an extended geometric optical and radiative transfer model

In order to prioritize the measurement requirements and accuracies of the two new lidar missions, a physical model is required for a fundamental understanding of the impact of surface topography, footprint size and off-nadir pointing on vegetation lidar waveforms and vegetation height retrieval. In this study, we extended a well developed Geometric Optical and Radiative Transfer (GORT) vegetation lidar model to take into account for the impacts of surface topography and off-nadir pointing on vegetation lidar waveforms and vegetation height retrieval and applied this extended model to assess the aforementioned impacts on vegetation lidar waveforms and height retrieval.Model simulation shows that surface topography and off-nadir pointing angle stretch waveforms and the stretching effect magnifies with footprint size, slope and off-nadir pointing angle. For an off-nadir pointing laser penetrating vegetation over a slope terrain, the waveform is either stretched or compressed based on the relative angle. The stretching effect also results in a disappearing ground peak return when slope or off-nadir pointing angle is larger than the “critical slope angle”, which is closely related to various vegetation structures and footprint size. Model simulation indicates that waveform shapes are affected by surface topography, off-nadir pointing angle and vegetation structure and it is difficult to remove topography effects from waveform extent based only on the shapes of waveform without knowing any surface topography information.Height error without correction of surface topography and off-nadir pointing angle is the smallest when the laser beams at the toward-slope direction and the largest from the opposite direction. Further simulation reveals within 20° of slope and off-nadir pointing angle, given the canopy height as roughly 25 m and the footprint size as 25 m, the error for vegetation height (RH100) ranges from − 2 m to greater than 12 m, and the error for the height at the medium energy return (RH50) from − 1 m to 4 m. The RH100 error caused by unknown surface topography and without correction of off-nadir pointing effect can be explained by an analytical formula as a function of vegetation height, surface topography, off-nadir pointing angle and footprint size as a first order approximation. RH50 is not much affected by topography, off-nadir pointing and footprint size. This forward model simulation can provide scientific guidance on prioritizing future lidar mission measurement requirements and accuracies.     

Long-term evaluation of three satellite ocean color algorithms for identifying harmful algal blooms (Karenia brevis) along the west coast of Florida: A matchup assessment

We present a simple algorithm to identify Karenia brevis blooms in the Gulf of Mexico along the west coast of Florida in satellite imagery. It is based on an empirical analysis of collocated matchups of satellite and in situ measurements. The results of this Empirical Approach is compared to those of a Bio-optical Technique - taken from the published literature - and the Operational Method currently implemented by the NOAA Harmful Algal Bloom Forecasting System for K. brevis blooms. These three algorithms are evaluated using a multi-year MODIS data set (from July, 2002 to October, 2006) and a long-term in situ database. Matchup pairs, consisting of remotely-sensed ocean color parameters and near-coincident field measurements of K. brevis concentration, are used to assess the accuracy of the algorithms. Fair evaluation of the algorithms was only possible in the central west Florida shelf (i.e. between 25.75°N and 28.25°N) during the boreal Summer and Fall months (i.e. July to December) due to the availability of valid cloud-free matchups. Even though the predictive values of the three algorithms are similar, the statistical measure of success in red tide identification (defined as cell counts in excess of 1.5 × 10⁴ cells L⁻¹) varied considerably (sensitivity—Empirical: 86%; Bio-optical: 77%; Operational: 26%), as did their effectiveness in identifying non-bloom cases (specificity—Empirical: 53%; Bio-optical: 65%; Operational: 84%). As the Operational Method had an elevated frequency of false-negative cases (i.e. presented low accuracy in detecting known red tides), and because of the considerable overlap between the optical characteristics of the red tide and non-bloom population, only the other two algorithms underwent a procedure for further inspecting possible detection improvements. Both optimized versions of the Empirical and Bio-optical algorithms performed similarly, being equally specific and sensitive (~ 70% for both) and showing low levels of uncertainties (i.e. few cases of false-negatives and false-positives: ~ 30%)—improved positive predictive values (~ 60%) were also observed along with good negative predictive values (~ 80%).     

Note on the quality of the (A)ATSR land surface temperature record from 1991 to 2009

Land-surface temperature (LST) is a key parameter in the physical study of atmosphere–land interactions as well as for global warming and climate change monitoring on a longer timescale. Remote sensing provides an excellent way to measure LST at a global scale with appropriate spatial and temporal resolution. To retrieve LSTs, measured radiances at the sensor have to be corrected for surface emissivity, atmospheric effects and contaminating clouds. This study is based on the LST data provided by the Along-Track Scanning Radiometers (ATSR-1 and ATSR-2) and the Advanced ATSR (AATSR) on board the three European Space Agency (ESA) satellites, European Remote-Sensing Satellite 1 (ERS-1), ERS-2 and Enviromental Satellite (ENVISAT). The analysis covers data from August 1991 up to December 2009 and contains detailed investigations on global as well as on regional scales with a temporal resolution of 1 month, outlining problems and restrictions within the time series mainly due to cloud contamination and effects of calibration drifts on the cloud detection tests. It is demonstrated that the trends observed, for cooling as well as for warming, are likely to be related to the trends in cloud contamination rather than the trends in the actual LST.     

Observations of Lake Baikal ice from satellite altimetry and radiometry

We demonstrate the potential of combining satellite altimetry and radiometry for lake ice studies using the example of Lake Baikal in Siberia. We show the synergy using active and passive microwave observations available from the recent satellite altimetry missions (TOPEX/Poseidon, Jason-1, ENVISAT and Geosat Follow-On), complemented by the SSM/I passive data. We assess the applicability of altimetric and radiometric data for ice/water discrimination, and discuss the drawbacks and benefits of each type of sensor. An ice discrimination method, based on the combined use of the data from the four altimetric missions and SSM/I, is proposed and validated using available in situ observations and MODIS imagery. The method is applied to the entire satellite data set and used to define specific dates of ice events (first appearance of ice, formation of stable ice cover, first appearance of open water, complete disappearance of ice) and associated uncertainties. Using these satellite-derived estimates, we can extend the existing time series of ice events in the Southern Baikal up to 2004 and provide new information on the Middle and Northern Baikal, regions where no recent in situ ice cover observations are available. Our data show that over the last 10-15 years, trends towards earlier ice formation and later ice break-up result in a tendency for longer fast ice duration over the whole Lake Baikal. The methods proposed here have been tested for Lake Baikal, but they are applicable for other lakes and water bodies with seasonal ice cover.     

Influence of phytoplankton pigment composition on remote sensing of cyanobacterial biomass

An extensive field campaign was carried out for the validation of a previously published reflectance ratio-based algorithm for quantification of the cyanobacterial pigment phycocyanin (PC). The algorithm uses band settings of the Medium Resolution Imaging Spectrometer (MERIS) onboard ENVISAT, and should accurately retrieve the PC concentration in turbid, cyanobacteria-dominated waters. As algae and cyanobacteria often co-occur, the algorithm response to varying phytoplankton composition was explored. Remote sensing reflectance and reference pigment measurements were obtained in the period 2001-2005 in Spain and the Netherlands using field spectroradiometry and various pigment extraction methods. Additional field data was collected in Spain in May 2005 to allow intercalibration of spectroradiometry and pigment assessment methods. Two methods for extraction of PC from concentrated water samples, and in situ measured PC fluorescence, compared well. Reflectance measurements with different field spectroradiometers used in Spain and the Netherlands also gave similar results. Residual analysis of PC predicted by the algorithm showed that overestimation of PC mainly occurred in the presence of chlorophylls b and c, and phaeophytin. The errors were strongest at low PC relative to Chl a concentrations. A correction applied for absorption by Chl b markedly improved the prediction. Without such a correction, the quality of the PC prediction still increased markedly with estimates > 50 mg PC m^− 3, allowing monitoring of the cyanobacterial status of eutrophic waters. The threshold concentration may be lowered when a high intracellular PC:Chl a ratio or cyanobacterial dominance is expected. Below the limit, predicted PC concentrations should be considered as the highest estimate. We evaluated that remote sensing of both PC and Chl a would allow assessment of cyanobacterial risk to water quality and public health in over 70% of our cases.     

Ice sheet survey over Antarctica using satellite altimetry: ERS-2, Envisat,SARAL/AltiKa,the key importance of continuous observations along the same repeat orbit

From September 2002 to October 2010, the Envisat radar altimeter surveyed Greenland and Antarctica ice sheets on a 35 day repeat orbit, providing a unique data set for ice sheet mass balance studies. Up to 85 repeat cycles are available and the whole Envisat data set may be along-track processed in order to provide height variability and trend with a good spatial resolution for the objectives of ice sheet survey.

Soon, a joint Centre National d’Etudes Spatiales/Indian Space Research Organisation mission, SARAL (Satellite with Argos and AltiKa), with the AltiKa payload on board, will be launched on exactly the same orbit (less than 1 km of the nomimal orbit in the across-track direction). This will allow an extension of previous European Remote Sensing (ERS) satellite, ERS-1 and ERS-2, and Envisat missions of the European Space Agency (ESA), in particular from the point of view of ice altimetry. However, AltiKa operates in the Ka band (36.8 GHz), a higher frequency than the classical Ku band (13.6 GHz), leading to important modifications and potential improvements in the interactions between radar wave and snow-pack.

In this paper, a synthesis is presented of all available information relevant to ice altimetry scientific purposes as derived from the Envisat mission: mean and temporal derivatives of the height ? but also of the backscatter and of the two waveform parameters ? snow-pack change corrections, across-track surface slope at 1 km scale, etc. The spatial and temporal variability of ice sheet surface elevation is investigated with the help of the high-resolution Envisat along-track observations. We show that at least 50 repeat cycles are needed to reach the required accuracy for the elevation trend. It is thus advocated as potentially highly beneficial for SARAL/AltiKa as a follow-on mission. Moreover, the novel characteristics of SARAL/AltiKa are promising in improving our understanding of snow penetration impact.     

Comparison and evaluation of Medium Resolution Imaging Spectrometer leaf area index products across a range of land use

Leaf area index (LAI) is a commonly required parameter when modelling land surface fluxes. Satellite based imagers, such as the 300 m full resolution (FR) Medium Spectral Resolution Imaging Spectrometer (MERIS), offer the potential for timely LAI mapping. The availability of multiple MERIS LAI algorithms prompts the need for an evaluation of their performance, especially over a range of land use conditions. Four current methods for deriving LAI from MERIS FR data were compared to estimates from in-situ measurements over a 3 km × 3 km region near Ottawa, Canada. The LAI of deciduous dominant forest stands and corn, soybean and pasture fields was measured in-situ using digital hemispherical photography and processed using the CANEYE software. MERIS LAI estimates were derived using the MERIS Top of Atmosphere (TOA) algorithm, MERIS Top of Canopy (TOC) algorithm, the Canada Centre for Remote Sensing (CCRS) Empirical algorithm and the University of Toronto (UofT) GLOBCARBON algorithm. Results show that TOA and TOC LAI estimates were nearly identical (R² > 0.98) with underestimation of LAI when it is larger than 4 and overestimation when smaller than 2 over the study region. The UofT and CCRS LAI estimates had root mean square errors over 1.4 units with large (∼ 25%) relative residuals over forests and consistent underestimates over corn fields. Both algorithms were correlated (R² > 0.8) possibly due to their use of the same spectral bands derived vegetation index for retrieving LAI. LAI time series from TOA, TOC and CCRS algorithms showed smooth growth trajectories however similar errors were found when the values were compared with the in-situ LAI. In summary, none of the MERIS LAI algorithms currently meet performance requirements from the Global Climate Observing System.     

Inter-comparison study of water level estimates derived from hydrodynamic-hydrologic model and satellite altimetry for a complex deltaic environment

Riverine deltas are hydrologically one of the most active terrestrial bodies supporting an intricate network of rivers, a highly unsteady flow regime, high agricultural productivity and large population centers. Understanding the complex hydrology of riverine deltas is challenging due to the paucity of conventional ground-based measurements on river water levels and flows that result in large spatial and temporal sampling gaps. One way to bridge this sampling issue is to employ hydrodynamic models in combination with remotely-sensed water level elevation data from satellite altimetry in a data assimilation framework. However, a good understanding of the performance of models and altimetry is required beforehand. Using Bangladesh as an example of a complex delta, an inter-comparison study was therefore performed for water level estimates derived from the two methods: 1) satellite altimetry and 2) hydrodynamic-hydrologic modeling framework. The Envisat mission was selected for satellite altimetry-based water level data. For the modeling framework, a calibrated 1-D hydrodynamic model, HEC-RAS, was set up for the major rivers of Bangladesh using in-situ river bathymetry, gaged stream flow and water level data. Envisat water level estimates were generally found to be exceeded by the model-based values by 0.20 m and 1.90 m for Monsoon and dry seasons, respectively. In general, the average RMSE between Envisat and modeled estimates is more than 2.0 m. The closest agreement with altimetry was observed during the high flow Monsoon season over the Brahmaputra river. Envisat estimates are found to disagree most with model-based estimates for small to medium-sized river basins that are mountainous and flashy. This inter-comparison study provides preliminary guidance on the relative weights to assign for each type of estimate when designing a data assimilation scheme for optimal water level prediction in ungaged basins.