Assessing the relationship between weather conditions and rockfall using terrestrial laser scanner to improve risk management

Since 1987, more than 13,200 rockfalls have been inventoried by the Québec Ministry of Transport (MTQ) as having impacted the National Road 132 in northern Gaspésie. This natural hazard represents a nearly permanent danger for users. Traditional mitigation measures can be ineffective on poorly consolidated, deformed and highly fractured rockwall such as those found in northern Gaspésie. To address this issue, implementing a preventive risk management based on the factors that 10 trigger rock instabilities could be the most effective method. Earthquake, rainfall and freeze-thaw cycles are commonly considered to be the main rockfall triggering factors. This study aims to better understand the climatic conditions conducive to rockfalls in northern Gaspésie in order to provide knowledge to implement an appropriate risk management strategy. Three rockwalls were scanned with a terrestrial laser scanner (TLS) during specific pre-targeted weather conditions. Over a period of 18 months, 17 surveys have allowed to identify 1,287 rockfalls with a magnitude above 0.005 m³ on a scanned surface of 15 12,056 m². In addition, meteorological instruments and a 550 cm thermistor string have been installed directly on a vertical rockwall. It appears that some weather conditions influence occurrence, frequency, and magnitude of rockfalls. In winter, rockfall frequency is 12 times higher during a superficial thaw than during a cold period in which temperature remains below 0°C. In summer, rockfall frequency is 22 times higher during a heavy rainfall event than during a period mainly dry. Superficial freeze-thaw cycle (< 50 cm) causes mostly a high frequency of small magnitude events while deeper spring thaw (> 100 cm) 20 results in a high frequency of large magnitude events. Influence of weather conditions on rockfall frequency and magnitude is crucial in order to improve risk management since large magnitude events represent higher potential hazards. This study provides a classification of weather conditions based on their ability to trigger rockfalls of different magnitudes. This knowledge could be used to implement a risk management strategy


Introduction 25
Rockfall is hillslope movement in which blocks detach from the surface of rocky escarpments (Budetta, 2004;Michoud et al., 2012;Piteau and Peckover, 1978;Selby, 1993). Rock mass properties, such as lithology, degree of alteration, discontinuity network characteristics and slope, are predisposing factors in the occurrence, magnitude and failure mode (e.g. planar, wedge or toppling) of rock instabilities (e.g. Selby, 1993;Turner and Schustler, 1996). Although they are spontaneous, rockfalls result from the long term interaction of a series of processes Schovanec, 2020). They never result solely 30 in Haute-Gaspésie, implementing hazard mitigation methods based on the factors that trigger rock instabilities could therefore be the most effective method of limiting hazards associated with rockfall (Laliberté et al., 2022). 65 Recent literature illustrates significant advances in our understanding of the impact of weather variables on rockfall (D'Amato et al., 2016;Delonca et al., 2014;Macciotta et al., 2015;Matsuoka, 2019;Pratt et al., 2019;Ravanel and Deline, 2011;Weidner and Walton, 2021). The primary objective of our study is to strengthen this knowledge for a geological context neglected by the literature and fill in its main gaps: 70

•
LiDAR data makes it possible to accurately quantify rock instabilities, but they are not always linked to suitable weather data. For example, freeze-thaw cycles are often derived from air temperature without consideration of the influence of solar radiation, even on south-facing rockwalls (e.g. Weidner and Walton, 2021). Furthermore, air temperature provides no insight on freezing front depth in the rock. This makes it difficult to relate spring thaw to the rockfalls that occur (e.g. Macciotta et al., 2015). 75

•
The relationship between rock instabilities and weather conditions is often studied at the monthly or seasonal scale (e.g. Kromer et al., 2018;Macciotta et al., 2017). While this temporality can reveal interesting trends, it cannot differentiate the respective impacts of different meteorological events on rockfall dynamics. For example, on a monthly scale, the respective roles of precipitation events and freeze-thaw cycles cannot be easily distinguished since their occurrences can overlap. Consequently, this temporality is not useful in the preventive hazard management of 80 rock instability triggers.

•
Many studies highlight the strong relationship between rockfall frequency and magnitude (e.g. Guerin et al., 2014Guerin et al., , 2020Rosser et al., 2005;Santana et al., 2012;van Veen et al., 2017;Williams et al., 2018). However, studies into the influence of weather variables on rockwall dynamics focus almost exclusively on rockfall frequency. One exception is Matsuoka (2019), who investigated the relationship between weather conditions and rockfall magnitude. 85 In order to mitigate the risk, the relationship between the magnitude of rock instabilities and their triggering factors is an essential area of study since the damage caused by instabilities is proportional to their magnitude.
Our study therefore aims to link rockfall dynamics to specific adequately measured weather conditions. We use TLS to quantify the frequency and magnitude of rock instabilities for each targeted weather condition. We also propose a classification of 90 weather conditions based on their level of concern for preventive hazard management.

Study sites
In Haute-Gaspésie, Route 132 is the sole transportation corridor linking a number of villages to the rest of the region. For the local population, it is critical infrastructure for accessing essential services. The national road runs along the Gulf of St.
Lawrence shoreline and is subject to coastal flooding and erosion (e.g. Drejza et al., 2014). The rockwalls that tower above the 95 road expose it to other natural hazards: snow avalanches (Fortin et al., 2011;Hétu, 2007), falling ice blocks (Gauthier et al., 2017), debris flows (Fortin et al., 2015), and rockfall. Those same rockwalls make moving the roadway inland from the shore unfeasible. To limit risk to road users, the Ministère des Transports du Québec (department of transportation) created a 24/7 patrol to monitor the 70 km of national road that run along the rockwall of Gaspé's uplands. The patrol is also responsible for clearing the road when it is obstructed by rock, snow or ice from hillslope movements. Between 1987 and 2020, the patrol 100 responded to 13,261 rockfall events along the 25 km where rockwalls overhang Route 132 or 16 rockfallsyear -1 km -1 (Ministère du Transport du Québec, 2021).
LiDAR was used to monitor three rockwalls (Figure 1, A1 and A2) that overhang Route 132 in Haute-Gaspésie. From east to west, the sites are near the villages of Manche-d'Épée (MAE), Gros-Morne (GMO) and Marsoui (MAR). The total area 105 scanned was 12,056 m 2 . Rockwalls were selected using five criteria: 1) structural and lithologic characteristics representative of Haute-Gaspésie rockwalls; 2) frequent rockfall incursions onto the national road; 3) year-round accessibility; 4) general morphology (verticality, height) and microtopography (roughness) compatible with LiDAR acquisition from the national road; and 5) absence of vegetation (occlusion) and drainage (signal reflection). The three sites are sedimentary rockwalls (flysch) with horizontal or subhorizontal bedding planes. The easternmost study site, MAE, covers an area of 3,154 m 2 . It has a vertical 110 (>80°) natural rockwall that is 35 m high and is oriented north (350°). It rises over an unvegetated 30 m long scree slope with an average gradient of 38°. The rockwall is primarily composed of siltstone (50%), shale (30%) and greywacke (20%). Rockfall of greywacke blocks is most likely to reach the national road, which runs about 10 m from the base of the scree slope. GMO is a rockwall with alternating spurs and re-entrants. On that site, three rock spurs with a total area of 5,370 m 2 were scanned, most of which have been dynamited. The orientation of the spurs ranges from 330° to 20°, and their gradients range between 115 70° and 90°. Their bases do not have scree slopes, and rockfall regularly reaches the national road. Only the lower 25 m of the three spurs were scanned to limit occlusion as much as possible. The rock spurs are composed of greywacke (50%), siltstone (40%) and shale (10%). Lastly, MAR is a vertical (90°) dynamited rockwall composed primarily of sandstone (70%), siltstone (20%) and shale (10%). The thickness of the siltstone strata tends to result in decimeter sized instabilities. Since the national road runs right along the base of the rockwall (no scree slope is present), these instabilities regularly reach the roadway. 120 However, its low height (<30 m) and verticality limit the distance of rockfall travel from the rockwall.
Haute-Gaspésie has a humid continental climate with short cool summers, according to the Köppen climate classification system (Beck et al., 2018). During the period 1991-2020, it had an annual mean temperature of 3.9 °C The mean temperatures of the warmest (July) and coldest (January) months were 16.3 °C and -9.2 °C, respectively. Over those same 30 years, average 125 annual precipitation was 888 mm, 33% of which fell as snow (Environment Canada, 2021). that measured the temperature to a depth of 550 cm (white circle). The red line is for the exposed to rockfalls Route 132.

Rockfall detection using a TLS
Over the 554-day period from 14 June 2019 to 18 December 2020, 17 LiDAR surveys were performed. Two TLS were used for data acquisition. The four first surveys used a Leica Geosystems ScanStation C10. A RIEGL VZ-400i was used for the 13 135 subsequent surveys. To obtain dense point clouds with the fewest possible occlusions, multiple surveys were conducted to adequately cover the scanned surfaces at each site (Figure 2a and 2b). Depending on site morphology and surface area, four to twelve 360° surveys were needed. The resulting point clouds were processed in the Cyclone 2.5 software. For each date, all these point clouds have been assembled together by a visual alignment using the "2D Scan Thumbnails" tool. Then, the point clouds have been optimized using the "point clouds autocorrelation" tool. This method allows to minimize the number of 140 occlusions and to produce point clouds with a centimeter-level density for each survey (Figure 2c and 2d).The same strategy was used to align the point clouds corresponding to the different dates of acquisition. For each of these dates, we have measured the surface differences between the cloud points using the M3C2 plugin of the open source free software CloudCompare v2.12.4 (Lague et al., 2013). The M3C2 algorithm operates directly on point clouds without meshing or gridding and computes the local distance between cloud points along the normal surface direction which tracks 3D variations in surface orientation. 145 Because it considers the surface roughness, this algorithm is particularly accurate to measure small surface changes and so, it is suitable to compare rockwall surfaces (Lague et al., 2013).
Once the surface differences between each survey have been calculated, we quantified the rockfall volumes between these periods ( Figure 2e). Each significative negative change was isolated from the whole point clouds and then, the 2.5D volume 150 calculation method was applied to compute its volume. Štroner et al. (2019) have compared this method, also available in CloudCompare, with software solutions using point clouds capable of calculating volume (Atlas DMT, 3D Reshaper, Leica Cyclone and Trimble RealWorks). They have proven than the 2.5D volume calculation method provides results identical to those computed by the best commercial solutions and has a much better time demand/accuracy ratio (Štroner et al., 2019).Finally, to avoid misinterpretation of rockfall volumes with, for example, wind deflation of snow or the presence of 155 vegetation, we have focused the analyses on vegetation-free cliff sections. A visual examination of the shape of each of the quantified volumes also helped to minimise interpretation errors (van Veen et al., 2017).

Protocol to monitor the influence of weather conditions on rockfalls 160
To study the influence of weather conditions on rockfall frequency and magnitude, sufficiently high spatial resolution is needed to identify all events at the lowest quantified magnitude (Abellán et al., 2014). Nevertheless, small-scale rockfalls can be difficult to distinguish from other processes such as wind deflation of snow over the rock slope surfaces. In order to avoid misinterpretation, a threshold of 10 cm for negative change has been determined and the scarce areas where vegetation is present were excluded of the studied areas. At the end, this method allows exhaustive detection of rockfalls larger than 165 0.005 m³.
Regarding the temporal resolution, Barlow et al. (2012) showed that a 19 month interval between surveys leads to an underestimation of small scale events. With overly long intervals, several distinct but spatially overlapping events can be erroneously interpreted as a single larger scale event. Williams et al. (2018) compared rockfall frequencies detected at temporal 170 resolutions of one hour and of 30 days. They found that the number of low magnitude events (<0.1 m 3 ) detected was three orders of magnitude greater at a temporal resolution of one hour. However, to achieve this high temporal resolution, a fixed scanner is required. This leads to other issues, such as greater occlusion on the rockwall surface and more limited scanned areas. Barlow et al. (2012) showed that while environmental factors influence rockfall, defining a constant frequency is not necessary for data acquisition. LiDAR surveys must be more frequent during periods when a greater rock instability frequency 175 is expected. Conversely, intervals between surveys can be longer during periods with a lower rockfall frequency. The protocol we implemented was designed to limit the scanning effort by targeting selected periods rather than performing regular scans.
This resulted in an inconsistent survey frequency that is justified by the influence of weather conditions on rockfalls. The scanning periods were selected based on knowledge of the weather conditions that cause major rock surface deformations and are likely to result in rockfall . Target periods could run over several dozen days (dry period, rainy 180 period, autumn period with near-daily freeze-thaw cycles and cold winter period) or monitor specific meteorological events (heavy and high intensity rainfall, winter freeze-thaw cycles and spring surface, 1 m deep and full rockwall thaws).
The periods to scan were determined over the course of the study using Environment Canada weather forecasts and then validated with local weather measurements. A Campbell TE525WS-L rain gauge was used to measure precipitation (mm). In 185 winter, it was equipped with a CS705 adapter to measure the water equivalent of solid precipitation. A thermistor string (with a GeoPrecision data logger operating at 915 MHz) inserted into a horizontal borehole measured the temperature (°C) every 30 cm to a depth of 550 cm. Using linear interpolation between the thermistors, we produced a temperature profile at different depths over time, quantified the number and depth of freeze-thaw cycles and tracked the depth of the freezing front in winter and the thawing front in spring. All weather instruments took data readings every 15 minutes. The instrumented rockwall is 190 located 200 m west of the study site GMO (49°15'21.0"N 65°33'52.2"W) and the precipitation data comes from a weather station located in a roadside rest area along the Gulf of the St-Lawrence (49°13'49.92"N 65°51'3.53"O) ( Figure 1).
Finally, the rockfall spatial distribution as well as its frequency and the rock slopes erosion rates are presented individually for the three study sites but the influence of weather conditions is based on the whole rockfall database. This approach has proved 195 necessary to avoid misinterpretation of the occurrence or non-occurrence of events during short microclimatic periods (e.g. heavy and high intensity rainfall or winter freeze-thaw cycles). Monitoring over a longer period would make it possible to establish links between the meteorological conditions and the geological and structural context of each of the studied rock slopes.

Hazard assessment 200
LiDAR surveys make it possible to quantify rockfall frequency and to calculate their individual volume. Based on this data, we can study the influence of different weather scenarios on rockfall magnitude. Because a high frequency of high magnitude events does not have the same significance as a high frequency of low magnitude events, the magnitude of rockfalls is a major issue for risk management. Following this precept, we have grouped and classified the 9 main weather scenarios in a matrix according to their ability to trigger rockfalls of different magnitudes. This matrix is composed of three magnitude classes 205 ("low", "intermediate" and "high") and five frequency classes from "low" to "extreme". This table provides a classification of weather conditions based on their ability to trigger rockfalls of different magnitudes. This knowledge could be used to implement a risk management strategy based on triggering factors.
In this case, it can be defined as: where f(V) is the mean daily frequency of rockfalls with a volume greater than or equal to V over the entire study period for the 12,056 m² of scanned rockwalls. This frequency is the equivalent of 2.2 rockfallsd -1 for events greater than or equal to 0.01 m³, 0.2 rockfallsd -1 for events greater than or equal to 0.1 m³, and 0.02 rockfallsd -1 for events greater than or equal to 1.0 m³. 69 rockfalls have a volume higher than 0.15 m³, (i.e. 5.3% of the total) of which 17 have higher than 0.5m 3 (i.e. 1.3% of the total) and 7 higher than 1 m 3 (i.e. 0.5% of the total) ( Figure 3). Where the frequency or spatial resolution of the LiDAR 220 surveys is insufficient, low magnitude instabilities are underrepresented (Barlow et al., 2012;Guerin et al., 2014;Malamud et al., 2004). The absence of a rollover in the logarithmic frequency-magnitude distribution obtained in this study makes it possible to validate whether a suitable instability detection threshold was chosen ( Figure 3b). This result confirms that the selected survey frequency and spatial resolution were appropriate for individually identifying the large majority of events with a volume exceeding 0.005 m³. 225

Global rockwall erosion rate 230
The annual rockwall erosion rates and the rockfall frequency were calculated using the entire 554-day time series of two summers and two autumns but only one winter and one spring. The 1,287 recorded rockfalls were distributed across the three study sites as follows: 363 at MAE, 486 at GMO and 438 at MAR. Figure 4, A1 and A2 show the surface differential at the three study sites between the first survey in June 2019 and the last survey in December 2020. On these figures, we have filtered out any change less than +/-1 cm and green clusters of points correspond to positive change and red clusters of points 235 correspond to negative changes. Over the 554-day period, many changes were visible on the rockwall surfaces. Most of the positive changes (in green) to the rockwall and scree slope are explained by low snow accumulation during the last LiDAR surveys on 18 December 2020. These positive changes can also be attributed to the development of rock instabilities large enough to be detected by LiDAR (Kromer et al., 2018;Oppikofer et al., 2008Oppikofer et al., , 2009Royán et al., 2014). For example, measurements taken at MAE show that between the months of June 2019 and December 2020, a rock mass measuring 3 m 240 high separated from the rockwall by 4 to 5 cm at its top and 1 to 2 cm at its base ( Figure 4a). The increased rate of displacement at the top of the rock mass suggests a toppling failure. This type of pre-failure deformation was not considered in the rockfall volume calculation. Only rockfalls that produced a negative change of more than 0.005 m 3 in the surface differential were considered in the erosion rate calculation. Because rockfalls of very small magnitude (< 0.005 m 3 ) were not considered, the erosion rates are slightly underestimated. It is important to recall that what appears to be a high magnitude event over the entire 245 study period may correspond to several successive smaller-scale events (Figure 4b) (Abellán et al., 2014;Barlow et al., 2012;van Veen et al., 2017;Williams et al., 2018). The change detection of +/-1 cm applied in Figure 4, A1 and A2 was only used to visualize the overall picture of the rock slope dynamic (e.g. snow accumulation, pre-failure deformation) but as described before, a filter of +/-10 cm was applied to quantify rockfalls. Considering the scanned surfaces separately, all three sites had fairly similar annual rockfall frequencies: 0.08 rockfallsm -² 255 for MAE and MAR, and 0.06 rockfallsm -² for GMO (Figure 5a). In contrast, the differences between the erosion rates of the three rockwalls were more pronounced. Over the study period, the erosion rate was higher for MAR (5.4 mmy -1 ) than for MAE (3.0 mmy -1 ) or GMO (2.8 mmy -1 ) ( Figure 5b). Therefore, despite comparable rockfall frequencies, MAR had a higher erosion rate than the other sites. This distinction shows that on average, MAR experienced higher magnitude events. This is borne out by the fact that in the entire time series, 5 of the 7 recorded events with a volume of over 1 m 3 occurred at MAR. 260 These events account for 33% of the total volume lost by the three rockwalls in the study period. At MAR, the eight largest magnitude events contributed as much to the erosion rate of the studied surface as the 430 lower magnitude events. Meanwhile, at GMO and MAE, the 5 and 20 highest magnitude events represented the volume of the 481 and 343 smallest events, respectively. Clearly, high magnitude events affect more substantially the erosion rate of the rockwalls (Figure 5c).

Weather conditions related to rockfalls 270
In this study, LiDAR surveys made it possible to quantify instabilities over 16 distinct periods ( Figure 6). Three of them totaling 104 days with a mean precipitation of 1.4 mmd -1 were defined as dry periods (DRY). Two periods totaling 100 days and a mean precipitation of 4.1 mm d -1 were defined as rainy periods (RAINY). Two periods totaling 93 d with a mean air temperature below the freezing point (-6.2 °C) and very low winter freeze-thaw frequency (0.02 thawsd -1 ) were defined as cold winter periods (WIN. <0 °C). The autumn periods (AUT. FT) experienced the first freeze-thaw cycles after summer. Two 275 autumn periods totaling 79 days were scanned. The specific meteorological events that were monitored included a heavy (54 mm) and high intensity (4 mm h -1 for 12 h 45 min) rainfall event in September 2019 (HEAVY RF) and two winter freezethaw cycles (WIN. FT) in winter 2020. One of these two winter freeze-thaw cycles was accompanied by rainfall. Since this kind of liquid precipitation event is only possible because of the thaw occurrence, we made the choice to consider this period as a WIN. FT. Four LiDAR surveys were used to segment the spring thaw (SPR. MELT) into three periods (Figure 7).. The 280 first phase corresponds to the initiation of the thawing season including two freeze-thaw events reaching a depth of 50 cm from the rockwall surface (SPR.0:50). The second phase is characterized by a thaw that penetrates to a depth of 95 cm and by a thaw of 77 cm from depth (SPR.50:100). The third phase corresponds to the complete thaw of the rockwall (SPR.100:350) ( Figure 7). Lastly, the period between 25 May and 16 September 2020 experienced a range of weather conditions too varied for categorization ( Figure 6). 285 In decreasing order, 241 rockfalls were associated with RAINY periods, 178 with AUT. FT periods, 145 with DRY periods, 120 with SPR. MELT periods, 97 with WIN. <0 °C periods and 91 with the HEAVY RF period and 44 with WIN. FT periods 295 ( Figure 8a). However, rockfall occurrence under these meteorological conditions is not representative of the dynamics of rockwall instability since some of those conditions occur over a much longer period than others. For example, DRY periods total 104 days, while WIN. FT periods total just four. Considering rockfall frequency rather than occurrence revealed different meteorological conditions as triggering factors of instabilities. WIN. <0 °C and DRY were associated with the lowest daily rockfall frequencies, respectively 0.9 and 1.4 events d -1 across the full scanned surface area (12,056 m²). This frequency tended 300 to double under AUT. FT (2.2 eventsd -1 ), RAINY (2.4 eventsd -1 ) and SPR. MELT (2.6 eventsd -1 ) conditions. Winter freezethaw cycles (WIN. FT) and heavy rainfall events (HEAVY RF) produced the highest rockfall frequencies. Frequencies for WIN. FT and HEAVY RF were 11.0 and 30.3 events d -1 , respectively (Figure 8b). Compared to a cold winter period (WIN. <0 °C), rockfall frequency was multiplied by a factor of 1.5 during DRY periods, 2.3 in AUT. FT periods, 2.8 in SPR. MELT periods, 11.7 in WIN. FT periods and 32.3 in HEAVY RF periods. 305

310
For each scanned survey period, precipitation intensity and freeze-thaw cycles were compared with rockwall erosion rates and rockfall frequency (Figure 9). The erosion rate and frequency shown for each period are relative to their respective mean values over the entire time series (standard score). Where values for a period were under one, the rockwall erosion rate and rockfall frequency were less than the time series mean. On 7 September 2019, Category 5 hurricane Dorian arrived in the Gulf of St.
Lawrence as a post-tropical storm. It dropped 54 mm of heavy rain (4 mmh -1 for 12 h 45 min) on the study site. Between 1987 315 and 2021, only 16 precipitation events of over 50 mm were recorded in Haute-Gaspésie, which represents a return period of 776 days (Environment Canada, 2021). The heavy rain event resulted in a rockwall erosion rate and rockfall frequencies that were 9.3 and 13.8 times the mean, respectively (Figure 9). Calculated for the full 72 h period between the two surveys (7 and 10 September 2019), those values would undoubtedly have been higher had they had been based solely on the actual duration of the event (14 h). 320 Between 23 and 26 February 2020, a series of three freeze-thaw cycles occurred (Figure 9). They followed a long period (70 days) of below freezing point temperatures. During these three freeze-thaw cycles, the air temperature rose to 3.1 °C, 8.3 °C and 2.1 °C, respectively, separated by 8 h periods around -2 °C. They resulted in a thaw 10 cm depth from rockwall surface.
The rockwall erosion rate and rockfall frequency were 6.3 and 8.5 above the mean, respectively. The absence of rain during 325 this period highlights the leading role of surface thaws in rockslide dynamics (Figure 9). The air temperature remained below the freezing point until 20 March 2020 before rising above zero for 24 hours, briefly reaching a high of 8.1 °C. This produced a surface thaw to a depth 14 cm accompanied by rockwall retreat and rockfall frequency rates of 2.3 and 4.4 times the mean, respectively (Figure 9). In the spring, the frequency and magnitude of freeze-thaw cycles increased (Figure 7 and Figure 9).
Rockfall frequency remained fairly low and lower than during the last winter freeze-thaw event (1.7 times more than the mean). 330 In contrast, erosion rates were higher than during the freeze-thaw event on 23 and 26 February and remained high for a period of 24 days (2.5 times more than the mean). This combination of a high erosion rate and fairly low rockfall frequency indicates an increased mean rockfall volume during the spring thaw.
Over the study period, the weather conditions that most contributed to rock instability dynamics were the winter freeze-thaw 335 cycles, the spring thaw and the heavy and high intensity rainfall event (Figure 9). In contrast, other weather conditions were associated with high stability in the studied rockwalls. In the period 5 December 2019-23 February 2020, both the rockwall erosion rate and rockfall frequency were very low, 0.6 and 0.3 times less than the mean, respectively (Figure 9). Temperatures remained below freezing point on all but two days of this 81-day period. The 19-day period 6-25 May 2020 had very low rockwall erosion rates and rockfall frequency (0.2 and 0.6 times less than the mean, respectively). During that period, rockwall 340 temperature remained constantly below freezing and only 5.6 mm of rain was recorded. In the study period, the cold winter and dry summer periods experienced the weather conditions associated with the highest stability in the scanned rockwalls ( Figure 9).

Hazard assessment 350
LiDAR surveys make it possible to quantify rockfall frequency and calculate the individual volume of each rockfall. Using this data, we can study the influence of different weather conditions on rockfall magnitude. Figure 10a presents a matrix of the weather conditions associated with different rockfall frequencies and magnitudes. We defined three classes of "low", "intermediate" and "high" magnitudes. The first one includes 1,100 low magnitude (0.005:0.05 m³) events. The second class comprises 119 events of intermediate magnitude (0.05:0.15 m³), while the third includes 68 events of the highest magnitude 355 (≥ 0.15 m³). We also defined five frequency classes for the matrix: low (<75% mean frequency for the study period), moderate (75:125%), high (125:300%), very high (300:500%) and extreme (> 500%). In that year, HEAVY RF and SPR.100:350 periods accounted for just 11 days. There was therefore a low occurrence of periods with weather conditions favorable to extreme frequency high magnitude events. The WIN. FT, SPR.0:50, SPR.50:100 and AUT. FT periods were also associated with high frequency or high magnitude rockfall and represent a cumulative occurrence of 44 days between 14 June 2019 and 13 June 2020 (Figure 10b). Lastly, the DRY, WIN. <0 °C and RAINY periods, which 375 do not coincide with periods of high instability, accounted for 310 days, or 85% of the year (Figure 10b). Thus, on an annual scale, the weather conditions that resulted in high frequency or magnitude rockfall are in a small minority.

Flysch rockwall erosion rate
The annual erosion rates measured over 554 days for the three flysch rockwalls in Haute-Gaspésie are very high (between 2.8 385 and 5.4 mmy -1 ) compared to reported rates for various lithologies and climate patterns (e.g. André, 1997;Barsch, 1977;Beylich, 2000;Curry and Morris, 2004;Hinchliffe and Ballantyne, 1999;Höllerman, 1983;Humlum, 2000;Olyphant, 1983;Sass, 1998). They are nonetheless perfectly consistent with those reported by Hétu and Gray (2000) for a 40 m flysch rockwall in the same study area (3.2 mmy -1 ). Those rates were measured by weighing the sediment accumulated on plastic tarps placed at the foot of rockwalls (Hétu and Gray, 2000). Direct measurements by TLS are extremely accurate, making it possible to 390 quantify the erosion rates of rockwalls in their morphoclimatic environments with a greater degree of certainty. The agreement between the results from the two methods demonstrates that they are both valid. However, surveys over very short periods clearly tend to underestimate actual rates by limiting the likelihood of observing very high magnitude events that would have a significant influence on erosion rate (e.g. Korup and Clague, 2009;Krautblatter and Dikau, 2007). Continuing to survey these rockwalls would confirm the accuracy of the erosion rates presented here. TLS is a modern method that merits more 395 widespread use in monitoring projects of this type and to build on our understanding of rockwall retreat and erosional feature dynamics (e.g. Ballantyne and Kirkbride, 1987;Duszyński et al., 2019;Higgins and Coates, 1990;Rapp, 1960).

Relationship between meteorological conditions and rockfall
Rockfall frequency during RAINY periods was 1.7 times higher than the mean relative to DRY periods. The highest magnitude and strongest intensity rain event (HEAVY RF) was associated with rockfall frequency 9.3 times higher than the study period 400 mean, or 21.6 times higher than during a dry period. Only RAINY conditions led to an increased frequency of low magnitude instabilities, while HEAVY RF conditions resulted in a significant increase in the frequency of high magnitude instabilities ( Figure 10). These results agree with those obtained in other geological and climatic contexts. For the limestone rockwalls of Isère, France, D' Amato et al. (2016) showed that rockfall frequency was 2.5 to 3 times higher during rain events than during the dry summer period, up to 7.5 higher during high magnitude rain (> 30 mm) and 27 higher during high intensity rain (> 405 5 mm h -1 ). For the rockwalls of Réunion Island, Delonca et al. (2014) showed that high magnitude rain (> 120 mm) resulted in an eightfold increase in the probability of daily rockfall. The lack of a time lag between rain events and instabilities shows that this precipitation has a nearly immediate (<24 h) triggering effect on rock instabilities. In the Japanese Alps, Matsuoka (2019) demonstrated the impact of liquid precipitation on increased water content and subsequent rockfall triggering. Certain high intensity rain events were associated with instabilities of one to multiple meters in size. These studies, like our own, 410 clearly show that rain events are aggravating and triggering factors of rockfall. However, they do not detail the process by which rainwater contributes to the formation of instabilities or improve our understanding of how water moves through the highly fractured rock layer near rockwall surfaces (Magnin et al., 2021;Stoll et al., 2020). It is well known that an increase in hydrostatic pressure in discontinuities reduces shear strength (Selby, 1993;Wieczorek and Jäger, 1996;Wyllie and Mah, 2004).
Groundwater recharge and hydrostatic pressure at depth have been put forward as kinematic factors in large rockslides 415 (Cloutier et al., 2015;Crosta et al., 2013). However, how rainwater influences the development of superficial instabilities is less well documented. In clay-rich rock, it has been shown that water content variations at the rockwall surface lead to irreversible mechanical deformations ) that could trigger rockfall. Rockwall surfaces can experience wide fluctuations in water content. Rain events result in subsurface runoff in the altered surface layer of rockwalls and represent a significant water input. This input can be amplified when wind spreads a layer of rain water over rockwall surfaces (Sass, 420 2005). This type of subsurface runoff can clearly contribute to reducing shear strength and trigger rockfall during high intensity rain events (Selby, 1993). Conversely, direct solar radiation on rockwall surfaces can lead to rapid and significant drying of the first decimeters of rock (Burnett et al., 2008). Limited access to rockwalls and the challenge of directly measuring water pressure limit our ability to measure moisture dynamics in this surface layer of rock and interpret its influence on the development of rock instabilities. 425 In this study, rockfall frequency was 11.7 times higher during winter freeze-thaw cycles (WIN. FT) than during cold periods when the temperature remains below the freezing point (WIN. <0 °C). We therefore observe that high rockfall frequencies occur during thawing periods (WIN. FT), while rockwalls remain very stable during the phases of the freezing process (WIN. <0 °C). The freezing phase can cause fractures to open (Coutard and Francou, 1989;Matsuoka, 2008;Matsuoka and Sakai, 430 1999). However, cohesion at the ice-rock interface is generally sufficient to hold newly unstable blocks of rock in place (Fiorio et al., 2002;Krautblatter et al., 2013). Consequently, rockfalls occur when the ice melts, such as during winter freeze-thaw cycles. This finding supports those of D' Amato et al. (2016), who showed that in limestone rockwalls in Isère, France, the frequency of rock instabilities is seven times higher during freeze-thaw cycles. Similarly, Delonca et al. (2014) found that rockfall frequency was doubled on basalt and granite rockwalls in Auvergne, France, for freeze-thaw cycles with a minimum 435 temperature of -5 °C. Our results also show that freeze-thaw cycles are associated with a very high frequency of low and moderate magnitude events but that their impact is limited when it comes to higher magnitude events (Figure 10a). This is consistent with the fact that winter freeze-thaw cycles influenced only the first 15 cm from surface (Figure 9a). Winter thaws can be accompanied by rainfall ( Figure 9) and a compounding effect of these conditions probably occurs when they are combined. Autumn freeze-thaw cycles (AUT. FT) do not appear to be as effective at developing rock instabilities as winter 440 freeze-thaw cycles. Their associated frequency of rockfall is near the annual mean for both low and high magnitude instabilities ( Figure 10a). This is likely partly due to the fact that autumn freeze-thaw cycles were not studied individually. The two AUT.
FT periods were 28 and 51 days long, respectively. They also coincided with weather conditions unmarked by high rockfall frequency that decreased the average frequency of rockfall specific to autumn freeze-thaw cycles. Matsuoka (2019) also notes that even if the temperature fluctuations around the freezing point penetrate to 40 cm in depth, the water in the pores and 445 discontinuities does not necessarily freezeor freeze completelydue to its high heat capacity. Lastly, the high intensity autumn rain events that preceded the period of many autumn freeze-thaw cycles (AUT. FT) may have purged the rockwall surface of its most unstable rock, leaving very little rock in a nearly unbalanced state.
The superficial spring thawing phase (SPR.0:50) does not appear to be a particularly unstable period, even for low magnitude 450 events ( Figure 10a). However, the thaw in the first meter of rockwall (SPR.50:100) results in very high frequency of moderate magnitude instabilities. When thawing at depth (SPR.100:350) occurs, the frequency of high magnitude instabilities increases ( Figure 10a). This relationship between seasonal freezing depth and rockfall magnitude is widely recognized (e.g. Dramis et al., 1995;Matsuoka, 2019;Matsuoka and Murton, 2008), but for the first time, this study compares the magnitude of rock instabilities with measurements of the spring thawing front at depth. The winter freeze-thaw cycles (WIN. FT) that occurred 455 before the first thawing phase (SPR.0:50) may have been responsible for purging the first decimeters of rockwall surfaces and explain why SPR.0:50 conditions are not associated with an unstable period. This surface thaw phase lasted 23 days. During that period, negative temperatures continued to dominate (conditions equivalent to WIN. <0 °C) and probably tended to attenuate rockfall frequency throughout the period. At greater depth, the spring thaw is the only thaw of the year (Figure 7 and Figure 9). Consequently, the mechanical action of this high amplitude freeze-thaw cycle has strong potential to destabilize 460 rock portions (Dramis et al., 1995;Matsuoka, 2019;Matsuoka and Murton, 2008). The impact of the spring thaw on high magnitude rockfall is amplified by water advection associated with snow cover melting, groundwater recharge and ice melting in rock discontinuities and pores. As spring progresses, rain precipitation becomes predominant, further increasing water inputs (Hasler et al., 2011). Large magnitude rockfall is likely to result from the reduced cohesion associated with ice melt and increases in hydrostatic and interstitial pressure in the altered layer of rockwall surfaces (Matsuoka, 2019;Schovanec, 2020;465 Wieczorek and Jäger, 1996).
Among the contributing factors to rock instability development, the roles of precipitation and freeze-thaw cycles are predominant (e.g. Collins and Stock, 2016;Coutard and Francou, 1989;D'Amato et al., 2016;Hungr et al., 1999;Matsuoka and Sakai, 1999;Rapp, 1960;Wieczorek and Jäger, 1996). To a lesser degree, other weather phenomena can trigger rockfall. 470 Large temperature fluctuations cause expansion and contraction within rock masses that can eventually result in rockfall (Collins and Stock, 2016;Eppes et al., 2016;Matsuoka, 2019). The impact of thermal cycles on the rockfall dynamics of the highly altered north-facing rockwalls of Haute-Gaspésie was considered to be minor and was not quantified in this study.

Rockfall hazard management
The magnitude of rockfalls is a major issue in risk management. Large dimension blocks tend to travel longer distances (e.g. 475 Dorren, 2003). If they reach part of the natural or built environment that is exposed, they have a larger probability of causing damage. Therefore, in risk management, a high frequency of high magnitude events does not have the same significance as a high frequency of low magnitude events. The rockfall on 16 April 2019 that occurred 1.3 km west of the MAE site resulted in the temporary closure of Route 132 ( Figure 11). It took place 18 hours after the first thaw to 60 cm in depth and coincided with the spring thaw (SPR.50:100). It also occurred immediately after a 20.1 mm rainfall event (HEAVY RF). In retrospect, 480 the SPR.50:100 and HEAVY RF conditions meant that extremely high frequency high magnitude rockfall was likely ( Figure 10). From a risk management standpoint, forecasting a major event such as the one on 16 April 2019 is possible using the risk assessment matrix developed in this study. High intensity rain (HEAVY RF) and spring thaw SPR.100:350 were identified as the weather conditions of greatest concern for public safety since they result in extreme frequency high magnitude events (Figure 10a). In our first full year of data collection (14 June 2019-13 June 2020), these conditions only occurred on 485 11 days (Figure 10b). The other weather conditions that lead to increased rockfall frequency and magnitude (WIN. FT, SPR.0:50, SPR.50:100 and AUT. FT) add 44 more days of concern. From the perspective of rockfall forecasting, it is noteworthy that potentially hazardous conditions are present only 55 days annually for users of Route 132. In the year under study, heightened vigilance was needed only 15% of the days on an annual scale (Figure 10b) that has wide seasonal disparities ( Figure 9 and Figure 10). The Ministère des Transports du Québec could implement different risk mitigation measures targeted 490 for the weather conditions forecast by the weather models. While RAINY conditions do not stand out in this study for causing high frequency and/or magnitude rock instabilities, rain events should be given special attention. Our methodology did not make it possible to clearly encompass all rain events on an annual scale. Given that HEAVY RF conditions are associated with extreme frequency high magnitude events, moderate intensity rain events must also feed back into frequencies and magnitudes of concern for the safety of Route 132 users. The same reasoning applies to the AUT. FT period. Had it been possible to better monitor fall freeze-thaw cycles, they would likely 500 have shown a frequency of associated rock instabilities more similar to that of winter freeze-thaw cycles (WIN. FT).
This study highlights that rockfall dynamic is strongly controlled by the weather conditions and by the thermal regime of the rock slopes. Some of these conditions are associated with rockwall stability while others significantly increase rockfall frequency and/or magnitude (Figure 8, Figure 9 and Figure 10). Because of the inconsistency of weather conditions in term of 505 duration (e.g. a cold winter period versus a winter thaw), we propose to focus on an event-based monitoring approach rather than a period-based approach to survey rockfalls. This approach allows to differentiate the respective impacts of different meteorological events that cannot be easily distinguished on a periodic scale since their occurrences overlap. For an equivalent number of surveys and therefore for the same effort, this event-based monitoring approach is more relevant to identify weather conditions conducive to rockfalls of different magnitudes. 510

Conclusion
The influence of weather conditions on the occurrence of rock instabilities is now widely recognized in the scientific community. This study proposed to quantify rockfall volumes using LiDAR during previously selected weather conditions rather than taking regular surveys (weekly, monthly, seasonal). For an equivalent number of surveys and therefore for the same effort, this protocol made it possible to efficiently monitor the influence of weather conditions on rockfall frequency and 515 magnitude. Over 18 months, 17 surveys were carried out on three rockwalls with a total area of 12,056 m 2 . They made it possible to monitor a range of weather conditions and identify 1,287 rockfalls. The rockwall erosion rate during that period was 2.8-5.4 mmy -1 , which is among the highest in the world but comparable to other observations in this type of lithology (e.g. André, 1997;Barsch, 1977;Beylich, 2000;Curry and Morris, 2004;Hinchliffe and Ballantyne, 1999;Höllerman, 1983;Humlum, 2000;Olyphant, 1983;Sass, 1998). 520 Some of the weather conditions are associated with periods of rockwall stability (dry summer and cold winter periods) while others significantly increase rockfall frequency and/or magnitude (winter freeze-thaw cycles, spring thaw, high intensity rainfall). Winter freeze-thaw cycles and moderate rainfall lead to a increase in the frequency of low magnitude rock instabilities. High intensity rainfall and spring thaw coincide with a significant increase in the frequency of large dimension 525 rockfall. In the realm of rockfall prevention, this study showed the importance of quantifying rockfall magnitude since the relationship between weather conditions and rockfall frequency is not necessarily applicable to rockfall magnitude. The results allow weather conditions to be categorized by their ability to trigger rock instabilities of varying magnitudes. Using this combined with current or forecast weather conditions, different risk mitigation measures to limit safety hazards can be implemented.

Data and code availability 540
All the data collected during this study are available on request. Codes were only used to realize the Figure 7 and can also be send on request.

Author contribution
Tom Birien made the field work, developed the data processing, carried out the data analysis and wrote the manuscript. Francis Gauthier helped to develop the initial idea and to prepare the manuscript. 545

Competing interests
The authors declare that they have no conflict of interest.

Acknowledgments
We thank the ministère des Transports du Québec (MTQ) as well the Natural Sciences and Engineering Research Council of Canada (NSERC) for funding this research. We acknowledge the LGGRM members who help us with the instrumentation 550 process of the study sites, the maintenance, and data collection. We also would like to thank the anonymous referees of Natural Hazards and Earth System Sciences for valuable comments that improved the quality of the paper.  Geomorphology, 139-140, 416-424, doi:10.1016Geomorphology, 139-140, 416-424, doi:10. /j.geomorph.2011Geomorphology, 139-140, 416-424, doi:10. .11.006, 2012