Multi-Criteria GIS-Based Analysis for Mapping Suitable Sites for Onshore Wind Farms in Southeast France

: Wind energy is critical to traditional energy sources replacement in France and throughout the world. Wind energy generation in France is quite unevenly spread across the country. Despite its considerable wind potential, the research region is among the least productive. The region is a very complicated location where socio-environmental, technological, and topographical restrictions intersect, which is why energy production planning studies in this area have been delayed. In this research, the methodology used for identifying appropriate sites for future wind farms in this region combines GIS with MCDA approaches such as AHP. Six determining factors are selected: the average wind speed, which has a weight of 38%; the protected areas, which have a relative weight of 26%; the distance to electrical substations and road networks, both of which have a signiﬁcant inﬂuence on relative weights of 13%; and ﬁnally, the slope and elevation, which have weights of 5% and 3%, respectively. Only one alternative was investigated (suitable and unsuitable). The spatial database was generated using ArcGIS and QGIS software; the AHP was computed using Excel; and several treatments, such as raster data categorization and weighted overlay, were automated using the Python programming language. The regions identiﬁed for wind turbines installation are deﬁned by a total of 962,612 pixels, which cover a total of 651 km 2 and represent around 6.98% of the research area. The theoretical wind potential calculation results suggest that for at least one site with an area bigger than 400 ha, the energy output ranges between 182.60 and 280.20 MW. The planned sites appear to be suitable; each site can support an average installed capacity of 45 MW. This energy beneﬁt will fulﬁll the region’s population’s transportation, heating, and electrical demands.


