Criteria-Based Fuzzy Logic Risk Analysis of Wind Farms Operation in Cold Climate Regions

: Different risks are associated with the operation and maintenance of wind farms in cold climate regions, mainly due to the harsh weather conditions that wind farms experience in that region such as the (i) increased stoppage rate of wind turbines due to harsh weather conditions, (ii) limited accessibility to wind farms due to snow cover on roads, and (iii) cold stress to workers at wind farms. In addition, there are risks that are caused by wind farms during their operation, which impact the surrounding environment and community such as the (iv) risk of ice throw from wind turbines, (v) environmental risks caused by the wind farms, and (vi) social opposition risk to installing wind farms in cold climate regions, such as the Arctic. The analysis of these six risks provides an overall view of the potential risks encountered by designers, operators, and decision makers at wind farms. This paper presents a methodology to quantify the aforementioned risks using fuzzy logic method. At ﬁrst, two criteria were established for the probability and the consequences of each risk; with the use of experts’ judgments, membership functions were graphed to reﬂect the two established criteria, which represented the input to the risk analysis process. Furthermore, membership functions were created for the risk levels, which represented the output. To test the proposed methodology, a wind farm in Arctic Norway was selected as a case study to quantify its risks. Experts provided their assessments of the probability and consequences of each risk on a scale from 0–10, depending on the description of the wind farm provided to them. Risk levels were calculated using MATLAB fuzzy logic toolbox and ranked accordingly. Limited accessibility to the wind farm was ranked as the highest risk, while the social opposition to the wind farm was ranked as the lowest. In addition, to demonstrate the effects of the Arctic operating conditions on performance and safety of the wind farm, the same methodology was applied to a wind farm located in a non-cold-climate region, which showed that the risks ranked differently.


Introduction
Wind energy applications in cold climate regions (CCRs) have gained more attention recently, and are growing at a rapid rate of approximately 20% per year according to the Global Wind Energy Council [1]. There are various encouraging reasons for installing wind farms (WFs) in CCRs: among others, the fact that the density of air in such regions is 10% higher than other non-cold climate regions, which results in higher availability of wind power resources [2]. In addition, abundant air resources exist in mountainous CCRs such as the Arctic region of Norway, where relatively steady winds with higher velocities [3]. Moreover, CCRs are known to be less and sparsely inhabited than other regions in the same country. Therefore, having WFs installed there will have less perceived impact on people than WFs built in large cities for example, and will likely encounter less opposition to installing WFs from residents. The paper is organized as follows: Section 2 presents the methodology followed to analyze the risks. Section 3 explains the fuzzy logic process. Section 4 identifies six risks to WFs in CCRs. Section 5 defines the criteria for the risks, considering five levels of probabilities (very low, low, medium, high, very high) and four for the severity of consequences (low moderate, high, very high). Section 6 considers a WF in Arctic Norway as a case study to demonstrate the proposed methodology and ranks the six risks; finally, conclusions are presented in Section 7.

Methodology
The methodology adopted in this work, shown in Figure 1, starts with identifying the potential risks usually WFs in CCRs are subjected to. The risks-relevant literature and research are being reviewed to define criteria for the inputs to the risk analysis, which are the probabilities of occurrence and the severity of consequences of the identified risks. Afterwards, the defined criteria are sent to a selected group of experts who will provide estimated values (between 0 and 10) for the different levels of the probabilities and consequences, which represent the input to the risk analysis and for the risks' levels, which represents the output. Based on the data collected from the experts, and by using MATLAB fuzzy logic toolbox, membership functions are graphed to represent the levels of the inputs and the output. The paper is organized as follows: Section 2 presents the methodology followed to analyze the risks. Section 3 explains the fuzzy logic process. Section 4 identifies six risks to WFs in CCRs. Section 5 defines the criteria for the risks, considering five levels of probabilities (very low, low, medium, high, very high) and four for the severity of consequences (low moderate, high, very high). Section 6 considers a WF in Arctic Norway as a case study to demonstrate the proposed methodology and ranks the six risks; finally, conclusions are presented in Section 7.

Methodology
The methodology adopted in this work, shown in Figure 1, starts with identifying the potential risks usually WFs in CCRs are subjected to. The risks-relevant literature and research are being reviewed to define criteria for the inputs to the risk analysis, which are the probabilities of occurrence and the severity of consequences of the identified risks. Afterwards, the defined criteria are sent to a selected group of experts who will provide estimated values (between 0 and 10) for the different levels of the probabilities and consequences, which represent the input to the risk analysis and for the risks' levels, which represents the output. Based on the data collected from the experts, and by using MATLAB fuzzy logic toolbox, membership functions are graphed to represent the levels of the inputs and the output.  This methodology is implemented to a WF in Arctic Norway as a case study for demonstration. Initially, data regarding the WF are gathered from the WF operator, research articles, and weather stations that publish their data online. Afterwards, the collected data conceded to a selected group of experts, who are asked to quantitatively evaluate the corresponding probability of occurrence and severity of consequences of each of the identified risks on a scale from 0 to 10. Experts' judgments were then fed into MATLAB fuzzy logic toolbox, which calculated the risk level, using specifically defined rules. Eventually, the risks are ranked depending on the resulting risk level, the highest risk level was assigned a rank of (1) and the lowest risk was assigned a rank of (6).

