The Possibility of Generating Electricity Using Small-Scale Wind Turbines and Solar Photovoltaic Systems for Households in Northern Cyprus: A Comparative Study

The increased energy demand and related environmental problems caused by burning fossil fuels have raised interest in alternative energy sources. This study investigated the wind characteristics and available wind energy for three urban regions in Northern Cyprus using the Weibull distribution function. The results illustrate that Gazimağusa is the most applicable location for harvesting the kinetic energy of the wind compared to Lefkoşa and Girne. Moreover, the solar potential at a specific location can be analyzed using a different simulation tool. In the present paper, the performance of a rooftop Photovoltaic (PV) system for household buildings in three selected is assessed. Three types of simulation software (PVGIS, PV*SOL, and PVWatts) are used to evaluate the performance of the 6.4 kWp grid-connected rooftop PV system. This study assessed the energy generation, performance ratio and capacity factor for this PV system. The results concluded that PVGIS is an easy, fast, and reliable software tool that can be used for the simulation of a solar PV system in the studied regions. Furthermore, an economic evaluation of renewable systems in the three urban regions is presented. As a result, a small-scale grid-connected solar/wind system that is able to generate electricity with an excellent percentage of clean energy was proposed and developed. The analysis indicates that the proposed PV projects showed significant potential in the studied locations. In addition, the proposed PV system is the most economical option for generating electricity compared to wind systems due to the low electricity prices and of the ability to recover the initial investment. Consequently, it is expected that the simulation results will help in demonstrating the advantages and challenges of installing grid-connected PV systems for households in Northern Cyprus in order to reduce the electricity consumption produced by fossil fuels.


Introduction
In Northern Cyprus, the electricity is currently produced using diesel generator power stations and PV power plant, which have been installed in Serhatköy with capacities of 212 MW and 1.27 MW, respectively [1,2]. Moreover, the growth of the population has led to an increase in energy demand, where nearly all of the energy production is currently dependent on fossil fuels. The increasing demand placed on conventional sources has encouraged the authors to investigate the field of renewable energy sources for generating electricity in the Northern part of Cyprus, especially wind and solar energy.
The rest of the paper is structured as follows: Section 2 presents overall information about the collected wind data, wind data adjustment, and analysis procedure. Section 3 describes the wind speed characteristics at the studied locations and analyzes the wind power densities at different heights to evaluate the wind energy potential in detail. It also discusses the economic evaluation and the performance of small-scale vertical and horizontal axis wind turbines. Furthermore, the solar resource potential, energy generation and performance comparison between the three different simulations tools are presented in Section 3. Section 4 presents the discussions, and Section 5 provides significant conclusions.

Materials and Methods
This section is divided into two parts. In the first part, the statistical analysis of wind speed measured at a height of 10 m at three urban regions in Northern Cyprus is discussed. The Weibull distribution function is used to determine the wind power density at the studied regions. The power law method is utilized to estimate the wind speed at various hub heights. The annual energy outputs, capacity factors and electricity-generated cost were derived for small-scale wind turbines of various sizes and types. In the second part, the solar potential at a specific location is analyzed using a different simulations tool. Therefore, the performance of a rooftop grid-connected PV system for small household buildings in selected regions in the Northern part of Cyprus is assessed. Three kinds of simulation software (PVGIS, PV*SOl and PVWatts) are used to evaluate the performance of the 6.4 kWp grid-connected rooftop PV system. Additionally, the energy generation, performance ratio capacity factor for this PV system and energy cost are calculated. Figure 1 illustrates the procedure analysis of the current study.
of PV systems and wind turbines that can be used to generate electricity at three different locations in Northern Cyprus.
The rest of the paper is structured as follows: Section 2 presents overall information about the collected wind data, wind data adjustment, and analysis procedure. Section 3 describes the wind speed characteristics at the studied locations and analyzes the wind power densities at different heights to evaluate the wind energy potential in detail. It also discusses the economic evaluation and the performance of small-scale vertical and horizontal axis wind turbines. Furthermore, the solar resource potential, energy generation and performance comparison between the three different simulations tools are presented in Section 3. Section 4 presents the discussions, and Section 5 provides significant conclusions.

Materials and Methods
This section is divided into two parts. In the first part, the statistical analysis of wind speed measured at a height of 10m at three urban regions in Northern Cyprus is discussed. The Weibull distribution function is used to determine the wind power density at the studied regions. The power law method is utilized to estimate the wind speed at various hub heights. The annual energy outputs, capacity factors and electricity-generated cost were derived for small-scale wind turbines of various sizes and types. In the second part, the solar potential at a specific location is analyzed using a different simulations tool. Therefore, the performance of a rooftop grid-connected PV system for small household buildings in selected regions in the Northern part of Cyprus is assessed. Three kinds of simulation software (PVGIS, PV*SOl and PVWatts) are used to evaluate the performance of the 6.4 kWp grid-connected rooftop PV system. Additionally, the energy generation, performance ratio capacity factor for this PV system and energy cost are calculated. Figure 1 illustrates the procedure analysis of the current study.

