Exploring Wind Energy Potential as a Driver of Sustainable Development in the Southern Coasts of Iran: The Importance of Wind Speed Statistical Distribution Model

: Wind energy as a clean and inexhaustible source of renewable energy can be a key element of sustainable development that decreases dependence of countries on fossil fuels. Therefore, implementing accurate and comprehensive feasibility studies in countries with a high level of consumption of traditional energy resources is vital; an approach encouraged and supported by green funds and climate change action. It is also crucial to helping spur economic and sustainable growth of these countries. In this regard, this study aims at accurate evaluation of onshore wind energy potential in seven coastal cities in the south of Iran. Six Probability Distribution Functions (PDFs) were examined over representative stations. It was deduced that the Weibull function, which is the most used PDF in similar studies, was only applicable to one station. Here, Gamma distribution offered the best ﬁt for three stations and for the other ones, Generalized Extreme Value (GEV) performed better. Considering the ranking of six examined PDFs and the simplicity of Gamma, it was identiﬁed as the effective function in the southern coasts of Iran bearing in mind the geographic distribution of stations. Moreover, six wind energy converter power curve functions contributed to investigating the capacity factor. It is found that, using only one function could cause under- or over-estimation. Then, stations were classiﬁed based on the National Renewable Energy Laboratory system. Last but not least, examining a range of wind energy converters enabled scholars to extend this study into practice and prioritize the development of stations considering budget limits.


Introduction
Renewable energies are harnessed in various types such as wind power, solar power, biopower, geothermal power, and ocean power. All these types, except geothermal and ocean energy, originate from the infinite energy of the sun, which emits the power about 1.74 × 10 17 w [1]. Considering 1-2 percent of this energy is transformed into wind energy, it is an interminable, environmentally friendly, clean, and reliable source, which is 50-100 times higher than energy conversion from all the plants in the earth combined [2]. It is estimated that the global wind energy potential is about 10 million MW, which could fulfill 35% of the total demand for world energy [3]. The cumulative global installed capacity of wind energy was 539,581 MW in 2017, and it is forecasted that the global wind capacity will reach 800,000 MW by 2021 [4]. In another study, it was estimated that the wind share  [15] and updated by authors based on [16]).