Fuzzy Logic Process
Fuzzy logic is based on fuzzy set theory developed by Zadeh [25]. Fuzzy sets are a generalization of the classical set theory, indicating that the classical set theory is a special case of the fuzzy set theory [26]. Fuzzy logic is an efficient tool in risk assessment as it compensates for the lack of knowledge and vagueness encountered when assessing the risks related to complex technological systems, and can be very helpful when dealing with fuzzy linguistic terms such as low, medium, high, etc., to describe the risks, as these terms do not have sharp boundaries for their definitions and can hold different interpretations when interpreted by different experts [23].
For example, taking X as a universal set that contains all objects used in the risk analysis process. Probability, consequences, and risk levels are the input and output variables used in the risk analysis. Each one of these variables contains a number of objects (fuzzy terms) that were previously defined in X as follows: X = {very low, low, medium, moderate, high, very high, moderate-high, extremely high} Input variable (probability) = {very low, low, medium, high, very high} Input variable (consequences) = {low, moderate, high, very high} Output variable (risk level) = {very low, low, moderate, moderate-high, high, very high, and extremely high} Each fuzzy term in the universal set X is described as a fuzzy subset (A), characterized by a membership function µ(x), which assigns to each object a degree of membership that has values between zero (no-membership) and one (complete membership). Based on that, a fuzzy subset A can be written as a set of pair: A = {(x,µ A (x)); x∈X}, where x is a numbering value provided by the experts to describe the input variable (i.e., the probability or the consequences) [23].
The fuzzy logic process followed in this paper is based on Mamdani method [27], which is the most commonly used method in fuzzy logic, using the center of gravity method to calculate the output value of the risk level during the defuzzification step, unlike the Sugeno method, which uses the weighted average method to calculate the risk level [28]. Figure 2 shows the three main steps (fuzzification, fuzzy logic inference, defuzzification) followed in applying the fuzzy logic process to calculate the risk level and ranking the risks [19] Probability of a risk is Low and the Consequences are High, Then the Risk level is Moderate.

•
Defuzzification: This is a counter step to the fuzzification step, where the resulted fuzzy risk levels are converted, using MATLAB fuzzy logic toolbox, into numbers reflecting how high or low the risk level is, where higher number reflects higher risk level and vice versa. Following this step, the risks to WFs can be ranked.

Risks Identification
The risks to WFs in CCRs are mainly caused by the harsh weather conditions they experience, which degrade their overall performance and resilience. In addition, WFs impose risks on their surroundings, these risks can be associated with environmental concerns and social opposition to WFs being installed in CCRs, such as the Arctic region. The following description of the six risks will be used in Section 5 to define the criteria for the probabilities and consequences of the risks.
Risk1. Increased WT stoppages due to harsh weather conditions (WT stoppage): Low temperatures affect the physical properties of materials and the normal operation of electronic devices [29], leading to increased failure rate in WT components. Lubrication oil viscosity, for example, changes under low temperatures and has an impact on the dimensions and mechanical properties of different components in the WT. This results in possible overheating and higher fatigue charges on components, with one of the most affected being the gearbox, as its lifetime is considerably reduced [30].
Moreover, ice accretion on WTs, which mostly takes place on the leading edge of the blades, increases mass and aerodynamic imbalances, and might render the operation of the WT unsafe, leading to the shutting down of the WT and the loss of power production until the ice is removed or melted [31]. The probability P of stoppages of WTs, due to failures and ice accretion, can be calculated using the Poisson distribution in Equation (1) [32]: (1) where λ denotes the rate of stoppages per WT per a specific period t (e.g., a month or a year), and k is the number of WT stoppages the distribution calculates the probability of. Risk 2. Cold stress to workers (Cold stress): Cold temperatures cause cold stress to crew workers, and limit their dexterity [33]. Serious cold-related illnesses and injuries, caused by trench foot, frostbite, and hypothermia, may occur in case of extreme cold temperatures, in addition to permanent tissue damage, and death that may result as a consequence to major cold-related injuries.
High wind speeds and cold temperatures are the two main factors contributing to cold stress to workers [34]. Wind Chill Temperature (WCT) is a measure that determines the likelihood that workers are subjected to the risk of frostbite, which can be calculated Fuzzification •Input from experts •Converting the numbering values for Prob. and Cons. into fuzzy terms.

Fuzzy logic inference
•Establishing the rules for combining probabilities and consequences.

Figure 2.
Overview of the fuzzy logic process.

