Tropical Deforestation and Habitat Fragmentation in the Amazon: Satellite Data from 1978 to 1988

David Skole and Compton Tucker

Landsat satellite imagery covering the entire forested portion of the Brazilian Amazon Basin was used to measure, for 1978 and 1988, deforestation, fragmented forest, defined as areas less than 100 square kilometers surrounded by deforestation, and edge effects of 1 kilometer into forest from adjacent areas of deforestation. Tropical deforestation increased from 78,000 square kilometers in 1978 to 230,000 square kilometers in 1988 while tropical forest habitat, severely affected with respect to biological diversity, increased from 208,000 to 588,000 square kilometers. Although this rate of deforestation is lower than previous estimates, the effect on biological diversity is greater.


D. Skole is with the Institute for the Study of Earth, Oceans, and Space, University of New Hampshire, Durham, NH 03824. C. Tucker is with the Laboratory for Terrestrial Physics, NASA Goddard Space Flight Center, Greenbelt, MD 20771.


Deforestation has been occurring in temperate and tropical regions throughout history (1). In recent years, much attention has focused on tropical forests, where as much as 50% of the original extent may have been lost to deforestation in the last two decades, primarily as a result of agricultural expansion (2). Global estimates of tropical deforestation range from 69,000 km2 year-1 in 1980 (3) to 100,000 to 165,000 km2 year-1 in the late 1980s; 50 to 70% of the more recent estimates have been attributed to deforestation in the Brazilian Amazon, the largest continuous region of tropical forest in the world (2, 4, 5).

The area and rate of deforestation in Amazonia are not well known, nor are there quantitative measurements of the effect of deforestation on habitat degradation. We used 1:500,000 scale photographic imagery from Landsat Thematic Mapper data and a geographic information system (GIS) to create a computerized map of deforestation and evaluate its influence on forest fragmentation and habitat degradation. Areas of deforestation were digitized into the GIS and the forest fragments and edge effects that result from the spatial pattern of forest conversion were determined.


Fig.1.Landsat Thematic Mapper color composite image of southern Rondonia state, Brazil, for path 230 and row 69 acquired on 5 June 1988. Areas of tropical forest, deforestation, regrowth, and isolated forest are labeled. The area identified as isolated forest is about 3 km by 15 km in size.


Background

Tropical deforestation is a major component of the carbon cycle and has profound implications for biological diversity. Deforestation increases atmospheric CO2 and other trace gases, possibly affecting climate (6,7). Conversion of forests to cropland and pasture results in a net flux of carbon to the atmosphere because the concentration of carbon in forests is higher than that in the agricultural areas that replace them. The paucity of data on tropical deforestation limits our understanding of the carbon cycle and possible climate change (8). Furthermore, while occupying less than 7% of the terrestrial surface, tropical forests are the home to half or more of all plant and animal species (9). The primary adverse effect of tropical deforestation is massive extinction of species including, for the first time, large numbers of vascular plant species (10).

Deforestation affects biological diversity in three ways: destruction of habitat, isolation of fragments of formerly contiguous habitat, and edge effects within a boundary zone between forest and deforested areas. This boundary zone extends some distance into the remaining forest. In this zone there are greater exposure to winds; dramatic micrometeorological differences over short distances; easier access for livestock, other nonforest animals, and hunters; and a range of other biological and physical effects. The result is a net loss of plant and animal species in the edge areas (11).

There is a wide range in current estimates of the area and rate of deforestation in Amazonia. Scientists at the Instituto Nacional de Pesquisas Espaciais (12-15) estimated a total deforested area of 280,000 km2 as of 1988 and an average annual rate of 21,000 km2 year-1 from 1978 to 1988. Other studies (2, 4, 5) have reported rates that range from 50,000 to 80,000 km2 year-1 (Table 1). Additional deforestation estimates have been made for geographically limited study areas in the southern Amazon Basin of Brazil with Landsat and meteorological satellite data (16-20).


