Data and code from: Climatic variation allows montane willows to escape an adaptive tradeoff

This repository contains the data and R code needed to reproduce all results and figures from Rosenblad & Ackerly (2024) Climatic variation allows montane willows to escape an adaptive tradeoff. File descriptions are below.

analysis.R contains all R code.

cutting_data.csv contains the following fields: 1- plant_id, a unique identifier for each genetically distinct in situ parent plant, from which cuttings were collected and propagated; 2- cutting_id, which specifies replicate cuttings from each in situ parent plant; 3- diameter_mm, the initial diameter at the base of each cutting (in millimeters) immediately following collection; and 4- height_cm, the initial height of each cutting.

freezing_resistance_data.csv contains the following fields: 1- plant_id (see above); 2- cutting_id_1 (see cutting_id above); 3- cutting_id_2, which specifies replicate greenhouse-grown twigs sampled from each greenhouse plant (specified by cutting_id_1); and 4- new_growth, a binary variable indicating whether each twig survived and grew roots following the freezing treatment.

leaf_and_shoot_mass_data.csv contains the following fields: 1- plant_id (see above); 2- cutting_id (see above); 3- shoots_g, the dry mass of all leaves harvested from the greenhouse plant at the end of the growth experiment; and 4- leaves_g, the dry mass of all shoots harvested from the greenhouse plant at the end of the growth experiment.

plant_id.csv contains the following fields: 1- plant_id (see above); and 2- site_id, a unique identifier for the provenance site of each plant.

psi_tlp_data.csv contains the following fields: 1- plant_id (see above); 2- cutting_id (see above); 3- psi_tlp_mpa, the measured turgor loss point (in megapascals); and 4- collection_environment, a binary variable distinguishing greenhouse plants from in situ parent plants.

root_mass_data.csv contains the following fields: 1- plant_id (see above); 2- cutting_id (see above); and 3- root_mass_g, the dry mass of roots for the given measurement. Some plants had too much root material to measure all at once, so these plants were split into multiple measurement batches, each of which receives a row in this CSV.

S_lemmonii_GBIF.csv contains occurrence data of S. lemmonii from the Global Biodiversity Information Facility. There are many fields. The important fields for this use case are decimalLatitude and decimalLongitude, the WGS 84 coordinates (EPSG 4326) of each observation. Information on the other fields is available at gbif.org.

treatment.csv contains the following fields: 1- plant_id (see above); 2- cutting_id (see above); and 3-treatment, a binary variable indicating whether each greenhouse plant received the drought or control treatment during the growth experiment.

fig1.png and other similarly named .png files are the figures in the paper. They can be re-generated by running analysis.R.

aet_1981_2010_bcm.tif contains mean annual actual evapotranspiration data for the period 1981-2010 from the Basin Characterization Model.

aprpck1981_2010_ave_HST_1641848708.tif, and other files with the same basename but different extensions, contain mean April 1 snow water equivalent data for the period 1981-2010 from the Basin Characterization Model.

ds542.shp, and other files with the same basename but different extensions, contains a polygon outlining the Sierra Nevada ecoregion.

KR_sites_WGS84.shp, and other files with the same basename but different extensions, contain polygons of the provenance sites.

tmn1981_2010may_ave_HST_1689800754.tif, and other files with the same basename but different extensions, contain mean May minimum temperature data fro the period 1981-2010 from the Basin Characterization Model.