This set of files includes downscaled projections of decadal means of annual day of freeze or thaw (ordinal day of the year), and length of growing season (numbers of days, 0-365) for each decade from 2010 - 2100 at 2km x 2km meter spatial resolution. Each file represents a decadal mean of an annual mean calculated from mean monthly data. ---- The spatial extent includes Alaska, the Yukon Territory, British Columbia, Alberta, Saskatchewan, and Manitoba. Each set of files originates from one of five top ranked global circulation models from the CMIP5/AR5 models and RPCs, or is calculated as a 5 Model Average. Day of Freeze, Day of Thaw, Length of Growing Season calculations: Estimated ordinal days of freeze and thaw are calculated by assuming a linear change in temperature between consecutive months. Mean monthly temperatures are used to represent daily temperature on the 15th day of each month. When consecutive monthly midpoints have opposite sign temperatures, the day of transition (freeze or thaw) is the day between them on which temperature crosses zero degrees C. The length of growing season refers to the number of days between the days of thaw and freeze. This amounts to connecting temperature values (y-axis) for each month (x-axis) by line segments and solving for the x-intercepts. Calculating a day of freeze or thaw is simple. However, transitions may occur several times in a year, or not at all. The choice of transition points to use as the thaw and freeze dates which best represent realistic bounds on a growing season is more complex. Rather than iteratively looping over months one at a time, searching from January forward to determine thaw day and from December backward to determine freeze day, stopping as soon as a sign change between two months is identified, the algorithm looks at a snapshot of the signs of all twelve mean monthly temperatures at once, which enables identification of multiple discrete periods of positive and negative temperatures. As a result more realistic days of freeze and thaw and length of growing season can be calculated when there are idiosyncrasies in the data. Please note that these maps represent climatic estimates only. While we have based our work on scientifically accepted data and methods, uncertainty is always present . Uncertainty in model outputs tends to increase for more distant climatic estimates from present day for both historical summaries and future projections.

This set of files includes downscaled projections of monthly totals, and derived annual, seasonal, and decadal means of monthly average temperature (in degrees Celsius, no unit conversion necessary) from 1901 - 2006 (CRU TS 3.0) or 2009 (CRU TS 3.1) at 771 x 771 meter spatial resolution.

This set of files includes downscaled future projections of vapor pressure (units=hPa) at a 1km spatial scale. This data has been prepared as model input for the Integrated Ecosystem Model (IEM). There can be errors or serious limitations to the application of this data to other analyses. The data constitute the result of a downscaling procedure using 2 General Circulation Models (GCM) from the Coupled Model Intercomparison Project 5 (CMIP5) for RCP 8.5 scenario (2006-2100) monthly time series and Climatic Research Unit (CRU) TS2.0 (1961-1990,10 min spatial resolution) global climatology data. Please note that this data is used to fill in a gap in available data for the Integrated Ecosystem Model (IEM) and does not constitute a complete or precise measurement of this variable in all locations. RCPs: 8.5 Centers, Model Names, Versions, and Acronyms: National Center for Atmospheric Research,Community Earth System Model 4,NCAR-CCSM4 Meteorological Research Institute,Coupled General Circulation Model v3.0,MRI-CGCM3 Methods of creating downscaled relative humidity data: 1. The GCM input data are distributed as relative humidity along with the CRU CL 2.0, therefore no conversion procedure was necessary before beginning the downscaling procedure. 2. Proportional Anomalies generated using the 20c3m Historical relative humidity data 1961-1990 climatology and the projected relative humidity data (2006-2100). 3. These proportional anomalies are interpolated using a spline interpolation to a 10min resolution grid for downscaling with the CRU CL 2.0 Relative Humidity Data. 4. The GCM proportional anomalies are multiplied by month to the baseline CRU CL 2.0 10min relative humidity climatology for the period 1961-1990. Creating a downscaled relative humidity projected time series 2006-2100. 5. Due to the conversion procedure and the low quality of the input data to begin with, there were values that fell well outside of the range of acceptable relative humidity (meaning that there were values >100 percent), these values were re-set to a relative humidity of 95 at the suggestion of the researchers involved in the project. It is well known that the CRU data is spotty for Alaska and the Circumpolar North, due to a lack of weather stations and poor temporal coverage for those stations that exist. 6. The desired output resolution for the AIEM modeling project is 1km, so the newly created downscaled time series is resampled to this resolution using a standard bilinear interpolation resampling procedure. 7. The final step was to convert the downscaled relative humidity data to vapor pressure using the calculation below, which uses a downscaled temperature data set created utilizing the same downscaling procedure. EQUATION: saturated vapor pressure = 6.112 x exp(17.62 x temperature/(243.12+temperature)) vapor pressure = (relative humidity x saturated vapor pressure)/100

Atmospheric rivers (ARs) were detected from ERA5 6hr pressure level data, using a detection algorithm adapted from Guan & Waliser (2015). The algorithm uses a combination of vertically integrated water vapor transport (IVT), geometric shape, and directional criteria to define ARs. See the sources listed below and the GitHub repository for more detail and other references. The AR database is a zipped archive containing multiple attributed shapefiles. Polygon data includes individual timestep ARs, ARs making landfall in Alaska, and aggregated landfalling AR events. Point data includes coastal impact points landfalling AR events.