Risks Identification
The risks to WFs in CCRs are mainly caused by the harsh weather conditions they experience, which degrade their overall performance and resilience. In addition, WFs impose risks on their surroundings, these risks can be associated with environmental concerns and social opposition to WFs being installed in CCRs, such as the Arctic region. The following description of the six risks will be used in Section 5 to define the criteria for the probabilities and consequences of the risks.
Risk 1. Increased WT stoppages due to harsh weather conditions (WT stoppage): Low temperatures affect the physical properties of materials and the normal operation of electronic devices [29], leading to increased failure rate in WT components. Lubrication oil viscosity, for example, changes under low temperatures and has an impact on the dimensions and mechanical properties of different components in the WT. This results in possible overheating and higher fatigue charges on components, with one of the most affected being the gearbox, as its lifetime is considerably reduced [30].
Moreover, ice accretion on WTs, which mostly takes place on the leading edge of the blades, increases mass and aerodynamic imbalances, and might render the operation of the WT unsafe, leading to the shutting down of the WT and the loss of power production until the ice is removed or melted [31]. The probability P of stoppages of WTs, due to failures and ice accretion, can be calculated using the Poisson distribution in Equation (1) [32]: where λ denotes the rate of stoppages per WT per a specific period t (e.g., a month or a year), and k is the number of WT stoppages the distribution calculates the probability of. Risk 2. Cold stress to workers (Cold stress): Cold temperatures cause cold stress to crew workers, and limit their dexterity [33]. Serious cold-related illnesses and injuries, caused by trench foot, frostbite, and hypothermia, may occur in case of extreme cold temperatures, in addition to permanent tissue damage, and death that may result as a consequence to major cold-related injuries.
High wind speeds and cold temperatures are the two main factors contributing to cold stress to workers [34]. Wind Chill Temperature (WCT) is a measure that determines the likelihood that workers are subjected to the risk of frostbite, which can be calculated using Equation (2), where V is the wind speed (km/h) 10 m above the surface and T is the air temperature ( • C) [34]: WCT[ • C] = 13.12 + 0.621T − 11.37V 0.16 + 0.3965TV 0. 16 (2)  Table 1 has been generated using Equation (2). The table is used to determine whether the workers at WFs in CCRs are subjected to the risk of frostbite or not, where the shaded region indicates an increasing risk of frostbite [34]. Risk 3. Limited accessibility to wind farms due to snow cover: CCRs are known for their diverse landscape, especially the Arctic region, including tundra, glaciers and steep mountains. The Arctic terrain can be challenging to move around. Snow accumulation on WFs roads and pathways will reduce accessibility to the WTs, this is especially significant when it comes to the maintainability of WTs, which might be reduced under such conditions [4]. The severity of snow accumulation determines whether the WF needs to implement snow removal strategies, using special snow removal vehicles, or if it would be enough to use specially equipped vehicles to access the WTs, such as snow mobiles and snowcats, in case normal vehicles were not useable.
Risk 4. Thrown ice pieces from operational wind turbines (Ice throw): This phenomenon can occur when pieces of ice are either thrown away from an operational WT, see Figure 3, due to the aerodynamic and centrifugal forces, or dropped down if the WT was idle. In both cases, ice pieces landing on the ground will represent a hazard to the safety of the WF including the WTs, WF facilities, crew personnel, and animals [8]. A simple equation (Equation (3)) can be used to measure the distance of thrown ice pieces from an operational WT as follows [35]: where (d) is the throwing distance, (D) is the rotor blade diameter, and (H) is the hub height. The probability of ice throw from operational WTs depends on the probability of ice formation on the WT blades, the probability of being thrown away to a location where ice pieces may represent a hazard to WF surroundings, and the probability of members of the public, crew personnel, and animals being present within the range of landing ice pieces [37]. The probability of ice accretion on WT blades depends on many factors such The probability of ice throw from operational WTs depends on the probability of ice formation on the WT blades, the probability of being thrown away to a location where ice pieces may represent a hazard to WF surroundings, and the probability of members of the public, crew personnel, and animals being present within the range of landing ice pieces [37]. The probability of ice accretion on WT blades depends on many factors such as the air temperature, wind speed, liquid water content (LWC), median volume diameter (MVD), and the elevation of WT from the sea level [38].
Rime ice and glaze ice are the two most common types of ice to accrete on the blades of WTs. Rime ice forms when supercooled water droplets freeze immediately upon impacting the surface of the blade, while glaze ice forms when the liquid water freezes shortly after impacting the surface of the blade [39]. Glaze ice accretion forms near the freezing point (0 • C) and has strong adhesion to the surface; it is transparent and has a higher density than rime ice. On the other hand, rime ice has lower adhesion to the surface and has a white or opaque color [40].
The probability of ice formation on WTs can be reflected by estimating the daily intensity of the icing events in kg/m 2 . Table 2 summarizes five site icing index categories that can be used to determine the intensity of icing on WTs in the WF location on a daily basis [41]. Risk 5. Environmental risks: The Arctic, as an example of CCR, is known for its sensitive environment. Locating WFs in the Arctic will lead to possible impacts on the critical habitats and threatened species. For example, bird mortalities caused by WTs have been debatable [42] for several years, even though it is stated that wind energy killed 20 times fewer birds compared to fossil fuels, and the number of birds killed by WTs can be negligible compared to some other human activities [43]. In addition, the construction phase of wind farms might result in pollution of nearby surface or underground water [44]. Similarly, the use of liquids such as the gearbox lubricating oil might result in pollution in case it leaked from the WT. Moreover, the Arctic area is known for reindeer grazing, therefore, having WFs built on winter grazing areas for reindeer might affect their density, especially during the construction phase or even after it. Risk 6. Social opposition: The visual presence of WTs can be annoying, especially to residents living nearby WFs. The presence of WFs might stop or limit the ability of local communities to utilize the surrounding lands and might affect its economy [45]. In addition, the generated noise by WTs might be annoying to residents living nearby WFs. The sources of the generated noise by WTs are the mechanical components such as the gearbox and control mechanisms, and the rotation of the WT blades through the air. Noise levels are measured by decibels (dB(A)), which is a scale designed to measure how the human ear perceives the sound frequencies. The day-evening-night noise level (L den ) is a European standard to express the noise levels from machines throughout an entire day [46]. The institution of occupational safety and health (IOSH) designed a scale for classifying noise levels [47], which can be used to assess the severity of noise emitted by WTs, as shown in Table 3.   Table 4 determines the criteria for estimating the probability of each of the identified risks. The criteria are based on reviewed research studies, measured data, and human evidence. Selecting the probability level for each risk type is primarily dependent on the WT or WF under the study. Table 4. Criteria for the probabilities of risks experienced by WFs in CCRs.

