Characterizing a cyanobacterial bloom in western Lake Erie using satellite imagery and meteorological data

The distribution and intensity of a bloom of the toxic cyanobacterium, Microcystis aeruginosa, in western Lake Erie was characterized using a combination of satellite ocean-color imagery, field data, and meteorological observations. The bloom was first identified by satellite on 14 August 2008 and persisted for more than 2 months. The distribution and intensity of the bloom was estimated using a satellite algorithm that is sensitive to near-surface concentrations of M. aeruginosa. Increases in both area and intensity were most pronounced for wind stress less than 0.05 Pa. Area increased while intensity did not change for wind stresses of 0.05–0.1 Pa, and both decreased for wind stress greater than 0.1 Pa. The recovery in intensity at the surface after strong wind events indicated that high wind stress mixed the bloom through the water column and that it returned to the surface once mixing stopped. This interaction is consistent with the understanding of the buoyancy of these blooms. Cloud cover (reduced light) may have a weak influence on intensity during calm conditions. While water temperature remained greater than 15°C, the bloom intensified if there were calm conditions. For water temperature less than 15°C, the bloom subsided under similar conditions. As a result, wind stress needs to be considered when interpreting satellite imagery of these blooms.

Modelling global pyrogeography using data derived from satellite imagery

Doutoramento em Engenharia Florestal e dos Recursos Naturais - Instituto Superior de Agronomia - UL

Integrating field survey data with satellite image data to improve shallow water seagrass maps : the role of AUV and snorkeller surveys

Repeatable and accurate seagrass mapping is required for understanding seagrass ecology and supporting management decisions. For shallow (< 5 m) seagrass habitats, these maps can be created by integrating high spatial resolution imagery with field survey data. Field survey data for seagrass is often collected via snorkelling or diving. However, these methods are limited by environmental and safety considerations. Autonomous Underwater Vehicles (AUVs) are used increasingly to collect field data for habitat mapping, albeit mostly in deeper waters (>20 m). Here we demonstrate and evaluate the use and potential advantages of AUV field data collection for calibration and validation of seagrass habitat mapping of shallow waters (< 5 m), from multispectral satellite imagery. The study was conducted in the seagrass habitats of the Eastern Banks (142 km2), Moreton Bay, Australia. In the field, georeferenced photos of the seagrass were collected along transects via snorkelling or an AUV. Photos from both collection methods were analysed manually for seagrass species composition and then used as calibration and validation data to map seagrass using an established semi-automated object based mapping routine. A comparison of the relative advantages and disadvantages of AUV and snorkeller collected field data sets and their influence on the mapping routine was conducted. AUV data collection was more consistent, repeatable and safer in comparison to snorkeller transects. Inclusion of deeper water AUV data resulted in mapping of a larger extent of seagrass (~7 km2, 5 % of study area) in the deeper waters of the site. Although overall map accuracies did not differ considerably, inclusion of the AUV data from deeper water transects corrected errors in seagrass mapped at depths to 5 m, but where the bottom is visible on satellite imagery. Our results demonstrate that further development of AUV technology is justified for the monitoring of seagrass habitats in ongoing management programs.

Estimating crop area using seasonal time series of enhanced vegetation index from MODIS satellite imagery

