Spanish Extreme Winds and Their Relationships with Atlantic Large-Scale Atmospheric Patterns

Alvaro Pascual, Francisco Valero, Maria Luisa Martín, Carlos García-Legaz

Agencia Estatal de Meteorología, C/Leonardo Prieto Castro, Madrid, Spain.
Dpto. Astrofísica y CC. Física de la Atmósfera, Facultad de CC Físicas, Universidad Complutense de Madrid, Ciudad Universitaria, Madrid, Spain.
Dpto. Matemática Aplicada, Escuela Universitaria de Informática, Universidad de Valladolid, Segovia, Spain.

DOI: 10.4236/ajcc.2013.23A003                

Abstract

The purpose of this work is to review procedures to obtain relationships between wind and large-scale atmospheric fields, with special emphasis on extreme situation results. Such relationships are obtained by using different methods and techniques such as wind cumulative probability functions and composite maps. The analyses showed different mean atmospheric situations associated with the different wind patterns, in which strong atmospheric gradients can be related to moderate to strong winds in Spain. Additionally, a statistical downscaling analog model, developed by the authors, is used for diagnosing large-scale atmospheric circulation patterns and subsequently estimating extreme wind probabilities. From an atmospheric circulation pattern set obtained by multivariate methodology applied to a large-scale atmospheric circulation field, estimations of wind fields, particularly extreme winds, are obtained by means of the analogs methodology. Deterministic and probabilistic results show that gust behaviour is quite better approximated than mean wind speed, in general. The model presents some underestimations except for strong winds. Moreover, the model shows better probabilistic wind results over the Spanish northern area, highlighting that the atmospheric situations coming from the Atlantic Ocean are better recovered to predict mean wind and gusts in the Northern Peninsula.

Keywords

Extreme Winds; Gusts; Atmospheric Circulation Variability; Analogs; Probabilistic Results

1. Introduction

Winter storms are responsible for more than 50% of the total economic loss in central Europe, due to natural hazards [1,2] and a single extreme storm event can cause economic losses exceeding 10 billion euros. A rise in storm-related monetary losses for Europe in the course of the 20th century has been observed by Barredo [3], explained principally by changes in economic and demographic conditions, with much of the recent infrastructure in various parts of the world increasingly constructed in zones at risk from severe weather [2]. Therefore, the knowledge of atmospheric circulation patterns, particularly the one dealing with atmospheric patterns conducive to risky meteorological situations related to extreme wind events, is especially important for wind energy applications [4-6]. Forecasters, energy producers and grid operators have different views on what extremes related to wind generation are. Extreme events have been categorized taken into account damages and economic loss

[7,8]. Extremes have been also identified as their occurrence probabilities [9]; they have been analyzed from the spatial-temporal characteristics of prediction errors [10], or taking into account their probabilistic forecasts by statistical scenarios [11,12], by ensemble predictions [13]. The final analysis can be used for nowcasting wind power over the whole area, and for data assimilation purposes (in order to update and improve wind power predictions) for better understanding the, or for issuing “global” warnings related to expected accuracy of weather and wind power forecasts over the area considered.

In determining temporal-spatial distribution changes of wind and other climatological elements, it is necessary to take into account the atmospheric circulation variability. The Western European climate steps are necessary on the available knowledge of natural variability in regional scales and its relationship to large-scale circulation [14- 20]. The relative location of different pressure centers over the North Atlantic area influence different air masses with distinct physical characteristics over Iberia to produce a wide range of differentiated regional climates, playing the topography a leading role. In fact, at local scales the development of cloud systems or the enhancement of wind speed over different areas can be especially affected due to the topography [21-23]; at largescale domains, topography can generate o redirected synoptic and mesoscale flows [24]. The present study is firstly focused on showing relationships between wind and large-scale atmospheric fields over an Atlantic area, with special emphasis on results involving extreme situations. These connections are attained by using different methods and procedures, such us cumulative probability curves and composite maps. Composites have already been used by the authors in several studies in order to analyse different fields, obtaining relationships between them, so that maximum and minimum intensity phases of a field can be related to the other one [16,18,19].

