EFFECTS OF URBAN SPRAWL ON FOREST CONSERVATION IN A METROPOLITAN WATER SOURCE AREA

EFEITOS DA EXPANSÃO URBANA NA CONSERVAÇÃO DE FLORESTAS EM UMA ÁREA DE MANANCIAL METROPOLITANO

ABSTRACT

The clearing of native vegetation in water source areas and its effects on water supply availability have become a major concern in large metropolises. This study examined the spatial and temporal changes in land use and cover and their effects on the conservation of Atlantic Forest in a water source area in the Metropolitan Region of São Paulo that is the second most densely populated territory in Brazil, the municipality of Diadema. Historical maps and trend scenarios were developed by applying the Markov chain in association with cellular automata. Forest conservation was assessed based on a set of landscape metrics of resource availability, fragment isolation, and trends of forest fragmentation. Our findings revealed an increase in urban growth of 37% in the last five decades directly associated to the clearing of 50% of forests in intermediate/advanced stages and 40% in early stages. The main drivers of landscape change were human-modified fields (41%) and net roads (54%). Optimal resource availability was drastically reduced (75%) while isolation increased (361.5%) among existing forest fragments. Therefore, future trends indicate that, without changes in the management of the area, deforestation may occur rapidly and continuously until 10% of forest remains, leading to potential adverse effects on water availability and quality of the Billings Reservoir and the local biodiversity.

Keywords:
Landscape metrics; Land use and cover; Billings reservoir

RESUMO

A supressão da cobertura vegetal em áreas de mananciais e seus efeitos sobre a disponibilidade de água para o abastecimento público têm se tornado cada vez mais uma preocupação das grandes metrópoles. Este estudo avaliou as mudanças espaço-temporais de uso e cobertura da terra e seus efeitos na conservação de florestas de Mata Atlântica em uma área de manancial de restrição hídrica para o abastecimento da Região Metropolitana de São Paulo, inserida na segunda maior densidade demográfica do Brasil, no município de Diadema. Foram elaborados mapas históricos e cenários de tendência futura por meio da aplicação da cadeia de Markov em associação com autômatos celulares. A conservação de florestas foi avaliada a partir de um conjunto de métricas da paisagem embasadas na disponibilidade de recursos, no isolamento e na tendência de fragmentação florestal. Os resultados confirmaram que houve incremento do uso urbano de 37% nas últimas cinco décadas, que está diretamente relacionado com o desmatamento de 50% de floresta em estágio médio/avançado e 40% em estágio inicial, sendo os campos antrópicos (41%) e as vias de acesso (54%) as principais forças motrizes de mudança da paisagem. Houve uma drástica redução na disponibilidade de recurso ótimo (75%) e no aumento do isolamento (361,5%) dos fragmentos florestais atuais. Como consequência, as tendências indicam que, se não houver mudanças na gestão da área, o desmatamento pode acontecer de forma acelerada e ininterrupta, restando apenas 10% de floresta o que pode, potencialmente, prejudicar a disponibilidade e qualidade hídrica do reservatório e a biodiversidade local.

Palavras-Chave:
Métrica da paisagem; Uso e ocupação da terra; Represa Billings

1. INTRODUCTION

Urban sprawls in water source areas in large cities, such as São Paulo, have caused concern regarding the availability and quality of public water supply, requiring increasingly farther sources (ANA, 2015; Alvim et al., 2015). Historically, the occupation of environmentally sensitive sites, such as water sources, occurs as a result of urban expansion from central areas to the outskirts of cities, where low-income populations settle in areas devoid of public services, often in a disorderly and irregular manner (Duarte and Malheiros, 2012). One of the main issues caused by these settlements is the clearing of the vegetation cover, which can result in siltation of rivers and reservoirs, hinder water infiltration into the soil and pollutant removal, compromising the water quality of the catchment area of water sources (Liu, 2015).

Given the effects of vegetation clearing, the quantity and quality of remaining forest fragments in water source areas are expected to reflect their environmental conservation status (Biao et al., 2010). According to principles of landscape ecology, the processes associated with environmental conservation can be inferred from patterns of vegetation composition and distribution and land uses (Turner et al., 2001). Based on this approach, structural patterns can be assessed by using landscape metrics or indexes (McGarigal and Marks, 1995), such as quality of the forests based on size, shape, isolation, succession stage, and anthropic boundaries (Zeng and Wu, 2005; Hardt et al., 2013a, 2014).

