Movimentos de massa podem afetar a integridade e o valor de uso de patrimônios culturais localizados em áreas suscetíveis à sua ocorrência ou expostas à trajetória de deslocamento destes processos. Diversos estudos foram realizados, sobretudo na Europa, com o objetivo de caracterizar o perigo a que estão expostos patrimônios culturais utilizando métodos quantitativos. No entanto, abordagens baseadas em modelos de base física para esta finalidade ainda são escassas em contextos tropicais úmidos e no Brasil. O objetivo desta pesquisa foi identificar o perigo a movimentos de massa em sítios arqueológicos pré-coloniais do Arquipélago de São Sebastião (SP), previamente identificados como potencialmente ameaçados em uma análise de suscetibilidade em escala regional (1:100.000). Foram aplicados o modelo SHALSTAB para identificar limiares de precipitação diários para deflagração de escorregamentos e o modelo de trajetória de movimentos de massa Gravitational Process Path (GPP). Dados de elevação foram derivados do Modelo Digital de Elevação TanDEM-X (12m), dados de propriedades geotécnicas foram obtidos na literatura e um inventário de cicatrizes foi utilizado para otimizar e validar ambos os modelos em uma área com características geológicas similares na Ilha de São Sebastião. Ambos os modelos aplicados foram submetidos a um experimento fatorial completo para otimização dos parâmetros, utilizando o inventário de cicatrizes e a área abaixo da curva ROC (AAC-ROC) como métrica de qualidade. Os modelos otimizados foram então replicados nas ilhas de Búzios, Vitória e Pescadores, onde se localizam os sítios potencialmente ameaçados. Para o modelo SHALSTAB, foram realizadas 405 combinações diferentes de parâmetros, sendo que o modelo otimizado apresentou AAC-ROC igual a 0.71 quando comparado com cicatrizes não utilizadas no processo de otimização. A otimização dos componentes de trajetória e de fricção do GPP foram realizadas, respectivamente, com 54 e 60 combinações possíveis de parâmetros, resultando em um modelo final com ambas as componentes acopladas possuindo AAC-ROC igual a 0.80. A reprodução dos modelos com parâmetros otimizados nas áreas com sítios potencialmente ameaçados resultou em dois cenários: a) um com as áreas de iniciação dos escorregamentos para o modelo GPP sendo aquelas classificadas pelo SHALSTAB como incondicionalmente instáveis; e b) áreas de iniciação associadas a limiares de precipitação diária do SHALSTAB entre 50 e 100mm. No primeiro cenário, dois dos 12 sítios são atravessados por 55% e 62% das simulações de trajetória, com velocidades de deslocamento do material entre 1.23 e 9.7 m/s. No segundo cenário, 8 sítios são afetados pelos movimentos de massa, todos eles atravessados por 60% ou mais das simulações de trajetória e com velocidades de deslocamento da massa deslizada entre 5.8 e 9.7 m/s. Conclui-se que os sítios identificados em escala regional estão de fato ameaçados por movimentos de massa, com probabilidade de afetarem sítios em ao menos 50% dos casos de precipitações abaixo de 50mm que deflagrem movimentos de massa e superior a 60% para precipitações entre 50 e 100mm. As velocidades estimadas indicam potencial de danos severos ao patrimônio cultural, com potencial soterramento dos sítios ou mobilização do material arqueológico.
