How Unusual Was 2015 in the 1984–2015 Period of the North Cascade Glacier Annual Mass Balance?

In 1983, the North Cascade Glacier Climate Project (NCGCP) began the annual monitoring of the mass balance on 10 glaciers throughout the range, in order to identify their response to climate change. Annual mass balance (Ba) measurements have continued on seven original glaciers, with an additional two glaciers being added in 1990. The measurements were discontinued on two glaciers that had disappeared and one was that had separated into several sections. This comparatively long record from nine glaciers in one region, using the same methods, offers some useful comparative data in order to place the impact of the regional climate warmth of 2015 in perspective. The mean annual balance of the NCGCP glaciers is reported to the World Glacier Monitoring Service (WGMS), with two glaciers, Columbia and Rainbow Glacier, being reference glaciers. The mean Ba of the NCGCP glaciers from 1984 to 2015, was −0.54 m w.e.a−1 (water equivalent per year), ranging from −0.44 to −0.67 m w.e.a−1 for individual glaciers. In 2015, the mean Ba of nine North Cascade glaciers was −3.10 m w.e., the most negative result in the 32-year record. The correlation coefficient of Ba was above 0.80 between all North Cascade glaciers, indicating that the response was regional and not controlled by local factors. The probability of achieving the observed 2015 Ba of −3.10 is 0.34%.


Introduction
The annual glacier mass balance measurements are the most accurate indicators of short-term glacier response to climate change [1,2]. The World Glacier Monitoring Service (WGMS) [2,3] has recognized that the changes in glacier mass are a key aspect of glacier monitoring, providing important information for assessing climatic changes, water resources, and sea level rise, and they maintain the most extensive data set on global glacier mass balance. The WGMS has recognized the continuous long term programs with consistent measurement programs as reference glaciers [3], there were 42 present world-wide in 2017. This network has proven valuable, but in many areas the number of glaciers that are monitored are limited [3]. There are three reference glaciers in the conterminous United States, and all are found in the North Cascades. The South Cascade Glacier (Figure 1) is monitored by the USGS (United States Geological Survey), and the Columbia Glacier and Rainbow Glacier are monitored by the North Cascade Glacier Climate Project (NCGCP) [3].
The North Cascade region in the United States extends from the Snoqualmie Pass to the Canadian Border and contains more than 700 glaciers. Glaciers in the North Cascade Range are important to water resources in many of the watersheds in the range [4,5]. The 2015 hydrologic year in the Pacific Northwest of North America was exceptional for its warmth and negative glacier mass balance.
Here we utilize the breadth of the mass balance record to examine the impact of climate on the glacier mass balance, and to place the exceptional 2015 mass balance year in context. NCGCP was founded in 1983 in order to monitor 10 glaciers throughout the range and identify the response of the North Cascade Range (Washington) glaciers to regional climate change [6,7]. The annual observations include the mass balance, terminus behavior, glacier surface area, and accumulation The NCGCP has continued Ba measurements on seven original glaciers that still exist, namely: Columbia, Daniels, Ice Worm, Lower Curtis, Lynch, Rainbow, and Yawning, as well as on Easton Glacier and Sholes Glacier, which were added in 1990. The annual measurement have been discontinued on three glaciers, namely: Lewis Glacier, Spider Glacier [7], and Foss Glacier, after it separated into several distinct sections. The Foss Glacier is becoming increasingly difficult to access and time spent doing so is not of a sufficient value, given the rapid recent demise and imminent loss of the glacier. The mass balance records provide both a direct measure of glacier runoff to rivers in the region and an understanding of that relationship with climate change [5,7,8]. No single glacier is representative of all of the others. In order to understand the causes and nature of changes in the glacier surface mass balance throughout a mountain range, it is necessary to monitor a significant number of glaciers [9]. The glacier mass balance in the North Cascades varies because of the geographic characteristics, including the aspect, elevation range, accumulation sources, and distance from crest of the range. When combined with the South Cascade Glacier [9], the NCGCP network represents the highest density of long term direct measurement, more than 25 years, of glacier mass balance over a given mountain range. The National Park Service is monitoring the mass balance of four glaciers in the North Cascades-namely, Silver Creek, North Klawatti, Noisy Creek, and Sandalee-starting the program in 1993, which are not used in this study because the duration of study has been less than 25 years [8].

