DOI 10.14295/derb.v43.772 Derbyana, São Paulo, 43: e772, 2022.

POTENCIALIDADES DOS LIMIARES DE CHUVA INTENSIDADE-DURAÇÃO NA

PREVISÃO DE DESLIZAMENTOS EM JARAGUÁ DO SUL (SC)

Rodrigo Augusto STABILE

Rafael Alexandre Ferreira LUIZ

Andrezza Marques FERREIRA

Harideva Maturano EGAS

Márcio Roberto Magalhães de ANDRADE

RESUMO

Um dos grandes desafios dos sistemas de monitoramento e alerta de desastres

naturais é prever em que momento e sob quais condições os deslizamentos ocorrem.
Os limiares intensidade-duração, ao utilizar a duração e a intensidade média dos even-
tos de chuva para estimar quais deles são capazes de deflagrar deslizamentos, têm
potencial prático e têm sido extensivamente aplicados mundialmente na análise de
deslizamentos deflagrados pela chuva. Consequentemente, buscando ampliar as pos-
sibilidades de previsão de deslizamentos no Brasil e explorar a boa resolução tem-
poral e espacial das informações de chuva dos pluviômetros do Centro Nacional de
Monitoramento e Alertas de Desastres Naturais (Cemaden), esta pesquisa teve como
objetivo definir limiares de chuva do tipo intensidade-duração para o município de
Jaraguá do Sul (SC). Para isso, foram utilizados dez eventos de chuva com registro de
deslizamentos, definidos por chuvas com um intervalo de até 12 h com ausência de
precipitação, para os quais foram calculadas a intensidade (I), medida em mm/h, e a
duração (D), medida em horas. Para a construção dos limiares os dados de I e D dos
eventos de chuva com deslizamentos foram plotados em gráfico log-log e ajustados a
partir de uma Função Potência com o mesmo expoente da linha de regressão dos da-
dos. Foi possível elaborar dois limiares de chuva: o Limiar A, que abrange quaisquer
ocorrências de deslizamento e o Limiar B, que delimita os eventos de precipitação
capazes de deflagrar dez ou mais deslizamentos. Com taxas de acerto médias de até
29,8% e 62,5%, respectivamente, estimou-se que 3,2% dos eventos de chuva ocorri-
dos entre 2016 e maio/2022 estão acima do Limiar A e que 0,4% superam o Limiar
B. Os resultados encorajam a expansão de estudos de limiares intensidade-duração
para outros municípios e regiões de maneira a aperfeiçoar o potencial de previsão de
desastres causados por deslizamentos no território brasileiro.

Palavras-chave: Desastres naturais; Escorregamentos; Movimentos de massa; Eventos

extremos; Geografia física.

ABSTRACT

POTENTIALITIES OF INTENSITY-DURATION RAINFALL THRESHOLDS

IN THE FORECAST OF LANDSLIDES IN JARAGUÁ DO SUL (SC). One of
the great challenges of natural disaster monitoring and early warning systems is to
predict when and under what conditions landslides are triggered. Intensity-duration
rainfall thresholds use rainfall events’ duration and intensity to estimate which of them
are capable of triggering landslide and have practical applicability. They have been
extensively applied worldwide in analyzing landslides triggered by rain. Thus, seeking
to expand the possibilities of forecasting landslides in Brazil and taking advantage of

1
Stabile et al.

temporal and spatial resolution of rainfall data of rain gauges from Brazil’s National
Center for Monitoring and Early Warning of Natural Disasters (Cemaden), this paper
aimed to define intensity-duration rainfall thresholds for the municipality of Jaraguá
do Sul (Santa Catarina state, Brazil). Ten events with records of landslides were used
to calculate rainfall intensity — I (mm/h) and duration — D (h). A rainfall event was
defined by rains with an interval of up to 12 hours with no precipitation. Rainfall
thresholds were created plotting I and D data in a log-log graph and adjusted using
a Power Function with the same exponent as the regression line of the data. Two
rainfall thresholds were developed: A, which comprises any landslide occurrences,
and B, which embraces rainfall events capable of triggering ten or more landslides.
With average success rates of up to 29.8% and 62.5%, respectively, it was estimated
that 3.2% of the rainfall events that occurred between 2016 and May 2022 overcame
threshold A and that 0.4% overcame threshold B. Results encourage an increase in
research on intensity-duration rainfall thresholds in other Brazilian regions in order to
improve the capacity of predicting disasters caused by landslides.

Keywords: Natural disasters; Landslides; Mass movements; Extreme events; Physical

geography.

1 INTRODUÇÃO deslizamentos (ALEOTTI 2004, WIECZOREK &

GLADE 2005, GUZZETTI et al. 2007, 2008).
Um sistema de monitoramento para emissão Os limiares do tipo intensidade-duração (li-
de alertas de desastres naturais eficaz deve identifi- miares ID), que utilizam a duração e a intensida-
car as condições para a deflagração dos processos de “média” do evento de chuva responsável por
em tempo real. No caso dos deslizamentos condi- rupturas nas encostas, são os mais comuns na li-
cionados pela chuva, além de equipamentos para
teratura internacional (LARSEN & SIMON 1993,
mensuração e estimativa da chuva e modelos de
WIECZOREK & GLADE 2005, GODT et al.
previsão do tempo adequados, um dos grandes de-
2006, GUZZETTI et al. 2007), desde o trabalho
safios está em prever em que momento e sob quais
pioneiro de CAINE (1980). A principal justifica-
condições as rupturas ocorrem.
tiva geotécnica para o uso e funcionamento des-
Atualmente, os limiares empíricos (GU-
te tipo de limiar é que os deslizamentos (sobre-
ZZETTI et al. 2007, GUZZETTI et al. 2008, SAI-
tudo os escorregamentos rasos) são comumente
TO et al. 2010) são os modelos mais utilizados em
deflagrados por eventos de precipitação de forte
sistemas de monitoramento e alerta para prever
quais eventos de chuva têm potencial para defla- intensidade (GUIDICINI & IWASA 1977). Estão,
grar movimentos de massa, especialmente desli- portanto, vinculados aos fenômenos instantâneos
zamentos. Para o escopo desta pesquisa, o limiar que detonam os processos de deslizamento, como
pode ser definido como a quantidade mínima de linhas de fluxos que geram forças de percolação,
chuva para a possível ocorrência de um movimento erosão e saturação do solo em descontinuidades
de massa (deslizamento) em um determinado perí- hidráulicas (TATIZANA et al. 1987, FERNAN-
odo (REICHENBACH et al. 1998, GUZZETTI et DES et al. 2001).
al. 2007, SEGONI et al. 2018). No Brasil, porém, desconhecemos estudos
Diferentemente dos modelos em bases físicas, de limiares do tipo intensidade-duração para mo-
que acoplam modelos de estabilidade e hidrológi- vimentos de massa. Em geral, os limiares de chuva
cos para determinar a ocorrência de deslizamentos definidos no território brasileiro utilizam acumu-
e, por isso, necessitam de complexos parâmetros lados de chuva de 24, 48 ou 72 h (GUIDICINI &
de entrada, os limiares empíricos são relativamente IWASA 1977, SANTORO et al. 2010, MOLINA
simples de serem implementados. Para sua cons- et al. 2015) ou associações entre a intensidade da
trução são necessárias informações sobre eventos chuva (algumas vezes apenas com dados de 24 h)
de deslizamento pretéritos e os dados de chuva. Es- e a chuva antecedente (TATIZANA et al. 1987,
tes últimos podem ser empregados das mais diver- D’ORSI 2011, D’ORSI et al. 2012, MENDES et
sas maneiras, fazendo com que diferentes tipos de al. 2015, METODIEV et al. 2018). De fato, com
limiares de chuva sejam utilizados na previsão dos exceções muito pontuais, o Brasil só passou a ter

2
 Derbyana, São Paulo, 43: e772, 2022.