These files include historical downscaled estimates of decadal average monthly snow-day fraction ("fs", units = percent probability from 1 – 100) for each month of the decades from 1900-1909 to 2000-2009 at 771 x 771 m spatial resolution. Each file represents a decadal average monthly mean. Version 1.0 was completed in 2015 using CMIP3. Version 2.0 was completed in 2018 using CMIP5. For more information on the methodology used to create this dataset, and guidelines for appropriate usage of the dataset, please see the data user's guide here: http://data.snap.uaf.edu/data/Base/AK_771m/historical/CRU_TS/snow_day_fraction/snow_fraction_data_users_guide.pdf

This dataset consists of single band GeoTIFFs containing total annual counts of wet days for each year from 1980-2100 for one downscaled reanalysis (ERA-Interim, 1980-2015) and two downscaled CMIP5 global climate models driven under the RCP 8.5 baseline emissions scenario (NCAR-CCSM4 and GFDL-CM3, 2006-2100), all derived from the same dynamical downscaling effort using the Weather Research and Forecasting (WRF) model (Version 3.5). A day is counted as a "wet day" if the total precipitation for that day is 1 mm or greater.

This dataset consists of four different sub-datasets: degree days below 65°F (or "heating degree days"), degree days below 0°F, degree days below 32°F (or "freezing index"), and degree days above 32°F (or "thawing index"). All were derived from the same dataset of outputs from dynamically downscaling one reanalysis (ERA-Interim) and two CMIP5 GCMs (GFDL-CM3, NCAR-CCSM4) over Alaska using the Weather Research and Forecasting model (WRF). Data from the GCMs are driven exclusively by the RCP 8.5 emissions scenario. Heating degree days, degree days below 0°F, and freezing index were computed in the following way: subtract the daily mean temperature values from the threshold value and compute the sum of this time series for the given calendar year. Thawing index is instead computed as the annual sum of the quantities resulting from subtracting the threshold (32°F) from the daily mean temperature values.

This set of files includes downscaled historical estimates of monthly total precipitation (in millimeters, no unit conversion necessary, rainwater equivalent) from 1901 - 2013 (CRU TS 3.22) at 10 min x 10 min spatial resolution with global coverage. The downscaling process utilizes CRU CL v. 2.1 climatological datasets from 1961-1990.

This set of files includes downscaled historical estimates of decadal means of annual day of freeze or thaw (ordinal day of the year), and length of growing season (numbers of days, 0-365) for each decade from 1910 - 2006 (CRU TS 3.0) or 2009 (CRU TS 3.1) at 2x2 kilometer spatial resolution. Each file represents a decadal mean of an annual mean calculated from mean monthly data. **Day of freeze or thaw units are ordinal day 15-350 with the below special cases.** *Day of Freeze (DOF)* `0` = Primarily Frozen `365` = Rarely Freezes *Day of Thaw (DOT)* `0` = Rarely Freezes `365` = Primarily Frozen *Length of Growing Season (LOGS)* is simply the number of days between the DOT and DOF. ---- The spatial extent includes Alaska, the Yukon Territories, British Columbia, Alberta, Saskatchewan, and Manitoba. Each set of files originates from the Climatic Research Unit (CRU, http://www.cru.uea.ac.uk/) TS 3.0 or 3.1 dataset. TS 3.0 extends through December 2006 while 3.1 extends to December 2009. **Day of Freeze, Day of Thaw, Length of Growing Season calculations:** Estimated ordinal days of freeze and thaw are calculated by assuming a linear change in temperature between consecutive months. Mean monthly temperatures are used to represent daily temperature on the 15th day of each month. When consecutive monthly midpoints have opposite sign temperatures, the day of transition (freeze or thaw) is the day between them on which temperature crosses zero degrees C. The length of growing season refers to the number of days between the days of thaw and freeze. This amounts to connecting temperature values (y-axis) for each month (x-axis) by line segments and solving for the x-intercepts. Calculating a day of freeze or thaw is simple. However, transitions may occur several times in a year, or not at all. The choice of transition points to use as the thaw and freeze dates which best represent realistic bounds on a growing season is more complex. Rather than iteratively looping over months one at a time, searching from January forward to determine thaw day and from December backward to determine freeze day, stopping as soon as a sign change between two months is identified, the algorithm looks at a snapshot of the signs of all twelve mean monthly temperatures at once, which enables identification of multiple discrete periods of positive and negative temperatures. As a result more realistic days of freeze and thaw and length of growing season can be calculated when there are idiosyncrasies in the data.

This set of files includes downscaled projected estimates of monthly temperature (in degrees Celsius, no unit conversion necessary) from 2006-2300* at 15km x 15km spatial resolution. They include data for Alaska and Western Canada. Each set of files originates from one of five top ranked global circulation models from the CMIP5/AR5 models and RCPs, or is calculated as a 5 Model Average. *Some datasets from the five models used in modeling work by SNAP only have data going out to 2100. This metadata record serves to describe all of these models outputs for the full length of future time available. The downscaling process utilizes CRU CL v. 2.1 climatological datasets from 1961-1990 as the baseline for the Delta Downscaling method.