Mapping sediment-laden sea ice in the Arctic using AVHRR remote-sensing data: Atmospheric correction and determination of reflectances as a function of ice type and sediment load

Exploiting the fact that the spectral characteristics of light backscattered from sediment-laden ice differ substantially from those of clean ice and that sediment tends to accumulate at the ice surface during the first melt season, remote-sensing techniques provide a valuable tool for mapping the extent of particle-laden ice in the Arctic basin and assessing its particulate loading. This study considers two fundamental problems that still need to be addressed in order to make full use of satellite observations for this type of assessment: (i) the effects of the atmosphere on surface reflectances derived from radiances measured by the satellite sensor need to be quantified and ultimately corrected for, and (ii) the spectral reflectance of the ice surface as a function of particle loading and sub-pixel distribution needs to be determined in order to derive quantitative estimates from the at-sensor satellite signal. Here, spectral albedos have been computed for different ice surfaces of variable sediment load with a radiative transfer model for sea ice coupled with an optical model for particulates included in sea ice. In a second step, the role of the atmosphere in modulating the surface reflectance signal is assessed with the aid of an atmospheric radiative transfer model applied to a “standard” Arctic atmosphere and surface boundary conditions as prescribed by the sea ice radiative transfer model. A series of sensitivity studies helps assess differences between top-of-the-atmosphere and true surface reflectance and has been utilized to derive a look-up table for atmospheric correction of Advanced Very High Resolution Radiometer (AVHRR) data over sediment-laden sea ice surfaces. In particular, the effects of solar elevation, viewing geometry, and atmospheric properties are considered. The atmospheric corrections are necessary for certain geometries and surface types. Large discrepancies between raw and corrected data are particularly evident in the derived coverage of clean ice and ice with small sediment loading.     

Submarine groundwater discharge in tidal flats revealed by space-borne synthetic aperture radar

Most space-borne sensors cannot detect subsurface features. Groundwater is a typical subsurface feature, and its discharge to coastal ocean waters plays an important role in transporting terrestrial chemical constituents and providing habitats for various species of fauna and flora. This is the first paper to report observational evidence for submarine groundwater discharge (SGD) in tidal flats using space-borne synthetic aperture radar (SAR). Tidal flats are composed of high-moisture-saturated sediments and water puddles. These shallow water puddles were imaged effectively by using SAR systems. The presence of water puddles is usually indicated by low radar backscatter in SAR images due to specular reflections on the water surface. This effect was proved by comparing radar backscattering coefficients obtained from two space-borne SAR systems, TerraSAR-X and RADARSAT-2, with those obtained from two theoretical scattering models, IEM and Oh model. We observed relatively large, widely distributed water puddles in belt shape along the upper parts of the tidal flat, which were confirmed to be related to the discharge of groundwater. The results of this research suggest that SAR can be a powerful tool for observing and determining the areal distributions of possible groundwater discharge in large tidal flats, which is normally difficult to detect with traditional measurement tools or survey techniques for groundwater discharge. We firmly believe that this technique can reduce significantly the efforts of field work to confirm SGD in tidal flats.     

Analysis of climate change impacts on lake ice phenology in Canada using the historical satellite data record

Variability and trends in lake ice dynamics (i.e. lake ice phenology) are related to climate conditions. Climate influences the timing of lake ice melt and freeze onset, ice duration, and lake thermal dynamics that feedback to the climate system initiating further change. Phenology records acquired in a consistent manner and over long time periods are required to better understand variability and change in climate conditions and how changes impact lake processes. In this study, we present a new technique for extracting lake ice phenology events from historical satellite records acquired by the series of Advanced Very High Resolution Radiometer (AVHRR) sensors. The technique was used to extend existing in-situ measurements for 36 Canadian lakes and to develop records for 6 lakes in Canada's far north. Comparison of phenology events obtained from the AVHRR record and in-situ measurements show strong agreement (20 lakes, 180 cases) suggesting, with high confidence especially in the case of break-up dates, the use of these data as a complement to ground observations. Trend analysis performed using the combined in-situ and AVHRR record ∼ 1950-2004 shows earlier break-up (average — 0.18 days/year) and later freeze-up (average 0.12 days/year) for the majority of lakes analyzed. Less confidence is given to freeze-up date results due to lower sun elevation during this period making extraction more difficult. Trends for the 20 year record in the far north showed earlier break-up (average 0.99 days/year) and later freeze-up (average 0.76 days/year). The established lake ice phenology database from the historical AVHRR image archive for the period from 1985 to 2004 will to a certain degree fill data gaps in the Canadian in-situ observation network. Furthermore, the presented extraction procedure is not sensor specific and will enable continual data update using all available satellite data provided from sensors such as NOAA/AVHRR, MetOp/AVHRR, MODIS, MERIS and SPOT/VGT.     

Preserving boundaries for image texture segmentation using grey level co-occurring probabilities

Texture analysis has been used extensively in the computer-assisted interpretation of digital imagery. A popular texture feature extraction approach is the grey level co-occurrence probability (GLCP) method. Most investigations consider the use of the GLCP texture features for classification purposes only, and do not address segmentation performance. Specifically, for segmentation, the pixels in an image located near texture boundaries have a tendency to be misclassified. Boundary preservation when using the GLCP texture features for image segmentation is important. An advancement which exploits spatial relationships has been implemented. The generated features are referred to as weighted GLCP (WGLCP) texture features. In addition, an investigation for selecting suitable GLCP parameters for improved boundary preservation is presented. From the tests, WGLCP features provide improved boundary preservation and segmentation accuracy at a computational cost. As well, the GLCP correlation statistical parameter should not be used when segmenting images with high contrast texture boundaries.     

Integration of radar and Landsat-derived foliage projected cover for woody regrowth mapping, Queensland, Australia

In Queensland, Australia, forest areas are discriminated from non-forest by applying a threshold (∼ 12%) to Landsat-derived Foliage Projected Cover (FPC) layers (equating to ∼ 20% canopy cover), which are produced routinely for the State. However, separation of woody regrowth following agricultural clearing cannot be undertaken with confidence, and is therefore not mapped routinely by State Agencies. Using fully polarimetric C-, L- and P-band NASA AIRSAR and Landsat FPC data for forests and agricultural land near Injune, central Queensland, we corroborate that woody regrowth dominated by Brigalow (Acacia harpophylla) cannot be discriminated using either FPC or indeed C-band data alone, because the rapid attainment of a canopy cover leads to similarities in both reflectance and backscatter with remnant forest. We also show that regrowth cannot be discriminated from non-forest areas using either L-band or P-band data alone. However, mapping can be achieved by thresholding and intersecting these layers, as regrowth is unique in supporting both a high FPC (> ∼ 12%) and C-band SAR backscatter (> ~ − 18 dB at HV polarisation) and low L-band and P-band SAR backscatter (e.g. < =∼ 14 dB at L-band HH polarisation). To provide a theoretical explanation, a wave scattering model based on that of Durden et al. [Durden, S.L., Van Zyl, J.J. & Zebker, H.A. (1989). Modelling and observation of radar polarization signature of forested areas. IEEE Trans. Geoscience and Remote Sensing, 27, 290-301.] was used to demonstrate that volume scattering from leaves and small branches in the upper canopy leads to increases in C-band backscattering (particularly HV polarisations) from regrowth, which increases proportionally with FPC. By contrast, low L-band and P-band backscatter occurs because of the lack of double bounce interactions at co-polarisations (particularly HH) and volume scattering at HV polarisation from the stems and branches, respectively, when their dimensions are smaller than the wavelength. Regrowth maps generated by applying simple thresholds to both FPC and AIRSAR L-band data showed a very close correspondence with those mapped using same-date 2.5 m Hymap data and an average 73.7% overlap with those mapped through time-series comparison of Landsat-derived land cover classifications. Regrowth mapped using Landsat-derived FPC from 1995 and JER-1 SAR data from 1994-1995 also corresponded with areas identified within the time-series classification and true colour stereo photographs for the same period. The integration of Landsat FPC and L-band SAR data is therefore expected to facilitate regrowth mapping across Queensland and other regions of Australia, particularly as Japan's Advanced Land Observing System (ALOS) Phase Arrayed L-band SAR (PALSAR), to be launched in 2006, will observe at both L-band HH and HV polarisations.     