Risk Very Low (Vl) Low Medium High Very High (Vh)
Increased WT stoppage rate [48] The probability of stoppage using Equation (1) is between 0-20% The probability of stoppage using Equation (1) is between 21-40% The probability of stoppage using Equation (1) is between 41-60% The probability of stoppage using Equation (1) is between 61-80% The probability of stoppage using Equation (1) is between 81-100% Cold stress [34] Mild wind chill conditions. The wind chill temperature can be larger or equal to Limited accessibility [4] No snow cover on the roads. The WF is easily accessible.
The roads of the WF are covered with snow but is still accessible with normal cars.
Accessing the WF requires the use of snowcats and snow mobiles due to snow cover.
There is a need to remove the snow off the road using special vehicles and equipment such as snowplows, blowers, loaders, and deicer trucks, etc.
The accessibility is very low due to extreme weather conditions and excessive snow cover on the roads.
Ice throw [41] The site icing index according to Table 2 is S5, indicating occasional icing. No roads, residential areas, or facilities are in the range of thrown ice pieces.
The site icing index according to Table 2 is S4, indicating light icing. Most roads residential areas, and facilities are not in the range of thrown ice pieces.
The site icing index according to Table 2 is S3, indicating moderate icing. Roads and facilities in the surroundings are in the range of thrown ice pieces.
The site icing index according to Table 2 is S2, indicating strong icing.
The probability of being hit by ice pieces is high.
Excessive ice accretion on the WT blades, S1. the main road is very close to the WF site; therefore, surroundings are in great danger of being struck by ice pieces thrown from the WTs.

Environmental risks
The WF is not built on a migration route for birds and is not built on winter grazing area for reindeer. No records of water or environmental pollution by the WF exist.
The WF is built on a migration route for birds and on a winter grazing area for reindeer, but the effects are not significant. No records of water or environmental pollution by the WF exist.
The WF is built on a migration route for birds and on a winter grazing area for reindeer and affect their existence. No records of water or environmental pollution by the WF exist.
The WF is built on a migration route for birds and on a winter grazing area for reindeer and affect their existence significantly high. There is a record of water and environmental pollution by the WF.
The WF affects the existence of migrating birds and reindeer density in the area very significantly, with significant water and environmental pollution record by the WF.
Social Opposition [45,49] The WF is located far from residential areas, does not have an impact on the livelihood of local communities, and the noise level is very low, L den = 0-40 dB(A).
The WF is located far from residential areas, does not have an impact on the livelihood of local communities, and the noise level is low, L den = 41-70 dB(A).
The WF is located near residential areas, with bearable effects on the livelihood of local communities, and the noise level is moderate, L den = 71-100 dB(A).
The WF is located near residential areas, with high effects on the livelihood of local communities, and the noise level is high, L den = 101-140 dB(A).
The WF is located close to residential areas, with very high effects on the livelihood of local communities, and the noise level is very high, L den > 140 dB(A). Table 5 shows the criteria selected for measuring the consequences of the identified six risks. The consequences of risks are defined differently based on the type of risk being assessed. The consequences can be evaluated depending on the resulting system deterioration, injuries or loss of lives, maintenance delays, and short-or long-term effects.  The WT stoppage did not cause deterioration in the WF operation and was slightly noticed by the operator.
The WT stoppage caused slight deterioration in the WF performance and was highly noticeable by the operator.
The WT stoppage was caused by a failure that significantly deteriorated the WF performance or led to minor injuries to humans nearby.
The WT stoppage would seriously affect the ability of the WF to continue operating, or cause damage, serious injury or death.
Cold stress [50] No injury or illness. Minor injury or minor occupational illness.
Medium injury or medium occupational illness.
Serious injury or death of humans.
Limited accessibility [51] No delay in carrying out maintenance activities to the failed WTs.
Maintenance is slightly delayed, with slight loss of power production Maintenance is significantly delayed, with significant loss of power production.
Maintenance of the failed WT is highly delayed, with so highly increased power losses.
Ice throw [50] No injury or illness. Minor injury or minor occupational illness.
Medium injury or medium occupational illness.
Serious injury or death of humans.
Environmental risks [50] Minor environmental damage, readily repaired and/or might incur slight costs to correct and/or in penalties.
Short-term environmental damage, with slight costs to correct and/or in penalties.
Medium-term environmental damage, with significant costs to correct and/or in penalties.
Long-term environmental damage, with very high costs to correct and/or in penalties.