Introduction
Increased energy consumption in developed and developing countries as a result of prolonged economic growth [1] may lead to fast resource depletion, environmental degradation, biodiversity loss, and climate change [2][3][4]. Therefore, governments are required to focus their efforts on reducing greenhouse gas emissions and other environmental, social, and economic problems [5][6][7], as well as converting energy supply to green energy production methods.
Indeed, France generates electricity, heat, and transportation using a variety of energy sources. This energy mix includes nuclear, fossil fuels, and renewable sources. According to the Minister of Environmental Transition's key energy numbers for 2021/2020, nuclear accounts for 40% of energy, oil for 28%, and natural gas for 16%. However, in France, renewable energy accounts for just 13% of total energy consumption [8]. Due to its support Indeed, France generates electricity, heat, and transportation using a variety of energy sources. This energy mix includes nuclear, fossil fuels, and renewable sources. According to the Minister of Environmental Transition's key energy numbers for 2021/2020, nuclear accounts for 40% of energy, oil for 28%, and natural gas for 16%. However, in France, renewable energy accounts for just 13% of total energy consumption [8]. Due to its support for activities concerned with the energy and ecological transition, as well as its obligations to reduce the dangers associated with global warming, France would want to see the 30% renewable energy target reached by 2030 [8].
Previous investigations have demonstrated that wind power is one of the most promising renewable energy sources [9][10][11][12]. It is becoming increasingly popular worldwide due to its several benefits, including simple access to efficient multi-megawatt wind turbines [13]. Furthermore, wind energy sources will supply more power than any other form of energy source by 2050 in the European Union's renewable energy decarbonization scenario [14]. In addition, due to variables such as the availability of stronger and longerlasting winds and land for installation, wind energy has recently become an important component of France's increasing renewable energy sector. Currently, there are already 11,625 onshore wind turbines in France. Thus, wind power has increased its proportion of the country's energy output from 2.2% 10 years ago to 7.9% in 2020, up from 6.1% in 2019 (https://www.revolution-energetique.com/, accessed on 2 July 2022). In addition, most wind turbines are unevenly spread across the Hauts-de-France and Grand-Est regions (as illustrated in Figure 1). Nonetheless, some regions, such as Aquitaine, Auvergne-Rhône-Alpes, and Provence-Alpes-Côte d'Azur, have insufficient infrastructures. Many studies have reported significant environmental, social, economic, political, legal, and technological issues associated with wind farms sitting around the world [16]. To address these limits, geographic information systems (GIS) and hierarchical multi-criteria analysis (AHP) methodologies have been frequently adopted. These methods have recently been applied in various studies, including the identification of suitable sites for sitting solar farms [17][18][19][20][21] and suitable sites for marine wind farms [16,[22][23][24][25][26][27].
Unfortunately, no study has been conducted in France on the application of multicriteria analysis approaches combined with GIS. However, there are studies that have been conducted on territories in Europe such as Greece [28,29] or southern Spain [30].
Different scenarios were considered. The result is a site suitability index map ranging from inadequate to highly suitable, which seems impractical for a sensitivity analysis of the overall site suitability index. Uncertainty is often inherent in such cases and leads to decision-making problems and inconsistency between decision makers' preferences [31]. Many studies have reported significant environmental, social, economic, political, legal, and technological issues associated with wind farms sitting around the world [16]. To address these limits, geographic information systems (GIS) and hierarchical multi-criteria analysis (AHP) methodologies have been frequently adopted. These methods have recently been applied in various studies, including the identification of suitable sites for sitting solar farms [17][18][19][20][21] and suitable sites for marine wind farms [16,[22][23][24][25][26][27].
Unfortunately, no study has been conducted in France on the application of multicriteria analysis approaches combined with GIS. However, there are studies that have been conducted on territories in Europe such as Greece [28,29] or southern Spain [30].
Different scenarios were considered. The result is a site suitability index map ranging from inadequate to highly suitable, which seems impractical for a sensitivity analysis of the overall site suitability index. Uncertainty is often inherent in such cases and leads to decision-making problems and inconsistency between decision makers' preferences [31]. All this research has highlighted the capabilities of GIS-based multi-criteria analysis approaches to site selection for onshore wind farms while considering regulations, legislation, and other constraints.
The results vary from one study area to another depending on the area of the study area, its topography, its natural resources (wind, temperature, etc.), and the criteria chosen for the study and the weights assigned to the criteria. We conducted this study with the aim of identifying suitable sites for the planning of future wind farm construction projects in an area that is extremely complex due to the existing constraints in the region, considering environmental or topographical factors. The proposed methodological approach can be applied to any region of the world by adapting the characteristics considered. The implementation of the proposed methodology could facilitate the achievement of national objectives in the energy sector and encourage energy interdependence between many geographical areas in France.
Thus, our study area choice is influenced by previously mentioned reasons, such as the lack of wind farms in this region, which has a high population and high-power consumption (including heating), as well as environmental and relief limits that make it difficult to find suitable sites. The research aims are as follows: (1) to promote the use of GIS-based multi-criteria analysis methodologies in decision-making processes; (2) to contribute to the country's growth by providing cartographic and documentary materials related to wind projects; (3) to provide users with a Python code template that combines each component of the multi-criteria analysis technique for choosing potential onshore wind farm locations. By replacing the criteria, this code can be used for any research that involves decision making. Some of the pre-processing activities must be performed using GIS software such as QGIS or ArcGIS because they are not included in the source code.
It is hoped that this research will contribute to France's efforts in spatial energy planning. Effective wind farm siting options, as outlined in this study, could help the state meet its energy goals and policies.

Study Area
The research area is in southeast France (Figure 2a), which is part of the Provence-Alpes-Côte d'Azur region and covers 80% of the Var department. Moreover, it is bounded to the west by the Bouches-du-Rhône department, to the north by the Alpes-de-Haute-Provence department, to the east by the Alpes-Maritimes department, and to the south by the Mediterranean Sea ( Figure 2b). It has an area of 11,208 km 2 with a perimeter of around 424 km (Figure 2c). The study area is known for its Mediterranean climate, of the Trewartha Cs or Köppen Csa type on the coast. Although it is a maritime climate, the annual temperature range is between 11 and 14 °C (Figure 3). The average annual rainfall is between 45 and 95 mm (Figure 3). The dominant winds are the Mistral (especially in Provence) and the Tramontane (especially in Languedoc) whose power comes from the channeling effect of  The study area is known for its Mediterranean climate, of the Trewartha Cs or Köppen Csa type on the coast. Although it is a maritime climate, the annual temperature range is between 11 and 14 • C ( Figure 3). The average annual rainfall is between 45 and 95 mm ( Figure 3). The dominant winds are the Mistral (especially in Provence) and the Tramontane (especially in Languedoc) whose power comes from the channeling effect of the surrounding massifs to the north and west (Alps, Pyrenees, and Massif Central). Generally, these winds dry the air and clear the sky, and their intensity is very variable from one place to another, depending strongly on the sheltering or accelerating effect of the neighboring massifs [32]. In recent years, the average annual wind speed has been between 5.2 and 5.9 m/s ( Figure 3). Interannual mean temperature ( • C) and precipitation (mm) data were collected from the Climatic Research Unit Time Series (CRUTS) database at the University of East Anglia (CRU TS v. 4.01, https://www.cru.uea.ac.uk/, accessed on 10 July 2022, Harris et al., (2014)). Wind speed data were downloaded from the "Power Data Access" site via the link (https://www.uea.ac.uk/web/groups-and-centres/climatic-research-unit/ data, accessed on 15 July 2022).  France's southeast is one of the country's most heavily populated regions. According to INSEE's 2019 census statistics (Figure 4), the research's area population distribution is heterogeneous, reaching around 5200 inhabitants/km 2 in the south (along the coastline), particularly in important towns such as Marseille, Toulon, and Aix-en-province. In contrast, population density in the region's north is modest, with values lower than 100 inhabitants/km 2 . Consequently, the high population density in southeastern France leads to high energy and electricity usage. However, with only one wind farm, as indicated in Figure 1, satisfying electricity needs using wind is not possible, prompting us to identify potential places for further wind project execution. France's southeast is one of the country's most heavily populated regions. According to INSEE's 2019 census statistics (Figure 4), the research's area population distribution is heterogeneous, reaching around 5200 inhabitants/km 2 in the south (along the coastline), particularly in important towns such as Marseille, Toulon, and Aix-en-province. In contrast, population density in the region's north is modest, with values lower than 100 inhabitants/km 2 . Consequently, the high population density in southeastern France leads to high energy and electricity usage. However, with only one wind farm, as indicated in Figure 1, satisfying electricity needs using wind is not possible, prompting us to identify potential places for further wind project execution.

Methodology
To map suitable sites for wind farm construction in southeast France, we adopted the methodology presented in Figure 5.

Methodology
To map suitable sites for wind farm construction in southeast France, we adopted the methodology presented in Figure 5.

Data Collection
This research's approach mobilizes the whole set of data that determines onshore wind project placement planning. The data collected covers socioeconomic, environmental, and technical parameters (Table 1).

Data Collection
This research's approach mobilizes the whole set of data that determines onshore wind project placement planning. The data collected covers socioeconomic, environmental, and technical parameters (Table 1).  The average wind speed raster data with a spatial resolution of 300 m was obtained from the "Global wind speed" website (https://globalwindatlas.info/, accessed on 15 July 2022). These data are based on ten years of hourly measurements recorded at a height of 100 m (2001-2010). Subsequently, the digital terrain model (DTM) retrieved from the USGS website was employed to generate a mosaic of two DTM rasters covering the whole research region with a spatial resolution of 30 m. Moreover, the IGN 2021 topo database [35] was used to collect information on protected areas (urbanized areas, industrial or commercial areas, infrastructures and equipment, continental waters), road networks (departmental, national, highways, railroads), and electrical substations. All data are resampled in 26 m.

GIS-Based Spatial Database Creation
Elevation is an important criterion; however, in numerous studies, high altitudes have not been indicated for wind projects [36][37][38]. The researchers mentioned have proposed that locations below 1000 m be considered extremely appropriate for wind projects.
This is due to access issues and a lack of basic infrastructure in higher places. Thus, to reduce the high expenses connected with construction, regions lower than 1000 m in height appear to be the most efficient and appropriate. As shown in Figure 6a, elevations above 1000 m account for less than 30% of the research region's total area.  Wind speed: Average wind speed has been the most essential and weighted parameter in wind farm location evaluation studies, as reported in most previous studies [24,28,39,40]. This parameter is directly related to the project's profitability [41]. In our wind farm siting analysis, areas, sites having an average annual wind speed of less than 5 m/s at a height of 100 m above mean sea level were considered inappropriate for wind farm sitings, as recommended by [28,42]. Nonetheless, several studies suggest that an annual average wind speed of more than 6 m/s is required for a functional wind farm installation [43,44]. Conversely, extremely high wind speeds can damage the wind turbines and the project execution in general.
Slope: A slope map in degrees is produced by combining two SRTM rasters acquired from the USGS website. This criterion can be applied to exclude areas with steep slopes of greater than 15 degrees and high relief. These are typically inaccessible and so unsuitable for wind turbines. The highest slopes, as illustrated in Figure 6c, are in the research's area northeast, towards the province of Alpes de Haute, and also surround the shoreline in the south.
Indeed, our choice of 15% (maximum limit of suitable slopes) has already been defined by research works [45,46]; others have adopted a constraint of 25% [47], while some [48,49] have raised the constraint threshold to 30%. In addition, some researchers have considered areas with slopes greater than 10% as infeasible areas for wind turbine installation [50,51]. Selection of land having a slope of less than 15% is planned to facilitate crane and truck accessibility to sites and to reduce installation and maintenance costs due to turbulence.
Protected area: Wind farm construction is controlled by various laws, most notably the French energy code, the urban planning code, and the environmental code. Any prospective wind project needs to evaluate its environmental impact by including parameters such as landscape impact, biodiversity, noise, and dangers to nearby inhabitants. The protected areas in this study, which include urban areas, wetlands, biodiversity parks, and water surfaces (Figure 7), were gathered from the IGN's BD 2021. Slope: A slope map in degrees is produced by combining two SRTM rasters acquired from the USGS website. This criterion can be applied to exclude areas with steep slopes of greater than 15 degrees and high relief. These are typically inaccessible and so unsuitable for wind turbines. The highest slopes, as illustrated in Figure 6c, are in the research's area northeast, towards the province of Alpes de Haute, and also surround the shoreline in the south.
Indeed, our choice of 15% (maximum limit of suitable slopes) has already been defined by research works [45,46]; others have adopted a constraint of 25% [47], while some [48,49] have raised the constraint threshold to 30%. In addition, some researchers have considered areas with slopes greater than 10% as infeasible areas for wind turbine installation [50,51]. Selection of land having a slope of less than 15% is planned to facilitate crane and truck accessibility to sites and to reduce installation and maintenance costs due to turbulence.
Protected area: Wind farm construction is controlled by various laws, most notably the French energy code, the urban planning code, and the environmental code. Any prospective wind project needs to evaluate its environmental impact by including parameters such as landscape impact, biodiversity, noise, and dangers to nearby inhabitants. The protected areas in this study, which include urban areas, wetlands, biodiversity parks, and water surfaces (Figure 7), were gathered from the IGN's BD 2021.
The study area is characterized by the presence of forests, pastures, beautiful landscapes, biodiversity parks (fauna and flora), NATURA 2000 protected sites [52,53], and the most important large urban agglomerations. These areas of environmental interest were not absolutely excluded according to the literature but also according to the national legislation (the minimum distances were determined after the decision approving the environmental conditions ("DAEC")). To avoid the destruction of these spaces and the negative impact of wind farms on the nature of these areas, a minimum distance of 2000 m is required [54,55].  The study area is characterized by the presence of forests, pastures, beautiful landscapes, biodiversity parks (fauna and flora), NATURA 2000 protected sites [52,53], and the most important large urban agglomerations. These areas of environmental interest were not absolutely excluded according to the literature but also according to the national legislation (the minimum distances were determined after the decision approving the environmental conditions ("DAEC")). To avoid the destruction of these spaces and the negative impact of wind farms on the nature of these areas, a minimum distance of 2000 m is required [54,55].
Road network: This was generated using data from the IGN 2021 database (Figure 6d). The road's proximity is a critical parameter in various studies. It is particularly relevant for studies related to the search for suitable sites for a large project implementation requiring massive equipment to keep transportation costs, as well as construction and maintenance costs, low [24,39,40,56].
Electrical substations: Close proximity to electrical substations minimizes wire costs, prevents power losses, and simplifies installation and maintenance processes [57]. Figure 6e represents electrical substations in the research area.

MCDM Using an AHP Approach
Suitable site selection for implementation of a sensitive project such as a wind farm is always difficult since it requires a combination of various parameters and criteria defining the project location. Therefore, decision-making solutions to overcome these obstacles have been developed by integrating all of these determining criteria. Generally, the multi-criteria decision-making (MCDM) approach is always applied to address problems with many stakeholders, criteria, and objectives [58]. Moreover, this approach has been widely applied in various fields, including the energy sector to plan renewable energy projects [59][60][61][62]. The analytic hierarchy process (AHP) is a well-known MCDA approach that was initially proposed by [63] and has subsequently been greatly improved.
The methods for weighting the criteria in the MCDA are diverse. Some of these methods include AHP, fuzzy measures [64], Analytic Network Process (ANP) [65], Swara [66], entropy [67], Dematel [68], and standard deviation [69]. Although these methods are quite limited, AHP is one of the most essential and widely used methods in MCDA. The AHP method is similar to Swara's in that the expert's opinion specifies the importance and prioritization of alternatives. As for the entropy method, there are two different views of this method. According to some studies, entropy is reliable and effective [70]. However, from another point of view, entropy results do not always take into account the importance of the indices [71]. Dematel is similar to the Swara method, except that the Dematel approach is used to solve extremely complex problems. In the Dematel decision process, the expert opinion is used to develop the pairwise comparison matrix, and it has three main characteristics. The attributes are compensatory and independent of each other. Qualitative attributes are transformed into quantitative attributes [66]. The Swara and Dematel methods have been widely used in MCDA problems, especially in the renewable energy sector [72][73][74]. In this study, the AHP method was employed to address site selection problems for several reasons: It is commonly used for its ease of design and implementation. It is highly compatible with GIS, which is widely used for planning and spatial analysis of site selection problems. The consideration of the consistency and inconsistency of alternatives is one of the main advantages of this method [60].
AHP can be combined with other methods of multicriteria analysis, genetic algorithms, neural networks, etc. [60]. It also takes into account quantitative and qualitative criteria to interpret the problem [75].
The AHP method can apply various sensitivity analyses to the criteria. AHP facilitates the decision-making process, using pairwise comparison between criteria [60]. For site selection problems, in which the main objective is to select the best locations, simple approaches such as AHP are satisfactory, and more complex approaches such as Fuzzy-AHP do not necessarily lead to distinct results [76]. AHP is a structured decision support method that is primarily focused on sophisticated computations with matrix algebra [77,78]. Through this approach, a decomposition of a complicated decision-making issue into a top-down hierarchical structure can be carried out in most cases. In recent years, as geographic information system technologies have improved, GIS integration with MCDA approaches has become increasingly popular. This integration is adaptable and suited to the qualitative and quantitative investigation of multi-criteria issues with a geographical component.
In this research, we developed the decision process required for the usage of the AHP approach. This approach is provided in four steps, each of which requires clear problem identification or study's objective.
Step 1: Deconstruct the decision-making problem and explain its main characteristics or components (criteria, sub-criteria, options, etc.). Then, using a limited number of levels, create a linear hierarchy of concerns ( Figure 8). Each level has a set number of selection criteria. The aim is expressed at the most fundamental level. Subsequently, the second and third layers comprise the criterion and sub-criteria. The bottom of the hierarchy is allotted to alternatives. popular. This integration is adaptable and suited to the qualitative and quantitative investigation of multi-criteria issues with a geographical component.
In this research, we developed the decision process required for the usage of the AHP approach. This approach is provided in four steps, each of which requires clear problem identification or study's objective.
Step 1: Deconstruct the decision-making problem and explain its main characteristics or components (criteria, sub-criteria, options, etc.). Then, using a limited number of levels, create a linear hierarchy of concerns ( Figure 8). Each level has a set number of selection criteria. The aim is expressed at the most fundamental level. Subsequently, the second and third layers comprise the criterion and sub-criteria. The bottom of the hierarchy is allotted to alternatives. Step 2: Design the judgment matrix and pairwise comparison matrices for each criterion. Based on the Saaty scale (Table 2), the pairwise comparisons are grouped into a matrix using the following criteria: .. 1 : In relation to the comparisons of two criteria C1 and C2, we designate an important value of the evaluation element "a". We place the "a" value in the cell column "i" and line "j" of an important criterion. Then, we need to place the value ratio "1/a" in the cell considered less important of the comparison. C1, C2, and Cn are the comparison criteria in row "i" and column "j", which correspond to the comparison values Ci and Cj. The entries aij are often taken from the ratio scale (1/9-9) [79]. The matrix's element semantic description is provided in Table 3.  Step 2: Design the judgment matrix and pairwise comparison matrices for each criterion. Based on the Saaty scale (Table 2), the pairwise comparisons are grouped into a matrix using the following criteria: In relation to the comparisons of two criteria C1 and C2, we designate an important value of the evaluation element "a". We place the "a" value in the cell column "i" and line "j" of an important criterion. Then, we need to place the value ratio "1/a" in the cell considered less important of the comparison. C1, C2, and Cn are the comparison criteria in row "i" and column "j", which correspond to the comparison values Ci and Cj. The entries aij are often taken from the ratio scale (1/9-9) [79]. The matrix's element semantic description is provided in Table 3.  Using the evaluations provided in the previous step, each hierarchy element's relative relevance was determined. Furthermore, the eigenvector problem is addressed to establish each matrix's element priority.
First, compute the sum of each jth column value as follows: Subsequently, a normalized comparison matrix n × n aij* is generated, in which each aij in the matrix is divided by the sum of its jth column, as expressed in Equation (3): The weights' ith criterion is then computed as follows: Step 3: The individual criteria weights are calculated using the eigenvalue procedure's pairwise comparison matrices. The eigenvalue λmax is calculated by multiplying each column value by the criteria weight as follows: Then, using the following equation, we determine the weighted sum value Sw by adding the sum of each preceding matrix's rows ai: Eventually, for each row, the ratio between the weighted value sum Sw and the weighting criterion is calculated as follows: By averaging the ratio i we obtain the highest eigenvalue max.
Step 4: Calculate the consistency ratio CR (Equation (8)). The final criteria weights are validated using this ratio. Discrepancies in the comparison matrix are identified at this stage: where the consistency index CI is calculated as follows: The value of the RI varies with the size of the matrix. Table 4 shows the RI values according to the number of criteria chosen. The RC should be lower than 10% to determine that the pairwise comparison evaluations are consistent. If this is not the case, the matrix should be updated, and the element values re-evaluated.
The weighted findings (Table 5) indicate the most important parameters in wind farm development. The wind's existence is the greatest driver of wind energy, with a relative weight of 38% (Table 4), followed in second place by protected regions with a relative weight of 26%. However, distances to power plants and road networks have a significant influence, with relative weights of 13% each, while slope and elevation have the lowest relative weights of 5% and 3%, respectively. Despite their low weights, these criteria should be considered in all wind projects to avoid or minimize potential negative impacts. This weighting choice is based on our study area's good knowledge. After completing all of the AHP calculation processes, the following step is to normalize the criteria ( Table 6). The vector data are then converted to a raster format, and the matrices are reclassified into two groups (adequate: code 1; and inadequate: code 0).

Weighted Superposition
The weighted overlay tool is one of the most frequently used methods for solving multi-criteria problems, such as site selection and suitability models. For instance, users can use this functionality to combine several spatial layers with varied weights to produce a final result. Each raster layer is assigned a weight in the suitability analysis. The raster layer values are re-ranked on a scale (two classes in our case). In this study, the weighted overlay analysis was utilized to identify the most suitable and appropriate sites for future wind farm siting based on the AHP-derived weights assigned to each evaluated parameter. According to Equation (10), all selected criteria in raster format that have been reclassified to equal size (number of columns equal to the number of rows) (Figure 9) are combined into a single raster layer ( Figure 10). Weighted overlay is defined as follows: where Wi is the weight of a specific choice criterion, Ri is the criterion's matrix layer, and n is the number of decision criteria. In total, 962,612 pixels define the wind turbine installation sites. The research area has a total size of 9319 km 2 . The detected locations cover an area of 650,725,712 m 2 , or 651 km 2 . This accounts for roughly 6.98% of the study area's surface.    In total, 962,612 pixels define the wind turbine installation sites. The research area has a total size of 9319 km 2 . The detected locations cover an area of 650,725,712 m 2 , or 651 km 2 . This accounts for roughly 6.98% of the study area's surface.

Results
AHP factor weights were computed using technical, environmental, and economic requirements for wind turbines in France. Factor weights used to evaluate appropriate sites for wind farm installation are shown in Figure 11. As can be seen in Figure 9, wind speed is the most important factor, with a weight of 38%. It is followed by the respect for buffer zones around protected areas (urban areas, wetlands, biodiversity parks, etc.) with a weight of 26%, and the proximity of electrical substations and the road network with a weight of 13% each. Slope and elevation are ranked last, with weights of 5% and 4%, respectively. It should be noted that the eigenvalue max (max = 6.26) is calculated after computing criteria's weights. CI and CR values are 0.05 and 0.04, respectively. The CR value is 10%, suggesting that the research was satisfactory.

Results
AHP factor weights were computed using technical, environmental, and economic requirements for wind turbines in France. Factor weights used to evaluate appropriate sites for wind farm installation are shown in Figure 11. As can be seen in Figure 9, wind speed is the most important factor, with a weight of 38%. It is followed by the respect for buffer zones around protected areas (urban areas, wetlands, biodiversity parks, etc.) with a weight of 26%, and the proximity of electrical substations and the road network with a weight of 13% each. Slope and elevation are ranked last, with weights of 5% and 4%, respectively. It should be noted that the eigenvalue max (max = 6.26) is calculated after computing criteria's weights. CI and CR values are 0.05 and 0.04, respectively. The CR value is 10%, suggesting that the research was satisfactory.
Land 2022, 11, x FOR PEER REVIEW 18 of 29 Figure 11. Decision criteria priority weights for selecting suitable sites for future wind projects. Figure 12 depicts the appropriate location distribution for planning future onshore wind farms in various research regions' departments. Figure 10 only shows locations greater in size than 400 hectares and the road network and the electricity substation's locations. Furthermore, calculating the eligible site's surface shows that 74.62%, or 35,127.92 hectares, is in the Var department, which controls more than 80% of the study area's surface. Only 10%, or 4962 hectares, of Alpes-de-Haute-Provence department is suitable for future wind farm development, compared to 13.70%, or 6449 hectares, in Bouches-de-  in size than 400 hectares and the road network and the electricity substation's locations. Furthermore, calculating the eligible site's surface shows that 74.62%, or 35,127.92 hectares, is in the Var department, which controls more than 80% of the study area's surface. Only 10%, or 4962 hectares, of Alpes-de-Haute-Provence department is suitable for future wind farm development, compared to 13.70%, or 6449 hectares, in Bouches-de-Rhones. In the Alpes-Maritimes, however, 1.45%, or 535 hectares, is protected ( Figure 13). Figure 11. Decision criteria priority weights for selecting suitable sites for future wind projects. Figure 12 depicts the appropriate location distribution for planning future onshore wind farms in various research regions' departments. Figure 10 only shows locations greater in size than 400 hectares and the road network and the electricity substation's locations. Furthermore, calculating the eligible site's surface shows that 74.62%, or 35,127.92 hectares, is in the Var department, which controls more than 80% of the study area's surface. Only 10%, or 4962 hectares, of Alpes-de-Haute-Provence department is suitable for future wind farm development, compared to 13.70%, or 6449 hectares, in Bouches-de-Rhones. In the Alpes-Maritimes, however, 1.45%, or 535 hectares, is protected ( Figure 13).   Figure 14 demonstrates that 1121 hectares, or 12.07% of the area suitable for future wind farm construction, have an average wind speed greater than 5 m/s, which is required for wind turbine development. Furthermore, 11.24% of the area, or 1044 hectares, is located lower than 1500 m above sea level, while just 8% is on slopes less than 15 degrees; Figure 13. Suitable site percentage distribution (%) for future wind project implementation by department. Figure 14 demonstrates that 1121 hectares, or 12.07% of the area suitable for future wind farm construction, have an average wind speed greater than 5 m/s, which is required for wind turbine development. Furthermore, 11.24% of the area, or 1044 hectares, is located lower than 1500 m above sea level, while just 8% is on slopes less than 15 degrees; 5.92% is next to roadways, 2.34% is near electrical substations, and 7.91%, or 735 hectares, is outside of protected areas. Figure 13. Suitable site percentage distribution (%) for future wind project implementation by department. Figure 14 demonstrates that 1121 hectares, or 12.07% of the area suitable for future wind farm construction, have an average wind speed greater than 5 m/s, which is required for wind turbine development. Furthermore, 11.24% of the area, or 1044 hectares, is located lower than 1500 m above sea level, while just 8% is on slopes less than 15 degrees; 5.92% is next to roadways, 2.34% is near electrical substations, and 7.91%, or 735 hectares, is outside of protected areas. Figure 14. Representation of the percentage regions in the "appropriate = score 1" class for each criterion for potential wind farm locations.
On Google Earth imagery, the selected appropriate site locations for wind farm development were projected ( Figure 15). Four locations were recommended, and their selection was based on their unique characteristics (area, location, elevation, slope, wind speed, accessibility, closeness to electrical substations, etc.) as well as their proximity to populous regions while respecting buffer zones relative to protected areas. Onshore wind turbines in France typically have a power range of 1.8 to 3 MW, with rotor diameters ranging from 80 to 110 m and total heights ranging from 80 to 155 m. In fact, a 2 MW wind turbine generates 4200 MWh per year, which is roughly equivalent to the average electricity consumption of over 800 French households [8]. France is classified by the International Electrotechnical Commission (IEC) as having strong winds with high average turbulence intensity. Some wind turbine types that are easily useable in the French market have been Figure 14. Representation of the percentage regions in the "appropriate = score 1" class for each criterion for potential wind farm locations.
On Google Earth imagery, the selected appropriate site locations for wind farm development were projected ( Figure 15). Four locations were recommended, and their selection was based on their unique characteristics (area, location, elevation, slope, wind speed, accessibility, closeness to electrical substations, etc.) as well as their proximity to populous regions while respecting buffer zones relative to protected areas. Onshore wind turbines in France typically have a power range of 1.8 to 3 MW, with rotor diameters ranging from 80 to 110 m and total heights ranging from 80 to 155 m. In fact, a 2 MW wind turbine generates 4200 MWh per year, which is roughly equivalent to the average electricity consumption of over 800 French households [8]. France is classified by the International Electrotechnical Commission (IEC) as having strong winds with high average turbulence intensity. Some wind turbine types that are easily useable in the French market have been chosen in accordance with IEC design criteria. Table 7 contains detailed information about the wind turbine types and their attributes. chosen in accordance with IEC design criteria. Table 7 contains detailed information about the wind turbine types and their attributes. Figure 15. Example of the most suitable sites selected and their characteristics (geographical coordinates of their centroids, area in hectares, perimeter in km, average wind speed, and average altitude) for the development of onshore wind farms. Theoretical wind power potential may be evaluated using Equation (11) based on wind turbine output capacity, rotor diameter, and total area of appropriate land [81][82][83]. Theoretical wind power potential may be evaluated using Equation (11) based on wind turbine output capacity, rotor diameter, and total area of appropriate land [81][82][83].
TWPP is theoretical wind power potential (MW), TA is the total area of the four appropriate locations (km 2 ) (Figure 13), and AF is the area factor (MW/km 2 ). Our computations were performed on wind turbines that were situated 7d × 5d apart, where d is the rotor diameter.
Based on the theoretical wind potential calculations, the four proposed sites for future wind turbine installations may generate between 182.60 and 280.20 MW of electricity. This energy benefit will suit the study region's population demands in terms of power consumption, heating, and transportation.

Discussion
In 2021, the wind sector in France grew in relevance, accounting for 7% of the country's net power generation. Furthermore, wind power now accounts for 7.7% of total consumption [32]. More than half of France's wind farms are concentrated in two regions: Hauts-de-France and Grand Est (Figure 16), with an almost complete absence in the country's southeast.

Discussion
In 2021, the wind sector in France grew in relevance, accounting for 7% of the country's net power generation. Furthermore, wind power now accounts for 7.7% of total consumption [32]. More than half of France's wind farms are concentrated in two regions: Hauts-de-France and Grand Est (Figure 16), with an almost complete absence in the country's southeast. Wind power production (wind farms), according to the International Energy Agency, is very unevenly distributed among areas. Despite the study region possessing the country's most populated cities, it has the lowest energy productivity (201 GWh in 2021). However, the study region (Provence-Alpes-Côte d'Azur, region code: 13) ranks second to last Wind power production (wind farms), according to the International Energy Agency, is very unevenly distributed among areas. Despite the study region possessing the country's most populated cities, it has the lowest energy productivity (201 GWh in 2021). However, the study region (Provence-Alpes-Côte d'Azur, region code: 13) ranks second to last in terms of wind energy generation [60], with 77 wind turbines and a very low installed capacity of 99 MW, compared to demand ( Figure 17). in terms of wind energy generation [60], with 77 wind turbines and a very low installed capacity of 99 MW, compared to demand ( Figure 17). The study area's onshore potential is greatly limited by certain constraints, including the requirement to avoid exclusion zones imposed by environmental protection areas, historical perimeters, and the requirement to build more than 500 m from homes, as well as habitat dispersion, which reduces the percentage of territory eligible for wind power. This research was conducted for all of these reasons. In addition, it may assist various governmental agencies, policymakers, researchers, and investors in planning and developing wind energy projects in this difficult location. Regarding criteria, the Ministers of Ecological Transition, Territorial Cohesion, and Energy Transition previously investigated a set of documents and reports on technical requirements, regulations, and environmental and urban planning issues related to wind turbine development in France, which served as the foundation for our criteria and con-  The study area's onshore potential is greatly limited by certain constraints, including the requirement to avoid exclusion zones imposed by environmental protection areas, historical perimeters, and the requirement to build more than 500 m from homes, as well as habitat dispersion, which reduces the percentage of territory eligible for wind power. This research was conducted for all of these reasons. In addition, it may assist various governmental agencies, policymakers, researchers, and investors in planning and developing wind energy projects in this difficult location.
Regarding criteria, the Ministers of Ecological Transition, Territorial Cohesion, and Energy Transition previously investigated a set of documents and reports on technical requirements, regulations, and environmental and urban planning issues related to wind turbine development in France, which served as the foundation for our criteria and constraints. The thresholds, on the other hand, must be closely tied to certain location characteristics. As a result, the criteria suggested in this study were used with considerable caution. They are currently being explored by wind planning professionals. Practical experience in the subject field is also advantageous for the assessment of visual findings [30].
In addition to the French regulations concerning the determining criteria for land use planning of future wind farm projects, a consultation of confidential reports and a discussion with experts and former researchers in the field of wind farm planning was carried out. The values for each criterion were selected according to French legislation. As no studies have been published on this topic in France, we also based our selection on research undertaken in Europe (e.g., Greece) as mentioned above. Other means can be used to define the important criteria, such as filling in questionnaires by experts in the field. Interviews with experts could also be an effective solution for the determination of criteria.
Furthermore, restriction criteria for onshore wind farm planning in France are well defined; only one scenario is required in this case: average wind speed greater than 5 m/s, altitude less than 1500 m, slope less than 15%, proximity to roads (2 km) and electrical substations (1.5 km), and at least 2 km from protected areas.
Unsurprisingly, as shown in Figure 18a, there is a substantial correlation between areas of high average wind speed and selected site locations as suitable for wind farms. Indeed, the project limits were previously chosen based on a wind speed map in France, which is a region where the average yearly wind speed (at least 50 m above ground level) exceeds 5 m/s. These locations are all accessible by national, regional, or occasionally freeway roads, and the majority of them are near electricity substations (Figure 18d). Furthermore, 30% of the locations are at high altitudes (over 1500 m) (Figure 18c). This is due to the high wind speed in these high-altitude areas, as well as high weight assigned to the wind parameter, and the low weight assigned to the elevation parameter. The selected sites are more than 2 km away from the protected regions (Figure 18e). Since the Southern Alps surround the northern and northwestern parts of the research region, various sites are on steep slopes. Nonetheless, we were able to choose really good locations on moderate slopes ( Figure 15). Other criteria, such as acceptance of these installations by populations and associations; administrative procedures and their validation by local authorities; energy demand in these territories; and pre-existing installation replacement, can all have an impact on decisions to install wind turbines. Indeed, the criteria required for wind farm construction, namely minimum average wind speed required, minimum acreage required, closeness of highways and electrical substations, and distance from protected areas, are all gathered at the identified sites. Some research employed the same choice criteria with various limitation levels based on the area and state restrictions. For example, [30] used the same criteria but with different limitation values because the study was conducted in Spain. Nonetheless, three situations were investigated, each with distinct limitation values and weights for the criterion. This scenario-based method is useful when there is little information or rigid wind-level limitations. However, [61] conducted a study to identify viable places for developing onshore wind projects in a rural zone using only four characteristics (urban area or habitat, vegetation, slope, and wind speed). Furthermore, they investigated three scenarios. The multi-criteria analysis technique used, however, is fuzzy rather than AHP. In summary, the multi-criteria technique used in the research varied, as did the number of scenarios examined, the number of criteria addressed, and weights assigned to the criteria. According to the findings of all studies on decision making for project implementation or multi-factor problem solving using GIS-based multi-criteria analysis approaches [62,84,85], when no restriction values are well defined by the state or agreed upon by experts in the field in question, it is recommended that several alternatives be implemented. The number and total area of suitable sites vary by region. For example, in our case, only 7% of the total area of the study area is suitable for the development of future wind projects. Even though it is a difficult location with various environmental, topographical, and urban restrictions, research and planning for wind farms in this area with great wind potential are still on hold. However, based on the theoretical wind potential calculation equation (Equation (10)), each of the sites depicted in Figure 18 produces an average of 45 MW, which is more than adequate to fulfill the demands of the local population. The number and total area of suitable sites vary by region. For example, in our case, only 7% of the total area of the study area is suitable for the development of future wind projects. Even though it is a difficult location with various environmental, topographical, and urban restrictions, research and planning for wind farms in this area with great wind potential are still on hold. However, based on the theoretical wind potential calculation equation (Equation (10)), each of the sites depicted in Figure 18 produces an average of 45 MW, which is more than adequate to fulfill the demands of the local population.

Conclusions
This research aimed to offer a method for identifying potential sites for future wind energy projects based on geographic information systems (GIS) and multi-criteria decisionmaking (AHP), as well as to contribute to the literature on renewable energy planning. To the best of the author's knowledge, this is the first research of its kind in France. Thus, this research was carried out in a region of France's southeast that has high wind energy potential but is also the least productive. Another reason for choosing this region was to overcome the numerous constraints that limited the region's energy output. For wind farm siting, six criteria were adopted, practically addressing in full the economic, technical, and socio-environmental challenges associated with these facilities and uses. Most of the criteria were based on worldwide literature, in addition to French wind turbine legislation.
According to the findings, many sites were identified as suitable for wind farms. Visual and manual analyses were performed on these sites to choose those with an area larger than 400 ha, a high average wind speed, accessibility by roads, proximity to electrical substations, and a distance from protected areas. Four sites with an average installed capacity of 45 MW were selected and must be confirmed by the appropriate state authorities. The decision tool provided in this article may be utilized in any part of the world by adapting it to the specific characteristics of each territory, as well as the distinct needs and policies.
Although the results presented in this paper are specific to France, the methodology presented provides an interesting reference model that can be transposed and adapted with relative ease. This assumes that the different constraints and criteria are adapted to the specific needs of energy planners and to the particularities of each study area.
Despite the quality and reliability of the IGN database used, the methodology followed (including the choice of the MCDA method), and the analysis performed in this work, certain aspects are to be recommended for future projects, such as taking into account the knowledge of the study area and the regulations put in place by the government concerned, and the relevant choice of determination (decision) criteria and their weightings, which often vary between experts' opinions and from one country to another.
The spatial resolution of the data used is also an important element, especially for topographic data (DEM) and wind data.
In addition, a validation of the identified suitable sites could lead to a more robust and real interpretation of the results, either through aerial photography (drone) or a field visit to these sites.
Future studies could consider extending the proposed method to investigate the theoretical energy potential of wind generation in order to benefit from their complementarity and overcome the inherent intermittency of renewable energies. This study can also be the starting point for a project to install wind turbines or solar panels in the study area by simply replacing the wind variable with the temperature variable. Furthermore, this study can also contribute to the creation of new investments and, consequently, new jobs in the region.