DOI: 10.5281/zenodo.17469867
title: Winter warming from climate change endangers the global hotspot of salamander diversity

description: Data and code supporting: Winter warming from climate change endangers the global hotspot of salamander diversity. 

Important READ ME instructions:

Code descriptions:

Dormancy_inactive_cody.py is a Python script that calculates spatial patterns of lipid consumption during winter dormancy across the species’ range. The code was written and tested in Python 3.9.14, and we only verified functionality with this version; we therefore recommend running it with Python 3.9.x to avoid compatibility issues.

Installing Python

To run the script, you first need a working Python installation:

Go to the official Python website:https://www.python.org/downloads/

Download the installer for your operating system (Windows, macOS, or Linux) for Python 3.9 (or the closest available 3.9.x release).

Run the installer and, on Windows, ensure the option “Add Python to PATH” is checked.

Complete the installation following the on-screen instructions. On most modern computers, this is a straightforward process and finishes quickly.

Alternatively, you may install Python via a scientific distribution such as Anaconda (https://www.anaconda.com) and create a dedicated environment with Python 3.9.

Required input data

The script requires gridded climate data in ASCII raster (.asc) format as described in the Methods. Use the climate .asc files provided in the associated data repository and ensure that:

The .asc files are downloaded and stored locally (e.g. in a data/ or climate/ directory).

The file paths in Dormancy_inactive_cody.py point to the correct location of these .asc files (edit the input paths near the top of the script if necessary).

How to run the script

You can run the script either from a command-line interface (Terminal, PowerShell) or from an IDE such as Spyder:

Option 1 – Terminal / command line

Open a Terminal (macOS/Linux) or Command Prompt / PowerShell (Windows).

Navigate to the directory containing Dormancy_inactive_cody.py, for example:

cd /path/to/your/code

Run the script with Python 3.9:

python Dormancy_inactive_cody.py

(If you have multiple Python versions installed, you may need python3 instead of python.)

Option 2 – Spyder or similar IDE

Open Spyder (e.g. from an Anaconda installation).

Open the file Dormancy_inactive_cody.py in the editor.

Make sure the working directory in Spyder is set to the folder containing the script and input .asc files.

Click Run to execute the script.

Expected output and runtime

When executed successfully, Dormancy_inactive_cody.py produces ASCII raster (.asc) files estimating the amount of lipids consumed across the range map over winter dormancy. These output rasters are written to the output directory specified in the script (see comments in the code for the output path). On a typical workstation, the full run generally completes in under one hour, depending on hardware and the spatial extent and resolution of the input climate data.

Reproducing figures

Instructions for reproducing the figures in the manuscript from the model outputs are provided in the Methods section of the main text (and/or the accompanying analysis scripts, as referenced there).

lipids_analysis_summary_data.py: This script processes ASCII raster maps representing the proportion of lipids lost during winter dormancy across the Plethodon jordani species complex. It imports climate-derived lipid depletion maps, reads the species’ range map, and calculates summary statistics for each scenario, including (1) the number and proportion of grid cells with fully depleted lipid stores, (2) the mean proportion of lipids lost across the range, and (3) the standard deviation of lipid depletion. The script outputs a CSV file (winter_dormancy_summary_stats.csv) containing these summary metrics. It was written in Python and uses standard scientific libraries (numpy, pandas, glob, and matplotlib).

lipid_analysis_ms.R: This script implements a logistic regression analysis to test whether the proportion of lipids consumed during winter dormancy predicts the probability of species presence across the landscape. The workflow includes data import, cleaning, model fitting, result export, and visualization of predicted relationships.vco2.R: This script conducts mixed-effects modeling and visualization of salamander metabolic rates (VCO₂) across temperature treatments, seasons, and elevations. It imports metabolic data, fits linear mixed models using lme4, performs post-hoc comparisons with emmeans, and generates predicted effects using ggeffects. The script produces publication-quality figures showing seasonal and elevational patterns in log-transformed VCO₂, as well as analyses of winter mass retention relative to metabolic rates. Output files include ANOVA tables (m1.csv, m2.csv) and several PNG figures saved to user-specified directories.