Investigation of Wind Power Potential in Mthatha, Eastern Cape Province, South Africa

: South Africa is currently grappling with a national energy crisis and the high infrastructure costs associated with expanding the national grid to remote areas. Simultaneously, the government has made substantial efforts to harness renewable energy technologies, particularly wind energy. The average wind speed in a speciﬁc region signiﬁcantly inﬂuences the energy yield from wind turbines. The vast open inland terrains, mountainous regions, and coastal areas in the Northern Cape, Eastern Cape, and Western Cape provinces of South Africa possess the most substantial wind potential. It is imperative to initiate wind energy projects in these provinces to cater to a signiﬁcant portion of the local electricity demand, especially in remote areas disconnected from the national grid. Wind energy generation is inherently stochastic, subject to variations in both time and space. Consequently, it is essential to gain a comprehensive understanding of the local wind patterns to assess the feasibility of utilizing wind resources. In the Eastern Cape Province, the Mthatha area still lags in household electriﬁcation, presenting an opportunity to electrify some households using wind energy. This study aimed to evaluate the wind resource potential for Mthatha area, utilizing data spanning from 2018 to 2023, provided by the South African Weather Services. Two distribution models, the two-parameter Weibull and three-parameter Weibull, were employed to characterize the provided wind data. To determine the parameters associated with each distribution model, two estimation methods, the Maximum Likelihood Method (MLM) and the Method of Moments (MOM), were utilized. The performance of these distribution models was assessed using the Root Mean Square Error (RMSE) statistical indicator. The results showed that Mthatha area predominantly experiences low wind speeds, with an annual average wind speed of 3.30 m/s and an overall wind power density of approximately 48.48 W/m 2 . The prevailing winds predominantly originate from the south and east–southeast directions. Consequently, Mthatha is recommended for stand-alone applications, with the added suggestion of augmented wind turbines for the area.


Introduction
The increase in energy demand due to growing population, a decrease in finite fossil fuels, and an increase in energy-linked climate change has led to an increased need to turn to renewable energy resources, such as solar, biomass, and wind, for electricity generation [1].Wind, as one of the fastest growing renewable energy technologies, witnessed the most significant year-on-year increase of 53% in 2020 [2].Wind is inherently intermittent and random in nature, varying across different time scales, including diurnal patterns, with its magnitude influenced by spatial location and site-specific characteristics [3].Establishing a successful wind turbine that offers affordable energy to end-users while ensuring a competitive return on investment for the investor requires an in-depth understanding of the local wind regime [4].It is important to understand wind resources as they play a pivotal role in turbine site selection and facilitate long-term planning.Understanding the wind speed distribution for a specific site serves as the foundation for estimating the operational parameters of a wind turbine, including the capacity factor and energy output [5].Therefore, the first step when planning to use wind resources for electrifying a particular area is to conduct a wind resource assessment to understand the fluctuating nature of the wind resources [6].
The Eastern Cape Province in South Africa has the lowest household access to electricity, with 82.64% in 2018 [7].Expanding the national grid can be challenging in certain areas of the province because of the remoteness of the areas and the lack of supportive infrastructure.This presents a significant challenge to the government's commitment to providing basic services to all, as outlined in the Reconstruction and Development Programme (RDP) implemented since 1994 [7].Wind energy can substantially contribute to electrifying remote rural areas in a province, such as Mthatha.To use wind resources for electrifying the Mthatha area, it is important to understand the wind regime for the area so that the project developers can make an informed decision in terms of turbine selection and project planning and implementation [8].
Simple measurement of wind speed at a specific site like Mthatha is insufficient for understanding the wind regime for predicting the site's energy potential due to the random nature of wind speed.The fluctuating nature of wind needs an appropriate way for its description, which is provided by probability distributions that are used for modelling the wind speed [9].This provides insight into the likelihood of specific wind speeds occurring at the site and identifies the most frequent wind speed, which is essential for determining the required wind power output.Different methods, such as the Johnson, logistics, lognormal, normal, Rayleigh, Weibull Rayleigh, and Weibull probability distribution functions (PDFs), have been proposed and used [4,10].In most of the literature, the two-parameter Weibull PDF (2WPDF) is preferred due to its flexibility, simplicity, and adaptability to a wide range of data [10,11].However, it proves less effective in regions with a high likelihood of minimal wind.In such scenarios, the three-parameter Weibull PDF (3WPDF) outperforms the 2WPDF [12].
The effectiveness of the WPDF in fitting the data depends on the method used for estimating its scale (c) and shape (k) parameters [11].Various methods, including the Maximum Likelihood Method (MLM), Empirical Method of Lysen (EML), Method of Moments (MOM), Modified MLM (MMLM), Least Squares Method (LSM), WAsP method, Openwind Method (OWM), Empirical Method of Justus (EMJ), and Standard Deviation Method (SDM), have been employed for estimating both the scale and shape parameters [6,10,11,[13][14][15][16]. MMLM and Openwind methods have been found to outweigh the other methods [11,13].The MLM, on the other hand, is typically employed when dealing with data featuring wide variations and frequent low wind speeds [14].
Various performance indicators, including Kolmogorov-Smirnov, the coefficient of determination (R 2 ), Mean Absolute Bias Error (MABE), Root Mean Square Error (RMSE), Relative Root Mean Square (RRMS), Correlation Coefficient (R), Mean Absolute Percentage Error (MAPE), and the Index of Agreement (IOA), are employed to evaluate the efficiency of both scale and shape parameter estimation methods.The selection of the appropriate performance indicator enhances our understanding of the wind regime at a given location.
The significance of wind resource assessment is exemplified by the extensive research conducted both locally and globally.For instance, [11] conducted a wind resource assessment for Upper Blinkwater, a small village in the Eastern Cape Province of South Africa.The researchers found that the average wind speeds at heights of 11, 20, and 30 m above ground level (AGL) are 4.36, 4.61, and 4.78 m/s, respectively, with corresponding wind power densities of 144.99, 171.52, and 192.19 W/m 2 .Predominant winds in the area are from the north-west direction.These results indicated that Upper Blinkwater was wellsuited for standalone applications utilizing small-scale wind turbines.Ref. [17] conducted wind potential assessment at six sites in the Eastern Cape Province and recommended small wind turbine installations in Fort Beaufort, Graaff-Reinet, Bisho, and Grahamstown, while Port Elizabeth was recommended for a large-scale wind turbine installation.However, in another study, Port Elizabeth experienced a mean wind speed exceeding 5 m/s, and thus, small-scale electricity generation projects were recommended [18].
The WPDF is the most used model to predict the wind resource potential.The Weibull parameters used for the prediction of the wind resource potential are estimated using various established methods and compared rather than using just one method [13].Thus, this study used the MLM and MOM methods to predict the 2WPDF and 3WPDF parameters for the Mthatha area to establish the wind resource potential to determine the feasibility of utilising wind power for off-grid household electrification.
The article is arranged into the following major sections: Section 2.1 briefly describes the location, size, and conditions of the Mthatha area.Section 2.2 describes the source of data, the size, and the characteristics of the dataset.Section 2.3 provides the wind speed descriptive statistics.Section 2.4 summarises the modelled PDFs, methods for predicting the 2WPDF and 3WPDF parameters, and estimation methods for performance indicators.Section 3 describes and explains the results in terms of the descriptive statistics of the wind speed, the Weibull PDF, and the wind direction.Section 4 summarises the research findings and provides some recommendations on the use of wind energy in the Mthatha area.

Study Area Description
Mthatha, shown in Figure 1, with an area of 54.97 km 2 with 538.30 households/km 2 , is the main city of the King Sabata Dalindyebo Local Municipality in the Eastern Cape Province of South Africa [19].It is among the most disadvantaged areas in South Africa, experiencing continuous load shedding for over a decade [20].
small wind turbine installations in Fort Beaufort, Graaff-Reinet, Bisho, and Grahamstown, while Port Elizabeth was recommended for a large-scale wind turbine installation.However, in another study, Port Elizabeth experienced a mean wind speed exceeding 5 m/s, and thus, small-scale electricity generation projects were recommended [18].
The WPDF is the most used model to predict the wind resource potential.The Weibull parameters used for the prediction of the wind resource potential are estimated using various established methods and compared rather than using just one method [13].Thus, this study used the MLM and MOM methods to predict the 2WPDF and 3WPDF parameters for the Mthatha area to establish the wind resource potential to determine the feasibility of utilising wind power for off-grid household electrification.
The article is arranged into the following major sections: The "Study area description" section briefly describes the location, size, and conditions of the Mthatha area.The "Meteorological data" section describes the source of data, the size, and the characteristics of the dataset.The "Wind Speed Analysis" section provides the wind speed descriptive statistics.The "Weibull probability density function" section summarises the modelled PDFs, methods for predicting the 2WPDF and 3WPDF parameters, and estimation methods for performance indicators.The "Results and discussion" section describes and explains the results in terms of the descriptive statistics of the wind speed, the Weibull PDF, and the wind direction.The "Conclusions" section summarises the research findings and provides some recommendations on the use of wind energy in the Mthatha area.

Study Area Description
Mthatha, shown in Figure 1, with an area of 54.97 km 2 with 538.30 households/km 2 , is the main city of the King Sabata Dalindyebo Local Municipality in the Eastern Cape Province of South Africa [19].It is among the most disadvantaged areas in South Africa, experiencing continuous load shedding for over a decade [20].

Meteorological Data
For reliable wind resource estimates, long-term observations spanning over 3 decades are considered necessary.In the absence of such long-term data, shorter periods of approximately 10 years can be used to produce acceptable estimates [8].For this study, hourly mean wind speed and direction data from January 2013 to 23 August 2023, measured at 10 m AGL, were provided by the South African Weather Services (SAWS).The

Meteorological Data
For reliable wind resource estimates, long-term observations spanning over 3 decades are considered necessary.In the absence of such long-term data, shorter periods of approximately 10 years can be used to produce acceptable estimates [8].For this study, hourly mean wind speed and direction data from January 2013 to 23 August 2023, measured at 10 m AGL, were provided by the South African Weather Services (SAWS).The data were recorded at Umthatha Wo Meteorological Station (at Lat: −31.5490,Lon: 28.6730) and an elevation of 745 m above sea level.The dataset of 93,288 entries contained 7307 h with no records, replaced by not a number (NAN) and removed before the analysis, which corresponded to 7.83% missing data.Using the hypothesis that the removal of data causes downward biases of annual energy production (AEP), the AEP downward bias was estimated to be approximately 3.35%, a value that remains within the range of data validity [21].

The Wind Speed Analysis Wind Speed Statistics
The data were analysed for the whole period, monthly and yearly.Average wind speed (v) was determined from (1), and its standard deviation (σ) was calculated using (2).Here, v i is the wind speed measurement at hourly time step i, and N is the total number of hours with non-zero wind speed per period.
The most probable wind speed (v mp ), which is the peak of a specific PDF, is the most frequent wind speed given by (3) for the WPDF [22].
where k is a dimensionless Weibull shape parameter depicting the dissymmetry, and c is the scale parameter with dimensions of wind speed.

Weibull Probability Density Function
The PDF (f(v)) provides an indication of the likelihood of finding wind of velocity (v) at a given site.The wind has a 2WPDF, and (f(v|k, c )) is given by (4) [23] for k > 0 and c > 0.
The frequency distribution becomes narrow with high values of k indicating low wind speed variation [24].For k = 2, (4) becomes the Rayleigh PDF [11].The value of c is directly proportional to average wind.The associated cumulative (C) PDF (F(v)) provides the possibility of detecting a wind speed less or equal to v.For 2WPDF, (F(v|k, c )) is given by (5).
The F(v) for the 3WPDF, (f(v|k, c, m )), is given by ( 6) [9].The parameters k and c have the same meaning with 2WPDF, and m ≥ 0 is the location parameter, which defines the starting point of the WPDF curve.This provides the 3WPDF a superior estimation accuracy when compared with the 2WPDF.For c < 1, the probability density of the WPDF approaches infinity as v approaches m.
The corresponding F(v|k, c, m ) is calculated using (7).

Approximating the Weibull Parameters
Different methods are used to estimate WPDF parameters.The MLM is the most commonly used method due to its strong consistency and efficiency.However, some researchers argue that the Method of Moments (MOM) is an equally competent competitor to MLM, particularly in terms of asymptotic properties.For those reasons, both MLM and MOM equations listed in Table 1 have been utilised for estimating the Weibull parameters.

Performance Indicators
To assess efficiency (goodness-of-fit) of the MOM and MLM in estimating the Weibull parameters, different statistical indicators may be used.The aim of this assessment is to show the effectiveness of each WPDF model in fitting the measured wind data.In this study, the popularly used RMSE given in (8) was used [11].
where WPd md,i is the wind power density calculated from the measured data, while WPd WPDF.i represents the wind power density predicted from the Weibull distribution.

Wind Energy Estimation
Wind power density (WPd(v)), which depends on the wind speed frequency distribution and the air density (ρ), was used as the performance indicator of wind resource potential and is calculated using (9) [1].
The mean measured wind power density (WPd m,md (v)) was calculated using (10).
WPd(v) is also calculated using WPDF parameters to provide Weibull wind power density (WPd WPDF (v)), and for the 2WPDF, (11) is used.

Statistical Analysis of the Wind Speed
Table 2 lists the wind monthly statistics for Mthatha in terms of the average wind speed (v); standard deviation (σ); coefficient of variance (CoV); minimum (min), median (med), and maximum (max) wind speeds; range; and kurtosis (kurt) and skewness (skw) of the measured wind speed for Mthatha.The area is dominated by low wind speed.The monthly v shown in Figure 2 varied from 2.6180 m/s in May to 3.8905 m/s in November.The σ varied between 1.4381 m/s in March and 2.1947 m/s in August.The CoV is an important indicator of the variation in wind at a particular site.According to [9] the wind speed at a site is highly variable if the CoV is between 40% and 70%.Mthatha winds are highly variable for all months, except for June, when the winds are extremely variable with a CoV of 71.4614%.The kurt is classified into three categories: normal kurtosis (kut = 3) for standard normal distribution, high kurtosis (kurt > 3), and low kurtosis (kurt < 3) [27].For Mthatha, kurt ranges from 0.9017 (low kurtosis) in February to 10.6159 in April, with the rest of the months lying in the high kurtosis category.The wind speeds distribution has visible tails compared to a normally distributed wind speed.This implies that there are more extreme wind speeds, providing the distribution a "pointier" appearance.The positive skw shown by all the months means that the majority of the values of measured wind speed are greater than the average wind speed, therefore depicting a better wind performance at Mthatha [11].The seasonal wind analysis, provided in Table 3, shows that the v varied from 2.8478 m/s in Autum to 3.6791 m/s in spring.The σ ranged from 1.5251 m/s in Autumn to 2.1210 m/s in Winter.The CoV was between 49.5716 in Summer and 67.3772 in Winter and showed that all the seasons have a highly variable wind speed.All the kurt values are above three, and thus, the wind speed distributions are in the high kurtosis region.The seasonal wind analysis, provided in Table 3, shows that the v varied from 2.8478 m/s in Autum to 3.6791 m/s in spring.The σ ranged from 1.5251 m/s in Autumn to 2.1210 m/s in Winter.The CoV was between 49.5716 in Summer and 67.3772 in Winter and showed that all the seasons have a highly variable wind speed.All the kurt values are above three, and thus, the wind speed distributions are in the high kurtosis region.The daily v for the entire period under observation, provided in Figure 3, shows that the maximum daily average wind speed of 4.6184 m/s occurred on the 1st of November.This is also the month with the maximum monthly average wind speed.The wind speed statistics for the entire period, provided in Table 4, indicates that the v of 3.3009 m/s was observed in Mthatha with a CoV of 56.7342% and a skw of 1.2870, implying that the wind speed is highly variable and right skewed.The observed kurt of 4.9809 implies that wind speed distribution is of high kurtosis.The European Wind Energy Association (EWEA) classifies areas for the installation of large wind turbines according to the wind speed, such as fairly good (6.5 m/s), good (7.5 m/s), and very good (8.5 m/s) [17].Thus, Mthatha area with the v of 3.3009 m/s is not suitable for the installation of large wind turbines.Similar results were found in other areas in the provinces of Upper Blinkwater (v between 4.36 and 4.78 m/s at 11, 20, and 30 m AGL), Bisho, Fort Beaufort, Graaff-Reinet, Grahamstown, Port Elizabeth, and Queenstown (3.1-5.6 m/s at 10 m AGL) where the use of smallscale standalone wind projects was recommended [11,17].This study also recommends small-scale standalone wind applications for the Mthatha area.Internationally, Puná Island's v of 4.8 m/s was found to be marginally good for a wind turbine installation.Wind turbines with large rotor diameters of greater or equal to 100 m were recommended for electricity production.
Mean wind speed (m/s) The daily v for the entire period under observation, provided in Figure 3, shows that the maximum daily average wind speed of 4.6184 m/s occurred on the 1st of November.This is also the month with the maximum monthly average wind speed.The wind speed statistics for the entire period, provided in Table 4, indicates that the v of 3.3009 m/s was observed in Mthatha with a CoV of 56.7342% and a skw of 1.2870, implying that the wind speed is highly variable and right skewed.The observed kurt of 4.9809 implies that wind speed distribution is of high kurtosis.The European Wind Energy Association (EWEA) classifies areas for the installation of large wind turbines according to the wind speed, such as fairly good (6.5 m/s), good (7.5 m/s), and very good (8.5 m/s) [17].Thus, Mthatha area with the v of 3.3009 m/s is not suitable for the installation of large wind turbines.Similar results were found in other areas in the provinces of Upper Blinkwater (v between 4.36 and 4.78 m/s at 11, 20, and 30 m AGL), Bisho, Fort Beaufort, Graaff-Reinet, Grahamstown, Port Elizabeth, and Queenstown (3.1-5.6 m/s at 10 m AGL) where the use of small-scale standalone wind projects was recommended [11,17].This study also recommends smallscale standalone wind applications for the Mthatha area.Internationally, Puná Island's v of 4.8 m/s was found to be marginally good for a wind turbine installation.Wind turbines with large rotor diameters of greater or equal to 100 m were recommended for electricity production.Figure 4 shows the daily cycle of wind speed averaged over a single day.The average hourly wind speed curve forms a bell shape.The wind speeds are low during the night, with v ranging from 2.8148 m/s at 22:00 hours to 2.4434 m/s around 07:00 h.During the day, the wind speed increases sharply to attain mean wind speed values of greater than 3.500 m/s between 13:00 and 19:00 h.From the figure, it is deduced that the Mthatha area is windy from mid-morning to early evening.4 shows the daily cycle of wind speed averaged over a single day.The average hourly wind speed curve forms a bell shape.The wind speeds are low during the night, with v ranging from 2.8148 m/s at 22:00 hours to 2.4434 m/s around 07: 00 h.During the day, the wind speed increases sharply to attain mean wind speed values of greater than 3.500 m/s between 13:00 and 19:00 h.From the figure, it is deduced that the Mthatha area is windy from mid-morning to early evening.

Weibull Probability Density Function
The 2WPDF and the 3WPDF were used to fit the Mthatha SAWS wind speed data.The MOM estimated parameters kMOM2 and cMOM2 and kMOM3, cMOM3, and mMOM3 for the 2WPDF and the 3WPDF, respectively, and are listed in Table 5 for the entire period under observation and are categorised as per its seasons.Also listed are corresponding MLM estimated parameters kMLM2 and cMLM2 and kMLM3, cMLM3, and mMLM3 for the 2WPDF and the 3WPDF, respectively.The measured (WPd m,md ) and estimated wind power densities WPd 2WPDFMOM , WPd 3WPDFMOM , WPd 2WPDFMLM , and WPd 3WPDFMLM are provided for the MOM and MLM estimation methods.In both the estimation methods, the 3WPDF achieves superior results than the 2WPDF.The 3WPDF's wind power density estimated using the MLM method is higher than the observed value.This could be due to an error in the parameter estimation.Further studies are recommended to investigate this anomaly.The WPd m,md ranged from 0.0297 kW/m 2 in Autumn to 0.0608 kW/m 2 in Spring.This corresponds to the wind resource temporal distribution, which shows the lowest and the highest mean wind speeds during these seasons.The United States classifies wind resources as fair (WPd < 0.1 kW/m 2 ), moderate (0.1 ≤ WPd < 0.3 kW/m 2 ), good (0.3 ≤ WPd < 0.7 kW/m 2 ), and excellent (WPd > 0.7 kW/m 2 ) [28].Also, the Battelle-Pacific Northwest Laboratory (Battelle-PNL) classifies areas according to their power density at 10 m above the ground, with class 1 (WPd ≤ 0.1 kW/m 2 ) being the lowest and class 7 (WPd ≤ 1 kW/m 2 ) being the highest [29].The Mthatha area, with a WPdof < 0.1 kW/m 2 in all seasons, has a fairly good wind resource and is recommended for small-scale wind turbines in standalone applications.From the conducted wind resource assessments in the Eastern Cape Province, small wind turbines were rec-ommended for measured wind speeds of up to 30 m AGL [11,17,18,30].More wind energy can be harnessed using wind turbines with large rotor diameters or by increasing the hub height.The performance of the parameter estimation methodology was compared in terms of how well the 2WPDF and 3WPDF distributions matched the SAWS measured wind speed data using RMSE.The RMSE for the distributions is provided in Table 5.For the whole period and its seasons, all distributions have an RMSE close to 0, which shows that they are good.The MOM 3WPDF methods provided the best estimates in all cases.This agrees with the PDF fits shown in Figures 5 and 6 for the whole period and its seasons, respectively.
Appl.Sci.2023, 13, x FOR PEER REVIEW 10 of 14 whole period and its seasons, all distributions have an RMSE close to 0, which shows that they are good.The MOM 3WPDF methods provided the best estimates in all cases.This agrees with the PDF fits shown in Figure 5 and Figure 6 for the whole period and its seasons, respectively.From the cumulative probability density plots, for the whole observation period, 57.5707% of the period had a wind speed of more than 3 m/s in Mthatha.In most cases, micro and small wind turbines have a cut-in wind speed of 3 m/s [2].This implies that, for 57.5707% of the observed time, some wind energy was produced in the Mthatha area.From the cumulative probability density plots, for the whole observation period, 57.5707% of the period had a wind speed of more than 3 m/s in Mthatha.In most cases, micro and small wind turbines have a cut-in wind speed of 3 m/s [2].This implies that, for 57.5707% of the observed time, some wind energy was produced in the Mthatha area.

Wind Direction for the Mthatha Area
The wind directions were analysed for the whole period and its seasons.The variation in the annual wind direction for the period from January 2013 to August 2023 is shown in Figure 7.The most prevailing direction is east-southeast (ESE) with a frequency of 12.0877%.For 36.5201% of the time, the wind speed was between 2 and 4 ms −1 in varying directions, while for 13.3845% of the time, null wind speed was observed.Only 0.0035% of the period had a wind speed of above 16 m/s, with 0.0011% in the south-southwest (SSW) direction and 0.0023% in the westerly (W) direction.The wind direction statistics are shown in Table 6 in terms of the period, most prevailing wind direction (MPWD), MPWD total frequency (MPWD tf ), most frequent wind speed range (MFWSR), MFWSR frequency (MFWSR f ) null wind frequency (NWF), lowest speed range (LSR), LSR frequency (LSR f ), highest speed range (HSR), and HSR frequency (HSR f ).

Conclusions
South Africa's Eastern Cape Province has the lowest household access to electricity, with 82.64% in 2018, and has experienced a prolonged load shedding for more than a decade.Grid extension in some remote areas of the province is hampered by the remoteness of certain areas and the lack of supportive infrastructure.Thus, there is a need to promote renewable energy technologies like wind to facilitate standalone and mini grid systems to alleviate load shedding and provide electricity to remote areas.This study investigated the wind potential for Mthatha, which is an important step in determining the country's wind energy potential.The 2WPDF and 3WPDF were used for modelling the wind speed for the period from January 2013 to August 2023.The MLM and MOM were applied to obtain the Weibull parameters, and their performance for fitting the data was compared using the RMSE.It was established that the Mthatha area is dominated with low wind speeds, with an annual wind speed average of 3.3009 m/s.The probability of receiving wind speeds of above 3 m/s was found to be 57.5707%.The area is not suitable for most wind turbines that are available on the wind energy market, which are designed for wind speeds of greater than 5 m/s [31,32].It is recommended that wind speed augmentation mechanisms be employed in the area, if wind turbines are to be used for providing energy.

Nomenclature c
Weibull

Figure 2 .
Figure 2. Monthly average speed variation for Mthatha at 10 m AGL.

Figure 2 .
Figure 2. Monthly average speed variation for Mthatha at 10 m AGL.

14 Figure 3 .Table 4 .Figure 3 .
Figure 3. Average daily wind speed for the entire period at 10 m AGL.Table 4. The wind speed statistics for the period from January 2013 to August 2023 at 10 m AGL.  CoV Min Med Max Range kurt skw 3.3009 1.8727 56.7342 0.1000 2.8000 17.2000 17.7000 4.9809 1.2870

Figure 3 .
Figure 3. Average daily wind speed for the entire period at 10 m AGL.
Mean wind spMean wind speed (m/s)

Figure 5 .Figure 5 .
Figure 5. WPDFs superimposed on to the measured data histogram for the period of 2013-2023.

Figure 5 .
Figure 5. WPDFs superimposed on to the measured data histogram for the period of 2013-2023.

Figure 6 .
Figure 6.WPDFs superimposed on to the measured data histogram for the four seasons of the year.

Table 1 .
Methods for estimating Weibull parameters.

Table 2 .
The monthly statistical wind data for Mthatha at 10 m AGL.

Table 3 .
The seasonal wind characteristics at 10 m AGL.

Table 4 .
The wind speed statistics for the period from January 2013 to August 2023 at 10 m AGL.

Table 4 .
The wind speed statistics for the period from January 2013 to August 2023 at 10 m AGL.

Table 5 .
Estimated Weibull parameters and wind power density for the whole period and its seasons at 10 m AGL.

Table 6 .
Wind direction statistics of Mthatha for the whole period and its seasons.