Review of the Literature
Wind power is mainly dependent on wind speed. Therefore, a reliable and accurate model for wind-speed data is a solid pavement for investigating wind energy potential. Numerous models have been used in the scientific literature (Section 3.1) and among them, Weibull is the most widely used distribution. However, all regions and their wind patterns cannot be modeled accurately by Weibull. There are numerous alternative distributions namely Gamma, Lognormal, Rayleigh, Nakagami, Gumbel, Burr, Generalized Extreme Value (GEV), Inverse Gaussian (IG), etc. [17,18]. Hence, in this study, different distributions (Weibull, Gamma, Rayleigh, GEV, IG, and Lognormal) are employed in order to find the most suitable model for wind-speed data based on the goodness of fit criteria.
There are many studies for Wind Resource Assessment (WRA) in Iran. In the past decade (from 2010 to 2019), 28 types of studies related to wind energy assessment in Iran gained more attention and citation. Table 1 presents the distribution function, method of estimation of parameters and the location in each study. Weibull was the most commonly used statistical distribution for describing wind speed data in such studies. In most of them (about 80%), the authors did not consider any further statistical investigations for selecting a suitable wind speed distribution and just used Weibull based on frequent usage while enjoying its simplicity. Moreover, there are a few studies of the coastal regions of Iran and most of them are for inland stations. Thus, a comparative study on the effectiveness of different wind speed distributions in the coastal stations is of great importance. In some of these studies performed statistical investigations to select the distribution or the method of estimation of parameters. For example, Mohammadi and Mostafaeipour used Weibull distribution for their study in Zarrineh city and calculated its parameters using two methods; standard deviation method and power density method [19]. They reported that the power density method is a better estimation method than the standard deviation method in monthly, seasonal, and yearly wind patterns.
In another study, Nedaei et al. compared the performance of four distribution functions (Weibull, Lognormal, Rayleigh, and Logistic) in Mahshahr station [20]. They found that the Weibull distribution is the best function to model wind data in Mahshahr, at two heights of 40 m and 10 m. They also used three different methods (graphical method, maximum likelihood, method of moments) for calculating Weibull parameters. The results show that the Graphical method at 10 m and method of moments at 40 m height are the best methods for calculating Weibull parameters in this city.
Furthermore, Alavi et al. investigated four different distribution functions (gamma, lognormal, Rayleigh, and Weibull) for five cities in Kerman province (Bam, Bardsir, Arzuiyeh, Rafsanjan, and Shahrbabak) [21]. They also used two methods for calculating Weibull parameters; maximum likelihood and method of moments. The experimental results show that the lognormal function produces better estimations for the actual data, while the Weibull model gives a better fit for the truncated wind speed data. They also concluded that calculating Weibull parameters using MLE perform perfectly well in comparison with MM estimations. Alavi et al. used eight PDFs (exponential, Weibull, gamma, lognormal, log-logistic, inverse-Gaussian, generalized extreme value, Nakagami) in their study and also computed the parameters using MLE for wind speed distribution in cities of Chabahar, Khaf, Lutak, Rafsanjan, and Zabol [22]. They found that Nakagami performs better than other distributions; however, because of the closed performance of the Nakagami and Weibull distributions, they finally proposed Weibull due to its flexibility and widely used. In [23], the graphical method, Maximum likelihood, and Method of moments were employed to estimate Weibull coefficients in a coastal area in the south of Iran in the Gulf of Oman. The authors asserted that the maximum likelihood method is the best algorithm to estimate Weibull parameters. Although the resulted output for the designated site is not considerable, more studies should be conducted for onshore wind power assessment in Iranian waters, the Persian Gulf and the Caspian Sea.
An alternative approach was proposed by Faghani et al. [24] used Weibull because of its popularity for WRA in 35 wind stations in nine provinces. Nevertheless, in order to calculate the Weibull parameters three different methods were used; including standard deviation method, empirical method of Lysen, power density method. Based on their investigations, authors concluded that the power density method is the most accurate method for extrapolating wind characteristics.
In a more comprehensive study, Nedaei et al. [25] fitted 46 different PDFs and suggested that the Wakeby distribution performs better than other distribution functions. They did not mention the method of estimating parameters. The distribution analysis is provided for only the ten most efficient PDFs, but details about the methodology for selecting these functions among others are not presented in the context of the paper.

Area of Interest
For the purpose of the current study, seven coastal stations are selected and investigated ( Figure 2). The reason for selecting these stations is that in the first place, the Renewable Energy and Energy Efficiency Organization of Iran (SATBA) has provided public access to validated wind speed and direction data for several stations in Iran [46]. Some of these stations are nearshore and some are not. Since the case study has been on the southern coasts of Iran, which comprises four provinces: Sistan and Baluchestan, Hormozgan, Bushehr, and Khuzestan, in each of these provinces, the nearshore stations with available validated information are selected to cover the whole southern shoreline of Iran from east to west. Hence, these stations have been selected. Table 2 shows some descriptive statistics of wind speed data in each of the seven stations including Latitude (Lat.), Longitude (Long.), data period, time interval, total recorded data, mean wind speed, standard deviation (Sd), maximum wind speed (Max.), locations and the data specifications for each station. According to Table 2, Chabahar (S7 station) has the highest mean wind speed with a value of 4.97 m/s, and therefore it has a great potential for harnessing wind energy.  One of the most effective factors in wind energy studies is wind direction distribution. Wind direction determines the optimal position of wind turbines and the whole array of wind farms. A wind rose is the best indicator for wind direction distribution. This graphical representative plays a vital role in choosing turbines and their alignments in the wind farm. Figure 3 shows wind roses of seven selected stations. As seen in these figures, dominant wind directions for stations S1 to S7 are NW, NW, N, NNE, W, W and SSE, respectively.