um registro sistemático da chuva em escala tempo- ção cálcio-alcalina geralmente básica, com porções
ral adequada para este tipo de análise (ex. a cada restritas de formações ferríferas, paragnaisses e ou-
10, 15 ou 60 min) a partir da instalação da rede tros tipos de rochas, e por Sedimentos Quaterná-
de pluviômetros automáticos do Centro Nacional rios, compostos por sedimentos aluvionares-cas-
de Monitoramento e Alertas de Desastres Naturais calheiras, areias e sedimentos siltico-argilosos. Na
(Cemaden), iniciada em 2011. porção norte do município, destacam-se os depósi-
Em vista disso, buscando ampliar as tos vulcânico-sedimentares do Grupo Campo Ale-
possibilidades de previsão de deslizamentos no gre e os maciços graníticos da Suíte Intrusiva da
Brasil e explorar o potencial das informações de Serra do Mar (TABALIPA 2019). Acerca da pedo-
chuva dos pluviômetros do Cemaden, esta pesquisa logia, verifica-se a predominância de cambissolos
teve como objetivo definir limiares de chuva do do tipo álico, distrófico e húmico álico, bem como
tipo intensidade-duração no município de Jaraguá de solos podzólicos vermelho-amarelo, latossóli-
do Sul (SC). cos álicos e podzólicos vermelho-amarelo álicos.
A oeste do município identifica-se menor porção
2 ÁREA DE ESTUDO de solos do tipo litólicos álicos (TABALIPA 2019).
Em relação aos aspectos geomorfológicos,
O município de Jaraguá do Sul (Figura 1), o município está localizado em área de planície
situado no nordeste de Santa Catarina, possui ex- formada na confluência de dois rios (Itapocu e Ja-
tensão territorial de 530,894 km2 e população esti- raguá), cercado por serras cristalinas subparalelas
mada em 184.579 habitantes (IBGE 2022). entre si, com disposição predominante NE-SW,
O substrato do município é formado prin- que se iniciam ao sul da Serra do Mar e do Pla-
cipalmente pelo Complexo Luiz Alves, que inclui nalto de São Bento do Sul, margeando todo lado
gnaisses granulíticos ortoderivados, de composi- leste do estado de Santa Catarina até a altura do

FIGURA 1 – Localização da área de estudo e setorização das áreas de risco. Adaptado de CPRM (2011, 2019).
FIGURE 1 – Location of the study area and sectorization of the risk areas. Adapted from CPRM (2011, 2019).

3
Stabile et al.

município de Jaguaruna, no sul catarinense. A Geológico do Brasil – CPRM em dois momentos

geomorfologia regional é caracterizada por vales distintos: o primeiro realizado em 2011, e o se-
profundos e encostas íngremes, o que, quando as- gundo realizado em 2018, como complemento ao
sociado à precipitação intensa, favorece a atuação primeiro. Ao analisar a tabela, verifica-se que o
de processos erosivos e inclusive a ocorrência de município possui 3.184 edificações (nas quais vi-
movimentos de massa (BERTOLI 2006, SANTA vem 13.226 pessoas), instaladas em 19 setores clas-
CATARINA 2016). sificados como risco muito alto, e 17 classificados
Segundo a classificação climática de Köppen, como risco alto, todos relacionados a pelo menos
o clima de Jaraguá do Sul faz parte do grupo C um ou mais processos em que ocorrem movimentos
(mesotérmico), úmido (f) e está incluído no sub- de massa (CPRM 2011, 2019).
tipo Cfa, que caracteriza a região como de clima As áreas mapeadas estão localizadas geral-
subtropical, com temperaturas médias no mês mais mente na porção basal e média das encostas, locais
frio inferior a 18 ºC e no mês mais quente acima para os quais a urbanização vem avançando, por
de 22 ºC, com verões quentes, pouca frequência de vezes de maneira desordenada (Figura 1). Observa-
geadas e tendência de concentração das chuvas nos -se ainda que o principal fator de risco no município
meses de verão, porém sem estação seca definida corresponde à ocupação de terrenos desfavoráveis,
(PANDOLFO et al. 2002). De acordo com a figura estando os setores de risco localizados em encostas
2, janeiro é o mês mais chuvoso, com média de 261 íngremes associados a taludes de corte, onde geral-
mm de precipitação, e o mês menos chuvoso é ju- mente é constatado elevado grau de umidade. De
nho, com média de 102 mm. maneira geral, o substrato nestes setores é composto
São variados os sistemas atmosféricos que por uma camada fina de solo sobre maciço rochoso
atuam na região e influenciam desastres relacio- (colúvio ou tálus), que podem conter fragmentos de
nados a deslizamentos, destacando-se: os Sistemas rocha de tamanhos diversos. (CPRM 2011, 2019).
Frontais, que geram instabilidades e produzem
muita chuva resultante do encontro de duas mas- 2.2 Histórico de desastres
sas de ar distintas; os Ciclones extratropicais, cuja
Devido sua localização, o município de Ja-
convergência de massas de ar devido a sistemas de
raguá do Sul apresenta suscetibilidade a processos
baixa pressão comuns na costa atlântica catarinense
naturais de denudação que, dadas as condições cli-
se propagam junto com frentes polares; os Sistemas
máticas, geológicas, geomorfológicas e de ocupação
Convectivos, comuns no verão, com a formação de
antrópica, aumentam a vulnerabilidade e o risco de
nuvens profundas; e, os Complexos Convectivos de
inundações e movimentos de massa. Os maiores
Mesoescala, com intensidade suficiente para gerar
eventos que geraram desastres ocorreram em 2008,
chuvas fortes (KOBIYAMA et al. 2006).
2011 e 2014 (SILVA & PINTO 2013, TABALIPA
2.1 Setores de risco 2019).
Embora os eventos de 2011 e 2014 tenham
A tabela 1 apresenta uma síntese dos setores resultado em severos danos à população, sobretudo
de risco mapeados em Jaraguá do Sul pelo Serviço devido aos processos hidrológicos registrados (en-
xurradas e inundações), os impactos mais significa-
tivos associados a movimentos de massa ocorreram
nos eventos de 2008 (SILVA & PINTO 2013, TA-
BALIPA 2019). GOUVEIA et al. (2018) ressaltam
que as condições meteorológicas anômalas obser-
vadas no último trimestre de 2008, associadas a um
sistema de bloqueio atmosférico, contribuíram para
que os maiores volumes de chuva no Sul do país
tenham se concentrado no estado de Santa Catari-
na, principalmente na região litorânea e no Vale do
FIGURA 2 – Média da precipitação mensal de Jara- Itajaí, ultrapassando os 1.000 mm de chuva somen-
guá do Sul/SC. Fonte: Atlas Pluviométrico do Brasil te no mês de novembro em muitos municípios da
(PINTO et al. 2011).
região.
O município de Jaraguá do Sul declarou situ-
FIGURE 2 – Average monthly precipitation in Jaraguá do Sul/ ação de emergência em 20 de novembro de 2008 e,
SC. Source: Pluviometric Atlas of Brazil (PINTO et al. 2011). posteriormente Estado de Calamidade Pública em

4
 Derbyana, São Paulo, 43: e772, 2022.

TABELA 1 – Síntese dos setores de risco em Jaraguá do Sul, SC. Adaptado de CPRM (2011, 2019).
TABLE 1 – Synthesis of the risk sectors in Jaraguá do Sul, SC. Adapted from CPRM (2011, 2019).

N.º do Setor Tipologia* Edificações Pessoas Grau de Risco Ano

1 D, C 101 404 Muito Alto 2011
2 D, C 98 320 Muito Alto 2011
3 D, C, I 220 480 Muito Alto 2011
4 D, C, I 162 1198 Muito Alto 2011
5 D, C 246 972 Muito Alto 2011
6 D, C, I 67 268 Muito Alto 2011
7 D, I 144 576 Muito Alto 2011
8 D 25 200 Muito Alto 2011
9 D 46 184 Muito Alto 2011
10 D, Q, C, Ex 195 780 Muito Alto 2011
11 D, I 258 1032 Muito Alto 2011
12 D, I 636 2758 Muito Alto 2011
13 D, Q 61 324 Muito Alto 2011
14 D, Q, C, E 137 548 Muito Alto 2011
15 D, I 177 708 Muito Alto 2011
16 D, E 107 428 Muito Alto 2011
17 D, C, E 55 250 Muito Alto 2011
18 D, E 167 668 Muito Alto 2011
19 D 211 844 Muito Alto 2011
20 D, Ex 11 44 Alto 2018
21 C 5 20 Alto 2018
22 C 6 24 Alto 2018
23 D, Er 5 20 Alto 2018
24 Q 1 4 Alto 2018
25 D 3 12 Alto 2018
26 D, Q 6 24 Alto 2018
27 D, Q 1 4 Alto 2018
28 D 1 4 Alto 2018
29 D, I 4 16 Alto 2018
30 D 3 12 Alto 2018
31 D 2 8 Alto 2018
32 D, Er 1 4 Alto 2018
33 C 2 8 Alto 2018
34 D 2 8 Alto 2018
35 C 3 12 Alto 2018
36 C 15 60 Alto 2018
*D = Deslizamento; C = Corrida de massa, I = Inundação; Q = Queda de bloco;
Ex = Enxurrada; Er = Erosão.

