An Assessment of Onshore and Offshore Wind Energy Potential in India Using Moth Flame Optimization

Wind energy is one of the supremely renewable energy sources and has been widely established worldwide. Due to strong seasonal variations in the wind resource, accurate predictions of wind resource assessment and appropriate wind speed distribution models (for any location) are the significant facets for planning and commissioning wind farms. In this work, the wind characteristics and wind potential assessment of onshore, offshore, and nearshore locations of India—particularly Kayathar in Tamilnadu, the Gulf of Khambhat, and Jafrabad in Gujarat—are statistically analyzed with wind distribution methods. Further, the resource assessments are carried out using Weibull, Rayleigh, gamma, Nakagami, generalized extreme value (GEV), lognormal, inverse Gaussian, Rician, Birnbaum–Sandras, and Bimodal–Weibull distribution methods. Additionally, the advent of artificial intelligence and soft computing techniques with the moth flame optimization (MFO) method leads to superior results in solving complex problems and parameter estimations. The data analytics are carried out in the MATLAB platform, with in-house coding developed for MFO parameters estimated through optimization and other wind distribution parameters using the maximum likelihood method. The observed outcomes show that the MFO method performed well on parameter estimation. Correspondingly, wind power generation was shown to peak at the South West Monsoon periods from June to September, with mean wind speeds ranging from 9 to 12 m/s. Furthermore, the wind speed distribution method of mixed Weibull, Nakagami, and Rician methods performed well in calculating potential assessments for the targeted locations. Likewise, the Gulf of Khambhat (offshore) area has steady wind speeds ranging from 7 to 10 m/s with less turbulence intensity and the highest wind power density of 431 watts/m2. The proposed optimization method proves its potential for accurate assessment of Indian wind conditions in selected locations.


Introduction
As populations and technological developments grow, the need for electricity consumption also increases, which impacts the environmental pollution conditions based on the type of fuels used for power generation. Remarkably, the fossil fuel-based generation from coal and oil threatens the climatic conditions due to its pollution particles (suspended in the air) which increase carbon emission greatly. As per the Paris Agreement commitments, the rising magnitude of average global temperature should be maintained at about 1.5 • C. Therefore, the Indian government has set an ambitious target to achieve • Average, standard deviation, minimum and maximum wind speed (m/s) are measured for wind generation suitability assessment, and turbine selection. • Average, standard deviation and maximum gust direction (degrees) are estimated for optimizing wind turbines and understand the spatial distribution of wind • Temperature ( • C) and vertical wind speed (m/s) are used for the air density and turbulence application respectively to measure the average and minimum/maximum value. • Average, standard deviation and minimum/maximum value of barometric pressure (kPa) are measured for air density applications • Relative humidity (%) and solar radiation (W/m 2 ) (average, minimum/maximum value) are measured for the icing and atmospheric analysis respectively.

Wind Speed Distribution Methods
The wind speed and its frequency of occurrence to form a wind distribution pattern is a critical piece of information. They are used to evaluate the power output of the wind turbine directly. The frequency of wind distribution denotes the number of intervals during period of data collected when the perceived wind speed falls within specific bins. They are normally occurred about 0.5 or 1 m/s space. They cover a minimum range of speeds described for the turbine power curve, i.e., 0 to 25 m/s and even more. The usual reports are presented in the bar chart/histogram and input as a table for software analysis. The probability density function is a mathematical one that visualizes the probable pattern of wind speed, which is random by nature in a continuous period. The detailed wind distribution models for the calculation of probability density function and cumulative distribution function [14] are presented in Appendix B.

Goodness of Fit
The goodness of fit is a statistical analysis method for qualifying the ten wind probability distributions considered for evaluation to match with actual data measured. The best scored root mean square error (RMSE) and correlation coefficient (R 2 ) [11] finds the best fit wind distribution method for the defined location. They are used to determine the fitness of ten distributions. These can be calculated as follows: 2.3.1. Root Mean Square Value (RMSE) RMSE estimate the residuals of the frequency of examined PDF (Probability density function) and the measured data.
where y i is the observed data plotted in a histogram, n stands for the number of bins and x i is an estimated PDF of wind distribution considered for evaluation. The evaluated RMSE value should be close to zero to form a best fit.

The Coefficient of Determination (R 2 )
The R 2 evaluates the linear correlation between predicted values from PDF functions and measured data.
where y defines the mean of the total observed value and the value of R 2 lies between 0 to 1 (closer to 1 is the best fit).

Wind Site and Measurement Details
The wind data collected from two Indian states namely Gujarat and Tamilnadu are considered for analysis. These two states are pioneers in wind energy generation contributing to the Indian power grid in a higher stake. The wind data segregated into three categories such as onshore from Kayathar (Tamilnadu, India), offshore from Gulf of Khambhat (Gujarat, India), and nearshore from Jafrabad (Gujarat, India). The statistics of each location are presented in Table 1.

Kayathar (Tamilnadu)-Onshore Location
The Tamilnadu State, situated at the southeastern end of the Indian peninsula, comprises of a higher rate of humidity and temperature throughout the year. The annual rainfall received from the South-West monsoon (June to September) and the North-East monsoon (October to December) is 911.6 mm and temperature ranges between 20 • C and 38 • C in the plains [28]. Kayathar is a small town located in the southern part of the Tamilnadu. The mapping coordinates of Kayathar are 08 • 57 N 77 • 48 E and it is situated at 78 m elevation (altitude) above sea level. The geolocation of the targeted site is shown in Figure 1a. The wind data collected from the mast installed at a height of 120 m by the National Institute of Wind Energy (NIWE) (the central agency for wind energy developments in India). Wind data measured at Kayathar were recorded with an average of 10 min along with different altitudes, directions of wind flow, temperature, and pressure for three years between the years 2014 and 2016. This region experiences a template of monsoon variations from South-West and North-East monsoons. Therefore, the summer season is very hot with high humidity, and hence this site is prone to high windy during monsoons. The wind data collected from this site is classified into three categories likely seasonal, yearly, and monthly to evaluate the wind characteristics [29].

Wind Site and Measurement Details
The wind data collected from two Indian states namely Gujarat and Tamilnadu are considered for analysis. These two states are pioneers in wind energy generation contributing to the Indian power grid in a higher stake. The wind data segregated into three categories such as onshore from Kayathar (Tamilnadu, India), offshore from Gulf of Khambhat (Gujarat, India), and nearshore from Jafrabad (Gujarat, India). The statistics of each location are presented in Table 1.

Kayathar (Tamilnadu)-Onshore Location
The Tamilnadu State, situated at the southeastern end of the Indian peninsula, comprises of a higher rate of humidity and temperature throughout the year. The annual rainfall received from the South-West monsoon (June to September) and the North-East monsoon (October to December) is 911.6 mm and temperature ranges between 20 °C and 38 °C in the plains [28]. Kayathar is a small town located in the southern part of the Tamilnadu. The mapping coordinates of Kayathar are 08°57′ N 77°48′ E and it is situated at 78 m elevation (altitude) above sea level. The geolocation of the targeted site is shown in Figure 1a. The wind data collected from the mast installed at a height of 120 m by the National Institute of Wind Energy (NIWE) (the central agency for wind energy developments in India). Wind data measured at Kayathar were recorded with an average of 10 min along with different altitudes, directions of wind flow, temperature, and pressure for three years between the years 2014 and 2016. This region experiences a template of monsoon variations from South-West and North-East monsoons. Therefore, the summer season is very hot with high humidity, and hence this site is prone to high windy during monsoons. The wind data collected from this site is classified into three categories likely seasonal, yearly, and monthly to evaluate the wind characteristics [29].