Analysis
A plausible resource assessment must be able to identify appropriate locations with strong, moderate, and weak winds. In the past, the assessment was performed through "wind deformed conifer trees". Here, trees and vegetation were the natural indicators for both wind speed and direction. By observing the local place for a long time, a rough idea about the wind richness of the place could be achieved [47]. Various techniques are used to show the effect of wind on trees in terms of numbers such as the Griggs-Putnam index [48]. A comprehensive review of these techniques was performed by Murthy et al. [49].
Although natural indicators are inexpensive, ubiquitous, simple, and conceivable to identify eligible locations, more detailed measurements must be performed to achieve better insight into station characteristics. With the advent of new technologies to measure wind properties, new horizons were broadened for the energy industry. To establish a meaningful analysis, a statistical analysis must be conducted. This procedure demands a time series record of wind speeds. Originally, measurement devices were used for airports across the world. Moreover, there are other technologies for on-station measurements. Therefore, the wind characteristics data pave the way for more precious and detailed researches in the energy industry.
Wind energy is the kinetic energy of air in motion. Total wind energy through an assumed surface area could be calculated through Equation (1): ρ is the air density, v is wind speed and Avt combination is equal to the volume of air passing through the surface at the given time. Air density could be analyzed through statistical distributions to gain more concrete results [50]. In this study, air density is assumed to be constant. Therefore, the term Avtρ is equal to the mass of air flowing through assumed surface area (e.g., Rotor swept area of wind turbine). Power is energy per time unit; therefore, the output power of wind is [51]: Figure 4 shows the flowchart of this study. The initial step for any wind power feasibility study is the technical assessment. The vital part of the technical assessment is to study and analyze the statistical characteristics of wind speeds because the energy output of the field considerably depends on wind speeds and their patterns. Various studies asserted that the main part of the assessment is finding the PDF through wind speeds data as the main input into wind turbine design, station planning and operational plan [52]. After finding a suitable distribution, wind energy potential can be determined. The next step is to calculate the wind turbine capacity factor for all available options. If the capacity factor is less than 0.25, the selected turbine model is rejected. Else, it is selected as a potential option and annual energy output is calculated for it [53].

Wind Speed Distribution Models
In statistics, a unimodal distribution function refers to a continuous probability distribution that contains a single clear peak or the most frequent value. A diverse range of unimodal distribution functions was used in previous studies for modeling the PDF of wind speed including Weibull function, Gamma function, Rayleigh function, logistic function, Log-normal function, etc. Equations of these PDFs are listed in Table 3. Sometimes, wind speed data do not have a single clear peak. In these cases, multimodal probability distributions may provide a better fit for the data. Among these functions, the Weibull function is most used because of its simplicity and flexible parameters [28,40,54]. Although the Weibull is capable of analyzing very low wind speeds close to 0 m/s, it cannot represent the wind structure nature thoroughly [55]. One of the main limitations of the Weibull function is that it cannot accurately capture the effect of calm winds. A solution for this issue is to use some mixture functions like two-component mixture Weibull function (Weibull-Weibull) or truncated normal Weibull mixture (Normal-Weibull). This approach has been comprehensively studied by Tian et al. [56]. Table 3. Investigated PDFs and their parameters.

