Dust impact on concentrated solar power: A review
Article information
Abstract
Many sites with high solar radiation face high dust loads that reduce energy generation by concentrated solar power plants. This review presents the attenuative impacts of atmospheric aerosols, as well as reflectivity losses due to soiling of solar reflectors, by covering both experimental investigations and numerical studies; along with presenting the theoretical background. The chemical nature of aerosols, and the physics of soiling and atmospheric extinction phenomena (scattering and absorption) are also reviewed. Suspended particles like aerosols result in atmospheric extinction of the solar radiation that reaches the concentrators, and the deposition of these particles on the solar reflectors provokes decreases up to 80% in their reflectivity, and thus enhances the cumulus of optical losses and the reduction of energy production. Even though dust affects both CSP and photovoltaics, CSP technologies suffer more losses. The impact of dust should be particularly considered during the planning phase of solar thermal plants, since its consequent reduction in energy output can be severe. While there have been multiple papers to review dust-related problems for PV, the present paper is the first literature review dedicated to the impact of soiling on concentrated solar power.
1. Introduction
During the last century, human activities resulted in a drastic increase in the amount of the atmospheric carbon dioxide [1–3]. To prevent global warming from rising above 1.5°C, 100% of the energy must be produced by transition to zero-emission in the upcoming years. Most of the energy should be generated using renewable energy [4–9]. Among the resources of renewable energy, solar energy is the most available one, and a feasible option for thermal energy applications [10, 11].
Concentrated Solar Power (CSP) technologies are very promising among the renewable energies and are growing fast, as they provide solutions to the decade’s most alarming issues of global warming combined with shortage of energy, water and food [12]. The fuel of the solar thermal energy is the solar irradiance that penetrates the atmosphere and is concentrated by CSP plants. However, the presence of water vapor and suspended atmospheric particles (aerosols) attenuates the solar radiation, and affects both PV and CSP technologies [13–17]. Also, the deposition of these particles on the collecting surfaces of the solar field provokes decreases in the reflectivity of solar mirrors and reduces the efficiency of thermal solar systems [18–20]. Therefore, being the common cause of the optical losses (extinction of irradiance and loss of reflectivity), aerosols are the most important factor that influences the exploitation of solar thermal energy (under clear-sky conditions). Not considering aerosols and dust will cause a significant uncertainty in the estimation of electricity production, and in the cost of large scale CSP projects [19, 21]. Moreover, since DNI is the fuel of CSP plants and small reflectivity losses result in serious drops in the plants’ overall efficiency, losses of produced energy due to soiling are more important for CSP projects than other solar energy applications [21, 22].
Forecasting the presence of aerosols and the consequent attenuation will help in reducing extra-costs and optimizing the performance of CSP plants [23, 24]. For the solar tower system, losses due to attenuation can reach over 10% [24]. The optical losses are enhanced further by the soiling of the solar reflectors [25].
Dust problems persist in the absence of better solutions, making research into aerosols-related issues critical for solar energy harvesting initiatives, particularly locations where there is no natural cleaning due to a lack of rain events. The aim of this paper is to review the studies that evaluated the effect of aerosols (suspended particles) and dust (deposited particles) on the performance of CSP technologies. The paper is organized as follows: in Section 2, previous informative literature reviews are presented, while Section 3, the theoretical background, explains the physics of the different phenomena and the chemical nature of aerosols. In Section 4, recent research articles that treated the impact of dust on the exploitation of solar thermal energy are analyzed and reviewed, followed, in Section 5, by the review of the wide literature of the effect of different climatic parameters on the performance of CSP plants. The main effects are -as already mentioned- the attenuation of the solar irradiance and the drop in the solar reflector’s optical efficiency. The range of the losses varies widely and depends on the sites surrounding CSP plants. In Section 6, the different mitigation approaches for soiling problems are presented and analyzed. The literature review is finalized with a summary (Section 7) and a conclusion on the losses caused by suspended and deposited dust particles. No previous studies reviewed the impact of dust on concentrated solar power, even though many papers reviewed the same issue for PV, which makes the present work the first literature review of this issue for CSP. Also, no previous review is available considering both of the optical losses, even though the extinction of solar irradiance and the soiling of solar reflectors are both caused by particles present in the vicinity of CSP plants. Since the phenomena are of common origin and occurrence, their negative impacts shouldn’t be separated.
2. Literature Reviews of Dust Impact on CSP
Previous informative papers have reviewed the literature of the issue at hand; some were interested in the attenuative effect of the suspended aerosols, while others focused on the impact of soiling on mirrors. Hereafter, we present the focuses and conclusions of these reviews:
Sarver et al. [21] presented a comprehensive state of the art covering studies of dust impact on solar energy (PV and CSP) between 1942 and 2013 and different measures to reduce the losses. The authors concluded that even though energy yield is reduced by dust for both PV and CSP, CSP is more affected since small reflectivity losses would result in major thermal losses (loss range: CSP, 20–50% compared to PV, 15–30%).
Ghazi et al. [26] reviewed the soiling of flat surfaces and noted that despite the high solar intensity available in many regions of the world, these regions host high aerosols loads, especially the MENA region. It was reported that solar devices should be inclined at an angle to the horizontal, no matter how small the angle is, in order to avoid the maximum soiling characterizing the horizontal position. As mitigation, it was recommended that CSP mirrors should be set to zero pitch angle early in the morning or late in the evening, so that the dust would slip off.
Costa et al. [27] authored a complete review of the literature of soiling-related solar energy issues for the period between 2012 and 2015, and updated it for 2016, but only few research papers were concerned by CSP soiling [28]. They highlighted the growing number of studies concerning the issue and presented a reference-listing table of the papers during the period of interest, many of the papers were concerned with PV soiling, compared to the fewer number of papers studying CSP soiling.
Hanrieder et al. [29] reviewed the approaches of estimating the atmospheric extinction in solar tower plants, inputs were obtained from on-site measures or satellite data. It was reported that power losses due to atmospheric extinction by scattering and absorption are considerable and highly dependent on the geometry and size of the heliostats field, the receiver, and the local atmospheric conditions (especially sites with high aerosols or water vapor presence). The authors recommended plant optimization design using time series of the atmospheric extinction (based on in-situ measurements) instead of the usual models.
Bouaddi et al. [30] reviewed the methods of cleaning soiled CSP reflectors like water-based methods, ultrasonic cleaning, dry cleaning methods, etc. Dry cleaning methods include coatings, electrodynamic screens (EDS), shaking, vibrating and applying high voltage to the reflectors. It was reported that ultrasonic cleaning uses 7 times less water than pressurized water cleaning.
Antonanzas-Torres et al. [31] authored a review of seventy clear-sky irradiance (CSI) models. The input data were obtained from the Baseline Surface Radiation Network (BSRN) and Aerosols Robotic Network (AERONET). Only 8 models got Root Mean Square Error (RMSE) less than 5%, 7.5% and 10% for Global Horizontal Irradiance (GHI), Diffused Irradiance (DI), and Direct Normal Irradiance (DNI) resp; and thus, were considered best models in Cabauw (Netherlands) and Carpentras (France).
In the present paper, we review the literature of dust impact on CSP. Even though this issue was reviewed for PV (like the recent review paper by Kazem et al. [11]), no previous literature review was dedicated to the consequences of soiling on CSP plants, thus the present up-to-date paper. Also, no review (and few research papers like [18]) linked the impact of atmospheric extinction to soiling problems despite the proved relationship between the two phenomena. Finally, we highlight in our paper insights and recommendations about research necessities regarding the issue at hand.
3. Theory
3.1. Chemical Nature of Aerosols
An aerosol is generally defined as a suspension of liquid or solid particles in a gas, with particle diameters in the range of 10−9 to 10−4 m. In scientific research, the term Particulate Matter (PM) refers to particles with aerodynamic diameters ≤ 1 μm (PM1), ≤ 2 μm (PM2) or ≤ 10 μm (PM10). In the troposphere (lower atmosphere), particles’ mass and number concentrations vary from 1 to 100 μg/m3 and from 102 to 105 particle/cm3 resp [32, 33]. Atmospheric aerosols differ widely in composition as a result of the diverse natural and anthropogenic origins: Carbonaceous Aerosols (CA) are constituted of Black Carbon (BC) and Organic Carbon (OC). OC is emitted from biomass burning (eg: forest fires) and transported biological materials (eg: wastes of animals & plants). BC is the inorganic black carbon (also known as elemental carbon), its main source is fossil fuels’ combustion [32]. Similarly, sulfur is a consequence -and an indicator- of hydrocarbon combustion and industrial activity, it is the result of long-range transport processes and continuous conversion to secondary particulate sulfur [32, 34].
Dust and sea salt are common species of atmospheric aerosols that are considered in simulation models (eg: [35]), dust presence in the atmosphere attenuates the solar energy and reduces the efficiency of CSP plants [14, 20], the main components of desert dust are quartz and calcite rich in silicon and calcium resp (eg: Egypt [34], Morocco [36]), while the existence of silicon is due to disturbances of desert soil by wind or human activity, calcium originates from industrial and/or building activities [34].
In the present paper, we are concerned with the attenuation of solar irradiance by natural aerosols (also referred to as: dust/soiling). Generally, the physical properties of the atmospheric aerosols are determined using the following techniques [32]: Differential Mobility Analysis (DMA), inertial separation, Scanning/Transmission Electron Microscopy (SEM/TEM), light scattering (Mie), Spectrophotometry, Photoacoustic Spectroscopy, and Nephelometry for absorption and scattering coefficients.
3.2. Physics of Aerosols’ Behavior
3.2.1. Extinction of solar irradiance (scattering and absorption)
When the solar energy penetrates Earth’s atmosphere, a part of the incident energy is lost by scattering and absorption, especially when the sky is turbid (Fig. 1) [37]. Scattering of irradiance is primarily due to aerosols [38], and happens at all wavelengths, a part of the scattered radiation is sent back to the outer space while the other part, the diffuse radiation, is scattered towards the ground. The solar radiation -unaffected by the scattering phenomenon- that reaches directly from the solar disk is the direct (or beam) normal irradiance DNI [37, 39]. Absorption happens only at selective wavelengths, the short-wave UV radiation below 0.29 μm is absorbed by atmospheric Ozone, while CA Carbon dioxide absorbs long-wave Infra-Red radiation [37, 39, 40]. Aerosols’ optical effect is much more scattering than absorbent, according to Pettit and Freese [41].
3.2.2. Soiling of solar reflectors
Dust is lifted off the ground; by the wind, and suspended in the air when the momentum of aerodynamic drag and lift forces is more important than the momentum of gravity and inter-particle forces [20], the vertical diffusion of aerosols in the air is caused mainly by the atmospheric turbulence, Brownian motion and gravity (see Fig. 2) [35]. When gravity becomes important, dust particles fall on concentrators of the solar field [42], another situation during which important soiling occurs is the windy (hazy) weather, not because higher wind speeds entrain more dust into the air, rather, the increase of soiling is because stronger wind exposes the reflectors to air flow with greater concentration of aerosols per unit of time [43]. Only dry dust deposition is considered since wet deposition happens during precipitations, which have a predominant cleaning effect [44].
Dust deposition is described by its most representative physical quantity: dust deposition flux Fd
where Cd is the dust concentration in the air (μg/m3), νd is the deposition velocity of particles and αt is the tilt angle of the solar collector [44], Fd represents the amount of dust that vertically falls on the reflector. νd is the sum of νg and νt, νg is the velocity of aerosols falling freely in the fluid (air in this case) under the equilibrium of gravity and dragging forces, νt is the velocity resulting from turbulence, diffusion and impaction, it depends mainly on the Brownian motion [44]:
In Eq. (3) ra and rb are respectively the aerodynamic and the quasi-laminar Brownian diffusion resistances. The detailed theory of dust deposition flux is given by Picotti et al. [44]
Aerosols particles are attracted and held to the mirrors’ surfaces by the adhesive forces: Van der Waals, electrostatic and capillary forces [42, 44]; even though the surface is tilted, the particles are not detached because the aforementioned adhesive forces are thousands of times stronger than the gravity component which is parallel to the surface, see Fig. 3 (the perpendicular component contributes to the adhesion) [45].
The removal of aerosols is generally due to natural washing by precipitations [36, 44], or to the action of the drag forces increased by increasing wind speed which, furthermore, generates turbulence eddies that contribute to the removal as well [43]. Humidity, however, prevents the removal of soiling through capillary forces, above certain values of relative humidity RH (70% to 80%) it’s very difficult for the removal by wind to happen [43, 46, 47].
4. Recent Research and Advances in Dust Impact on CSP
4.1. Effect of Soiling on Solar Reflectors
Since the most recent review paper reviewed the up-to-2016 literature of the impact of dust on solar energy (including thermal solar energy) [28], we review in this subsection the papers that studied dust deposition on solar reflectors and were published in 2017 or later.
To assess the impact of dust deposition on the CSP plant’s energy output, Alami Merrouni et al. [48] used Ebsilon software. Without consideration of soiling, the energy yield in Morocco and Portugal was 1.4 GWeh and 1.3 GWeh resp, but with consideration of soiling, both sites produce comparable electricity supply of 1.25 GWeh during the exposure period, this is expected given the higher soiling in Morocco compared to Portugal [48, 49]. Azouzoute et al. [50] evaluated the impact of soiling on energy production of a 5 MW plant. The cleanliness index dropped by up to 30% in just 8 days. Also, different meteorological parameters were measured for one year, the data were implanted in System Advisor Model (SAM), two weeks of soiling caused average daily production loss of 22%. While Hachicha et al. [51] reported 36% decrease in the thermal efficiency of parabolic trough collectors in UAE due to a reflectivity drop of 63% (3 months of exposure). Also, dust caused more loss in CSP than PV by a factor of 3 to 5, and was correlated to wind speed and direction. Whilst in Portugal, soiling rate was reported to be 8 to 14 times higher for CSP than for PV for an identical quantity of soiling according to Bellmann et al. [52]. Based on Mie-scattering and the various cleanliness and incidence angle measurements (Fig. 4), a model that evaluates the optical losses due to soiling for both CSP and PV was developed.
Bouaddi et al. [53] analyzed soiling’s effect on different types of CSP mirrors; glass reflectors recover perfectly their reflectivity after rain, better than aluminum ones; however, glass mirrors tend to lose their reflectivity faster. Also, clean mirrors were soiled faster than unclean ones. The highest soiling rate was recorded for glass mirrors in August: −2.3% per day, by the summer’s end, the mirrors lose up to 73% of cleanliness. A similar study was done by Wiesinger et al. [54] and the mirrors’ average reflectance losses were 36.9% for aluminum and 10.8% for glass. Azouzoute et al. [55] confirmed that aluminum mirrors were less soiled than glass mirrors, their maximum soiling rate was 25% in comparison with 35% for glass ones. Authors recommended aluminum mirrors for CSP plants, especially in arid regions for industries with low temperature needs.
The severity of dust deposition does not vary only with the solar mirror’s material type, but also with the different climates and regions hosting the CSP plants: Sansom et al. [56] run a comparative study in Algeria, Iran and Libya. Only aerosols smaller than 250 μm adhered to mirrors. Dust particles up to 1 mm can be lifted by strong wind and damage the surface, thus particles larger than 250 microns should be used in laboratory erosion simulations. Guerguer et al. [57] (Morocco, 3 years) compared soiling between seaside and desert sites: mirrors exposed in the seaside were more soiled than those in the desert. The reflectivity losses were up to 80% (August) for seaside site and 23% (March) for desert site. More than 80% of the deposited particles were smaller than 30 microns. In a similar coast/desert comparison, Endaya et al. [58] simulated -by accelerated tests- the impact of dust on a CSP reflector in Libya. The fine coastal dust caused more optical losses than the desert’s one: the reflectivity losses were 8% and 4%, respectively.
Long-range aerosols-transport phenomena can also cause soiling, Conceição et al. [49] compared soiling of mirrors in Portugal and Morocco, the daily dirtiness rates were 0.013 in Morocco and 0.004 in Portugal. Red rain in Morocco caused reflectance drop by 19.7% in one day. Similarly, Alami Merrouni et al. [48] reported high soiling rates of 10 to 21% in Morocco compared to 3 to 7% in Portugal. In a wider comparison, Wiesinger et al. [59] run 1 year of exposure in six CSP-potential sites in Spain, Morocco and Chile, as well as accelerated aging tests with an open-loop setup. An erosion potential risk matrix is proposed and validated in 5 sites, the matrix has 6 meteorological and geological parameters as inputs. Likewise, Wiesinger et al. [60] studied the erosion risk of reflectors in Morocco and Kuwait. The maximum loss of the specular reflectance of mirrors was in Kuwait: 42.5% in just 9 months (in Morocco only 5.9% after 25 months). Matal et al. [61] performed erosion simulations. The reflectivity loss is proportional to air velocity and tilt angle, the reflectivity loss due to erosion at 90° was greater than that at 45° (for 25 m/s, the optical losses were 11.03% and 5.31% resp). Recently, in 2021, Wiesinger et al. [62] investigated the effect of height and orientation on the erosion of solar mirrors. The erosion intensity decreased with increasing height and increased with the increasing tilt angle. Also, short, strong wind events caused more intense erosion than the more frequent, moderate events. Buendía-Martínez et al. [63] established a lifetime prediction model for solar reflectors, based on accelerated aging tests and outdoor exposure at 10 sites. The model considers the effect of corrosion, erosion and soiling by combination of the results of both tests.
Some studies were interested in the relationship between the deposited dust’s weight, and the consequent reflectance drop: Zhao et al. [64] studied soiling of Fresnel linear reflectors in China. After 48 days, dust density of 2.5 g/m2 caused 9.4% loss in average relative reflectivity. The soiling was maximal for the horizontal position (9.19%) and minimal for the vertical one (0.72%). A model for the cleanliness factor was developed with about 1% standard deviation. Another study in China was done by Wu et al. [65] to assess the soiling of a parabolic trough collector. After a month, the collector’s bottom edge was more soiled and lost more reflectivity (1.03 g/m2, 15%) than the top edge (0.83 g/m2, 12%). More severe losses were reported by Azouzoute et al. [66] in Morocco. The highest soiling loss rate was after two weeks: 39%, correlated with a deposition density of 0.8 g/m2. Usamentiaga et al. [67] studied the impact of soiling on the reflectivity in the visible and infrared spectra, with 6 different pollutants. The pollutants weight varied from 0.1 g to 5 g. For some pollutants, the reflectivity was higher than the clean surface, this due to the opposite optical behaviors of these pollutants between infrared and visible domains.
An efficient way to bypass the aforementioned experimental studies is by developing models -on physical bases- that evaluate the soiling of CSP reflectors, like the one validated by Picotti et al. [44] which considers: dust deposition rate, gravity, diffusion and turbulence, Brownian motion, aerodynamic drag and van der Waals forces, the model considers also dust shading and blockage effects; the model’s details can be found in ref. [44]. The inputs are: aerosols concentration and size distribution, mirrors’ position, wind speed and ambient temperature. The average relative error is 14% for the 45° tilt. A model considering incidence angles’ impact on the optical efficiency of CSP plants was developed and validated by Heimsath & Nitz [68]. It showed proportionality between incidence angles and the intensity of scattering by aerosols, and inverse proportionality to specular reflectivity. The optical efficiency of parabolic trough reflectors in Spain was found to be overestimated by 1.5%. Wu et al. [65] proposed a model to predict dust characteristics and tilt angles’ influence on the soiling rate; the model is validated, with a less than 3% standard deviation, by the measured cleanliness of a parabolic trough collector in China. Similarly, Heimsath et al. [69] proposed a model that considers the effects of variable reflectance angles and soiling rates.
In addition to natural dust, bird droppings contribute to the soiling of solar mirrors as well and result in serious reflectivity drops [64, 65], the losses in the overall performance might even reach 60% [67]. The severity of the losses due to bird droppings makes them comparable to major soiling losses such as localized dust deposition and red rain events [48, 49].
Finally, according to Vicente et al. [70] (Spain, 2 years), the application of anti-soiling coatings to the solar mirrors is of positive effect. They confirmed tilt angle and soiling to be inversely proportional.
In Table 1, we present studies that proposed interesting methods and instruments to measure the soiling rates of CSP reflectors.
Summary of Recent Studies Proposing Different Methods and Devices to Measure CSP Reflectors’ Soiling Rates
4.2. Aerosols’ Attenuation of Solar Irradiance
The optical losses in CSP plants result from the cumulus of two overlapping phenomena: the reflectivity loss due to dust deposition and the attenuation (extinction) of DNI by suspended aerosols, aerosols’ suspension results from the uplift of dust particles by wind, and dust deposition results from aerosols’ settling, thus the atmospheric turbidity should also be evaluated when studying the risk of soiling in CSP sites because reflectivity and AOD are inversely proportional as reported by Raillani et al. [18] (Morocco, 1 year) who confirmed decreasing reflectivity when AOD increases. This supports the proposition of Griffith et al. [46] who recommended the monitoring of AOD around CSP plants to assess the risk of soiling, especially for plants near industrial zones.
AOD is important for extinction assessment, especially in solar tower plants: Polo et al. [76] (PSA, 2 years) used the data of two high-resolution cameras to improve Polo extinction model [77] for central tower plants. The model was validated with AOD data of the nearby AERONET station. The improved model has 6.1% of RMSE (previously 10.9%). Similarly, Carra et al. [78] developed and validated an Extinction AOD method using extinction data at PSA. The data were used to validate AERONET and MERRA2 datasets, as well as the data of MODIS Aqua and Terra satellites. Hanrieder et al. [79] (Spain and Morocco, 35 months) also studied the transmittance in solar tower plants and validated a model [80] for 1 km slant range. Aerosols were described by three approaches: Homogenous distribution, Vertical aerosols profiles of LIVAS database [81] and The Boundary Layer Height [82]. The average broadband transmittance was 0.89 for PSA (Spain), 0.87, and 0.86 for the Moroccan sites. Statistical parameters were MBE < 5% and RMSE < 8%, which reflects the good model fitting.
Turbidity parameters (mainly AOD) are crucial for the prediction of DNI by CSI models as the models’ performance varies globally due to AOD variations. Ruiz-Arias and Gueymard [22] compared 15 frequently cited CSI models. AOD’s effect on DNI was more important than Pw’s. DNI inter-model discords are higher and broader than GHI, especially in Asia, MENA and Central Africa. Large DNI discords are more apparent for extreme solar zenith angles, these deficiencies are related to extreme solar altitudes and would not affect most solar applications. Benkciali et al. [83] (Algeria, 3 years) studied 18 CSI models under different climates. The best models were: ESRA, Dogniaux, MAC and Yang. Turbidity parameters (TL, AOD, α and β) were found to have the most important impact on the models’ accuracy. Simple models over-estimated the DNI, while others underestimated it. Grosjean et al. [84] (10 sites worldwide, 1 year) reported that ASTM spectrum [85] has a tendency of underestimating AOD. The difference between the solar spectrum calculated with ASTM and the data from 10 AERONET stations is less than 0.1%.
Satellites-derived AOD data might become dependable inputs for CSI models: Goto et al. [35] studied the data measured by two next-generation satellites. The agreements with AERONET were high to moderate, thus the data are reliable. The satellites could detect a large load of aerosols which no other device could detect. Yet satellites cannot retrieve AOD of areas with thick clouds or high surface-albedo. Bright and Gueymard [86] (452 AERONET stations, 18 years) validated the AOD data monitored by MODIS Aqua and Terra satellites. Their combination offered more precise AOD values than those given separately by Aqua or Terra. Both satellites lack the accuracy for solar energy applications, notably for a decisive parameter like AOD.
Many authors -worldwide- developed local CSI models with high to acceptable accuracy. Behar et al. [87] (Chile, 1 year) proposed and validated a method to predict solar energy and atmospheric turbidity using meteorological inputs. The respective MBE of DNI and TL are 0.91% and 1.31%, which shows the accuracy of the model despite its simplicity. In Northern Africa, Marif et al. [88] (Algeria, 31 months) studied the turbidity and developed a CSI model, they confirmed that TL and β are highly affected by wind speed, ambient temperature, and humidity. Aerosols’ loads were high in summer and low in winter. The mean TL and β were 3.63 and 0.094 resp, while the highest values were 5.07 and 0.159. Capderou’s Linke turbidity formula [89] was validated, and was judged suitable for low altitudes. In Asia, Rathore et al. [39] (India, 8 years) generated a clear-sky solar map of India; they validated the r.sun model. The GHI, DNI and DI were calculated with low MBE under the value of ±10% recommended by [90], thus the r.sun model is qualified to compute the irradiation in India. Yu et al. [91] (China, 84 sites) studied the radiative effect of aerosols and water vapor, and used datasets to simulate the radiative effect with a RMSE of 9.00%. A consistency is observed between the 15 years averaged distribution of AOD and that of aerosols’ radiative effect; the yearly radiative effect was 0.51 W/m2 (3.35% over 15 years). In Europe, Abreu et al. [92] (Portugal, 10 years) validated a CSI model with AERONET and BSRN datasets. Yang [93] (SURFRAD network, USA, 4 years) proposed a new factor to evaluate CSI models, namely: the mean square error scaling.
Bright et al. [94] (5 BSRN sites, 3 years) developed a global clear sky detection model. The model performed well under five different climatic zones and overcame the limitations observed in other models (eg: poor performance under high zenith angles). The model is freely available and coded in Matlab. Also, Bright et al. [95] provided access to reanalysis data for CSI modelling through a free python package of MERRA-2.
5. Effect of Dust Deposition and Climate Parameters on CSP Performancce
5.1. Effect of Soiling on Reflectors
To evaluate the impact of dust on solar tower plants, Singh et al. [96] investigated the soiling of heliostats and Open Volumetric Air Receivers (OVAR), the tests were done with single and double heliostat set-ups. Dust blocked 20% of the OVAR’s pores, which caused reduction of the convective heat transfer and the overall efficiency. For the single heliostat and the first one of the double-heliostat set-up, the deposition of dust was uniform, for the second heliostat, the deposition was localized due to the turbulent flow caused by the first heliostat. The effect of soiling on a parabolic trough power plant was studied by Niknia et al. [97], 44% and 60% reflectance losses were observed after 44 days and two months resp, the same losses were correlated with dust weights of 1 g/m2 and 1.5 g/m2 resp. The following empirical energy loss correlation was given as a function of the deposited dust thickness:
Factors such as the tilt angles of solar reflectors or their exposure direction have an impact on the intensity of the soiling. Alami Merrouni et al. [36] evaluated this effect for different materials of solar reflectors. The monthly cleanliness drops of the horizontal mirrors were 45% and 33% for glass and aluminum resp, while for +45° mirrors the drop was 14% for both materials, the 0° (vertical) and −45° (ground-facing) mirrors remained clean (3% drop). An important soiling factor highlighted by the authors was the wind direction. Similarly, Griffith et al. [46] studied soiling of tilted mirrors at a candidate CSP site. The soiling of the horizontal mirrors was twice the +45° ones and the mean reflectivity loss was 0.5% per day. Other similar studies on the effect of reflectors’ material, tilt and direction on the intensity of soiling can be found in table II [98–101].
Bouaddi et al. [102] modeled the cleanliness of CSP reflectors with Dynamic Linear Models, then compared it with that of exposed 45° tilted mirrors; the dust distribution was homogenous (as noted also by [96, 103]); the lowest measured cleanliness was around 20% after 3 months. The local linear trend was the best to describe the cleanliness factor. Other models are proposed in Table 2 [104, 105].
Summary of Selected Experimental and Computational Studies Treating Aerosols’ Role in Both the Reflectivity Loss of Solar Mirrors and the Attenuation of DNI
Finally, even though most of the soiling is due to the deposition of dust, many authors highlighted other sources of optical losses due to soiling, mainly bird droppings (eg: [46, 106]).
To sum up, dust deposition is non-negligible and can cause energy losses up to 88% [107] and reflectivity losses up to 80% [102], the soiling rate depends upon the wind direction and the exposure direction of the reflectors [36, 46, 108], and their placement within the solar field, in particular, first-row reflectors suffer from uniform soiling and are soiled more than the rest of the solar field [96, 102, 103]. These are important factors to consider when planning optimal cleaning schedules that keep the efficient reflectors’ performance at good levels. To assess the risk of soiling, monitoring climatic parameters at CSP sites is crucial, namely: AOD [46], precipitations [36, 108], wind speed and direction [108], and relative humidity [53]. Tilted reflectors are less soiled, and glass mirrors -compared to other materials- lose their reflectivity faster and restore it better after cleaning, these properties are helpful in sites where frequent cleaning is possible [36, 46, 98]. Erosion experiments should include dust particles larger than 250 μm, since these particles can be lifted and deposited on the surfaces of reflectors or cause their erosion [43, 56].
5.2. Aerosols’ Effect on Solar Irradiance
Clear-sky irradiance (CSI) is affected by various parameters and CSI modeling is crucial for proper exploitation of solar energy. Ineichen [109] validated seven CSI models over eight years at 22 locations in Europe, Africa and Asia. Parameters other than AOD and Pw have minor effects on irradiation. The three best models were: Solis, REST2 and McClear. Ineichen and Perez [110] provided an air-mass-independent formulation for TL in order to bypass the TL’s dependence upon solar geometry. The developed TL formulation is independent of altitude and AM (Eq. (5)), and is coherent with the previous studies for AM = 2.
I0 is the solar constant, b is a multiplicative factor and AM = 2.
Several papers studied the attenuation of solar energy due to aerosols, like the study of Singh et al. [111] in India. Hazy skies caused very high mean AOD of 1.17; for very hazy days AOD was up to 3.08; high AOD values are associated with close to zero α values, which indicated the dominance of coarse particles; the average α was 0.328 and the lowest α was −0.06 (negative α was correlated with the presence of water insoluble dust). The attenuation of solar irradiance by dust had an average value of −13.6 ±1.4 W/m2 for every 0.1 increase in τ500. In Malta, Bilbao et al. [112] quantified the attenuation of the solar shortwave radiation SW (especially UV, the most energetic fraction of SW radiation), they found that an AOD unit causes UV radiation attenuation by −28.12% to −52.42%, while the absorption of the global SW by the water vapor column was between −2.44 and −4.53 (%/cm). Haywood et al. [40] evaluated the top of atmosphere (TOA) radiative effect of the Saharan dust in West Africa, the attenuation of DNI during dust events was as much as −60±5 W/m2; also, measurements of aerosols’ optical properties were in agreement with Mie’s scattering theory.
Tahboub et al. [113] elaborated a model that evaluates heliostats-receiver irradiance attenuation in UAE, four pyrheliometers recorded one year of data (see Fig. 5), the data were filtered for various reasons (eg: night data, shading) and 55.8% of the recorded data were selected, the model is based on the Beer-Lambert equation and the attenuation of irradiance is as follows:
DNIn is measured by stations 1, 2, or 3, “X” and “d” are, resp, the extinction coefficient and the distance between station n and 4. This model is similar to the one of Sengupta and Wagner [24] although the later was not validated at a specific site according to [113].
Cardemil et al. [38] compared the sizing of a solar tower plant in Brazil by different CSI models, while Pitman & Vant-Hull was in agreement with Sengupta & Wagner, DELSOL didn’t agree with both, the Power Tower Generator (PTGen) program that was used for plant sizing uses DELSOL3 code, which resulted in a 4% underestimation of the attenuation of irradiance due to aerosols and water vapor, thus, the solar field should be 4% larger to avoid long term losses. In fact, DELSOL was developed for the specific elevation of Barstow (0.7 km), CA, USA, corrections should be introduced before applying it to different elevations.
To conclude this section, it is clear from the review of the literature that the radiative effect of aerosols must be accurately evaluated for better long-term exploitation of the solar energy. The performance of the models varies spatially, thus, appropriate corrections should be introduced for better accuracy (as noted by many authors, eg: [38]). The models seem to underestimate high DNI and perform poorly for low solar zenith angles [22, 109], this should be taken into account when developing or applying CSI models to areas with either of the aforementioned climate conditions. Additionally, the complex models are not always more accurate than simple ones (because of the error propagation) [22]. All these factors must be considered to maximize the accuracy of DNI estimation, especially during the planning phase of CSP plants (more similar papers can be found in Table 2). It might also be of interest to rely on next-generation satellites’ data, as they can effectively detect the atmospheric disturbances that happen at high altitudes and are not detected from the ground (eg: large TOA dust event reported by [35]), after all, these disturbances cause non-negligible DNI attenuation and can result in considerable losses.
6. Mitigation Approaches
6.1. Washing
The cleaning of reflectors by rain can be very efficient if it precedes natural adhesive factors, namely high humidity and dew [103]. Niknia et al. [97] studied the impact of cleaning frequencies on the energy yield of a parabolic trough collectors’ plant, the extremes of cleaning frequencies were 7 and 42 days, the respective yearly water consumption was 398.9 and 69.0 m3 and the reflectors’ efficiency losses were 11% and 38%; monthly cleaning consumed 99.7 m3 per year and resulted in a yearly efficiency loss of 27%. Periodic cleaning seems to be the most optimal approach according to Ashley et al. [150], who applied a heuristic approach to define the optimal sub-routing for cleaning the field of heliostats. Picotti et al. [151] developed a model to compute economically optimized cleaning schedules in tower plants; the model is based on the physical behavior of soiling and is applied to solar tower plants in Australia and UAE, the results indicate possible up to 15% gain of the total cleaning cost. Fernández-García et al. [152] tested a variety of washing methods for 2 years under a semi-desert climate, the most effective method was using a brush with demineralized water, the reflectivity was restored to an average of 98.8%; adding detergents to this method did not result in big gains. Bourdon et al. [153] proposed water-saving solutions for CSP plants; first with dust-barriers which stopped 47% of dust particles larger than 25μm. An ultrasonic device (reported to be 4 times more water-economic than the usual cleaning) was used and resulted in a cleanliness factor up to 99%.
6.2. Electrodynamic Screens “EDS”
EDS [154, 155] are an emerging method to prevent CSP mirrors’ soiling, EDS consist of a transparent dielectric film in which alternated electrodes are deposed; when activated with phased voltage pulses, the electrodynamic field prevents dust particles from depositing on the solar mirror and removes the already deposited particles by charging them, lifting them, and propelling them off the reflector [156, 157]. Stark et al. [158] investigated the efficiency of EDS in preventing CSP mirrors’ soiling. Even though EDS application caused an initial reflectivity loss of 3%, EDS-integrated mirrors kept > 90% reflectivity without any cleaning operation. Mazumder et al. [159] reported similar experimental results and highlighted the low cost and energy consumption of EDS, namely $0.005 and 0.2 Wh per m2 per cleaning cycle. After 3 years of study (2017–2020), Mazumder et al. [160] reported that the dust removal efficiency and reflectivity restauration of EDS are both high and that EDS can operate under various temperatures and relative humidity; the authors concluded after an economic analysis that EDS technology is on a commercial level. However, other studies confirmed EDS to be greatly affected by humidity [161, 162], Javed and Guo [163] reported a decreasing dust removal efficiency from 22% at RH of 10% to 9% at RH of 80%, and suggested to increase the frequency of EDS device activation (hourly or bi-hourly) in order to improve its cleaning efficiency. Kawamoto [162] recommended to activate the EDS before the ambient temperature gets lower than the dew point at night. While Bernard et al. [164] proposed to coat the EDS with zinc oxide to shield it from the environmental conditions. Further field studies remain necessary to prove the efficiency of this method for CSP plants.
6.3. Anti-Soiling Coatings
The idea of using coatings to protect the solar reflectors is more mature than that of using EDS (see for example [165]). Polizos et al. [166] performed durability tests on anti-soiling coatings (spray) composed from multifunctional silica particles and polymeric binders, which are hydrophobic coatings with low surface energy and a high contact angle. The initial reflectivity was higher than 90%, accelerated aging showed that anti-soiling spray could maintain good performances for a period of 18 months to two years. Similar investigations by Hunter et al. [25] reported better results: 99% of initial solar weighted reflectance, and two to three years of outdoor durability.
While most anti-soiling coatings are hydrophobic, Aranzabe et al. [167] tested a titania-based hydrophilic coating for CSP anti-soiling purposes in Spain for 43 months. Hydrophilic coatings have a high surface energy and a low contact angle, with a totally wet surface being the extreme case. The mirrors that were coated showed a higher average reflectance difference of 3.3% when compared to the uncoated ones, moreover, the highest reflectance difference was 9.5%. Lorenz et al. [168] studied the economic feasibility of a hydrophilic coating with photo-catalytical functions, a 3% yearly energy yield gain is expected if the coating is applied.
Wiesinger et al. [169] tested different anti-reflective coatings as an erosion mitigation; the coatings improved the optical performance of the tested glass by around 4.9%; however, if the CSP site often suffers extreme weather conditions, anti-reflective coatings would cause faster optical losses. Even though adding a hydrophobic coating to an anti-reflective one slightly reduced optical performance, it was very effective against sand erosion.
Since 2019, more recent papers report innovative methods of water-saving mitigation by coatings. Wette et al. [170] tested a new hydrophilic coating under real outdoor conditions in PSA for 22 months; the effect of the coating was positive, lower soiling rates and higher cleanliness were observed for the coated mirrors in contrast to the uncoated ones. The optical properties and durability of this coating were studied by Fernández-García et al. [171], the coating did not cause significant change in the initial reflectance (−0.001±0.003). The durability was assessed via accelerated aging tests and outdoor exposure, the same previous negligible values were observed, which denotes the durability of the anti-soiling coating and its applicability to CSP reflectors. Matal and Naamane [172] conducted accelerated erosion tests of two types of coated solar reflectors: one with a hard coating, the other with an anti-soiling coating; the results suggest that the hard coating is better in terms of optical performance and resistance to mechanical stress. Wette et al. [173] conducted 6 years of outdoor exposure, as well as laboratory abrasion test, to evaluate anti-soiling coatings. One of the tested coatings improved the initial reflectance while the other one did not change it much; this advantageous behavior diminishes for biweekly-cleaned coatings (2 years) faster than the monthly-cleaned ones (4 years), which confirms that coatings are water saving. Lopes et al. [174] studied horizontal and tilted coated mirrors in Portugal. The coatings’ effect was positive, especially for the tilted samples, while the horizontal ones presented similar soiling rates, coated or not. Pescheux et al. [175] performed extreme accelerated aging tests on different coatings including a commercialized one, the tests simulated severe conditions of salination, humidity, temperature, UV radiations and solar radiation. SiO2 nanoparticles-based hydrophobic coatings performed better than the others, including the commercial one. Sutter et al. [176] studied the impact of dust on tube receivers of the parabolic trough technology for 5 years in Morocco and Spain. Anti-reflective coatings reduced the impact of dust to only −0.014 transmittance loss in Spain, while in Morocco, severe dust events resulted in a transmittance loss of 0.05; accelerated aging tests were also conducted and the highest degradation rate was −0.041. Wette et al. [177] tested two newly developed anti-soiling coatings during 2 years of exposure in Spain, the coatings were cleaned with pressurized water and a brush; the results suggest that the coatings perform differently and dependently on the cleaning technique: for one type, the pressurized water resulted in a cleanliness gain, for the other type it was vice-versa (poor washability). Yilbas et al. [178] presented recently, in 2021, a comprehensive analysis of dust characterization, soiling mechanisms, and dust mitigation by hydrophobic surfaces for solar energy applications.
7. Summary
With the increasing environmental issues caused by greenhouse gases, the transition towards cleaner energy is imperative. One of the most abundant renewable energies is the solar energy, thus its utilization is of prime benefit, especially concentrated solar power (CSP) that is getting increasing interest worldwide. However, the presence of dust in the atmosphere induces an attenuation of the solar energy and a reduction in the efficiency of concentrated solar power (CSP) plants [14, 20]. The globality of this issue is evident by the widespread of the reviewed studies (Fig. 6(a)).
Studies of soiling problems and mitigation approaches in CSP plants: (a) The repartition of the reviewed studies over the continents. (b) Comparison between the MENA region and the rest of the world.
Generally, the regions blessed with high DNI suffer from high aerosols loads, this is especially the case of the MENA region [26, 179, 180], this explains the higher number of the -reviewed- papers studying CSP-related soiling problems in the MENA region (Fig. 6(b)). Soiling is an issue for solar energy applications in general, however, CSP technologies suffer from soiling problems more than PV [51]. The review of the literature shows that no previous review was dedicated to the effect of dust on CSP (unlike PV), and that very few research papers -and no review papers- were interested in both effects of suspended aerosols and deposited dust, or linked turbidity to reflectivity loss (like the studies of [18, 46]), although both effects ought to be simultaneously investigated, since turbid skies accompany the deposition of dust on solar mirrors. Shortage in long term energy production results from the underestimation of both DNI extinction (eg: [38]) and reflectivity loss of the solar field (eg: [50]). Aerosols’ attenuation of DNI is reported to be up to 25% [122], and the soiling of solar mirrors causes up to 80% reflectivity loss [57, 102], while the resulting diminution in energy production can be more than 88% [107].
Attenuation of the solar energy by aerosols is primarily via scattering more than absorption [41]. The major parameter identifying the attenuation is the Aerosols Optical Depth (AOD) [123], it is found to particularly be the input parameter with crucial impact on CSI models’ outputs [83], as well as on extinction models used in solar tower plants [76]. Introducing measured real-time AOD data into models does not improve their accuracy [139]. High AOD values are correlated with low Ångström exponent α [111, 181], thus large aerosols particles result in very hazy atmospheres and are responsible for most of the extinction [18]. The highest measured AOD range from 0.30 (China [133]) to 3.08 (India [111]). To consider the particularities of sites, and since extinction levels in CSP plants are specific to the hosting site [29], it is better to build models for sites with potential for CSP plants (eg: [65, 105]), or introduce the necessary corrections (eg: [38, 120]). The models underestimate high DNI and perform poorly for low solar zenith angles [22, 109]. The complexity of models is not always an assurance of better performance, sometimes, simple models perform better than complex ones [22], this is because complex models have more error propagation than the simple ones. It might be of interest to depend on data of next-generation satellites, which can detect important atmospheric disturbances that take place at high altitudes and escape ground measurements [35]. The attenuation of solar irradiance by atmospheric aerosols occurs via two optical phenomena: scattering and absorption [37], even though other atmosphere constituents contribute in scattering, most of it is primarily due to aerosols [38]. Theoretically, the optical phenomena happening in the lower atmosphere are best described by the well-known Mie theory [40].
The deposition of dust on solar reflectors is mainly due to gravity more than diffusion [46]. Since AOD is proportional to the soiling rate of solar mirrors [18], monitoring AOD in the surroundings of CSP plants ought to be done, in order to gauge the risk of soiling [46]. Intense, rare soiling events cause much more important reflectivity drops than the frequent, moderate dust deposition [49, 62]. Many authors developed innovative methods to measure soiling rates [71–75, 115]. Soiling of pyrheliometers causes erroneous evaluation of DNI in potential CSP sites [119]. First row mirrors are soiled fast and uniformly, and suffer more dust deposition than the rest of the solar field [53], this should be considered when scheduling cleaning operations. The tilt angle and material of solar reflectors have a crucial impact on the soiling rate: horizontal mirrors are always more soiled than tilted ones, thus, a tilt angle, no matter how slight it is, can reduce the deposition of dust; also, glass mirrors are soiled faster than aluminum ones, and restore their reflectivity faster and better after being washed [36, 46]. Soiling rates can reach high values (eg: −2.3% per day [53]), it would be of benefit to keep the factors controlling dust deposition under observation; other than the exposure direction, multiple meteorological factors are directly related to soiling events, namely: wind speed and direction, relative humidity and atmospheric pressure; while high relative humidity values intensify dust deposition, low pressure is associated with rainy days [46, 47]. Monitoring these factors would help program efficient cleaning operations.
Many authors developed and validated soiling models for CSP applications (eg: [65, 69]). The consideration of soiling when modeling the performance of CSP plants resulted in considerable losses of the yielded energy [48]. The weight of the deposited dust can also be used to develop correlations of reflectivity losses [64–67].
Other than the deposition of dust on the reflective surfaces, the optical losses are also caused by the erosion of the solar reflectors, mainly due to sandstorms [60–62]. Erosion intensity varies from site to site, as a result of the variation of particles’ shapes that cause different microscopic cracks. The erosion’s intensity is affected by the shape of particles more than their hardness [56–58]. The general trend of erosion simulations is to use particles smaller than 250 microns, but larger particles should also be used to avoid erosion underestimation [56].
The primary effect of rain is cleaning the reflecting mirrors more than enhancing their soiling, the latter is mainly due to very light rain [29, 44]. The highest possibility of soiling is during the early morning, which makes that timing favorable for maximum profit of cleaning operations [118]. Down-facing mirrors remain clean, which makes stowing mobile solar reflectors a good preventive approach [36, 107]. As for materials, the reflectivity of glass mirrors is perfectly restored after cleaning operations [53]. Most mitigation approaches are restorative and water-based, despite the shortage of water in the regions blessed with high DNI. Dry preventive methods include mainly coatings and electrodynamic screens “EDS” [30]. EDS are a water-saving approach. With their high dust removal efficiency and reflectivity restauration, EDS are becoming ready for large-scale use [160]. Coatings of solar mirrors enhance water saving since less cleaning frequencies are required; and the coated reflectors restore their reflectance better than the uncoated ones [173]. Furthermore, the coatings’ anti-soiling effect will be more significant when the soiling is more intense [170]. Coatings can be associated with cleaning methods to enhance the optical performance of the reflectors; however, this association of mitigation approaches must be done dependently on the type of coating, different anti-soiling coatings perform differently for cleaning methods, the wrong association can lead to losses in the mirrors’ reflectance instead of restoring it [177].
To finalize this section, aerosols and soiling problems still remain in need for further research. The resulted decrease in the solar power plant’s overall performance might increase the cost of cleaning operations and therefore increasing the cost of the produced kWh.
The amount and composition of dust -and the consequent losses-relate to the geography and climate of each region in the world, thus the appropriate mitigation approaches are site-dependent, this explains the recent R&D trend of evaluating experimentally the factors controlling soiling losses in different regions of the world rather than standardizing those losses. The evaluation of the amount of dust is very important to assess the cost of cleaning tasks, and thus assess the viability of CSP projects around the world. The cost of cleaning solar fields must be optimized by adapting the cleaning method to the factors governing the deposition of dust at the plant’s site, the careful association of different mitigation approaches can lead to more efficient cleaning operations. The usual scarcity of water in areas with high DNI potential urges the necessity for more research on water-economic solutions. Future research should also investigate more the relationship between AOD and reflectivity drop, this can help in the early prediction and prevention of solar mirrors’ soiling (by stowing mobile reflectors for example). It is recommended to build site-dependent soiling and extinction models for precise estimation of the optical losses.
8. Conclusions
CSP plants suffer production losses due to aerosols-related problems, more than any of the other solar energy technologies; mid-air suspended aerosols result in attenuation of the incident direct normal irradiance, furthermore, the deposition of dust on the reflecting mirrors causes serious reflectivity drops. A literature review of the studies concerned with both DNI attenuation and reflectivity loss is presented, indeed, investigations reported aerosols-related problems to be of serious impact on the production of energy in CSP plants. The main parameters characterizing the atmospheric turbidity are: Aerosols Optical Depth, Linke turbidity factor, and Ångström coefficient and exponent. The state of solar mirrors can be expressed through the soiling rate or the cleanliness factor, these parameters are influenced by many factors, mainly tilt angle, wind speed and direction, and relative humidity.
Soiling issues and mitigation strategies still require more research and development. The research, development, and challenges connected to dust issues for CSP plants are examined and summarized in this paper. The review has offered a thorough examination and overview of studies on DNI attenuation and reflectivity losses, as well as the main characteristics of air turbidity, factors of soiling, and mitigation strategies. The following main conclusions can be drawn:
Dust and soiling are still issues that need to be addressed. This is particularly true in the world’s desert-sand regions, which, ironically, have some of the best solar conditions as well as some of the "best" dust conditions. The lack of natural washing from rain and the small amount of moisture that does fall might result in even more severe soiling.
Concentrating solar systems’ energy delivery can be considerably reduced when dust is present. Small reductions in a mirror’s or heliostat’s reflectivity result in significant reductions in CSP performance. Concentrating systems, in general, are more susceptible to dust collection and require greater maintenance in terms of dust mitigation.
Mitigation strategies have been discussed including restorative approaches and preventative approaches. Passive methods (which use coatings to prevent dust attachment) and active methods (which actively oppose the charged dust particles) were used as preventative measures. The best mitigation strategy is determined by the geographic location and climate zone. Further research on water-economic mitigation approaches (restorative and preventive approaches) is needed
It is best to build site-specific extinction models to consider the particularities of the sites hosting CSP plants.
More papers should investigate the relationship between AOD and reflectivity drop, this can lead to the early prediction and even prevention of the reflector’s soiling.
Glossary
α
Ångström turbidity exponent
β
Ångström turbidity coefficient
Δz
Height difference
δ
Declination angle
δcda
Optical thickness of the clear-dry atmosphere
λ
Wavelength
vd
Deposition velocity
vg
Free fall velocity
vt
Turbulence velocity
τ
Optical Depth
Cd
Dust concentration
D
Particle diameter
d
Distance
F
Aerosol Forcing
Fd
Dust deposition flux
I0
Normal incident extraterrestrial irradiance
Nj
Julian day number
n
Refraction index
P
Pressure
ra
Aerodynamic resistance
rb
Brownian diffusion resistance
Sfc
Surface level
TL
Linke’s turbidity factor
X
Extinction coefficient
Z
Solar Zenith Angle
Supplementary Information
Notes
Author Contributions
K.Z. (Ph.D student) collected data and wrote the manuscript. A.G. (Assistant Professor) supervised, wrote sections 4 and 6 and revised the manuscript. M.A. (Professor) supervised and wrote sections 4 and 6. N.R. (Ph.D) wrote sections 1–3, 7, 8. F.Y. (Assistant Professor) wrote sections 1–3, 7, 8. N.L.P. (Assistant Professor) wrote section 5, revised the manuscript and checked the language.