Evolução da vulnerabilidade à erosão costeira na Praia de Massaguaçú (SP), Brasil*

Evolution of vulnerability to coastal erosion at Massaguaçú Beach, Brazil

Juliana dos Santos Ribeiro ^{@, I}; Paulo Henrique Gomes de Oliveira Sousa ^I; Danilo Rodrigues Vieira ^I; Eduardo Siegle ^I

^@Corresponding author: juliana.santos.ribeiro@usp.br, jusantosribeiro@gmail.com

^IInstituto Oceanográfico - Universidade de São Paulo, Praça do Oceanográfico, 191. São Paulo, Brasil

RESUMO

A erosão costeira pode causar grandes prejuízos tanto ambientais quanto sociais nas áreas atingidas e a ocupação humana mal planejada nesses locais pode agravar o processo erosivo. A praia de Massaguaçú, localizada no litoral norte do Estado de São Paulo, vem sofrendo processos erosivos intensos nos últimos anos. O presente estudo teve como objetivo analisar a vulnerabilidade à erosão da praia em questão nos últimos 40 anos a partir de nove indicadores de vulnerabilidade: largura da praia, posição da linha de costa, configurações ao largo, presença de rios e/ou desembocaduras, elevação do terreno, vegetação, obras de engenharia costeira, taxa de ocupação e permeabilidade do solo. Os resultados mostraram que a vulnerabilidade à erosão não teve um aumento constante ao longo do período analisado. A parte central de Massaguaçú, hoje alvo de um processo erosivo intenso, foi a que apresentou a maior vulnerabilidade à erosão ao longo do tempo. Apesar de a ocupação humana aparentemente ter tido participação nos resultados, concluiu-se que a mesma pode não ser a causa direta do atual processo erosivo em Massaguaçú, uma vez que houve predominância de vulnerabilidades baixas mesmo em anos em que a ocupação esteve crescente. O método utilizado no estudo se mostrou eficiente para responder as questões propostas; é um procedimento rápido e de baixo custo para a avaliação do ambiente, tornando-o acessível para os gestores e tomadores de decisões realizarem avaliações rápidas do meio.

Palavras-chave: Indicadores de vulnerabilidade, Sistemas de Informação Geográfica (SIG), planejamento urbano.

ABSTRACT

Coastal erosion can cause major environmental and social damage to coastal areas, and disorderly and inadequate human occupation at these sites can lead to erosion.
Massaguaçú Beach, on the northern coast of São Paulo State, has suffered severe erosion in recent years. In the most severe cases, erosion had already begun to affect the road way Dr. Manuel Hyppolito Rego (SP-55) that runs next to the beach. Several efforts have been implemented to contain this problem, but have proved to be unsuccessful. The cause of these erosions remain unclear, however, studies suggest that these changes in the coast line may be related to the increasing urbanization around the beach.
Because of this phenomenon, this study aimed to analyze the coastal erosion vulnerability at the beach over a period of 40 years, with the objective of evaluating if human occupation had a significant impact in the local erosion. Data was gathered using nine vulnerability indicators: beach width, shoreline position, offshore settings, presence of rivers and/or inlets, terrain elevation, vegetation, coastal engineering structures, occupation rate and soil permeability.

Keywords: Indicators of vulnerability, Geographic Information System (GIS), urban planning.

1. INTRODUÇÃO

1.1.Motivação

A avaliação da vulnerabilidade à erosão costeira tem sido uma ferramenta importante para identificar e prevenir impactos socioeconômicos nas cidades litorâneas. Entretanto, a lacuna oriunda da escassez de dados de monitoramento de longo prazo da zona costeira, sobretudo no Brasil, é um obstáculo difícil de ser vencido. O esforço realizado nesse sentido tem partido de pesquisadores com interesse fundamentalmente científico.

Ao longo de praticamente todo o litoral brasileiro são verificados pontos de erosão costeira, resultantes de processos que, em muitos casos, são fruto de ações humanas que causam interferências na morfodinâmica e no aporte sedimentar (Muehe, 2005; 2006). Muehe (2005) não considera esses processos erosivos como uma ameaça para a totalidade da orla, contudo esses processos têm se traduzido em danos socioeconômicos, ambientais e culturais significativos na medida em que casas são instaladas na pós-praia, a infraestrutura urbana é afetada e obras de engenharia costeira são construídas.

De forma sumarizada, segundo Cai et al. (2009), as causas da erosão costeira de escala global ou regional podem ser naturais: aumento do nível do mar, intensificação de tempestades, subsidência tectônica, alterações nas bacias hidrográficas; e antrópicas: subsidência do solo e construção de barragens. À escala local, podem-se considerar relevantes as construções que interferem na deriva litoral, e a retirada de areia no sistema fluvial-litoral. Os processos naturais e antrópicos que resultam na erosão costeira inter-relacionam-se ao reduzir o balanço sedimentar ou alterar a dinâmica costeira. Além da perda de área emersa, outros impactos são, sobretudo, a intrusão marinha nas águas superficiais, a curto prazo, e a intrusão de água salgada no lençol freático e a perda da biodiversidade a longo prazo (Nicholls 2004; Nicholls & Cazenave, 2010).