14 de dezembro do mesmo ano. Segundo a defesa áreas não legalizadas pela administração municipal.
civil municipal, entre os meses de novembro e de- Enfatiza ainda que a conformação do relevo de Ja-
zembro de 2008 foram registradas 13 vítimas fatais raguá do Sul, com declive acentuado, determina que
(12 em decorrência de deslizamentos), 147 feridos e apenas 42% de sua área total não seja constituída de
40 famílias desalojadas devido a sucessivos desliza- encostas e que somente 21,3 % de sua área urbana,
mentos, corridas de massa e inundações decorrentes até meados da década de 90, fossem considerados
da elevada e contínua precipitação (SILVA & PIN- urbanizáveis. Com o escasso espaço superficial pla-
TO 2013). no e livre de inundações, o município possui forte
Sobre este aspecto, SILVA (2014) chama crescimento populacional e consequente demanda
atenção para a insegurança nas áreas de instabilida- por moradia, o que tem levado à ocupação das en-
de natural do município, agravadas pela ocupação costas tanto pelas classes de baixa como de alta ren-
irregular de encostas de morros, fundos de vale e da (BUTZKE 1995 apud SILVA 2014).

5
Stabile et al.

3 MATERIAIS E MÉTODOS descartados por motivos diversos, entre os quais: pro-

cesso predominante ser diferente de deslizamento,
3.1 Análise dos deslizamentos a chuva não constituir o principal fator deflagrador
(ex. deslizamentos em obras ativas com alterações
Os deslizamentos foram levantados no sistema significativas do equilíbrio da encosta), ausência de
de Registros de Eventos de Inundação e Desliza- pluviômetro automático num raio menor que 5 km no
mentos do Cemaden (REINDESC), compreenden- caso de eventos com um único deslizamento. O raio
do eventos de deslizamento ocorridos entre 2016 e de 5 km foi definido empiricamente considerando a
maio/2022. Esta base de dados é fomentada pelos variabilidade de um mesmo evento de chuva dentre
especialistas do Cemaden a partir de informações dis- os diversos pluviômetros utilizados na análise. No
ponibilizadas pelas Defesas Civis, outros sistemas de caso de deslizamentos múltiplos sem identificação da
registro nacionais, assim como por notícias da mídia localização, os eventos foram considerados, uma vez
em geral. Portanto, a qualidade dos dados – locali- que seria altamente improvável que ao menos um de-
zação, data e horário, caracterização geotécnica/hi- les não estivesse próximo a um pluviômetro.
drológica e classificação do processo etc. – apresenta Os dez eventos de deslizamento utilizados no
variações significativas. Dentre os catorze registros estudo (Figura 3 e Tabela 2) constituem, de modo
para o município de Jaraguá do Sul, quatro foram geral, deslizamentos rasos em áreas urbanas e pe-

FIGURA 3 – Localização dos eventos com deslizamentos. A classificação dos eventos quanto à precisão es-
pacial segue os seguintes critérios: Ausente - informação fonte faz referência à ocorrência em município sem
detalhamento da localização; Baixa - informação fonte refere-se à ocorrência em nível de bairro ou região; Mé-
dia - informação fonte refere-se à ocorrência em nível de localização ou rua; Alta - informação fonte identifica
local específico com possibilidade de identificação de coordenada geográfica. Elaborado pelos autores.
FIGURE 3 – Location of events with landslides. The classification of events in terms of spatial precision follows the criteria: Absent
- source information refers to the occurrence in the municipality without detailing the location; Low - source information refers to
the occurrence at the neighborhood or region level; Average - source information refers to occurrence at location or street level; High
- source information identifies specific location with possibility of geographic coordinate identification. Prepared by the authors.

6
 Derbyana, São Paulo, 43: e772, 2022.

riurbanas (Figura 4). A maior parte deles têm infor- cia (madrugada/manhã/tarde/noite). Cabe ressaltar
mações espaciais com precisão de bairro ou região que um dos grandes desafios para a construção de
e precisão temporal apenas do período da ocorrên- limiares confiáveis é a escassez de informações

TABELA 2 – Eventos com deslizamento em Jaraguá do Sul entre Janeiro/2016 e Maio/2022 utilizados para a
elaboração dos limiares.
TABLE 2 – Rainfall events with landslides in Jaraguá do Sul between January/2016 and May/2022 used for the elaboration of the
thresholds.

Código Ti* Classificação Pluviômetro ΔS (km)* I (mm/h)* D (h)*

01 20/02/2016 18:30 1 deslizamento 420890607A 1 14,11 4,17

02 22/12/2018 17:40 1 deslizamento 420890603A 3 5,97 12,00

03 26/05/2017 23:10 2 a 9 deslizamentos 420890605A 1 2,48 39,67

04 19/01/2019 16:00 2 a 9 deslizamentos 420890605A 2,5 8,61 8,50

05 16/02/2019 15:50 2 a 9 deslizamentos 420890602A 3 4,78 42,5

06 30/05/2019 03:20 2 a 9 deslizamentos 420890605A 1,5 2,17 40,83

07 12/12/2020 20:50 2 a 9 deslizamentos 420890603A N/A 25,53 1,83

08 03/11/2019 19:40 10 ou mais deslizamentos 420890606A N/A 29,18 3,67

09 18/01/2021 15:10 10 ou mais deslizamentos 420890607A 2 1,88 105,33

10 21/02/2022 18:40 10 ou mais deslizamentos 420890605A N/A 41,08 2,17

* Ti = início do evento de chuva em que ocorreu o(s) deslizamento(s); ΔS = distância aproximada entre o evento de desliza-
mento e o pluviômetro (quando não foi possível identificar a localização do deslizamento = N/A); I = intensidade média do
evento de chuva (mm/h); D = duração do evento de chuva (h). A definição de evento de chuva e o método para escolha do
pluviômetro de referência de cada evento estão no texto.

FIGURA 4 – Deslizamentos registrados em Jaraguá do Sul nos eventos 09 (a) e 05 (b), conforme códigos
informados na tabela 1. Fonte: (a) Prefeitura do Município de Jaraguá do Sul, disponível em https://scc10.
com.br/cotidiano/tempo/tres-imoveis-afetados-por-deslizamentos-em-jaragua-do-sul-passam-por-avaliacao-
-da-defesa-civil/, acesso em: 02/09/2022; (b) Prefeitura do Município de Jaraguá do Sul, disponível em <ht-
tps://g1.globo.com/sc/santa-catarina/verao/2019/noticia/2019/02/18/chuva-em-sc-causa-deslizamentos-deixa-
-familias-desabrigadas-e-vias-interditadas-nesta-segunda-feira.ghtml>, acesso em: 02/09/2022.
FIGURE 4 – Landslides recorded in Jaraguá do Sul in events 09 (a) and 05 (b), according to the codes reported in table 1. Source:
(a) Municipality of Jaraguá do Sul, available at https://scc10.com.br/cotidiano/tempo/tres-imoveis-afetados-por-deslizamentos-em-
-jaragua-do-sul-passam-por-avaliacao-da-defesa-civil/, accessed on: 02/09/2022; (b) Municipality of Jaraguá do Sul, available at
https://g1.globo.com/sc/santa-catarina/verao/2019/noticia/2019/02/18/chuva-em-sc-causa-deslizamentos-deixa-familias-desabriga-
das-e-vias-interditadas-nesta-segunda-feira.ghtml, accessed on: 02/09/2022.