Wind Data Measurement
In the investigation of the wind energy potential and characteristics, three studied locations; namely, Lefkoşa, Girne, and Gazimağusa in Northern Cyprus, hourly wind speed data was collected from the Meteorology Department located in Lefkoşa for various periods (2010-2016). The obtained

Wind Data Measurement
In the investigation of the wind energy potential and characteristics, three studied locations; namely, Lefkoşa, Girne, and Gazimagusa in Northern Cyprus, hourly wind speed data was collected from the Meteorology Department located in Lefkoşa for various periods (2010-2016). The obtained wind data are of the surface type obtained at a height of 10 m. These have been extrapolated to various heights using the Weibull distribution technique and the locations' wind shear exponents. The meteorological information and locations of the studied regions are given in Table 1 and Figure 2, respectively. wind data are of the surface type obtained at a height of 10 m. These have been extrapolated to various heights using the Weibull distribution technique and the locations' wind shear exponents. The meteorological information and locations of the studied regions are given in Table 1 and Figure  2, respectively.

Wind Characteristics Model
Several methods are used for the assessment of wind resources for a specific site based on data measured by meteorological stations. Direct methods such as Weibull and Rayleigh and indirect methods like atmospheric boundary layer wind tunnel testing and numerical simulation with Computational Fluid Dynamics (CFD) are widely used for the assessment of wind resources for a specific region [20].
In general, estimating the wind speed in an urban environment is difficult due to the varying roughness, the drag exerted by surface-mounted obstacles on the flow and the presence of adjacent buildings [20]. Hence, it is conceivable that the lack of accurate approaches for the assessment of wind speed in urban areas is a major impediment to the successful development of micro-scale energy generation [21]. The most dependable method for wind assessment in the urban environment is to directly measure the wind speed in the region, ideally at the position and the height of the proposed wind turbine [20].
In the current study, two parameters Weibull distribution is utilized to evaluate the wind speed characteristics of the studied locations. When analyzing the wind energy potential of the interested region, a function, which gives the probability density ( ) and cumulative distribution ( ) functions for the Weibull distribution, is expressed in Equations (1) and (2)

Wind Characteristics Model
Several methods are used for the assessment of wind resources for a specific site based on data measured by meteorological stations. Direct methods such as Weibull and Rayleigh and indirect methods like atmospheric boundary layer wind tunnel testing and numerical simulation with Computational Fluid Dynamics (CFD) are widely used for the assessment of wind resources for a specific region [20].
In general, estimating the wind speed in an urban environment is difficult due to the varying roughness, the drag exerted by surface-mounted obstacles on the flow and the presence of adjacent buildings [20]. Hence, it is conceivable that the lack of accurate approaches for the assessment of wind speed in urban areas is a major impediment to the successful development of micro-scale energy generation [21]. The most dependable method for wind assessment in the urban environment is to directly measure the wind speed in the region, ideally at the position and the height of the proposed wind turbine [20]. In the current study, two parameters Weibull distribution is utilized to evaluate the wind speed characteristics of the studied locations. When analyzing the wind energy potential of the interested region, a function, which gives the probability density ( f (v)) and cumulative distribution (F(v)) functions for the Weibull distribution, is expressed in Equations (1) and (2) [22][23][24][25]: where c is the scale parameter in m/s and k is the shape factor of distribution. Both equations depend on wind speed, scale parameter and shape factor. Determination of these parameters plays an important role when analyzing the wind characteristics. The maximum likelihood method (MLM) is used to estimate the parameters of the distribution function. This method evaluates the parameters of the distribution function in a time-series format. Generally, this method may require a higher number of iterations when compared to other methods. However, it is very popular among similar studies due to its simplicity. The equations that are used to determine the scale (c) and shape (k) parameters of the distribution are given as [25,26]: where n is used to represent the total number of data points in a particular period of time, and v i is the speed of the wind measured at the interval i. In most cases, an accurate value of the surface roughness coefficient is not readily available or ascertained. Therefore, another approach is to use the Weibull probability function parameter values which are determined at the measured height and extrapolate them to the hub heights using the following expressions [27,28] as: k(z) = k 0 1 − 0.088 ln where c 0 and k 0 are the scale and shape factors, respectively, at the measurement height z 0 , while z is the hub height. The exponent n is defined as [27,28]:

