Next Article in Journal
Production of Microalgal Biomass Using Aquaculture Wastewater as Growth Medium
Next Article in Special Issue
A Decrease in the Regulatory Effect of Snow-Related Phenomena in Spanish Mountain Areas Due to Climate Change
Previous Article in Journal
Encounter Probability and Risk of Flood and Drought under Future Climate Change in the Two Tributaries of the Rao River Basin, China
Previous Article in Special Issue
Hydro-Meteorological Characterization of Major Floods in Spanish Mountain Rivers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of North Atlantic Oscillation on the Snowpack in Iberian Peninsula Mountains

Instituto Pirenaico de Ecología, CSIC. Campus de Aula Dei, Av. Montañana 1005, 50059 Zaragoza, Spain
*
Author to whom correspondence should be addressed.
Water 2020, 12(1), 105; https://doi.org/10.3390/w12010105
Submission received: 11 November 2019 / Revised: 19 December 2019 / Accepted: 20 December 2019 / Published: 28 December 2019

Abstract

:
The North Atlantic Oscillation (NAO) is considered to be the main atmospheric factor explaining the winter climate and snow evolution over much of the Northern Hemisphere. However, the absence of long-term snow data in mountain regions has prevented full assessment of the impact of the NAO at the regional scales, where data are limited. In this study, we assessed the relationship between the NAO of the winter months (DJFM-NAO) and the snowpack of the Iberian Peninsula. We simulated temperature, precipitation, and snow data for the period 1979–2014 by dynamic downscaling of ERA-Interim reanalysis data, and correlated this with the DJFM-NAO for the five main mountain ranges of the Iberian Peninsula (Cantabrian Range, Central Range, Iberian Range, the Pyrenees, and the Sierra Nevada). The results confirmed that negative DJFM-NAO values generally occur during wet and mild conditions over most of the Iberian Peninsula. Due to the direction of the wet air masses, the NAO has a large influence on snow duration and the annual peak snow water equivalent (peak SWE) in most of the mountain ranges in the study, mostly on the slopes south of the main axis of the ranges. In contrast, the impact of NAO variability is limited on north-facing slopes. Negative (positive) DJFM-NAO values were associated with longer (shorter) duration and higher (lower) peak SWEs in all mountains analyzed in the study. We found marked variability in correlations of the DJFM-NAO with snow indices within each mountain range, even when only the south-facing slopes were considered. The correlations were stronger for higher elevations in the mountain ranges, but geographical longitude also explained the intra-range variability in the majority of the studied mountains.

1. Introduction

The North Atlantic Oscillation (NAO) is a mode of atmospheric circulation, and its index value is calculated as the difference in atmospheric pressure for the dipole centered on Iceland and the Azores Islands [1]. The NAO is considered to be the main pattern of climate variability in the Northern Hemisphere [2,3], especially during winter. Negative (positive) NAO phases generally are related to weak (strong) westerly winds to Europe, leading to positive (negative) anomalies of temperature and precipitation in Southern (Northern) Europe [4,5,6]. It is also clear that NAO anomalies are more related with the winter precipitation variability than with the interannual winter temperature fluctuations [7,8].
As NAO fluctuations are related with the winter temperature and precipitation over much of the Northern Hemisphere, it is unsurprising that many studies have reported that the NAO is a good indicator of the duration and depth of the snowpack. This is the case for Mediterranean mountains [7], the Pyrenees [9,10], the Alps [11,12,13], Nordland [14], Western Poland [15], Bulgaria [16], and the Russian Arctic [17]. At wider spatial scales, Henderson and Leathers [18] associated negative (positive) phases of the NAO with expansion (reduction) of the snow cover in Europe. Keylock [19], Jomelli et al. [20], and García et al. [21] also found significant relationships between NAO phases and the frequency of avalanches in Iceland, the French Alps, and the Pyrenees, respectively. The main value of understanding the relationship between snow and interannual fluctuations of the snowpack is that, despite the stochastic nature of the NAO, the statistic and dynamic seasonal forecasting of this phenomenon has been very promising, to the point of achieving useful levels of predictability [22,23,24]. Improving predictions of snowpack anomalies has a clear applied importance because of the direct impact of the snowpack on water availability, plant and animal phenology, and a wide range of economic activities in mountain and downstream areas [25,26]. In addition, a good understanding of the relationship between the NAO and the snowpack is necessary to accurately interpret the reported trends in the snowpack over recent decades, as a strong correlation has been observed between trends in the snowpack and the duration of dominant NAO anomalies [27].
Although studies have provided clear insights into how the NAO is related with the snowpack in many mountains and cold regions of Europe, at the local scale, many uncertainties about its effects remain. This is because mountain ranges present a topographical barrier to main wind flows, causing shadow effects in the main atmospheric circulation patterns [28,29]. In addition, elevation also affects the impact relationship of NAO with snow in mountain areas because of its greater relationship with precipitation relative to temperature. Thus, at high elevations, interannual variability of the snowpack is mainly driven by precipitation, so such areas are expected to exhibit a stronger response to the NAO than lower elevation sites, where temperature has a major role in snow variability [30]. López-Moreno and Vicente-Serrano [9] suggested that this phenomenon is the main explanation for the control of NAO on the snowpack in areas above 2000 m a.s.l. in the Spanish Pyrenees. The spatial variability in the relationship between the NAO and snow cover makes it difficult to assess NAO effects in many areas, because of the limited data on snowpacks, especially in mountain areas at the highest elevations.
The aim of this study is to extend knowledge of the NAO and snowpack relationship in the main mountainous regions of the Iberian Peninsula. This is one of the highest elevation regions of Europe, where the accumulation of snow is of major importance to water resource availability [31,32], and on which many economically important activities, including agriculture, forestry, and winter tourism, depend [33,34]. Because of the limited data on snowpacks in the Iberian Peninsula, we used a recently developed database on snow water equivalent (SWE) and snow depth (SD) for the main mountain areas of the Iberian Peninsula [35], to enable us to undertake the first detailed assessment of the relationship between the winter (December, January, February, and March) NAO (DJFM-NAO) on the interannual variability in the peak SWE and the duration of snow on the ground.