Social
Opposition [50] The WF has minor effects on the touristic activities in the area. The WF noise levels do not cause hearing impairments.
The WF has short-term effects on the touristic activities in the area. The WF noise levels cause minor hearing impairments.
The WF has medium-term effects on the touristic activities in the area. The WF noise levels cause severe hearing impairments.
The WF has long-term effects on the touristic activities in the area. The WF noise levels might cause permanent hearing loss.

Experts' Judgments
The preceding criteria were sent to seven experts, who were asked to provide their range of values for each fuzzy linguistic term used to assess the probability (Very low (Vl), Low, Medium, high, Very high (Vh)) and the consequences (Low, Moderate, High, Very high) term, as well as the risk levels, which are described as Very low (Vl), Low (L), Moderate (M), Moderate-high(MH), High(H), Very high (Vh), and Extremely high (Eh). Afterwards, the data gathered from the experts were used to design the membership functions that reflected these fuzzy terms.
The selected experts have backgrounds ranging from university professors working on wind energy to experienced staff at wind farms in Arctic Norway. Experts were assumed to have equal weights for their answers. Based on the average values for the probabilities, consequences, and risk levels collected from the experts, and by using MATLAB fuzzy logic toolbox, the triangular membership functions for the input variables (i.e., probabilities and consequences) and the output variable (the risk levels) were defined, as shown in Figures 4 and 5, where the x-axis represents the input values provided by the experts (from 0 to 10), and the y-axis represents the degree of membership (from 0 to 1) for each membership function. The combination between these three variables can be represented by a 3-dimensional surface plot, in Figure 6, which shows a fuzzy risk matrix. illness.
Environmental risks [50] Minor environmental damage, readily repaired and/or might incur slight costs to correct and/or in penalties.
Short-term environmental damage, with slight costs to correct and/or in penalties.
Medium-term environmental damage, with significant costs to correct and/or in penalties.
Long-term environmental damage, with very high costs to correct and/or in penalties.
Social Opposition [50] The WF has minor effects on the touristic activities in the area. The WF noise levels do not cause hearing impairments.
The WF has short-term effects on the touristic activities in the area. The WF noise levels cause minor hearing impairments.
The WF has medium-term effects on the touristic activities in the area. The WF noise levels cause severe hearing impairments.
The WF has long-term effects on the touristic activities in the area. The WF noise levels might cause permanent hearing loss.

Experts' Judgments
The preceding criteria were sent to seven experts, who were asked to provide their range of values for each fuzzy linguistic term used to assess the probability (Very low (Vl), Low, Medium, high, Very high (Vh)) and the consequences (Low, Moderate, High, Very high) term, as well as the risk levels, which are described as Very low (Vl), Low (L), Moderate (M), Moderate-high(MH), High(H), Very high (Vh), and Extremely high (Eh). Afterwards, the data gathered from the experts were used to design the membership functions that reflected these fuzzy terms.
The selected experts have backgrounds ranging from university professors working on wind energy to experienced staff at wind farms in Arctic Norway. Experts were assumed to have equal weights for their answers. Based on the average values for the probabilities, consequences, and risk levels collected from the experts, and by using MATLAB fuzzy logic toolbox, the triangular membership functions for the input variables (i.e., probabilities and consequences) and the output variable (the risk levels) were defined, as shown in Figures 4 and 5, where the x-axis represents the input values provided by the experts (from 0 to 10), and the y-axis represents the degree of membership (from 0 to 1) for each membership function. The combination between these three variables can be represented by a 3-dimensional surface plot, in Figure 6, which shows a fuzzy risk matrix.     For more illustration, the y-axis in Figure 4a,b represents the degree to which a certain input (a probability or a consequence) can be described as low, medium, high, etc. For example, if the average value of the probability of a certain risk, given by experts, was 5, that means that the degree of membership of the probability of that risk is 80% medium, as per Figure 4a. Similarly, if the severity of a specific risk was determined by experts to be 9, that means that the severity of that risk is 100% very high, as per Figure 4b. Following that, the defined inference rules will determine the level of the risk in fuzzy terms. A logical inference rule that applies to this example can be: if the probability is medium and severity is very high, then risk level is high. The degree of membership of the risk level is determined using the minimum operator as in Equation (4) [20]: μ (Risk is high) = min (0.8; 1) = 0.8 (4) where μ (Risk is high) is the degree of membership of the risk level as a high risk. Afterwards, the risk level is determined by referring to Figure 5, where the x-axis value for the risk level that corresponds to 80% high risk is equal to 7.5. The fuzzy inference functions in MATLAB used in this risk analysis application are shown in Table 6. The membership functions defined by the experts are used to generate the fuzzy rules that will be used to rank the risks. A total of 5×4 = 20 rules were generated. Examples of these rules are as follows: • Rule 1: If probability is very low and consequence is low, then the risk level is very low. • Rule 11: If probability is medium and consequence is high, then the risk level is moderate-high. • Rule 19: If probability is very high and consequence is high, then the risk level is very high. A detailed description of the WF was sent to the experts, who were asked to provide their numeric values for the probabilities and consequences of each risk. Afterwards, the For more illustration, the y-axis in Figure 4a,b represents the degree to which a certain input (a probability or a consequence) can be described as low, medium, high, etc. For example, if the average value of the probability of a certain risk, given by experts, was 5, that means that the degree of membership of the probability of that risk is 80% medium, as per Figure 4a. Similarly, if the severity of a specific risk was determined by experts to be 9, that means that the severity of that risk is 100% very high, as per Figure 4b. Following that, the defined inference rules will determine the level of the risk in fuzzy terms. A logical inference rule that applies to this example can be: if the probability is medium and severity is very high, then risk level is high. The degree of membership of the risk level is determined using the minimum operator as in Equation (4) [20]: µ (Risk is high) = min (0.8; 1) = 0.8 (4) where µ (Risk is high) is the degree of membership of the risk level as a high risk. Afterwards, the risk level is determined by referring to Figure 5, where the x-axis value for the risk level that corresponds to 80% high risk is equal to 7.5. The fuzzy inference functions in MATLAB used in this risk analysis application are shown in Table 6. The membership functions defined by the experts are used to generate the fuzzy rules that will be used to rank the risks. A total of 5×4 = 20 rules were generated. Examples of these rules are as follows:

•
Rule 1: If probability is very low and consequence is low, then the risk level is very low. • Rule 11: If probability is medium and consequence is high, then the risk level is moderate-high. • Rule 19: If probability is very high and consequence is high, then the risk level is very high. A detailed description of the WF was sent to the experts, who were asked to provide their numeric values for the probabilities and consequences of each risk. Afterwards, the fuzzy logic toolbox in MATLAB was used to calculate the risk level using the centroid method [52]. Figure 7 illustrates the different probability and consequence combinations, based on the defined 20 fuzzy rules, to calculate the risk level values.

Case Study: Wind farm in Arctic Norway
A WF in Arctic Norway, with a layout shown in Figure 8, was selected to strate the proposed methodology. The selected WF is located in a valley at aroun above the sea level. The WF consists of fourteen 2.3 MW WTs. In order to sup experts in their estimations, operational and site specifications data about the W acquired from different sources. For example, data regarding failure rates, icin maintenance reports, and WTs performance were acquired from the WF opera data covered two years of WF operation, from 2019 to 2020. In addition, data re the WF's site specifications and weather data were acquired from published rese ticles and manufacturers' technical publications, as well as weather stations that their data online.

Analysis
MATLAB fuzzy logic toolbox was utilized to calculate the level of each ris basis of the average values of the probabilities and consequences of each risk prov the experts. Table 7 shows the average probability and consequence values det by the experts as well as the resulting risk level calculated by MATLAB. In addi

Case Study: Wind farm in Arctic Norway
A WF in Arctic Norway, with a layout shown in Figure 8, was selected to demonstrate the proposed methodology. The selected WF is located in a valley at around 420 m above the sea level. The WF consists of fourteen 2.3 MW WTs. In order to support the experts in their estimations, operational and site specifications data about the WF were acquired from different sources. For example, data regarding failure rates, icing rates, maintenance reports, and WTs performance were acquired from the WF operator; this data covered two years of WF operation, from 2019 to 2020. In addition, data regarding the WF's site specifications and weather data were acquired from published research articles and manufacturers' technical publications, as well as weather stations that publish their data online.

Case Study: Wind farm in Arctic Norway
A WF in Arctic Norway, with a layout shown in Figure 8, was strate the proposed methodology. The selected WF is located in a val above the sea level. The WF consists of fourteen 2.3 MW WTs. In o experts in their estimations, operational and site specifications data acquired from different sources. For example, data regarding failur maintenance reports, and WTs performance were acquired from th data covered two years of WF operation, from 2019 to 2020. In addit the WF's site specifications and weather data were acquired from pu ticles and manufacturers' technical publications, as well as weather s their data online.