Wind Power Density (WPD)
The wind energy potential at a specific region is evaluated by wind power density (WPD). The air density and wind speed are the major parameters that affect the value of WPD. The WPD can be calculated by measured values (Equation (8)) and PDF values (Equation (9)). It can be estimated as [29] P A = 1 2 ρv 3 (8) Moreover, for a period measurement, the average WPD can be determined using Equation (10) [30].
where P is wind power density in W, A is swept area in m 2 , ρ is air density (ρ = 1.225 kg/m 3 ), f (v) is the probability density function (PDF), P the mean wind power density in W and v is the mean wind speed in m/s.

Most Probable Wind Speed and Wind Speed Carrying Maximum Energy
The wind speeds that are most possible or probable (v mp ) and carry the highest (maximum) energy (v maxE ) are necessary for approximating wind power. The two wind speeds are obtained from the scale and shapes factors, as expressed in Equations (11) and (12) [30].

Wind Speed Variation
In order to determine the energy produced by a wind turbine, the power law model is used to estimate the wind speed at different hub heights [30,31].
where v is the wind speed at the wind turbine hub height z, v 10 is the wind speed at original height z 10 (original height is the measured height, which is 10 m height), and α is the surface roughness coefficient (Equation (14)) [30,31].

Wind Turbine Energy Output
The output power of any turbine can be determined power curve of turbine and wind speed of the specific locations [32]. The power curve of wind turbines can be approximated with a parabolic law as given by [32][33][34] where v i is the vector of possible wind speed at a given site, P(v i ) is the vector of corresponding power of the wind turbine (W), P r is the rated power of the turbine (W), v ci is the cut-in wind speed (m/s), v r is the rated wind speed (m/s) and v co is the cut-off wind speed (m/s) of the wind turbine. Suitable to simulate the power curve of a pitch-controlled wind turbine and to a lesser extent a stall-or a yaw-controlled wind turbine, which do not have a constant power range and thus neglects the power output exceeding P r .
The coefficient of performance (C p ) can be calculated as [32] The total energy generated (E wt ) by the operation of the wind turbine over a period (t) can be determined as [32] Lastly, the capacity factor (CF) of a wind turbine can be estimated as [32]

Economic Analysis
The capital cost of the project, cost of operation and maintenance system and the economic design life of the turbine are the main factors that govern the wind power costs for the specific region [35,36]. The present value of costs (PVC) method is widely used to determine the wind energy cost [37]. It can be expressed as where r is the discount rate, i is the inflation rate, n is the machine life as designed by the manufacturer, C omr is the cost of operation and maintenance, I is the investment summation of turbine price and other initial costs, including provisions for civil work, land, infrastructure, installation and grid integration and S is the scrap value of the turbine price and civil work. The electricity generated cost per kWh (EGC) can be estimated as [37]:

PV System Description
A residential building with a small space available on the rooftop area (roughly 70 m 2 ) is used in this study. Table 2 shows the description of the 6.4 kWp rooftop system used. The proposed PV panel for the 6.4 kWp system is composed of Mono-crystalline cell material. The system is of a fixed stand type and can sufficiently power a household comprised of a small family. The components of the grid-connected solar PV plants are The mounting of PV panels is a major consideration; it is important to ensure that they are mounted at optimal angles according to the site conditions. • Grid connection: includes a sub-station and its components like transformers, net metering systems, protection systems, DC/AC inverter, DC/AC breakers fuses. • DC/AC cables: cables are required for connecting panels, inverter and to the grid.

Solar PV Simulation Software
In this study, the solar resource potentials for the selected towns in Northern Cyprus are taken from the radiation databases available from various software. Table 3 shows the various system software used for assessing the performance analysis of the rooftop solar PV systems.

PVWatts
An open source research tool for performance assessment of PV technology in geographical regions and designed by the National Renewable Energy Laboratory (NREL) allowing a user to perform a simulation with geographical data of the location Location, system, and auxiliary devices requirements.

Wind Speed Characteristics at 10 m Height
Variations of the average monthly wind speeds for each region during the investigation period (2010-2016) are illustrated in Figure 3. In addition, monthly wind speeds for the entire measurement period are shown in Figure 3. It is found that Gazimagusa has the highest wind speed values compared to other regions. The mean monthly wind speeds for Gazimagusa range from 3.7 m/s to 7.2 m/s and the general trend is that the mean wind speed decreases from March to August and then starts to increase afterward for the rest of the year. Furthermore, it is observed that the lowest wind speeds for Lefkoşa are around 1.5 m/s in November and highest values appear in June with 3.5 m/s. The monthly wind speed values illustrate that the minimum and maximum average wind speeds vary between approximately 1.5 m/s and 3.4 m/s for Girne.
Moreover, the mean hourly wind speeds that vary within a 24-h period in the selected regions in Northern Cyprus, are shown in Figure 4. It is evident from the charts in Figure 4 that similar patterns within the 24-h period are observed during the investigation period. The hourly average wind speeds slowly decrease early in the mornings and then start to increase until they reach a peak. After the highest values of the period, wind speeds are observed to decrease in Lefkoşa and Girne, but not in Gazimagusa. The average wind speed in the Gazimagusa region decreases from 1 a.m. to 8 a.m. and shows a sharp increase afterward, where it reaches its maximum value at around 1 p.m. The wind speeds decrease after 1 p.m. in Famagusta until 8 p.m. and the mean values show a marginal increase through the night. Overall, it can be determined from the data in Figure 3 that the coastal areas record maximum average wind speeds late in the afternoon and the minimum value occurs between 4 a.m. and 6 a.m. In contrast, the maximum wind speeds were observed at 2 p.m. and the minimum speeds were between 3 a.m. and 4 a.m. in Lefkoşa, which is the capital city of the Northern part of Cyprus and has the highest building density. Moreover, the mean hourly wind speeds that vary within a 24-h period in the selected regions in Northern Cyprus, are shown in Figure 4. It is evident from the charts in Figure 4 that similar patterns within the 24-h period are observed during the investigation period. The hourly average wind speeds slowly decrease early in the mornings and then start to increase until they reach a peak. After the highest values of the period, wind speeds are observed to decrease in Lefkoşa and Girne, but not in Gazimağusa. The average wind speed in the Gazimağusa region decreases from 1 a.m. to 8 a.m. and shows a sharp increase afterward, where it reaches its maximum value at around 1 p.m. The wind speeds decrease after 1 p.m. in Famagusta until 8 p.m. and the mean values show a marginal increase through the night. Overall, it can be determined from the data in Figure 3 that the coastal areas record maximum average wind speeds late in the afternoon and the minimum value occurs between 4 a.m. and 6 a.m. In contrast, the maximum wind speeds were observed at 2 p.m. and the minimum speeds were between 3 a.m. and 4 a.m. in Lefkoşa, which is the capital city of the Northern part of Cyprus and has the highest building density.  Moreover, the mean hourly wind speeds that vary within a 24-h period in the selected regions in Northern Cyprus, are shown in Figure 4. It is evident from the charts in Figure 4 that similar patterns within the 24-h period are observed during the investigation period. The hourly average wind speeds slowly decrease early in the mornings and then start to increase until they reach a peak. After the highest values of the period, wind speeds are observed to decrease in Lefkoşa and Girne, but not in Gazimağusa. The average wind speed in the Gazimağusa region decreases from 1 a.m. to 8 a.m. and shows a sharp increase afterward, where it reaches its maximum value at around 1 p.m. The wind speeds decrease after 1 p.m. in Famagusta until 8 p.m. and the mean values show a marginal increase through the night. Overall, it can be determined from the data in Figure 3 that the coastal areas record maximum average wind speeds late in the afternoon and the minimum value occurs between 4 a.m. and 6 a.m. In contrast, the maximum wind speeds were observed at 2 p.m. and the minimum speeds were between 3 a.m. and 4 a.m. in Lefkoşa, which is the capital city of the Northern part of Cyprus and has the highest building density.  As can be seen from Figures 3 and 4, the characteristics of the mean wind speed v i may be somewhat different from year to year. It is generally believed that in order to obtain a quantitative, representative and persuasive interpretation of wind characteristics, longer periods of wind data are highly preferable.
Unfortunately, on most occasions, long period wind measurements are often unavailable [38]. Therefore, this study introduces the term relative error (ξ) to account for any possible bias originating from insufficient wind data [39]. The relative error (ξ) for each year is calculated through Equation (21) and plotted in Figure 5.
where v i denotes the annual mean wind speed of the concerned year and v is the mean wind speed during the whole period from 2010 to 2016. Figure 4 shows that during the seven-year period, the maximum relative error occurs in the year 2013 (6.18%) for Gazimagusa. For Lefkoşa and Girne, the maximum relative errors are found to be, 6.55% (2013) and 5.95% (2013), respectively.
Environments 2019, 6, 47 10 of 22 As can be seen from Figures 3 and 4, the characteristics of the mean wind speed vi may be somewhat different from year to year. It is generally believed that in order to obtain a quantitative, representative and persuasive interpretation of wind characteristics, longer periods of wind data are highly preferable.
Unfortunately, on most occasions, long period wind measurements are often unavailable [38]. Therefore, this study introduces the term relative error ( ) to account for any possible bias originating from insufficient wind data [39]. The relative error ( ) for each year is calculated through Equation (21) and plotted in Figure 5.
where vi denotes the annual mean wind speed of the concerned year and ̅ is the mean wind speed during the whole period from 2010 to 2016. Figure 4 shows that during the seven-year period, the maximum relative error occurs in the year 2013 (6.18%) for Gazimağusa. For Lefkoşa and Girne, the maximum relative errors are found to be, 6.55% (2013) and 5.95% (2013), respectively.  Table 4 illustrates the shape (k) and scale (c) values evaluated at each location for the entire wind data obtained for the seven-year period between 2010 and 2016. Additionally, Figure 6 presents the yearly variations of the shape and scale parameters for the seven-year period at the three studied locations. It can be noted that the calculated yearly k parameter does not show significant differences throughout the years for all measurement locations. However, the c parameter shows profound changes during the measurement period. Moreover, the annual wind speed frequency distribution data for each region are presented in Figure 7 for the period 2010 to 2016.   Table 4 illustrates the shape (k) and scale (c) values evaluated at each location for the entire wind data obtained for the seven-year period between 2010 and 2016. Additionally, Figure 6 presents the yearly variations of the shape and scale parameters for the seven-year period at the three studied locations. It can be noted that the calculated yearly k parameter does not show significant differences throughout the years for all measurement locations. However, the c parameter shows profound changes during the measurement period. Moreover, the annual wind speed frequency distribution data for each region are presented in Figure 7 for the period 2010 to 2016.