Cereal grain is one of the main export commodities of Australian agriculture. Over the past decade, crop yield forecasts for wheat and sorghum have shown appreciable utility for industry planning at shire, state, and national scales. There is now an increasing drive from industry for more accurate and cost-effective crop production forecasts. In order to generate production estimates, accurate crop area estimates are needed by the end of the cropping season. Multivariate methods for analysing remotely sensed Enhanced Vegetation Index (EVI) from 16-day Moderate Resolution Imaging Spectroradiometer (MODIS) satellite imagery within the cropping period (i.e. April-November) were investigated to estimate crop area for wheat, barley, chickpea, and total winter cropped area for a case study region in NE Australia. Each pixel classification method was trained on ground truth data collected from the study region. Three approaches to pixel classification were examined: (i) cluster analysis of trajectories of EVI values from consecutive multi-date imagery during the crop growth period; (ii) harmonic analysis of the time series (HANTS) of the EVI values; and (iii) principal component analysis (PCA) of the time series of EVI values. Images classified using these three approaches were compared with each other, and with a classification based on the single MODIS image taken at peak EVI. Imagery for the 2003 and 2004 seasons was used to assess the ability of the methods to determine wheat, barley, chickpea, and total cropped area estimates. The accuracy at pixel scale was determined by the percent correct classification metric by contrasting all pixel scale samples with independent pixel observations. At a shire level, aggregated total crop area estimates were compared with surveyed estimates. All multi-temporal methods showed significant overall capability to estimate total winter crop area. There was high accuracy at pixel scale (>98% correct classification) for identifying overall winter cropping. However, discrimination among crops was less accurate. Although the use of single-date EVI data produced high accuracy for estimates of wheat area at shire scale, the result contradicted the poor pixel-scale accuracy associated with this approach, due to fortuitous compensating errors. Further studies are needed to extrapolate the multi-temporal approaches to other geographical areas and to improve the lead time for deriving cropped-area estimates before harvest.

High resolution satellite imagery and GPS collars can help assist in the assessment of patch selection by grazing cattle in semi-arid savannas.

The selection of different patch types for grazing by cattle in tropical savannas is well documented. Advances in high resolution satellite imagery and computing power now allow us to identify patch types over an entire paddock, combined with GPS collars as a non instrusive method of capturing positional data, an accurate and comprehensive picture of landscape use by cattle can be quantified.

Patch selection by cattle can be quantified using satellite imagery and GPS in extensive, semi-arid savannas.

Patch selection by grazing animals is difficult to quantify, particularly in large, extensive paddocks like those in northern Australia. However, advances in high resolution satellite imagery now allow identification of patch types over an entire paddock which combined with GPS collars to capture positional data, can give an accurate and comprehensive picture of landscape use by cattle.

Burnt area mapping in insular Southeast Asia using medium resolution satellite imagery

Burnt area mapping in humid tropical insular Southeast Asia using medium resolution (250-500m) satellite imagery is characterized by persisting cloud cover, wide range of land cover types, vast amount of wetland areas and highly varying fire regimes. The objective of this study was to deepen understanding of three major aspects affecting the implementation and limits of medium resolution burnt area mapping in insular Southeast Asia: 1) fire-induced spectral changes, 2) most suitable multitemporal compositing methods and 3) burn scars patterns and size distribution. The results revealed a high variation in fire-induced spectral changes depending on the pre-fire greenness of burnt area. It was concluded that this variation needs to be taken into account in change detection based burnt area mapping algorithms in order to maximize the potential of medium resolution satellite data. Minimum near infrared (MODIS band 2, 0.86μm) compositing method was found to be the most suitable for burnt area mapping purposes using Moderate Resolution Imaging Spectroradiometer (MODIS) data. In general, medium resolution burnt area mapping was found to be usable in the wetlands of insular Southeast Asia, whereas in other areas the usability was seriously jeopardized by the small size of burn scars. The suitability of medium resolution data for burnt area mapping in wetlands is important since recently Southeast Asian wetlands have become a major point of interest in many fields of science due to yearly occurring wild fires that not only degrade these unique ecosystems but also create regional haze problem and release globally significant amounts of carbon into the atmosphere due to burning peat. Finally, super-resolution MODIS images were tested but the test failed to improve the detection of small scars. Therefore, super-resolution technique was not considered to be applicable to regional level burnt area mapping in insular Southeast Asia.

Enhanced processing of SPOT multispectral satellite imagery for environmental monitoring and modelling