On the other hand, the improvement of meteorological forecasts of wind by means of dynamic modelling has been progressing by means of limited area models or ensemble prediction systems in several research projects (ANEMOS, ANEMOS.plus). However, this methodology bears high computational costs. In order to overcome this problem, the analog method for predicting time series can be used [25]. With this method, local prediction models are obtained finding in a set of historic data similar situations to a particular situation [26-28]. This technique has been implemented for both climatic anomaly predictions [29,30] and short-range prediction [31], revealing as an alternative to other more complex models with high computational cost. In the framework of the European Project SafeWind, the authors have been developing several works based on multivariate methodologies for obtaining atmospheric situations analog to a situation associated with extreme winds [32]. One of the final purposes of this European Project is to develop a statistical downscaling model (ANPAF: ANalog PAttern Finder) for diagnosing large-scale atmospheric circulation patterns and subsequently estimating extreme wind probabilities. In the present paper, from an atmospheric circulation pattern set obtained by multivariate methodology [18,19,33,34] applied to a large-scale atmospheric circulation field, estimations of wind fields, particularly extreme winds, are obtained by means of the analogs methodology.

The study is organized as follows. In Section 2, data used in the study are described. In Section 3, the connections between wind speeds and large-scale atmospheric patterns are shown, presenting the relationships between large-scale atmospheric patterns and wind speed patterns statistically obtained. Moreover, the interactions between observational winds and large-scale atmospheric circulation statistical modes are provided and analyzed. Section 4 is devoted to analyze the analog results for both the large-scale atmospheric field and the Spanish mean wind speed and wind gust, in terms of some deterministic and probabilistic tools. The main conclusions are drawn in Section 5.

2. Data

In order to analyze the relationships between wind speeds and large-scale atmospheric fields and to extract information about extreme situations it is very important to select the appropriate datasets. In this work, in order to characterize the atmospheric circulation, 1000 hPa daily geopotential heights at 12:00 UTC (Z1000) for 36 winters from 1971 to 2007 covering from 51.5˚W to 15.5˚E and 20˚ to 60˚N have been used. Z1000 data are a product of the ERA40 Reanalysis [35]. Concerning the wind speed, firstly daily mean wind speed (MWS) data for 21 stations distributed over Spain (Figure 1) during the winter (D-J-F) season from 1970 to 2002 have been considered. These wind data come from in-situ measurements of the station network of the Spanish Meteorological Service (Agencia Estatal de Meteorología, AEMET). In order to analyze the relationships between wind speeds and large-scale atmospheric fields and to extract information about extreme situations it is very important to select the appropriate datasets. In this work, in order to characterize the atmospheric circulation, 1000 hPa daily geopotential heights at 12:00 UTC (Z1000) for 36 winters from 1971 to 2007 covering from 51.5˚W to 15.5˚E and 20˚ to 60˚N have been used. Z1000 data are a product of the ERA40 Reanalysis [35]. Concerning the wind speed, firstly daily mean wind speed (MWS) data for 21 stations distributed over Spain (Figure 1) during the winter (D-J-F) season from 1970 to 2002 have been considered. These wind data come from in-situ measurements of the station network of the Spanish Meteorological Service (Agencia Estatal de Meteorología, AEMET).

On the other hand, some techniques for estimating and forecasting wind speeds are reviewed, with special emphasis in extreme winds. To do this wind speed and wind gust estimations in Spain, with special emphasis in extreme values, are obtained using the analog methodology applied to the Z1000 data base. To do this, a additional data set is considered, the daily wind gust (WGU) data over Spain. The WGU data used in this paper consist of 73 time series of daily gusts in Spain (Figure 1). Taking into account the observational data quality and the methodology employed in this contribution, in this part of study described in Section 4, the three datasets finally cover the common period from 1971 to 2002.

Figure 1. The Iberian Peninsula with its orography detailed with the mean wind speed stations in red crosses and the wind gust stations in black crosses.

3. Extreme Wind Speeds-Large Scale Atmospheric Patterns Connections

3.1. Large-Scale Atmospheric Patterns—Wind Speed Statistical Mode Relationships

