Deadly disasters in southeastern South America: flash floods and landslides of February 2022 in Petrópolis, Rio de Janeiro
←
→
Page content transcription
If your browser does not render page correctly, please read the page content below
Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023 https://doi.org/10.5194/nhess-23-1157-2023 © Author(s) 2023. This work is distributed under the Creative Commons Attribution 4.0 License. Deadly disasters in southeastern South America: flash floods and landslides of February 2022 in Petrópolis, Rio de Janeiro Enner Alcântara1 , José A. Marengo1,2 , José Mantovani1 , Luciana R. Londe1,2 , Rachel Lau Yu San3 , Edward Park3 , Yunung Nina Lin4 , Jingyu Wang3 , Tatiana Mendes1,5 , Ana Paula Cunha1,2 , Luana Pampuch1,5 , Marcelo Seluchi2 , Silvio Simões1 , Luz Adriana Cuartas1,2 , Demerval Goncalves2 , Klécia Massi1,5 , Regina Alvalá1,2 , Osvaldo Moraes1,2 , Carlos Souza Filho6 , Rodolfo Mendes1,2 , and Carlos Nobre1,7 1 Graduate Program in Natural Disasters, Unesp/CEMADEN, São José dos Campos, Brazil 2 National Center for Monitoring and Early Warning of Natural Disasters (CEMADEN), São José dos Campos, Brazil 3 National Institute of Education, Earth Observatory of Singapore and Asian School of the Environment, Nanyang Technological University (NTU), Singapore 4 Institute of Earth Sciences, Academia Sinica, Taipei, Taiwan 5 Department of Environmental Engineering, Institute of Science and Technology, São Paulo State University (Unesp), São José dos Campos, Brazil 6 Institute of Geosciences (IG/Unicamp), University of Campinas, Campinas, Brazil 7 Institute of Advanced Studies, University of São Paulo (IEA/USP), São Paulo, Brazil Correspondence: Enner Alcântara (enner.alcantara@unesp.br) Received: 2 June 2022 – Discussion started: 10 June 2022 Revised: 15 February 2023 – Accepted: 5 March 2023 – Published: 21 March 2023 Abstract. On 15 February 2022, the city of Petrópolis in ducive for landslides as they recorded landslide occurrences the highlands of the state of Rio de Janeiro, Brazil, received of about 9 % to 11 %. Regarding the soil moisture, higher an unusually high volume of rain within 3 h (258 mm), gen- variability was found in the lower altitude (842 m) where the erated by a strongly invigorated mesoscale convective sys- residential area is concentrated. Based on our land deforma- tem. It resulted in flash floods and subsequent landslides that tion assessment, the area is geologically stable, and the land- caused the deadliest landslide disaster recorded in Petrópo- slide occurred only in the thin layer at the surface. Out of the lis, with 231 fatalities. In this paper, we analyzed the root 1700 buildings found in the region of interest, 1021 are on causes and the key triggering factors of this landslide dis- the slope between 20 to 45◦ and about 60 houses were di- aster by assessing the spatial relationship of landslide oc- rectly affected by the landslides. As such, we conclude that currence with various environmental factors. Rainfall data the heavy rainfall was not the only cause responsible for the were retrieved from 1977 to 2022 (a combination of ground catastrophic event of 15 February 2022; a combination of un- weather stations and the Climate Hazards Group InfraRed planned urban growth on slopes between 45–60◦ , removal of Precipitation – CHIRPS). Remotely sensed data were used to vegetation, and the absence of inspection were also expres- map the landslide scars, soil moisture, terrain attributes, line- sive driving forces of this disaster. of-sight displacement (land surface deformation), and urban sprawling (1985–2020). The results showed that the average monthly rainfall for February 2022 was 200 mm, the heavi- est recorded in Petrópolis since 1932. Heavy rainfall was also 1 Introduction and background recorded mostly in regions where the landslide occurred, ac- cording to analyses of the rainfall spatial distribution. As for The municipality of Petrópolis is nestled in the moun- terrain, 23 % of slopes between 45–60◦ had landslide occur- tains, 68 km from the city of Rio de Janeiro. It presents a rences and east-facing slopes appeared to be the most con- rugged relief with numerous cliffs, and it is populated by ap- proximately 305 687 inhabitants. The city is predominantly Published by Copernicus Publications on behalf of the European Geosciences Union.
1158 E. Alcântara et al.: Deadly disasters in southeastern South America urban (IBGE: https://cidades.ibge.gov.br/brasil/rj/petropolis/ count in the state was over 200. Approximately 160 others panorama, last access: 23 May 2022) (95.1 %) (Fig. 1). It has were injured and 15 000 had to leave their homes across the been ravaged yet again by mudslides and floods on February, state. By the time the rain stopped, several days later, 12 000 2022, when heavy rains triggered landslides that left hun- people were homeless (https://www.france24.com/en/ dreds dead or missing. On that day, there was an unusually 20100406-rio-de-janeiro-flooding-death-toll-passes-100, high amount of rain within 3 h, 258 mm, as measured by rain last access: 23 May 2022). gauges in the city. This was more than the prior 30 d accu- The January 2011 disaster was one with the highest im- mulated rainfall and triggered the worst case of heavy rain- pact in the history of disasters in Brazil (Rosi et al., 2019). fall in the city since 1932; these data were obtained by IN- With over 260 mm pouring down in less than 24 h, the rain MET (Brazilian Institute of Meteorology) using data from triggered a series of mudslides throughout the região ser- a rain gauge in Petrópolis (number 83804), which operates rana (mountainous region) northwest of the state of Rio de from 1 October 1912, to 1 August 1960. According to the Janeiro, causing death and destruction in five nearby cities: National Center for Monitoring and Early Warning of Natu- Teresópolis, Petrópolis, Nova Friburgo, Sumidouro, and São ral Disasters (Cemaden, http://www.cemaden.gov.br, last ac- José do Vale do Rio Preto. Nova Friburgo and Teresópolis cess: 7 June 2022), 250 mm of rain was recorded between each counted over 250 deaths with the final death toll of 16:20 and 19:20 UTC (February climatology is 185 mm). 900 and 300 still missing in the região serrana. In Teresópo- The previous record occurred in January 2011, with accu- lis alone, 1000 people were left without homes (Rosi et mulated rainfall on the order of 241.8 mm in 24 h and a peak al., 2019) and claimed the lives of 73 people in Petrópolis of 61.8 mm in 1 h (Dourado et al., 2012). (Marengo and Alves, 2012). Such precipitation occurred be- The abundant precipitation caused flash floods and cause of the formation of an intense episode of rainfall gen- mudslides in various sectors of the city, tearing down erated by the action of the South Atlantic Convergence Zone dozens of homes on the hillsides and causing floods in (SACZ) that was intensified by a deep trough on the south- the streets. Images and videos on social media showed eastern coast of Brazil. The SACZ is a typical meteorological rivers of mud rushing through the city’s streets, sweep- system of the rainy season in the central area of Brazil, which ing away cars, trees, and people in its way. According normally develops between November and March. The scale to CNN (https://edition.cnn.com/2022/02/17/world/gallery/ of the disaster brought attention to the rapid urbanization, in- brazil-landslides/index.html, last access: 23 May 2022) and frastructure, and housing rights issues within the context of CPRM (2022), Brazil’s Civil Defense Secretariat reported disasters in Rio de Janeiro and Brazil, leading to the modern- 269 landslides during this event. By the end of February, ization of the civil defense structure in Brazil and the creation the death toll had risen to 231 (https://floodlist.com/america/ of CEMADEN. brazil-floods-landslides-petropolis-february-2022, last ac- Once again, on 20 March 2022, torrential rains caused cess: 23 May 2022). flooding and landslides in the metropolitan area of Petrópo- While this is the deadliest flood and mudslide in the history lis and surrounding areas of the mountainous region of of Petrópolis, heavy rains are not uncommon during Brazil’s the state of Rio de Janeiro (https://floodlist.com/america/ austral summer months (November to March; warm/rainy brazil-floods-landslides-petropolis-march-2022, last access: season). 20 April 2022). Petrópolis Civil Defense reported 415 mm of According to Guerra (1995), Petrópolis was impacted by rain in the São Sebastião district in just 10 h, with 217.4 mm 1161 catastrophic events between 1940 and 1990, which in- of that total falling in a 4 h period between 14:00 and clude landslides, mudslides, rockfalls, and floods. Most of 18:00 UTC. Around 100 incidents of landslides and floods these events were caused by torrential precipitation. The have been reported across the city, including flooded roads number of deaths appears to increase over time, and nearly and buildings damaged or destroyed, as well as seven deaths. 90 % of the events occurred within urban areas. Moreover, The above examples demonstrated the adverse impacts of the human fatalities are also found to directly correlate with extreme weather and subsequent landslides on the people the spread of new urban settlements onto steep deforested of Petrópolis and highlight the importance of assessing the hillsides. These disasters were responsible for the deaths of causes of these landslides for future disaster management. 526 people in the last 50 years, though seemingly worse over Therefore, this study aims to better understand the main trig- the last 20 years with 300 deaths. gers for landslides and flash floods in Petrópolis, specifically On 5 February 1988, it poured rain in Petrópolis, which for the disaster in February 2022. To achieve this objective, partially destroyed the city. This tragedy caused 171 deaths, we performed an integrated analysis of the urban sprawl over leaving more than 4000 people homeless and thousands exposed areas, removal of vegetation, and shallow saturated injured (Guerra, 1995). On 5 April 2010, precipitation in soil. Rio de Janeiro was the highest recorded in the past 30 years: weather stations recorded 288 mm within 24 h, well over the average rainfall amount for April (140 mm). It caused 52 fatalities in the city of Rio de Janeiro, and the total death Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023 https://doi.org/10.5194/nhess-23-1157-2023
E. Alcântara et al.: Deadly disasters in southeastern South America 1159 Figure 1. Petrópolis and the study site. (a) Location of the state of Rio de Janeiro in Brazil, (b) location of the study site in the state of Rio de Janeiro. Panels (a)–(b) were based on the Brazilian Institute of Geography and Statistics (IBGE). (c) Location of the study site in the city of Petrópolis and the position of weather stations and (d) elevation map of the study site (ALOS PALSAR, https://asf.alaska.edu/data-sets/ sar-data-sets/alos-palsar/, last access: 10 May 2022). (e) Planet image natural color from 14 February 2022 and (f) Planet image natural color from 22 February 2022. 2 Data and methods tion on topography since the entire region was submitted to tectonic events during the Precambrian period (Penha et al., 2.1 Study area 1981; Fonseca et al., 1998; Gonçalves, 1998; IBGE 2018; Rosi et al., 2019). The drainage network of the region is Petrópolis is located within the Atlantic Forest biome, a bio- strongly influenced by brittle regional structures, which play diversity hotspot (Myers et al., 2000) (Fig. 1). Currently, only an essential role in its organization and the relief pattern 13 % of the original cover remains (Fundação SOS Mata and modeling. This set of geological characteristics, such Atlântica, 2018), due to intense deforestation and human dis- as highly foliated and fractured rocks, triggers mass move- turbance that mostly occurred in the first half of the 19th ments, mainly shallow landslides, and debris flow. Processes century (Dean, 1996). Historically, land degradation in the of folding, fault reactivation, and block remobilization re- region is associated with a combination of different geomor- sulted in Petrópolis’s terrain (Gonçalves, 1998), which pro- phic processes, deforestation overexploitation (Nehren et al., duced braided drainage pattern, with elongated hills which 2018), and urban expansion (Guerra, 1995; Rosi et al., 2019). separate valleys. The soil is mostly occupied by clay with Petrópolis is geologically located in the Rio Negro low fertility, is well drained, and has high aluminum satura- Complex, of Paleoproterozoic origin, formed mainly by tion (Carvalho-Filho et al., 2000). migmatites and granitoids. This part of the state of Rio de Janeiro suffered the effects of several regional metamorphic 2.2 Rainfall data and analysis phases, resulting in highly foliated rocks cut by large duc- tile shear zones. These rocks are severely sectioned by frac- To identify the extreme rainfall event in the municipality tures and faults of regional extension, with a strong reflec- of Petrópolis on 15 February 2022, different precipitation https://doi.org/10.5194/nhess-23-1157-2023 Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023
1160 E. Alcântara et al.: Deadly disasters in southeastern South America datasets were retrieved from INMET and CEMADEN station Table 1. Categorization of slope aspect into classes of 45◦ . data. CEMADEN has 17 rain gauges installed in Petrópo- lis, with measurements available from 2015 onwards, al- Aspect (◦ ) Direction though some stations have a large amount of missing data −1 Flat for long periods, especially in the early years of the series. 0–22.5 and 337.5–360◦ N CEMADEN also has five geotechnical stations, installed at 22.5–67.5◦ NE the end of 2021, consisting of sensors that monitor rainfall 67.5–112.5◦ E and soil moisture at different depths (0.5 to 3.0 m). 112.5–157.5◦ SE At the local scale, ground rain gauges provide direct-point 157.5–202.5◦ S estimations necessary for extreme rainfall analysis. However, 202.5–247.5◦ SW a long precipitation time series is required for studying the 247.5–292.5◦ W extreme rainfall event in Petrópolis. Since there is no single 292.5–337.5◦ NW weather station with enough temporal coverage, rainfall data from five different weather stations located up to 15 km from the highest accumulated rainfall location (Fig. 1c) were con- sidered. The hourly rain gauge observations from June 1976 2.3 Terrain analysis using DEM to February 2022 were obtained from the Rio de Janeiro State Environmental Institute – Inea (http://www.inea.rj.gov.br/ The terrain morphology of Petrópolis could have also in- ar-agua-e-solo/monitoramento-hidrometeorologico/, last ac- fluenced the slope stability (Costanzo et al., 2012), causing cess: 10 May 2022) and the Cemaden (https://www.gov.br/ certain areas to be more susceptible to landslides than oth- cemaden/pt-br, last access: 10 May 2022) databases. ers. Terrain attributes widely used in landslide susceptibil- While we initially considered using satellite-based prod- ity assessments include elevation, slope angle, aspect, and ucts for studying rainfall, it proved to be inaccurate; curvature (Catani et al., 2013; Chen et al., 2017; Reichen- thus, this approach was not carried on. Gridded data for bach et al., 2018; Dias et al., 2021). Different elevation cre- global (Climate Hazards Group InfraRed Precipitation – ates diverse environmental conditions in temperature, rainfall CHIRPS; Funk et al., 2015) and regional domains (Satellite- regimes, and vegetation (Dai and Lee, 2002; Costanzo et al., based Global Precipitation Measurement (GPM) – Integrated 2012; Catani et al., 2013), which also influence human de- Multi-satellite Retrievals for GPM (IMERG) combined with velopment (Chau and Chan, 2005). Slope angle affects the data from surface observations – MERGE, from the Cen- shear stress acting on the slope and has been considered one ter for Weather Forecasts and Climate Studies – CPTEC; of the biggest determinants for landslide occurrences (Catani Rozante et al., 2010) were found to exhibit poor capabilities et al., 2013; Kanwal et al., 2017). Slope aspect refers to the in capturing the extreme rainfall event in Petrópolis (Huff- direction of a slope and determines sunlight exposure, often man et al., 2020; Reis et al., 2020). Both products underesti- correlating with moisture retention and vegetation cover and mated the storm-accumulated rainfall by 90 % on 15 Febru- consequently landslide initiation (Guzzetti et al., 1999; Dai ary 2022, throughout the study site. Previous studies have and Lee, 2002). Slope curvature is the rate of change of the also shown that retrieving rainfall extremes from satellite- slope, and this controls the direction of landslide motion by based products can present certain uncertainties concerning concentrating or dispersing surface runoff and gravitational accuracy probability due to instrument issues and retrieval stresses (Ohlmacher, 2007; Costanzo et al., 2012). algorithms or the interpolation process (Dembélé and Zwart, To assess the effects of terrain morphology on this land- 2016; Hermance and Sulieman, 2018; Jiang and Bauer- slide event, the digital elevation model (DEM) was obtained Gottwein, 2019). from ALOS PALSAR, available at the Alaska Satellite Fa- To examine the synoptic weather pattern that spawned the cility (ASF). ALOS PALSAR images were resampled (cubic mesoscale convective system (MCS) with unprecedented tor- interpolation) from 30 to 12.5 m pixel size with orthometric rential rainfall, the fifth generation of atmospheric reanal- altitude (EGM96 geoid model) before being converted to a ysis (ERA5) produced by European Centre for Medium- geometrical altitude (ellipsoidal). The following terrain at- Range Weather Forecasts (ECMWF; Hersbach et al., 2020) tributes were derived from the DEM: slope angle, slope as- is adopted for the analysis of the atmospheric environment. pect, and slope curvature. The spatial variation of the three ERA5 features a horizonal resolution of 31 km and 137 slope attributes within this landslide extent was assessed. The vertical levels with hourly output frequency, which pro- slope angle was classified into five classes: 0–15, 15–25, 25– vides a detailed picture of the MCS-hatching synoptic back- 35, 35–45, and 45–60◦ . Slope aspect was categorized into ground. Meanwhile, the morphology and organization of classes of 45◦ (Table 1). high-precipitation MCS is analyzed using the National Aero- Slope curvature was reclassed into three categories – up- nautics and Space Administration (NASA) merged geosta- wardly convex (negative), flat (zero), and upwardly concave tionary satellite half-hourly 4 km resolution infrared bright- (positive). To account for the slope morphology of the en- ness temperature data (merged IR; Janowiak et al., 2017). tire study area, the landslide extent within each class of slope Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023 https://doi.org/10.5194/nhess-23-1157-2023
E. Alcântara et al.: Deadly disasters in southeastern South America 1161 Table 2. Satellite imagery used in this study. et al. (2004): Sensor/satellite Date LSTmax − LST SMI = , (1) LSTmax − LSTmin OLI and TIRS/Landsat 8 28 January 2021 OLI and TIRS/Landsat 8 17 March 2021 where LSTmax and LSTmin are, respectively, the maximum OLI and TIRS/Landsat 8 2 April 2021 and minimum values of LST within the image for a given OLI and TIRS/Landsat 8 7 July 2021 NDVI, expressed as OLI and TIRS/Landsat 8 23 July 2021 OLI and TIRS/Landsat 8 24 August 2021 LSTmax = a1 × NDVI + b1 (2) OLI and TIRS/Landsat 8 9 September 2021 OLI-2 and TIRS-2/Landsat 9 23 January 2022 LSTmin = a2 × NDVI + b2 , (3) PlanetScope 14 February 2022 PlanetScope 22 February 2022 where a and b are empirical parameters defining the dry and wet edges modeled as a linear fit to the data (Parida et al., 2008). angle, aspect, and curvature was calculated in terms of pro- 2.5 Interferometric synthetic aperture radar (InSAR) portion to the study area. analysis of ground deformation 2.4 Optical remote sensing data and soil moisture Ground deformation was studied using a multi-temporal pattern interferometric synthetic aperture radar (MTInSAR) time- series analysis. C-band Sentinel-1B SAR data of descend- Soil moisture was analyzed through images from the Oper- ing path 155 acquired between 8 May 2015 and 20 Decem- ational Land Imager (OLI) and the Thermal Infrared Sensor ber 2021 were collected and processed in InSAR Scientific (TIRS) on board of the Landsat 8 and 9 satellites. Images Computing Environment (ISCE) to form a coregistered stack from the Planet satellites were also selected (see Table 2). of single look complex (SLC) images. No new images were Images from Landsat satellites are surface reflectance prod- acquired between 21 December 2021 and 22 February 2022, ucts at 30 m resolution, which were obtained from the USGS due to the Sentinel-1B anomaly, meaning that there are no Earth Explorer site (https://earthexplorer.usgs.gov, last ac- data acquired right before and after the landslides in mid- cess: 1 May 2022). The Planet constellation consists of more February 2022. The coregistered SLC stack is then processed than 130 orbital Earth observation satellites acquiring daily in Fine Resolution InSAR using Generalized Eigenvectors images of the Earth on four bands in the visible and near- (FRInGE; Fattahi et al., 2019) to form a wrapped phase time infrared wavelength range at 3–0 m spatial resolution. series via a phase linking approach (Ansari et al., 2018), The soil moisture index (SMI; Lambin and Ehrlich, 1996; which is then unwrapped epoch by epoch using SNAPHU to Zhan et al., 2004) was used to reconstruct the dynamics of generate the final displacement time series (Chen and Zebker, soil moisture from 2021 to the beginning of 2022. This in- 2001). Tropospheric correction is carried out on the displace- dex can be estimated from land surface temperature (LST) ment time series using the ERA5 weather model (Jolivet et and normalized difference vegetation index (NDVI). Accord- al., 2014). Line-of-sight (LOS) velocities and velocity errors ing to Lambin and Ehrlich (1996), the scatterplot of LST vs. are then estimated using the L1-norm iteratively reweighted NDVI results in a trapezoidal shape, and all types of land least squares algorithm (Schlossmacher, 1973). cover fall within the trapezoid of the LST–NDVI space. The upper envelope of the trapezoid (upper limit of surface tem- 2.6 Mapping urban sprawling and forest losses perature for a given vegetation cover) represents the dry con- dition (warm edge), while the lower limit represents the wet The urban sprawl analysis was performed using data avail- condition (cold edge) (Parida et al., 2008). able at the MapBiomas project – Collection 6 (http:// Landsat 8 and 9 images were used to calculate the SMI. mapbiomas.org, last access: 1 May 2022), with spatial reso- OLI/Landsat 8 and 9 Level-2 satellite images were obtained lution at 30 m. The MapBiomas project provides a historical at the USGS Earth Explorer website (https://earthexplorer. series of land use and land cover (LULC) information and the usgs.gov) from 2021 to 2022. Level-2 data are atmospheri- transition data between 1985 and 2020 to the whole coun- cally corrected and generated from the LASRC (Landsat 8 try, based on a random forest algorithm applied to Landsat Surface Reflectance Code), which yields surface reflectance archive using Google Earth Engine (Souza et al., 2020). The at a 30 m spatial resolution suitable for studying the dynam- LULC data of eight years (1985, 1990, 1995, 2000, 2005, ics of our study site (Silva et al., 2017). OLI images were 2010, 2015, and 2020) were analyzed focusing on the urban used to calculate the NDVI, whereas TIRS/Landsat 8 and 9 area and forest formation classes to observe the growth of were used to calculate LST. SMI values, at a range between 0 occupation in the Petrópolis municipality and the area of in- (drier soil) and 1 (wet soil), are calculated according to Zhan terest (AOI, see Fig. 1f). Urban sprawl and loss of vegetation https://doi.org/10.5194/nhess-23-1157-2023 Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023
1162 E. Alcântara et al.: Deadly disasters in southeastern South America (suppressed forest) were visual and quantitatively analyzed currence of high-precipitation MCS, which later produced through the transition data for seven 5-year periods, from the unprecedented torrential rainfall that caused the flood and 1985 to 2020. According to Pisano et al. (2017) land cover landslide. can be considered the most important conditioning factor in landslide susceptibility analysis. 3.2 Rainfall variability 2.7 Deriving the landslide scar map Figure 3 shows the January–February daily rainfall distribu- tion from 1977 to 2022. By comparing the shape of the box- The building footprints were compared to landslide scars plot for 2022 with the other years, the third quartile of 2022 to quantify either impacted or destroyed buildings because was observed to be the biggest of the time series (22.35 mm). of this event. Moreover, these building footprints were con- It indicates that in January and February 2022, 25 % of the fronted with the slope data to obtain the slope range they days accumulated daily precipitation higher than 22.35 mm; occupy on the hill. The building footprints were obtained 2022 also registered three outliers (circles). Compared with by the OpenStreetMap (OSM) database. The OSM is a the other years, this number is not so expressive, but the crowdsourced mapping project that aims to create and pro- records stand out for being abnormally high (rainfall accu- vide a freely available geographical information database of mulated of 88.8, 219.6, and 260 mm in one day). In 2022, the world (Minghini and Frassinelli, 2019; Goldblatt et al., the most intense outliers were also higher than in all other 2020). The mapping of this area was performed as an ur- years of the series (only four events greater than 189.6 mm) gent project after the event on 17 February, which resulted in (Fig. 3). The value registered on 15 February 2022 (260 mm) flooding and landslides in several parts of the municipality. is the second largest value for that month in the 46-year se- ries. 3 Results An intense rainfall, starting on the evening of 15 Febru- ary 2022, caused mudslides on the hillsides above downtown 3.1 Synoptic background analysis Petrópolis and produced flooding that did more damage in the streets below. Images and videos on social media showed Right before the onset of record-breaking rainfall, central rivers of mud rushing through the city’s streets, sweeping ev- Brazil was dominated by a widespread subtropical high- erything along the way: cars, buses, trees, and sometimes pressure center at 500 hPa level (Fig. 2a), whose southern people. Petrópolis recorded 252.8 mm of rain in just 3 h on boundary pushed to 22◦ S parallel above Petrópolis. As re- 15 February (São Sebastião station – 330390604G) (loca- vealed in many previous studies (e.g., Maddox et al., 1978; tions 1 and 2 in Fig. 5). The average for the entire month Mitchell et al., 1995; Wang et al., 2019), the atmosphere di- of February is approximately 200 mm. Landslides and flood- rectly under the subtropical high-pressure center is most sta- ing events in Rio de Janeiro happen most frequently during ble because a subsidence inversion layer (capping inversion) the rainy season (from December to April austral summer). is formed as a result of widespread descending air. However, Such hourly rates and amounts have not been recorded in the atmosphere becomes more unstable toward the edge of Brazil and are rare even for other parts of the world. the high pressure because of the weakening inversion, which From our results, the mean precipitation for the region in promotes the occurrence of convection. Consequently, a ring January is 304.19 mm, while in February it is 229.28 mm of precipitation area commonly forms at the periphery of the (Fig. 4). January 2022 registered an accumulated rainfall of high pressure, known as a “ring of fire” pattern (Galarneau 581.4 mm, the second biggest value in entire series. In Febru- and Bosart, 2006). The 850 hPa specific-humidity distribu- ary, the accumulated rainfall was 650 mm, the largest value tion (Fig. 2b) echoes this pattern, as the warm and moist air is for the 46-year series and considered an outlier for the month being pushed southeastward towards the edge of the subtrop- (circle). ical high, and the moisture is further advected to Petrópolis, For the 15 February 2022, event, Fig. 4a depicts the which fuels the convection there. Figure 2c exhibits the con- station-based 24 h accumulated rainfall over Petrópolis. The current cloud distribution, where the cold colors (< 241 K) highest rainfall accumulation was observed at São Sebastião represent the presence of deep clouds whose cloud top tem- station (330390604G), flagged as “1” in Fig. 5. perature is cooler than the surrounding area without clouds, It recorded about 260 mm in 24 h and approximately and an MCS is formed over Petrópolis. In addition to the 230 mm in two 3 h intervals from 16:00 to 21:00 LT (local influence of the high-pressure center, the MCS that affects time). Another nearby station – Dr. Thouzet (330390603G) Petrópolis is further intensified by a long-wave trough as in- – also registered 221 mm in 24 h (flagged as “2” in dicated in Fig. 2a. The city is located at the bottom of the Fig. 4). Other stations registered the following: Indepen- 500 hPa trough where there is a strong lifting of air that facil- dencia 2 (330390612A), 147 mm (“3” in Fig. 4), Qui- itates the intensification of convection. In summary, the edge tandinha (330390601G), 143 mm (“4” in Fig. 4), Bin- of a subtropical high intersects with a long-wave trough over gen (30390605G), 142 mm (“5” in Fig. 4) and Rua Ama- Petrópolis, forming a very favorable environment for the oc- zonas/Quitandinha (330390618A), 131 mm (“6” in Fig. 4). Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023 https://doi.org/10.5194/nhess-23-1157-2023
E. Alcântara et al.: Deadly disasters in southeastern South America 1163 Figure 2. Large-scale environmental variables and satellite observation at 16:00 LT. (a) Wind speed (shading), geopotential height (contour), and winds (arrow) at 500 hPa, where the trough line is indicated using a dashed brown line and the location of Petrópolis is highlighted by a yellow triangle. (b) Specific humidity (shading) and winds (arrow) at 850 hPa. (c) Merged IR brightness temperature, where cold color represents the presence of deep clouds. Figure 3. The long-term precipitation time series was created from two weather stations from INMET (A603 – Duque de Caxias and A610 – Pico do Couto), two weather stations from INEA (2243238 – Xeren and 2243235 – Andorinhas), and one from Cemaden (São Sebastião). The INEA hydrological station, Alto da Serra (2243315), was the only one with available data for this event. The peak level was 3.59 m (Fig. 4b) for an accumulated rainfall of 207 mm until the peak (the total rainfall was 223 mm). The maximum peak previously registered was 2.58 m on 18 March 2013 (data period 2012–2021), for accumulated rainfall of 186 mm. It is noteworthy that a heavy precipita- tion particularly affected the urban part of the city, where risk areas are located. Further away, over the rural part of the municipality, no rainfall was recorded. 3.3 Characterization of the terrain to assess landslide Figure 4. January and February rainfall (mm) distribution (boxplot) susceptibility for 1977–2022. The blue line indicates the mean value for the 1977– 2022 period (for the months of January and February), and the red line represents the accumulated rainfall for 2022. Circles represent Located in the mountainous regions, Petrópolis is at a rela- outliers. tively high elevation, 953.5 m from 821 to 1086 m (Fig. 6a). https://doi.org/10.5194/nhess-23-1157-2023 Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023
1164 E. Alcântara et al.: Deadly disasters in southeastern South America Figure 5. (a) Accumulated precipitation in 24 h on 15 February 2022 at CEMADEN rain gauges. Red contours correspond to areas of risk for landslides at Petrópolis. The black line indicates the municipality limit. The cartographic base was obtained from © OpenStreetMap contributors 2022. Distributed under the Open Data Commons Open Database License (ODbL) v1.0. (b) Rainfall (blue) and river level (red) time series for 15 February 2022 time event, at Alto da Serra hydrological station (INEA). Although the landslide scars were found all elevations, ceptibility (Costanzo et al., 2012; Reichenbach et al., 2018). most of them were above 950 m elevation. As for the slope However, complementary analysis is required for the occur- morphology, the highest slope angle recorded in this study rence of landslides in slopes lower than 25◦ , as they may be area is approximately 60◦ (Fig. 6b), the east-facing slopes associated with the execution of slopes of cuts and/or fills accounts for the largest area (Fig. 6c), and most of the slopes (Mendes et al., 2018a, b; Ávila et al., 2021). were either concave or convex (Fig. 6d). Spatial analysis re- For the aspect factor, the east-facing slopes appear to be vealed that, out of the five categories of slope angles, slope the most conducive for landslides, as they account for more angles of 45–60◦ had the highest percentage of landslide oc- than 30 % of landslide occurrences (Fig. 7b). Conversely, currence with 23 % (Fig. 7a). There was also a distinct in- west-facing slopes had the lowest landslide occurrence with creasing trend in the proportion of landslide occurrences as only 1 % experiencing landslides. This could be attributed the slope becomes steeper. This corresponds to numerous to multiple reasons such as the exposure of the surface to previous findings that steeper slopes increase slope instabil- varying extents of wind and solar radiation prior to the land- ity due to greater resistance required to maintain stability slide event (Guzzetti et al., 1999; Dai and Lee, 2002). This (Catani et al., 2013; Reichenbach et al., 2018). It is often one could have conditioned the soil moisture, humidity, as well of the most important factors when assessing landslide sus- as vegetation growth on the slopes (Catani et al., 2013; Re- Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023 https://doi.org/10.5194/nhess-23-1157-2023
E. Alcântara et al.: Deadly disasters in southeastern South America 1165 Figure 6. Distribution map of terrain attributes in the study area of Petrópolis. (a) Elevation, (b) slope, (c) aspect, and (d) curvature of Petrópolis. Light-yellow polygons represent the extent of the February 2022 landslide patches. Figure 7. Percentage of the landslide area in relation to the total study area, in each classification of (a) slope (blue trend line), (b) aspect, and (c) curvature. ichenbach et al., 2018). Vegetation growth in Petrópolis was (Fig. 7c). While upwardly convex (negative curvature) and observed to be more abundant on west-facing slopes, which concave (positive curvature) slopes (6 %) had a higher per- could have increased the slope stability in that direction due centage of landslide occurrences than flat curvatures (4 %), to their extensive root systems. Additionally, the rainfall di- the difference is small at about 2 %. Based on the spatial rection during the storm could have also favored landslides analysis of the terrain morphology, this Petrópolis landslide in the eastward orientation. event was likely associated with the terrain morphology of In terms of slope curvature, these landslides did not ap- the mountain even though the trigger was induced by the in- pear to vary substantially between the different curvatures tensive rainfall during the storm. https://doi.org/10.5194/nhess-23-1157-2023 Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023
1166 E. Alcântara et al.: Deadly disasters in southeastern South America 3.4 Soil moisture trend prior to the landslides Table 3. Percentage of change in LULC classes over 5-year periods for the Petrópolis municipality. The estimations of SMI from January 2021 to January 2022 are shown in Fig. 8. Values near zero represent dry soil condi- Period Forest Urban Rocky Pasture and Other tions, and values near 1 mean wet soil conditions. The results area outcrop agriculture show that the wet soil conditions were more frequently in the 1985–1990 3.54 22.45 −0.37 −6.39 −13.53 highest altitude and dryer lower altitude during the summer 1990–1995 3.00 7.04 0.70 −4.97 −13.05 from January 2021 to April 2021 (compare Fig. 6a for eleva- 1995–2000 −0.72 4.83 −0.35 0.43 −9.74 tion). The soil moisture becomes higher in the south portion 2000–2005 0.84 4.16 1.81 −2.06 10.64 during the winter (July to September 2021), mainly lower 2005–2010 1.81 2.48 2.80 −3.47 23.03 2010–2015 −1.46 4.38 2.14 1.22 28.63 altitude, with one exception in August, when it resembles 2015–2020 −0.04 3.69 −2.92 −0.26 21.74 the summer pattern. In January 2022, the moisture condition spread throughout the area, just like in January 2021. The spatial standard deviation image showed a higher variability Table 4. Urban sprawl and forest loss classes in the area of interest. in lower altitude than in higher altitude. The main concentra- tion of buildings within our study site was also found within Period Urban sprawl Forest loss the region with higher moisture variability. (m2 ) (m2 ) 1985–1990 76 301.57 27 977.24 3.5 Interferometric synthetic aperture radar time 1990–1995 22 042.68 18 651.49 series 1995–2000 12 716.93 25 433.86 2000–2005 40 694.17 33 911.81 InSAR time-series analysis revealed a narrow zone on the 2005–2010 23 738.27 17 803.70 eastern slope of the study site with a continuous motion away 2010–2015 9325.75 22 890.47 from the satellite at a rate of 2–3 mm yr−1 (Fig. 9c). Given the 2015–2020 10 173.54 18 651.49 line-of-sight (LOS) direction, the motion can be either west- ward or downward. Since westward motion is nearly impos- sible on an east-facing slope, we interpreted the motion as In Petrópolis municipality, forest cover has historically a continuous down-slope movement with part of the down- been replaced mostly by pasture and agriculture (3.38 %), ward motion (negative LOS velocity) canceled by the east- followed by urbanization (0.77 %). Looking closer to the ward motion (positive LOS velocity). study area, forest loss was constant over the entire study in- This interpretation means that the actual down-slope ve- terval, reaching about 16.6 ha by 2020 (Fig. 10 and Table 4). locity needs to be larger than the LOS velocity. Considering When the urban expansion data for the study area are ob- the slope angle (about 50–60◦ ; see Fig. 6) and the satellite served (Fig. 10 and Table 3), it is verified that this study area angle (about 38◦ ), the actual down-slope velocity within this reflects the data presented for the municipality, with an in- zone can reach as high as 10–15 mm yr−1 . The rest of the crease in the urban area over the periods and a higher urban study site showed little ground deformation between May expansion for the period of 1985 to 1990. The urban sprawl 2015 and December 2021. can be seen in Fig. 10, which presents the transition informa- 3.6 Urban sprawling, forest losses and landslide scars tion between classes for the study area. For the study area, approximately 1700 buildings that oc- During the last 35 years, forest formation in the munic- cupy the area around the hill were counted. When the build- ipality of Petrópolis has been increasing (7.08 %: from ings were confronted with landslide scars, it was found that 410.55 to 439.64 km2 ) and urban areas have largely ex- at least 60 houses were affected by landslides. Due to the panded (58.78 %: from 32.07 to 50.93 km2 ), while pasture characteristics of the relief in the region, occupation on hill- and agriculture have lost cover (14.73 %: from 336.70 to side areas is very common. Most of them (1021 buildings) 287.09 km2 ). However, when analyzing during a 5-year inter- are on the slope range of 20 % to 45 % (Fig. 11). It was ob- val, forest cover mostly increased from 1985–2010 (Table 2: served that 343 buildings occupy areas with slopes above there was a small decrease from 1995–2000). Since 2010, 45 %, which are considered “Permanent Preservation Ar- though, there has been a loss of 1.50 % in vegetation cover eas”, according to Brazilian Environmental Legislation (Law (6.72 km2 ). Pasture and agriculture consistently decreased no. 12.651/2012; Brasil, 2012). Considering the use and oc- during the study interval (exceptions were 1995–2000 and cupation of the study area, further stability analysis is rec- 2010–2015; Table 3). Regarding urban areas, 1985 to 1990 ommended to verify the possibility of variation in the safety was the period with the highest increase (about 22 %), and factor of the slopes due to the execution of slopes of cuts and the urban expansion has stayed at a constant rate of 2 %–7 % fills (Mendes et al., 2018a, b). since then, accruing a total increase of 58 % in 35 years. The landslide area coinciding with the urban area totaled about 17 000 m2 (Fig. 12a). The scars occurred in areas of Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023 https://doi.org/10.5194/nhess-23-1157-2023
E. Alcântara et al.: Deadly disasters in southeastern South America 1167 Figure 8. Soil moisture index from 28 January 2021 to 23 January 2022. Figure 9. (a) Optical image before February 2022 landslides (image from Google Earth Pro). (b) Digital elevation map (ALOS PALSAR). (c) InSAR line-of-sight (LOS) velocity between 8 May 2015 and 20 December 2021. Note that most areas have nearly zero deformation except for the eastern slope (pointed by a white arrow). (d) Standard deviation of LOS velocity. forested area, urban area, pasture, and agriculture, common in agriculture and pasture. This finding may be an indica- classes in the urban region of the municipality. From the tion that urban sprawl at the beginning of the analysis period zonal statistics considering the landslide scars of this event (1985) is sufficient for landslides to cause destruction, ma- and the LULC class as the majority in its interior, it was ver- terial damage, and loss of life in the region. The graph in ified that there were comparatively fewer landslide scars in Fig. 12b shows the zonal statistic result for the 17 landslide urban areas in 1985. Over the years, more scars were found scars (2022-year event) and LULC classes for the years of with urban area predominance, and fewer scars were found analysis. https://doi.org/10.5194/nhess-23-1157-2023 Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023
1168 E. Alcântara et al.: Deadly disasters in southeastern South America Figure 10. Transition classes for the study area: (a) 1985–1990, (b) 1990–1995, (c) 1995–2000, (d) 2000–2005, (e) 2005–2010, (f) 2010– 2015, (g) 2015–2020, and (h) area (m2 ) of urban sprawl and forest loss over the analyzed period. Figure 11. Building footprints and slope in the area of interest (a) and (b) number of buildings per slope interval. 4 Discussion Rio de Janeiro and created a favorable environment for at- mospheric convection. A similar analysis on meteorological 4.1 Origin/cause of the extreme rainfall, subsequent evolution leading to landslides can be found in Martinotti et flood and landslide al. (2017). The authors studied a period of intense rainfall from 1 to 6 September 2014 in the Gargano Promontory, a The weather forecasts issued by INMET and CEMADEN for karst area in Puglia, southern Italy; heavy rainfall is known the mountainous region of Rio de Janeiro, released earlier to cause multiple geohydrological hazards, including inun- on 14 February, warned about isolated convective rainfall, dations, flash floods, landslides, and sinkholes. which could occur in some areas of the city. However, no The other key element that led to the extraordinary rains in meteorological model predicted such significant amounts of the center of Petrópolis was the passage of a cold front, with rainfall over the region. The heavy rainfall that occurred in very particular characteristics, which occurred at the exact the city of Petrópolis on 15 February was caused by the ac- moment the rain showers began to form over the city. This tion of a meteorological phenomenon known as mesoscale cold front, on the one hand, was weak to dissipate the in- convective cell, with extraordinary characteristics without stability related to the storm clouds, and, on the other hand, known recorded antecedents. The situation was influenced it was strong enough to change the wind direction, which by the presence of the South Atlantic Convergence Zone came from the south, exactly perpendicular to the moun- (SACZ), which at the time was positioned over the state of Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023 https://doi.org/10.5194/nhess-23-1157-2023
E. Alcântara et al.: Deadly disasters in southeastern South America 1169 Figure 12. (a) Landslide scars superimposed in the urban area classes (the background image is from Planet) and (b) zonal statistics based on LULC classes (for every 5 years, from 1985 to 2020). tainous region of Petrópolis. As a result of this combination of factors, the mesoscale convective cell cloud (technically called of cumulonimbus), which should have lasted a few minutes, lasted several hours due to the interaction of the southerly winds associated with the cold front in the moun- tain. Figure 13 shows the surface winds and temperature for the 18:00 UTC on 15 February. It is noticed that there is a cold front with a temperature gradient and changes in wind direction. Southerly winds on the mountain region where the urban part of the city of Petrópolis (region prone to land- slides) turned the “orographic cloud” (cloud that positions on the top of mountains for hours and that normally do not precipitate) into a convective cell, which is very rare. Due to the sudden formation and the null displacement of the storm, the weather radars also did not allow its anticipated track- ing. It is noteworthy that the current state of knowledge and Figure 13. Global Forecast System (GFS) analysis at 18:00 UTC meteorological forecast do not allow predicting where each on 15 February 2022, for the Petrópolis region. The shaded areas individual cloud will form, with which this event could not represent surface temperature, the isolines depict sea level pressure, be predicted in advance. and the wind speed is indicated by barbs (in knots). Source: http: Figure 14 (lower side) shows the radar images of Pico do //atmosferalivre.com/ (last access: 5 May 2022). Couto site, where the formation of the cloud can be seen ex- actly above the center of the city of Petrópolis. It should be noted that only the residential area of the municipality was quent landslides in Petrópolis are translational of layers, with affected by the rain, which lasted more than 3 h. The accumu- varying thickness from 0.5 to 2 m. The municipality’s hy- lated rainfall over the Petrópolis station (Fig. 14, upper side) drography indicates the convergence of rivers and streams, shows the most intense rain between 19:00 and 21:00 UTC. which are in an anthropized hydrographic basin, while its to- The highest record of 260 mm in just 4 h, which occurred be- pographic characteristics resulted in high speed and energy tween the afternoon and evening of 15 February, is unprece- surface runoff. Due to these characteristics and the meteoro- dented in the city. The curvature of the storm produced by the logical event that hit the municipality, there was an increase southerly winds (represented by blue arrows in the figure) is in the levels of rivers and streams during the intense and con- also noted, which resulted in its long persistence. centrated rains. The drainage systems were overloaded, and, According to CEMADEN, from a hydrological point of as a result, the rainwater runs off the surface, causing flash view, the events recorded on 15 February 2022 in the city of floods and floods. Petrópolis were characterized as landslide, flood, and flash Regarding the occurrence of landslides, the geotechni- flood typologies. Flash floods in Petrópolis occur mostly af- cal stations of CEMADEN installed in the municipality of ter heavy rainfall associated with unplanned land and poor Petrópolis indicate that there was a significant increase in uses of land. According to Avelar et al. (2013) the most fre- soil moisture during the rainfall event that occurred between https://doi.org/10.5194/nhess-23-1157-2023 Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023
1170 E. Alcântara et al.: Deadly disasters in southeastern South America Figure 14. Hourly rainfall for various CEMADEN’s weather station in the city of Petrópolis during 13–16 February 2022 (a); radar images at the Pico do Couto site, taken at 19:00 UTC (b); 19:30 UTC (c); 20:00 UTC (d); and 20:30 UTC for 15 February 2022 (source: CEMADEN). 15:30 and 19:00 UTC on 15 February 2022. The monitoring study has shown (Rosi et al., 2019). Petrópolis is in a moun- station also indicated that prior to the event, soil moisture tainous region with lower deforestation, compared to less was already high since the previous rains exceeded 220 mm hilly sites (Silva et al., 2017), and natural regeneration may per 14 d and 350 mm per 21 d (recorded in some of the be happening there (especially replacing old, abandoned pas- CEMADEN stations). The preceding rise in soil moisture tures), as observed in other parts of the southeastern Atlantic was an inducing, preparatory factor for the occurrence of Forest (Silva et al., 2017). In addition, the municipality is part landslides. The abrupt elevation of soil moisture in a short of a protected area (sustainable use: APA Federal da Região time interval, due to intense and concentrated precipitation Serrana de Petrópolis), which may help forest conservation and, consequently, the oversaturation of the ground, provided and restoration with its management actions. However, the the triggering of mass movement processes. On the other forest located in places of interest, within the urbanized re- hand, it was also verified that very high soil moisture was gion, has substantially decreased in extent (Rosi et al., 2019). not recorded at deep levels (sensors in depths 2.0, 2.5, and At the same time, irregular and unplanned urbanization has 3.0 m). In other words, a priori, the processes were closely been increasing (Guerra, 1995). Thus, public policies, such related to intense surface runoff and percolation with very as a master plan, are needed to prevent land use changes that high positive pressure in fractures. The high fracture density might consequently lead to disasters. favors the formation of large blocks of rocks on the slopes. In Petrópolis, urbanization has been intensifying since the Blocks of rocks and colluvial soil deposits were incorporated 1930s, increasing the rural-to-urban migration and unequal into the mass of debris and deposited in the valley bottom as urban development. The population increased rapidly, reach- a poorly sorted material. ing 75 000 inhabitants in 1940 and 255 000 in 1990. The growth of tourism and the development of a land market 4.2 Tackle urban sprawling for disaster risk for second homes for citizens from Rio de Janeiro increased management the real estate speculation. This speculation excluded low- income households from the formal land market, driving Other than rainfall being the contributing factor, land use them to occupy hazardous zones. The rate and form of ur- change likely worsened the effect of the landslide event. banization expansion overcame precipitation as the most im- While we verified an increase in forest cover and urban area portant driver of landslides. Slope areas were unsuitably in- and a decrease in agriculture and pasture along the study in- corporated into the urban network. In the 1980s, the number terval in the municipality of Petrópolis, the forest located of mass movements was larger than those of the 1960s and in urbanized region has substantially decreased, as another Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023 https://doi.org/10.5194/nhess-23-1157-2023
E. Alcântara et al.: Deadly disasters in southeastern South America 1171 1970s, in contrast to rainfall values which were lower than environmental education, community meetings, and preven- the previous decade (Guerra, 1995). Our study revealed that tion advice and are as important as the structural ones. Early most of the buildings are located in areas of high slope (above warnings of disasters are essential to significantly reduce 20◦ ), which is not suitable for human settlements. As such, death toll and injuries. For instance, the 20 March event had better management is required to ensure that people inhabit no more than 4 % of the death toll compared to the 15 Febru- safer areas. ary event, despite the meteorological drivers of the heavy Despite not being an official data source, the Open- rainfall in March being of high magnitude. A detailed study StreetMap platform – OSM – is able to generate geographic is currently underway to go into depth to explain the differ- information about human settlements with high quality (Al- ence between the two events and the significance of early buquerque et al., 2016). With the need for rapid disaster re- warnings in saving lives. sponse, the crowdsourced data have made it possible to esti- In Petrópolis, there are ongoing preventive works, such as mate the impact caused by human settlements in areas with- the NUPDEC Vale do Cuiabá: a group of people who live out cadastral data available. Examples include the Wenchuan in risk areas and are in contact with the local civil defense (2008) and Haiti (2010) earthquakes (2008) and Cyclone to have coordinated actions during a disaster. Complement- Idai and Cyclone Kenneth (2019) in Mozambique (Li et al., ing the early warnings of disasters, issued by federal, state, 2020). and local authorities, Petrópolis also counts on sound warn- Dias et al. (2018) created a methodology to associate de- ings (sirens). The National Centre for Monitoring and Early mographic census data with disaster-prone areas with risk Warnings of Natural Disasters (CEMADEN) is responsible of landslides and floods, from 1999 to 2012. They found for issuing the warnings, which are forwarded to the federal that, out of the total population exposed to risk in Petrópo- civil defense, which distributes the early warnings to state lis, 8.19 % were children, 11.24 % were elderly, 48.20 % and municipal civil defense agencies. The local civil defense male, and 51.80 % female. In addition, an estimated 26 % of is responsible for evacuating risk areas and addressing rescue the at-risk population was living in subnormal agglomerates activities and shelter establishment. Concerning the civil de- (slums), from which 54 % had no water supply and 14 % had fense activities, the main needs, according to a recent survey, inadequate sanitation services (Dias et al., 2018). This im- are related to financial support, structure, and capacity build- portant information about the conditions of use and occupa- ing. Another recurrent theme among the civil defense pro- tion of the studied area reinforces the potential contribution fessionals was the lack of acknowledgement of their work. of anthropic inducing factors (leaks in pipes, for example) The actions to improve prevention and response to disasters, in the deflagration of landslides, and its detailed investiga- therefore, should necessarily consider the improvement of tion may even be possible from stability modeling coupled these mentioned deficiencies. with transient flow analyses (Mendes et al., 2018a, b). Cur- rently there are no data for 2022 using the same methodol- ogy, nor updated data from IBGE – the Brazilian Institute for 5 Conclusion Geography and Statistics. We retrieved some numbers from the local press and from the municipality’s official commu- Contrary to what the media reported, the seriousness of the nications. Beyond the impacts on human lives and injuries, Petrópolis disaster under analysis in the current study was not the economic consequences are also serious, with estimates exclusively due to the large volume of rain. The large death around BRL 78 000 000 (FIRJAN, 2022). toll had to do with the lack of a precise early warning due Marchezini and Wisner (2017) described a series of dy- to the unpredictable nature of the meteorological event that namic pressures when analyzing disaster risks in Petrópo- hit the region and due to the large number of areas of risk. lis: societal deficiencies (lack of planning and investments), A set of drivers were responsible for the deadliest disaster in business cycles, dense urbanization with population change, Petrópolis: the heavy rainfall, combined with the already sat- deforestation, poor governance, weakening of environmen- urated soil, increased unplanned urban sprawl, replacement tal legislation, real estate speculation, and land use change. of vegetation for surfaces with lower capacity of infiltration, They also described root causes (social and economic struc- and the lack of early warnings, was responsible for the disas- tures, history and culture heritage and ideologies) and factors ter that occurred on 15 February 2022. related to unsafe locations, such as sewage leaking into soil, By utilizing available precipitation and human settlement limited skills and formal education, marginalized groups and datasets, as well as employing multiple remote sensing data individuals, lack of access to formal credit, and lack of dis- in our analysis, including optical and radar images, we aster preparedness. showed that urban sprawling could have a significant ef- Disasters have serious consequences, like deaths, injuries, fect on slope stability and consequently landslide suscepti- and economic losses. Preventive actions to tackle risks and bility. Several spatial patterns of the 15 February 2022 land- disasters may be structural and non-structural. Structural as- slide event were identified. Other than the average rainfall pects refer to drainage, slope contention, urban services, and for February 2022 being the heaviest recorded since 1932, engineering works in general. Non-structural actions refer to the urban part of the city was particularly affected by heavy https://doi.org/10.5194/nhess-23-1157-2023 Nat. Hazards Earth Syst. Sci., 23, 1157–1175, 2023
You can also read