Wind Speed and Wind Power Density at Various Heights
The optimum wind speed for a typical wind turbine should be equal to or higher than 6.7 m/s. At the same time, it is important to note that wind speeds higher than 11 m/s can be dangerous; therefore, it is not safe to invest in a wind turbine in regions that have a wind speed of more than 11 m/s wind speed during the year [31]. The roughness coefficient is expressed by the exponent α,

Wind Speed and Wind Power Density at Various Heights
The optimum wind speed for a typical wind turbine should be equal to or higher than 6.7 m/s. At the same time, it is important to note that wind speeds higher than 11 m/s can be dangerous; therefore, it is not safe to invest in a wind turbine in regions that have a wind speed of more than 11 m/s wind speed during the year [31]. The roughness coefficient is expressed by the exponent α,

Wind Speed and Wind Power Density at Various Heights
The optimum wind speed for a typical wind turbine should be equal to or higher than 6.7 m/s. At the same time, it is important to note that wind speeds higher than 11 m/s can be dangerous; therefore, it is not safe to invest in a wind turbine in regions that have a wind speed of more than 11 m/s wind speed during the year [31]. The roughness coefficient is expressed by the exponent α, which is associated with the characteristics of the land surface and its value varies between 0.05 and 0.5 [25,31]. The surface roughness values (α) determined by using Equation (14) for the different locations are given in Table 5. In this paper, annual mean wind speeds are estimated for different heights of 30, 50, 80, and 90 m by using the roughness coefficients listed in Table 3. The wind speed increases as one moves higher above the ground and this variation is called the wind shear profile. Figure 8 presents the wind shear profiles in the three different regions that are included in this study. which is associated with the characteristics of the land surface and its value varies between 0.05 and 0.5 [25,31]. The surface roughness values (α) determined by using Equation (14) for the different locations are given in Table 5. In this paper, annual mean wind speeds are estimated for different heights of 30, 50, 80, and 90m by using the roughness coefficients listed in Table 3. The wind speed increases as one moves higher above the ground and this variation is called the wind shear profile. Figure 8 presents the wind shear profiles in the three different regions that are included in this study.  Data collected from each region at a height of 10 m have been extrapolated to the various heights, as shown in Table 6. Among the regions investigated in this study, the maximum estimated power density became prominent in the Gazimağusa region, where the highest density is 288.96 W/m 2 at a height of 90 m (Table 6). According to the results of the current study and the mean Data collected from each region at a height of 10 m have been extrapolated to the various heights, as shown in Table 6. Among the regions investigated in this study, the maximum estimated power density became prominent in the Gazimagusa region, where the highest density is 288.96 W/m 2 at a height of 90 m (Table 6). According to the results of the current study and the mean power density ranges found in the literature [40], all of the locations chosen for investigation indicate poor wind energy potential. Therefore, high capacity wind turbines (MWs) are not feasible to be investigated in these areas. Nevertheless, small-scale wind turbines can be used to gather wind energy potential in these regions.

