NAAPS-TimeSeries-2023

NAAPS-TimeSeries-2023 This is a repository containing data and code used in "Two decades of wildfire smoke deposition on the Greenland Ice Sheet: Spatiotemporal trends, the extreme 2023 season, and implications for surface albedo" (Hall et al., in preparation). Requirements Only R is needed for this analysis. The .csv files needed are in the Data folder. The analysis was performed using R-Markdown files, which are in the Scripts folder. Usage The data required for all of the analysis is in the Data folder. Download the data, then download the R-Markdown scripts in the Scripts folder to the same location. Then, create a new R project using the folder where the data and scripts were downloaded. The scripts will run with the data provided. DATA NAAPS Data Contains daily smoke deposition fluxes from the Navy Aerosol Analysis and Prediction System (NAAPS) reanalysis model from each transect in individual 1º x 1º tiles from 2003 - 2023. The AllTransects_XDeposition.csv files contain NAAPS daily smoke deposition flux data from all transects together NAAPS-AWS-Joined Contains daily NAAPS data and daily (averaged) measurements of albedo and other meteorological variables from PROMICE/GC-Net stations that intersect the given NAAPS tiles PROMICE_GCNet_Hourly_RawData Contains the hourly raw data for each PROMICE / GC-Net weather station used in analysis Data originally downloaded from: https://dataverse.geus.dk/dataset.xhtml?persistentId=doi:10.22008/FK2/IW73UU SCRIPTS HourlyAlbedoAveraging Use this script to aggregate the hourly PROMICE/GC-Net data to a daily resolution using a mid-day window for the most accurate albedo measurements AWS_NAAPS_Combining Use this script to combine the NAAPS and daily averaged AWS measurements into harmonized .csv files BaselineThresholding_LongTermTrends Use this script to define the ice-sheet-wide baseline deposition values for setting a deposition threshold for substantial deposition for each deposition mechanism (wet, dry, total) Then assess long-term trends in smoke deposition, smoke days, smoke events, and smoke event intensity Also used to compare smoke deposition in 2023 to the previous 20-year average DepositionTrends_Strucchange Use this script to assess trends in smoke deposition and identify the point in which there is a significant structural change in the slope of the deposition trend (2019) 2019_Trends_Pre_PostComparison Use this script for assessing smoke deposition trends before and after the identified breakpoint in 2019. ExtremeEvents Use this script to assess extreme smoke events (>= 95th percentile of all days over the baseline values) over time, pre/post-2019 and in 2023 compared to average FullTS_GAMs Use this script to run generalized additive models (GAMs) and relative importance analyses to compare the impacts of wildfire smoke deposition on albedo reduction relative to other variables measured at PROMICE/GC-Net stations