A Principal Component Analysis (PCA) is applied to the MWS and Z1000 fields in order to know its general behaviour and to extract the most significant patterns from the original data [36]. However beyond mere data compression, a PCA is a very useful tool for exploring large multivariate data sets because of its potential for yielding substantial insights into both the spatial and temporal variations of the analysed fields. This methodology applied to spatial data enables patterns to be identified that can be attributed to specific physical processes by statistical assessment. The new uncorrelated variables are called principal components (PCs) and consist of linear combinations of the original variables derived from the diagonalization of the covariance/correlation matrix. The coefficients of the linear combinations represent the weight of the original variables in the PCs and they are named loadings or PC patterns. The PCs indicate modes of variation of the original field and are numbered according with their related variance. Thus, the first PC is the linear combination with the maximum possible variance; the second one is the linear combination with the maximum possible variance which is uncorrelated with the first PC and so on. The projection of the original series onto each eigenvector gives as result the time-dependence coefficient named scores or PC time series. In our case, the PCA was applied to the correlation matrices of both data sets, the MWS and Z1000 fields, being a set of eigenvalues and eigenvectors produced for each data set. Generally, the most important (the first ones) eigenvectors tend to describe regions with largest fluctuations. Thus, most relevant information from the data can be represented using fewer numbers of the principal components and a much smaller data set. Five leading modes for both datasets have been selected (not all shown). They account for more than 66% and 77% of the total variability for MWS and Z1000, respectively.

For reasons of brevity only the first mode is shown. In Figure 2(a), the eigenvector or spatial pattern of the retained MWS PCs is shown which helps highlighting diverse areas of different wind behaviour over Spain. The leading wind PC pattern (Figure 2(a)) accounts for the most important percentage of variance in the original data (37.9%). In Figure 2(a) the spatial pattern shows homogeneous wind behaviour in inner Iberia, and also underlines the area to the North Iberian Plateau with high correlation values. This conduct in the wind field could be related to the predominant westerly circulation regime (Poniente) in the Iberian Peninsula. The time variability of the spatial pattern above described is depicted showing the evolution of its PC time series obtained by applying the PCA over the MWS data in wintertime (Figure 2(b)). Significant trends are not found after applying a Mann-Kendall test and a spectral analysis of the PC time series. As stated previously, the first spatial pattern of Figure 2a showed homogeneous wind pattern over Iberia, underlying areas corresponding to the North Iberian Plateau. This behaviour can also be represented in the corresponding time series (Figure 2(b)) with mostly positive and high score values over the selected period (1970-2002).

However, the derived modes are statistically obtained. To analyze the extreme situations is needed to find connections between wind speed and the atmospheric field.

(a)(b)(c)

Figure 2. (a) Spatial patterns of the first PC of wind speed. The positive (negative) correlations are solid (dashed); (b) Time series of the first PC of the wind speed. Units are in standard deviations in the y-axis and the x-axis corresponds to the time period; (c) Illustration of picked up dates from the PC time series. Red rows indicate both the 5% high positive and negative scores used to build the composite maps.

Thus, to examine the real atmospheric circulation features associated with the winter wind speed patterns a set of positive and negative composite plots (of Z1000 and MWS) was constructed from the dates associated with 5 and 95 percentiles of the scores of the time series obtained of the PCA (Figure 2(c)). The composite maps represent configurations of the variable which are comparable to observations. Composites are defined here as the averaged ensemble of sets of maps of the large-scale atmospheric variable and the wind speeds [37]. Physical distinctive features in the composite plots are achieved through obtaining additional information to the statistical meaning of the derived spatial modes. Here, the anomaly composites of large-scale atmospheric variables have been built for those weather configurations associated with the highest and lowest PC scores of the wind speed. This way, the composites represent the atmospheric state associated with particular extreme wind characteristics. Positive (negative) composites are constructed directly from a number of configurations with high (low) scores of the PC time series because they indicate situations in which the corresponding PC mode is dominant in its positive (negative) phase. The selected number of configurations represents 5% of the total number of cases in the dataset.