Nesse contexto, a projeção de aumento do nível do mar (NM) feita pelo IPCC (2007) para 2100 foi de 0,2 a 0,6 m. Adicionalmente, existem projeções como a de Rahmstorf (2007) onde o cenário mais otimista para o mesmo período é de um aumento da ordem de 0,6 m e o mais pessimista é de 1,4 m. Esses estudos são importantes para fins de planejamento e adaptação. No caso de algumas praias onde um rápido aumento do NM foi observado, foi constatado que as ações de adaptação partem primeiramente da população afetada e pouco ou quase nada é feito pelo governo (Linnekamp et al., 2011).

Certamente as ações com vista à adaptação ao aumento do NM afetam diretamente a vida da população que habita as proximidades do litoral. No caso do delta de Ebro, Eha, Fatorić & Chelleri (2012) constataram que a comunidade mostrou preocupação em relação a como a adaptação seria feita. Esta optou por um processo suave, com ação realizada em consonância com a manutenção dos recursos naturais da região.

Sabe-se, no entanto, que os impactos humanos podem somar-se àqueles causados pela natureza. Além disso, sobretudo às escalas local e regional, podem ser mais danosos do que processos naturais que atuam em grande escala temporal e espacial, como as mudanças climáticas. A melhor forma de lidar com prejuízos causados por processos erosivos é tomar medidas preventivas antes que eles atinjam níveis mais críticos. Desta maneira, estudos que contemplem a avaliação da vulnerabilidade à erosão costeira assumem papel primordial para o planejamento urbano e regional em áreas litorâneas.

A vulnerabilidade da zona costeira pode ser avaliada considerando diferentes fatores, como o aumento do nível do mar, erosão costeira e tempestades. Vários trabalhos foram feitos considerando esses fatores (Pendleton et al., 2010; Bush et al., 1999; Bosom & Jiménez, 2011), no entanto, todos foram realizados para áreas de grande extensão geográfica e utilizando bancos de dados com vários anos de monitoramento. Alguns estudos de vulnerabilidade costeira já foram realizados na praia de Massaguaçú (SP) (Figura 1) (Sousa et al., 2011; Sousa et al., 2013). Entretanto, esses trabalhos fizeram uma análise instantânea da vulnerabilidade com base em informações recentes. Diferentemente, o objetivo deste trabalho é avaliar a evolução da vulnerabilidade à erosão costeira local ao longo de 40 anos observando como se deu o crescimento urbano e como este influenciou os processos erosivos atuais na região. O estudo tem como base uma série de indicadores de ordem, sobretudo, qualitativa, adaptados para a área em apreço: largura da praia, posição da linha de costa, configurações ao largo, presença de rios e/ou desembocaduras, elevação do terreno, vegetação, obras de engenharia costeira, taxa de ocupação e permeabilidade do solo. Através destes, foram estabelecidos graus de vulnerabilidade nos diferentes anos e setores do arco praial.

A praia de Massaguaçú (Figura 1), pertencente ao município de Caraguatatuba (litoral norte de São Paulo), está localizada próximo ao Trópico de Capricórnio, possui 7,5 km de extensão e orientação NE-SE. A praia é composta por areia grossa e média, com grau de selecionamento moderado (Rogacheski, 2010). A maior parte da praia apresenta morfologia com características de praia refletiva de alta energia, com tendências intermediárias (Sousa & Luna, 2010). A linha praial é interrompida por dois pequenos rios que alcançam o oceano nas extremidades da praia: rio Capricórnio ao sul e rio Bracuí ao norte.

A área está sujeita a um regime de micro-marés (Sousa et al, 2013). Quanto ao clima de ondas, a direção dominante da onda é leste no verão (com período de 6-8 s e altura de onda de 1-2 m) e sul nas demais estações (com períodos de 10-12 s e alturas de onda de 1-2 m, 1-3 m e 2-3 m para verão, outono e inverno, respectivamente), de acordo com os resultados de Pianca et al. (2010) para a costa sudeste brasileira.

2. MATERIAL E MÉTODOS

Foram utilizadas técnicas de geoprocessamento baseadas em Sistemas de Informações Geográficas (SIG). Para isso, a praia de Massaguaçú, com extensão de 7,5 km, foi dividida em três setores: norte, centro e sul. As técnicas utilizadas para este trabalho estão descritas a seguir.

2.1. Imagens aéreas

A maioria dos indicadores de vulnerabilidade utilizados para estimar a evolução da vulnerabilidade à erosão foram analisados com base em imagens digitais. Para tanto, utilizaram-se conjuntos de imagens aéreas abrangendo um período de aproximadamente 40 anos (Tabela 1).

As imagens foram utilizadas para determinar o grau de vulnerabilidade da área em cada ano e gerar mapas temáticos para uma melhor visualização dos indicadores analisados. O material foi fornecido pelo Instituto Oceanográfico da Universidade de São Paulo (fotografias aéreas dos anos de 1962, 1973, 1977 e 1987) e pelo INPE (Instituto Nacional de Pesquisas Espaciais; fotografia digital do ano de 2001).

A imagem de 2001 foi fornecida georreferenciada no sistema de coordenadas UTM (Universal Transverse Mercator), fuso 23S e datum planimétrico SAD69 (South American Datum 1969). Desta forma, essa imagem serviu como base para o georreferenciamento das demais que foram projetadas através dos programas IDRISI 15.0 – Andes Edition e ArcMAP 9.3 da ESRI (Environmental Systems Research Institute).

Durante o processo, o erro encontrado no georreferenciamento (RMS - Root Mean Square) foi sempre menor que o erro máximo permitido, que corresponde ao erro de graficismo da imagem (Júnior et al., 2006). A resolução das imagens após digitalização foi de 300 DPIs.

2.2. Delimitação dos setores e extensão da faixa costeira analisada

Para facilitar a análise, o arco praial de Massaguaçú, com extensão de 7,5 km, foi dividido em três setores distintos: setor norte, setor centro e setor sul. Essa divisão se deu em virtude das diferentes feições encontradas ao longo da área de estudo (Figura 1). No setor norte há a presença das Ilhas Cocanha, que servem de barreira natural para a praia, além da presença do rio Bracuí. O setor centro é caracterizado por conter a maior extensão da rodovia, sendo que esta margeia a praia por todo o comprimento do setor. Além disso, nessa área há ausência de rios e barreiras naturais ao largo. Já o setor sul é caracterizado pela quase ausência da rodovia, e pela presença do rio Capricórnio. A extensão da faixa costeira analisada foi de 500 m, partindo do limite da linha d’água (obtida através da imagem do ano de 2001) em direção ao interior.

2.3. Indicadores de Vulnerabilidade

Para se obter o grau de vulnerabilidade à erosão da área de estudo, nove indicadores foram analisados. Foram utilizados parâmetros qualitativos e quantitativos para classificar os indicadores como sendo de baixa, média ou alta vulnerabilidade (Tabela 2).

Os indicadores analisados foram adaptados de Sousa et al. (2013) e são brevemente descritos abaixo:

Posição da linha de costa: linhas de costa estáveis ou em progradação representam um cenário de baixa vulnerabilidade, enquanto que linhas de costa que sofrem retrogradação evidenciam um cenário de alta vulnerabilidade.

Largura da praia (ao longo do arco praial): sistemas praiais largos são mais propensos à estabilidade.

Configurações ao largo: a área de mar aberto que se encontra adjacente à linha de costa determina, parcialmente, a altura das ondas (de tempestade ou não) que podem surgir. Este indicador tem relação com feições migratórias ou fixas que, em geral, funcionam como barreiras que atenuam o ataque das ondas à zona costeira.

Elevação do terreno: a elevação do sistema praia-duna e pós-praia é um fator primordial quando se trata da suscetibilidade à inundação. Áreas baixas e planas estão expostas a ataques de ondas, sobrelavagem, inundação por tempestade (incluindo a ondulação e a sobreelevação metereológica), e às vezes transporte eólico de areia.

Vegetação: a presença de vegetação bem desenvolvida, seja ela rasteira, arbustiva ou arbórea na zona de pós-praia, sugere um cenário de baixa erosão e rara intrusão de água salina.

Presença de rios e/ou desembocaduras: essas feições podem migrar por centenas de metros, sugerindo instabilidade em suas adjacências.

Permeabilidade do solo: a permeabilidade dos solos está intimamente relacionada com a ocupação e com o tipo de uso do solo. Solos pouco permeáveis são mais susceptíveis à inundação.

Taxa de ocupação: a ocupação da zona costeira altera o ambiente natural, seja com o desmatamento de vegetação nativa, com a impermeabilização do solo, a criação de corredores de penetração de inundações ou com outras ações. Ademais, as estruturas construídas são muito afetadas por processos erosivos. Quanto maior a taxa de ocupação do local, maior a vulnerabilidade à erosão costeira.

2.3.1. Posição da linha de costa

A linha de costa (LC), do ponto de vista físico, corresponde ao ponto de encontro entre o continente e o oceano. Deve-se entender que a LC é uma linha móvel, em escalas espacial e temporal, devido aos inúmeros fatores que atuam nesse ambiente dinâmico (processos erosivos e deposicionais, marés, ondas, dentre outros; Araujo et al., 2009). Tendo em vista essa mobilidade, é importante que haja a definição de um critério para delimitação da LC, para que os resultados obtidos possam ser comparados.