The Taita Hills in southeastern Kenya form the northernmost part of Africa’s Eastern Arc Mountains, which have been identified by Conservation International as one of the top ten biodiversity hotspots on Earth. As with many areas of the developing world, over recent decades the Taita Hills have experienced significant population growth leading to associated major changes in land use and land cover (LULC), as well as escalating land degradation, particularly soil erosion. Multi-temporal medium resolution multispectral optical satellite data, such as imagery from the SPOT HRV, HRVIR, and HRG sensors, provides a valuable source of information for environmental monitoring and modelling at a landscape level at local and regional scales. However, utilization of multi-temporal SPOT data in quantitative remote sensing studies requires the removal of atmospheric effects and the derivation of surface reflectance factor. Furthermore, for areas of rugged terrain, such as the Taita Hills, topographic correction is necessary to derive comparable reflectance throughout a SPOT scene. Reliable monitoring of LULC change over time and modelling of land degradation and human population distribution and abundance are of crucial importance to sustainable development, natural resource management, biodiversity conservation, and understanding and mitigating climate change and its impacts. The main purpose of this thesis was to develop and validate enhanced processing of SPOT satellite imagery for use in environmental monitoring and modelling at a landscape level, in regions of the developing world with limited ancillary data availability. The Taita Hills formed the application study site, whilst the Helsinki metropolitan region was used as a control site for validation and assessment of the applied atmospheric correction techniques, where multiangular reflectance field measurements were taken and where horizontal visibility meteorological data concurrent with image acquisition were available. The proposed historical empirical line method (HELM) for absolute atmospheric correction was found to be the only applied technique that could derive surface reflectance factor within an RMSE of < 0.02 ps in the SPOT visible and near-infrared bands; an accuracy level identified as a benchmark for successful atmospheric correction. A multi-scale segmentation/object relationship modelling (MSS/ORM) approach was applied to map LULC in the Taita Hills from the multi-temporal SPOT imagery. This object-based procedure was shown to derive significant improvements over a uni-scale maximum-likelihood technique. The derived LULC data was used in combination with low cost GIS geospatial layers describing elevation, rainfall and soil type, to model degradation in the Taita Hills in the form of potential soil loss, utilizing the simple universal soil loss equation (USLE). Furthermore, human population distribution and abundance were modelled with satisfactory results using only SPOT and GIS derived data and non-Gaussian predictive modelling techniques. The SPOT derived LULC data was found to be unnecessary as a predictor because the first and second order image texture measurements had greater power to explain variation in dwelling unit occurrence and abundance. The ability of the procedures to be implemented locally in the developing world using low-cost or freely available data and software was considered. The techniques discussed in this thesis are considered equally applicable to other medium- and high-resolution optical satellite imagery, as well the utilized SPOT data.

Assessment of insect occurrence in boreal forests based on satellite imagery and field measurements.

The presence/absence data of twenty-seven forest insect taxa (e.g. Retinia resinella, Formica spp., Pissodes spp., several scolytids) and recorded environmental variation were used to investigate the applicability of modelling insect occurrence based on satellite imagery. The sampling was based on 1800 sample plots (25 m by 25 m) placed along the sides of 30 equilateral triangles (side 1 km) in a fragmented forest area (approximately 100 km2) in Evo, S Finland. The triangles were overlaid on land use maps interpreted from satellite images (Landsat TM 30 m multispectral scanner imagery 1991) and digitized geological maps. Insect occurrence was explained using either environmental variables measured in the field or those interpreted from the land use and geological maps. The fit of logistic regression models varied between species, possibly because some species may be associated with the characteristics of single trees while other species with stand characteristics. The occurrence of certain insect species at least, especially those associated with Scots pine, could be relatively accurately assessed indirectly on the basis of satellite imagery and geological maps. Models based on both remotely sensed and geological data better predicted the distribution of forest insects except in the case of Xylechinus pilosus, Dryocoetes sp. and Trypodendron lineatum, where the differences were relatively small in favour of the models based on field measurements. The number of species was related to habitat compartment size and distance from the habitat edge calculated from the land use maps, but logistic regressions suggested that other environmental variables in general masked the effect of these variables in species occurrence at the present scale.

Monitoring changes in phytoplankton abundance and composition in the Northwest Atlantic: a comparison of results obtained by Continuous Plankton Recorder sampling and colour satellite imagery

Phytoplankton abundance in the NW Atlantic was measured by continuous plankton recorder (CPR) sampling along tracks between Iceland and the western Scotian Shelf from 1998 to 2006, when sea-surface chlorophyll (SSChl) measurements were also being made by ocean colour satellite imagery using the SeaWiFS sensor. Seasonal and inter-annual changes in phytoplankton abundance were examined using data collected by both techniques, averaged over each of four shelf regions and four deep ocean regions. CPR sampling had gaps (missing months) in all regions and in the four deep ocean regions satellite observations were too sparse between November and February to be of use. Average seasonal cycles of SSChl were similar to those of total diatom abundance in seven regions, to those of the phytoplankton colour index in six regions, but were not similar to those of total dinoflagellate abundance anywhere. Large inter-annual changes in spring bloom dynamics were captured by both samplers in shelf regions. Changes in annual (or 8 months) averages of SSChl did not generally follow those of the CPR indices within regions and multi-year averages of SSChl, and the three CPR indices were generally higher in shelf than in deep ocean regions. Remote sensing and CPR sampling provide complementary ways of monitoring phytoplankton in the ocean: the former has superior temporal and spatial coverage and temporal resolution, and the latter provides better taxonomic information.

Satellite imagery as a single source of predictor variables for habitat suitability modelling: how Landsat can inform the conservation of a critically endangered lemur

1. Statistical modelling of habitat suitability is an important tool for planning conservation interventions, particularly for areas where species distribution data are expensive or hard to collect. Sometimes however the predictor variables typically used in habitat suitability modelling are themselves difficult to obtain or not meaningful at the geographical extent of the study, as is the case for the Alaotran gentle lemur Hapalemur alaotrensis, a critically endangered lemur confined to the marshes of Lake Alaotra in Madagascar.2. We developed a habitat suitability model where all predictor variables, including vegetation indices and image texture measures at different scales (as surrogates for habitat structure), were derived from Landsat7 satellite imagery. Using relatively few presence records, the maximum entropy (Maxent) approach and AUC were used to assess the performance of candidate predictor variables, for studying the effect of scale, model selection and mapping suitable habitat.3. This study demonstrated the utility of satellite imagery as a single source of predictor variables for a Maxent habitat suitability model at the landscape level, within a restricted geographical extent and with a fine grain, in a case where predictor variables typically used at the macro-scale level (e.g. climatic and topographic) were not applicable.4. In the case of H. alaotrensis, the methodology generated a habitat suitability map to inform conservation management in Lake Alaotra and a replicable protocol to allow rapid updates to habitat suitability maps in the future. The exploration of candidate predictor variables allowed the identification of scales that appear ecologically relevant for the species.5. Synthesis and applications. This study presents a cost-effective combination of maximum entropy habitat suitability modelling and satellite imagery, where all predictor variables are derived solely from Landsat7 images. With a habitat modelling method like Maxent that shows good performance with few presence samples and Landsat images now freely available, the methodology can play an important role in rapid assessments of the status of species at the landscape level in data-poor regions, when typical macro-scale environmental predictors are of little use or difficult to obtain.

Temporal analysis of the reduction in gas emission in areas of mechanically-harvested sugarcane using satellite imagery

