This repository consists of NCL and MATLAB code for processing, analyzing, and plotting LOcalized Constructed Analogs (LOCA; Pierce et al. 2014) statistically downscaled, gridded climate data. The gridded climate data includes daily maximum surface air temperature (tasmax), surface air minimum temperature (tasmin), and daily precipitation (pr) at 1/16th degree resolution over the Northeastern United States domain for 29 Coupled Model Intercomparison Project 5 (CMIP5) models. Simulations include historical (1980-2005), lower emissions rcp4.5 (2006-2099), and higher emissions rcp8.5 (2006-2099).
This is a repository in progress and updates can be expected throughout Summer 2018. Code will not run "as-is", but will require customization (i.e., paths, sub-folders).
Below is the list of scripts and general order of usage in my own processing.
ann_hist.ncl
- Calculates the historical (1980-2005) climatological average & linear trend in annual temperature (min & max) andprecipitation (total).
- Inputs include daily LOCA data for three variables:
- tasmin (surface minimum air temperature)
- tasmax (surface maximum air temperature)
- precipitation (daily total)
- Output includes network Common Data Format (netCDF, or .nc) files for climatological mean and annual trend in eachvariable for each model:
- ann_pr_historical_1980-2005.nc
- ann_tasmax_historical_1980-2005.nc
- ann_tasmin_historical_1980-2005.nc
ann_future.ncl
- calculates the future projected (2006-2099) climatological average annual temperature (min & max) andprecipitation (total) for early (2010-2039), mid (2040-2069), and late (2070-2099) century.
- Inputs include daily LOCA data for three variables:
- tasmin (surface minimum air temperature)
- tasmax (surface maximum air temperature)
- precipitation (daily total)
- Output includes network Common Data Format (netCDF, or .nc) files for annual climatological mean and trend for eachscenario (rcp$ = rcp45 or rcp85) and within the file for each variable for each model:
- ann_pr_rcp$_2006-2099.nc
- ann_tasmax_rcp$_2006-2099.nc
- ann_tasmin_rcp$_2006-2099.nc
seas_hist.ncl
- calculates the historical (1980-2005) climatological average and linear trend in seasonal temperature (min & max) andprecipitation (total).
- Seasons include:
- Spring (March, April, May; MAM),
- Summer (June, July, and August; JJA)
- Fall (September, October, and November; SON)
- Winter (December, January, February; DJF). Winter years are associated with the January of a given year (i.e., Winter2018 is December 2017, January and Feburary 2018).
- Inputs include daily LOCA data for three variables:
- tasmin (surface minimum air temperature)
- tasmax (surface maximum air temperature)
- precipitation (daily total)
- Output includes network Common Data Format (netCDF, or .nc) files for climatological mean and seasonal ($SEAS$ = MAM, JJA,SON, or DJF) trend in each variable for each model:
- $SEAS$_pr_historical_1980-2005.nc
- $SEAS$_tasmax_historical_1980-2005.nc
- $SEAS$_tasmin_historical_1980-2005.nc
seas_future.ncl
- calculates the future projected (2006-2099) climatological average seasonal temperature (min & max) andprecipitation (total) for early (2010-2039), mid (2040-2069), and late (2070-2099) century.
- Seasons include:
- Spring (March, April, May; MAM),
- Summer (June, July, and August; JJA)
- Fall (September, October, and November; SON)
- Winter (December, January, February; DJF). Winter years are associated with the January of a given year (i.e., Winter2018 is December 2017, January and Feburary 2018).
- Inputs include daily LOCA data for three variables:
- tasmin (surface minimum air temperature)
- tasmax (surface maximum air temperature)
- precipitation (daily total)
- Output includes network Common Data Format (netCDF, or .nc) files for climatological mean and seasonal ($SEAS$ = MAM, JJA,SON, or DJF) trend for each scenario (rcp$ = rcp45 or rcp85) and within the file for each variable for each model:
- $SEAS$_pr_rcp$_2006-2099.nc
- $SEAS$_tasmax_rcp$_2006-2099.nc
- $SEAS$_tasmin_rcp$_2006-2099.nc
ExtremeTemp_hist.ncl
- calculates the historical (1980-2005) climatological average and linear trend in extreme temperature metrics:
- days90F = number of days per year when daily tasmax is greater than 90F (32.2C)
- days95F = number of days per year when daily tasmax is greater than 95F (35C)
- days32F = number of days per year when daily tasmin is less than 32F (0C)
- days0F = number of days per year when daily tasmin is less than 0F (-17.8C)
- txHottest = tasmax on the hottest day of the year
- txColdest = tasmax on coldest day of the year
- tnColdest = tasmin on coldest night of the year
- tnHottest = tasmin on hottest night of the year
- Inputs include daily LOCA data for two variables:
- tasmin (surface minimum air temperature)
- tasmax (surface maximum air temperature)
- Output includes a single network Common Data Format (netCDF, or .nc) file for climatological mean and seasonal trend ineach extreme temperature metric for each model:
ExtremeTemp_future.ncl
- calculates the historical (2006-2099) climatological average for early (2010-2039), mid (2040-2069), and late (2070-2099)century in extreme temperature metrics:
- days90F = number of days per year when daily tasmax is greater than 90F (32.2C)
- days95F = number of days per year when daily tasmax is greater than 95F (35C)
- days32F = number of days per year when daily tasmin is less than 32F (0C)
- days0F = number of days per year when daily tasmin is less than 0F (-17.8C)
- txHottest = tasmax on the hottest day of the year
- txColdest = tasmax on coldest day of the year
- tnColdest = tasmin on coldest night of the year
- tnHottest = tasmin on hottest night of the year
- Inputs include daily LOCA data for two variables:
- tasmin (surface minimum air temperature)
- tasmax (surface maximum air temperature)
- Output includes network Common Data Format (netCDF, or .nc) files for climatological mean and seasonal trend ineach extreme temperature metric for each model:
- ExtremeTemp_rcp45_2006-2099.nc
- ExtremeTemp_rcp85_2006_2099.nc
Time series files are named similarily to the .ncl scripts above except with *_ts.ncl appeneded to the end. An example is given for ann_hist_ts.ncl.
ann_hist_ts.ncl
- Calculates the historical (1980-2005) time series in annual temperature (min & max) and precipitation (total).
- Inputs include daily LOCA data for three variables:
- tasmin (surface minimum air temperature)
- tasmax (surface maximum air temperature)
- precipitation (daily total)
- Output includes network Common Data Format (netCDF, or .nc) files of the annual time series for eachvariable and for each model:
- ann_pr_historical_1980-2005_timeseries.nc
- ann_tasmax_historical_1980-2005_timeseries.nc
- ann_tasmin_historical_1980-2005_timeseries.nc
References:
Pierce, DW, DR Cayan, and B Thrasher. 2014. Statistical Downscaling Using Localized Constructed Analogs (LOCA). J Hydromet. 15: 2558-2585. doi: 10.1175/JHM-D-14-0082.1.