São diversos os indicadores que podem ser utilizados para se delimitar a linha de costa, como o limite da vegetação costeira, a linha d’água e a interface da zona entremarés (correspondente ao limite máximo da preamar). A delimitação deve ser feita pelo indicador mais discernível e constante no material analisado (Araujo et al., 2009; Boak & Turner, 2005), que no caso das imagens aéreas utilizadas no presente estudo é a interface água - sedimento. Na área de estudo, em função de mudanças na linha de vegetação causadas pela ação antrópica (plantio nas margens da estrada), optou-se pela utilização da linha d’água como indicador de linha de costa. Apesar das possíveis limitações envolvidas na análise desse parâmetro, ele apresentou-se como melhor opção para a área de interesse.

As linhas de costa foram digitalizadas com o uso do ArcGIS 9.3. Para determinar a variação dessas linhas nos vários anos, uma rotina em ambiente MATLAB foi criada. A função consiste em analisar a variação de duas linhas de costa a partir de “n” transectos traçados perpendicularmente a uma linha de base suavizada. O retorno da função é uma matriz de duas colunas, sendo a primeira a coordenada do meridiano UTM do transecto e a segunda um valor positivo ou negativo representando a distância de progradação ou retrogradação da costa, respectivamente.

2.3.2. Largura da praia

A partir da imagem aérea de 2001, foram traçados 37 transectos ao longo da área de estudo, perpendiculares à linha de costa e distantes um dos outros cerca de 200 metros. A largura da praia de cada ano corresponde à medida do comprimento do segmento de cada transecto que tinha como limites a linha d’água e a faixa de vegetação. Para determinar o grau de vulnerabilidade de cada setor, fez-se uma divisão do comprimento dos transectos em três intervalos iguais, considerando como extremos o maior e menor valor encontrado ao longo dos anos. Assim, os setores que apresentaram um maior número de transectos estreitos (de 8,78 a 37,17 metros) foram classificados como de alta vulnerabilidade, os que apresentaram transectos de média extensão (37,18 a 65,56 metros) foram considerados de média vulnerabilidade e os setores que continham em sua maioria transectos largos (65,57 a 93,96 metros) foram classificados como sendo de baixa vulnerabilidade.

2.3.3. Presença de rios e / ou desembocaduras

A vulnerabilidade quanto a esse indicador foi analisada a partir da distância de cada setor aos rios presentes na praia. Nos setores em que havia a presença de rios ou desembocaduras, a distância considerada foi zero (alta vulnerabilidade). Nos outros setores, a vulnerabilidade foi considerada baixa se a distância dos rios era maior que 100 metros e considerada média se a distância era entre 50 e 100 metros. Para esse cálculo, foram traçados polígonos em volta dos rios que tivessem uma distância fixa dos mesmos (função buffer do ArcGIS 9.3), criando “áreas de influência” dos rios de 50 e 100 metros. A partir disso, pôde-se ter certeza do quão distante cada setor estava dos rios.

2.3.4.Elevação do terreno

Em virtude da inexistência de dados para a região costeira na praia de Massaguaçú, foram utilizados dados de altitude extraídos do projeto “TOPODATA - Banco de Dados Geomorfométricos do Brasil” (INPE, 2008). A partir das curvas de nível obtidas, criou-se uma rede triangular irregular (TIN - Triangular Irregular Network) para facilitar a análise altimétrica. O método de interpolação empregado foi a triangulação de Delaunay.

A TIN foi dividida em três intervalos de altura: menor que 3 metros (alta vulnerabilidade), entre 3 e 6 metros (média vulnerabilidade) e maior que 6 metros (baixa vulnerabilidade). A partir desse mapa, foram traçados polígonos ao redor das regiões de mesma altura, e suas áreas foram calculadas para se obter a vulnerabilidade predominante em cada setor.

Como as imagens aéreas antigas não possuem informações altimétricas e o período analisado é muito curto para se considerar variações expressivas da elevação do terreno, considerou-se a mesma altura para todas as fotografias aéreas analisadas.

2.3.5.Taxa de ocupação

A taxa de ocupação em porcentagem (T) foi obtida, para cada ano e setor analisados, através da Equação 1:

onde A_o representa a área de ocupação e A_u representa a área útil.

A área útil (A_u, representando aqui a área passível de ocupação) foi determinada a partir da Equação 2:

É importante ressaltar que a área útil dos setores foi diferente para cada ano, uma vez que o curso dos rios e a posição da linha de costa mudaram durante o período estudado.

2.3.6.Permeabilidade do solo

A permeabilidade do solo foi classificada de acordo com o resultado da taxa de ocupação, uma vez que ela é relativa à pavimentação e à urbanização do local analisado. Assim, nos setores que apresentaram alta taxa de ocupação o indicador foi considerado como de alta vulnerabilidade, para taxa média de ocupação o indicador foi considerado de média vulnerabilidade e para baixa taxa de ocupação, baixa vulnerabilidade.

2.3.7.Demais indicadores

Os indicadores configurações ao largo, vegetação e obras de engenharia costeira (vide Tabela 2) foram classificados a partir de dados da literatura e observação visual das fotografias aéreas.

3. RESULTADOS

3.1. Indicadores de vulnerabilidade

Os indicadores de vulnerabilidade foram analisados e classificados de acordo com sua respectiva vulnerabilidade (alta, média ou baixa) para cada ano analisado. A compilação dos resultados da análise da vulnerabilidade dos indicadores está apresentada na Tabela 3, e alguns dos resultados estão descritos em detalhes a seguir.

Através da função de MATLAB já explicitada no item 2.1.1 determinou-se se a linha de costa progradou ou retrogradou nos intervalos de tempo especificados. A rotina gerou gráficos com informações quantitativas da variação da linha de costa. Quando os resultados apresentaram um maior número de perfis em progradação, a vulnerabilidade correspondente ao setor foi baixa. Já nos setores em que os resultados eram de retrogradação em sua maioria, a vulnerabilidade foi alta. O ano de 1962 não apresentou resultados para esse indicador porque não havia linha de costa anterior para se comparar.

A Figura 2 apresenta os gráficos da variação da posição da LC entre os anos subsequentes e setores analisados. O setor norte apresentou retrogradação no período de 1962 a 1977 e progradação no período de 1977 a 2001. O setor centro apresentou progradação nos períodos de 1973 a 1987 e retrogradação nos períodos de 1962 a 1973 e 1987 a 2001. Já o setor sul apresentou retrogradação no período de 1962 a 1977 e progradação no período de 1987 a 2001. Analisando os quase 40 anos como um todo, ou seja, o intervalo entre 1962 e 2001, observa-se que houve progradação no setor norte e retrogradação nos setores centro e sul.

3.1.2.Largura da praia

Na medição da extensão praial ao longo da área de estudo, o maior comprimento obtido dentre todos os anos foi 93,96 metros, e o menor foi 8,78 metros. Com a divisão dos resultados em três intervalos iguais, ficou estabelecido que perfis praiais de 8,78 a 37,17 metros seriam considerados de alto risco, de 37,18 a 65,56 metros seriam considerados de médio risco e de 65,57 a 93,96 metros seriam considerados de baixo risco. Os resultados da classificação dos transectos quanto à vulnerabilidade estão apresentados na Tabela 4.

3.1.3.Presença de rios e / ou desembocaduras

O arco praial da área de estudo apresenta dois rios ao longo de seu comprimento (presentes em todas as imagens aéreas analisadas), o rio Bracuí no setor norte e o rio Capricórnio no setor sul. O setor centro dista mais de 100 metros desses rios.

Desse modo, a vulnerabilidade à erosão das imagens aéreas quanto a esse indicador foi classificada como alta nos setores norte e sul e como baixa no setor centro, em todos os anos analisados.

3.1.4.Elevação do terreno

A Tabela 5 apresenta as áreas correspondentes a cada tipo de vulnerabilidade, nos diferentes setores. A vulnerabilidade média prevaleceu em todos os setores analisados, representando 47,8%, 67,1% e 44,4% da área total dos setores norte, centro e sul, respectivamente. Dessa forma, a vulnerabilidade à erosão das imagens aéreas analisadas foi classificada como média em todos os setores e anos analisados em relação a esse indicador.

Ao se analisar a evolução da taxa de ocupação, percebe-se que a mesma cresceu em todos os setores analisados, especialmente no setor centro, que teve um crescimento praticamente exponencial a partir de 1977 (Figura 3).

A evolução da área ocupada pela urbanização fica ainda mais visível na Figura 4, onde também se pode observar o consequente aumento da impermeabilização do solo e diminuição da área ocupada pela vegetação. Na figura, a área representada como “continente” representa o solo exposto, aquele não ocupado nem por vegetação nem por ocupação antrópica.

3.1.6.Permeabilidade do solo

A permeabilidade do solo foi classificada como média no ano de 2001, para todos os setores, e como baixa nos demais anos e setores. Essa classificação seguiu a classificação da taxa de ocupação.

3.1.7.Configurações ao largo

No setor norte da praia encontram-se as Ilhas Cocanha, que protegem a costa fazendo com que as ondas incidentes percam energia antes de atingir a praia. Essa configuração caracteriza uma vulnerabilidade baixa.

Essa configuração foi observada para todos os anos analisados. Deste modo, a vulnerabilidade à erosão das imagens aéreas, quanto a esse indicador, foi classificada como baixa no setor norte e alta nos demais setores.

3.1.8.Vegetação