Gulf of Khambhat (Gujarat)-Offshore Location
Gujarat state has an international boundary with common border at the North-Western (NW) fringe [30]. It has a long coastline of about 1600 km, which is the longest among all states of the country . India earmarked the states Gujarat and Tamilnadu for pacing up rapidly towards offshore wind energy development through the supervision of the nodal ministry MNRE-Ministry of New and Renewable Energy and nodal agency NIWE-National Institute of Wind Energy). The offshore development activities are started with preliminary site assessments in Gujarat coastal area. The proposed location of the offshore wind farm is about 23 to 40 km off the coast from the Pipavav port, Gulf of Khambhat, Gujarat. It covers about 400 km 2 and the wind measurements and other data collections are through LiDAR. Figure 2 displays the placement of the LiDAR instruments, showing a monopole structure erected in the Gulf of Khambhat at 20 • 45 19.10" N, 71 • 41 10.93" E with 25 km distance from the coast and 15 m water depth. The nearest port is Pipavav located at approximately 23 km from the Gujarat coast in the South-East direction. The monopole structure with a platform comprises of an anemometer, wind vane, temperature monitor, and pressure instruments are also commissioned. The LiDAR instrument is a wind cube configured for 12 various heights with a minimum height of 40 m and maximum height of 200 m. The overall data collected from LiDAR for the period from December 2018 to November 2019 is about 75.85% of the total, and the remaining 24% of data are not successful due to technical issues.

Gulf of Khambhat (Gujarat)-Offshore Location
Gujarat state has an international boundary with common border at the North-Western (NW) fringe [30]. It has a long coastline of about 1600 km, which is the longest among all states of the country with a geographical area of 1.96 lakh sq. km. It is situated between 20°1′ and 24°7′ North latitudes and 68°4′ and 74°4′ East longitudes with an elevation of about 9 m [30]. India earmarked the states Gujarat and Tamilnadu for pacing up rapidly towards offshore wind energy development through the supervision of the nodal ministry MNRE-Ministry of New and Renewable Energy and nodal agency NIWE-National Institute of Wind Energy). The offshore development activities are started with preliminary site assessments in Gujarat coastal area. The proposed location of the offshore wind farm is about 23 to 40 km off the coast from the Pipavav port, Gulf of Khambhat, Gujarat. It covers about 400 km 2 and the wind measurements and other data collections are through LiDAR. Figure 2 displays the placement of the LiDAR instruments, showing a monopole structure erected in the Gulf of Khambhat at 20°45′19.10″ N, 71°41′10.93″ E with 25 km distance from the coast and 15 m water depth. The nearest port is Pipavav located at approximately 23 km from the Gujarat coast in the South-East direction. The monopole structure with a platform comprises of an anemometer, wind vane, temperature monitor, and pressure instruments are also commissioned. The LiDAR instrument is a wind cube configured for 12 various heights with a minimum height of 40 m and maximum height of 200 m. The overall data collected from LiDAR for the period from December 2018 to November 2019 is about 75.85% of the total, and the remaining 24% of data are not successful due to technical issues.

Jafrabad (Gujarat) Nearshore Location
The meteorological mast tower is installed with 100 m height at the Jafrabad coast in line of sight with the LiDAR location, at the distance of 25 km (Figure 1b). The main purpose of this installation is to correlate and validate with the LiDAR measurement for assessing offshore Gulf of Khambhat wind potential. The wind speed is measured from anemometers and temperature sensors placed at various levels such as 100, 80, 50, and 20 m above ground level. The wind direction sensors are positioned at 98, 78, and 48 m. This mast is installed at 20 53′29.81″ N, 71 27′35.68″ E with an altitude of 9 m above the ground level, with a mast height of 100 m. The wind data was collected with a 99.89% recovery rate for the period from December 2018 to November 2019.

Jafrabad (Gujarat) Nearshore Location
The meteorological mast tower is installed with 100 m height at the Jafrabad coast in line of sight with the LiDAR location, at the distance of 25 km (Figure 1b). The main purpose of this installation is to correlate and validate with the LiDAR measurement for assessing offshore Gulf of Khambhat wind potential. The wind speed is measured from anemometers and temperature sensors placed at various levels such as 100, 80, 50, and 20 m above ground level. The wind direction sensors are positioned at 98, 78, and 48 m. This mast is installed at 20 53 29.81" N, 71 27 35.68" E with an altitude of 9 m above the ground level, with a mast height of 100 m. The wind data was collected with a 99.89% recovery rate for the period from December 2018 to November 2019.

Seasonal Wind Periods
The highest wind power potential in India concentrated on two major states, namely Tamilnadu [28] and Gujarat [29]. Concerning the classification of seasons for Gujarat and Tamilnadu states, the influence of seasonal winds from North-East monsoon (NEM) and South-West monsoon (SWM) are taken into consideration for analysis [30]. These seasons are grouped based on seasonal months, as follows: Winter: January and February (2 months)

Moth Flame Optimization (MFO) Method
The application of artificial intelligence and soft computing techniques along with optimization is adopted in various activities. In this work, this method is implemented for estimating the parameters of the bimodal Weibull-Weibull method. The MFO algorithm [31] is a naturally inspired algorithm where moths are fancy insects, which are highly similar to the family of butterflies. The main inspiration of MFO is the navigation method of moths in nature, called transverse orientation. Moths navigate at night by maintaining a fixed angle to the moon. This is an effective mechanism for traveling long distances in a straight line [32]. It is considered as one of the promising metaheuristic algorithms and successfully applied in various optimization problems in a wide range of fields, such as power and energy systems, economic dispatch, engineering design, image processing, and medical applications. The working principle of MFO is based on the individuals/set of moths (M) called population and flame (F) is said to be the best solution for each moth. The moth-flame combination contains one flame per moth, which is considered to be as its best position. During the iteration, this flame will get updated if any better solution found. The matrix OM stores the corresponding fitness (objective) value, which can be written as follows: where n is the number of moths, I is the objective function, and d is the number of variables. Notably, each moth flies around its corresponding flame. Therefore, the flames matrix is the same size as the moths matrix [30]. The set of flames can be represented in the matrix F, while the matrix OF indicates the corresponding fitness value and can be written as follows.
The main difference between moths and flames is the way they are treated and updated in each iteration. The moths are search agents moving around the search space, while flames are the best position of moths obtained by the current iteration [30]. The equation to find the best optimal value is given as: where M i indicates the i-th moth, Fj indicates the j-th flame, and S is the spiral function. The spiral function of the moth is given as: where D i indicates the distance of the i-th moth for the j-th flame, b is a constant for defining the shape of the logarithmic spiral, and t is a random number. D i is calculated as follows: Energies 2020, 13, 3063 where F i indicates the i-th flame. The moth can converge or exploit in the given search space by changing the value 't' when the position of a moth gets changed. A mechanism to solve the number of flames during each iteration is mathematically defined as follows: where N is the maximum number of flames, l is the number of current iteration and T is the maximum number of iterations. The objective function to be considered for minimization of estimated Bimodal Weibull distribution parameters by MFO and observed actual wind speed distribution is derived as: where f est is the estimated value from bimodal-MFO and f o is the observed data from the histogram: where, the parameters k1 and k2 are shapes, c1 and c2 stand for scales, w1 and w2 define weights and v terms the wind speed. The MFO pseudocode algorithm with program logic is derived in Appendix A and the flow chart for the MFO algorithm is described in Figure 3.
Energies 2020, 13, 3063 8 of 41 where Fi indicates the i-th flame. The moth can converge or exploit in the given search space by changing the value 't' when the position of a moth gets changed. A mechanism to solve the number of flames during each iteration is mathematically defined as follows: where N is the maximum number of flames, l is the number of current iteration and T is the maximum number of iterations. The objective function to be considered for minimization of estimated Bimodal Weibull distribution parameters by MFO and observed actual wind speed distribution is derived as: where is the estimated value from bimodal-MFO and is the observed data from the histogram: where, the parameters k1 and k2 are shapes, c1 and c2 stand for scales, w1 and w2 define weights and v terms the wind speed. The MFO pseudocode algorithm with program logic is derived in Appendix A and the flow chart for the MFO algorithm is described in Figure 3.