2. Study Area

The Iberian Peninsula is located at the southwest edge of Europe, between latitudes 36° N and 44° N. It is an area of contrasting climate, as a consequence of Atlantic and Mediterranean influences. The topography is characterized by two main plateaus surrounded by mountainous areas and comprises five main mountain ranges that have seasonal snow cover, including the Cantabrian Range, the Central Range, the Iberian Range, the Pyrenees, and the Sierra Nevada. The highest elevation (3478 m a.s.l.) is Mulhacén Peak in the Sierra Nevada, followed by Aneto Peak in the Pyrenees (3404 m a.s.l.). The other mountainranges exceed 2000 m a.s.l., reaching 2648 m a.s.l. (Torrecerredo Peak) in the Cantabrian Range, 2596 m a.s.l. (Almanzor Peak) in the Central Range, and 2314 m a.s.l. in the Iberian Range.
The spatial distribution of these mountain ranges (Figure 1), which extend within Spain from west to east at differing latitudes and distances to the sea, results in wide snowpack variability at similar elevations, even within the same mountain range [36].

3. Methods

3.1. NAO Index and Temperature and Precipitation Patterns

The NAO index, which is defined as the normalized sea surface pressure (SLP) difference between the Icelandic low pressure and the Azores high pressure regions [1], is a widely used indicator of the atmospheric situation over the North Atlantic. In this study the NAO index was calculated from the SLP difference between the positions 65° N 20° W (near Reykjavik, Iceland) and 35° N 5° W (near Gibraltar Straight), following the procedure and recommendations of Osborn et al. [37], Trigo et al. [38], and Fernández-González et al. [39], among others. The SLP data were downloaded from the NCAR SLP reanalysis dataset (https://climatedataguide.ucar.edu/climate-data/ncar-sea-level-pressure; last accessed 20 October 2019). The normalization was performed in relation to normal values for the study period 1979–2015.
The accumulated winter precipitation and the minimum and maximum averaged daily temperatures were linearly correlated with the DJFM-NAO index. The temperature and precipitation meteorological variables were obtained from the regional atmospheric simulation described below. In this study the winter period was considered to encompass the months December, January, February, and March.

3.2. Snowpack Database and Statistical Analysis

The SWE and SD database used in this study was developed using the weather research and forecast (WRF) model [40], forced using the ERA-Interim reanalysis [41,42]. The spatial resolution of the atmospheric simulations was ~10 km, and the simulations covered the period 1979–2014 [43]. The surface meteorological outputs of the simulations were subsequently projected to various elevation bands that differed by 100 m, using an array of psychrometric and radiative formulae and lapse rates [35]. The new meteorological data were used as forcing data for the physically based energy and mass balance snowpack factorial snow model (FSM) [44]. The final product was a semi-distributed daily SWE and SD database. Alonso-González et al. [35] validated the database, using in situ observations and remote sensing information, and were able to reproduce the spatial and temporal patterns of the snowpack over the entire Iberian Peninsula.
We computed the annual peak SWE and the annual snowpack duration for each elevation band from the snowpack database. These values were linearly correlated with the averaged DJFM-NAO, along with the Pearson’s R value. To identify different responses of the snowpack on the windward and leeward sides of each mountain range, we separated the correlation values for each mountain range along the principal longitudinal axis, which runs generally from west to east. The Wilcoxon–Mann–Whitney test [45] was used to identify statistical differences between the groupings of north- and south-facing slopes.
Following calculation of the Pearson’s R values between the DJFM-NAO and the snow indices, we correlated the value with the elevation and geographical longitude, to investigate the spatial patterns of DJFM-NAO influence on the snowpack within each mountain range.
Finally, to highlight the spatial patterns of the correlation between the DJFM-NAO and the snow indices, 500 random samples of 100 values of the Pearson’s R values were calculated. We performed partial correlations between snow indices and longitude (removing the effect of the elevation), and elevation (removing the effect of the longitude). This recursive resampling of the dataset was performed to avoid using excessively big samples, which can lead to overestimation of statistically significant correlations and has associated problems caused by spatial autocorrelation [46,47].

4. Results and Discussion

4.1. Relationship between the DJFM-NAO and Temperature and Precipitation

Except for small areas on the Mediterranean coast, most of the Iberian Peninsula showed positive correlations between the NAO and the maximum temperature (Figure 2), being consistent with previous studies in the Mediterranean mountain ranges [48]. A clear difference was evident between the north- and south-facing slopes along the main longitudinal axis of the mountain ranges, including in the Central Range (north approximately 0.5, south approximately 0.8) and Cantabrian Range (north approximately 0.4, south approximately 0.8). This result is consistent with previous studies based on meteorological observations [49], and can be explained by the Foehn effect on the northern slopes during negative NAO phases [50]. This phenomenon is very common in the Cantabrian and Pyrenean ranges during DJFM-NAO negative phases.

4.2. Relationship between the Winter NAO and Temperature and Precipitation

High elevation areas generally showed strong positive correlations, which we attribute to the fact that mountain climates are strongly related to the movement of air masses [51]. Thus, during positive NAO phases, the maximum air temperature on the slopes and summits increase, and during negative NAO phases, the lapse rates increase and the summits cool, although the existence of southern fluxes [28].
The correlations between the DJFM-NAO and minimum temperatures tended to be positive, but showed a progressive decrease towards the southwest and lower areas of the main hydrological basins. In the valley bottoms the correlations became negative, as reported by del Río et al. [52] for the Duero basin. This progressive change is because of the increased cloudiness during negative NAO phases in this Atlantic southwest area, related to the associated southern and western cyclonic flows. The cloudiness slows the nighttime reduction in the air temperature [53], but this effect does not penetrate fully inland to the north, making the signal less obvious.
The DJFM-NAO and temperature were generally positively correlated. During negative DJFM-NAO phases the occurrence of more frequent cyclonic circulations brings lower temperatures than during anticyclonic periods, but the most frequent wind directions during negative DJFM-NAO phases has a clear Atlantic influence, which generally brings mild temperatures during snowfall events. This causes a marked increase of the elevation where precipitation changes from liquid to solid phase compared to other synoptic situations over Iberia [28,54].
The correlation between monthly accumulated precipitation and the NAO index was negative for the entire Iberian Peninsula, as previously reported by other authors over central Portugal [38] or in the northwestern part of Spain [39], as well as in the main mountain ranges of the Iberian Peninsula [48] (Figure 2); although, for small Mediterranean and Cantabrian areas, the correlations were negligible [55,56]. Southwest areas showed the greatest negative correlations (values of approximately 0.8), which is consistent with the findings of Muñoz-Díaz and Rodrigo [57], where they found different areas of correlation of the precipitation with the NAO index over the Iberian Peninsula and, more specifically, Queralt et al. [58] about extreme precipitation events associated with DJFM-NAO phases.
There was differentiation between the highly correlated southern slopes and the less correlated northern slopes, as previously reported for the Spanish Pyrenees by Vada et al. [59] and Buisan et al. [60], and for the Alps by Stefanicki et al. [51]. This differentiation was especially evident for the mountain ranges located at higher elevations (i.e., the Cantabrian Range and the Pyrenees).
The differences between the southward and northward slopes of the mountain ranges were clear for the Pyrenees (north approximately 0.1; south approximately 0.7), the Cantabrian Range (north approximately 0.3; south approximately 0.6), and the Iberian Range (north approximately 0.2; south approximately 0.6). The differences were less marked for the Central Range (north approximately 0.4; south approximately 0.5), and for the Sierra Nevada were almost negligible. This effect is caused likely by the geographical position of Sierra Nevada and Central Range, and its limited potential to be an effective topographical barrier against the westerly and southwesterly wind directions associated to negative DJFM-NAO phases.

4.3. NAO Spatial Influence on Peak SWE and Snow Season Duration

Figure 3 shows that the DJFM-NAO has a clear correlation with the interannual variability of the snow duration on the southern slopes of the Pyrenees and the Cantabrian, Central, and Iberian ranges. The correlation of the DJFM-NAO with snow duration was high and statistically significant at the highest elevations of the Sierra Nevada, but with no statistically significant differences between the north- and south-facing slopes. However, the other mountain areas showed statistical significance (p-value < 0.05; Wilcoxon–Mann–Whitney test) between north- and south-facing slopes. This finding is consistent with the spatial patterns shown in Section 3.2, between the NAO and the interannual variability of winter temperature and precipitation.
The spatial patterns and magnitude of the correlations for peak SWE were very similar to those for snow duration (Supplementary Figure S1). However, there were differences between these snow variables. Figure 4 shows the absolute difference in the Pearson’s correlation coefficients of the DJFM-NAO with the snow duration and the peak SWE for each pixel for the south-facing slopes, where high correlation values were found (except for the Sierra Nevada, for which all pixels were considered). The greatest differences were found for the higher elevations in the Pyrenees and the Sierra Nevada, where the correlations for peak SWE were higher than those for snow duration. These two mountain ranges exhibit large areas above 2000 m a.s.l, in contrast with the other mountain ranges. We hypothesize that, at high elevations, the correlation of the DJFM-NAO with the temperature and precipitation leads to high levels of snow accumulation during the DJFM-NAO phases, as a consequence of the driving mechanism that controls the DJFM-NAO; however, this is not reflected in the duration of the snow cover. This is mainly because snow duration also depends on solid precipitation occurring during April and May, which is not related to the DJFM-NAO. In addition, the very late and rapid melting that occurs in these high elevation areas [61] is mostly due to the marked increase of the incoming solar radiation at the end of the snow season period.
Figure 4 shows that there was large spatial variability in the correlation coefficients between the DJFM-NAO and the snow indices within each of the analyzed mountain ranges. The influence of elevation on the distribution of the DJFM-NAO correlation values between snow duration and the peak SWE (on south-facing slopes, except in the Sierra Nevada, where both slopes were considered) is shown in Figure 5, and the comparable influence of geographical longitude is shown in Figure 6.
For the Cantabrian Range, most of the correlations of the DJFM-NAO with snow duration ranged from 0.25 to 0.75, and with peak SWE ranged from 0 to 0.75. This variability showed no relationship to elevation, but a clear link with longitude was evident, with an increase in correlation westward. For the Pyrenees, most correlations ranged from 0 to 0.75 for snow duration and peak SWE. While the response of snow duration to elevation was unclear, there was a strong negative relationship for peak SWE. The increase in the correlation of the DJFM-NAO on snow magnitude at high elevations has previously been reported by López-Moreno and Vicente-Serrano [9], where they found statistically significant correlations of the DJFM-NAO with snowpack observations above 1700 m a.s.l., being consistent with our results. These findings suggest that the frequent westerly and southwesterly air fluxes during negative DJFM-NAO years bring mild temperatures during snowfall events, which may cause rain precipitation at low and mid elevations, but heavy snowfalls over 1800–2000 m a.s.l. The absence of a vertical gradient on the relationship of the DJFM-NAO on snow duration is explained by the presence of snow cover and snowfalls until late spring (end of May to early June); these are also highly dependent on spring precipitation, which is unrelated to the DJFM-NAO. Figure 6 does not show statistically significant correlations with longitude but reveals a clear pattern for both indices, characterized by the highest correlations (values reaching 0.75) for the Central Pyrenees, and a decrease westward and eastward. This effect is likely a consequence of the small surface with elevation above 1800 m a.s.l. in the western part compared to the Central Pyrenees, causing low correlation between DJFM-NAO and snow indices (Figure 5). The also-low correlations at the easternmost part are explained because southwestern air fluxes, which are frequent during negative DJFM-NAO phases and generally affect only the Western and Central Pyrenees. Thus, DJFM-NAO are not clearly related with both temperature and precipitation in this sector (Figure 2). In the Central Range, the correlation with the DJFM-NAO ranged from 0 to 0.75 for snow duration and 0.1 to 0.6 for peak SWE. This variability showed a clear negative trend with elevation for both indices. There was no strong relationship with longitude for the correlations between the snow indices and the DJFM-NAO for the Central Range. For the Iberian Range, most of the correlations with the DJFM-NAO ranged from 0.1 to 0.6 for snow duration, and 0 to 0.75 for peak SWE. The reduced area at high elevations and their limited geographical extent made it difficult to distinguish clear spatial patterns in the relationship of the DJFM-NAO and the snow indices. The Sierra Nevada showed the largest variability in the correlation values for both indices, ranging from slightly positive values to 0.75. For both indices the impact of the DJFM-NAO was clearer at higher elevations and eastward in this southern mountain range.
This study shows that the relationship of the DJFM-NAO and the snowpack can be seen in long-term trends of snow cover over wide areas of the Iberian Peninsula. Thus, very different trends in snow accumulation have been found in the Pyrenees when different periods have been analyzed. López-Moreno [54] found a statistically significant decrease in snow accumulation in the central Spanish Pyrenees for the period 1950–1999, but Buisan et al. [60] found no significant trends for the period 1985–2014. Such disagreements were explained by strong decadal fluctuations in the snow data, which correlated with the DJFM-NAO. The marked spatial differences observed in the correlation between the DJFM-NAO and snow indices may explain the marked contrast in snow trends over very short distances, as has been observed in the Pyrenees [10,62]. The presence of areas showing strong correlations between the snowpack and the DJFM-NAO, independent of the study period (as demonstrated by Buisan et al. 2015 for the Pyrenees), provides the promise of improving the seasonal forecasting of snowpack based on short- and medium-term projections for the DJFM-NAO from climate models [63]. This could improve optimization of water management and benefit the skiing industry, both of which are highly affected by interannual fluctuations in snow accumulation and the duration of the snowpack.

5. Conclusions

The DJFM-NAO is related to the atmospheric circulation patterns over the Iberian Peninsula during the snow-accumulation period. The NAO is related to contrasts in precipitation and temperature over the Iberian Peninsula, and hence related to the snow-accumulation and melting processes.
The mountain ranges analyzed in this study are topographic barriers to the main westerly and southwesterly air fluxes related to the DJFM-NAO and explain their variable impacts on the snowpack. Thus, the snowpack on north-facing slopes along the main mountain range axes did not show statistically significant correlations in most of its surface, but there were significant correlations in south-facing slopes. The Sierra Nevada was the only mountain range where this north south difference was not observed. For all Iberian mountain areas, the statistically significant correlations between the DJFM-NAO and snow indices were negative, and, in some cases, R values exceeded 0.7. Marked spatial differences were found, even for those slopes where the DJFM-NAO generally affects the temporal evolution of the snow. In most of the mountain areas, the strength of the relationship between the DJFM-NAO and peak SWE increased with elevation. This also occurred for snow duration, except in the Pyrenees, where spring snowfall can have a strong influence on the snow duration at high elevations, and snowmelt occurs later and more rapidly. This explains why the snowpack duration and peak SWE are not always related, and explains the weaker correlations between the DJFM-NAO and snow duration compared with peak SWE.
In several mountain areas, we found spatial differences that were related to geographical longitude. This was evident for the Cantabrian Range and the Sierra Nevada, where the correlation values decreased from west to east. For the Pyrenees, the maximum correlations were found for the central valleys and decreased toward its western and eastern edges. The results of this study have clear implications for understanding the effect of the length of the study period on the detection of robust long-term trends in snow data, and also spatial differences in trend analyses over short distances.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4441/12/1/105/s1, Figure S1: Pearson’s R values between the DJFM-NAO and the peak SWE. The lines represent the boundaries of the main catchments in the Iberian Peninsula.

Author Contributions

E.A.-G. designs and conducted the analysis and wrote the paper. J.I.L.-M. and J.R. supervise the interpretations of the results. F.M.N.-S. contributes in the interpretation of the relashionship between the NAO index and the temperature and precipitation. All authors have read and agreed to the published version of the manuscript.

Funding

Esteban Alonso-González is the recipient of a pre-doctoral FPI grant by the Spanish Ministry of Economy and Competitiveness (BES-2015-071466). This study was funded by the Spanish Ministry of Economy and Competitiveness project CGL2017-82216-R (HIDROIBERNIEVE). Francisco Navarro-Serrano is the recipient of a pre-doctoral FPU grant (Spanish Ministry of Education, Culture and Sports). All the simulation data are freely downloadable at https://zenodo.org/record/854619.

Conflicts of Interest

The authors declare no conflict of interest in this article

References

  1. Hurrell, J.W. Decadal trends in the North Atlantic oscillation: Regional temperatures and precipitation. Science 1995, 269, 676–679. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Hurrell, J.W.; Deser, C. North Atlantic climate variability: The role of the North Atlantic Oscillation. J. Mar. Syst. 2010, 79, 231–244. [Google Scholar] [CrossRef]
  3. Kjellström, E. North atlantic oscillation (Nao). In Encyclopedia of Earth Sciences Series; Springer Dordrecht: Heidelberg, Germany, 2016. [Google Scholar]
  4. Wanner, H.; Brönnimann, S.; Casty, C.; Gyalistras, D.; Luterbacher, J.; Schmutz, C.; Stephenson, D.B.; Xoplaki, E. North Atlantic Oscillation—Concepts and Studies. Surv. Geophys. 2001, 22, 321–381. [Google Scholar] [CrossRef]
  5. Trigo, R.; Osborn, T.; Corte-Real, J. The North Atlantic Oscillation influence on Europe: Climate impacts and associated physical mechanisms. Clim. Res. Clim. Res. 2002, 20, 9–17. [Google Scholar] [CrossRef]
  6. Sánchez-López, G.; Hernández, A.; Pla-Rabes, S.; Toro, M.; Granados, I.; Sigró, J.; Trigo, R.M.; Rubio-Inglés, M.J.; Camarero, L.; Valero-Garcés, B.; et al. The effects of the NAO on the ice phenology of Spanish alpine lakes. Clim. Chang. 2015, 130, 101–113. [Google Scholar] [CrossRef] [Green Version]
  7. López-Moreno, J.I.; Vicente-Serrano, S.M.; Morán-Tejeda, E.; Lorenzo-Lacruz, J.; Kenawy, A.; Beniston, M. Effects of the North Atlantic Oscillation (NAO) on combined temperature and precipitation winter modes in the Mediterranean mountains: Observed relationships and projections for the 21st century. Glob. Planet. Chang. 2011, 77, 62–76. [Google Scholar] [CrossRef]
  8. Durán, L.; Rodríguez-Fonseca, B.; Yagüe, C.; Sánchez, E. Water vapour flux patterns and precipitation at Sierra de Guadarrama mountain range (Spain). Int. J. Climatol. 2015, 35, 1593–1610. [Google Scholar] [CrossRef]
  9. López-Moreno, J.I.; Vicente-Serrano, S.M. Atmospheric circulation influence on the interannual variability of snow pack in the Spanish Pyrenees during the second half of the 20th century. Nord. Hydrol. 2007, 38, 33–44. [Google Scholar] [CrossRef] [Green Version]
  10. Revuelto, J.; López-Moreno, J.I.; Morán-Tejeda, E.; Fassnacht, S.R.; Vicente Serrano, S.M. Variabilidad Interanual del Manto de Nieve en el Pirineo: Tendencias Observadas y su Relación con Índices de Teleconexión Durante el Periodo 1985–2011; Asociación Española de Climatología: Barcelona, Spain, 2012; Available online: https://fundacion.usal.es/es/congresos-y-reuniones-cientificas1111/45-congresos-y-jornadas-actuales/472-8o-congreso-internacional-cambio-climatico-extremos-e-impactos (accessed on 28 December 2019).
  11. Bensiton, M. Variations of snow depth and duration in the Swiss Alps over the last 50 years: Links to changes in large-scale climatic forcings. Clim. Chang. 1997, 36, 281–300. [Google Scholar] [CrossRef]
  12. Scherrer, S.C.; Appenzeller, C. Swiss Alpine snow pack variability. Clim. Res. 2006, 32, 187–199. [Google Scholar] [CrossRef] [Green Version]
  13. Schöner, W.; Koch, R.; Matulla, C.; Marty, C.; Tilg, A.-M. Spatiotemporal patterns of snow depth within the Swiss-Austrian Alps for the past half century (1961 to 2012) and linkages to climate change. Int. J. Climatol. 2019, 39, 1589–1603. [Google Scholar] [CrossRef]
  14. Theakstone, W.H. Long-term variations of the seasonal snow cover in Nordland, Norway: The influence of the North Atlantic Oscillation. Ann. Glaciol. 2013, 54, 25–34. [Google Scholar] [CrossRef] [Green Version]
  15. Bednorz, E. Snow cover in western Poland and macro-scale circulation conditions. Int. J. Climatol. 2002, 22, 533–541. [Google Scholar] [CrossRef]
  16. Brown, R.D.; Petkova, N. Snow cover variability in Bulgarian mountainous regions, 1931–2000. Int. J. Climatol. 2007, 27, 1215–1229. [Google Scholar] [CrossRef]
  17. Bednorz, E.; Wibig, J. Spatial distribution and synoptic conditions of snow accumulation in the Russian Arctic. Polar Res. 2016, 35, 25916. [Google Scholar] [CrossRef] [Green Version]
  18. Henderson, G.R.; Leathers, D.J. European snow cover extent variability and associations with atmospheric forcings. Int. J. Climatol. 2009, 30, 1440–1451. [Google Scholar] [CrossRef]
  19. Keylock, C.J. The North Atlantic Oscillation and snow avalanching in Iceland. Geophys. Res. Lett. 2003, 30. [Google Scholar] [CrossRef]
  20. Jomelli, V.; Delval, C.; Grancher, D.; Escande, S.; Brunstein, D.; Hétu, B.; Filion, L.; Pech, P. Probabilistic analysis of recent snow avalanche activity and weather in the French Alps. Cold Reg. Sci. Technol. 2007, 47, 180–192. [Google Scholar] [CrossRef]
  21. García, C.; Martí, G.; Oller, P.; Moner, I.; Gavaldà, J.; Martínez, P.; Carlos Peña, J. Major avalanches occurrence at regional scale and related atmospheric circulation patterns in the Eastern Pyrenees. Cold Reg. Sci. Technol. 2009, 59, 106–118. [Google Scholar] [CrossRef]
  22. Cohen, J.; Fletcher, C. Improved Skill of Northern Hemisphere Winter Surface Temperature Predictions Based on Land–Atmosphere Fall Anomalies. J. Clim. 2007, 20, 4118–4132. [Google Scholar] [CrossRef]
  23. Scaife, A.A.; Arribas, A.; Blockley, E.; Brookshaw, A.; Clark, R.T.; Dunstone, N.; Eade, R.; Fereday, D.; Folland, C.K.; Gordon, M.; et al. Skillful long-range prediction of European and North American winters. Geophys. Res. Lett. 2014, 41, 2514–2519. [Google Scholar] [CrossRef] [Green Version]
  24. Wang, L.; Ting, M.; Kushner, P.J. A robust empirical seasonal prediction of winter NAO and surface climate. Sci. Rep. 2017, 7, 279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Beniston, M. Climatic Change in Mountain Regions: A Review of Possible Impacts. Clim. Chang. 2003, 59, 5–31. [Google Scholar] [CrossRef]
  26. Fayad, A.; Gascoin, S.; Faour, G.; López-Moreno, J.I.; Drapeau, L.; Le Page, M.; Escadafal, R. Snow hydrology in Mediterranean mountain regions: A review. J. Hydrol. 2017, 551, 374–396. [Google Scholar] [CrossRef]
  27. Buisan, S.T.; López-Moreno, J.I.; Saz, M.A.; Kochendorfer, J. Impact of weather type variability on winter precipitation, temperature and annual snowpack in the Spanish Pyrenees. Clim. Res. 2016, 69, 79–92. [Google Scholar] [CrossRef]
  28. Navarro-Serrano, F.; López-Moreno, J.I. Spatio-temporal analysis of snowfall events in the Spanish pyrenees and their relationship to atmospheric circulation | Análisis espacio-temporal de los eventos de nevadas en el pirineo Español y su relación con la circulación atmosférica. Cuad. Investig. Geogr. 2017, 43, 233–254. [Google Scholar] [CrossRef] [Green Version]
  29. Riaz, S.M.F.; Iqbal, M.J.; Hameed, S. Impact of the North Atlantic Oscillation on winter climate of Germany. Tellus A Dyn. Meteorol. Oceanogr. 2017, 69, 1406263. [Google Scholar] [CrossRef] [Green Version]
  30. Morán-Tejeda, E.; López-Moreno, J.I.; Beniston, M. The changing roles of temperature and precipitation on snowpack variability in Switzerland as a function of altitude. Geophys. Res. Lett. 2013, 40, 2131–2136. [Google Scholar] [CrossRef] [Green Version]
  31. Lorenzo-Lacruz, J.; Vicente-Serrano, S.M.; López-Moreno, J.I.; González-Hidalgo, J.C.; Morán-Tejeda, E. The response of Iberian rivers to the North Atlantic Oscillation. Hydrol. Earth Syst. Sci. 2011, 15, 2581–2597. [Google Scholar] [CrossRef] [Green Version]
  32. Morán-Tejeda, E.; Lorenzo-Lacruz, J.; López-Moreno, J.I.; Rahman, K.; Beniston, M. Streamflow timing of mountain rivers in Spain: Recent changes and future projections. J. Hydrol. 2014, 517, 1114–1127. [Google Scholar] [CrossRef] [Green Version]
  33. López-Moreno, J.I.; Beniston, M.; García-Ruiz, J.M. Environmental change and water management in the Pyrenees: Facts and future perspectives for Mediterranean mountains. Glob. Planet. Chang. 2008, 61, 300–312. [Google Scholar] [CrossRef] [Green Version]
  34. Gilaberte-Búrdalo, M.; López-Martín, F.; Pino-Otín, M.R.; López-Moreno, J.I. Impacts of climate change on ski industry. Environ. Sci. Policy 2014, 44, 51–61. [Google Scholar] [CrossRef] [Green Version]
  35. Alonso-González, E.; López-Moreno, J.I.; Gascoin, S.; García-Valdecasas Ojeda, M.; Sanmiguel-Vallelado, A.; Navarro-Serrano, F.; Revuelto, J.; Ceballos, A.; Esteban-Parra, M.J.; Essery, R. Daily gridded datasets of snow depth and snow water equivalent for the Iberian Peninsula from 1980 to 2014. Earth Syst. Sci. Data 2018, 10, 303–315. [Google Scholar] [CrossRef] [Green Version]
  36. Alonso-González, E.; López-Moreno, J.I.; Navarro-Serrano, F.; Sanmiguel-Vallelado, A.; Revuelto, J.; Domínguez-Castro, F.; Ceballos, A. Snow climatology for the mountains in the Iberian Peninsula using satellite imagery and simulations with dynamically downscaled reanalysis data. Int. J. Climatol. 2019. [Google Scholar] [CrossRef]
  37. Osborn, T.J.; Briffa, K.R.; Tett, S.F.B.; Jones, P.D.; Trigo, R.M. Evaluation of the North Atlantic Oscillation as simulated by a coupled climate model. Clim. Dyn. 1999, 15, 685–702. [Google Scholar] [CrossRef]
  38. Trigo, R.M.; Zêzere, J.L.; Rodrigues, M.L.; Trigo, I.F. The Influence of the North Atlantic Oscillation on Rainfall Triggering of Landslides near Lisbon. Nat. Hazards 2005, 36, 331–354. [Google Scholar] [CrossRef]
  39. Fernández-González, S.; Del Río, S.; Castro, A.; Penas, A.; Fernández-Raga, M.; Calvo, A.I.; Fraile, R. Connection between NAO, weather types and precipitation in León, Spain (1948–2008). Int. J. Climatol. 2012, 32, 2181–2196. [Google Scholar] [CrossRef] [Green Version]
  40. Skamarock, W.C.; Klemp, J.B.; Dudhia, J.; Gill, D.O.; Barker, D.M.; Dudha, M.G.; Huang, X.; Wang, W.; Powers, Y. A Description of the Advanced Research WRF Version 3. NCAR Tech. Note 2008. [Google Scholar] [CrossRef]
  41. Berrisford, P.; Dee, D.P.; Poli, P.; Brugge, R.; Fielding, K.; Fuentes, M.; Kallberg, P.W.; Kobayashi, S.; Uppala, S.; Simmons, A. The ERA-Interim archive Version 2.0. ERA Rep. Ser. 2011, 23, 13177. [Google Scholar]
  42. Dee, D.P.; Uppala, S.M.; Simmons, A.J.; Berrisford, P.; Poli, P.; Kobayashi, S.; Andrae, U.; Balmaseda, M.A.; Balsamo, G.; Bauer, P.; et al. The ERA-Interim reanalysis: Configuration and performance of the data assimilation system. Q. J. R. Meteorol. Soc. 2011, 137, 553–597. [Google Scholar] [CrossRef]
  43. García-Valdecasas Ojeda, M.; Gámiz-Fortis, S.R.; Castro-Díez, Y.; Esteban-Parra, M.J. Evaluation of WRF capability to detect dry and wet periods in Spain using drought indices. J. Geophys. Res. Atmos. 2017, 122, 1569–1594. [Google Scholar] [CrossRef]
  44. Essery, R. A factorial snowpack model (FSM 1.0). Geosci. Model Dev. 2015, 8, 3867–3876. [Google Scholar] [CrossRef] [Green Version]
  45. Siegel, S.; Castellan, N.J. Nonparametric Statistics for the Behavioral Sciences; McGraw-Hill: New York, NY, USA, 1988. [Google Scholar]
  46. Revuelto, J.; López-Moreno, J.I.; Azorin-Molina, C.; Vicente-Serrano, S.M. Topographic control of snowpack distribution in a small catchment in the central Spanish Pyrenees: Intra- and inter-annual persistence. Cryosphere 2014, 8, 1989–2006. [Google Scholar] [CrossRef] [Green Version]
  47. Koenig, W.D. Spatial autocorrelation of ecological phenomena. Trends Ecol. Evol. 1999, 14, 22–26. [Google Scholar] [CrossRef]
  48. López-Moreno, J.I.; Vicente-Serrano, S.M.; Morán-Tejeda, E.; Lorenzo-Lacruz, J.; Zabalza, J.; El Kenawy, A.; Beniston, M. Influence of Winter North Atlantic Oscillation Index (NAO) on Climate and Snow Accumulation in the Mediterranean Mountains. In Hydrological, Socioeconomic and Ecological Impacts of the North Atlantic Oscillation in the Mediterranean Region; Springer: Dordrecht, The Netherlands, 2011; pp. 73–89. [Google Scholar]
  49. Sáenz, J.; Zubillaga, J.; Rodríguez-Puebla, C. Interannual winter temperature variability in the north of the Iberian Peninsula. Clim. Res. 2001, 16, 169–179. [Google Scholar] [CrossRef] [Green Version]
  50. Lorente-Plazas, R.; Montávez, J.P.; Jimenez, P.A.; Jerez, S.; Gómez-Navarro, J.J.; García-Valero, J.A.; Jimenez-Guerrero, P. Characterization of surface winds over the Iberian Peninsula. Int. J. Climatol. 2015, 35, 1007–1026. [Google Scholar] [CrossRef]
  51. Stefanicki, G.; Talkner, P.; Weber, R.O. Frequency Changes of Weather Types in the Alpine Region since 1945. Theor. Appl. Climatol. 1998, 60, 47–61. [Google Scholar] [CrossRef]
  52. Del Río, S.; Fraile, R.; Herrero, L.; Penas, A. Analysis of recent trends in mean maximum and minimum temperatures in a region of the NW of Spain (Castilla y León). Theor. Appl. Climatol. 2007, 90, 1–12. [Google Scholar] [CrossRef]
  53. Esteban-Parra, M.J.; Pozo-Vázquez, D.; Rodrigo, F.S.; Castro-Díez, Y. Temperature and Precipitation Variability and Trends in Northern Spain in the Context of the Iberian Peninsula Climate. In Mediterranean Climate; Springer Berlin Heidelberg: Berlin/Heidelberg, Germany, 2003; pp. 259–276. [Google Scholar]
  54. López-Moreno, J.I. Recent variations of snowpack depth in the central Spanish Pyrenees. Arct. Antarct. Alp. Res. 2005, 37, 253–260. [Google Scholar] [CrossRef]
  55. Ríos-Cornejo, D.; Penas, Á.; Álvarez-Esteban, R.; del Río, S. Links between teleconnection patterns and precipitation in Spain. Atmos. Res. 2015, 156, 14–28. [Google Scholar] [CrossRef]
  56. Rodrigo, F.S.; Esteban-Parra, M.J.; Pozo-Vázquez, D.; Castro-Díez, Y. Rainfall variability in southern Spain on decadal to centennial time scales. Int. J. Climatol. 2000, 20, 721–732. [Google Scholar] [CrossRef] [Green Version]
  57. Muñoz-Díaz, D.; Rodrigo, F.S. Impacts of the North Atlantic Oscillation on the probability of dry and wet winters in Spain. Clim. Res. 2004, 27, 33–43. [Google Scholar] [CrossRef] [Green Version]
  58. Queralt, S.; Hernández, E.; Barriopedro, D.; Gallego, D.; Ribera, P.; Casanova, C. North Atlantic Oscillation influence and weather types associated with winter total and extreme precipitation events in Spain. Atmos. Res. 2009, 94, 675–683. [Google Scholar] [CrossRef]
  59. Vada, J.A.; Rodríguez-Marcos, J.; Buisán, S.; San Ambrosio, I. Climatological comparison of 2011–2012 and 2012–2013 snow seasons in Central and Western Spanish Pyrenees and its relationship with the North Atlantic Oscillation (NAO). In Proceedings of the International Snow Science Workshop, Chamonix Mont-Blanc, France, 7–11 October 2013. [Google Scholar]
  60. Buisan, S.T.; Saz, M.A.; López-Moreno, J.I. Spatial and temporal variability of winter snow and precipitation days in the western and central Spanish Pyrenees. Int. J. Climatol. 2015, 35, 259–274. [Google Scholar] [CrossRef]
  61. Musselman, K.N.; Clark, M.P.; Liu, C.; Ikeda, K.; Rasmussen, R. Slower snowmelt in a warmer world. Nat. Clim. Chang. 2017, 7, 214–219. [Google Scholar] [CrossRef]
  62. Morán-Tejeda, E.; López-Moreno, J.I.; Sanmiguel-Vallelado, A. Changes in Climate, Snow and Water Resources in the Spanish Pyrenees: Observations and Projections in a Warming Climate. In High Mountain Conservation in a Changing World; Springer: Cham, Switzerland, 2017; pp. 305–323. [Google Scholar]
  63. Dunstone, N.; Smith, D.; Scaife, A.; Hermanson, L.; Eade, R.; Robinson, N.; Andrews, M.; Knight, J. Skilful predictions of the winter North Atlantic Oscillation one year ahead. Nat. Geosci. 2016, 9, 809–814. [Google Scholar] [CrossRef]
Figure 1. Geographical locations of the main mountain ranges of the Iberian Peninsula.
Figure 1. Geographical locations of the main mountain ranges of the Iberian Peninsula.
Water 12 00105 g001
Figure 2. Pearson’s R values between the DJFM-NAO and minimum temperatures (A), maximum temperatures (B), and precipitation (C). The lines represent the boundaries of the main catchments in the Iberian Peninsula.
Figure 2. Pearson’s R values between the DJFM-NAO and minimum temperatures (A), maximum temperatures (B), and precipitation (C). The lines represent the boundaries of the main catchments in the Iberian Peninsula.
Water 12 00105 g002
Figure 3. Pearson’s R values between the DJFM-NAO and the snowpack duration. The lines represent the boundaries of the main catchments in the Iberian Peninsula.
Figure 3. Pearson’s R values between the DJFM-NAO and the snowpack duration. The lines represent the boundaries of the main catchments in the Iberian Peninsula.
Water 12 00105 g003
Figure 4. Absolute differences in the Pearson’s R statistic for correlations of the DJFM-NAO with the snow duration and the peak SWE. The colors in the scatterplots indicate the relative density of points, ranging from blue (lower relative density) to red (higher relative density).
Figure 4. Absolute differences in the Pearson’s R statistic for correlations of the DJFM-NAO with the snow duration and the peak SWE. The colors in the scatterplots indicate the relative density of points, ranging from blue (lower relative density) to red (higher relative density).
Water 12 00105 g004
Figure 5. Influence of elevation on the correlation of the DJFM-NAO with snow duration and the peak SWE. Red lines indicate no statistical significance of the correlations of the random samples (p-value > 0.05), while green lines show the statistically significant correlations.
Figure 5. Influence of elevation on the correlation of the DJFM-NAO with snow duration and the peak SWE. Red lines indicate no statistical significance of the correlations of the random samples (p-value > 0.05), while green lines show the statistically significant correlations.
Water 12 00105 g005
Figure 6. Influence of geographical longitude in the distribution of correlation values for the relationship between the DJFM-NAO and the snow duration and the peak SWE. Red lines indicate no statistical significance of the correlations of the random samples (p-value > 0.05), while green lines show the statistically significant correlations.
Figure 6. Influence of geographical longitude in the distribution of correlation values for the relationship between the DJFM-NAO and the snow duration and the peak SWE. Red lines indicate no statistical significance of the correlations of the random samples (p-value > 0.05), while green lines show the statistically significant correlations.
Water 12 00105 g006

Share and Cite

MDPI and ACS Style

Alonso-González, E.; López-Moreno, J.I.; Navarro-Serrano, F.M.; Revuelto, J. Impact of North Atlantic Oscillation on the Snowpack in Iberian Peninsula Mountains. Water 2020, 12, 105. https://doi.org/10.3390/w12010105

AMA Style

Alonso-González E, López-Moreno JI, Navarro-Serrano FM, Revuelto J. Impact of North Atlantic Oscillation on the Snowpack in Iberian Peninsula Mountains. Water. 2020; 12(1):105. https://doi.org/10.3390/w12010105

Chicago/Turabian Style

Alonso-González, Esteban, Juan I. López-Moreno, Francisco M. Navarro-Serrano, and Jesús Revuelto. 2020. "Impact of North Atlantic Oscillation on the Snowpack in Iberian Peninsula Mountains" Water 12, no. 1: 105. https://doi.org/10.3390/w12010105

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop