Wind speed investigations in future period
1. Introduction
Green energy and renewable energy are one of the most essential fundamentals for future developments of countries all over the world. Running out of fossil fuels in near future, makes the use of renewable energies almost inevitable. Evaluation of the capability of using these types of energies is an essential issue at present. Extreme dependency of mankind to the limited energy resources is both dangers and unsustainable and finding a way out is one of the most important challenges that we face. Achieving some unlimited energy sources has been human’s dream. Daily increase of energy demands and limitation of fossil energy resources from one side, and increase of environmental pollutants caused by using these resources in the other side, has made application of renewable energies more essential and widespread. The wind, as one of the climatological factors has a wide effects on agriculture, transport, pollution, energy, manufacturing and industrial plans. Wind power is one of the first energy sources discovered by humans. It is applied for different purposes such as power source for ship movement, irrigation and milling in some countries such as Iran and China. However, construction of wind power plants would make this type of energy more applicable. Nowadays, coal, oil and gas are the main resources to provide energy. In recent years, an increasing trend in price of the mentioned materials has become vivid due to the globalization and political events and mostly because of water crises. These increasing prices make us to find an appropriate solution to decrease the expenses and increasing the stability. Furthermore, renewable energies are considered as available, exchangeable and inexhaustible resources. New energies are available until these resources exist. Wind energy is used in two ways, first in direct way in which the wind is applied for drying and ventilation, and secondly, the wind is utilized indirectly for milling the grains, to pump water to the fields and to generate electricity. The evidence shows that in some countries such as Iran, Iraq, Egypt, China, Italy, Spain, the wind energy has been used for milling and irrigation. According to IPCC special report on renewable energy sources and decreasing the climate change damages, there is an increasing trend in the magnitude of the wind power plants from 17 meters high and 75 kilowatts in 1980 decade to 80 meters high and 1800 kilowatts in 20052010. It is predicted that the size of these generators should be 250 meters high and 20000 kilowatts in future years which is clarifying the increasing trend of application of this green energy. An important challenge in front of this kind of power plants is the uncertainty in accessible capacity of electrical power. This problem has been caused by random nature of effective factors such as random variation in mechanical forces generating wind power. In other words, due to the continuous variations in meteorological and climatological conditions, the wind speed, duration, density and power are randomly changing. Thus, for using its power, it is necessary to study the windy conditions in the area and statistical data reported by meteorology centers. However, the analysis of such large amount of data recorded in meteorology centers to estimate the mechanical power input of wind power plants is not possible except using some applied methods. Obviously, the study of the behavior and speed of winds will lead to more accurate estimation of accessible capacity in wind power plants. Moreover, spatial analysis of this climatological phenomenon will provide some essential knowledge about the areas with potential capacity of constructing wind power plants. The frequency analysis is an operative tool in assessing this aim. This study is going to describe these materials and their effects on regionalization of wind speed under effect of climate change.
2. Previous investigations on wind prediction
There are some studies reported in the literature about wind speed in future periods which some of them are summarized in Table 1.






Historical wind time series for future application  Harmsen et al. (2009)  
Historical wind time series for future application  Doria et al. (2006)  
Wind speed change is stronger than the changes obtained for the same region on daily precipitation and temperature  Case Study  Nebraska  2 scenarios  MPI  Bogardi and Matyasovszky (1996) 
Dynamical Downscaling  Case study  German  2 scenarios  ECHAM  Hoyme and Zieleke (2001) 
The GCM models of climate change become more reliable and as tools are refined for improving results at regional scales, it will be desirable to include improved estimates of vulnerabilities in the wind power site selection decision process.  Case study  USA  2 scenarios  HADCM2, CGCM1  Breslow and Sailor (2002) 
ANN and statistical downscaling A2 is preferred to A1b changes in daily mean wind speeds at each location and they are presented and discussed with respect to potential implications for wind power generation. 
Case study  USA  A1b, A2  4 GCM  Sailor et al. (2000) 
Impact of the climate change scenarios on wind power may be as high as a 40% reduction in summer time generation potential.  Case study  Texas and California  2 scenarios  CCM  Sailor et al. (2008) 
In order to create the scenarios of changes in power density under climate change, a method is developed for mapping the dailyresolution downscaled GCM output to the hourly level. The model results for either of the SRES scenarios were similar. In most cases, the impacts of the A2 scenario were slightly larger than those for the A1B scenario, but these interscenario differences were smaller than the intermodel differences, even after downscaling. 

The average of wind speed will increase significantly towards the end of this century under the changing climate impacts.  Case study  Finland  A2  HWIND and SIMA  Peltola et al. (2010) 
 Use of NCEP and ANN for downscaling.  For the two locations considered, the increase in the 100year wind was found to be varying from 44% to 74%. 
Case study  Indian  A2  CGCM3  Deepthi and Deo (2010) 
By the end of the twentyfirst century, there is an evidence for small magnitude changes in the wind resource and increases in extreme wind speeds, and the declines in sea ice and icing frequencies.  Review Paper  Some changes associated with climate evolution will likely benefit the wind energy industry while other changes may negatively have an impact on wind energy developments, with such ‘gains and losses’ depending on the region under consideration.  Gumbel is presented for the probability distribution of extreme wind speeds 
Pryor and Barthelmie (2010) 

15 – 30 % wind speed growth  Case study  Brazil  A1B  HADCM3  Pereira et al. (2012) 
3. Climate change impacts
Fast development of industries and its outcome as increasing of the emission of greenhouse gases, has led to destruction of climatic equilibrium of the earth. This phenomenon is called “Climate Change” (IPCC 2007, Leander et al. 2006). The research is indicating the negative impacts of this phenomenon on different systems such as water resources, agriculture, environment, health, industry, and economy. The importance and hazardous of climate change has been emphasized in different international communities such as the group of eight (G8) which is a forum for the governments of eight of the world's largest economies and some of its facing solutions to save water resources, agriculture, and environmental resources have been suggested. As the water is an important resource, which is extremely under effect of climate change, the analysis of its changes in future years can provide a very useful key for future droughts, floods, evapotranspiration and etc.
The first step in the study of climate change impacts on future resources is to simulate the behavior of climatological factors under the effect of greenhouse gases. A general circulation model (GCM) is a three dimensional mathematical models of the general circulation of a planetary atmosphere or ocean. Atmospheric and oceanic GCMs (AGCM and OGCM) are key components of global climate models, which are systems of differential equations. Using such models, scientists divide the atmosphere, hydrosphere, geosphere, cryosphere, and biosphere of the planet into a 3dimensional grid, apply the basic equations, and evaluate the results. Atmospheric models calculate winds, heat transfer, radiation, relative humidity, and surface hydrology within each grid and evaluate the interactions with neighboring points. Different greenhouse gases emission scenarios such as A1, B1, A2, and B2 are going to be used during the simulation process.
4. Generation of climatic scenarios
The aforementioned 3dimensional joint atmosphericoceanic general circulation models (AOGCM) are used in this study among different available methods for the generation of climatic scenarios. The GCM models have a physical basis presented by mathematical relations. They are going to be solved in a 3dimensional grid all over the planet. In order to simulate the climate of the planet, the fundamental climatic processes in the atmosphere, hydrosphere, geosphere, cryosphere, and biosphere will be simulated in separate secondary models. Then, these atmospheric and oceanic secondary models are joining together to form AOGCMs. To study the condition of climate in the past periods, the observed values of greenhouse gases, solar radiation changes, and volcanic eruption aerosols until the 2000 are entered as input to the GCM models and the climatic variables are simulated as time series. After simulation of these variables in the past periods by using the GCM models, the introduction of future greenhouse gas conditions is necessary for simulation of these variables in future periods. For this purpose, at first the amounts of emitted greenhouse gases presented by emission scenarios (which are always until 2100) are transformed to concentrate and then to the amount of solar radiation and these values would be the input of the GCM models. The results obtained from the GCM models under emission scenarios will form the time series of climatic variables up to 2100.
5. Downscaling
One of the main challenges using the output of the AOGCM models is the spatial scale of their calculation cell and the downscaling method is used to solve this challenge. These methods are generally consists of two main groups of dynamic and statistical ones. In these methods, the downscaling procedure is done by using the observed meteorological data. A considerable point in application of the final outputs is the different sources of uncertainty, which can be evaluated using the Bootstrap method (Efron, 1993) in each confidence level.
6. Frequency analysis
The magnitude of an extreme event has an inverse relation with its frequency. In other words, the higher magnitude is the event, the less is its frequency of occurrence. The primary objective of frequency analysis is to relate the magnitude of extreme events to their frequency of occurrence through application of probability distributions (Chow et al., 1988). The first assumption in this manner is that under study data are independent and identically distributed and their underlying system is random and is spatially and temporally independent. This would be available when there is no correlation between observations. In application, these conditions can be achieved by using annual maximum values noting the independency of events between years. However, the wind speed parameter has been rarely examined by this method among other meteorological parameters; therefore, the studies in this field are at the beginning.
7. Probability distribution function and frequency formula
To describe the probability distribution of a random variable
This is the probability of the random variable
The occurrence of extreme events is not according to a constant regime or with a fixed magnitude and the time interval between two such events is variable. Thus, the return period defined as the average interarrival time between two extreme events is an applicable tool in such cases. An extreme event
Thus, the above mentioned nonexceedance probability can be presented as follows:
Equation (3) is presenting the magnitude of an extreme event correspond to a return period
which gives the value of
8. Selection of probability distribution
The probability distribution is a function for describing the probability of occurrence of a random event. Large amount of statistical information will be summarized in the distribution and its parameters by fitting a probability distribution on a set of hydrologic data. The most important and widely used methods for estimation of distribution parameters from data samples are method of moments, linear moments, and maximum likelihood which will be described more in next paragraph.
9. Method of moments
The method of moments was first introduced by Pearson (1902). He found that the appropriate estimations of the parameters of a probability distribution are those which their moments match with corresponding sample moments in the best way. In this method, general formula for calculation of moments of order
The method of moments is describing the relation between moments and distribution parameters. The most important moments around the mean are the mean, variance, skewness, and kurtosis, which are the one to four order moments, respectively.
10. Maximum likelihood method
Maximumlikelihood estimation was recommended, analyzed and vastly popularized by R. A. Fisher between 1912 and 1922 (Aldrich, John 1997). He argued that the best value of a parameter of a probability distribution should be one, which maximizes the likelihood or joint probability of occurrence of the sample. Assume that the sample space is divided into parts with length
As the distance
The maximum likelihood method is theoretically the most accurate method in estimation of the parameters of probability distributions. In fact, it is estimating those parameters with minimum average error with respect to correct parameters.
11. Lmoments method
Probability weighted moments (PWMs) are defined by Greenwood et al. (1979) as
Which can be rewritten as:
Where
The first four order Lmoments can be calculated as bellow:
Lmoment ratios, which are analogous to conventional moment ratios, are defined by Hosking (1990) as:
where
However, Lmoments have significant advantages over PWM's, specially their ability to summarize a statistical distribution in a more meaningful way. Since Lmoment estimators are linear functions of the ordered data values, they are virtually unbiased and less influenced by outliers. Also they have relatively small sampling variance and the bias of their small sample estimates remains quite small. Lmoments have become popular tools for solving various problems related to parameter estimation, distribution identification, and regionalization in different fields such as hydrology, water resources, and especially in regional analysis of rainfall and floods.
12. Atstation goodness of fit tests
After estimation of the parameter of the prescribed distributions, the usual question is the selection of the best fitted distribution to the observation sample. For this aim, the goodnessoffit tests are used to compare fitted theoretical distributions and observations. Two very common goodnessoffit tests with wide application in the literature are Chisquare and KolmogorovSmirnov tests. Another extensively used test is root mean square error (RMSE) which estimates the root of the square differences between observed values and calculated ones divided by the sample size. Among theoretical probability distributions, some of them are selected by researchers for describing wind speed, which are presented in Table 2.



1940 to 1945  Pearson Type III  Putnum^{*} 
1951  Pearson Type III  Sherlock^{*} 
1950 to 1970  Bivariate Distributions of Two Components  Essenwanger^{ *} 
TwoParameter Normal  Crutcher and Bear^{ *}  
1970s  Isotropic Gaussian Model of McWilliams et al.  Justus and Koeppl^{ *} 
1974, 1976, 1977  ThreeParameter Lognormal  Luna and Church^{ *} Kaminsky, Justus et al. 
1976 to 1977  Squareroot Normal Model  Winger^{ *} 
1978  Threeparameter Weibull  Stewart and Essenvanger^{ *} 
1980  ThreeParameter Generalized Gamma  Auwera^{ *} 
1980  Inverse Gaussian  Bardsley^{ *} 
1983  Pearson Type I (Beta)  Lavagnini et al. ^{*} 
1994  Weibull  Stelios Pasardes^{ *} 
1996  LogNormal  Bogardi and Matyasovski 
2009  Beta  Carta et al. 
2009  TwoComponent Weibull  Akdog et al.^{ *} 
2010  Gumbel and Weibull  Deepthi and Deo 
13. Regional frequency analysis
One of the main problems in frequency analysis is the lack of adequate long time data in locations under study, which beside insufficient accuracy in data recording in atsite estimates, has caused regional frequency analysis to be more applicable in such studies. Steps of regional frequency analysis are presented below.
14. Discordancy test
For the screening of similar sites, the discordancy measure, in terms of the sample Lmoment ratios (LCV, Lskew, and Lkurt.) of the gauging sites’ observed data is suggested by Hosking and Wallis (1997). The aim of the data screening performed using Lmoments based on a discordancy measure
Where N is the number of stations,
15. Heterogeneity test
Lmoment heterogeneity tests allow assessing whether a group of sites may reasonably be treated as a homogeneous region. The heterogeneity measure compares the betweensite variations in sample Lmoments for a group of sites with what would be expected for a homogeneous region (Hosking and Wallis, 1993).
The homogeneity test used in this study is the test that is proposed by Hosking and Wallis (1993) and is based on various orders of sample Lmoment ratios. It is particularly based on the variability of three different levels of tests: a test based on the LCV only; a test based on the LCV and Lskew; and a test based on the Lskew and Lkurt. These tests are called the Vstatistics and are respectively defined as:
The heterogeneity measure is then defined as:
where
where
Declare the region to be heterogeneous if
In order to achieve reliable estimates
Hosking and Wallis (1993) revealed that
16. Estimation of the parameters of regional frequency distribution
Four first orders Lmoments for each site inside a homogeneous region is making dimensionless by dividing them by the average of the data. Weighted values of dimensionless Lmoments are used to calculate standardized regional Lmoments.
where
Where
17. Spatial interpolation
In classic statistics, the samples derived from a population are usually considered as random sets and the recorded value of a particular variable in an individual sample cannot present any details about the value of that variable in another sample with specific distance. In geostatistics, it is possible to link the values of one variable in a population and distance and direction of samples relating to each other. Furthermore, in classic statistics, it is assuming that the variables are randomly changing, while in geostatistics, some parts of the variable is random and some other parts have structure and is a function of distance and direction. Thus, using geostatistics, first the existence or the absence of a spatial structure between data is considered and then, in presence of a spatial structure, the data will be analyzed. It is possible to adjacent data to be spatially dependent together in a certain distance. In such cases, as in the presence of spatial structure, the variations in a certain space have more chances to be effective on near spaces with regard to more far ones, it is clear that variables are maybe more similar in closer samples.
Geostatistics is a field of statistics with the base of “local variables theory”. Any variable distributed in 3d space with spatial dependence is called local variable and can be studied and analyzed in geostatistical studies. Some method of geostatistical studies are Inverse Distance Weighting (IDW), Global Polynomial (GP), Local Polynomial (LP), Radial Basis Functions (RBF), Kriging (Simple, Ordinary, Universal, Disjunctive and CoKriging).
Applying these methods needs to parallel application of spatial and statistical analysis, which is possible only in some environments like ArcMap. To visualize the above process, the steps are summarized as follows in Figure 1.
18. Comparison of different methods
In above sections, all methods are separately presented but their comparison is very important to choose the most appropriate method for analysis. However, a wrong selection in this step will lead to large amount of uncertainty in output results.
IDW  Deterministic  Prediction  Fast  Yes  Little flexibility, few parameter decision  Few decisions  NO assessment of prediction errors, bull'seyes around data location  None 
GP  Deterministic  Prediction  Fast  No  Little flexibility, few parameter decision  Few decisions  NO assessment of prediction errors, may be too smooth, edge points have large influence  None 
LP  Deterministic  Prediction  Fairly fast  No  Some flexibility, more parameter decision  Flexible  NO assessment of prediction errors, may be hard to choose a good local neighborhood  None 
RBF  Deterministic  Prediction  Fairly fast  Yes  Flexibility, more parameter decision  Flexible  NO assessment of prediction errors, may be too automatic  None 
Kriging  Stochastic  Prediction. Standard Error, Probability, Quintile  Fairly fast  Yes without measurement error, No with measurement error  Very flexibility, assess spatial autocorrelation, obtain standard errors, many decisions more parameter decision  Flexible with modeling tools; prediction standard errors  Many decision on transformations trends, models, parameters, and neighborhoods  Stationary, some methods require a normal data distribution 
Cokriging  Stochastic  Prediction. Standard Error, Probability, Quintile  Fairly fast  Yes without measurement error, No with measurement error  Very flexibility, assess spatial autocorrelation, obtain standard errors, very many decisions  Flexible with modeling tools; prediction standard errors  Many decision on transformations trends, models, parameters, and neighborhoods  Stationary, some methods require a normal data distribution 
*by compus.esri.com 
19. A case study
To make the above discussions more clear, a case study including all steps needed for assessment of climate change effects on positioning of wind power plant station is briefly presented. Figure 2 shows these steps applied in the assessment.
For this purpose, a region in Southern Khorasan, Iran is chosen. Five synoptic stations are considered as reference stations as displayed in Figure 3.
Considering high level of sensitivity of the results to different sources of uncertainty in this study, uncertainty analysis was applied with bootstrap method at 95% confidence interval on the results after downscaling steps (Efron, 1993; Khan, 2006; Fakhry, 2012a; Fakhry, 2012b). Figure 3 shows the results of uncertainty analysis of downscaled data.
After preparing average wind speed parameter for each station, frequency analysis is performed and the Weibull distribution is selected using Lmoments method. Finally, the quantile values corresponding to each return period were derived. Then, geostatistics method was applied for local interpolation and final map for the historical period was prepared.
As the average wind speed was the single used parameter, only one map was drawn. In the case of many existed maps, final positioning could be possible by weighting, according to importance of the maps. Figure 4 shows precedence of potential locations to install wind power station in the region according to long duration of historical records. It is clear from the Figure that the regions with the highest and lowest potential of wind power plant construction are respectively located in the south eastern and north western parts.
To investigate the effect of climate change on positioning process, first the output of HADCM3 model under A2 emission scenario is derived from IPCC website and downscaled by using statistical downscaling techniques. In this study, linear regression method is used and by making wind speed time series of under study stations for near future period (20102039), the frequency analysis is performed again and the future map is presented.
As presented in Figure 6, climate change impacts on wind power plant construction in the regions with low priority will be negligible but in high priority ones, the best points are concentrated in the north eastern part.
20. Conclusions
This chapter made an attempt to initiate discussions on the impact of the climate change on construction of wind farms through various cases that have been presented. Climate change impact studies on wind farms should consider the following principle steps.
This study suggests that the main focus should be on using observational data, time series and frequency analysis and utilizing spatial interpolation methods to create the initial maps, which depict the wind velocity, direction and the average power generation could be calculated. The decision making of the wind farm location could start based on these initial maps and power analysis, which takes various factors into account. The following step would involve downscaling process of the GCM model for initially selected sites. Comparing the GCM model output for future conditions to the historical maps drawn based on existing data. This process will be instrumental in helping to create a path to coupe with a changing climate and its impact on the wind farms.
References
 1.
Ackermann Th 2005 Wind power in power system, Royal Institute of Technology, Stockholm, Sweden.  2.
Akday S. A Bagiorgas H. S Mihalakakou G 2010 Use of Twocomponent weibull mixtures in the analysis of wind speed in the eastern, Applied Energy:2560 2573  3.
Aldrich, John ( 1997 R. A. Fisher and the making of maximum likelihood 19121922. Statistical Science,12 3 162 176  4.
Asif M Muneer T 2007 Energy supply, its demand and security issues for developed and emerging economies, Renewable and Sustainable Energy Reviews,11 1388 413  5.
Bogardi I Matyasovszky I 1996 Estimating daily wind speed under climate change, Solar Energy57 3 239 248  6.
Bonhamcarter G. F 1994 Geographic Information Systems for Geoscientists: Modelling with GIS, 1^{st} Ed. Pergamon Press, Oxford, UK.  7.
Breslow P. B Sailor D. J 2002 Vulnerability of wind power resources to climate change in the continental United States, Renewable Energy:27 585 598  8.
Vela´zquez, S. (Carta J. A Ramirez P 2009 A review of wind speed probability distributions used in wind energy analysis Case studies in the Canary Islands. Renewable and Sustainable Energy Reviews.13.933 955  9.
Chow V. T andD. R Maidment L. W Mays 1998 Applied Hydrology, McGrawHill.  10.
Cunnane C 1989 Statistical distributions for flood frequency analysis, world meteorological organization operational hydrology, Report33 WMONO.718, Geneva, Swiczerland.  11.
Deepthi R Deo M. C 2010 Effect of climate change on design wind at the Indian offshore locations, Ocean Engineering37 1061 1069  12.
Diakoulaki D Karangelis F 2007 Multicriteria decision analysis and costbenefit analysis of alternative scenarios for the power generation sector in Greece, Renewable and Sustainable Energy Reviews;11 716 27  13.
Doria R Madramooto C Mehdi B. B 2006 Estimation of Future Crop Water Requirements for 2020 and 2050, Using CROPWAT. IEEE.  14.
Efron B Tibshrani R. J 1993 An Introduction to the Bootstrap, Monographs on Statistics 57 & Applied Probability, Chapman & Hall/CRC.  15.
Pereira, Fernando R. Martins, Marcelo P. Pes, Eliude I. da Cruz Segundo, André de A. Lyra. (Enio B 2012 The impacts of global climate changes on the wind power density in Brazil, Renewable Energy :1 4  16.
Eslamian S. S H Feizi 2007 Maximum Monthly Rainfall Analysis Using LMoments for an Arid Region in Isfahan Province, Iran. J. Appl. Meteorol. Clim.46 494 503  17.
M, J. (Fakhry M Farzaneh M R Eslamian S. S Khordadi 2012a Uncertainty assessment of downscaled rainfall to investigate impact of climate change on the probability of flood, Journal of Flood Engineering,3 1  18.
M. J. (Fakhry M Farzaneh M R Eslamian S. S Khordadi 2012b Confidence interval assessment to estimate dry and wet spells under climate change in Shahrekord Station, Iran, ASCE, Journal of Hydrologic Engineering, (In Press).  19.
Fisher R. A 1912 On an absolute criterion for fitting frequency curves. Messenger of Mathematics 41 155 160. CP1 Z. in Bennet 1971,1 Z.  20.
Fisher P. F 1989 Knowledgebased Approaches to Determining and Correcting Areas of Unreliabilityin Geographic Databases, Goodchild, M., Gopal, S.: The Accuracy of Spatial Databases,45 54  21.
Greenwood J. A Landwehr J. M andMatalas N. C Wallis J. R 1979 Probability Weighted Moments: definition and relation to parameters of several distributions expressible in inverse form, Water Resources Research,15 5 1049 1054  22.
Harmsen E Miller N. L andSchlegel N. J Gonzalez J. E 2009 Seasonal climate change impacts on evapotranspiration, precipitation deficit and crop yield in Puerto Rico. Agricultural Water Management96 1085 1095  23.
Hiremath R. B Shikha S Ravindranath N. H 2007 Decentralized energy planning; modeling and application review, Renewable and Sustainable Energy Reviews;11 729 52  24.
Hosking J. R 1986 The theory of probability weighted moments. Res. Rep. RC12210, IBM Research Division, Yorktown Heights, New York 10598.  25.
and estimation of distributions using linear combinations of order statistics. J. Royal Stat. Soc. BHosking J. R 1990 Lm. o. m. e. n. t. s Analysis 52 105 124  26.
Hosking J. R 1991 Approximations for Use in Constructing LMoment Ratio Diagrams. Res. Rep. RC 16635, IBM Research Division, Yorktown Heights, New York 10598.  27.
andHosking J. R J. R Wallis 1991 Some Statistics Useful in Regional Frequency Analysis. Res. Rep. RC 17096, IBM Research Division, Yorktown Heights, New York 10598.  28.
andHosking J. R J. R Wallis 1993 Some statistics useful in regional frequency analysis. Wat. Resour. Res.29 271 281  29.
andHosking J. R J. R Wallis 1997 Regional Frequency Analysis: An Approach based on Lmoments. Cambridge University Press, Cambridge, UK.  30.
Hoyme H Zielke W 2001 Impact of climate changes on wind behaviour and water levels at the German North Sea coast. Estuarine, Coastal & Shelf Science53 4 451 458  31.
IPCC 2012 The IPCC Special Report on Renewable Energy Sources and Climate Change Mitigation.  32.
IPCC 2007 in: Alfsen, K., Barrow, Bass, E.B, Dai, Desanker, X.P., Gaffin, Giorgi, S.R.F., Hulme, M.M., Lal, L.J., Mata, L.O., Mearns, J.F.B., Mitchell, T., Morita, R., Moss, D., Murdiyarso, J.D., PabonCaicedo, J., Palutikof, M.L., Parry, C., Rosenzweig, B., Seguin, R.J., Scholes, D., and Whetton, P.H., General guidlines on the use of scenario data for climate impact and adaptation assessment. Cambridge University Press, UK.  33.
IPCC 2001 in: Watson, R.T., Zinyowera, M.C., Moss, R.H., Dokken, D.J. (Eds.)., Sepecial Report on The Regional Impacts of Climate Change, An Assessment of Vulnerability. Cambridge University Press, UK.  34.
cart, Penelope Ramirez, Celia Bueno. (Jose A 2008 A joint probability density function of wind speed and direction for wind energy analysis, Energy Conversion and Managemant,1309 1320  35.
Karvanen J 2006 Estimation of quantile mixtures via Lmoments and trimmed Lmoments, Computational Statistics & Data Analysis51 2 947 959  36.
Khan M. S Coulibaly P Dibike Y 2006 uncertainity analysis of statistical downscaling methods,. journal of Hydrology,319 357 382  37.
Leander R Bouishand T. A 2006 Resampling of regional climate model output for the simulation of extreme river flows. Journal of Hydrology,332 487 496  38.
Pearson K 1902 On the systematic fitting of curves to observations and measurements, Biometrika,1 265 303  39.
S. (Peltola H Ikonen V P Gregow H Strandman H Kilpelainena A Venalainenb A Kellomaki 2010 Impacts of climate change on timber production and regional risks of windinduced damage to forests in Finland, Forest Ecology and Management260 833 845  40.
Pereira Enio B Martins Fernando R., Marcelo P. Pes, Eliude I. da Cruz Segundo, André de A. Lyra. (2012 The impacts of global climate changes on the wind power density in Brazil, Renewable Energy :1 4  41.
Pryor S. C R. J Barthelmie 2010 Climate change impacts on wind energy: A review, Renewable and Sustainable Energy Reviews14 430 437  42.
Sailor D J Smith M Hart M 2008 Climate change implications for wind power resources in the Northwest United States, Renewable Energy33 2393 2406  43.
Sailor D. J Hu T Li X Rosen J. N 2000 A neural network approach to local downscaling of GCM output for assessing wind power implications of climate change, Renewable Energy19 359 378  44.
Sharifi Forood. Bihamta, M.R )Samadi S. Z Mahdavi M 2009 Methodology for selecting the best predictor for climate change impact assessment in Karkheh Basin, Iran, Journal of Environment Engineering and Science, 51(4), 249256.  45.
Valle Costa C La Roverea E., Dirk A. (2008 Technological innovation policies to promote renewable energies: Lessons from the European experience for the Brazilian Case, Renewable and Sustainable Energy Reviews;12 65 90