Results and Discussions
This segment outlines the statistical investigation of the wind data measured in the three different landscapes. A total of ten wind distribution functions are adopted to estimate the correctness of the distribution that best fits the wind speeds of the targeted location. Later, the optimization method, i.e., the MFO method is applied for parameter estimation and compared with ten wind distribution parameters. The MATLAB R2018b software with in-house developed code is used for estimating the parameters of ten wind distribution functions along with the MLM. Then, it is compared with the estimated parameters through MFO method using bimodal (Weibull-Weibull) distribution. As stated earlier, the goodness of fit is to evaluate the best wind distribution method with a lower rate of RMSE value nearer to zero. Additionally, the higher rate of R 2 should be nearer to one with reduced error fitting for evaluating the wind power density. Further, the results are grouped for the best wind power density with the selected distribution. Later, the annual mean wind speed along with mean turbulence intensity (15 m/s) are considered for categorizing the turbulence class for the turbine selection to form a standalone wind turbine/wind farm. The complete steps/procedures for the assessment of wind resources are given in Figure 4.

Results and Discussions
This segment outlines the statistical investigation of the wind data measured in the three different landscapes. A total of ten wind distribution functions are adopted to estimate the correctness of the distribution that best fits the wind speeds of the targeted location. Later, the optimization method, i.e., the MFO method is applied for parameter estimation and compared with ten wind distribution parameters. The MATLAB R2018b software with in-house developed code is used for estimating the parameters of ten wind distribution functions along with the MLM. Then, it is compared with the estimated parameters through MFO method using bimodal (Weibull-Weibull) distribution. As stated earlier, the goodness of fit is to evaluate the best wind distribution method with a lower rate of RMSE value nearer to zero. Additionally, the higher rate of R 2 should be nearer to one with reduced error fitting for evaluating the wind power density. Further, the results are grouped for the best wind power density with the selected distribution. Later, the annual mean wind speed along with mean turbulence intensity (15 m/s) are considered for categorizing the turbulence class for the turbine selection to form a standalone wind turbine/wind farm. The complete steps/procedures for the assessment of wind resources are given in Figure 4.

Wind Characteristics
The detailed wind statistical analysis for the Kayathar station (onshore), Jafrabad station (Mastnearshore), and Gulf of Khambhat station (offshore) are presented in this section. The key characteristics of any wind assessment are to ascertain the wind resource feasibility such as mean wind

Wind Characteristics
The detailed wind statistical analysis for the Kayathar station (onshore), Jafrabad station (Mast-nearshore), and Gulf of Khambhat station (offshore) are presented in this section. The key characteristics of any wind assessment are to ascertain the wind resource feasibility such as mean wind speed, standard deviations of wind, turbulence intensity, and wind directions. The collected statistical analysis is used to determine the site eligibility for wind turbines, selection and energy yield to make the wind project technically viable. The statistical analysis of wind speed, wind direction, and its relation is represented by the wind rose plot. The mean wind speed over the periods for seasonal, monthly, and annual analysis represent the energy contained in the wind, and for the analysis of turbulence intensity to classify turbines. The PDF for fitting the wind speed distribution, wind power potential analysis, and statistical characteristics are presented in the following sections.   Table 3 illustrates the seasonal wind speed peaks on the South-West monsoon period (SWM) of about 10.03 m/s and maximum seasonal wind speed observed as 18.52 m/s. Comparing the monsoon periods, the North-East monsoon (NEM) fetches a low wind mean speed of 4.5 m/s. The monthly mean wind speed is observed as having a maximum during June of 10.29 m/s and the maximum wind speed during July is 19.81 m/s. The monthly standard deviation is stable from January to March with an average value of 1.8 m/s. Further, it is observed as 3.99, 3.5, and 3.01 m/s for May, July, and October respectively during the peaks in monsoon periods due to seasonal winds. The mean turbulence intensity (MTI) is estimated as a ratio between the standard deviation of the experiential wind speed every 10 min and the mean observed wind speed for the same period. Additionally, the turbulence is higher at a low wind speed as shown in Figure 5c. The defined location, namely the Kayathar wind station, falls on turbulence category 'B' (as per International Electrotechnical Commission (IEC) standard) with MTI at 15 m/s of 0.132, i.e., 13.2% (Table 3).    Figure 5 shows the annual wind statistics at 100 m height, wind speed at different altitude and intensity of turbulence. It is observed that the wind speed increases on higher altitudes. When considering the 10 m, the wind speed attained a scale of 2 m/s and further reaches about 4 m/s (100 m). It is known fact that the wind power is the cubic of wind speed; it is suggested to erect the turbines at Energies 2020, 13, 3063 12 of 41 higher altitude and hub height to obtain higher wind energy extraction. The wind speed observed in the year 2016 was not great, and particularly relatively lower than 2015 and 2014 due to climatic changes. Table 3 shows the observed wind speed of the Kayathar station for seasonal case and it recorded better characteristics. It is situated in the southern part of Tamilnadu and produces more wind power during the SWM, which stakes the major portion of state demand.
Furthermore, the annual maximum wind speeds hinge on wind direction; the West direction observed a maximum wind speed of 10.26 m/s (38.66%) followed by WSW and WNW, i.e., 11.27% and 7.57% respectively, as shown in Table 4. The second maximum wind generation of NEM in NE wind direction is observed as 6.51%, with a mean wind speed of 4.766 m/s followed by NNE (North-Northern east) of 6.50%. The wind rose plots exposed in Figure 6 indicate the influence of South-West and North-East monsoons. The North-East monsoon influenced the wind direction during October to December, and the West direction shifts to N-NNE (Northern NNE). During SWM periods, (i.e., June-September) wind directions are found to be West and WSW (West southern west). It also represents the main wind direction from the West and ranges between 258.75 • -281.25 • during the South-West. The winter season experiences NEM winds and influences at the end of the South-West monsoon period. The summer wind rose plot experiences wind from SWM and the wind direction from SE-SSE (Southeast-South Southeast) and Western direction. It indicates that the Kayathar region experiences maximum wind power from the Southwest direction. During low-wind periods from NEM and winter seasons, the wind speed reduced relatively which deteriorate to low wind generation.

