Climate warming drives rockfall from an increasingly unstable mountain slope

This readme file provides all data and R codes used to perform the analyses presented in Figs. 2-4 of the main text and Supplementary Information Figures S1-S2-S3.

FIGURE 2
-  Seasonally_dated_GDs.txt: Contains information on the timing (Season) of rockfall (GD) in a given tree (Id) and a given year (yr) over the past 100 years. Inv refers to the operators which analyzed growth disturbances in the tree-ring series. Lat / Long refers to the position of the tree in CH1903/ Swiss Grid projection. Intensity (1-4) refers to (1), intermediate (2)  and strong (3) GD. Intensity 4 was attributed to injuries (I). Only the 408 GD rated 3 (strong TRD) and 4 (injuries) were used in Fig. 2. Acronyms used for Response_type read as follows: TRD: Tangential rows of traumatic resin ducts; I: Injuries. Acronyms used for Season refer to Dormancy (1_D), early (2_EE), middle (3_ME) and late (4_LE) earlywood, whereas a GD found in the latewood was attributed to either the early (5_EL) or late (6_LL) latewood.

- Trends_in_seasonality_R1.R: The data contained in "Seasonally_dated_GDs" were processed with the R script "Trends_in_Seasonality.R". This seasonal trend analysis code is inspired by work published by Schlögl et al. (2021; https://doi.org/10.1016/j.crm.2021.100294) and Heiser et al. (2022; https://doi.org/10.1029/2011JF002262).

FIGURE 3-4-S1
- Tasch_GD.txt: Contains the raw data on rockfall impacts (GD) in a given year (yr) as found in all trees available in that same year (Sample_depth) as well as the cumulated diameter at breast height (cumulated_DBH) of all trees present in that same year.
- Rockfall_frequency_climate.R: The data contained in "Tasch_GD.txt" were processed with the R script "Rockfall_frequency_climate.R".  
- The temperature (Imfeld23_tmp.txt) and precipitation (Imfeld23_prc.txt) data used in Fig. 3 are from the Imfeld et al. 2023 (10.5194/cp-19-703-2023) gridded dataset (1x1 km lat/long) and were extracted at the grid point centered on the Täschgufer site.
- The script set with temperature series enables to compute Fig. 4 (l.149:216) and Fig. 3 (l. 216:330); the script set with precipitation series enables to compute Fig. S1

FIGURE S2
- Tasch_GD.txt: Contains the raw data on rockfall impacts (GD) at the Täschgufer site in a given year (yr) as found in all trees available in that same year (Sample_depth) as well as the cumulated diameter at breast height (cumulated_DBH) of all trees present in that same year.
- Rockfall_frequency_borehole.R: is adapted from "Rockfall_frequency_climate.R" to work with the borehole dates.  
- Corvatsch0_6R1: Contains the Corvatsch borehole temperature series (2000-2020, 0.6m depth) (Hoelzle, M. et al. https://doi.org/10.5194/essd-14-1531-2022, 2022).

FIGURE S3
- Plattje_GD.txt: Contains the raw data on rockfall impacts (GD) at the Plattje site in a given year (yr) as found all trees available in that same year (Sample_depth) as well as the cumulated diameter at breast height (cumulated_DBH) of all trees present in that same year.
- - Rockfall_frequency_climate_Plattje.R: The data contained in "Plattje_GD.txt" were processed with the R script "Rockfall_frequency_climate_Plattje.R".  
- The temperature (Imfeld23_tmp_Plattje.txt) and precipitation (Imfeld23_prc_Plattje.txt) data used in Fig. 3 are from Imfeld et al. 2023 (10.5194/cp-19-703-2023) gridded dataset (1x1 km lat/long) and were extracted at the grid point centered on the Plattje site.