Name Probability Distribution Functions Parameters
Weibull [57] In more complex and accurate studies, air density could be related to the height from the ground or sea surface and model the relation between them for calculating Equation (2). Therefore, there must be two variations for estimating power output, wind speed, and air density. A Bivariate Probability Model (BPM) was proposed by Carta et al. to model wind speed and air density interactions in Spain [59].
As stated before, Weibull distribution is the most used PDF for the sake of statistical analysis of wind speed data. The first and prior stage to use Weibull distribution for investigating wind speed patterns, is to compute its distribution parameters. There are several methods in scientific contexts about this subject. Some of them are: Three of the above-mentioned methods were most used in various studies: maximum likelihood [67][68][69][70], the empirical method [26,53,71,72], and the graphical method [73,74]. Moreover, the method of moments [64,75,76], energy pattern factor, equivalent energy method, and WAsP [61] are the other methods used for the calculation of Weibull coefficients. Different studies have compared the effectiveness of these methods. As mentioned in Shoaib et al.'s study, using the energy pattern factor method, maximum likelihood, and modified maximum likelihood results stated that maximum likelihood method is the most reliable method [77]. In one of the most comprehensive studies among different algorithms for determining Weibull coefficients, Rocha et al. attempted to evaluate and compare graphical, empirical, moment, Energy pattern factor, maximum likelihood, modified maximum likelihood, and equivalent energy methods in Brazil. It was observed that the equivalent energy method is the best method to calculate coefficients and graphical and energy pattern factor methods are the least effective ones to find coefficient values [63]. Solyali et al., used WAsP (Wind Atlas Analysis and Application), maximum likelihood, and graphical method algorithms [78]. Results indicate that the WAsP algorithm gains the highest correlation with the actual data. Allouhi et al. chose maximum likelihood, graphical, and WAsP methods to compute coefficients in Morocco [62]. Maximum likelihood presents the best fit with the actual data. Masseran, comprehensively compared Weibull, Rayleigh, Lognormal, Burr, Exponential, Inverse Gaussian, and Inverse Gamma in two stations in Malaysia [79]. Surprisingly, the Gamma function performed the best accuracy for the two stations, based on a combination of four goodness of fit indicators. It is worth mentioning here that there are uncertainties and variabilities related to renewable energy sources, e.g., wind. Therefore, it is important to address this challenge in feasibility studies. Using stochastic optimization approaches may be very helpful to address this challenge. For further information, references are recommended [80][81][82][83].

Wind Speed Extrapolation
The data record used in this study represents the wind specifications at the height of 10m above ground. To find wind speed at the height of the turbine hub, power law has been used to convert wind speed. The wind speed fluctuations near the ground surface are referred to as wind shear. In the atmospheric surface boundary layer in which height does not exceed 150 m above the ground surface, the power-law rule is a reliable tool to extrapolate wind speeds at different heights [84]. It can be expressed by Equation (3): where V H is the wind speed at height H, and V R is the reference wind speed at the reference height R which is 10m in the wind speed data, α coefficient is the Hellmann exponent, also named the wind shear coefficient (WSC). The WSC varies regarding ground level height, time of day, atmospheric stability, humidity, and roughness of terrain. The general values of WSCs are summarized in Table 4 [85]. In this study, the WSC coefficient is assumed to be equal to 0.2 due to a lack of reliable data for the stations.

Goodness of Fit Tests
Now it is time to assess the goodness of fitted PDFs for modeling the wind speed probability distribution. There are various Goodness-of-Fit (GoF) tests to assess the quality of fitted distributions. The coefficient of determination (R 2 ) is used to measure the linear relationship between the observed and predicted probabilities. Additionally, root-meansquare error (RMSE) is used to show the level of concentration of data around the fitted distribution. Moreover, because of using the MLE method for parameter estimation Akaike information criterion (AIC) and Bayesian information criterion (BIC) is used to assess the accuracy of the fitted distribution. Table 5 presents the formulae and definitions of parameters for each of these four statistical indicators. Lower values for RMSE, AIC, and BIC indicate higher goodness of fit, while on the contrary, a larger value for R 2 shows better effectiveness of the fitted distribution. Table 5. Goodness-of-Fit (GoF) tests in this study.

Indicator
Formula Parameters y i : Observed data y ic : fitted data n : Number of data samples.

Wind Power and Energy Density
After calculating Weibull distribution function parameters, the next crucial step will determine wind power density through using Weibull distribution parameters. Wind power density is one of the most meaningful indicators to depict how powerful the winds are in a region within the defined time period. By contemplating the wind power density, quantification of the potential wind energy electricity is possible. Wind power is directly related to the area of the turbine (i.e., swept area) and the cube of wind velocity at the station. Wind power can be calculated as Equation (4): Note that in the above equation, the area of turbine and air density is considered as constant values. The more tangible indicator for wind energy is the power per unit area. It can be written as Equation (5): The wind power density using Weibull distribution for wind speed patterns will be stated as [53,62]: While using the Rayleigh probability distribution function as a distribution of wind speeds, the wind power density function will be written as [86]: For another distributions, there are not any analytic answers for the integral in Equation (4). Therefore, the wind power density is calculated using numerical integration.
There is a reliable classification for every station based on National Renewable Energy Laboratory (NREL) to determine the relative wind power potential. This classification is based on observed and recorded wind speeds at the height of 50 m. The wind power classification is given in Table 6 [87]. Wind energy density is the amount of energy produced for any arbitrary period of time. It can be written as Equation (8):