Gulf of Khambhat (Gujarat Offshore) Station
The mean wind speed (MWS) collected from LiDAR measurements at various heights (40 to 200 m) are presented in Figure 7 for the period of December 2018 to November 2019. It is observed that the wind speed attained an increased rate for higher altitudes with higher hub height (100 to 160 m) which fetches more wind power. The wind speed of the Gulf of Khambhat during prevailing months is high

Gulf of Khambhat (Gujarat Offshore) Station
The mean wind speed (MWS) collected from LiDAR measurements at various heights (40 to 200 m) are presented in Figure 7 for the period of December 2018 to November 2019. It is observed that the wind speed attained an increased rate for higher altitudes with higher hub height (100 to 160 m) which fetches more wind power. The wind speed of the Gulf of Khambhat during prevailing months is high compared with the Kayathar region, particularly during winter and NEM periods, which are suitable for sustained wind generation.  Table 5 shows the monthly-maximum mean wind speed in July (10.13 m/s) and low mean wind speed in November (4.81 m/s). Considering the 100 m measurements for analysis, the annual mean wind speed and maximum wind speed during June are observed as 7.59 and 22.99 m/s respectively. The shear analysis power-law coefficient is perceived as 0.0782 based on the measured wind speeds at various altitudes. In comparison with the Kayathar onshore shear coefficient, the Gulf of Khambhat is less due to its offshore landscape by nature. The wind rose plot for seasonal periods is plotted and presented in Figure 8. The major wind directions during SW (Southwest) monsoon are observed in the SW-South direction and during the NE (Northeastern) monsoon on the North-NE direction. The annual and seasonal wind rose plots are shown in Figure 8. The maximum range of annual wind speed occurred in SSW and SW direction. Regarding seasonal changes, the winter season took the maximum wind that occurs in the North to NE direction.   Table 5 shows the monthly-maximum mean wind speed in July (10.13 m/s) and low mean wind speed in November (4.81 m/s). Considering the 100 m measurements for analysis, the annual mean wind speed and maximum wind speed during June are observed as 7.59 and 22.99 m/s respectively. The shear analysis power-law coefficient is perceived as 0.0782 based on the measured wind speeds at various altitudes. In comparison with the Kayathar onshore shear coefficient, the Gulf of Khambhat is less due to its offshore landscape by nature. The wind rose plot for seasonal periods is plotted and presented in Figure 8. The major wind directions during SW (Southwest) monsoon are observed in the SW-South direction and during the NE (Northeastern) monsoon on the North-NE direction. The annual and seasonal wind rose plots are shown in Figure 8. The maximum range of annual wind speed occurred in SSW and SW direction. Regarding seasonal changes, the winter season took the maximum wind that occurs in the North to NE direction.    The month-based wind speed shown peak on July (10.13 m/s) and maximum wind speed occurred during June (22.99 m/s) due to SEM. The low wind speed occurred during October and November, i.e., 5.4 and 4.8 m/s respectively, due to NEM. The skewness factor of the wind speed obtained a maximum value of 1.07 during October. The annual skewness is recorded as 0.25. It indicates that the wind distribution is positive with a moderate skew, and its range falls between 0.5 to 1. The annual kurtosis rate (0.57) indicates that the wind distribution is shorter; tails are thinner than the normal distribution.
The annual wind direction obtained a maximum value at the SW (Southwest) direction of about 15.8% with a mean wind speed of 9.8501 m/s. Further, SSW (South Southwest) and WSW (West Southwest) obtained about 15.8% and 8.7% respectively, due to SWM as shown in Table 7. The next maximum wind direction occurred in the NNE (North Northeast) of about 11.07%, with a mean wind speed of 7.0987 m/s followed by North direction (10.17%) with a mean wind speed of 7.83%. The wind direction starts at N-NNE-NE during the winter season initially, and moves towards the South-Southwest direction during summer with the traces of the SWM season picking up. Later, the peaks have risen during the SWM with SW-S direction and finally through the NEM season in the N-NE direction. The major wind directions are obtained in the S-SW and N-NE direction during the seasonal periods. The turbulence intensity for offshore LiDAR data is plotted in Figure 7b, and its mean value took place at 15 m/s (5.9%). As per the IEC standard, it falls on the turbulence Ç category. On comparison with the onshore Kayathar station, the MTI value is low due to the offshore area. From the above discussed wind characteristics, the Gulf of Khambhat shows better results, with specifically steady wind speed and low turbulence intensity. These positive factors provide sustained wind generation and help the wind farm planners to invest in offshore wind power production in the region of Gulf of Khambhat. The wind data collected from mast anemometers with various altitudes (20-100 m) for the period of December 2018 to November 2019 is considered. The mean wind speed (MWS) data collected from mast measurements at various heights (10 to 100 m) for the same period is presented in Figure 9. The maximum mean wind speed is attained about 8.92 m/s at 100 m in June and a low mean wind speed of 5.11 m/s in November is perceived. Further, the comparison of mean wind speed gain is made while increasing the altitudes (June) at 20 and 100 m are 6.48 m/s, and 8.92 m/s, respectively. Wherein, the 100 m measurements are taken for analysis and the annual mean wind speed of 6.99 m/s at 100 m is observed. The shear analysis power-law coefficient is observed as 0.0228 based on the measured wind speeds at various altitudes. On comparison with the Kayathar onshore shear coefficient, the Jafrabad nearshore is less due to its coastal area. The distance between the Jafrabad wind station and Gulf of Khambhat (offshore) is 25 km. The wind speed behaves dynamically with different altitudes, and some steady wind speed patterns are observed at 100 and 80 m. The wind turbine with the hub height of 80 m shows steady wind generation, and is suitable for efficient power production.  The wind data collected from mast anemometers with various altitudes (20-100 m) for the period of December 2018 to November 2019 is considered. The mean wind speed (MWS) data collected from mast measurements at various heights (10 to 100 m) for the same period is presented in Figure 9. The maximum mean wind speed is attained about 8.92 m/s at 100 m in June and a low mean wind speed of 5.11 m/s in November is perceived. Further, the comparison of mean wind speed gain is made while increasing the altitudes (June) at 20 and 100 m are 6.48 m/s, and 8.92 m/s, respectively. Wherein, the 100 m measurements are taken for analysis and the annual mean wind speed of 6.99 m/s at 100 m is observed. The shear analysis power-law coefficient is observed as 0.0228 based on the measured wind speeds at various altitudes. On comparison with the Kayathar onshore shear coefficient, the Jafrabad nearshore is less due to its coastal area. The distance between the Jafrabad wind station and Gulf of Khambhat (offshore) is 25 km. The wind speed behaves dynamically with different altitudes, and some steady wind speed patterns are observed at 100 and 80 m. The wind turbine with the hub height of 80 m shows steady wind generation, and is suitable for efficient power production.  The annual wind direction from the wind rose plots are shown in Figure 10. The maximum wind direction occurs in SW to West direction. Regarding the winter season and NEM, the wind direction The annual wind direction from the wind rose plots are shown in Figure 10. The maximum wind direction occurs in SW to West direction. Regarding the winter season and NEM, the wind direction goes toward the NNE to East direction. During the SWM season and summer monsoon periods, the majority of wind occurs in the SW to WNW direction.   Tables 8 and 9 display the annual, seasonal, and monthly wind parameters. Table 10 shows the annual wind direction at the SW direction with maximum stake of 15.04% and mean wind speed of 7.991 m/s. Later, the WSW and WNW shares about 11.631% and 8.891% respectively, due to the South-West monsoon. The next maximum wind direction of 10.158% occurred in the NE with a mean wind speed of 7.335 m/s, followed by an NNE direction (9.42%) with a mean wind speed of 7.217%.