Analysis
MATLAB fuzzy logic toolbox was utilized to calculate the level basis of the average values of the probabilities and consequences of ea the experts. Table 7 shows the average probability and consequence the experts. Table 7 shows the average probability and consequence values determined by the experts as well as the resulting risk level calculated by MATLAB. In addition, the table shows the ranking of the risks, where limited accessibility to WF risk is assigned the highest rank (1), and social opposition risk is assigned the lowest rank (6). Risk 1 (WT stoppage). From the data gathered from the WF operator, there were 1993 stoppages experienced by the WTs during 2019, mainly due to maintenance. In addition, December 2019 was the month that witnessed the highest rate of WT stoppages due to icing, which was 65 stoppages. Equation (1) was applied separately to determine the probability of stoppage per WT per month due to failure and per month due to icing. It was found out that the probability of stoppage was 29% per WT per month. Those stoppages did not result in deterioration of the WF production and operation was resumed as regular, according to the WF operator.
Risk 2 (Cold stress). The coldest average ambient temperatures impacting the WF occur during January, February, and March, according to Figure 9a, with the coldest average temperature of as low as nearly −12 • C recorded in February. In addition, the average monthly wind speed during the same month is 13.6 km/h, as shown in Figure 9b [53]. By applying Equation (2), the WCT during February is calculated to be about −19 • C, which indicates no risk of frostbite to workers at the WF. In addition, according to the WF operator, there are no injuries recorded or illnesses caused to workers at the WF by very cold WCTs. Risk1 (WT stoppage). From the data gathered from the WF operator, there were 1993 stoppages experienced by the WTs during 2019, mainly due to maintenance. In addition, December 2019 was the month that witnessed the highest rate of WT stoppages due to icing, which was 65 stoppages. Equation (1) was applied separately to determine the probability of stoppage per WT per month due to failure and per month due to icing. It was found out that the probability of stoppage was 29% per WT per month. Those stoppages did not result in deterioration of the WF production and operation was resumed as regular, according to the WF operator.
Risk 2 (Cold stress). The coldest average ambient temperatures impacting the WF occur during January, February, and March, according to Figure 9a, with the coldest average temperature of as low as nearly −12 °C recorded in February. In addition, the average monthly wind speed during the same month is 13.6 km/h, as shown in Figure 9b [53]. By applying Equation (2), the WCT during February is calculated to be about −19 °C, which indicates no risk of frostbite to workers at the WF. In addition, according to the WF operator, there are no injuries recorded or illnesses caused to workers at the WF by very cold WCTs.  [54]. Data were gathered from Weather Atlas.

Risk 3 (Limited accessibility). The Norwegian Meteorological Institute [55]
provides information about the snow depth of specific sites in Norway through weather stations distributed in different areas of the country. The closest weather station to the WF is located in Straumsnes. Figure 10 shows that the highest recorded snow depth in 2020 occurred during February until May, with a maximum snow depth of 75 cm recorded in April. Such accumulation of snow requires specially equipped vehicles, such as snowcats and snowmobiles to be available all the time, to maintain access to the WF under severe weather conditions, which was confirmed by the WF operator during a visit to the WF by the authors. Besides that, during periods of peak snow accumulation, snow removal strategies are implemented, which can be costly [4]. Average wind speed (km/h) (b) Figure 9. Monthly average temperatures (a) and wind speeds (b) at the WF's location [54]. Data were gathered from Weather Atlas.

Risk 3 (Limited accessibility). The Norwegian Meteorological Institute
[55] provides information about the snow depth of specific sites in Norway through weather stations distributed in different areas of the country. The closest weather station to the WF is located in Straumsnes. Figure 10 shows that the highest recorded snow depth in 2020 occurred during February until May, with a maximum snow depth of 75 cm recorded in April. Such accumulation of snow requires specially equipped vehicles, such as snowcats and snowmobiles to be available all the time, to maintain access to the WF under severe weather conditions, which was confirmed by the WF operator during a visit to the WF by the authors. Besides that, during periods of peak snow accumulation, snow removal strategies are implemented, which can be costly [4].

Risk 4 (Ice throw).
A field study shows that the daily intensity of icing in Na city close to the WF, is approximately 14.5-18.5 kg/m 2 during January, Februa March, as shown in Figure 11 [56]; those measurements can be assumed close eno be applicable to the site under consideration. This indicates a light icing intensi site icing index for the WF site, according to Table 1. By applying Equation (3) to ca the throwing distance of ice pieces, it is found that the throwing distance is 255 closest residential area to the WF is kilometers away. Moreover, a main road pass to the WF, but it is around 500 m away from the closest WT, which means that n dents or personnel exist within the range of ice throw risk.

Risk 5 (Environmental risks).
A study by Jacobsen et.al. [57] revealed that the n of registered birds migrating through the WF area was low, where the average num

Risk 4 (Ice throw).
A field study shows that the daily intensity of icing in Narvik, a city close to the WF, is approximately 14.5-18.5 kg/m 2 during January, February, and March, as shown in Figure 11 [56]; those measurements can be assumed close enough to be applicable to the site under consideration. This indicates a light icing intensity with site icing index for the WF site, according to Table 1. By applying Equation (3) to calculate the throwing distance of ice pieces, it is found that the throwing distance is 255 m. The closest residential area to the WF is kilometers away. Moreover, a main road passes next to the WF, but it is around 500 m away from the closest WT, which means that no residents or personnel exist within the range of ice throw risk.

Risk 4 (Ice throw).
A field study shows that the daily intensity of icing in Narvik, a city close to the WF, is approximately 14.5-18.5 kg/m 2 during January, February, and March, as shown in Figure 11 [56]; those measurements can be assumed close enough to be applicable to the site under consideration. This indicates a light icing intensity with site icing index for the WF site, according to Table 1. By applying Equation (3) to calculate the throwing distance of ice pieces, it is found that the throwing distance is 255 m. The closest residential area to the WF is kilometers away. Moreover, a main road passes next to the WF, but it is around 500 m away from the closest WT, which means that no residents or personnel exist within the range of ice throw risk.