Table 1.Tropical forest area (3) and reported tropical deforestation rates by country. The deforestation rates from the 1970s are from the Food and Agriculture Organization (FAO) (3). The 1980s data are from Meyers (2) and the World Resources Institute (WRI) (5).


The Amazon Basin of Brazil has been defined by law to include the states of Acre, Amapa, Amazonas, Para, Rondonia, Roraima plus part of Mato Grosso, Maranh o, and Tocantins and is referred to as the Legal Amazon (21). It covers an area of ~5,000,000 km2, of which ~4,090,000 km2 is forested, ~850,000 km2 is cerrado or tropical savanna, and ~90,000 km2 is water (Table 2). Confusion has arisen among researchers regarding the stratification of the Brazilian Amazon into forest, cerrado, and water strata. A Food and Agriculture Organization (FAO)-United Nations Environmental Program (UNEP) study (3) found 3,562,800 km2 of forest, whereas Fearnside and co-workers claim there is 4,195,660 km2 of forest, 793,279 km2 of cerrado (17), and 4,906,784 km2 total (13). Meanwhile, an IBGE study (22) found 20,972 km2 of water, 3,793,664 km2 of forest, and 1,149,943 km2 of cerrado for a total of 4,964,920 km2. These differences prevent comparison of different deforestation studies.


Table 2.Predeforestation water, forest, and cerrado land cover for the Brazilian Amazon by state as used in this study. The values determined in this study were based on the IBGE vegetation map and by interpretation of satellite data (22). Areas obscured by clouds were excluded from deforestation and affected forest habitat analyses (97% of the cloud-affected data were over tropical forest).


The use of satellite data and the GIS make it possible to explicitly stratify Amazonia on the basis of cover types (22), thereby providing a means of comparison with other studies. This approach is also necessary for spatial analysis of habitat fragmentation and edge effects of deforestation. Finally, GIS provides a data management tool with which we could manage large amounts of spatial data and precisely merge and geocode information from the more than 200 satellite images used in this study.

Remote Sensing

The large area of the Brazilian Amazon necessitates a straightforward and accurate method of measurement. Landsat Thematic Mapper photo products are inexpensive and of sufficient spatial and spectral resolution for the determination of deforestation. Analysis with visual interpretation techniques produces quantitative results similar to digital processing of full- resolution, multispectral data from the Thematic Mapper and SPOT (23).

We acquired 210 black and white photographic images of the entire Brazilian Amazon. They were obtained with channel five of the Landsat Thematic Mapper (1.55 to 1.75 micrometers) at 1:500,000 scale and were primarily from 1988 (24). We digitized the deforested areas with visual deforestation interpretation and standard vector GIS techniques (Fig. 1). The digitized scenes were projected into equal-area geographic coordinates (latitude, longitude), edge matched, and merged in the computer to form a single, seamless data set for the entire Brazilian Amazon.


Fig.2.Representation of deforestation in the Amazon of Brazil from (A) 1978 and (B) 1988. THe deforestation represented in these figures is confined exclusively to the forest strata. The data were averaged into 16 km grid cells. Spatial analysis of the geometry of deforestation is critical to the estimation of forest fragmentation and the edge effect. If 100 km2 of tropical deforestation occurs as a 10 km by 10 km square and we assume that the edge effect is 1 km, the total area affected is ~143 km2. In contrast, if the 100 km2 of deforestation is distributed as ten strips, each 10 km by 1 km, the affected area is ~350 km2.


We extracted forest fragments <100 km2 that were isolated by deforestation and computed edge effects for a zone of 1 km along the boundaries. All areas of closed-canopy tropical forest deforested by 1988 were delineated, including areas of secondary growth on abandoned fields and pastures where visible (Fig. 1). Areas of long-term forest degradation along river margins in central Amazonia were also included, as were scattered small clearings associated with rubber tappers, mining operations, airfields, and other small disturbances. All visible roads, power line right of ways, pipelines, and similar human-made features were also digitized into the GIS and treated as deforestation. We used 50 digital Landsat Multispectral Scanner (MSS) scenes from 1986 and 15 digital Thematic Mapper images from 1988 for detailed examination of Acre, Amazonas, Mato Grosso, Para, and Rondonia.