Wind Distribution Fitting
For estimating the wind energy potential, accurate wind distribution modeling is essential to extract wind power meritoriously. The performance of targeted location wind potential assessment hinges on the selection of the PDF to illustrate the computed wind speed behavior plotted in the frequency distribution.
The collected wind speed data to fit into ten numbers of wind distribution such as Weibull, Rayleigh, Birnbaum-Sandras (BS), gamma (GM), inverse Gaussian (IG), Rician (RI), lognormal (LN), Nakagami (NK), GEV and bimodal Weibull-Weibull (BM) distribution methods. The bimodal WW method is adopted, since the visual traces of probability distribution function have similarity in the bimodal with two peaks. The bimodal nature is predicted by numerical calculations of mean and standard deviation [33,34]. The goodness of fit and best-fitted wind distribution is ascertained with RMSE and R 2 methods by comparing actual measured data with individual wind distribution methods through PDFs. In this work, three methods are adopted, namely the maximum likelihood method (MLM), expectation-maximization algorithm (EM) and moth flame optimization (MFO). The bimodal-Weibull-Weibull methods are used to estimate the parameters of wind speed distributions. The EM method consists of two steps, namely the expectation-step (E-step) and the maximization-step (M-step). The E-step computes the conditional expectation and the M-step maximizes the expectation. The parameters of the bimodal and multimodal distribution are estimated by an iterative procedure until convergence is reached. Additionally, the bimodal WW parameters (Appendix B) were estimated through the MFO with objective function by minimizing the error between estimated and observed PDF. The distribution of wind speed data (also known as wind frequency distribution) is presented by histogram plots which are a usual method of presenting wind data for a specific year. The histogram displays the percentage of time of occurrence of each wind speed ranges.

Kayathar Station (Onshore)
The annual wind speed observed at the 100 m hub height with estimated parameters using 10 distributions through probability distribution function (PDF) and cumulative probability distribution (CDF) plots are shown in Figure 11. The observed wind speed plotted into frequency distribution is fitting with other distributions to observe the difference in the fitting of each distribution with real histogram data.

Wind Distribution Fitting
For estimating the wind energy potential, accurate wind distribution modeling is essential to extract wind power meritoriously. The performance of targeted location wind potential assessment hinges on the selection of the PDF to illustrate the computed wind speed behavior plotted in the frequency distribution.
The collected wind speed data to fit into ten numbers of wind distribution such as Weibull, Rayleigh, Birnbaum-Sandras (BS), gamma (GM), inverse Gaussian (IG), Rician (RI), lognormal (LN), Nakagami (NK), GEV and bimodal Weibull-Weibull (BM) distribution methods. The bimodal WW method is adopted, since the visual traces of probability distribution function have similarity in the bimodal with two peaks. The bimodal nature is predicted by numerical calculations of mean and standard deviation [33,34]. The goodness of fit and best-fitted wind distribution is ascertained with RMSE and R 2 methods by comparing actual measured data with individual wind distribution methods through PDFs. In this work, three methods are adopted, namely the maximum likelihood method (MLM), expectation-maximization algorithm (EM) and moth flame optimization (MFO). The bimodal-Weibull-Weibull methods are used to estimate the parameters of wind speed distributions. The EM method consists of two steps, namely the expectation-step (E-step) and the maximization-step (M-step). The E-step computes the conditional expectation and the M-step maximizes the expectation. The parameters of the bimodal and multimodal distribution are estimated by an iterative procedure until convergence is reached. Additionally, the bimodal WW parameters (Appendix B) were estimated through the MFO with objective function by minimizing the error between estimated and observed PDF. The distribution of wind speed data (also known as wind frequency distribution) is presented by histogram plots which are a usual method of presenting wind data for a specific year. The histogram displays the percentage of time of occurrence of each wind speed ranges.

Kayathar Station (Onshore)
The annual wind speed observed at the 100 m hub height with estimated parameters using 10 distributions through probability distribution function (PDF) and cumulative probability distribution (CDF) plots are shown in Figure 11. The observed wind speed plotted into frequency distribution is fitting with other distributions to observe the difference in the fitting of each distribution with real histogram data.
The wind distribution from the bimodal Weibull-Weibull method closely follows the peaks of observed histogram wind patterns. The other distribution plots seem to be under fitting and close fit to The wind distribution from the bimodal Weibull-Weibull method closely follows the peaks of observed histogram wind patterns. The other distribution plots seem to be under fitting and close fit to the observed histogram as per the pictorial observation. Figure 12 shows the seasonal wind speed pattern by estimated parameters with distribution functions and plotted with PDF functions. The visual examination of PDF-specifically the Rayleigh distribution-behaves distinctly from other distributions during winter, summer, and SEM monsoon periods. During winter, Rayleigh performs well owing to the incidence of low wind speed. the observed histogram as per the pictorial observation. Figure 12 shows the seasonal wind speed pattern by estimated parameters with distribution functions and plotted with PDF functions. The visual examination of PDF-specifically the Rayleigh distribution-behaves distinctly from other distributions during winter, summer, and SEM monsoon periods. During winter, Rayleigh performs well owing to the incidence of low wind speed. It is also observed that the Rayleigh distribution suffers from histogram fitting during the winter and Southwest Monsoon periods. The bimodal behavior observed during the summer season and the curve follows the histogram variations.

Bi-Modal Behavior
The bi-modal WW method was applied to observe the data in two parts with different weightage, as shown in Table 11. It is estimated using bimodal methods through the EM (expectationmaximization) algorithm and the values are plotted on each part along with weightage, mean speed, standard deviation, shape, and scale parameters. The bmc1 and bmc2 are the bimodality numerical check apart from visual analysis. Where bmc1 = |μ2 − μ1|/(2 (σ1 × σ2) (1/2) ) when bmc1 > 1. Likewise, bmc2 == |μ2 − μ1|/(2* min (σ1, σ2)) when bmc2 > 2. (μ1 and μ2 are means, and σ1 and σ2 are standard deviations of 2 parts of the bimodal model.  It is also observed that the Rayleigh distribution suffers from histogram fitting during the winter and Southwest Monsoon periods. The bimodal behavior observed during the summer season and the curve follows the histogram variations.

Bi-Modal Behavior
The bi-modal WW method was applied to observe the data in two parts with different weightage, as shown in Table 11. It is estimated using bimodal methods through the EM (expectation-maximization) algorithm and the values are plotted on each part along with weightage, mean speed, standard deviation, shape, and scale parameters. The bmc1 and bmc2 are the bimodality numerical check apart from visual analysis. Where bmc1 = |µ2 − µ1|/(2 (σ1 × σ2) (1/2) ) when bmc1 > 1. Likewise, bmc2 == |µ2 − µ1|/(2* min (σ1, σ2)) when bmc2 > 2. (µ1 and µ2 are means, and σ1 and σ2 are standard deviations of 2 parts of the bimodal model. The wind distribution examination in the Kayathar region falls on multiple peaks and differs from normal bell-shaped histogram curves. The bimodal behavior is exclusively analyzed with the bimodal Weibull-Weibull method. Figure 13 shows the bimodal behavior of the annual and winter season. The annual season has more bimodality than winter through visual analysis, as given in Table 11, i.e., bmc1 and bmc2 estimation. Wherein, winter comprises a lesser range of bmc1 and bmc2 of 0.64 and −0.1511 respectively with the annual value of 1.43 and 3.05 for bmc1 and bmc2 respectively. The wind distribution examination in the Kayathar region falls on multiple peaks and differs from normal bell-shaped histogram curves. The bimodal behavior is exclusively analyzed with the bimodal Weibull-Weibull method. Figure 13 shows the bimodal behavior of the annual and winter season. The annual season has more bimodality than winter through visual analysis, as given in Table 11, i.e., bmc1 and bmc2 estimation. Wherein, winter comprises a lesser range of bmc1 and bmc2 of 0.64 and −0.1511 respectively with the annual value of 1.43 and 3.05 for bmc1 and bmc2 respectively.  Table 12 shows the yearly and seasonal parameters from the observed wind speed by ten wind distribution models with the parameter estimation from the maximum likelihood method. The MATLAB developed in-house software program codes are used for parameter estimation.   Table 12 shows the yearly and seasonal parameters from the observed wind speed by ten wind distribution models with the parameter estimation from the maximum likelihood method. The MATLAB developed in-house software program codes are used for parameter estimation.  Table 13 displays the goodness of fit methods of RMSE and R 2 between estimated wind distribution parameters and observed parameters by the ten distributions PDF. The annual wind distribution rate of the Kayathar region shows lower RMSE and R 2 value of 0.004 and 0.9878 respectively. The seasonal wind distribution fitness such as winter (best fit by gamma method), summer (inverse gamma), SWM (Rician method) and winter season (GEV method) are observed. From the Kayathar wind distribution fitting observations, the bimodal Weibull-Weibull wind distribution is best suited, compared to other wind speed distributions on annual wind flow analyses.

Gulf of Khambat (Offshore) Wind Distribution
The annual wind speed observed at 100 m with estimated parameters by 10 distributions through probability distribution function (PDF) and cumulative probability distribution (CDF) plots in the Gulf of Khambhat region is presented in Figure 14. The annual wind speed observed at 100 m with estimated parameters by 10 distributions through probability distribution function (PDF) and cumulative probability distribution (CDF) plots in the Gulf of Khambhat region is presented in Figure 14. The Gulf of Khambhat wind distribution presented in a frequency distribution histogram plot is shown in the above Figure. Some bimodal traces and the bimodal distribution follows the peaks of a histogram, and most of the wind distributions seem to be slightly over fit. As per the observations, all the wind distributions follow the histogram if there are no bimodal traces. On bimodal traces, all distributions are slightly deviating except bimodal wind distributions. Figure 15 shows the seasonal wind speed pattern by the estimated parameters using distribution functions and plotted with PDF functions. The Rician, Weibull and Nakagami seem to follow the histogram. Whereas, the bimodal distribution seems to be a close fit. The Rayleigh distribution shown in the orange line in the PDF plot during the Southwest monsoon deviates from the observed frequency histogram plot. The Rayleigh distribution looks to be fit on low wind speed seasons only. Table 14 shows the bi-modal WW method parameters. The bmc1 and bmc2 are the bimodality numerical check apart from visual analysis. The Gulf of Khambhat wind distribution presented in a frequency distribution histogram plot is shown in the above Figure. Some bimodal traces and the bimodal distribution follows the peaks of a histogram, and most of the wind distributions seem to be slightly over fit. As per the observations, all the wind distributions follow the histogram if there are no bimodal traces. On bimodal traces, all distributions are slightly deviating except bimodal wind distributions. Figure 15 shows the seasonal wind speed pattern by the estimated parameters using distribution functions and plotted with PDF functions. The Rician, Weibull and Nakagami seem to follow the histogram. Whereas, the bimodal distribution seems to be a close fit. The Rayleigh distribution shown in the orange line in the PDF plot during the Southwest monsoon deviates from the observed frequency histogram plot. The Rayleigh distribution looks to be fit on low wind speed seasons only. Table 14 shows the bi-modal WW method parameters. The bmc1 and bmc2 are the bimodality numerical check apart from visual analysis. The annual wind speed observed at 100 m with estimated parameters by 10 distributions through probability distribution function (PDF) and cumulative probability distribution (CDF) plots in the Gulf of Khambhat region is presented in Figure 14. The Gulf of Khambhat wind distribution presented in a frequency distribution histogram plot is shown in the above Figure. Some bimodal traces and the bimodal distribution follows the peaks of a histogram, and most of the wind distributions seem to be slightly over fit. As per the observations, all the wind distributions follow the histogram if there are no bimodal traces. On bimodal traces, all distributions are slightly deviating except bimodal wind distributions. Figure 15 shows the seasonal wind speed pattern by the estimated parameters using distribution functions and plotted with PDF functions. The Rician, Weibull and Nakagami seem to follow the histogram. Whereas, the bimodal distribution seems to be a close fit. The Rayleigh distribution shown in the orange line in the PDF plot during the Southwest monsoon deviates from the observed frequency histogram plot. The Rayleigh distribution looks to be fit on low wind speed seasons only. Table 14 shows the bi-modal WW method parameters. The bmc1 and bmc2 are the bimodality numerical check apart from visual analysis.  The wind distribution analysis in the Gulf of Khambhat region falls on multiple peaks. Figure 15 shows the bimodal behavior of the annual and winter season. The annual season has less bimodality than winter through visual analysis. The estimated value of bmc1 and bmc2 are observed as 1.29 and 2.35 respectively. The annual value has a lower score of bmc1 and bmc2 of about 1.09 and 1.50 respectively. Table 15 shows the yearly and seasonal parameters from the observed wind speed by ten wind distribution models with the parameter estimation from the maximum likelihood method. Table 16 shows the goodness of fit methods of RMSE and R 2 between estimated wind distribution parameters and observed parameters by ten distribution models. The annual wind distribution pattern follows the bimodal-WW wind distribution model with a lower rate of RMSE and R 2 values.   The wind distribution analysis in the Gulf of Khambhat region falls on multiple peaks. Figure 15 shows the bimodal behavior of the annual and winter season. The annual season has less bimodality than winter through visual analysis. The estimated value of bmc1 and bmc2 are observed as 1.29 and 2.35 respectively. The annual value has a lower score of bmc1 and bmc2 of about 1.09 and 1.50 respectively. Table 15 shows the yearly and seasonal parameters from the observed wind speed by ten wind distribution models with the parameter estimation from the maximum likelihood method. Table 16 shows the goodness of fit methods of RMSE and R 2 between estimated wind distribution parameters and observed parameters by ten distribution models. The annual wind distribution pattern follows the bimodal-WW wind distribution model with a lower rate of RMSE and R 2 values.  By comparing the annual best fit, the bimodal WW distribution performed well compared with the other distributions as indicated with the R 2 and RMSE values of 0.98 and 0.004, respectively. During seasonal periods, the bimodal WW performed well on summer and SWM periods. The Weibull and BM share the best fit during the winter season, and the Rayleigh distribution performed well during the North East monsoon (NEM) period. Figure 16 shows the annual wind speed observed at the 100 m hub height with computed parameters by 10 distributions of PDF and CDF plots in the Jafrabad nearshore region. The wind distribution representation of the Jafrabad plot such as inverse Gaussian, gamma, and lognormal distribution looks to over fit the histogram. Whereas, the Weibull, Rician, Nakagami and bimodal distributions appear to fit the frequency distribution plot. The seasonal wind speed pattern of the estimated parameters by the distribution functions and plotted with PDF is shown in Figure 17. The bimodality is observed in the winter season of the Jafrabad plot with peaks of histogram and summer season. The Rayleigh shows the usual trend, which has diverged during the Southwest monsoon. It behaves worse during the high-wind period of the SEM. The other wind distributions, such as the Rician, Weibull and Nakagami, are responding well during the high-wind periods.  Figure 16 shows the annual wind speed observed at the 100 m hub height with computed parameters by 10 distributions of PDF and CDF plots in the Jafrabad nearshore region. The wind distribution representation of the Jafrabad plot such as inverse Gaussian, gamma, and lognormal distribution looks to over fit the histogram. Whereas, the Weibull, Rician, Nakagami and bimodal distributions appear to fit the frequency distribution plot. The seasonal wind speed pattern of the estimated parameters by the distribution functions and plotted with PDF is shown in Figure 17. The bimodality is observed in the winter season of the Jafrabad plot with peaks of histogram and summer season. The Rayleigh shows the usual trend, which has diverged during the Southwest monsoon. It behaves worse during the high-wind period of the SEM. The other wind distributions, such as the Rician, Weibull and Nakagami, are responding well during the high-wind periods.  Figure 16 shows the annual wind speed observed at the 100 m hub height with computed parameters by 10 distributions of PDF and CDF plots in the Jafrabad nearshore region. The wind distribution representation of the Jafrabad plot such as inverse Gaussian, gamma, and lognormal distribution looks to over fit the histogram. Whereas, the Weibull, Rician, Nakagami and bimodal distributions appear to fit the frequency distribution plot. The seasonal wind speed pattern of the estimated parameters by the distribution functions and plotted with PDF is shown in Figure 17. The bimodality is observed in the winter season of the Jafrabad plot with peaks of histogram and summer season. The Rayleigh shows the usual trend, which has diverged during the Southwest monsoon. It behaves worse during the high-wind period of the SEM. The other wind distributions, such as the Rician, Weibull and Nakagami, are responding well during the high-wind periods.  Table 17 shows the bi-modal WW method parameters. The bmc1 and bmc2 are the bimodality numerical check, apart from visual analysis. There is no bimodality in comparison with other locations, due to the lower value of bmc1 and bmc2 of about 0.860 and 0.328, respectively. Figure 17 shows the bimodal behavior of winter season with the bmc1 and bmc2 scores of 1.321 and 1.088, respectively, with not much visual changes on bimodality during the summer season. Table 18 shows the yearly and seasonal estimated parameters from the observed wind speed by ten distribution models with the parameter estimation from the maximum likelihood method.   Table 17 shows the bi-modal WW method parameters. The bmc1 and bmc2 are the bimodality numerical check, apart from visual analysis. There is no bimodality in comparison with other locations, due to the lower value of bmc1 and bmc2 of about 0.860 and 0.328, respectively. Figure 17 shows the bimodal behavior of winter season with the bmc1 and bmc2 scores of 1.321 and 1.088, respectively, with not much visual changes on bimodality during the summer season. Table 18 shows the yearly and seasonal estimated parameters from the observed wind speed by ten distribution models with the parameter estimation from the maximum likelihood method. It is observed that the bimodal WW performed well during the winter periods particularly on the SWM and NEM periods ( Table 19). The value of R 2 (0.989) shared the best fit during summer, whereas the Weibull shares the best fit with RMSE (0.006). By summarizing all the locations, the bimodal WW stands at the top, and the Jafrabad Weibull distribution ranked as the best fit. The annual PDF parameters estimated from the ten distributions are compared with the proposed optimization method using the moth flame optimization algorithm are in the subsequent section.

Optimization Parameters Comparison
The parameters estimated from the wind distribution methods in the previous sections through unimodal and EM algorithm for bimodal (BM-EM) are compared with annual bimodal Weibull parameters through the moth flame optimization method. The estimation of the moth flame method parameter is presented in Table 20. The optimization parameter is projected through the MFO with the objective function for the bimodal Weibull distribution. The wind parameters from the MFO method, PDF fitness analysis with RMSE and R 2 method, and the snapshot of WPD results are presented in Table 21. The MFO parameters are estimated in two stages by satisfying the PDF objective function firstly with minimum RMSE/R 2 value. Then, the over fit and under fit of wind power density calculated from estimated parameters from optimization is tuned with MFO optimization again for higher accuracy. The observed results satisfy both the lower RMSE/R 2 value and the annual wind power density. It is accurately estimated through the MFO method. The results of MFO for Kayathar, Jafrabad, and the Gulf of Khambhat (GoK) outperform the bimodal parameters estimated through the EM algorithm and MFO stands top-graded in all the stations.

Wind Power Density Analysis (WPD)
The WPD value is the major factor for the wind energy potential assessment of a targeted location. It also helps to evaluate the economic feasibility to form a potential wind farm. The modelling of WPD outlines the distributions of wind energy at several wind speed and it is proportional to the cube of wind speed. Its value is represented in W/m 2 and depends on the wind location, air density, and wind speed. It is derived as follows: where P is the wind power measured in watts, A stands for a swept area in m 2 , ρ states the air density (ρ = 1.225 kg/m 3 ), and v represents the wind speed in m/s. Considering all these, the distribution of wind speed is taken into account and then WPD can be expressed as: For evaluating the WPD of a particular distribution, the wind power density distribution function f (v) can be expressed as: The best fit of WPD is calculated as follows: where WPDo, the observed wind power density, WPDest is the estimated values from distribution fitting, and eWPD is the error between observed and estimated values.

Comparison of Wind Power Density
The comparison of annual WPD with estimated MFO parameters and seasonable wind power density for the three locations are given in Tables 22 and 23 respectively. The different landscapes of three locations such as offshore, nearshore, and onshore are considered for the analysis. Wind power density is evaluated after carrying goodness of fit with RMSE and R 2 using probability distribution methods. The goodness of fit indicates the wind power density and demonstrates the best distribution suits for the selected location. Table 24 represent the wind power density calculated from different distributions along with measured value during seasonal periods for three selected locations. Additionally, the annual WPD shares (region wise) of Jafrabad, Gulf Khambhat, and Kayathar are shown in Figure 18. The percentage shares of the SW monsoon take the upper hand in all three sites compared to other seasonal shares. Moreover, Figure 19 shows the monthly wind power density statistics using monthly WPD analysis for Kayathar station. The bimodal method scores better during February, April, July, August, and September (five months), whereas GEV distribution scores best during March, October, November, and December (four months) followed by Weibull during May and June as the best fit. For monthly wind power density of Gulf of Khambhat, BM scores better except February, June, and November whereas Rician distribution scores best during February and November followed by Weibull during June. For Jafrabad region monthly WPD analysis, BM scores better except for in April, May, and August, whereas the Rician distribution scores best during April. The generalized extreme value (GEV) distribution scored better during May, and the Weibull during August was the best fit.  Table 24 represent the wind power density calculated from different distributions along with measured value during seasonal periods for three selected locations. Additionally, the annual WPD shares (region wise) of Jafrabad, Gulf Khambhat, and Kayathar are shown in Figure 18. The percentage shares of the SW monsoon take the upper hand in all three sites compared to other seasonal shares. Moreover, Figure 19 shows the monthly wind power density statistics using monthly WPD analysis for Kayathar station. The bimodal method scores better during February, April, July, August, and September (five months), whereas GEV distribution scores best during March, October, November, and December (four months) followed by Weibull during May and June as the best fit. For monthly wind power density of Gulf of Khambhat, BM scores better except February, June, and November whereas Rician distribution scores best during February and November followed by Weibull during June. For Jafrabad region monthly WPD analysis, BM scores better except for in April, May, and August, whereas the Rician distribution scores best during April. The generalized extreme value (GEV) distribution scored better during May, and the Weibull during August was the best fit.     Figure 19c represents the wind power density share in the different seasons by comparing each region. Gulf of Khambhat (offshore) shares about 37% of maximum wind power density during SWM and minimum during the NEM period with 15%. The winter and summer have a steady share of 22% and 26% respectively. The Gulf of Khambhat (GoK) has the advantage of steady wind power density, due to the advantage of offshore characteristics which is free from wind obstruction comparing with landmass. Jafrabad shows a similar pattern of WPD with peak WPD during the SWM of 40% and minimum on the North East period from October to December (about 14%). During winter and summer, WPD shares about 24% and 26% respectively. Kayathar displays wide marginal variations on WPD share with peaks during the SWM of about 68% and low WPD during North-East monsoon of 10% share. The winter stakes about 7% and summer share about 15% of the rest wind power density. Finally, the wind power class of Kayathar region falls on Class 2, Gulf of Khambhat falls on Class 3 and Jafrabad falls on Class 2 [35].

Research Findings
The wind energy potential assessments for selected locations of different landscape are evaluated using bimodal factors along with ten wind distribution methods and MFO optimization methods. The summary of the outcomes is represented in Table 25. accurate wind distribution fit and wind power density are achieved using the moth flame optimization method (on estimating wind distribution parameters) in comparison with the maximum likelihood method. The selection of technologies for better wind resource assessment is the prime motive for ascertaining the accurate wind potential in targeted location. The commercially available wind analysis software should be tuned by incorporating the latest developments and observation of concepts from these types of research studies to cope the practical field conditions. This research study demonstrates the comparison of different wind distribution methods using different parameter estimation methods. It can be used to select the best wind distribution method from available methodologies. Though the proposed method shows numerous advantages over other methods, limitations such as number of iterations, number of moths, and measured quality of wind speed data are considered before adopting this method. Moreover, the measured data from field instruments are properly calibrated for high accuracy in results. The data collection during wind resource assessment and data loggers from instruments should be closely monitored for ensuring the continuance of data. In nutshell, this research explored the multiple wind distributions and their behavior on different landscapes using bimodality influences. In India, Tamil Nadu and Gujarat stand at the top on wind energy generation. During data collection for this research, certain issues are observed in day to practice on wind energy management. The major hurdle is the low voltage ride through (LVRT) issue [36]. Because the power generating plants must continue to operate through short periods of low grid voltage, and should not get disconnected from the grid. As per the Indian Electricity Grid Code (IEGC), wind farms connected above 66 kV need to stay connected in the grid. During system faults, it acts as a fault recovery by providing reactive power compensation to support the grid. In some of the old wind turbines with stall regulation types, the provision of putting add on LVRT is technically not possible. It is noted that about 11,510 such turbines were installed till April 2014 in India, which have to be replaced with new wind turbines to make a smooth grid operation. The active network management describes the control systems that manage the generation and load for specific purposes to implement in India to control the wind turbines during high grid frequency. However, the SCADA (supervisory control and data acquisition) integration and network communication have to be strengthened. During a higher rate of renewable energy injection in the grid, the control of renewable energy should be adopted to maintain the grid by modernizing the network and information technology systems. The wind resource assessment is based on field data measurements, history of data available in the archive, and satellite data. However, raising standards of data quality through big data analytics [37] will help to make the wind and solar resource assessment with high accuracy.

Conclusions
The wind speed characteristics and distribution patterns in the offshore, onshore, and nearshore in parts of India are compared at different timescales such as monthly, seasonal, and annual in detail to predict the wind behavior. The following conclusions are drawn based on the analysis:

•
In the offshore coastal area of the Gulf of Khambhat, Gujarat, the turbulence intensity is at the low of 0.0782 due to low surface roughness, in comparison to the onshore Kayathar wind station. The prevailing wind direction for the Gulf of Khambhat is observed on the SW (15.8%) and SSW (13.05%) from the Southwest monsoon. The northeast monsoon fetches low wind and prevailing wind direction with a wind speed of about 11.07% (NNE) and North direction of 10.17% wind speed. The WPD density measured at 100 m of Gulf of Khambhat is highest with 431.53 watts/m 2 , in comparison with Kayathar and Jafrabad. The annual mean wind speed is better at about 7.51 m/s.

•
In the onshore area, the Kayathar region, the WPD obtained a maximum of 68% share during southwest monsoon which is highest compared with the Gulf of Khambhat and Jafrabad. During winter, it is perceived as 7%, and the bimodality nature in the annual wind pattern is achieved a great value comparing with the other two regions. The wind power class for Kayathar and Jafrabad belongs to Class 2 categories while the Gulf of Khambhat falls under the Class 3 category.

•
In the nearshore area, the Jafrabad region, the mast is installed at 100 m in line of sight with LiDAR (installed at the offshore location) for correlation and validation. The distance between the two locations is roughly about 25 km. However, the wind pattern has significant variations between the offshore and nearshore, with Jafrabad with a low mean speed of 6.99 m/s in comparison with offshore measurements of 7.51 m/s.

•
The influence of SW (South-west) and NEM (North east monsoon) taken for seasonal analysis. The SWM (South west monsoon) fetches more wind power generation than NE (North east).

•
The conventional unimodal wind distribution models are overridden by the bimodal (Weibull-Weibull) method with the moth flame optimization method. Comparatively, the bimodal WW performed well, followed by the Nakagami, Rician, and Weibull distributions.

•
This research outcome helps with investment in offshore wind farms in Gujarat, due to its higher power density obtained through the proposed method. Moreover, there would be a drop in pricing of offshore wind energy tariffs and land acquisition cost compared to onshore wind farm.

•
The proposed MFO method can be extended to short term wind forecasting by optimizing the weight and bias parameters of the artificial neural network (ANN) method. The optimizing parameters of window size and neurons of long short-term memory (LSTM) would be adopted for accurate predictions.
In future, the concept of big data analysis will be adopted by integrating historical data and atmospheric environmental parameters in the targeted location for prediction of detailed wind resource assessment accurately.