Economic Analysis of Wind Potential
In the wind turbine industry, the wind turbine is classified into two categories, which are horizontal and vertical axis wind turbines [41]. Horizontal axis wind turbines are the most common turbines used today [42]. The selection of wind turbine is a function of the wind power density of the region and class. It is essential that the wind resources are accurately modelled for site evaluation and sizing of the wind turbine. The amount of electricity that can be produced from the wind turbine depends on the wind speed of the specific region. Therefore, the wind speed measurements of the studied region and the power curve of the selected wind turbine are the most important factors for choosing the best wind turbine for the specific region. The selected wind turbines have been chosen after an overall comparison between different types of wind turbines. The specification of the selected wind turbines is tabulated in Table 7.
In this work, the capacity factor and generated electricity of different types of wind turbine systems that could be placed on the roof of a reference building with different heights are studied.  Table 8 presents the results of the capacity factor (CF) and annual generated electricity (GE) of 12 wind turbine models for all selected locations. It is observed that the values of CF of the wind systems in Gazimagusa are the highest compared to other locations and range from 19.25% to 56.31%, which depends on the type of wind turbine. Additionally, it is observed that the CF values of Alize, Iskra, Eurowind 1 and Eurowind 2 at all locations are very low, which indicates that these models are not suitable for generating electricity for all studied locations.
Moreover, the cost of unit energy per kWh (UCE) based on the PVC method for the wind turbine systems in thestudied locations is presented in Table 8. It is observed that the AWT(2)2000 model has lower values for electricity cost compared to other models. In addition, it is observed that Gazimagusa has the lowest values for electricity cost compared to other regions. Furthermore, the wind system of Gazimagusa is a more economical option for generating electricity because it has higher capacity factor values as well as lower UCE values.

Solar Resource Potential in Three Studied Cities
In general, solar potentials will vary from location to location based on weather conditions. For this study, the solar resource potentials for the selected locations are taken from the radiation databases available from the PVGIS and National Renewable Energy Laboratory (NREL). For the PV plant modeling in the PVGIS simulation tool, the solar radiation considered is sourced from the PVGIS-SARAH database. Figures 9-11 show a comparison of the solar irradiation potential using three different simulation tools for the selected regions. Considering the PVGIS Simulation, the maximum solar irradiation is recorded in July and the minimum value of solar irradiation is obtained in January, as shown in Figure 9. Furthermore, the maximum monthly solar radiation using PV*SOL is obtained in June (224 kWh/m 2 ). The maximum and minimum global irradiation using PVWatts are found in June (201.36 kWh/m 2 ) and in the December (26.25 kWh/m 2 ), respectively (Figure 9).  Based on Figure 10, the solar irradiation at Girne varies from 26.24 to 233 kWh/m 2 . In addition, the maximum and minimum solar irradiation using PVGIS are recorded in July and January, respectively. It is also found that the minimum value of solar irradiation is achieved in December, while the maximum is recorded in June using the PV*SOL and PVWatts simulation tools. Considering the PVGIS simulation tool, the lowest and highest solar irradiation values for Gazimağusa are achieved in July and January with values of 123 kWh/m 2 and 232 kWh/m 2 , respectively, as shown in Figure 11. Additionally, considering the PV*SOL and PVWatts simulation tools, the solar irradiation values are within the range of 72-224kWh/m 2 using PV*SOL and 26.28-201.33kWh/m 2 using PVWatts. Based on Figure 10, the solar irradiation at Girne varies from 26.24 to 233 kWh/m 2 . In addition, the maximum and minimum solar irradiation using PVGIS are recorded in July and January, respectively. It is also found that the minimum value of solar irradiation is achieved in December, while the maximum is recorded in June using the PV*SOL and PVWatts simulation tools.  Based on Figure 10, the solar irradiation at Girne varies from 26.24 to 233 kWh/m 2 . In addition, the maximum and minimum solar irradiation using PVGIS are recorded in July and January, respectively. It is also found that the minimum value of solar irradiation is achieved in December, while the maximum is recorded in June using the PV*SOL and PVWatts simulation tools. Considering the PVGIS simulation tool, the lowest and highest solar irradiation values for Gazimağusa are achieved in July and January with values of 123 kWh/m 2 and 232 kWh/m 2 , respectively, as shown in Figure 11. Additionally, considering the PV*SOL and PVWatts simulation tools, the solar irradiation values are within the range of 72-224kWh/m 2 using PV*SOL and 26.28-201.33kWh/m 2 using PVWatts. Considering the PVGIS simulation tool, the lowest and highest solar irradiation values for Gazimagusa are achieved in July and January with values of 123 kWh/m 2 and 232 kWh/m 2 , respectively, as shown in Figure 11. Additionally, considering the PV*SOL and PVWatts simulation tools, the solar irradiation values are within the range of 72-224 kWh/m 2 using PV*SOL and 26.28-201.33 kWh/m 2 using PVWatts.

Variation of Solar Energy Generation
Variations in annual energy generation are observed from the simulated studies conducted using three different software, namely PVGIS, PV*SOL and PVWatts. The proposed PV rooftop system for the household as per the simulated studies was done using various simulation tools in urban regions in order to reduce the amount of electricity generated using diesel generators and to compare the cost of electricity generation based on solar energy with wind and fossil fuels.
A residential building with a small space available on the rooftop area (roughly 70 m 2 ) is used in this study. This paper proposes a 6.4 kW solar PV rooftop system for the three urban regions in Northern Cyprus.
The monthly variations in energy production and energy consumption using the simulation tools are shown in Figures 12-14. Based on the PVGIS simulation results, it is found that the maximum and minimum energy production values are obtained in July and January, respectively for all selected regions. The annual energy production is found to be 10708 kWh for Lefkoşa, 10631 kWh for Girne and 10,936 kWh for Gazimagusa. The specific energy production values for Lefkoşa, Girne, and Gazimagusa are found to be 1673.125 kWh/kWp, 1661.094 kWh/kWp and 1708.75 kWh/kWp, respectively.

Performance Comparison
Based on PV*SOL, the annual energy consumption of the household is calculated as 3500 kWh. Based on the simulation result, PV supplies 1378 kWh (39.37%) and the remaining 2122 kWh (60.63%) is covered from the grid, as shown in Table 9. Moreover, Tab1e 9 shows a performance comparison of the 6.4 kWp rooftop PV system at Lefkoşa using three simulation software (PVGIS, Considering the PV*SOL and PVWatts Simulation, the maximum energy production at the selected regions is achieved in June followed by July. The minimum potential is observed in the month of December. The annual energy production using PV*SOL is found to be 8040 kWh for Lefkoşa, 6003 kWh for Girne and 5986 kWh for Gazimagusa. In addition, the specific energy production values for Lefkoşa, Girne and Gazimagusa are found to be 1256.25 kWh/kWp, 937.96 kWh/kWp and 935.31 kWh/kWp, respectively. Moreover, considering the PVWatts results, the annual energy production values are found to be 3861.812 kWh for Lefkoşa, 3795 kWh for Girne and 3864.375 kWh for Gazimagusa. Additionally, the specific energy production values for Lefkoşa, Girne, and Gazimagusa are found to be 603.41 kWh/kWp, 592.97 kWh/kWp and 603.81 kWh/kWp, respectively.

Performance Comparison
Based on PV*SOL, the annual energy consumption of the household is calculated as 3500 kWh. Based on the simulation result, PV supplies 1378 kWh (39.37%) and the remaining 2122 kWh (60.63%) is covered from the grid, as shown in Table 9. Moreover, Tab1e 9 shows a performance comparison of the 6.4 kWp rooftop PV system at Lefkoşa using three simulation software (PVGIS, PV*SOL, and PVWatts). The PVGIS software gives the highest annual energy output compared to PV*SOL and PVWatts. In terms of maximum energy output, PVGIS shows that it is in the month of July while PV*SOL and PVWatts show maximum energy output in June. The capacity factor and performance ratio obtained by PVGIS are significantly higher than PV*SOL and PVWatts.

Discussion
With rapid growth in populations, energy consumption is expected to increase in the future. Moreover, the increased energy demand and environmental problems caused by burning fossil fuels have raised interest in alternative energy sources, particularly in urban environments. The aim of this study is to investigate the wind and solar potential at three urban regions in the Northern part of Cyprus.
The results of the collected wind speed data and the analysis show that the annual actual wind speeds at the selected regions are 2.47 for Girne, 2.50 m/s for Lefkoşa and 4.65 m/s for Gazimagusa ( Figure 3). Thus, according to the results obtained at a height of 10 m, Gazimagusa has higher wind speeds compared to the other regions (Lefkoşa and Girne), which is confirmed by [43]. Based on the wind power density classes, the evaluation of the wind resource available in the three selected locations as class 1 wind power sites indicates their suitability for small-scale wind turbines (Table 6). From the perspective of the costs of generating electricity, the AWT(2)2000 model has lower values of electricity cost compared to the other models. Additionally, the wind system for Gazimagusa is a more economical option for generating electricity because it has higher capacity factor values as well as lower energy cost per kWh values.
Furthermore, the solar rooftop PV system is an attractive alternate electricity source for households. In addition, the potential of a solar PV at a given site can be evaluated through software simulation tools [44]. The performance ratio of PVGIS is the highest compared to PV*SOL and PVWatts; our results are thus in good agreement with those found by [44]. Based on the analysis, the maximum and minimum annual energy yields are estimated to be 1708.75 kWh/kWp for Gazimagusa and 1661.09 kWh/kWp for Girne. It is shown that the grid-connected system installations are technically viable energy solutions for the selected area.
Out of three simulation software tools, PVGIS is found to be the most widely used and most effective tool because of its ability to carry out multiple analyses, which makes it easier and faster to evaluate different system configurations. In fact, according to Ref. [45], the performance ratio of the Serhatköy PV Power Plant, which is located in the Northern part of Cyprus, is about 73.7%. Based on the results in Table 9, only PVGIS has a performance ration value close to that of Serhatköy. Therefore, PVGIS can be considered as a reliable software tool for the simulation of a solar PV system. Moreover, the authors performed an economic evaluation of developing gird-connected solar/wind power systems for the urban regions. It is found that the relevant authorities are not willing to invest in this type of technology due to high electricity prices and the poor recovery of initial investment. This type of technology is also not suitable for use in Lefkoşa and Girne due to the low capacity factor. In addition, it is observed that PV systems are the most economical option for generating electricity in comparison to producing electricity by wind turbines because of the lower value of the cost of energy per kWh.
According to Kibris Türk Elektrik Kurumu (2018) [46], the diesel electricity cost is about 0.15$/kWh (0-250 kWh) and 0.17$/kWh (250-500 kWh). The current results show that these systems are the most economical option for generating electricity in comparison to producing electricity by diesel generators because of the lower value of the cost of energy per kWh. The results show that using a solar system for producing electricity in Northern Cyprus could reduce the electricity costs by an average of 70% compared to diesel electricity costs. Moreover, the net fuel savings provided by electricity-generation system with solar systems are nearly three times greater than those obtained with wind systems.

Conclusions and Future Works
Urban wind and solar energy consist of the utilization of wind and solar energy technology in applications in urban and suburban built environments. In this assessment, 2010-2016 wind speed data for three urban regions in Northern Cyprus have been used to investigate the potential for generating electricity from wind energy. The Weibull distribution function was used to analyze the wind speed characteristics in the selected regions. The Power law method was used to estimate the wind speeds at various hub heights. A techno-economic assessment was made for the generation of electricity using 12 wind turbines with various characteristics and types in all the studied regions. Moreover, this paper analyzed the solar energy potential based on three simulation tools. The performance of a grid-connected rooftop solar photovoltaic system in terms of energy generation, performance ratio and capacity factor has been analyzed. The main outcomes of the study are:

•
All the considered regions have annual mean wind speeds above 2 m/s, and the wind power densities range between 11.72 W/m 2 and 68.49 W/m 2 at a height of 10 m. The wind power analysis shows that Gazimagusa is the best location for harvesting wind energy. • AWT(2)2000 model with a power rating of 4 kW has the lowest energy production cost among the considered wind turbine technologies.

•
Vertical axis wind turbines have significant potential to generate electricity for the urban environment.

•
Among the simulation software studied, PVGIS has been demonstrated to be an easy, fast, and reliable software tool for the simulation of solar PV systems in the studied towns.

•
In comparison with solar and wind systems, PV systems are the best option for producing electricity in the studied regions.

•
It is concluded that the PV system project is the most economical option for generating electricity in the studied locations compared with the wind system due to the low electricity prices and recovery of initial investment.
• Based on the overall performance of the rooftop solar PV power plant, it is evident that solar PV power is a feasible solution for supplying power and reducing the large volume of CO 2 emissions emitted into the atmosphere.
Generally, solar radiation can be considered constant in a wide area, but wind is a local resource, and it must be studied in the exact location. The prediction of wind speed in the built environment is difficult due to the varying roughness and the drag exerted by surface-mounted obstacles on the flow, which reduce the wind speed close to the ground. Therefore, it is conceivable that there is a lack of accurate approaches for the assessment of wind speed in urban areas. An interesting area for future study is Computational Fluid Dynamics (CFD), which can be used to predict the wind speed in urban environments (selected regions) for conventional buildings. In addition, the evaluation of potential locations for the installation of small-scale wind turbines in urban areas can be investigated using CFD.
Author Contributions: R.A.Z. and Y.K. analyzed the data. Y.K. and H.G. wrote the paper.
Funding: This research received no external funding.