Capacity Factor
The capacity factor (C f ) of a wind turbine is an indicator that defines the output viability of a wind turbine at a selected station. It determines the ratio of average power yield to the rated power of the turbine. C f is one the most reliable measures for choosing wind turbine because it inherently shows the performance of the wind turbine. C f can be expressed as [53]: Installed power (MW) . Annual hours (9) P r Represents the power output of wind turbine when the turbine runs as its rated power all the time in one year. Needless to say, the capacity factor is always less than 1. The P avg can be expressed as [62]: where f (v) is the probability distribution of wind speed and P(v) is the power curve of wind turbines. The cut-in wind speed is the minimum wind speed at which the turbine blades overcome the friction and begin to rotate. The cut-out wind speed is the speed at which turbine blades stopped rotating to prevent probable damages from high winds. It is notable that not all turbines have a well-defined cut-out speed. The schematic concept of power curve for a wind turbine is presented in Figure 5. For pitch controlled turbines, the generated power from wind turbines can be estimated from the power curve as [53]: The function c(v) denotes the wind turbine's output power in the interval between cut-in and cut-off speeds. In recent years, several studies have been conducted by scholars aiming to approximate this function [53,89,90]. The six most widely used functions for this purpose are listed as follows from Equations (12)- (17): The simultaneous power output of wind turbines in a wind farm is not the same because of the spatial distribution of turbines and the unpredictable and stochastic nature of wind speed distribution. Therefore, a more complicated approach is required to obtain the power output. Based on Wang et al., the average of the above six functions is used for estimating the capacity factor [52].
Therefore, the wind turbine will work within the speed period between cut-in and cut-off speeds.

Availability Factor
According to Section 2, wind turbine will work within the wind speed interval between the cut-in and cut-off speeds. Thus, the availability factor of a wind turbine can be expressed as Equation (19): (19) where f (v) is the wind speed probability distribution function.

Analysis of Distribution Functions
Six different probability distributions, which are most widely used in the scientific literature, are fitted to the wind speed data of each station. The Rayleigh PDF as a oneparameter function, Gamma, Log-normal, Weibull and Inverse Gaussian as two-parameter functions and finally, GEV as a three-parameter function are used for this purpose. MLE is used to estimate PDF parameters, since it is more efficient than other estimation methods, and generates a lower mean squared error [22]. The estimated parameters for each region are reported in Table 7. Note that the parameters are calculated from wind data at 10 m above the ground. The wind speed distribution diagrams (PDF plots) for seven selected stations are shown in Figure 6. The horizontal axis of these diagrams shows the range of wind speed and the vertical axis is the probability density. It varies from zero to the highest possible value for each curve. As seen in figures, all distributions are skewed to the right (i.e., positive Skewness), and therefore, the mean of wind speeds is bigger than the mode and median of data.
Four different GoF indicators are calculated for these functions and the results are reported in Table 8. Note that lower RMSE, AIC, and BIC show better fitness. On the other hand, higher R 2 shows a better correlation between observed data and fitted distribution. Note that different GoF indicators can yield different results. For example, in the S3 station, Gamma performs better in terms of R 2 , whereas Lognormal performs better in terms of RMSE. This paper assigns R 2 a greater weight for the assessment and selects it as the first reference index. Results show that Gamma is the best distribution for S1, S3, S5, and GEV has the best fit for S2, S4, and S6. Weibull is only suitable for S7. Results show that R 2 values for Weibull distribution are 1 to 7 % lower than that of the best distribution in stations S1 to S6.
As seen above, one cannot use Weibull distribution for all stations without further investigations. Nevertheless, for engineering applications, it is valuable to propose an effective function that has the best performance to be used for a specific geographical location. The results show that Gamma and GEV are the most effective function for the southern coasts of Iran since they have the best statistical performance. However, Gamma has two parameters, whereas GEV has three parameters. Therefore, Gamma is selected as the effective function for the southern coasts because of its more simplicity. It should be mentioned that by increasing the length of the data and the number of distributions and even the fitting method, one can clearly expect different fitted distributions. This originates from the nature of statistical research. Considerably more work will need to generalize the obtained distributions and to propose them as the optimal distributions for the selected sites.