The primary objective of this study was to estimate the amount of gas not emitted into the air in areas cultivated with sugarcane (Saccharum officinarum) that were mechanically harvested. Satellite images CBERS-2/CCD, from 08-13-2004, 08-14-2005, 08-15-2006 and 08-16-2007, of northwestern São Paulo State were processed using the Geographic Information System (GIS)-IDRISI 15.0. Areas of interest (the mechanically-harvested sugarcane fields) were identified and quantified based on the spectral response of the bands studied. Based on these data, the amount of gas that was not emitted was evaluated, according to the estimate equation proposed by the Intergovernmental Panel on Climate Change (IPCC). The results of 396.65 km(2) (5.91% for 2004); 447.56 km(2) (6.67% for 2005); 511.54 km(2) (7.62% in 2006); and 474.60 km(2) (7.07% for 2007), calculated from a total area of 6,710.89 km(2) with sugarcane, showed a significant increase of mechanical harvesting in the study area and a reduction of gas emissions of more than 300,000 t yr(-1).

Multitemporal analysis of satellite imagery: implementation of a land evolution model in ALCS (ATSR LAND CLASSIFICATION SYSTEM)

In this report it was designed an innovative satellite-based monitoring approach applied on the Iraqi Marshlands to survey the extent and distribution of marshland re-flooding and assess the development of wetland vegetation cover. The study, conducted in collaboration with MEEO Srl , makes use of images collected from the sensor (A)ATSR onboard ESA ENVISAT Satellite to collect data at multi-temporal scales and an analysis was adopted to observe the evolution of marshland re-flooding. The methodology uses a multi-temporal pixel-based approach based on classification maps produced by the classification tool SOIL MAPPER ®. The catalogue of the classification maps is available as web service through the Service Support Environment Portal (SSE, supported by ESA). The inundation of the Iraqi marshlands, which has been continuous since April 2003, is characterized by a high degree of variability, ad-hoc interventions and uncertainty. Given the security constraints and vastness of the Iraqi marshlands, as well as cost-effectiveness considerations, satellite remote sensing was the only viable tool to observe the changes taking place on a continuous basis. The proposed system (ALCS – AATSR LAND CLASSIFICATION SYSTEM) avoids the direct use of the (A)ATSR images and foresees the application of LULCC evolution models directly to „stock‟ of classified maps. This approach is made possible by the availability of a 13 year classified image database, conceived and implemented in the CARD project (http://earth.esa.int/rtd/Projects/#CARD).The approach here presented evolves toward an innovative, efficient and fast method to exploit the potentiality of multi-temporal LULCC analysis of (A)ATSR images. The two main objectives of this work are both linked to a sort of assessment: the first is to assessing the ability of modeling with the web-application ALCS using image-based AATSR classified with SOIL MAPPER ® and the second is to evaluate the magnitude, the character and the extension of wetland rehabilitation.

A new information fusion method for land-use classification using high resolution satellite imagery

Satellite image classification involves designing and developing efficient image classifiers. With satellite image data and image analysis methods multiplying rapidly, selecting the right mix of data sources and data analysis approaches has become critical to the generation of quality land-use maps. In this study, a new postprocessing information fusion algorithm for the extraction and representation of land-use information based on high-resolution satellite imagery is presented. This approach can produce land-use maps with sharp interregional boundaries and homogeneous regions. The proposed approach is conducted in five steps. First, a GIS layer - ATKIS data - was used to generate two coarse homogeneous regions, i.e. urban and rural areas. Second, a thematic (class) map was generated by use of a hybrid spectral classifier combining Gaussian Maximum Likelihood algorithm (GML) and ISODATA classifier. Third, a probabilistic relaxation algorithm was performed on the thematic map, resulting in a smoothed thematic map. Fourth, edge detection and edge thinning techniques were used to generate a contour map with pixel-width interclass boundaries. Fifth, the contour map was superimposed on the thematic map by use of a region-growing algorithm with the contour map and the smoothed thematic map as two constraints. For the operation of the proposed method, a software package is developed using programming language C. This software package comprises the GML algorithm, a probabilistic relaxation algorithm, TBL edge detector, an edge thresholding algorithm, a fast parallel thinning algorithm, and a region-growing information fusion algorithm. The county of Landau of the State Rheinland-Pfalz, Germany was selected as a test site. The high-resolution IRS-1C imagery was used as the principal input data.