To determine the extent of deforestation in 1978, we used the GIS to digitize maps of scale 1:500,000 from single-channel Landsat MSS data, produced jointly by the Instituto Brasiliero de Desenvolvimento Florestal (IBDF) and the Instituto de Pesquisas Espaciais (INPE) in the early 1980s (12, 23). These maps did not differentiate between forest and cerrado clearing. We compiled forest, cerrado, and water data by combining a vegetation map with analysis of Landsat images and meteorological satellite data (25). Our deforestation and affected habitat analyses for 1978 and 1988 were restricted to closed-canopy forest of the Brazilian Amazon.


Fig.3.Map of the Brazilian Amazon Basin showing where biological diversity was adversely affected in 1988 by deforestation, isolation of forest, and the 1-km edge effect of deforestation. The largest contributor to the area of negative effects on biological diversity was the 1-km edge effect from adjacent areas of deforestation. Isolation of forest patches was not a large contributor to this problem. The affected-habitat data were averaged into 16 km by 16 km grid cells.


Deforestation and Forest Fragmentation

Distribution of deforestation and affected habitat in the Brazilian Amazon for 1978 and 1988 (Figs. 2 and 3) was concentrated in a crescent along the southern and eastern fringe of the Amazon [a spatial pattern similar to the distribution of fires observed from thermal anomalies in data from Landsat's Advanced Very-High Resolution Radiometer (AVHRR) (20)] and along major transportation corridors in the interior of the Amazon. Deforestation increased between 1978 and 1988 (78,000 to 230,000 km2), while the total affected habitat increased (208,000 to 588,000 km2) (Table 3). The total area deforested increased by a factor of two to three or more in every state except Amapa; but it is likely that the deforested area in Amapa is higher than our assessment because excessive cloud cover in this region prevented complete analysis (Table 2). We found that 6% of closed-canopy forest had been cleared as of 1988 and ~15% of the forested Amazon was affected by deforestation-caused habitat destruction, habitat isolation, and edge effects (Fig. 2 and Table 3).


Table 3. Tropical deforestation, forest isolated or cut off by deforestation, and the area of forest adversely affected by a 1-km edge effect from adjacent areas of deforestation in the Brazilian Amazon. Areas that were obscured by clouds were omitted from this analysis. Parentheses following the edge effect entries contain the ratio between a 500-m buffer and a 1,000-m buffer.


Our analysis of the spatial pattern of deforestation found a strong tendency toward spatial concentration; areas of undisturbed tropical forest tended to be sizable (Table 4). This is more pronounced than Table 4 indicates because many of the large areas of undisturbed tropical forest are contiguous among states.

For the entire Brazilian Amazon, our deforestation estimate is close to, but lower than, the estimates of Fearnside et al. (13) and the INPE (15) of ~280,000 km2 as of 1988. The difference is a result of three factors: (i) different stratification of forest, cerrado, and water; (ii) slightly different estimates of secondary growth, which is spectrally similar to intact forest in channel five; and (iii) positional accuracy, interpretation, and boundary generalization. We estimate that ~30,000 km2 of the difference is from a different evaluation of the forest-cerrado boundaries in Mato Grosso and Tocantins. By comparison, our analysis suggests that deforestation estimates based on coarse-resolution meteorological satellite data in the southern Amazon of Brazil have overestimated deforestation by ~50% (18, 23).


Table 4.Spatial characteristics of isolated and remaining tropical forest within the \ Legal Amazon as determined by analysis of 1988 Landsat Thematic Mapper imagery. Isolated forest refers to areas of forest <100 km2 surrounded by deforestaion. Remaining forest refers to tropical forest that has not been deforested and includes both isolated and larger areas of forest. Many of the largest remaining areas of tropical forest are contiguous amoung states. Areas affected by clouds were omitted from this analysis.


The average deforestation rate in the closed-canopy forests from 1978 to 1988 (~15,000 km2 year-1) (Table 3) is higher than the rate from 1975 to 1978 (3) but considerably lower than recent estimates (2, 4, 5, 20). Our estimates can be used in assessments of net flux of carbon from land clearing and biomass burning in the Brazilian Amazon. Current estimates of these fluxes have largely been based on model calculations with deforestation values much higher than we report. In addition, many deforested areas are in stages of regrowth following abandonment (26). If regrowth is widespread, estimates of the net flux of carbon should be further reduced because carbon accumulates in regrowing biomass.

The preponderance of affected habitat results from proximity to areas of deforestation (~341,000 km2 for a 1-km edge effect) and not from isolation of forest (~16,000 km2) or deforestation per se (~230,000 km2). While the rate of deforestation averaged ~15,000 km2 year-1 in the Brazilian Amazonia from 1978 to 1988, the rate of habitat fragmentation and degradation was ~38,000 km2 year-1. Implications for biological diversity are not encouraging and provide added impetus for the minimization of tropical deforestation.

References and Notes

1. R. P. Tucker and J. F. Richards, Global Deforestation and the Nineteenth Century World Economy (Duke Univ. Press, Durham, NC, 1983); J. F. Richards, Environment 26, 6 (1984); M. Williams, Prog. Hum. Geogr. 13, 176 (1989); in The Earth as Transformed by Human Action, B. L. Turner II et al., Eds. (Cambridge Univ. Press, Cambridge, 1990), pp. 179-201.

2. N. Myers, Clim. Change 19, 3 (1991).

3. "Los Recursos Forestales de la America Tropical," 32/6.1301-78-04, Tech. Rep. No. 1 (Food and Agriculture Organization of the United Nations, Rome, 1981); "Forest Resources of Tropical Africa," 32/6.1301-78-04, Tech. Rep. No. 2 (Food and Agriculture Organization of the United Nations, Rome, 1981); "Forest Resources of Tropical Asia," 32/6.1301-78-04, Tech. Rep. No. 3 (Food and Agriculture Organization of the United Nations, Rome, 1981).

4. "The Forest Resources of the Tropical Zone by Main Ecological Regions," Report to the United Nations Conference on Environment and Development by the Forest Resource Assessment 1990 Project (Food and Agriculture Organization of the United Nations, Rome, 1992). The FAO Forest Assessment 1990 Project has produced several reports, and estimates from them have varied considerably. The recent release of another report [P. Aldhous, Science 259, 1390 (1993)] provides slightly different estimates than those reported in 1992.

5. World Resources 1990-91: A Report by the World Resources Institute in Collaboration with the United Nations Environment Program and The United Nations Development Program (Oxford Univ. Press, New York, 1990).

6. J. H. C. Gash and W. J. Shuttleworth, Clim. Change 19, 123 (1991); R. A. Houghton, et al., Nature 316, 617 (1985); R. A. Houghton and D. L. Skole, in The Earth as Transformed by Human Action, B. L. Turner II, et al., Eds. (Cambridge Univ. Press, Cambridge, 1990) pp. 393-408; _________, D. S. Lefkowitz, J. For. Ecol. Manage. 38, 173 (1991); M. Keller, D. J. Jacob, S. C. Wofsy, R. C. Harriss, Clim. Change 19, 139 (1991); E. Salati, in The Geophysiology of Amazonia: Vegetation and Climate Interaction, R. E. Dickinson, Ed. (Wiley Interscience, New York, 1987), pp. 273-296; _________ and C. A. Nobre, Clim. Change 19, 177 (1991); E. Salati and P. B. Vose, Ambio 12, 67 (1983); Science 225, 129 (1984); J. Shukla, C. Nobre, P. Sellers, ibid. 247, 1322 (1990).

7. R. A. Houghton, Clim. Change 19, 99 (1991).

8. Climate Change 1992: the Supplementary Report to the IPCC Scientific Assessment, J. T. Houghton, B. A. Callander, S. K. Varney, Eds. (Intergovernmental Panel on Climate Change, Cambridge Univ. Press, Cambridge, 1992).

9. E. O. Wilson, in Biodiversity, E. O. Wilson and F. M. Peters, Eds. (National Academy Press, Washington, DC, 1988), pp. 3-20; A. H. Gentry, Proc. Natl. Acad. Sci. U.S.A. 85, 156 (1988). The immense biological diversity of tropical forests is difficult to comprehend. For example, ten selected 1-ha plots in Borneo contained 700 species of trees and 1 ha of tropical Peru contained 300 tree species. By comparison, 700 tree species occur in all of North America. "Species" is used to mean organisms that can breed freely with each other and conversely, cannot breed freely with other species.

10. G. T. Prance and T. S. Elias, Eds., Extinction is Forever (New York Botanical Garden, New York, 1982); R. Lewin, Science 234, 14 (1986); T. L. Erwin, Coleopt. Bull. 36, 74 (1982); Bull. Entomol. Soc. Am. 30, 14 (1983); P. R. Ehrlich and E. O. Wilson, Science 253, 758 (1991); E. O. Wilson and F. M. Peters, Eds., Biodiversity (National Academy Press, Washington, DC, 1988); W. V. Reid and K. R. Miller, Keeping Options Alive: The Scientific Basis for Conserving Biodiversity (World Resources Institute, Washington, DC, 1989); E. O. Wilson, Issues Sci. Technol. 2, 20 (fall 1985); T. L. Erwin, in Biodiversity, E. O. Wilson and F. M. Peters, Eds. (National Academy Press, Washington, DC, 1988), pp. 123-129; R. M. May, Science 241, 1441 (1988); A. H. Knoll, in Extinctions, M. H. Nitecki, Ed. (Univ. of Chicago Press, Chicago,

1984), pp. 21-68. 11. J. F. Franklin and R. T. T. Forman, Landscape Ecol. 1, 5 (1987); L. D. Harris, The Fragmented Forest: Island Biogeographic Theory and the Preservation of Biotic Diversity (Univ. of Chicago Press, Chicago, 1984); L. D. Harris, Conserv. Biol. 2, 330 (1988); D. H. Janzen, Oikos 41, 402 (1983); T. E. Lovejoy et al., in Extinctions, M. H. Nitecki, Ed. (Univ. of Chicago Press, Chicago, 1984), pp. 295-325; T. E. Lovejoy et al., in Conservation Biology: The Science of Scarcity and Diversity, M. E. Soule, Ed. (Sinauer, Sunderland, MA, 1988), pp. 257-285; D. S. Wilcove, C. H. McLellan, and A. P. Dobson, ibid., pp. 273-256; D. A. Saunders, R. J. Hobbs, C. R. Margules, Conserv. Biol. 5, 18 (1991); J. Terborgh, Science 193, 1029 (1976); B. A. Wilcox, in Conservation Biology an Evolutionary-Ecological Perspective, M. E. Soule and B. A. Wilcox, Eds. (Sinauer, Sunderland, MA, 1980), pp. 95-117; K. H. Redford, BioScience 42, 412 (1992).

12. A. T. Tardin et al., Rep. 411-NTE/142 (Instituto Nacional de Pesquisas Espaciais, S o Jose dos Campos, Brazil, 1979); A. T. Tardin et al., "Subprojeto Desmatamento," IBDF/CNPq-INPE (Instituto de Pesquisas Espaciais, Sao Jose dos Campos, Brazil, 1980).

13. P. M. Fearnside, A. T. Tardin, L. G. M. Meira, Deforestation Rate in Brazilian Amazonia (National Secretariat of Science and Technology, Brasilia, Brazil, 1990).

14. A. T. Tardin and R. P. da Cunha, Report INPE-5015-RPE/609 (Instituto de Pesquisas Espaciais, Sao Jose dos Campos, Brazil 1990).

15. Deforestation in Brazilian Amazonia (Instituto Nacional de Pesquisas Espaciais, Sao Jose dos Campos, Brazil, 1992).

16. D. J. Mahar, Government Policies and Deforestation in Brazil's Amazon Region (World Bank, Washington, D.C., 1989).

17. P. M. Fearnside, Environ. Conserv. 17, 213 (1990).

18. A. M. Cross, J. J. Settle, N. A. Drake, and R. T. M. Paivinen, Int. J. Remote Sensing 12, 1119 (1991); J. P. Malingreau and C. J. Tucker, Ambio 17, 49 (1988).

19. R. F. Nelson and B. N. Holben, Int. J. Remote Sensing 7, 429 (1986); R. F. Nelson, N. Horning, T. A. Stone, ibid. 8, 1767 (1987); C. J. Tucker, B. N. Holben, T. E. Goff, Remote Sensing Environ. 15, 255 (1984); G. M. Woodwell, R. A. Houghton, T. A. Stone, R. F. Nelson, W. Kovalick, J. Geophys. Res. 92, 2157 (1987).

20. A. W. Setzer and M. C. Pereira, Ambio 20, 19 (1991); A. W. Setzer, Relatorio INPE- 4534-RPE/565, (Instituto Nacional de Pesquisas Espaciais, Sao Jose dos Campos, Brazil, 1988).

21. Anuario Estatistico do Brasil 1991 (Fundac o Instituto Brasileiro de Geografia e Estat¡stica, Rio de Janeiro, Brazil, 1991), vol. 51, pp. 1-1024.

22. Mapa de Vegetac o do Brasil (Fundac o Instituto Brasileiro de Geografia e Estat¡stica, Rio de Janeiro, Brazil, 1988).

23. D. L. Skole, thesis, University of New Hampshire (1992).

24. Images: 7 from 1989, 175 from 1988, 8 from 1987, and 20 from 1986. All data from the Brazilian Landsat receiving station. The exact boundary between intact forest and deforested land was digitized in the Universal Transverse Mercator projection and then edited and error-checked with use of clear velum plots of the line-work overlaid on each photographic image. Each Landsat scene contained coordinate control points in decimal degree units, such that each scene could be geographically registered within precise tolerances and mosaicked together. For digitization, vertices were placed approximately every 50 m of ground position. Tests of positional accuracy in digitizing followed those of R. Dunn, R. Harrison, and J. C. White [Int. J. Geograph. Inf. Syst. 4, 385 (1990)] and indicated encoding; hence, area-estimation errors were less than 3% (23). The variance associated with interpretation and delineation of boundaries between intact forest and deforested areas was less than 10% overall. Further accuracy assessment was made in test sites established in Rondonia, where fragmentation was very high. An explicit spatial comparison between our estimate of deforestation and the same derived from high-resolution (20-m resolution) SPOT satellite imagery was highly correlated (r2 = 0.98; y = 1.11x 57.358). Additional ground checking and verification was done in eastern Para state, north of Manaus, and along the Rio Negro, both in Amazonas.

25. Fundamental to our analysis was a specified representation for water, cerrado or savanna, and forest for the Brazilian Amazon. We used a vegetation map (23) that was augmented by Landsat Thematic Mapper and meteorological satellite imagery for more accurate depiction of cerrado and water. This GIS representation is available upon request.

26. R. B. Buschbacher, Bioscience 36, 22 (1986); C. Uhl, R. B. Buschbacher, E. A. S. Serrao, J. Ecol. 76, 663 (1988); R. B. Buschbacher, C. Uhl, E. A. S. Serrao, ibid., p. 682.

27. This work was supported by National Aeronautics and Space Administration's mission to planet Earth and the Eos Data Information System's Landsat Pathfinder Program. We acknowledge S. Tilford and W. Huntress for initiating this research, W. Chomentowski for assistance in developing the satellite and GIS database, and A. Nobre for his assistance in interpreting the satellite data. G. Batista, M. Heinicke, and T. Grant assisted with the GIS representation of forest, water, and cerrado.