Risk 5 (Environmental risks).
A study by Jacobsen et.al. [57] revealed that the number of registered birds migrating through the WF area was low, where the average number of observed birds passing through the WF was around five birds during 4 h of observation in a day. This is primarily due to the topography and local conditions of that area, which is not inviting to migrating birds. However, the area around the WF has great value for reindeer husbandry and is used for reindeer grazing during autumn and winter. In this respect, field results show that the density of reindeer in WF areas in Arctic Norway do not change significantly, and that reindeer can adapt and keep using WF areas for grazing [58]. With regards to environmental pollution, there is no evidence that the WF Risk 5 (Environmental risks). A study by Jacobsen et.al. [57] revealed that the number of registered birds migrating through the WF area was low, where the average number of observed birds passing through the WF was around five birds during 4 h of observation in a day. This is primarily due to the topography and local conditions of that area, which is not inviting to migrating birds. However, the area around the WF has great value for reindeer husbandry and is used for reindeer grazing during autumn and winter. In this respect, field results show that the density of reindeer in WF areas in Arctic Norway do not change significantly, and that reindeer can adapt and keep using WF areas for grazing [58]. With regards to environmental pollution, there is no evidence that the WF caused pollution to nearby waters or to the environment in general.
Risk 6 (Social opposition). The presence of the WTs and the noise they generate are not expected to cause annoyance to humans due to the fact that the WF is located away from residential areas, in a valley with surrounding mountains as high as 2000 m above sea level from the north and south. Besides, the location of the WF is not classified as a touristic area. However, there is a main road that passes next to the WF, which means that drivers and passengers passing by the WF will be subjected to this noise level for very short time. According to one study, the level of noise generated by a similar WT located at approximately 500 m far, which is nearly the distance between the closest WT and the main road, would be 48 dB(A) at a wind speed of 10 m/s [59], meaning that the noise severity can be described as low.

A Wind Farm under Normal Operating Conditions
In order to compare the effects of cold climate operating conditions on the calculated risk level and ranking of risks, the same methodology is applied to a WF located in a non-cold-climate region. The Kozbeyli WF [45], in Turkey has higher reliability, with lower rate of WT stoppages, no ice accretion on the blades of the WTs, less snow accumulation on the roads of the WF, and a 3.1 km access road built to guarantee access to the WF. However, the WF is close to an Environmentally Protected Area in Foça, with a bird migration route 4 km to the south of the WF [60]. Moreover, in the WF area, there are endangered species such as Passer Domesticus and Crocidura Russula, that were identified and listed by the Bern Convention [45].
Furthermore, the WF is located 1.3 km near to a touristic village, which has natural and historical values. Social acceptance of the WF was determined to be poor due to its impact on tourism in that area. Based on the preceding information, Table 8 demonstrates the ranking of risks using experts' judgments and MATLAB fuzzy logic toolbox. It can be seen from the table that the social opposition risk ranked the highest, followed by the environmental risks, while the risk of ice throw from WTs ranked as the lowest risk, flowed by the risk of WT stoppage.

Conclusions
In this paper, six types of risks that are related to the operation of WFs in CCRs. These risks were distinguished as being caused by the harsh weather conditions, and risks caused by the WFs on their surrounding environment and community. The identified risks were analyzed using expert judgment and MATLAB fuzzy logic toolbox. The identified risks are the increased WT stoppages risk, cold stress to workers risk, limited accessibility to WFs risk due to snow accumulations on the roads, ice throw from WTs risk, environmental risks, and social opposition risk.
Based on a research review, gathered data, and published data, and experts' reasoning, two criteria tables were defined for the probability of occurrence and the severity of consequences of each risk. Furthermore, experts' judgments and MATLAB fuzzy logic toolbox were used to graph the membership functions for the probabilities and consequences of each risk (the inputs), as well as risk levels (the output). The risk levels were calculated based on a set of 20 rules generated using the experts' data.
A WF in Arctic Norway was selected to illustrate the proposed methodology. Experts were provided with a description of the WF and were asked to deliver their assessed values for the probabilities and consequences of each risk. Through the calculated risk level, it was concluded that limited accessibility to the WF ranked as the highest risk, followed by WT stoppage. On the other hand, social opposition was ranked as the lowest risk followed by the ice throw.
In order to demonstrate the methodology further, a WF that is not subjected to cold climate operating conditions, located in Turkey, was selected. The social opposition to the WF was ranked as the highest risk followed by the environmental risks, where the ice throw risk and WT stoppage risk were deemed to be the lowest-ranked risks. This was due to the fact that the Turkish WF was installed close to a village with touristic value and in a area that is a home for endangered species.
Author Contributions: A.M.M. developed the methodology, conducted the literature review, and gathered the data to establish the probabilities and consequences criteria, gathered the data from the experts, performed the risk analysis using MATLAB fuzzy logic toolbox, and wrote the paper; A.B., as a main supervisor, followed up all study steps and gave helpful guidance and reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.

Data Availability Statement:
The data gathered from the wind farm are subjected to a non-disclosure agreement. Therefore, an approval from the wind farm is required in case the data are needed. The data gathered from the experts can be provided upon request from the editors.

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