Classification of sediments on exposed tidal flats in the German Bight using multi-frequency radar data

We present a new method for the extraction of roughness parameters of sand ripples on exposed tidal flats from multi-frequency synthetic aperture radar (SAR) data. The method is based on the Integral Equation Model (IEM) which predicts the normalized radar cross-section (NRCS) of randomly rough dielectric surfaces. The data used for this analysis were acquired in the German Bight of the North Sea by the Spaceborne Imaging Radar-C/X-Band SAR (SIR-C/X-SAR) in 1994. In-situ measurements of the root-mean-squared (rms) height and the correlation length of the sand ripples clearly demonstrate a relationship between these roughness parameters and the C-band NRCS determined from an ERS SAR image. Using the IEM we have calculated NRCS isolines for the three frequency bands deployed by SIR-C/X-SAR (L, C, and X band), as a function of the rms height and the correlation length of the sand ripples. For each SIR-C/X-SAR image pixel these two roughness parameters were determined from the intersections of the NRCS isolines at different radar bands, and they were used for a crude sediment classification for a small test area at the German North Sea coast. Comparing our results with available sediment maps, we conclude that the presented method is very promising for tidal flat classification by using data from presently existing airborne and future spaceborne multi-frequency SAR systems.     

Multi-temporal monitoring of wetland water levels in the Florida Everglades using interferometric synthetic aperture radar (InSAR)

Interferometric synthetic aperture radar (InSAR) techniques can successfully detect phase variations related to the water level changes in wetlands and produce spatially detailed high-resolution maps of water level changes. Despite the vast details, the usefulness of the wetland InSAR observations is rather limited, because hydrologists and water resources managers need information on absolute water level values and not on relative water level changes. We present an InSAR technique called Small Temporal Baseline Subset (STBAS) for monitoring absolute water level time series using radar interferograms acquired successively over wetlands. The method uses stage (water level) observation for calibrating the relative InSAR observations and tying them to the stage's vertical datum. We tested the STBAS technique with two-year long Radarsat-1 data acquired during 2006-2008 over the Water Conservation Area 1 (WCA1) in the Everglades wetlands, south Florida (USA). The InSAR-derived water level data were calibrated using 13 stage stations located in the study area to generate 28 successive high spatial resolution maps (50 m pixel resolution) of absolute water levels. We evaluate the quality of the STBAS technique using a root mean square error (RMSE) criterion of the difference between InSAR observations and stage measurements. The average RMSE is 6.6 cm, which provides an uncertainty estimation of the STBAS technique to monitor absolute water levels. About half of the uncertainties are attributed to the accuracy of the InSAR technique to detect relative water levels. The other half reflects uncertainties derived from tying the relative levels to the stage stations' datum.     

Unsupervised image segmentation using a simple MRF model with a new implementation scheme

A simple Markov random field model with a new implementation scheme is proposed for unsupervised image segmentation based on image features. The traditional two-component MRF model for segmentation requires training data to estimate necessary model parameters and is thus unsuitable for unsupervised segmentation. The new implementation scheme solves this problem by introducing a function-based weighting parameter between the two components. Using this method, the simple MRF model is able to automatically estimate model parameters and produce accurate unsupervised segmentation results. Experiments demonstrate that the proposed algorithm is able to segment various types of images (gray scale, color, texture) and achieves an improvement over the traditional method.