Analysis of Wind Power and Energy Density
The wind power density for each station is determined using the best-fitted probability distribution at 10, 30, and 50 m above the ground. Once the wind power density is calculated, wind energy density can be found. In this study, wind energy density is computed for one year. The results are gathered in Table 9. For the sake of comparison, wind powers for each height across the stations are shown in Figure 7. Table 9. Estimated wind power density and wind energy density and classification of the stations for the wind power based on NREL.

Wind Turbine Selection
Fifteen wind turbines with hub heights ranging 20-100 m and rated powers ranging 0.055-4.2 MW are exploited for assessing wind performance. This range of wind turbine power classes could be useful in the next steps of economic and feasibility analysis for establishing wind farms. Wind turbines characteristics are listed in Table 10.
To evaluate the energy output of each turbine, wind speed distribution parameters have been calculated at the hub-height. Consequently, wind turbine capacity factor and availability factor are determined using Equations (18) and (19), respectively. For cost-effective investment in the wind energy sector, it is asserted that the capacity factor of wind turbines should be between 0.25-0.45 to be economically efficient and feasible [92]. According to Ayodele et al., any wind turbine with a capacity factor lower than 0.25 will not be suitable to be integrated into the grid [53]. Estimated capacity factors for all wind turbines are listed in Table 11. Capacity factors higher than 0.25 could be plausible for a wind turbine [53]. Based on this, in Abadan, T11 to T15 are suitable for wind energy production. The table reveals that in Mahshahr again T11 to T15 are reliable for energy production. After that, in Delvar station, there is not any capacity factor higher than 0.25 but, T13 has the highest capacity factor. In Bordekhoon, one of the most suitable stations, T3 and T11 to T15 are suitable choices. Kish, the most important island of the Persian Gulf for tourism, again T11 to T15 are good for wind energy production. Jask has not any acceptable capacity factors due to poor wind power density. Eventually, in Chabahar station, T2, T3, and T11 to T15 are acceptable wind turbines with significant capacity factors.
The availability factor for each wind turbine in each station has been calculated and shown in Table 12. As shown in Table 12, Bordekhoon and Chabahar stations have the greatest availability factors among stations. This fact could be realized from their greater wind power density and their higher average wind speed. After them, Kish, Mahshahr, and Abadan have the greatest availability factors. Jask and Delvar are the least in availability factors due to their poor wind power density. Among wind turbines, T5 has the greatest average of availability factor because it has the least cut-in speed and this turbine cover the wider range of applicable wind speeds.