Regional Climate
The North Cascades have had a temperate maritime climate with mild year-round temperatures, abundant winter precipitation, and dry summers. The two key climate variables for the glaciers have The NCGCP has continued Ba measurements on seven original glaciers that still exist, namely: Columbia, Daniels, Ice Worm, Lower Curtis, Lynch, Rainbow, and Yawning, as well as on Easton Glacier and Sholes Glacier, which were added in 1990. The annual measurement have been discontinued on three glaciers, namely: Lewis Glacier, Spider Glacier [7], and Foss Glacier, after it separated into several distinct sections. The Foss Glacier is becoming increasingly difficult to access and time spent doing so is not of a sufficient value, given the rapid recent demise and imminent loss of the glacier. The mass balance records provide both a direct measure of glacier runoff to rivers in the region and an understanding of that relationship with climate change [5,7,8]. No single glacier is representative of all of the others. In order to understand the causes and nature of changes in the glacier surface mass balance throughout a mountain range, it is necessary to monitor a significant number of glaciers [9]. The glacier mass balance in the North Cascades varies because of the geographic characteristics, including the aspect, elevation range, accumulation sources, and distance from crest of the range. When combined with the South Cascade Glacier [9], the NCGCP network represents the highest density of long term direct measurement, more than 25 years, of glacier mass balance over a given mountain range. The National Park Service is monitoring the mass balance of four glaciers in the North Cascades-namely, Silver Creek, North Klawatti, Noisy Creek, and Sandalee-starting the program in 1993, which are not used in this study because the duration of study has been less than 25 years [8].