The current state of environmental conservation is also associated to land use history, which can be characterized based on the identification and understanding of the driving forces of changes observed over time (Bürgi et al., 2004; Plieninger et al., 2016). Understanding the dynamics of these changes can assist planning and assessing future conservation actions, such as the construction of trend scenarios to prevent or mitigate adverse impacts (Ferraz and Vettorazzi, 2003; Willis et al., 2007; Whitlock et al., 2017).

Given the increasing urban pressure of large metropolises, it is necessary to understand its effects on forest conservation. Therefore, studies are needed on the application of methods and techniques for the conservation of sensitive areas, such as water sources, in order to ensure the supply of water, especially in densely populated areas. In this sense, this study was aimed to evaluate the history and future trends of spacetime changes in land use and cover and their effects on the conservation of the Atlantic Forest in an area of water scarcity, in the second most demographically dense area of Brazil, the municipality of Diadema. Our findings may be used as a framework for planning and managing water sources, as it reveals the consequences of urban sprawling on the remaining forest fragments based on changes in landscape structural patterns.

2. MATERIALS AND METHODS

2.1 Study area

This study was conducted in the southern area of the municipality of Diadema in the Metropolitan Region of São Paulo (SP), Southeastern Brazil. The population projection of the municipality for 2017 was 418 thousand inhabitants in a small territory of 30.7 km² -13.1 hab./km² (IBGE, 2017). Of these, 7.2 km² are located over a water source protection area of the Billings Reservoir (São Paulo, 2009), where the population has increased considerably in recent decades and modified the native Atlantic Forest formed by remnants of Dense Ombrophilous Forest (São Paulo, 2010).

The study area encompasses 7.9 km² and consists of 12 micro-watersheds, tributaries of the Billings Reservoir. The limits of the area were determined with the software ArcGIS^® (version 10.2.2), using the municipality limits to the south and the watershed delimitation to the north (Figure 1).

Figure 1
Location and delimitation of the study site in the water source area of the Billings reservoir, in Diadema municipality, state of São Paulo - Brazil.
Figura 1
Localização e delimitação da área de estudo na região de manancial da represa Billings no município de Diadema, estado de São Paulo - Brasil.

2.2 Land use and cover maps and scenarios

Historical maps and future scenarios of land use and cover were elaborated for the area of the Billings Reservoir in Diadema.

We mapped the historical changes in ArcGIS^® using manual vectorization based on visual assessment of aerial orthophotos of 1962 and 2011, provided by the Municipality of Diadema, at a scale of 1:25,000 and 1:10,000, respectively. Twelve types of land use and cover were categorized based on region characteristics, using as criteria: pattern, shape, texture, granularity, and color tone of the photos.

Future scenarios wereconstructed from the projection of trends of land use and cover using the Markov chain, an empirical mathematical model that assumes that the relationship among present variables will be maintained in the future (Baker, 1989). According to this model, the condition of the system after time t+1 depends only on the state of the system at time t and those likely to occur, represented in transition matrices (Eastman, 2012). The future prediction was obtained using the software IDRISI^® (Taiga version) based on the probability of transition between the years 1962 and 2011, with the application of the Markov chain with a proportional error of 0.15, considering that the precision of the maps of use and coverage was always higher than 85%. Predictions were made at 10-year intervals (2021 to 2060), until the maximum interval of the historical analysis (1962-2011), which was 49 years (2060). Sincethe product generated by the Markovian analysis does not provide a single image of the representation of the future probabilities, the cellular automata technique was used to create trend scenarios with the CA_Markov tool of the same software with neighborhood rules “5x5”. In this application, also was used as criteria the number of iterations equal to the time interval of each scenario, ranging from 10 to 49 iterations, following Eastman (2012).

The historical and trend changes in land use and cover were calculated using IDRISI^® by overlaying maps to create cross-tabulation (CROSSTAB tool). In order to facilitate the visualization and interpretation of these changes, graphs were generated in Excel^® to show the percentage of losses and gains of a certain use/cover in relation to the other types.

2.3 Selection and application of landscape metrics