Através da análise visual das fotografias aéreas do período, observa-se que 1962 é o único ano em que há vegetação densa na região de estudo. A área coberta pela vegetação diminui nos demais anos, e observa-se pouca ou nenhuma vegetação no entorno da praia nos últimos anos analisados.

Desse modo, definiu-se a vulnerabilidade como sendo baixa em todos os setores do ano de 1962. Nos setores do ano de 1973 e nos setores norte e centro do ano de 1977 a vulnerabilidade foi classificada como média. Já para o setor sul do ano de 1977 e para todos os setores dos anos de 1987 e 2001, definiu-se a vulnerabilidade como sendo alta.

3.1.9.Obras de engenharia costeira

De acordo com Ceccarelli (2009), desde 2005, obras de engenharia têm sido feitas na praia de Massaguaçú. A presença dessas obras na área de estudo indica uma alta vulnerabilidade à erosão. Entretanto, em relação a esse indicador, o presente estudo não definiu a vulnerabilidade como alta. Isso porque o indicador foi analisado até o ano de 2001, e as obras só foram iniciadas em 2005. Sendo assim, houve ausência de estruturas de contenção costeira no período analisado (1962 - 2001), e a vulnerabilidade à erosão foi classificada como baixa em relação a esse indicador em todos os setores e anos analisados.

3.2.Integração dos resultados dos indicadores de vulnerabilidade

A análise de vulnerabilidade se deu através de porcentagens, utilizando como base a Tabela 3. Todos os indicadores tiveram o mesmo peso na análise, e a porcentagem se deu com relação ao número de indicadores analisados para cada ano (oito indicadores para 1962, já que para esse ano não houve análise de posição de linha de costa devido à falta de dados anteriores, e nove indicadores para os demais anos). A evolução das taxas de vulnerabilidades dos setores norte, centro e sul está explicitada na Figura 5.

No ano de 1962 houve predominância de vulnerabilidade baixa em todos os setores analisados (62,5% nos setores norte e centro e 50% no setor sul). No ano de 1973 houve predominância de vulnerabilidade baixa nos setores norte e centro (44,4% nos dois setores) e maior taxa de vulnerabilidade alta no setor sul (44,4%). Em 1977, os setores norte e centro apresentaram predominância de vulnerabilidade baixa (44,4% e 55,6%, respectivamente) e o setor sul teve maior porcentagem de vulnerabilidade alta (55,6%). O ano de 1987 apresentou predominância de vulnerabilidade baixa nos setores norte e centro (55,6%), e as mesmas porcentagens de vulnerabilidade alta e baixa foram observadas nos setores norte e sul (44.4%). Já o ano de 2001 apresentou maior porcentagem de vulnerabilidade média nos setores norte e sul (44,4%), e predominância de vulnerabilidade alta no setor centro (44,4%). Na análise do período completo (1962 – 2001), percebe-se que as vulnerabilidades alta e média aumentaram em todos os setores, sendo que a alta aumentou principalmente no setor centro (de 12,5% para 44,4%).

4. DISCUSSÃO

Analisando a evolução dos resultados para cada setor, percebe-se grande variação nas porcentagens da vulnerabilidade à erosão devido às oscilações da linha de costa e da extensão praial durante o período da avaliação. Nos setores norte e sul a vulnerabilidade aumentou até o ano de 1977, sofreu uma queda em 1987 e voltou a aumentar em 2001. No setor centro a vulnerabilidade também aumentou, mas teve queda nos anos de 1977 e 1987.

Dentre os nove indicadores empregados nesse estudo, quatro deles apresentaram resultados constantes ao longo do tempo (configurações ao largo, presença de rios e/ou desembocaduras, elevação do terreno e obras de engenharia costeira). Quanto aos outros cinco indicadores, três deles (vegetação, taxa de ocupação e permeabilidade do solo) estão intrinsecamente relacionados à urbanização da área de estudo, e apresentaram resultados proporcionais ao aumento da ocupação do local. Somente dois dos indicadores, posição da linha de costa e largura da praia, oscilaram ao longo do tempo, não apresentando um padrão regular de comportamento.

O aumento da ocupação humana teve impacto na região de estudo ao longo dos anos, mas não parece ser um fator determinante ao aumento da vulnerabilidade, uma vez que houve predominância de vulnerabilidades baixas mesmo em anos em que a ocupação esteve crescente. Entretanto, ao se analisar o período entre 1962 e 2001 percebe-se que a vulnerabilidade aumentou principalmente devido ao crescimento da taxa de ocupação e, consequentemente, dos demais indicadores relacionados a ela (vegetação e permeabilidade do solo). Os maiores índices de vulnerabilidade ocorreram especialmente no setor centro (a vulnerabilidade alta à erosão cresceu 31,9%), fato corroborado pelo intenso processo erosivo que tem afetado diretamente a rodovia presente neste setor. Com base nos estudos de Rogacheski (2010) e Luca (2011), acredita-se que a erosão severa em Massaguaçú tem suas causas advindas do vetor de transporte de deriva litorânea, uma vez que a corrente retira sedimentos da parte central da praia depositando-os nos limites norte e sul.

O método para a avaliação do ambiente utilizado no estudo é um procedimento rápido e de baixo custo para a avaliação do ambiente, tornando-o acessível para os gestores dos ambientes costeiros realizarem avaliações rápidas do meio. Conhecer as áreas mais vulneráveis à erosão auxilia os gestores nas tomadas de decisão e no melhor planejamento da orla costeira. Estudos como o proposto neste trabalho são importantes, principalmente no atual cenário de mudanças climáticas globais, uma vez que, com o provável aumento do nível médio do mar, as cidades litorâneas devem se planejar para mitigarem as consequências de alagamentos, intrusão de água salina e outros tantos problemas decorrentes da mudança do padrão climático.

A vulnerabilidade à erosão costeira apresentou um comportamento oscilatório devido às variações dos indicadores largura da praia e posição da linha de costa. Apesar disso, houve tendência de crescimento da vulnerabilidade durante o período analisado em razão dos indicadores taxa de ocupação, vegetação e permeabilidade do solo; esse aumento foi mais acentuado no setor central.

Apesar de a ocupação humana ser um influente fator na estrutura e no equilíbrio costeiro, a mesma não contribui efetivamente para o processo erosivo em Massaguaçú. A maior contribuição para o agravamento da erosão na área de estudo parece ser o comportamento da corrente de deriva litorânea, que transporta os sedimentos da parte central da praia para os setores vizinhos, como sugerido por Rogacheski (2010) e Luca (2011).

Para estudos futuros, sugere-se uma análise quantitativa da vulnerabilidade da região, para que seja possível estimar a magnitude das mudanças na costa ocasionadas pelo provável aumento do nível médio do mar.

AGRADECIMENTOS

Os autores agradecem à FAPESP (Fundação de Amparo à Pesquisa do Estado de São Paulo) e ao CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico) pelo suporte financeiro.

REFERÊNCIAS BIBLIOGRÁFICAS

Araujo, R.S.; Silva, G.V.; Freitas, D.; Klein, A.H.F. (2009) - Georreferenciamento de fotografias aéreas e análise da variação da linha de costa. In: J. Alcántara-Carrió, I.D. Correa, F. Ila, M. Alvarafo, A.H.F. Klein, J.A. Cabrera & R. Sandoval (eds.), Métodos en Teledetección Aplicada a la Prevención de Riesgos Naturales en el Litoral, pp.123-138, Cádiz, Eha. ISBN: 978-84-96023-67-3. Disponível em: http://www.rc.unesp.br/igce/planejamento/download/isabel/sist_inf_geogr_aula5/foto_aerea_cap8.pdf        [ Links ]

Boak, E.H.; Turner, I.L. (2005) - Shoreline Definition and Detection: A Review. Journal of Coastal Research (ISSN: 0749-0208), 21(4):688–703, West Palm Beach, FL, USA. Disponível em http://www6.hawaii.gov/dlnr/occl/files/Shoreline/JCR-VOL21-4.pdf

Bosom, E.; Jiménez, J.A. (2011) - Probabilistic coastal vulnerability assessment to storms at regional scale e application to Catalan beaches (NW Mediterranean). Natural Hazards and Earth System Science, 11:475-484. DOI: 10.5194/nhess-11-475-2011        [ Links ]

Bush, D.M.; Neal, W.J.; Young, R.S.; Pilkey, O.H. (1999) - Utilization of geoindicators for rapid assessment of coastal-hazard risk and mitigation. Ocean & Coastal Management, 42(8):647-670. DOI: 10.1016/S0964-5691(99)00027-7        [ Links ]

Cai, F.; Su, X.; Liu, J.; Li, B.; Lei, G. (2009) - Coastal erosion in China under the condition of global climate change and measures for its prevention. Progress in Natural Science, 19(4):415-426. DOI: 10.1016/j.pnsc.2008.05.034.         [ Links ]

Ceccarelli, T.S. (2009) - Paradigmas para os projetos de obras marítimas no contexto das mudanças climáticas. 125p., Dissertação de Mestrado, Universidade de São Paulo, São Paulo, SP, Brasil. Não publicado.         [ Links ]

Fatorić, S.; Chelleri, L. (2012) - Vulnerability to the effects of climate change and adaptation: The case of the ish Ebro Delta. Ocean & Coastal Management, 60: 1-10. DOI: 10.1016/j.ocecoaman.2011.12.015        [ Links ]

INPE (2008) - TOPODATA - Banco de Dados Geomorfométricos do Brasil (Portal internet). Instituto Nacional de Pesquisas Espaciais (INPE), São José dos Campos, SP, Brasil. In: http://www.dsr.inpe.br/topodata/index.php        [ Links ]