Regional Climate
The North Cascades have had a temperate maritime climate with mild year-round temperatures, abundant winter precipitation, and dry summers. The two key climate variables for the glaciers have been the accumulation-season precipitation (November-April) and ablation-season temperature (May-September) [7,9]. A warming of 0.8 • C in the mean annual temperature from 1900 to 2012 was observed from in the North Cascades [10]. The warming was accelerated to +0.20 C per decade, for the 1980-2012 period. Every season, except spring, had experienced warming, particularly during the 1980-2012 period [10].
Approximately 70% of the region's precipitation occurred during the wet season (October-April) when the North Cascades were on the receiving end of the Pacific storm track [11,12]. From late spring to early fall, high pressure to the west kept the Pacific Northwest comparatively dry. Occasionally in the winter, warm fronts elevated temperatures and freezing levels, which resulted in rainfall at the glacier elevations. Rain on snow events had increased in frequency, this led to an increase in the ratio of the winter precipitation falling as rain versus it falling as snow. Mote [11] noted a decline in the snowpack storage efficiency in the Pacific Northwest, with the ratio between the total accumulation season precipitation and April 1 having retained snowpack SWE (snow water equivalent). Pelto [8] used the Diablo Dam and concrete weather stations for the total accumulation season precipitation for November-March-the period of accumulating the snowpack at the six SNOTEL sites, with long term records in the North Cascades-in order to determine the snowpack storage efficiency (Fish Lake, Lyman Lake, Park Creek, Rainy Pass, Stampede Pass, and Stevens Pass). For these six USDA (United States Department of Agriculture) SNOTEL stations that were utilized in this study, the mean April 1 SWE declined by 29% from 1946 to 2014. During the same period, the winter precipitation had increased slightly. The change in SWE was comparable with the change in the snowpack storage efficiency, which indicated that this was the primary cause of the reduced April 1 SWE at the North Cascade SNOTEL stations. The freezing level was a key factor in determining the snowpack storage efficiency, and an online application the North American Freezing Level Tracker was developed by Abatzoglou [12] that was utilized for the comparison of the freezing level during the winter season for Mount Baker.
The other key factor in glacier mass balance is the magnitude of ablation, which is primarily controlled by air temperature [13]. That air temperature is the key is indicated by the success of degree day factors (DDF) for assessing glacier ablation on glaciers in this region [13,14]. Nearly all of the ablation occurred during the May-September period, with the majority of the ablation having occurred from June-September. The most reliable weather station in the region was Diablo Dam, which was used for DDF derivation on the South Cascade Glacier [13]. An examination of trends in the melt season temperature at this station indicated a nearly identical pattern for May-September and June-September, which indicated that either period could be used to identify the melt season climate change. We utilized June-September, in this study, as the melt season. Six of the ten warmest melt seasons during the 1946-2014 period occurred since 2003. The long term melt season warming was 0.7 • C at Diablo Dam. The average June-September temperature from 2003 to 2014 was 0.6 • C above the mean of that for the 1946-2002 period.
The Pacific Decadal Oscillation Index (PDO) has been the leading principal component of North Pacific monthly sea surface temperature variability, poleward of 20 N [15]. During the positive PDO phase, warm weather was favored in the Pacific along the Northwest Coast and over the Pacific Northwest. During the negative phase, cool ocean water was found off the Northwest Coast and cooler temperatures were found across the Pacific Northwest [15]. In the past century, "cool" PDO regimes prevailed from 1890 to 1924 and again from 1947 to 1976, while "warm" PDO regimes dominated from 1925 to 1946 and again from 1977 to 1998 [15].
The El Niño/Southern Oscillation (ENSO) phenomenon was the most observable of the atmospheric circulation indices that led to year-to-year climate variability. ENSO positive events (El Nino) heralded abnormally warm sea surface temperatures (SST) over the eastern half of the equatorial Pacific. La Niña, was the opposite phenomenon, which was indicative of abnormally cold SST in the eastern half of the equatorial Pacific [16]. ENSO was an east-west atmospheric pressure seesaw that directly affected the tropical weather around the globe and indirectly impacted a much larger area [16]. The ENSO multivariate index (MEI-ENSO) that was used was based on the principal observed climate variables over the tropical Pacific. The index was a weighted average of the main ENSO features that were contained in the following six variables, namely: the sea-level pressure, east-west and north-south components of the surface wind, SST, surface air temperature, and total amount of cloudiness [16]. Positive MEI-ENSO values were usually accompanied by a sustained warming of the central and eastern tropical Pacific Ocean. Negative values of the MEI-ENSO index were associated with stronger Pacific trade winds and warmer sea temperatures in the Western Pacific to the north of Australia [16].
Bitz and Battisti [17] noted the importance of PDO and MEI-ENSO to glaciers in the region and observed that PDO had a greater influence during the 1960-1995 period. Josberger et al. [8] indicated that the importance of PDO had declined recently. They used three time periods, 1966-2004, 1966 to 1988, and 1989-2004, and found a significant change in the relationship between the PDO and the winter balances of the Wolverine Glacier, however less so for the South Cascade Glacier. Pelto [18] utilized the indices during the accumulation season as a first order forecast for the annual glacier mass balance, and found that the impact was the strongest when PDO and ENSO were either both positive or both negative.

Surface Mass Balance
Annual surface mass balance (Ba) was the difference between the annual accumulation of snow/ice and the loss of snow/ice by ablation. It was typically measured on a water year basis, beginning approximately October 1 and ending September 30, in the Northern Hemisphere.
Since 1984, NCGCP monitored the Ba of 9-10 glaciers every year [18][19][20]. Seven glaciers had a 32-year record, namely, the Columbia, Daniels, Ice Worm, Lower Curtis, Lynch, Rainbow, and Yawning glaciers. The Foss Glacier had a 30 year record (1984-2013) and was discontinued because of the glacier separating into several individual bodies. Sholes and Easton Glacier had a 26-year record . The glaciers represented a range of geographic characteristics and spanned the North Cascade Range (Table 1 and Figure 1). The key geographic variables were the glacier orientation, elevation, accumulation sources and distance to the mountain range watershed, and climate divide. The Columbia Glacier and the Rainbow Glacier were part of the 42 reference glaciers of the WGMS data set. Table 1. The geographic characteristics of the nine glaciers where the annual balance has been monitored annually ( Figure 1). Accumulation sources are listed in order of significance: wind drifting = WD, avalanche accumulation = AV, direct snowfall = DS. NCGCP measured the conditions on a glacier near the time of minimal mass balance, at the end of the water year, using a fixed date method. NCGCP methods emphasized the surface mass balance measurements with a relatively high density of sites on each glacier (>100 sites km −2 ), consistent measurement methods that were applied on fixed dates and at fixed measurement locations with consistent supervision [18,19,21]. The use of a high measurement density and consistent methods generated errors, which resulted from an imperfectly representative measurement network that was largely consistent and correctable; the error range had been observed at ±0.10-0.15 ma −1 [21]. Fischer [22] examined the mass balance errors and observed that the error had declined with increased density of measurements from 0.33 ma −1 , with a lower density to 0.10 ma −1 on a glacier that had a high density.
Any additional ablation that occurred after the last visit to a glacier was measured during the subsequent hydrologic year. The methods were reviewed in detail by Pelto [18][19][20][21].