Figure 3 shows the anomaly composites for Z1000 displaying the positive and negative composite plots conditioned by the 5% highest and lowest PC scores of the MWS. Subsequently, mean maps of Z1000 anomalies are drawn up from these days, and highlight the mean atmospheric state conditioned by predominant oscillation of the selected wind speed PC mode. Additionally, maps of MWS, also corresponding to those days, are picked up to illustrate the behaviour of the wind speed field over Spain in such atmospheric situations. Thus, the Z1000 anomaly composites associated to the first wind speed PC (Figures 3(a) and (b) first positive and negative composites) highlight two different mean atmospheric situations associated with the wind behaviour. Thus, in the first positive anomaly composite (Figure 3(a)), a strong gradient of Z1000 is observed over the Iberian Peninsula, underlying strong winds over Iberia as it can be noted in Figure 3(c) with wind speeds exceeding 8 m∙s⁻¹ (30 km∙h⁻¹) in daily average. In contrast to this atmospheric situation, the first negative composite (Figure 3(b)) displays high anomaly pressure over Iberia with little gradient over it and a nucleus over northern France. This situation is indicative of low wind speed

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	“Swiss Re,” 2000. http://www.swissre.com/about_us/art_architecture/Swiss_Re_Next.html
[2]	U. Ulbrich, A. H. Fink, M. Klawa and J. G. Pinto, “Three extreme storms over Europe in December,” Weather, Vol. 56, No. 3, 2001, 1999, pp. 70-80. doi:10.1002/j.1477-8696.2001.tb06540.x
[3]	J. I. Barredo, “No Upward Trend in Normalised Windstorm Losses in Europe: 1970-2008,” Natural Hazards of Earth System Sciences, Vol. 10, 2010, pp. 97-104. doi:10.5194/nhess-10-97-2010
[4]	J. P. Palutikof, P. M. Kelly, T. D. Davies and J. A. Halliday, “Impacts of Spatial and Temporal Windspeed Variability on Wind Energy Output,” Journal of Applied Meteorology, Vol. 26, No. 9, 1987, pp. 1124-1133. doi:10.1175/1520-0450(1987)026<1124:IOSATW>2.0.CO;2
[5]	R. H. Thuilleier, “Real-Time of Local Wind Patterns for Application to Nuclear-Emergency Response,” Bulletin of the American Meteorological Society, Vol. 68, No. 9, 1987, pp. 1111-1115. doi:10.1175/1520-0477(1987)068<1111:RTAOLW>2.0.CO;2
[6]	J. A. Zuranski and B. Jaspinka, “Directional Analysis of Extreme Wind Speeds in Poland,” Journal of Wind Engineering and Industrial Aerodynamics, Vol. 65, No. 1-3, 1996, pp. 13-20. doi:10.1016/S0167-6105(97)00018-4
[7]	M. Gaya, J. Amaro, M. Aran and M. C. Llasat, “Preliminary Results of the Societal Impact Research Group of MEDEX: The Request Database (2000-2002) of Two Meteorological Services,” Proceedings of 9th EGS Plinius Conference, Nisosia, 2008, p. 12.
[8]	J. Amaro, M. Aran, L. Barberia and M. C. Llasat, “The Strong Wind Event of 24th January 2009 in Catalonia: A Social Impact Analysis,” Proceedings of 10th EGS Plinius Conference, Barcelona, 2009, p. 10.
[9]	M. J. OrtizBeviá, E. SánchezGómez and F. J. Alvarez-García, “North Atlantic Atmospheric Regimes and Winter Extremes,” Natural Hazards and Earth System Science, Vol. 11, 2011, pp. 971-980. doi:10.5194/nhess-11-971-2011
[10]	J. Tastu, P. Pinson, E. Kotwa, H. Aa. Nielsen and H. Madsen, “Spatio-Temporal Analysis and Modeling of Wind Power Forecast Errors,” Wind Energy, Vol. 14, No. 1, 2011, pp. 43-60. doi:10.1002/we.401
[11]	P. Pinson, G. Papaefthymiou, B. Klockl, H. Aa. Nielsen and H. Madsen, “From Probabilistic Forecasts to Statistical Scenarios of Short-Term Wind Power Production,” Wind Energy, Vol. 12, No. 1, 2009, pp. 51-62. doi:10.1002/we.284
[12]	P. Pinson and H. Madsen, “Adaptive Modeling and Forecasting of Wind Power Fluctuations with Markov-Switching Autoregressive Models,” Journal of Forecasting, Vol. 31, No. 4, 2012, pp. 281-313. doi:10.1002/for.1194
[13]	P. Pinson, H. Aa. Nielsen, H. Madsen and G. Kariniotakis, “Skill Forecasting from Ensemble Predictions of Wind Power,” Applied Energy, Vol. 86, No. 7-8, 2009, pp. 1326-1334. doi:10.1016/j.apenergy.2008.10.009
[14]	M. L. Martín, D. Santos-Munoz, F. Valero and A. Morata, “Evaluation of an Ensemble Precipitation Prediction System over the Western Mediterranean Area,” Atmospheric Research, Vol. 98, No. 1, 2010, pp. 163-175. doi:10.1016/j.atmosres.2010.07.002
[15]	M. Y. Luna, M. L. Martín, F. Valero and F. González-Rouco, “Wintertime Iberian Peninsula Precipitation Variability and Its Relation to North Atlantic Atmospheric Circulation,” In: M. Brunet and D. López, Eds., Detecting and Modelling Regional Climate Change and Associated Impacts, Springer-Verlag, Berlin, 2001, pp. 369-376. doi:10.1007/978-3-662-04313-4_31
[16]	F. Valero, M. Y. Luna, M. L. Martín, A. Morata and F. González-Rouco, “Coupled Modes of Large-Scale Climatic Variables and Regional Precipitation in the Western Mediterranean in Autumn,” Climate Dynamics, Vol. 22, No. 2-3, 2004, pp. 307-323. doi:10.1007/s00382-003-0382-9
[17]	M. L. Martín, M. Y. Luna, A. Morata and F. Valero, “North Atlantic Teleconnection Patterns of Low-Frequency Variability and Their Links with Springtime Precipitation in the Western Mediterranean,” International Journal of Climatology, Vol. 24, No. 2, 2004, pp. 213-230. doi:10.1002/joc.993
[18]	M. L. Martín, F. Valero, A. Morata, M. Y. Luna, A. Pascual and D. Santos-Munoz, “Springtime Coupled Modes of Regional Wind in the Iberian Peninsula and Large-Scale Variability Patterns,” International Journal of Climatology, Vol. 31, No. 6, 2011, pp. 880-895. doi:10.1002/joc.2127
[19]	M. L. Martín, F. Valero, A. Pascual, A. Morata and M. Y. Luna, “Springtime Connections between the Large-Scale Sea Level Pressure Field and Gust Wind Speed over Iberia,” Natural Hazards of Earth System Sciences, Vol. 11, 2011, pp. 191-203. doi:10.5194/nhess-11-191-2011
[20]	E. García-Ortega, L. López and J. L. Sánchez, “Atmospheric Patterns Associated with Hailstorm Days in the Ebro Valley, Spain,” Atmospheric Research, Vol. 100, No. 4, 2011, pp. 401-427. doi:10.1016/j.atmosres.2010.08.023
[21]	R. Frouin, A. F. Fiúza, I. Ambar and T. J. Boyd, “Observations of a Poleward Surface Current off the Coasts of Portugal and Spain during the Winter,” Journal of Geophysical Research, Vol. 95, No. C1, 1990, pp. 679-691. doi:10.1029/JC095iC01p00679
[22]	R. D. Haynes and E. D. Barton, “A Poleward Flow along the Atlantic Coast of the Iberian Peninsula,” Journal of Geophysical Research, Vol. 95, No. 1, 1990, pp. 11425-11442. doi:10.1029/JC095iC07p11425
[23]	P. Bougeault, B. Benech, P. Bessemoulin, B. Carissimo, A. Lar, J. Pelon, M. Petitdidier and E. Richard, “PYREX: A Summary of findings,” Bulleting American Meteorological Society, Vol. 78, No. 4, 1997, pp. 637-650. doi:10.1175/1520-0477(1997)078<0637:PASOF>2.0.CO;2
[24]	F. Valero, Y. Luna and M. L. Martín, “An Overview of a Heavy Rain Event at Southeastern Iberia: The Role of the Large-Scale Meteorological Conditions,” Annales Geophysicae, Vol. 15, 1997, pp. 494-502. doi:10.1007/s00585-997-0494-3
[25]	E. N. Lorenz, “Atmospheric Predictability as Revealed by Naturally Accourring Analogues,” Journal of Atmospheric Science, Vol. 26, No. 4, 1969, pp. 636-646. doi:10.1175/1520-0469(1969)26<636:APARBN>2.0.CO;2
[26]	T. Hastie, R. Tibshirani and J. Friedman, “The Elements of Statistical Learning,” Springer, New York, 2001. doi:10.1007/978-0-387-21606-5
[27]	J. Fernandez and J. Saenz, “Improved Field Reconstruction with the Analog Method: Searching the CCA Space,” Climate Research, Vol. 24, No. 3, 2003, pp. 199-213. doi:10.3354/cr024199
[28]	K. Fraedrich, C. C. Raible and F. Sielmann, “Analog Ensemble Forecasting of Tropical Cyclone Tracks in the Australian Region,” Weather and Forecasting, Vol. 18, 2003, pp. 3-11.
[29]	E. Y. Zorita and H. von Storch, “The Analog Method as a Simple Statistical Downscaling Technique: Comparison with More Complicated Methods,” Journal of Climate, Vol. 12, No. 8, 1999, pp. 2474-2489. doi:10.1175/1520-0442(1999)012<2474:TAMAAS>2.0.CO;2
[30]	R. L. Wilby and T. Wigley, “Downscalling General Circulation Model Output. A Review of Methods and Limitations,” Progress in Physical Geography, Vol. 21, No. 4, 1997, pp. 530-548. doi:10.1177/030913339702100403
[31]	H. M. Dool van den, “Searching for Analogs, How Long Must We Wait?” Tellus, Vol. 46A, 1994, pp. 314-324.
[32]	A. Pascual, M. L. Martin, F. Valero, D. Santos-Munoz,, A. Morata and M. Y. Luna, “Development of an Analogous Model for Wind Prediction Using Principal Components,” The SAFEWIND Workshop, Oldemburg, 2010, 22 p. www.safewind.eu
[33]	A. Morata, M. L. Martín, M. G. Sotillo, F. Valero and M. Y. Luna, “Iberian Autumn Precipitation Characterization through Observed, Simulated and Reanalysed Data,” Advances of Geosciences, Vol. 16, 2008, pp. 49-54. doi:10.5194/adgeo-16-49-2008
[34]	F. Valero, M. L. Martín, M. G. Sotillo, A. Morata and M. Y. Luna, “Characterization of the Autumn Iberian Precipitation from Long-Term Data Sets: Comparison between Observed and Hindcasted Data,” International Journal of Climatology, Vol. 29, No. 4, 2009, pp. 527-541. doi:10.1002/joc.1726
[35]	A. J. Simmons and J. K. Gibson, “The ERA-40 Project Plan,” ERA-40 Project Report Series No. 1. ECMWF, Reading, 2000.
[36]	R. W. Preisendorfer, “Principal Component Analysis in Meteorology and Oceanography,” Elsevier Science Publishers BV, Amsterdam, 1998.
[37]	A. Pascual, M. L. Martín, F. Valero, M. Y. Luna and A. Morata, “Wintertime Connections between Extreme Wind Patterns in Spain and Large-Scale Geopotential Heath Field,” Atmospheric Research, Vol. 122, 2013, pp. 213-228. doi:10.1016/j.atmosres.2012.10.033
[38]	H. B. Bluestein, “Synoptic Dynamic Meteorology in Midlatitudes. Vol. II. Observations and Theory of Weather Systems,” Oxford University Press, Oxford, 1993.
[39]	T. M. Hamill, “Interpretation of Rank Histograms for Verifying Ensemble Forecasts,” Monthly Weather Review, Vol. 129, 2001, pp. 550-560. doi:10.1175/1520-0493(2001)129<0550:IORHFV>2.0.CO;2
[40]	P. Pinson and R. Hagedorn, “Verification of the ECMWF Ensemble Forecasts of Wind Speed against Observations,” Meteorological Applications, Vol. 19, No. 4, 2012, pp. 484-500. doi:10.1002/met.283
[41]	S. Herrera, S. Pazo, J. Fernández and M. A. Rodríguez, “The Role of Large-Scale Spatial Patterns in the Chaotic Amplification of Perturbations in a Lorenz’96 Model,” Tellus A, Vol. 63, No. 5, 2011, pp. 978-990. doi:10.1111/j.1600-0870.2011.00545.x
[42]	A. Pascual, F. Valero, M. L. Martín, A. Morata and M. Y. Luna, “Probabilistic and Deterministic Results of the ANPAF Analog Model for Spanish Wind Field Estimations,” Atmospheric Research, Vol. 108, 2012, pp. 39-56. doi:10.1016/j.atmosres.2012.01.011
[43]	A. S. Cofino, “Técnicas Estadísticas y Neuronales de Agrupamiento Adaptativo Para la Predicción Probabilística de Fenómenos Meteorológicos Locales. Aplicación en el Corto Plazo y en la Predicción Estacional,” Tesis Doctoral, Universidad de Cantabria, Cantabria, 2004.
[44]	J. M. Gutiérrez, R. Cano, A. S. Cofino and M. A. Rodríguez, “Clustering Methods for Statistical Down-scaling in Short-Range Weather Forecast,” Monthly Weather Review, Vol. 132, No. 9, 2004, pp. 2169-2183. doi:10.1175/1520-0493(2004)132<2169:CMFSDI>2.0.CO;2
[45]	G. W. Brier, “Verification of Forecasts Expressed in Terms of Probabilities,” Monthly Weather Review, Vol. 78, No. 1, 1950, pp. 1-3. doi:10.1175/1520-0493(1950)078<0001:VOFEIT>2.0.CO;2