IPCC. (2007) - Summary for Policymakers. In: Palutikof, J., Van der Linden, P., Hanson, C. (eds.), Climate Change 2007: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, pp. 7–22, IPCC - Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge, U.K.

Júnior, J.B.T.; Ouverney, M.I.M.; Ubert, M.S.; Antunes, M.A.H. (2006) - Avaliações de Imagens Ikonos II e QuickBird para obtenção de Bases Cartográficas para cadastro técnico municipal. In: Congresso Brasileiro de Cadastro Técnico Multifinalitário. Anais. Florianópolis, SC, Brasil. Disponível em http://www.ufrrj.br/institutos/it/deng/marlene/downloads/Artigos/COBRAC%202006%202.pdf        [ Links ]

Linnekamp, F.; Koedam, A.; Baud, I.S.A. (2011) - Household vulnerability to climate change: Examining perceptions of households of flood risks in Georgetown and Paramaribo. Habitat International, 35(3):447-456. DOI: 10.1016/j.habitatint.2010.12.003        [ Links ]

Luca, C.B. (2011) - Implementação de ferramentas numéricas e bases de dados no SMC-Brasil e sua aplicação no estudo piloto da praia de Massaguaçu-Brasil. 213p., Dissertação de Mestrado, Universidade de Cantabria, Santander, Eha. Não publicado.         [ Links ]

Muehe, D. (2005) - Aspectos gerais da erosão costeira no brasil. Revista Mercator (ISSN: 1984-2201), 4(7):97-110, Fortaleza, CE, Brasil. Disponível em http://www.mercator.ufc.br/index.php/mercator/article/view/113/85        [ Links ]

Muehe, D. (org). (2006) - Erosão e progradação no litoral brasileiro. 476 p., Instituto Brasileiro do Meio Ambiente e dos Recursos Naturais Renováveis, Brasília, DF, Brasil. ISBN: 85-7738-028-9.         [ Links ]

Nicholls, R. (2004) - Coastal flooding and wetland loss in the 21st century: changes under the SRES climate and socio-economic scenarios. Global Environmental Change, 14(1):69-86. DOI: 10.1016/j.gloenvcha.2003.10.007        [ Links ]

Nicholls, R.J.; Cazenave, A. (2010) - Sea-level Rise and its impact on coastal zones. Science, 328(5985):1517-1520. DOI: 10.1126/science.1185782        [ Links ]

Pendleton, E.A.; Barras, J.A.; Williams, S.J.; Twichell, D.C. (2010) - Coastal Vulnerability Assessment of the Northern Gulf of Mexico to Sea-level rise and Coastal Change. 26p., U.S. Geological Survey, Open-File Report 2010-1146, Reston, Virginia, USA. http://pubs.usgs.gov/of/2010/1146/pdf/ofr2010-1146.pdf        [ Links ]

Pianca, C.; Mazzini, P.L.F.; Siegle, E. (2010) - Brazilian offshore wave climate based on NWW3 reanalysis. Brazilian Journal of Oceanography, 58(1):53-70. DOI: 10.1590/S1679-87592010000100006        [ Links ]

Rahmstorf, S. (2007) - A semi-empirical approach to projecting future sea-level rise. Science, 315(5810):368-370. DOI: 10.1126/science.1135456        [ Links ]

Rogacheski, C.E. (2010) - A Dinâmica Sedimentar e a Caracterização de Zonas de Erosão Acentuada (ZEA) ao Longo do Arco Praial de Massaguaçu, SP. 168p., Dissertação de Mestrado, Universidade de São Paulo, São Paulo, SP, Brasil. Não publicado.         [ Links ]

Sousa, P.H.G.O.; Siegle, E.; Tessler, M.G. (2013) - Vulnerability assessment of Massaguaçú Beach (SE Brazil). Ocean and Coastal Management, 77:23-40. DOI: 10.1016/j.ocecoaman.2012.03.003        [ Links ]

Sousa, P.H.G.O.; Siegle, E.; Tessler, M.G. (2011) - Environmental and Anthropogenic Indicators for Coastal Risk Assessment at Massaguaçú Beach (SP) Brazil. Journal of Coastal Research (ISSN: 0749-0208), SI64:319-323, Szczecin, Poland. Disponível em http://www.ics2011.pl/artic/SP64_319-323_P.H.G.O.Sousa.pdf        [ Links ]

Souza, C.R.G.; Luna, G.C. (2010) - Variação da linha de costa e balanço sedimentar de longo período em praias sob risco muito alto de erosão do município de Caraguatatuba (Litoral Norte de São Paulo, Brasil). Revista de Gestão Costeira Integrada, 10(2):179-199. DOI: 10.5894/rgci176        [ Links ]

*Submission: 29 October 2012; Evaluation: 4 December 2012; Accepted: 4 December 2012; Available on-line: 31 May 2013