7
Stabile et al.

satisfatórias (ALEOTTI 2004, GUZZETTI et al. Um dos grandes problemas metodológicos na

2008), levando em conta ainda que as informações construção e uso dos limiares ID consiste na defi-
da Defesa Civil e da mídia provavelmente só vão nição temporal do evento de chuva deflagrador. A
relatar ocorrências com algum tipo de impacto à escolha deste parâmetro interfere diretamente nos
população e à infraestrutura e não deslizamentos valores de intensidade (I) e duração (D) dos even-
isolados, ou cujo impacto tenha sido imperceptível tos de chuva analisados. GODT et al. (2006), por
àquela comunidade. exemplo, estabeleceram um período com ausência
de chuva de pelo menos 3 h para limitar diferen-
3.2 Coleta e tratamento dos dados de chuva
tes eventos. SAITO et al. (2010) e MATURIDI et
O Cemaden tem sete pluviômetros automá- al. (2021), por outro lado, separaram os eventos de
ticos no município de Jaraguá do Sul (Figura 5), chuva a partir de um período de ausência de chuva
distribuídos de 2 a 6 km entre si e que registram a de 24 h. No primeiro caso, os autores enfatizaram
precipitação a cada 10 min. Por conseguinte, para a intensidade da chuva como parâmetro de análise
cada um dos eventos de deslizamento listados foi do limiar ID, mas o associaram a um balanço hí-
atribuído um pluviômetro para ser utilizado como drico para explicar o papel da chuva antecedente
referência nos cálculos do evento de chuva defla- nos deslizamentos. Nos últimos, percebeu-se que
grador (Tabela 2). Como regra geral foi utilizado nas chuvas de longa duração apenas a intensidade
o pluviômetro mais próximo do evento quando média do evento de chuva é considerada, muitas
sua localização era conhecida. Nos casos de des- vezes desprezando o papel de picos de chuva in-
lizamentos múltiplos com localização conhecida tensos na deflagração dos deslizamentos (SAITO
foi escolhido um pluviômetro em um ponto médio et al. 2010), em conclusão similar à de GUZZET-
entre as ocorrências e nos casos de deslizamentos TI et al. (2007), que apontam que a intensidade da
múltiplos com localização imprecisa foi seleciona- chuva mensurada em períodos curtos ou longos
do o pluviômetro com as maiores intensidades de tem um significado físico diferente. Neste estudo,
chuva registradas. após a análise do comportamento de alguns even-

FIGURA 5 – Localização dos pluviômetros do Cemaden. Elaborado pelos autores.

FIGURE 5 – Location of Cemaden rain gauges. Prepared by the authors.

8
 Derbyana, São Paulo, 43: e772, 2022.

tos de chuva e levando em conta os estudos citados, Cabe notar que mesmo nos casos em que o
decidiu-se definir o evento de chuva a partir da au- horário aproximado da ocorrência de deslizamen-
sência de chuva por um período de pelo menos 12 to era conhecida, sendo estes exceções à precisão
h (antes ou depois do evento). temporal geral, foram utilizados os valores de I
A intensidade (I), medida em mm/h, e a du- e D do evento de chuva como um todo, ou seja,
ração (D), medida em horas, foram então estabele- mesmo após a virtual ocorrência do deslizamen-
cidas para cada um dos eventos de chuva com des- to. O uso da intensidade média e da duração do
lizamento. A figura 6 ilustra três eventos de chuva evento de chuva na sua totalidade é o método mais
utilizados neste estudo. A partir deles fica nítida a difundido na aplicação de limiares intensidade-
problemática da definição do evento de chuva: nas -duração (GUZZETTI et al. 2007, 2008), e dife-
figuras 6b e 6c um mesmo evento tem vários picos rem de métodos que consideram o horário preciso
de intensidade com períodos de ausência de chuva dos deslizamentos para calcular a chuva deflagra-
ou chuva fraca, de modo que uma especificação di- dora (ex. ALEOTTI 2004). Há duas justificativas
versa para o evento de chuva resultaria em valores principais para esta escolha. A primeira baseia-se
de I e D diferentes. na baixa precisão dos registros de deslizamento,

FIGURA 6 – Gráficos representando três eventos de chuva, códigos 10 (a), 09 (b) e 05 (c) conforme a tabela 2.
Cada barra azul representa a intensidade da chuva (mm/h) em um período de 10 min tal como mensurado pelo
pluviômetro automático.
FIGURE 6 – Graphs representing three rain events, codes 10 (a), 09 (b) and 05 (c) as shown in table 2. Each blue bar represents the
intensity of rain (mm/h) over a 10-minute period as measured by the automatic rain gauge.

9
Stabile et al.

de modo que o uso do horário do deslizamento deslizamento com o mesmo expoente β da linha de
na análise seria impossível na maioria dos casos regressão. Foi possível elaborar dois limiares de
e provavelmente adicionaria erros aos registros de chuva para o município de Jaraguá do Sul. O Li-
intensidade e duração da chuva. A segunda advém miar A (Equação 2) abrange quaisquer ocorrências
da possibilidade de comparação com os eventos de de deslizamento, enquanto o Limiar B (Equação
chuva que não geram deslizamentos. 3) delimita os eventos de precipitação capazes de
deflagrar dez ou mais deslizamentos. Nota-se que
não há diferença visível da relação entre I e D en-
4 LIMIARES INTENSIDADE-DURAÇÃO
tre aqueles eventos com apenas um deslizamento
Os limiares Intensidade-Duração (ID) são e aqueles que compreendem de dois a nove desli-
normalmente representados por uma Função Po- zamentos.
tência na seguinte forma (Equação 1): 𝐼 = 36,456𝐷 −0,76 (1,83 < 𝐷 < 105,3) (Equação 2)
𝐼 =∝ 𝐷 𝛽 (Equação 1) 𝐼 = 64,816𝐷 −0,76 (1,83 < 𝐷 < 105,3) (Equação 3)
Onde: I é a intensidade (média) da chuva Estes limiares delimitam eventos com
(mm/h), D é a duração da chuva (h), enquanto α e 1,83<D<105,33 e I variando entre 1,88 e 41,08
β são parâmetros da equação. mm/h. O Limiar B, particularmente, foi sistemati-
Para a construção dos limiares de Jaraguá do zado a partir de apenas três eventos com durações
Sul os dados de I e D dos eventos de chuva com curtas ou muito longas, sem dados intermediários,
deslizamentos (Tabela 2) foram plotados em um devendo ser observado com cautela.
gráfico log-log. A distribuição dos pontos mostrou Alguns estudos (GUZZETTI et al. 2007,
uma clara correlação negativa entre I e D (Figura SAITO et al. 2010) têm apontado que os mecanis-
7), levando à interpretação de que chuvas mais in- mos para deflagração de deslizamentos nas chuvas
tensas necessitam de menos tempo para deflagrar de forte intensidade e curta duração provavelmente
deslizamentos. não são exatamente os mesmos das chuvas mais
Após o ajuste da linha de regressão dos dados longas e de menor intensidade. Por outro lado, os
na forma de uma função potência, o limiar foi esta- dados da figura 6 sugerem que mesmo em chuvas
belecido no limite inferior de todos os eventos com de longa duração há picos de alta intensidade de
chuva que podem estar também associados à defla-
gração de deslizamentos, denotando que as ocor-
rências precisam de um gatilho, mesmo com altos
volumes de chuva antecedente.
Efetivamente, mesmo que os deslizamentos
sejam deflagrados pela água que imediatamente in-
filtra na encosta e gera forças de percolação, erosão
e saturação do solo (TATIZANA et al. 1987, FER-
NANDES et al. 2001), sabe-se que a precipitação
antecedente atua na redução da resistência ao cisa-
lhamento (GUIDICINI & IWASA 1977, TATIZA-
NA et al. 1987) e favorece a percolação da água no
solo devido à condutividade hidráulica na condição
saturada (GUIDICINI & IWASA 1977, WOLLE &
CARVALHO 1989, ALEOTTI 2004, GODT et al.
2006).

4.1 Parâmetros α e β e comparação com outros li-