Accumulation Area Ratio
At regional scales and on specific glaciers there were two common proxies for assessing mass balance without detailed observations. They were the AAR and equilibrium line altitude (ELA), both of which could be derived from satellite imagery or photographs [23,24]. The AAR was the ratio of a glacier that was in the accumulation zone. The ELA was the elevation at which the ablation equaled accumulation, on temperate alpine glaciers this was coincidental with the transient snow line (TSL) at the end of the melt season. The ELA was not typically an easily discernible line or elevation on the North Cascades glaciers, because of the variability of the snow accumulation from the impacts of wind and avalanche redistribution. The ELA could be calculated from the balance gradient, but in such cases it was not an independent variable. The AAR was a more accurately determined parameter and a better proxy in this case [20]. The accumulation zone was a patchwork of the retained snowpack and ablation areas. Each patch of the retained snowpack was mapped and included in the AAR determination. The AAR in the North Cascades was determined from either photographs or direct surface mapping, by measuring the GPS along the TSL around each patch of the retained snowpack. The AAR-Ba method was proven to be reliable for the annual balance estimates [20,24,25]. A comparison of the annual AAR and Ba observations in WGMS [26,27] indicated correlation coefficients (Pearson's r in all cases in the paper) ranging from 0.70 to 0.92 for fifteen glaciers, with at least 10 years of records. The World Glacier Monitoring Service (WGMS) had adopted the reporting of AAR with mass balance values [26,27] and was plotting the relationship for each glacier. A combination of the AAR observations and Ba measurements from 1984 to 2015 on the North Cascade glaciers provided an opportunity to assess the relationship for each glacier and the variability between glaciers.
The AAR0 value was the AAR for a glacier with an equilibrium mass balance [28]. Braithwaite and Muller [29] noted that the AAR0 for an alpine glacier with an equilibrium balance had averaged 0.67. The mean AAR0 that was reported for 89 temperate alpine glaciers to the WGMS was 0.57, the AAR0 value was determined from a regression of the observed Ba and AAR [26,27].

Annual Mass Balance
The glaciers in the North Cascades exhibited consistent Ba responses to climate from year to year [19,20]. Figure 2 illustrates the closely correlated pattern of the Ba fluctuations from 1985 to 2015.
In most years, all of the glaciers responded in step with each other to the variations in the winter retained snowpacks and summer temperature. There was a mean annual range between the North Cascade glaciers of 0.8 m w.e., from the maximum Ba to the minimum Ba in a given year, but the inter-annual trend was the same for each glacier ( Table 2). This regional response was indicated by the high cross-correlation values of the Ba between individual glaciers, ranging from an r 2 of 0.80 to 0.98, between each glacier pair, including the South Cascade Glacier that was monitored by the USGS (Table 3) [2,30]. The data for 2013-2014 had not been finalized for the South Cascade Glacier, hence, the comparison period ended in 2012 for the cross correlation. The correlation between the Ba and accumulation season precipitation was also quite similar, ranging from a low of 0.67 on the South Cascade Glacier to a high of 0.79 on the Columbia Glacier ( Table 4). The correlation between the Ba and ablation season temperature was also quite similar, ranging from a low of 0.67 on the Columbia Glacier to a high of 0.80 on the Sholes Glacier (Table 4).
The mean Ba had been −0.54 m w.e.a −1 for the 1984-2015 period, on the glaciers that were monitored annually, and ranged from a maximum of −0.44 m w.e.a −1 to a minimum of −0.67 m w.e.a −1 . The mean Ba, −0.54 m w.e.a −1 for the North Cascade glaciers during the 1984-2015 period, was quite similar to the WGMS global reference glacier mean Ba of 0.56 ma −1 [31]. The mean cumulative mass balance loss had been −17.2 m w.e., which was 19-20 m of glacier thickness that was lost. The mean thickness of the several glaciers that were investigated by the geophysical methods ranged from 30 to 60 m. Each glacier that was observed was larger than the mean of the North Cascade glaciers or the mean of the glaciers that were observed by NCGCP [32,33]. Thus, at least 30% of the volume of these glaciers had been lost since 1984. The observations by the USGS at the South Cascade Glacier indicated that the mean Ba from 1984 to 2012 was −0.75 m w.e.a −1 on the South Cascade Glacier, which was 0.21 m w.e.a −1 more negative than the glaciers that we monitored. The South Cascade Glacier had a much less negative mean Ba of −0.15 m w.e.a −1 from 1956 to 1975 [30].
Water 2018, 10, x FOR PEER REVIEW 6 of 15 thickness of the several glaciers that were investigated by the geophysical methods ranged from 30 to 60 m. Each glacier that was observed was larger than the mean of the North Cascade glaciers or the mean of the glaciers that were observed by NCGCP [32,33]. Thus, at least 30% of the volume of these glaciers had been lost since 1984. The observations by the USGS at the South Cascade Glacier indicated that the mean Ba from 1984 to 2012 was −0.75 m w.e.a −1 on the South Cascade Glacier, which was 0.21 m w.e.a −1 more negative than the glaciers that we monitored. The South Cascade Glacier had a much less negative mean Ba of −0.15 m w.e.a −1 from 1956 to 1975 [30].   thickness of the several glaciers that were investigated by the geophysical methods ranged from 30 to 60 m. Each glacier that was observed was larger than the mean of the North Cascade glaciers or the mean of the glaciers that were observed by NCGCP [32,33]. Thus, at least 30% of the volume of these glaciers had been lost since 1984. The observations by the USGS at the South Cascade Glacier indicated that the mean Ba from 1984 to 2012 was −0.75 m w.e.a −1 on the South Cascade Glacier, which was 0.21 m w.e.a −1 more negative than the glaciers that we monitored. The South Cascade Glacier had a much less negative mean Ba of −0.15 m w.e.a −1 from 1956 to 1975 [30].

Annual Mass Balance 2015
On June 15, when the automatic weather station and discharge station were installed adjacent to the Sholes Glacier, the snowpack was similar to a typical early August snow cover. On the Sholes Glacier, the AAR fell from 0.55 on 9 July to 0.00 on 9 September (Figure 4). This was the first year since the monitoring had begun in 1984 that the mean AAR in early August was below 0.25. The result was an exposure of the older firn layers and a general decrease in albedo. In early August, the AAR was below 0.1 for all of the glaciers, except for the Easton Glacier. On the Columbia Glacier, the AAR on August 1 was the lowest observed yet at 0.12, with six weeks remaining in the melt season ( Figure 5). The early exposure of glacier ice was important as the melt rate was faster, as was indicated by the greater melt factor [13,14].

Annual Mass Balance 2015
On June 15, when the automatic weather station and discharge station were installed adjacent to the Sholes Glacier, the snowpack was similar to a typical early August snow cover. On the Sholes Glacier, the AAR fell from 0.55 on 9 July to 0.00 on 9 September (Figure 4). This was the first year since the monitoring had begun in 1984 that the mean AAR in early August was below 0.25. The result was an exposure of the older firn layers and a general decrease in albedo. In early August, the AAR was below 0.1 for all of the glaciers, except for the Easton Glacier. On the Columbia Glacier, the AAR on August 1 was the lowest observed yet at 0.12, with six weeks remaining in the melt season ( Figure 5). The early exposure of glacier ice was important as the melt rate was faster, as was indicated by the greater melt factor [13,14]. In 2015, out of the nine glaciers where the Ba was examined, the AAR was 0.00 on seven of the glaciers, 0.05 on the Rainbow Glacier, and 0.26 on the Easton Glacier.     [30]. In 1958, a more negative mass balance was reported, although the winter and summer balances were not measured. This was the first full year of observations and the methods had not been standardized, hence, it was not considered part of the official USGS record [30]. In 2005, the South Cascade Glacier had a mass balance of −2.45 m w.e., the second most negative Ba since 1984.
The fact that the 10 glaciers in the North Cascades, with a Ba record of more than 25 years, had experienced their most negative mass balance in one year, indicated how unusual 2015 was. A probability function was used to determine the likelihood of the mean balance achieving a Ba of −3.19 m w.e. or less, given that the mean for the 1983-2014 period was −0.47 m w.e. and had a standard deviation of 0.97. The probability of achieving the observed 2015 Ba of −3.10 was 0.34%. The same approach was applied to individual glaciers that yielded a range of 0.1-2.5%.  [30]. In 1958, a more negative mass balance was reported, although the winter and summer balances were not measured. This was the first full year of observations and the methods had not been standardized, hence, it was not considered part of the official USGS record [30]. In 2005, the South Cascade Glacier had a mass balance of −2.45 m w.e., the second most negative Ba since 1984.
The fact that the 10 glaciers in the North Cascades, with a Ba record of more than 25 years, had experienced their most negative mass balance in one year, indicated how unusual 2015 was. A probability function was used to determine the likelihood of the mean balance achieving a Ba of −3.19 m w.e. or less, given that the mean for the 1983-2014 period was −0.47 m w.e. and had a standard deviation of 0.97. The probability of achieving the observed 2015 Ba of −3.10 was 0.34%. The same approach was applied to individual glaciers that yielded a range of 0.1-2.5%.

Climate in 2015
The 2015 winter accumulation season featured 51% of the mean  winter snow accumulation at six long-term USDA SNOTEL stations in the North Cascades, namely, Fish Creek, Lyman Lake, Park Creek, Rainy Pass, Stevens Pass, and Stampede Pass. This was exceptional as it was the second lowest out of the 32 years of the mass balance observation series. The winter season was exceptional for warmth, being the warmest winter season on record in the state of Washington [34]. The freezing level in 2015 averaged 1645 m in the Mount Baker region from November-March, compared with an average of 1077 m [12]. The previous record for the mean November-March freezing level, since the record began in 1948, was 1500 m. The result was the snowpack at the six SNOTEL sites averaging 0.52 m w.e., which was 51% of the 1984-2014 average. The minimum snowpack was in 2005, with 44% of the mean snowpack, 2005 had less precipitation. The limited snowpack was not due to reduced precipitation, as the November-April precipitation at the Diablo Dam and concrete weather stations was 5-10% above the 1984-2014 mean. In 2015, the snowpack storage efficiency was the lowest, with the mean 1 April retained SWE at the six long-term SNOTEL stations being 19% of the November-March precipitation. Compared with a 1984-2014 mean of 41%, the second lowest year was 2005 at 20%. The snowline on 1 May 2015, which was the end of the typical accumulation season in the Mount Baker region, was 1400 m versus the long term average of 800 m.
In 2015, the mean May-September temperature at Diablo Dam was 2.2 • C warmer than the long term mean, and it was the second warmest to 1958 in the 1950-2015 record. For June-September, the mean temperature was 2.0 • C warmer than the long term mean, and was also second to 1958 as the warmest. The combination of the warmest melt season in over 50 years and the second lowest accumulation season snowpack in the last 30 years was a good indication that the glacier mass balance would be quite negative.
In 2015, the sea surface temperature waters that had developed in the winter of 2013/14, persisted off the coast of the Pacific Northwest, with anomalies generally exceeding 2 • C [34,35]. Such conditions were associated with a positive PDO and warm air temperatures in the Pacific Northwest [34]. From the winter of 2013/14 to 2014/15, the Northeast Pacific experienced the largest marine heatwave ever recorded [35]. In 2015, the warmest sea surface temperature paralleled the entire Pacific North American coast, which exhibited a PDO pattern [35]. A comparison of the 2015 values for the aforementioned climate parameters and the longer term means is provided in Table 5. A comparison of annual, winter, and summer values for PDO from 1984 to 2015 indicated that 2015 ranked in the top five for each seasonal measure (Table 6). Examining the 1984-2015 period for ENSO indicated that, in 2015, the summer and annual values ranked in the top five. The 2015 hydrologic year had the most positive PDO of any year, however ENSO was not the most positive. The combined value of PDO-ENSO in 2015 ranked in the top five for both the annual and seasonal rankings, during the 1984-2015 period. The positive values for both ENSO and PDO were inversely related to the glacier mass balance and also reinforced each other [18]. This illustrated that 2015 was predisposed for a negative annual balance, based solely on the climate indices. The link between the PDO and ENSO forced the anomalous atmospheric and sea surface temperature event in 2015 [35].

Discussion
The examination of the Pacific Northwest glacier mass balance for the full period of research for the NCGCP and USGS benchmarked glaciers indicated that the winter balance fluctuations of the maritime glaciers dominated the Ba fluctuations [8]. Looking at the full period of record from 1984 to 2015, it was apparent that the summer balances of the benchmark glaciers had become more negative during recent years, as a result of the warmer and drier summers and, hence, the summer balance had a greater role in determining the Ba [8]. Pelto [18]  Here we used two periods, 1984-1998 and 1999-2014, with 1999 marking a transition from a period of primarily positive ENSO and PDO, to a period of negative ENSO and PDO, in order to compare the annual mass balance for each glacier to the winter PDO index (October-April) and the annual mean MEI-ENSO (October-September). During the first period, the mean correlation was −0.34 for the winter PDO and −0.59 for the annual MEI-ENSO. During the latter period, the mean correlation had increased to −0.54 for the winter PDO and −0.66 for the annual MEI-ENSO. The annual MEI-ENSO index remained a better indicator of the annual mass balance in the North Cascade Range than the winter PDO.
Overall, 2015 was exceptional in terms of the sea surface temperature in the Northeast Pacific, air temperature in the North Cascade region, 1 April snowpack, winter freezing levels, and in the positive values of PDO. ENSO was significantly positive but not exceptional. Given the similar correlation between the key climate variables of the ablation season air temperature and accumulation season precipitation and Ba, it was expected that the 2015 climate conditions would lead to record Ba losses in 2015. The mean Ba for the WGMS reference glacier network in 2015 was −1.12 m w.e., the second most negative after 2003 [31]. The mean Ba of all of the 149 WGMS reporting glaciers in 2015 was −1.04 m. Of these, 11 glaciers reported a loss of more than 3 m w.e., nine were in the state of Washington, in addition to the NCGCP glaciers-namely, the Eel Glacier in Olympic National Park and the Noisy Creek Glacier in the North Cascades National Park, as reported by the National Park Service [31]. Of the WGMS reference glacier network, the Columbia Glacier had the most negative balance of any glacier. The Ba continued to trend negatively, both globally and in the North Cascades, which indicated that, instead of approaching equilibrium as the glaciers retreated, they continued to be in disequilibrium with the current climate (Figures 4 and 5) [2,36].

Conclusions
The impact of the warmer temperatures and reduced winter snowpack on the North Cascade glaciers has been an average annual balance of −0.54 ma −1 over the past 32 years. The net loss of −17.3 m w.e. represents a significant portion (30% of the total glacier volume), resulting in a substantial retreat and thinning. The resultant retreat is ubiquitous, rapid, and increasing. There is no evidence that the North Cascade glaciers are close to equilibrium. Their ongoing thinning indicates that all of the glaciers will continue to retreat for the foreseeable future. In cases where the thinning is substantial along the entire length of the glacier, no point of equilibrium can be achieved with present climate, and the glacier is unlikely to survive [36].
Continued glacier retreat is inevitable; 75% of the North Cascade glaciers that we observed are in disequilibrium and will melt away during this century with the current climate [36]. The loss of glacier area will lead to further declines in summer runoff in glacier fed rivers, as the glacier area available for melting in the summer declines [37,38]. This will impact salmon in streams such as the Nooksack River [39]. The correlation between the cumulative mass balance records of glaciers in various alpine glaciated regions around the globe, suggest that the global climate signal is the dominant driver of mass balance losses [2].