The different maps and scenarios were compared by selecting a set of landscape indexes or metrics to assess the conservation of forests. Thus, quantity and quality of forests were estimated considering the composition and configuration of landscape elements from a set of metrics based on: i) Optimal Resource Availability (ORA), calculated as a function of the Forest Ratio (FR), Size (FS) Shape (FSh), and Successional Stage (FSS); ii) analysis of landscape connectivity using the Euclidean Nearest Neighbor distance (ENN) and Proximity Index (PROX) measurements and; (iii) analysis of forest fragmentation trends, based on frequency, measured as the Edge Segments Density (ESD), and length, calculated as the Edge Density (ED) for the total boundaries between forest and the different types of human-related uses (Table 1).

Table 1
Criteria used to assess the forest conservation in historic landscapes and future scenarios. Modified from Hardt et al. (2013a, 2014).
Tabela 1
Critérios de avaliação da conservação da floresta nas paisagens históricas e nos cenários futuros. Modificado de Hardt et al. (2013a, 2014).

	Metric	Formula	Description	Software
Forest resource availability variable from 0 to 1	Forest Ratio -FR	$\mathit{FR}=\frac{\mathit{FA}}{\mathit{TA}}$	Forest cover area (FA) per total landscape area (TA)	ArcGIS®
	Forest Size - FS	$\mathrm{FS}=\frac{{\Sigma }_{\mathit{i}=l}^{\mathit{n}}\left({\mathit{A}}_{\mathit{i}}.{\mathit{CS}}_{\mathit{i}}\right)}{\mathit{FA}.{\mathit{CS}}_{\mathit{max}}}$	Patch forest area (Ai) weighted by coefficient of patch size (CSi), ranging from 0 to 4: 0 (≤ 0,01 ha), 1 (≤ 10 ha), 2 (10 - 50 ha), 3 (50 e 4 (≥ 100 ha), in relation to the maximum size (CSmax), adapted from Hardt et al. (2014)	ArcGIS®
	Forest Shape - FSh	$\mathit{FSh}=\frac{{\Sigma }_{\mathit{i}=l}^{\mathit{n}}\left({\mathit{CA}}_{\mathit{i}}\right)}{\mathit{FA}}$	Proportion of forest core area (CAi), considering a 70 m edge width, based on data obtained by Hardt et al. (2013b) for Atlantic Forest fragments with urban boundary	ArcGIS® V-LATE extension
	Forest Successional Stage - FSS	$\mathit{FSS}=\frac{{\mathit{FA}}_{\mathit{adv}}}{\mathit{FA}}$	Proportion of forest in the most advanced successional stage (FAadv), identified during the mapping based on aerial photographs differences roughness	ArcGIS®
	Optimal Resource Availability - ORA	$\mathit{ORA}=\frac{{\Sigma }_{\mathit{i}=l}\left({\mathit{CA}}_{\mathit{i}}.{\mathit{CS}}_{\mathit{i}}.{\mathit{CSS}}_{\mathit{i}}\right)}{\mathit{TA}.{\mathit{CS}}_{\mathit{max}}.{\mathit{CSS}}_{\max }}$	Forest quality, weighted by core area (CAi), size coefficient (CSi) and successional stages (CSSi = 1 for initial stage and 2 for medium/ advanced stage) in relation to the maximum potential in the total area (TA)	ArcGIS®
Landscape connectivity (McGarigal and Marks, 1995)	Euclidean Nearest Neighbor distance - ENN	$\mathit{ENN}={\mathit{h}}_{\mathit{ij}}$	Measure of focal patchi isolation by the shortest straight-line distance (hij) to nearest neighboring patchj of the same type, based on patch edge-to-edge distance	FRAGSTATS®
	Proximity Index -PROX	$\mathit{PROX}=\sum _{\mathit{s}=l}^{\mathit{n}}\frac{{\mathit{a}}_{\mathit{ijc}}}{{\mathit{h}}_{\mathit{ijs}}^2}$	Sum of patch forest area (aijs) divided by the nearest edge-to-edge distance squared (hijs) between the focal patch and all patches of the corresponding patch type whose edges are within 500 m from the focal patch
Trends to forest fragmentation	Frequency/Edge Segments Density - ESD	$\mathit{ESD}=\sum _{\mathit{k}=l}^{\mathit{n}}{\mathit{s}}_{\mathit{ik}}$	Number of boundary segments (sik) between forest patchi and human-related neighboring patchk (Zeng and Wu, 2005), calculated by the sum for boundary type	ArcGIS®
	Length/Edge Density - ED	$\mathit{ED}=\sum _{\mathit{k}=l}^{\mathit{m}}{\mathit{e}}_{\mathit{ik}}$	Cumulative total boundaries length (ED) between the forest patchi and human-related neighboring patchk, calculated in meters by the sum for boundary type (Mcgarigal and Marks, 1995; Zeng and Wu, 2005)	FRAGSTATS®

3. RESULTS

In 49 years, more than half of the Atlantic Forest surrounding the Billings Reservoir was cleared in Diadema, totaling 182 ha (Figure 2). The increase in urban uses observed in 2011 was 37% and accounts for 50% of the clearing of forest in intermediate/advanced stages and 40% in early stages (Figure 3A). This change is concentrated in the northern part of the study area, where net roads (54%), grouping of trees (50%) and human-modified fields (41%) were also frequent changes observed in urban areas (Figures 2 and 3A).

Figure 2
Historic maps and trend scenarios of land use and cover in the water source area of the Billings reservoir in the municipality of Diadema -SP.
Figura 2
Mapas históricos e cenários de tendência de uso e cobertura da terra na região de manancial da represa Billings no município de Diadema-SP.

Figure 3
Percentage of change in land use and cover types in the water source area of the Billings reservoir in the municipality of Diadema-SP: A -between the years 1962 and 2011; B -between the years 2011 and 2060. Bars below and above the zero axis indicate, respectively, the losses and gains of land use indicated in the X-axis (acronym) for the other types of use.
Figura3
Porcentagem de mudança nos tipos de uso e cobertura da terra na região de manancial da represa Billings no município de Diadema-SP: A - entre os anos de 1962 e 2011; B - entre os anos de 2011 e 2060. Barras abaixo e acima do eixo zero indicam, respectivamente, as perdas e os ganhos do uso indicado no eixo X (sigla) para os outros tipos de uso.

The future scenarios present a pessimistic perspective for the conservation of the region. Trends indicate an increase in forest loss, mainly as a result of urban sprawling to the north and west, towards the municipality’s water source areas, in the southeast portion of the study area (Figure 2). By 2060, over half of the current native forest will potentially be replaced by a urban high density, from the same types of land use changes - human-modified fields (37%), grouping of trees (29%), and net roads (24%) (Figure 3B). Another important force driving deforestation may come from the expansion of Eucalyptus spp. or Pinus spp. reforestation with an increase in monoculture plantations to the detriment of understory crops.

Since 1962, human-modified fields and net roads have been the two types of land use most often observed bordering the forest, both in frequency and in length (Figure 4A). These boundaries predominate in 2011 but with less intensity, as the forest begins change with other kinds of fragmentation pressures from new anthropogenic uses, such as reforestation with understory, residential with green spaces (residential areas with gardens and grouping of trees) and highdensity urban areas. In future scenarios, the boundary between forest and human-altered field will continue to predominate, but as these fields are transformed into high-density urban areas (Figure 2 and 3B), the urban-forest boundary increases primarily in frequency and can further intensify beginning in 2041 (Figure 4A).

Figure 4
Length and frequency of boundaries between forest patches and other land uses and cover (A) and isolation measurements, Proximity Index - PROX and Euclidean Nearest Neighbor distance - ENN, among the forest fragments (B) in the water source area of the Billings reservoir in the municipality of Diadema in the different years of examined in the study.
Figura 4
Comprimento e frequência de fronteiras entre as manchas de floresta e os outros usos e coberturas da terra (A) e medidas de isolamento, Índice de Proximidade - PROX e Distância Euclidiana do Vizinho mais Próximo - ENN, entre os fragmentos de floresta (B) da região de manancial da represa Billings no município de Diadema nos diferentes anos de estudo.

Although in 1962 the forest covered nearly half (43%) of the study area, had high successional stage -FSS (0.91) and forest fragments were closer -PROX, the rural landscape of this decade had low optimal resources availability -ORA (0.02). This low ORA value was associated with the irregular shape -FSh (0.10) and size -FS (0.38) of forest fragments, most of them less than 10 ha. During the 49 years of analysis of historical changes (1962-2011), ORA further declined (0.005), representing a reduction of 75%, mostly as a result of the sharp decrease in the proportion of forest -FR (0.18), especially in fragments in more advanced successional stages -FSS (0.83). The distances from forest fragments -ENN increased dramatically (361%) between the years of 1962 and 2011 and will continue to increase until 2060. Another concerning situation for the maintenance of landscape connectivity are PROX values approaching zero in the 2021 scenario. Until this year, ORA values remain the same as those for 2011 (0.005) and then decline until 2051 (0.002). The expected future variations for ORA are explained by the further reductions in the proportion of forest -RF (0.14 to 0.10), successional stage -FSS (0.8 to 0.76) and size of fragments with areas less than 10 ha (Figure 2 and 4).

4. DISCUSSION

The urban high-density growth of Diadema occurred towards the outskirts of the city, mainly as a result of deforestation. Studies have shown that unplanned urban occupation is directly associated with clearing of the vegetation cover, and the processes of urbanization and deforestation are often simultaneous events in order to meet the demand for new spaces (Weber and Puissant, 2003; Atmis et al., 2007; Defries et al., 2010; Martins, 2011). In the case of Diadema, at the municipal and state levels, the protection of the water sources is regulated by the State Law 13.579/ 09 of APRM-B -Area of Protection and Recovery of water sources of the Billings Reservoir (São Paulo, 2009) and by the complementary Municipal Law 273/ 08 of the Master Plan (Diadema, 2008). According to the environmental zoning of these legal provisions, most of the areas cleared for urbanization was not officially intended for urban expansion and were occupied irregularly. In large Brazilian metropolises, the exponential increase of the urban population, unemployment, and real estate speculation do not provide adequate housing alternatives for the low-income population (Canettieri, 2015). As observed in the study area, this population tends to be restricted to the outskirts of cities, with lower-cost housing options, usually located in irregular lots in areas of environmental protection and water sources and devoid of urban infrastructure and environmental sanitation (Jacobi, 2004; Canettieri, 2015).

The analysis of changes in land use and cover revealed the driving forces responsible for shifting the rural landscape from the 1960s to the recent urbanization (2011). The presence of grouping of trees, human-modified fields, and net roads facilitate the emergence of new residential and urban uses (Hardt et al., 2013a). Studies show that access to previously unoccupied areas allows new opportunities for economic activities that, in turn, generate income and become places of interest to live (Chomitz and Gray, 1996; Freitas et al., 2010; Barber et al., 2014). Our finding about how grouping of trees and human-modified fields function as drivers of urbanization is very relevant for the management of green spaces in large urban centers, such as São Paulo. It means that squares and gardens, represented by these uses, can be facilitators for the occupation of new irregular uses (Ferraz and Machado, 2014).

Future projections indicate that, without improvements in the environmental management of the area, continuous and rapid deforestation can occur. Therefore, management changes are a priority, as urban expansion continues toward the water source of the Billings Reservoir, an environmentally sensitive area of the municipality. Pollution generated by the direct dumping of sanitary effluents in the water body is a concern, as well as and the clearing of riparian forest, increasing soil vulnerability and the risks of landslides and floods (Jacobi,2006;Anelli, 2015).Forests protect soils against erosion and rivers and reservoirs from silting, favoring infiltration and reducing runoff, which ensures the recharge of water sources and aquifers and helps regulate water flow during times of flood and drought (Neary et al., 2009; Liu, 2015). In addition, the conservation of riparian forests is able to filter the environmental pollutants, contributing to the maintenance of adequate conditions for aquatic biodiversity (Neary et al., 2009; Pereira-Silva et al., 2011). Thus, the trend of replacing the forest cover of water source areas for urban high-density uses and the lack of planning can adversely impact the maintenance of the water quality and quantity from the Billings Reservoir.

Boundary metrics are usually used as indicators of connectivity or landscape fragmentation (Metzger and Muller, 1996). Increases in frequency and length of edges between human-related activities and forest remnants promote landscape fragmentation (Zeng and Wu, 2005). Therefore, the dominance of forests boundaries with human-modified fields and net roads since 1962 supports the influence of these uses on the historical changes of forest loss and fragmentation and the emergence of the present urban area. A new contact zone between the forest and the urban matrix is created, increasing pressure on the forest due easier access and higher human circulation, such as trails, improper waste disposal, vandalism, and wood collecting (Matlack, 1993). The forest-urban contact zone also has the potential to increase the extent of the natural edge effect (Hardt et al., 2013b), which may result in changes in the composition and density of fauna and flora, as well as invasion or dominance of exotic species that are more competitive than less tolerant species (Fahrig, 2003; Ries et al., 2004; Villaseñor et al., 2014). As a result, urban areas become new agents of fragmentation and to a higher extent, responsible for the drastic reduction in the optimal resource availability - ORA observed in this study.

The results obtained for 1962 support the assumption that quantity of forest without habitat quality does not always promote environmental conservation (Fahrig, 2001; Pflüger and Balkenhol, 2014; Magura et al., 2017). Small, elongated, and irregularly shaped forest patches are the main limitations for the most demanding species that are restricted to the interior of forests because they are not able to survive under the environmental conditions at edge microclimates (Wiens et al. 1993; Fahrig, 2003; Kahilainen et al., 2014).

The current data (2011) show that the forest cover of the water source area of Diadema was considerably reduced, with dramatic qualitative changes. Future projection scenarios indicate that the process of landscape fragmentation will advance and as a result, forest fragments will become more disconnected and with lower environmental quality. The physical and functional reach between habitats, variable with home range, dispersal capacity, and species migration (Matos, et al., 2016; Pflüger and Balkenhol, 2014) represent another concerning factor for the survival of local species and ecologic groups.

Given the effects of conservation trends, the first improvement for the conservation of Diadema forests is an increase in forest cover in order to protect water sources and consequently the public water supply. However, our findings emphasize that special attention should be giving to the size and shape of fragments and the conduction of studies for the identification of priority areas for the restoration of the functional connectivity of the landscape and facilitation of the genetic flow of local species.

5. CONCLUSIONS

The analysis of the historical changes in land use and cover in the water source area of the Billings Reservoir in Diadema revealed that deforestation in the last five decades is directly associated to urban sprawling in the municipality, resulting in great damage to the local conservation of the Atlantic Forest. Consequently, there is a clear trend towards even further reductions in quantity, quality, and connectivity of forest fragments, which can have adverse effects on water availability and quality of the reservoir and be a potential threat to the survival of local native species. Reversing this scenario depends on changes in planning and management aimed at ensuring the maintenance of the remaining forest remnants, as well as incentives to viable alternatives for restoring the native vegetation.

The main scientific contribution of our study is to demonstrate the utility of assessment methods of forest conservation to understand the effects of the urban expansion densely populated in water source areas with irregular occupations, and their potential to be applied in landscape planning and management.

The main challenge of this methodological application, however, is the lack of historical geospatial data and in-depth studies about the effects of urbanization on the conservation of biological diversity, especially of the species and ecological groups that occur in this region.

6. ACKNOWLEDGMENTS

This study was supported by the Conselho Nacional de Desenvolvimento Científico e Tecnológico - CNPq - Universal Grant Program 01/2016 [proc. 426178/20165] and Science Intern Fellowship Program - PIBIC/ CNPq/Unifesp [proc. 152593/2015-2 and 101335/20172]. The authors thank the Geoprocessing Office of Diadema City for providing the original spatial dataset of the municipality.

7. REFERENCES

• Agência Nacional de Águas-ANA. Atlas Brasil (Abastecimento Urbano de Água). Apresenta o trabalho de diagnóstico e planejamento nas áreas de recursos hídricos e saneamento no Brasil. [acessado em: 03 set. 2015]. Disponível em: http://atlas.ana.gov.br/Atlas/forms/AglomeradosUrbanos.aspx
  » http://atlas.ana.gov.br/Atlas/forms/AglomeradosUrbanos.aspx
• Alvim ATB, Kato VRC, Rosin JR. A urgência das águas: intervenções urbanas em áreas de mananciais. Cadernos Metrópole. 2015;17(33):83-107.
• Anelli RLS. Uma nova cidade para as águas urbanas. Estudos Avançados. 2015;29(84):69-84.
• Atmis E, Özden S, Lise W. Urbanization pressures on the natural forests in Turkey: An overview. Urban Forestry & Urban Greening. 2007;6(2):83-92.
• Baker WL. A review of models of landscape change. Landscape Ecology. 1989;2(2):111-33.
• Barber CP, Cochrane MA, Souza CM, Laurance WF. Roads, deforestation, and the mitigating effect of protected areas in the Amazon. Biological Conservation. 2014;177:203-9.
• Biao Z, Wenhua L, Gaodi X, Yu X. Water conservation of forest ecosystem in Beijing and its value. Ecological Economics. 2010;69(7):1416-26.
• Bürgi M, Hersperger AM, Schneeberger N. Driving forces of landscape change - current and new directions. Landscape Ecology. 2004;19(8):857-68.
• Canettieri T. Uma nova segregação metropolitana: as periferias fractais - evidências encontradas na Região Metropolitana de Belo Horizonte. Revista Brasileira de Estudos Urbanos e Regionais. 2015;17(2):25-39.
• Chomitz MK, Gray AD. Roads, land use, and deforestation: A spatial model applied to belize. World Bank Economic Review. 1996;10:487-512.
• Defries R, Rudel T, Uriarte M, Hansen M. Deforestation driven by urban population growth and agricultural trade in the twenty-first century. Nature Geoscience. 2010;3:178-81.
• Diadema. Prefeitura Municipal. Lei complementar nº 273, de 08 de jul de 2008. Dispõe sobre o Plano Diretor do Município de Diadema. [acessado em: 23 out 2015]. Disponível em: http://www.cmdiadema.sp.gov.br/legislacao/leisintegra.php?chave=10027308
  » http://www.cmdiadema.sp.gov.br/legislacao/leisintegra.php?chave=10027308
• Duarte CG, Malheiros TF. Habitação e gestão ambiental em áreas de mananciais: o caso do município de Santo André (SP). Saúde e Sociedade. 2012;21(3):82-95.
• Eastman JR. IDRISI Selva: Tutorial. 17ªth. ed. Massachusetts: Clark University; 2012. p.182-5.
• Fahrig L. How much habitat is enough? Biological Conservation. 2001;100(1): 65-74.
• Fahrig L. Effects of habitat fragmentation on biodiversity. Annual Review of Ecology, Evolution and Systematics, 2003;34:487-515.
• Ferraz SFB, Vettorazzi CA. Identificação de áreas para recomposição florestal com base em princípios de ecologia de paisagem. Revista Árvore. 2003;27(4):575-83.
• Ferraz MTS, Machado BA. Eu não tenho onde morar, é por isso que moro na rua: Os “sem-teto”: moradores ou transgressores? Cadernos Metrópole. 2014;16(32):609-23.
• Freitas SR, Hawbaker JT, Metzger PJ. Effects of roads, topography, and land use on forest cover dynamics in the Brazilian Atlantic Forest. Forest Ecology and Management. 2010;259:410-7.
• Hardt E, Santos RF, Pablo CL, Agar PM, Pereira-Silva EFL. Utility of landscape mosaics and boundaries in forest conservation decision making in the Atlantic Forest of Brazil. Landscape Ecology. 2013a;28(3):385-99.
• Hardt E, Pereira-Silva EFL, Santos RF, Tamashiro JY, Ragazzi S, Lins D. The influence of natural and anthropogenic landscapes on edge effects. Landscape and Urban Planning. 2013b;120:59-69.
• Hardt E, Santos RF, Pereira-Silva EFL. Evaluating the ecological effects of the social agent scenarios for a housing development in the Atlantic Forest. Ecological Indicators. 2014;36:120-130.
• Instituto Brasileiro de Geografia e Estatística - IBGE. Apresenta estimativas populacionais dos municípios para 2017 [acessado em: 8 dez. 2017]. Disponível em: https://agenciadenoticias.ibge.gov.br/agencia-detalhe-demidia.html?view=mediaibge&catid=2103&id=1328
  » https://agenciadenoticias.ibge.gov.br/agencia-detalhe-demidia.html?view=mediaibge&catid=2103&id=1328
• Jacobi P. Cidade e meio ambiente: percepções e práticas em São Paulo. 2ª. ed. São Paulo: Annablume; 2006. p.46-79.
• Jacobi P. The challenges of multi-stakeholder management in the watersheds of São Paulo. Environment & Urbanization. 2004;16(2):199-211.
• Kahilainen A, Puurtinen M, Kotiaho JS. Conservation implications of species-genetic diversity correlations. Global Ecology and Conservation. 2014;2:314-23.
• Liu E. Hydrological processes of water resource conservation forest. In: Zhang J. editor. Environmental protection and sustainable ecological development. London: Taylor & Francis Group; 2015. p.175-9.
• Magura T, Lövei GL, Tóthmérész B. Edge responses are different in edges under natural versus anthropogenic influence: a meta-analysis using ground beetles. Ecology and Evolution. 2017;7(3):1009-17.
• Martins MLR. São Paulo, centro e periferia: a retórica ambiental e os limites da política urbana. Estudos Avançados. 2011;25(71):59-72.
• Matlack GR. Sociological edge effects: Spatial distribution of human impact in suburban forest fragments. Environmental Management. 1993;17(6):829-35.
• Matos FAR, Magnago LFS, Gastauer M, Carreiras JMB, Simonelli M, Meira-Neto JAA. et al. Effects of landscape configuration and composition on phylogenetic diversity of trees in a highly fragmented tropical forest. Journal of Ecology. 2016;105:265-76.
• McGarigal K, Marks B. Fragstats: spatial pattern analysis program for quantifying landscape structure. Corvallis: Oregon State University; 1995. 134p.
• Metzger JP, Muller E. Characterizing the complexity of landscape boundaries by remote sensing. Landscape Ecology. 1996;11(2):65-77.
• Neary DG, Ice GG, Jackson CR. Linkages between forest soils and water quality and quantity. Forest Ecology and Management. 2009;258(10):2269-81.
• Pereira-Silva EFL, Pires JSR, Hardt E, Santos JE, Ferreira WA. Avaliação da qualidade da água em microbacias hidrográficas de uma Unidade de Conservação do Nordeste do estado de São Paulo, Brasil [Internet]. Revista Brasileira de Biociências. 2011;9:371-81.
• Pflüger FJ, Balkenhol N. A plea for simultaneously considering matrix quality and local environmental conditions when analysing landscape impacts on effective dispersal. Molecular Ecology. 2014;23(9):2146-56.
• Plieninger T, Draux H, Fagerholm N, Bieling C, Bürgi M, Kizos T. et al. The driving forces of landscape change in Europe: A systematic review of the evidence. Land Use Policy. 2016;57:204-14.
• Ries L, Fletcher RJ, Battin J, Sisk TD. Ecological responses to habitat edges: mechanisms, models, and variability explained. Annual Review of Ecological, Evolution and Systematics. 2004;35:491-522.
• São Paulo (Estado). Lei nº 13.579, de 13 de julho de 2009. Define a Área de Proteção e Recuperação dos Mananciais da Bacia Hidrográfica do Reservatório Billings -APRM-B. [acessado em: 10 out. 2015]. Disponível em: http://www.al.sp.gov.br/repositorio/legislacao/lei/2009/lei-1357913.07.2009.html
  » http://www.al.sp.gov.br/repositorio/legislacao/lei/2009/lei-1357913.07.2009.html
• São Paulo (Estado). Infraestrutura de Dados Espaciais Ambientais do Estado de São Paulo. Base territorial ambiental unificada - Inventário Florestal 2010. [acessado em: 10 out. 2016. Disponível em: http://datageo.ambiente.sp.gov.br/app/ctx=DATAGEO#
  » http://datageo.ambiente.sp.gov.br/app/ctx=DATAGEO#
• Turner MG, Gardner RH, O’Neill RV. Landscape ecology in theory and practice: pattern and process. New York: Spring-Verlag; 2001.
• Villaseñor NR, Driscoll DA, Escobar MAH, Gibbons P, Lindenmayer DB. Urbanization impacts on mammals across urban-forest edges and a predictive model of edge effects. PLoS ONE. 2014;9(5):e97036.
• Weber C, Puissant A. Urbanization pressure and modeling of urban growth: example of the Tunis Metropolitan Area. Remote Sensing of Environment. 2003;86(3):341-52.
• Wiens JA, Stenseth NC, Horne BV, Ims RA. Ecological Mechanisms and Landscape Ecology. Oikos. 1993;66(3):369-80.
• Willis KJ, Araújo MB, Bennett KD, Figueroa-Rangel B, Froyd CA, Myers N. How can a knowledge of the past help to conserve the future? Biodiversity conservation and the relevance of long-term ecological studies. Philosophical Transactions of The Royal Society B. 2007;362(2367):175-86.
• Whitlock C, Colombaroli D, Conedera M, Tinner W. Land-use history as a guide for forest conservation and management. Conservation Biology. 2017:1-14.
• Zeng H, Wu XB. Utilities of edge-based metrics for studying landscape fragmentation. Computers, Environment and Urban Systems. 2005;29:159-78.

Publication Dates

• Publication in this collection
  13 Aug 2018
• Date of issue
  2018

History

• Received
  19 Oct 2017
• Accepted
  21 Mar 2018