FIGURA 7 – Gráfico log-log com a intensidade (eixo miares
vertical) e duração (eixo horizontal) dos eventos de Os limiares estabelecidos neste estudo pare-
chuva com deslizamento(s) e os limiares A e B (méto- cem ter parâmetro α relativamente alto (α=36,456
do de elaboração no texto). para o Limiar A e α=64,816 para o Limiar B) em
FIGURE 7 – Log-log graph with intensity (vertical axis) and comparação com outros estudos (Figura 8). Ain-
duration (horizontal axis) of rainfall events with landslide(s) da que não seja totalmente adequado interpretá-
and thresholds A and B (method of elaboration in the text). -lo isoladamente (sem β), dentre os limiares apre-

10
 Derbyana, São Paulo, 43: e772, 2022.

Acreditamos que as explicações para este resultado

inusitado baseiem-se em diferenças metodológicas
na definição do evento de chuva, mas, principal-
mente na forma de coleta dos dados, uma vez que o
autor utilizou dados de outros trabalhos publicados
internacionalmente que provavelmente relatavam
eventos de grande magnitude relativa.
O expoente β=-0,76 nesta pesquisa asseme-
lha-se ao de outros estudos realizados em áreas
tropicais (β=-0,82 em Porto Rico e β=-0,722 na
Malásia Peninsular), bem como a estudos regionais
que utilizaram eventos com períodos de tempo am-
plos (β=-0,7 para o sul e centro europeu). Valores
maiores de β: -0,39 em CAINE (1980), -0,44 em
GUZZETTI et al. (2008) e -0,26 em SAITO et al.
(2010), foram encontrados em estudos globais ou
em áreas muito amplas e com grande variação am-
biental, onde os dados utilizados para a construção
dos limiares normalmente têm uma maior disper-
são. Valores muito menores: β=-1,13 em GODT et
al. (2006) podem estar associados aos eventos com
amplitude de duração restrita utilizados por estes
autores. Outra explicação para estas variações po-
deria estar na própria escolha do método de defini-
FIGURA 8 – Comparação dos limiares A e B com ção das curvas (neste estudo utilizamos o mesmo
outros limiares intensidade-duração publicados. expoente da linha de regressão), mas as pesquisas
de GUZZETTI et al. (2007, 2008) empregaram a
FIGURE 8 – Comparison of thresholds A and B with other pu-
mesma técnica de estatística Bayesiana e encontra-
blished intensity-duration thresholds.
ram valores bastante diferentes de β, refutando, a
princípio, esta suposição.
sentados na figura 8, apenas os limiares definidos Assim, β ~ 0,7 ou 0,8 parece ser um valor
por LARSEN & SIMON (1993) para Porto Rico típico em limiares de chuva intensidade-duração
(α=91,46) e GODT et al. (2006) no NW dos Esta- para deslizamentos, quando o horizonte temporal
dos Unidos (α=82,73) apresentam um parâmetro α dos eventos é suficientemente extenso e quando
maior. Cabe ressaltar que na pesquisa de LARSEN são utilizados dados de uma região com relativa
& SIMON (1993), os autores utilizaram eventos homogeneidade geográfica, mas esta conjectura só
com dezenas a centenas de escorregamentos, en- pode ser confirmada a partir da análise de outros
quanto GODT et al. (2006) utilizaram eventos com limiares no Brasil e no mundo.
um período de duração restrito (entre 23 e 46 h).
Em geral, percebe-se que os limiares globais
5 TAXAS DE ACERTO
ou definidos para áreas muito grandes (ex. SAITO
et al. 2010, cujo limiar abrange todo o território Considerando que os limiares foram defini-
japonês) têm um parâmetro α bastante menor. O li- dos apenas a partir de eventos com deslizamentos e
miar global definido por GUZZETTI et al. (2008), objetivando seu uso prático (por exemplo, na emis-
por exemplo, teria pouco efeito prático na previsão são de alertas de desastres naturais), interessou-nos
de deslizamentos em Jaraguá do Sul, uma vez que saber qual seria a sua representatividade em rela-
ele abrangeria boa parte dos eventos de chuva re- ção aos eventos de chuva registrados em Jaraguá
gistrados nos pluviômetros. do Sul. Para fins metodológicos, todos os eventos
Neste sentido, chama a atenção o limiar glo- de chuva em que não constam registros de desli-
bal inicialmente proposto por CAINE (1980), onde zamentos no sistema de Registros de Eventos de
α=14,82, porém, β=-0,39, de modo que para chu- Inundação e Deslizamentos do Cemaden (REIN-
vas de maior duração (D>12h), seu limiar acaba DESC) foram considerados eventos com ausência
sendo mais restrito do que aqueles aqui propostos. de deslizamentos.

11
Stabile et al.

Calculamos os eventos de chuva acima dos para o Limiar A e entre 0% e 100% (=35%) para
limiares estabelecidos para cada um dos pluviôme- o Limiar B.
tros de Jaraguá do Sul (Tabela 3) e, em seguida, ava- Os pluviômetros 04A e 06A (Figura 5), situ-
liamos uma taxa de acerto, conforme a equação 4. ados mais a leste do município, apresentam taxas
𝐸𝑙𝑖𝑚𝑖𝑎𝑟;𝑑𝑒𝑠𝑙𝑖𝑧𝑎𝑚𝑒𝑛𝑡𝑜 de acerto notadamente mais baixas. Acreditamos
𝑇𝐴 (%) = × 100 (Equação 4) que isto ocorra por estarem mais distantes da
𝐸𝑙𝑖𝑚𝑖𝑎𝑟
maioria dos deslizamentos registrados no muni-
Onde: TA (%) é a taxa de acerto do limiar na- cípio, bem como por representarem muitas vezes
quele pluviômetro; Elimiar é o total de eventos acima um momento diferente da chuva daqueles pluviô-
do limiar registrados no pluviômetro; Elimiar;deslizamento metros situados mais a montante, entre as serras
é o total de eventos acima do limiar no pluviômetro e morros do vale do rio Itapocu. Já o pluviôme-
e com registro de deslizamentos (ou da quantidade tro 02A tem uma taxa de acerto particularmente
de deslizamentos prevista pelo limiar). baixa para o Limiar B. A explicação sustenta-se
Cabe notar que, uma vez que no nosso estu- no fato deste pluviômetro ter seu funcionamento
do o cálculo da intensidade (I) e da duração (D) suspenso em meados de 2019, antes que dois dos
dos eventos com deslizamento utilizou um pluvi- três registros de eventos com dez ou mais desli-
ômetro de referência para cada um dos eventos, é zamentos em Jaraguá do Sul tivessem ocorrido.
possível que um evento com deslizamento esteja Efetivamente, quando excluímos os pluviômetros
abaixo do limiar em pluviômetros diferentes da- 04A e 06A dos cálculos de taxa de acerto do Li-
quele, principalmente quando eles estão mais dis- miar A e excluímos os pluviômetros 02A, 04A e
tantes (Figura 9). Assim, o mesmo evento de chuva
06A dos cálculos de taxa de acerto do Limiar B,
que causou o(s) deslizamento(s) pode ter valores
as taxas médias aumentam para 29,8% no Limiar
de precipitação menores ou estar registrado em ou-
A e 62,5% no Limiar B.
tro horário em pluviômetros diferentes, devido ao
deslocamento e à dinâmica da tempestade. 5.1 Considerações sobre a intensidade e os acumu-
Sobre a frequência da ocorrência de eventos lados prévios de chuva
acima dos limiares (Tabela 3), nota-se que eventos
acima do Limiar A representam entre 2,6% e 3,6% A observação da figura 9 nos permite veri-
do total de eventos (=3,2%) e que eventos acima ficar também que os eventos de chuva acima do
do Limiar B constituem entre 0,1 e 0,9 (=0,4%) Limiar A e sem registros de deslizamentos se con-
do total de eventos de chuva registrados. As taxas centram tanto nos eventos de curta duração (até 6
de acerto variam entre 12,5% e 33,3% (=25,5%) h), como nos eventos de chuva de longa duração

TABELA 3 – Taxas de acerto dos limiares de chuva.

TABLE 3 – Success rates of the rainfall thresholds.

Pluviômetro* 01A 02A 03A 04A 05A 06A 07A

18 13 21 24 21 27 21
Elimiar
(2,6%) (3,4%) (3,6%) (3,3%) (2,9%) (3,9%) (2,9%)

Limiar A Elimiar;deslizamento** 5 4 6 3 7 6 6

TA 27,8% 30,8% 28,6% 12,5% 33,3% 22,2% 28,6%

2 2 2 4 1 6 3
Elimiar
(0,3%) (0,5%) (0,3%) (0,5%) (0,1%) (0,9%) (0,4%)

Limiar B Elimiar;deslizamento** 1 0 1 1 1 1 2

TA 50% 0% 50% 25% 100% 16,7% 66,7%

Total de eventos 684 380 582 728 734 699 722

* Os códigos dos pluviômetros foram resumidos, todos eles são precedidos pela sequência “4208906”.
** Para o cálculo de Elimiar;deslizamento do Limiar A foram utilizados todos os eventos com deslizamento e acima do limiar,
enquanto para o Limiar B foram utilizados apenas os eventos com dez ou mais deslizamentos e acima do limiar.

12
 Derbyana, São Paulo, 43: e772, 2022.

FIGURA 9 – Eventos de chuva registrados em cada um dos pluviômetros, com destaque para os even-
tos com ocorrência de deslizamentos. Pode acontecer de eventos com deslizamentos estarem plotados
abaixo dos limiares por não ser o seu pluviômetro de referência. Os eventos com ocorrência no plu-
viômetro de referência são aqueles em que os dados de I e D utilizados para a construção dos limiares
foram coletados naquele pluviômetro e estão destacados. Código dos pluviômetros representados em
cada gráfico: (a) 420890601A, (b) 420890602A, (c) 420890603A, (d) 420890604A, (e) 420890605A,
(f) 420890606A, (g) 420890607A.
FIGURE 9 – Rain events recorded in each of the pluviometers, highlighting the events with landslide occurrences. It may
happen that events with landslides are plotted below the thresholds because it is not their reference rain gauge. The events with
landslide occurrence in the reference rain gauge are those in which the I and D data used to build the thresholds were collected
in that rain gauge and are highlighted. Code of the rain gauges represented in each graph: (a) 420890601A, (b) 420890602A, (c)
420890603A, (d) 420890604A, (e) 420890605A, (f) 420890606A, (g) 420890607A.

13
Stabile et al.

(acima de 24 h). Sabe-se que os limiares inevitavel- mente, os principais indutores dos deslizamentos
mente representam uma simplificação da relação registrados entre 2016 e maio/2022 em Jaraguá do
entre a chuva e a deflagração de movimentos de Sul são os próprios eventos de chuva deflagradores.
massa (ALEOTTI 2004), de modo que não sendo
a chuva a causa direta das rupturas e sim a infiltra- 6 POSSIBILIDADES DE APLICAÇÃO EM
ção da água na encosta e todos os mecanismos ge- SISTEMAS DE ALERTA
otécnicos envolvidos, a ocorrência de deslizamen-
tos irá variar também em função de outros fatores Os limiares intensidade-duração trazem
controladores como as propriedades hidrológicas e novas perspectivas quanto às possibilidades de
físicas do solo, a geomorfologia do terreno, a estru- emissão de alertas de deslizamentos no território
tura geológica, a cobertura vegetal, assim como as brasileiro. Diferentemente dos limiares tradicio-
interferências antrópicas. Estas últimas, no caso de nalmente utilizados que empregam acumulados de
deslizamentos isolados reportados pelas Defesas chuva em 24, 48 ou 72 h, os limiares intensidade-
Civis, certamente têm influência muito grande na -duração são capazes de predizer os eventos de
instabilização das encostas. chuva de curta duração que potencialmente podem
Certamente esta justificativa associada à deflagrar deslizamentos.
subnotificação de ocorrências explica parte dos Contudo, no caso do uso desses limiares em
eventos de chuva acima dos limiares e sem des- um sistema de alertas, há de se ponderar um núme-
lizamentos. Contudo, interessou-nos saber se ro considerável de eventos falsos positivos, sobre-
outros aspectos da chuva, particularmente os pi- tudo para o Limiar A. Isso pode ocorrer, nas taxas
cos de intensidade e os acumulados prévios, não aqui reportadas, pela subnotificação das ocorrên-
poderiam esclarecer por que alguns eventos de cias de deslizamento pelos motivos já descritos no
chuva acima dos limiares não foram capazes de item 3.1, mas também pelos limiares de chuva tra-
deflagrar deslizamentos. Para isso foi calculada a tarem-se de modelos simples de uma realidade am-
intensidade máxima da chuva (em 10 min e em plamente complexa. Por isso, um sistema de alertas
1 h) de cada um dos eventos de chuva acima dos deve envolver necessariamente a conscientização
limiares (Figura 10a 10b) e o acumulado de chuva da população não só acerca dos riscos (INTRIERI
em até 96 h antes do evento considerado (Figura et al. 2013, WMO 2018), mas também das limi-
10c e 10d). tações metodológicas na previsão dos fenômenos.
As figuras 10a e 10b sugerem que nos even- A UNDRR (2015) define que um sistema
tos de chuva com duração acima de 24 h, os picos de alerta corresponde a um sistema integrado de
de intensidade parecem ter tido alguma importân- monitoramento de perigos, previsão e predição,
cia na deflagração de deslizamentos. Verifica-se avaliação de risco de desastres, assim como a siste-
que 80% destes eventos de chuva com deslizamen- mas para comunicação e preparação que permitem
tos tiveram um pico de chuva de pelo menos ~40 que indivíduos, comunidades, governos, empresas
mm/h em 10 min e ~20 mm/h em 1h. Da mesma e outros tomem medidas para reduzir os riscos de
forma, 80% dos eventos de chuva com D > 24h e desastres antes da ocorrência dos eventos perigo-
sem deslizamentos têm picos de chuva abaixo des- sos. Consequentemente, a antecedência é um com-
tes valores. ponente fundamental dos alertas.
De fato, estes resultados corroboram as análi- No exemplo dos limiares aqui definidos para
ses de GUIDICINI & IWASA (1977) e TATIZANA Jaraguá do Sul (SC), o ideal é que o alerta seja emi-
et al. (1987) de que os escorregamentos rasos são tido antes de atingir o Limiar A, que já prediz a
deflagrados por processos geomecânicos relacio- ocorrência de deslizamentos pontuais. Tendo em
nados aos eventos de precipitação de forte inten- conta que os deslizamentos são fenômenos rápi-
sidade. Mesmo eventos de precipitação que geram dos, que ocorrem no espaço temporal de alguns
altos acumulados, especialmente acima de 24 h, segundos ou minutos, a única maneira de prevê-los
precisam de picos de chuva intensa para deflagrar e enviar alertas com antecedência suficiente é utili-
deslizamentos, ainda que menores do que os picos zando a previsão do tempo.
de chuva que instabilizam as encostas em eventos Fatalmente, a previsão do tempo irá adicionar
de menor duração. incertezas ao alerta emitido, que se somam àquelas
Por outro lado, as figuras 10c e 10d não apre- já inerentes aos limiares de chuva, especialmente
sentam dados consistentes com qualquer conclusão para os sistemas meteorológicos convectivos que
a respeito dos acumulados prévios. Consequente- se desenvolvem rapidamente e causam chuvas que

14
 Derbyana, São Paulo, 43: e772, 2022.

(a) (b)

(c) (d)

FIGURA 10 – Gráficos mostrando as relações entre chuva e duração dos eventos acima do Limiar A (para os
eventos de chuva sem deslizamentos foram considerados os dados do pluviômetro 420890607A). Os gráficos
(a) e (b) apresentam a relação entre os picos de intensidade de chuva de cada evento (em 10 min e em 1 h) e a
sua duração. Os gráficos (c) e (d) representam a relação entre os acumulados de chuva prévios (desconsiderado
o evento deflagrador) em até 96 horas e a duração e intensidade dos eventos acima do Limiar A.
FIGURE 10 – Graphs showing the relationships between rainfall and duration of events above Threshold A (for rainfall events
without landslides, data from pluviometer 420890607A were considered). Graphs (a) and (b) show the relationship between the
rainfall intensity peaks of each event (in 10 minutes and in 1 hour) and their duration. Graphs (c) and (d) represent the relationship
between previous rainfall accumulated (disregarding the triggering event) in up to 96 hours and the duration and intensity of events
above Threshold A.

atingem os limiares em poucas horas. Portanto, a visão do tempo apropriados para emissão de alertas
compreensão dos sistemas meteorológicos capazes com a antecedência necessária. Neste sentido, os
de gerar chuvas que ultrapassam os limiares para limiares intensidade-duração fornecem um impor-
deflagração de deslizamentos, bem como sua sazo- tante ponto de partida, pois permitem discernir os
nalidade e frequência, é fundamental para se deter- eventos de curta duração daqueles eventos causa-
minar os equipamentos, modelos e métodos de pre- dos por fenômenos climáticos mais duradouros.

15
Stabile et al.

7 CONCLUSÕES tarina e à Defesa Civil do Município de Jaraguá

do Sul pelas informações acerca dos eventos des-
Este estudo demonstrou que os limiares in- lizamento no município, assim como aos colegas
tensidade-duração são capazes de auxiliar na pre- do Centro Nacional de Monitoramento e Alertas de
visão dos deslizamentos deflagrados pela chuva. Desastres Naturais (Cemaden) que contribuíram
Uma vez implementados em sistemas de alertas, na sistematização destas informações. Também
estes limiares têm o potencial de contribuir para a agradecem aos revisores anônimos cujas valiosas
prevenção de desastres naturais causados por des- contribuições e reflexões originais enriqueceram o
lizamentos. texto.
Praticamente inexistentes no Brasil, a cons-
trução destes limiares depende da existência de
dados sistematizados de chuva com resolução 9 REFERÊNCIAS BIBLIOGRÁFICAS
temporal adequada, bem como do registro de ocor- ALEOTTI, P. 2004. A warning system for
rências de deslizamentos com dados confiáveis de rainfall-induced shallow failures.
localização e horário da deflagração. Informações Engineering Geology, 73: 247-265. https://
como estas só passaram a ter abrangência nacional doi.org/10.1016/j.enggeo.2004.01.007
a partir da instalação dos pluviômetros do Cema-
den e do registro de ocorrências de desastres natu- BAUM, R.L.; GODT, J.W. 2010. Early warning of
rais pelo Cemaden e outros órgãos relacionados à rainfall-induced shallow landslides and debris
prevenção de desastres. flows in the USA. Landslides, 7: 259-272.
Foram elaborados dois limiares para o mu- https://doi.org/10.1007/s10346-009-0177-0
nicípio de Jaraguá do Sul: um para eventos com
pelo menos um deslizamento (Limiar A) e outro BERTOLI, D. 2006. Dinâmica da paisagem da
sub-bacia do ribeirão Chico de Paulo (Ja-
para eventos com dez ou mais deslizamentos (Li-
raguá do Sul - SC): urbanização e conflitos
miar B), com taxas de acerto médias de até 29,8%
decorrentes. Universidade Federal de Santa
e 62,5%, respectivamente. Estimou-se que 3,2%
Catarina, Florianópolis, Dissertação de Mes-
dos eventos de chuva ocorridos entre 2016 e
trado, 117 p.
maio/2022 estão acima do Limiar A e que 0,4% es-
tão acima do Limiar B. A previsão destes eventos CAINE, N. 1980. The rainfall intensity-duration
pode se beneficiar de futuros estudos que tenham control of shallow landslides and debris
como objetivo a compreensão da sua sazonalida- flows. Geografiska Annaler, Series A,
de, bem como dos sistemas meteorológicos a eles Physical Geography, 62(1/2): 23-27. https://
associados. doi.org/10.2307/520449
Os resultados aqui apresentados encorajam
a expansão dos estudos de limiares intensidade- CPRM – SERVIÇO GEOLÓGICO DO BRASIL.
-duração para outros municípios e regiões de 2011. Setorização de áreas de alto e muito
maneira a aperfeiçoar o potencial da previsão de alto risco a movimentos de massa, enchentes
desastres naturais causados por deslizamentos. e inundações: Jaraguá do Sul, SC. CPRM,
Análises posteriores podem explorar o aprimo- Ministério de Minas e Energia, Relatório Téc-
ramento da metodologia aqui proposta ao testar nico. Disponível em https://rigeo.cprm.gov.
outras definições de evento de chuva, assim como br/jspui/handle/doc/18719. Acessado em 01
incluir parâmetros ou modelos que descrevam a set. 2022.
precipitação acumulada ou a umidade anteceden-
CPRM – SERVIÇO GEOLÓGICO DO BRASIL.
te do solo (GODT et al. 2006, BAUM & GODT
2019. Setorização de áreas de alto e muito
2010), bem como explorar o papel dos picos de
alto risco a movimentos de massa, enchentes
intensidade da chuva na deflagração dos desliza-
e inundações: Jaraguá do Sul, SC. CPRM,
mentos.
Ministério de Minas e Energia, Relatório Téc-
nico. Disponível em https://rigeo.cprm.gov.
8 AGRADECIMENTOS br/jspui/handle/doc/18719. Acessado em 01
set. 2022.
Os autores agradecem aos colegas do Centro
Nacional de Gerenciamento de Riscos e Desastres D’ORSI, R.N. 2011. Correlação entre pluviome-
(CENAD), à Defesa Civil do Estado de Santa Ca- tria e escorregamentos no trecho da Serra

16
 Derbyana, São Paulo, 43: e772, 2022.

dos Órgãos da Rodovia Federal BR-116 RJ duration control of shallow landslides and
(Rio-Teresópolis). Universidade Federal do debris flows: an update. Landslides, 5: 3-17.
Rio de Janeiro, Rio de Janeiro, Tese de Dou- https://doi.org/10.1007/s10346-007-0112-1
torado, 287 p.
IBGE – INSTITUTO BRASILEIRO DE GEO-
D’ORSI, R.N.; PAES; N.M.; MAGALHÃES, GRAFIA E ESTATÍSTICA. 2022. Cidades.
M.A.; COELHO, R.S.; COSTA, L.C. 2012. Disponível em https://cidades.ibge.gov.br/
Relatório de correlação entre chuvas e es- brasil/sc/jaragua-do-sul/panorama. Aces-
corregamentos para a cidade do Rio de Ja- sado em 03 set. 2022.
neiro no ano de 2010. Fundação Instituto de
Geotécnica do Município do Rio de Janeiro, INTRIERI, E.; GIGLI, G.; CASAGLI, N.;
Relatório GEO-RIO/DEP/GPE n. 10/2012, NADIM, F. 2013. Landslide Early Warning
58 p. System: toolbox and general concepts.
Natural Hazards and Earth System
FERNANDES, N.F.; GUIMARÃES, R.G.; Sciences, 13: 85-90. https://doi.org/10.5194/
GOMES, R.A.T.; VIEIRA, B.C.; MON- nhess-13-85-2013
TGOMERY, D.; GREENBERG, H. 2001.
Condicionantes geomorfológicos dos desli- KOBIYAMA, M.; MENDONÇA, M.; MORENO,
zamentos nas encostas: avaliação de meto- D.A.; MARCELINO, I.P.V.O; MARCELI-
dologias e aplicação de modelo de previsão NO, E.V.; GONÇALVES, E.F.; BRAZETTI,
de áreas susceptíveis. Revista Brasileira L.L.P.; GOERL, R.F.; MOLLERI, G.S.F.;
de Geomorfologia, 2(1): 51-71. https://doi. RUDORFF, F.M. 2006. Prevenção de De-
org/10.20502/rbg.v2i1.8 sastres Naturais – Conceitos Básicos. Orga-
nic Trading, Curitiba, 109 p.
GODT, J.W.; BAUM, R.L.; CHLEBORAD, A.F.
2006. Rainfall characteristics for shallow LARSEN, M.C.; SIMON, A. 1993. A rainfall
landsliding in Seattle, Washington, USA. intensity-duration threshold for lanslides in
Earth Surface Processes and Landforms, a humid-tropical environment, Puerto Rico.
31(1): 97-110. https://doi.org/10.1002/ Geografiska Annaler, Series A, Physical
esp.1237 Geography, 75(1-2): 13-23. https://doi.org/1
0.1080/04353676.1993.11880379
GOUVEIA, R.L.; CAMPOS, C.C.; MENEZES,
J.T.; MOREIRA, G.F. 2018. Análise de fre- MATURIDI, A.M.A.M.; KASIM, N.; TAIB, K.A.;
quência de precipitação e caracterização de AZAHAR, W.N.A.W; TAJUDDIN, H.B.A.
anos secos e chuvosos para a Bacia do Rio 2021. Empirically based rainfall threshold for
Itajaí. Revista Brasileira de Climatologia, landslides occurrence in Peninsular Malaysia.
22: 309-323. http://dx.doi.org/10.5380/ab- KSCE Journal of Civil Engineering, 25:
clima.v22i0.55276 4552-4566. https://doi.org/10.1007/s12205-
021-1586-4
GUIDICINI, G.; IWASA, O.Y. 1977. Tentative
correlation between rainfall and landslides MENDES, R.M.; VALERIO FILHO, M.; BER-
in a humid tropical environment. Bulletin of TOLDO, M.A.; SILVA, M.F. 2015. Estudo
the International Association of Engineering de limiares críticos de chuva deflagradores
Geology, 16: 13-20. https://doi.org/10.1007/ de deslizamentos no município de São José
BF02591434 dos Campos/SP (Brasil). Territorium, 22:
119-129. https://doi.org/10.14195/1647-
GUZZETTI, F.; PERUCCACCI, S.; ROSSI, M.; 7723_22_8
STARK, C.P. 2007. Rainfall thresholds
for the initiation of landslides in central METODIEV, D.; ANDRADE, M.R.M.; MENDES,
and southern Europe. Meteorology and R.M.; MORAES, M.A.E; KONIG, T.;
Atmospheric Physics, 98: 239-267. https:// BORTOLOZO, C.A.; BERNARDES, T.;
doi.org/10.1007/s00703-007-0262-7 LUIZ, R.A.F.; COELHO, J.O.M. 2018.
Correlation between rainfall and mass
GUZZETTI, F.; PERUCCACCI, S.; ROSSI, M.; movements in North Coast Region of
STARK, C.P. 2008. The rainfall intensity- São Paulo State, Brazil for 2014-2018.

17
Stabile et al.

International Journal of Geosciences, 9: 669- rainfall thesholds for landslide occurence.

679. https://doi.org/10.4236/ijg.2018.912040 Landslides, 15: 1483-1501. https://doi.
org/10.1007/s10346-018-0966-4
MOLINA, E.A.N.; CARDOSO, A.O.; NOGUEI-
RA, F.R. 2015. Relação precipitação-desli- SILVA, F.A. 2014. Percepção dos riscos e desas-
zamento no município de São Bernardo do tres ambientais dos anos 2008 e 2011 em
Campo – SP. Ciência e Natura, 37: 46-54. Jaraguá do Sul – SC. Universidade Federal
https://doi.org/10.5902/2179460X16214 de Mato Grosso do Sul, Três Lagoas, Disser-
tação de Mestrado, 106 p.
PANDOLFO, C.; BRAGA, H.J.; SILVA JÚNIOR,
V.P.; MASSIGNAN, A.M.; PEREIRA, E.S.; SILVA, F.A.; PINTO, A.L. 2013. Vulnerabili-
THOMÉ, V.M.R; VALCI, F.V. 2002. Atlas dade no município de Jaraguá do Sul, SC:
climatológico do Estado de Santa Catarina. Estudo dos efeitos das chuvas de 2008.
Florianópolis, EPAGRI. Periódico Eletrônico Fórum Ambiental
Da Alta Paulista, 9(8): 62-70. https://doi.
PINTO, E.J.A.; AZAMBUJA, A.; FARIAS, org/10.17271/19800827982013466
J.A.M.; SALGUEIRO, J.P.B.; PICKBREN-
NER, K. 2011. Atlas pluviométrico do Bra- TABALIPA, L.K. 2019. Aplicação do Modelo
sil. Rio de Janeiro, Serviço Geológico do SHALSTAB para o mapeamento de áreas
Brasil (CPRM). suscetíveis a deslizamentos translacionais:
estudo de caso em Jaraguá do Sul, SC. Uni-
REICHENBACH, P.; CARDINALI, M.; DE versidade Federal de Santa Catarina, Floria-
VITA, P.; GUZZETTI, F. 1998. Regional nópolis, Dissertação de Mestrado, 150 p.
hydrological thresholds for landslides and
floods in the Tiber River Basin (Central TATIZANA, C.; OGURA, A.T.; CERRI, L.E.S.;
Italy). Environmental Geology, 35(2-3): 146- ROCHA, M.C.M. 1987. Análise de correla-
159. https://doi.org/10.1007/s002540050301 ção entre chuvas e escorregamentos – Serra
do Mar, município de Cubatão. In: ABGE,
SAITO, H.; NAKAYAMA, D.; MATSUYAMA, CONGRESSO BRASILEIRO DE GEO-
H. 2010 Relationship between the initiation LOGIA DE ENGENHARIA, 5, São Paulo,
of a shallow landslide and rainfall intensity- Anais, 2: 225-236.
duration thresholds in Japan. Geomorphology,
118: 167-175. https://doi.org/10.1016/j. UNDRR – UNITED NATIONS OFFICE FOR
geomorph.2009.12.016 DISASTER RISK REDUCTION. 2015.
Global assessment report on disaster risk
SANTA CATARINA. 2016. Atlas geográfico de reduction: making development sustainable
Santa Catarina: diversidade da natureza – The future of disaster risk reduction. UN,
- fascículo 2. Editora UDESC, Florianópo- New York, 316 p.
lis, Secretaria de Estado do Planejamento.
Diretoria de Estatística e Cartografia, 2.ed. WIECZOREK, G.F.; GLADE, T. 2005. Climatic
Disponível em https://www.udesc.br/faed/ factors influencing occurrence of debris
geografia/atlasgeografico. Acessado em 03 flows. In: M. Jakob & O. Hungr (eds.) Debris
set. 2022. flow hazards and related phenomena. Berlin,
Springer, p. 325-362.
SANTORO, J; MENDES, R.M.; PRESSINOTTI,
M.M.N.; MANOEL, G.R. 2010. Correla- WMO – WORLD METEOROLOGICAL
ção entre chuvas e deslizamentos ocorridos ORGANIZATION. 2018. Multi-hazard Warning
durante a operação do Plano Preventivo de Systems: A Checklist. Disponível em https://
Defesa Civil em São Paulo, SP. In: ABGE/ public.wmo.int/en/our-mandate/focus-areas/
CTCGG, SIMPÓSIO BRASILEIRO DE natural-hazards-and-disaster-risk-reduction/
CARTOGRAFIA GEOTÉCNICA E GEO- mhews-checklist. Acessado em 01 set. 2022.
AMBIENTAL, 7, Maringá, Anais, 1-15.
WOLLE, C.M.; CARVALHO, C.S. 1989. Desliza-
SEGONI, S.; PICIULLO, L.; GARIANO, S.L. mentos em encostas na Serra do Mar - Brasil.
2018. A review of the recent literature on Solos e Rochas, 12: 27-36.

18
 Derbyana, São Paulo, 43: e772, 2022.

Endereço dos autores:

Rodrigo Augusto Stabile*, Rafael Alexandre Ferreira Luiz, Andrezza Marques Ferreira, Harideva
Maturano Egas e Márcio Roberto Magalhães de Andrade – Centro Nacional de Monitoramento e Alertas
de Desastres Naturais (Cemaden), Estrada Doutor Altino Bondesan, 500, Distrito de Eugênio de Melo,
CEP: 12.247-016, São José dos Campos, SP, Brasil. E-mails: [email protected], rafael.
[email protected], [email protected], [email protected], marcio.
[email protected]

* Autor correspondente

Artigo submetido em 25 de setembro de 2022, aceito em 29 de novembro de 2022.

This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International License.

19