Comparison with Previous Studies
In Abadan (S1), there is one similar study for wind energy potential [30]. The data used in this study is the same as in the current research. Weibull distribution was deployed for characterizing wind speed patterns. Nevertheless, the Weibull distribution ranked fourth among distributions used in this study. The best PDF for S1 is Gamma distribution. A precise comparison for the goodness of fit is not possible because it was not reported in that paper.
For Mahshahr city (S2), Nedaei et al. studied data with 19 months duration for the region and used Weibull, Rayleigh, and Lognormal distributions [20]. Reported R 2 in that paper are 0.923, 0.910, and 0.901, respectively, which are consistent with R 2 of the current study with the values 0.926, 0.924, and 0.917. Furthermore, in the current study, GEV has the best performance in terms of R 2 value and has R 2 = 0.988.
In Delvar city (S3), Dabbaghiyan et al. studied wind speed data for 2011 [40]. They deployed Weibull distribution for wind speed, but R 2 is not published. While, in this study, Gamma function is selected for this location with R 2 = 0.969.
Two studies were conducted on Kish Island (S5). Nedaei studied wind data between 2006 and 2007 and also selected Weibull distribution [31]. R 2 for Weibull distribution was 0.986, while Gamma distribution has the highest R 2 equal to 0.993 in the current study. Although R 2 is slightly increased, different data sources should be considered. Another study in the region is performed by Mohammadi et al. [35] based on long-term data from 2002 to 2009. Again, Weibull distribution opted. R 2 was not calculated. Wind power density at the height of 10m reported 111 W/m 2 , while this research calculated 144 W/m 2 using Gamma density function with R 2 equal to 0.993.
In the Jask station (S6), Nedaei studied wind power density using data from 2006 to 2007 and Weibull distribution [31]. R 2 reported 0.9782 while in current research GEV function was the most suitable for the region with R 2 = 0.986. Weibull distribution ranked third among six PDFs with R 2 = 0.975, which is in line with the previous study.
Chabahar (S7) is the most studied location with three studies. Minaeian et al. studied wind speed data for 2007 and used Weibull distribution but R 2 was not reported [41]. Mohammadi et al. studied data from 2002 to 2009 and used Weibull distribution [35]. Again, R 2 was not announced. Wind power density at the height of 10 m is calculated 111 W/m 2 while in the current study, using Gamma function, wind power density is determined as 119 W/m 2 which is slightly higher. This distinction might cause due to different wind speed data and distribution used. Moreover, Alavi et al. studied the station with data from 2008 to 2009. They conducted analysis using Weibull, Gamma, Lognormal, and GEV functions with R 2 = 0.999, 0.999, 0.998, and 0.999 respectively which are significantly high. Additionally, Nakagami distribution function yielded the best fitness with R 2 = 0.9999. Accordingly, in the current study, R 2 for those distributions are 0.994, 0.984, 0.923, and 0.969, respectively. Apart from Nakagami distribution, Weibull shows the best fitness in both studies and R 2 of the two analyses are approximately equal. The negligible difference might occur because of different data. It should be noted that this region is very important in the development programs of Iran and therefore has great potential for the construction of coastal and offshore structures [93,94].
As mentioned above, current research enhanced goodness of fit for wind speed data according to R 2 in comparison with previous studies. Moreover, in other sections of the analysis, a more precise approach is conducted to compute capacity factors. Additionally, to increase the practicality of the article, a broad range of wind turbines are considered to analysis to obtain a more concrete insight toward wind energy capacity in the south coastal zone of Iran.

Conclusions
This paper presented a comprehensive and step-by-step methodology for the wind energy assessment along with a case study in seven different stations covering the southern coastal borders of Iran. Moreover, an extensive literature review is conducted to gain insight into the concept of wind energy assessment and efforts made in this field, their methodologies, and locations. The majority of previous studies in Iran used the Weibull function in order to estimate wind speed patterns in the selected area. Here, six different probability distribution functions ranging from one to three parameters were exploited in order to find the best fitness to the wind speed data. Furthermore, they exploited only a single function for determining the capacity factor of wind turbines. This study has shown that although the Weibull PDF is the most well-known function as wind speed distribution, in the area of interest, Gamma, GEV, and lognormal functions have better correlations with observed data. The second major finding was that based on the distribution ranks, the Gamma function is considered as the most suitable distribution at the initial step of wind power assessment in the southern coasts of Iran. Additionally, the annual wind power density and energy output are determined for three heights above the ground surface. Based on power densities, Bordekhoon, Kish, and Chabahar are more suitable for harnessing wind energy. Jask has the least power density and Abadan, Mahshahr, and Delvar are suitable for stand-alone wind power facilities. Moreover, the capacity factor of wind turbines is calculated through six different functions and the average of these six functions is used as a final capacity factor value to have more reliable results. This provides an accurate evaluation of wind energy capacity. To analyze practical aspects of the current study, fifteen wind turbines with hub heights ranging 20-100 m and rated powers ranging 0.055-4.2 MW are exploited for assessing wind performance. The achieved results show that as capacity factor increases with the hub height of wind turbine, T11 to T15 are the best wind turbines due to their capacity factors. The present study is the primary comprehensive practical assessment of wind energy along with assessing different actual turbine models.
Further research can be undertaken to explore economic aspects of implementing wind energy farms in the southern coasts of Iran by employing the levelized cost of energy (LCOE). This approach is used in similar research by the authors [10] for the tidal current energy in this area.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: