1. Introduction
In general, global horizontal irradiation (GHI) data can be obtained from routine ground-based measurements. However, as solar power plants are integrated into the electricity grids, it is important to forecast how much solar energy will be produced two or three days ahead. The management of total power production considering different sources and the decision making in the energy market depend on the forecasting of GHI. This information is helpful for grid operators in order to better accommodate the variable generation of electricity in scheduling and regulating power. Studies on the impact of clouds and aerosols on GHI are of great importance for reliable GHI forecasting. GHI depends on cloudiness, relative humidity and aerosols; therefore, it is highly variable (Ohtake et al., 2012, 2015; Pedro and Coimbra, 2012; Diagne et al., 2014; Neves and Silva, 2014; Zempila et al., 2016). For example, Zempila et al. (2016) evaluated four short wave (SW) radiation schemes of the Weather Research and Forecast (WRF) model in terms of GHI and found that for all-sky conditions hourly GHI-based on WRF is overestimated for all schemes. Furthermore, all the schemes are found to perform better for clear skies. Ohtake et al. (2015) showed that seasonal and regional variation of cloud types are related to relatively large GHI forecast errors. Thus, accurate short-term solar forecasts are critical for solar power planning and production, and for an efficient use of the fluctuating energy output. Inman et al. (2013) and Diagne et al. (2013) reviewed the current methods used to forecast solar irradiance.
Numerical weather prediction (NWP) models are widely used for weather forecasting and atmospheric research. Some NWP models combine numerous options for physical parameterizations of radiation, cloud microphysics, planetary boundary layer (PBL), surface layer fluxes, turbulence, cumulus convection and others for the representation of sub-grid scale processes. A meteorological model includes parameterizations of radiation and clouds. In general, there are several schemes available to represent each of these processes. NWP models have been proven to be powerful tools for solar radiation forecasting (Lorenz et al., 2009a, b; Rincón et al., 2011; Lara-Fanego et al., 2012; Diagne et al., 2013; Pérez et al., 2013). Early modeling attempts for estimating GHI were based on the fifth-generation Penn State/National Center for Atmospheric Research (NCAR) Mesoscale Model (MM5) (Grell et al., 1998; Zamora et al., 2003, 2005; Heinemann et al., 2006). In the study of Rincón et al. (2011), the WRF model was applied with 4 × 4 km horizontal resolution over the Iberian Peninsula. These models predict GHI by resolving physical processes and using radiation transfer models (RTM) (Morcrette et al., 2007). Studies were reported for various locations of the USA (Zamora et al., 2003, 2005; Mathiesen and Kleissl, 2011) and Europe (e.g., Germany and Spain) (Heinemann et al., 2006; Lorenz et al., 2009a, b). Mathiesen et al. (2013) mentioned that most operational NWP models consistently overpredicted irradiance. Moreover, Lara-Fanego et al. (2012) assessed the WRF-derived solar irradiation predictions in southern Spain and found a striking dependence between model skill in reproducing GHI and cloud cover. Recently Ruiz-Arias et al. (2013) used high quality ground data to test the radiative closure in WRF for USA with the RRTMG (a modified version of the well-known Rapid Radiative Transfer Model [RRTM] scheme in WRF) and Goddard SW radiation schemes in clear-sky conditions, while Zempila et al. (2016) evaluated the differences induced in the WRF model predictions in Greece when using different short wave radiation schemes.
There are several studies indicating the usage of solar energy resources as an alternative to conventional energy in Turkey with a promising prospect for the future (Taşdemiroǧlu and Sever, 1986; Topcu and Oney, 1994; Bulut and Büyükalaca, 2007; Bakirci, 2009; Alta et al., 2010; Incecik et al., 2012). Since Turkey is situated between 36º N and 42º N latitudes, it has been recognized as a region with abundant solar radiation in the Mediterranean belt. The southeastern Anatolia (SEA) region has a highly favorable geographical location for the utilization of solar energy, similar to southern Europe and Mediterranean countries such as Spain, Portugal, France, Italy and Greece. According to the values presented by the General Directorate of Renewable Energy in Turkey (GDRE, 2018), the SEA region in Turkey is the most significant one in terms of annual total solar radiation and sunshine duration.
The aim of this study is to analyze the forecasting skills of the WRF model with four different configurations of physical parameterizations and scheme combinations. Model evaluation is performed for 72-h GHI forecasts in the SEA region of Turkey under two different cloud conditions during May 8-10 and August 20-22, 2011. The evaluated forecasts of these two cases were based on hourly values measured on August 20-22, 2011 for clear days (the clear sky condition refers to the absence of clouds) and May 8-10, 2011 for partly cloudy days. The stations selected for this study were Kilis, Ceylanpınar, Bozova, Şırnak, and Mardin in SEA. The meteorological data for these stations were obtained from the Turkish State Meteorological Service (TSMS).
2. Methods
2.1. Study area
The SEA region is located at 36-38º N and covers 57 200 km2, approximately 7.5% of the Turkish territory. The region presents a diverse topography. Its southern part is an almost homogeneous, flat area with a mean elevation of around 744 masl. This area extends around the lower Tigris River basin. On the other hand, the eastern part of the region presents a very complex topography, with steep elevation gradients reaching altitudes over 3800 m.
In terms of GHI, the SEA region is the most significant in Turkey with annual total solar radiation of 1460 kWh m-2 and 2993 h of sunshine, while average values for Turkey are 1311 kWh m-2 y-1 and 2640 h, respectively (GDRE, 2018). Kilis, Ceylanpınar, Bozova, Şırnak, and Mardin stations have GHI measurements in SEA. Table I href="?a5-143-158"> shows the geographical characteristics of these stations.
2.2. Cloudy conditions
Clouds and their accompanying weather patterns are among the most important atmospheric phenomena limiting solar radiation at the Earth’s surface. GHI has a nonlinear dependency on cloud cover due to its dynamic chaotic characteristics (Paulescu and Badescu, 2011; Chiacchio and Vitolo, 2012). Therefore, clouds are an important source of uncertainty in the prediction of solar irradiance.
Cloudiness or cloud fraction is a meteorological measurement traditionally made by human observers. Cloud cover measurements are reported by meteorological services in SYNOP and METAR reports (González et al., 2012). Because of their subjectivity, these measurements are not always accurate.
The daily sunshine index (s/S) was considered as a criterion, in order to determine the cloud conditions of the study periods. s/S was obtained by normalizing s, the bright sunshine duration, with S, the astronomical day length. The values of daily index greater than 0.75 were considered descriptive of clear-sky conditions and values between 0.25 and 0.75 were considered as partly cloudy conditions. Simulating solar radiation in clear sky conditions (the absence of clouds) is generally easier than in partly-cloudy ones.
2.3. Model setup
The WRF model (Skamarock et al., 2008; Wei et al., 2008) used in the present study (v. 3.3) was integrated using three nested domains with grid spacing values of 36, 12, and 4 km. The SEA region was covered with 4-km grid resolution, as shown in Figure 1.
The simulations were driven based on the National Center for Environmental Prediction (NCEP) Global Forecast System (GFS) data. The temporal resolution of the GFS analysis was 6 h, while the spatial resolution was 1º. Model configurations are given in Table II. The atmospheric column was decomposed into 35 vertical levels. Only the innermost domain was used in the evaluation. The WRF model ran for 72 h in each case. The evaluations of the 72-h GHI forecasts were based on the hourly values measured on August 20-22, 2011 for clear days and on May 8-10, 2011 for partly cloudy days for the given stations in SEA. These dates are chosen after examining observations of all stations during 2011. The model uses two-way nested boundary conditions. No cumulus parameterization was used for domain 3.
Domain 1 | Domain 2 | Domain 3 (SEA) | |
Horizontal grid interval | 36 km | 12 km | 4 km |
Number of horizontal grids | 146 × 108 | 220 × 133 | 157 × 97 |
Vertical levels | 35 | 35 | 35 |
Initial conditions (IC) / boundary conditions (BC) | GFS 1.0, 6-h intervals | Domain 1 | Domain 2 |
Nesting | Two-way | Two-way | Two-way |
Cumulus | Kain-Fritsch | Kain-Fritsch | - |
The WRF model offers numerous physical options such as PBL, surface physics, radiation, microphysics and cumulus parameterization. The WRF model was operated with four different configurations defined as WRF0, WRF1, WRF2 and WRF3. The physical parameterizations and scheme choices used in this study are presented in Table III.
Parameterizations | WRF0 | WRF1 | WRF2 | WRF3 | |
Radiation | SW | RRTMG | RRTMG | FLG | Goddard |
LW | RRTM | RRTM | FLG | Goddard | |
Microphysics | WSM6 | Thomson | Thomson | Goddard | |
Surface layer | Monin-Obukhov with Carlson-Boland |
Monin-Obukhov with Carlson-Boland |
Monin-Obukhov with Carlson-Boland |
Monin-Obukhov with Carlson-Boland |
|
Surface physics | Noah | Noah | Noah | Noah | |
Boundary layer | YSU | YSU | YSU | YSU | |
Cumulus | Kain-Fritsch | Kain-Fritsch | Kain-Fritsch | Kain-Fritsch |
Radiation is the main force that regulates the surface energy budget. An atmospheric RTM calculates the radiative transfer of electromagnetic radiation through the atmosphere. Both SW and longwave (LW) radiation parameterizations in an atmospheric model calculate radiation fluxes and heating rates (Price et al., 2013). Specifically, radiation schemes resolve atmospheric heating from radiative flux divergence and surface downward LW and SW radiation. The Rapid Radiative Transfer Model (RRTM) for general circulation models (GCMs) (RRTMG) is the modified version of the RRTM for SW solar irradiance (Clough et al., 2005; Iacono et al., 2008). In the present study, RRTMG (Iacono et al., 2008), RRTM (Mlawer et al., 1997), Fu-Liou-Gu (FLG) (Gu et al., 2011), and Goddard radiation schemes (Chou and Suarez, 1994) were used for SW and LW with different microphysics groups in four configurations. Table IV presents a comparison of these four radiation schemes.
Characteristics | RRTMG | RRTM | FLG | Goddard | |
Number of bands |
SW | 14 | 14 | 6 | 11 |
LW | 16 | 16 | 12 | 10 | |
Gases (optical parameters) |
H2O, CO2,
O3, N2O, CH4, O2, and four types of halocarbons: CFC-11, CFC- 12, CFC-22, CCl4 |
H2O, CO2,
O3, N2O, CH4, O2, CO, and four types of halocarbons: CFC-11, CFC-12, CFC- 22, CCl4 |
H2O, CO2,
O3, N2O, CH4, O2, CO, NO, SO2, NO2, and three types of halocarbons: CFC-11, CFC-12, CH3Cl |
H2O, CO2,
O3, N2O, CH4, and two types of halocarbons: CFC-11, CFC-12 |
|
Gaseous absorption | Correlated k-distribution method |
Correlated k-distribution method |
Correlated k-distribution method |
Correlated k-distribution method |
|
Cloud overlap | Maximum-Random cloud overlapping |
Maximum-Random cloud overlapping |
Maximum-Random cloud overlapping |
Maximum-Random cloud overlapping |
|
Calculation of flux |
Two-stream algorithm |
Four-stream adding DISORT |
Delta 4-stream (SW) Delta 2-stream (LW) |
Two-stream algorithm |
For example, the number of bands for RRTM and RRTMG are equal to 14 for SW and 16 for LW, respectively. RRTM performs the calculations through the Discrete Ordinates Model (DISORT) for radiative transfer using four streams. As seen in Table IV, both models employ gaseous absorption coefficients extracted from the line-by-line radiation model and the computationally efficient correlated-k method for radiation transfer calculations (Goody et al., 1989). Furthermore, the RRTMG scheme uses a two-stream algorithm for multiple scattering and then uses reduced set of g-intervals for integration over the extinction in each band (Table IV). Moreover, the RRTMG-SW includes the Monte-Carlo Independent Column Approximation (McICA) method to represent sub-grid cloud variability (Barker et al., 2006, 2008). Ruiz-Arias et al. (2013) have shown that the RRTMG has high potential skill for long-term assessment of clear-sky GHI. Furthermore, Mlawer et al. (1997) showed that the RRTM-LW is an accurate and efficient correlated-k, LW RTM. FLG (Fu and Liou, 1993) uses 6 SW and 12 LW spectral bands. FLG uses Delta 4-stream approximation for SW and Delta 2-stream approximation for LW.
The Goddard-SW solar radiation model has been developed at the NASA (Chou and Suarez, 1999). Goddard-SW includes extinction by water vapor (H2O), ozone (O3), carbon dioxide (CO2), nitrous oxide (N2O), methane (CH4), and CFC-11, CFC-12 (Chou et al., 2001). The SW solar radiation model contains attenuation by H2O, O3, oxygen, CO2, aerosols and clouds and the fluxes are calculated at 11 spectral bands, using the two-stream added approximation. Its accuracy is expected to be within a few W m-2, whereas the atmospheric heating rate between 0.01 hPa and the surface is accurate to within 5%, relative to line-by-line calculations (Chou and Suarez, 1999). Iacono and Nehrkom (2010) showed that the Goddard-SW model produces downward SW surface fluxes that are about 30 W m-2 higher than RRTMG-SW. Moreover, the Goddard-LW uses 10 spectral bands. In this study, we employed the newest Goddard scheme under WRF v3.3.
A proper representation of clouds in solar irradiance prediction is rather problematic. Mathiesen et al. (2013) showed that the forecast accuracy is low in case of the regions with dynamic cloud systems. In the present study, the WRF Single-Moment 6-Class (WSM6) (Hong and Lim, 2006) and the Thompson (Thompson et al., 2004) microphysics schemes were used for representing microphysical processes in clouds while only the Kain-Fritsch scheme (Kain and Fritsch, 1990; Kain, 2004) was considered for cumulus parameterization.
In microphysics, the WSM6 scheme uses recent scientific developments regarding ice nucleation and includes a dependence on aerosol concentration (Skamarock et al., 2008). Besides, the scheme includes a parameterization for calculating the direct radiative effects of several types of aerosol, including maritime, continental, urban, mineral dust, soot (black carbon), and sea salt. The WSM6 has been one of the options of microphysical processes in the WRF model since 2004. It is based on the revised ice-microphysics suggested by Hong, Dudhia and Chen (HDC) (Hong et al., 2004; Hong and Lim, 2006). In addition to the characteristics mentioned above, the WSM6 behaves realistically in response to the appropriate grid resolvable forcing (Hong and Lim, 2006).
The Thompson scheme, developed for use in WRF and other mesoscale models, is another scheme used for parameterization of microphysics. The Thompson microphysics scheme includes improvements to the earlier bulk scheme of Reisner et al. (1998). It was tested and compared with idealistic case studies, specifically for mid-latitudes (Hall et al., 2005).
There are many different PBL parameterization schemes in the literature (Kleczek et al., 2014). Many of them are designed to perform well during specific weather conditions. The variation of the predicted parameters such as temperature, mixing height and GHI with different parameterization schemes were examined by several investigators. For example, Shin and Hong (2011) compared five PBL schemes in the WRF model for a single day from a field program. They reported a large variation in the mixing height values computed by the five schemes both in daytime and nighttime. Hu et al. (2010) compared the 3-month mean diurnal variations of the simulated mixing heights, using different combinations of parameterization schemes in the WRF model with the observed values for Southeast Texas. Similarly, Han et al. (2008) evaluated five PBL schemes in the MM5 model for the East Asian domain and reported a large difference in the mixing heights predicted by the various combinations of parameterization schemes. Ruiz-Arias et al. (2008) conducted an evaluation study for the performance of different parameterizations of the model. In this study, the model’s performance was found highly dependent on the sky conditions. For example, solar radiation RMSE values were about one order of magnitude higher during broken-clouds and overcast conditions compared to clear-sky ones. More recently, Zempila et al. (2016) evaluated the shortwave radiation schemes such as Dudhia, updated RRTMG, updated Goddard and the GFDL schemes based on the WRF-model derived GHI predictions for cloudy and cloudless conditions. All SW radiation schemes presented better GHI predictions during the warm period due to the reduced presence of clouds.
The Monin-Obukhov similarity theory with Carlson-Boland viscous sublayer was also used for representing surface-layer parameterization; the Noah model was selected as land surface model for surface physics. The Yonsei University (YSU) scheme for PBL (Hong and Lim, 2006) and the Kain-Fritsch scheme for cumulus parameterizations (Kain and Fritsch, 1990) were used.
2.4. Error metrics
In the present study, we employed two error metrics for the evaluation of the model simulations according to the observations: mean bias error (MBE) and root mean square error (RMSE).
where x e is the estimated data and x m is the measured data. In addition to these metrics, normalized values rMBE and rRMSE are defined as relative MBE and relative RMSE, respectively. Nearest grid-point is used for verification with the station observations.
3. Results and discussion
This article investigates the forecasting skills of solar radiation by the WRF model with different configurations of physical parameterization and scheme combinations, both in terms of partly-cloudy and clear-sky conditions at each site in the SEA region. This evaluation will improve our understating of the strengths and weaknesses of the four different configurations of the WRF model in forecasting hourly GHI values in the region.
In this study, the three days ahead hourly GHI simulations based on the four different configurations were evaluated for the two 3-day simulation periods representing different cloud conditions including May 8-10, 2011 (partly cloudy) and August 20-22, 2011 (clear sky) in the regions. Model predictions were compared to GHI measurements from installed pyranometers at five meteorological stations in SEA. The results are structured and analyzed in two parts as follows:
3.1. Partly cloudy case (May 8-10, 2011)
The period May 8-10, 2011 during which the s/S varied between 0.25 and 0.75, was selected as the partly-cloudy case. Seventy two-hour simulation results of the model based on the four different configurations named WRF0, WRF1, WRF2 and WRF3 for partly cloudy days at five stations are shown in Figure 2.
The MBE and RMSE are computed for all stations and on an average throughout the SEA region according to Eqs. (1) and (2); the results are shown in Table V. The rMBE and rRMSE values at each station for all configurations in partly cloudy conditions are calculated and shown in Figures 3 and 4, respectively. MBEs are found positive indicating overestimation of GHI for the partly cloudy conditions at five stations and for all configurations (Table V). The findings are supported by the results of Zempila et al. (2016). The overestimated WRF-derived hourly GHI values can also be clearly seen in Figure 2. As seen in Table V, MBE varies from 49.49 (WRF2) to 80.98 W m-2 (WRF3) depending on the configuration and site of the region for the partly cloudy situation.
Site | Statistics | WRF0 | WRF1 | WRF2 | WRF3 |
Bozova | MBE | 62.06 | 64.58 | 58.88 | 73.89 |
RMSE | 139.49 | 141.54 | 139.94 | 151.40 | |
Ceylanpınar | MBE | 63.76 | 68.40 | 67.88 | 80.98 |
RMSE | 177.23 | 183.16 | 183.30 | 191.17 | |
Kilis | MBE | 56.70 | 58.64 | 49.49 | 68.57 |
RMSE | 148.89 | 151.44 | 152.44 | 159.48 | |
Mardin | MBE | 53.61 | 63.12 | 53.11 | 70.13 |
RMSE | 173.52 | 179.68 | 174.35 | 175.98 | |
Şırnak | MBE | 61.09 | 59.72 | 54.24 | 68.63 |
RMSE | 122.32 | 121.87 | 121.07 | 131.07 | |
Average | MBE | 59.44 | 62.89 | 56.72 | 72.39 |
RMSE | 152.29 | 155.54 | 154.22 | 161.82 |
Note: minimum errors are shown in bold.
Moreover, the partly-cloudy days on May 8-10, 2011 are linked with high RMSEs. For example, the regional average RMSE for the scheme configurations varied between 152.29 and 161.82 W m-2 at SEA. These results are in agreement with the findings of Mathiesen and Kleissl (2011), who also presented higher biases during partly-cloudy conditions in the USA. The MBE values reported for NAM forecasts under moderately-cloudy conditions in the USA are positive and vary between 10.8 and 62.6 W m-2 while the RMSE values vary between 93.5 and 149 W m-2.
Furthermore, from the relative view of the errors, the rMBE varies between 17 and 30% depending on the configuration and site while the rRMSE varies between 41 and 70%, respectively (Figs. 3 and 4). The difference between the forecasting values of WRF-derived hourly GHI and the observations shows that the WRF3 combination has the highest error among the four configurations in the region for partly-cloudy situation (Fig. 2). The lowest RMSE was found in Bozova (in the range of 139.49 and 151.40 W m-2, depending on the configuration) while the highest was found at Ceylanpınar (177.23 and 191.17 W m-2) and Mardin (173.52 and 179.68 W m-2).
In summary, the WRF model with the four configurations used in the study did not adequately represent the fluctuations in the partly-cloudy case. This may be explained by the impact of clouds and aerosols (Sezen et al., 2013). Although the impacts of the aerosols were not considered in this study, spring is the favorable season for dust transport from the southern desert areas to many regions of Turkey including SEA. Specifically Mardin is very well known as a dusty city. As a final remark for the partly-cloudy, case the WRF0 configuration had a better performance compared to other configurations.
3.2. Clear-sky case (August 20-22, 2011)
In this study, the August 20-22, 2011 period has been selected to represent the clear sky situations using s/S values higher than 0.75. The modelled and measured GHI values for 72 h at SEA sites are shown in Figure 5. As seen from Table VI, MBEs are positive, indicating overestimation of GHI for the clear-sky days at the five stations and for all configurations, similar to the partly cloudy case in May. The results of the WRF model at each site in the region, in terms of performance of each statistical metric, are also shown in Figures 6 and 7. The comparison of the four model configurations presents significant differences in the relative MBE and RMSE at the five sites in the region. The relative MBE is between 10 and 18% and the relative RMSE is between 20 and 39% throughout the region depending on the WRF model configuration (Figs. 6 and 7). Furthermore, from a regional perspective, the average MBE and RMSE values varied in the ranges of 36.73-48.26 W m-2 and 71.24-81.13 W m-2, respectively, depending on the configuration.
Site | Statistics | WRF0 | WRF1 | WRF2 | WRF3 |
Bozova | MBE | 45.03 | 45.90 | 45.27 | 55.76 |
RMSE | 90.81 | 91.43 | 92.65 | 102.72 | |
Ceylanpınar | MBE | 46.97 | 46.28 | 44.86 | 55.32 |
RMSE | 71.75 | 69.10 | 69.23 | 82.04 | |
Kilis | MBE | 40.97 | 33.69 | 26.09 | 50.83 |
RMSE | 81.70 | 79.77 | 87.41 | 93.37 | |
Mardin | MBE | 34.68 | 34.72 | 31.47 | 41.33 |
RMSE | 57.71 | 57.63 | 57.43 | 70.19 | |
Şırnak | MBE | 39.89 | 39.54 | 35.98 | 47.94 |
RMSE | 58.69 | 58.26 | 55.42 | 70.11 | |
Average | MBE | 41.51 | 40.03 | 36.73 | 48.26 |
RMSE | 72.13 | 71.24 | 72.43 | 81.13 |
Note: minimum errors are shown in bold.
According to the statistical metrics, the model performance is better for the clear sky case than the partly cloudy conditions, as expected.
GHI forecast by the WRF2 model configuration is close to the observations particularly at Şırnak and Mardin, with the lowest RMSE values of 55.42 and 57.43 W m-2, respectively (Table VI). These sites are located at the rugged mountainous area with higher altitudes in SEA (Table I). RMSE values are supported by the findings of Ohtake et al. (2015) in Japan using the Japan Meteorological Agency (JMA) mesoscale model.
In the WRF2 model configuration, rMBE varies between 10 and 17% while rRMSE varies between 20 and 35% depending on the site within the study region. These results by the WRF2 configuration are in agreement with the results by Diagne et al. (2014), who reported rMBE of 6.34% and rRMSE of 37.18% at Saint Pierre together with an rMBE of 16.08% and rRMSE of 57.3% at Tampan, two sites on Reunion Island in the Indian Ocean. Generally, the statistics of WRF2 indicated slightly better performance than WRF0 and WRF1 and significantly better performance than WRF3 for clear sky conditions. WRF3 presents poor results compared to the other configurations. Hence, we can emphasize that it is a combined effect of all parameterization schemes which translates into a good simulation of GHI by the WRF model except for the WRF3 combination.
4. Conclusions
In the present study, the skill of the WRF model in forecasting hourly GHI values in the SEA region of Turkey was assessed at five sites and under two different cloud conditions with four different combinations of parameterization schemes to improve our understanding of which schemes perform better.
The skill of the WRF model was evaluated through the combination of MBE and RMSE at the five sites in the region. The performance of the GHI forecasting up to three days ahead based on the WRF mesoscale model applications including physical parameterizations (radiation, microphysics, PBL, surface layer, surface physics, and cumulus) with four different scheme configurations were examined in this region. The hourly values of solar radiation measured during the periods of May 8-10, 2011 (for partly cloudy conditions) and August 20-22, 2011 (for clear-sky conditions) were used to evaluate the performance of the various configurations.
The numerical experiments carried out to select the optimum combination of parameterization schemes for estimating GHI have shown that the configuration for physics option consisting of RRTMG as the SW radiation, RRTM as the LW radiation, WSM6 or Thomson as the microphysics, YSU as the PBL, Monin-Obukhov with Carlson-Boland as the surface layer, and Noah as the land-surface model performs reasonably well in reproducing the observed GHI, especially for clear-sky conditions in southern Turkey.
The role played by the radiation model for a better simulation of the SEA region under clear-sky conditions was also emphasized. However, no such conclusive statement can be made for the partly-cloudy cases. Studies for the proper choice of a parameterization scheme are useful when a prognostic weather model like WRF is selected to forecast the amount of solar energy that will be produced at a photovoltaic power station (PV plant) site. The results are as follows:
It was found that the forecasted hourly GHI values by all configurations of the WRF model tended to be overestimated for all sites and for both partly-cloudy and clear-sky cases. The regional average of the MBE values ranged from 56.7 to 72.4 W m-2 for the partly-cloudy case and from 36.7 to 48.3 W m-2 for the clear-sky period. There were no significant differences in the MBE values between the sites except with the WRF3 combination for both cases. The systematic overestimation found in both cloud-cover cases and with each configuration may be related to the solar radiation parameterizations of the WRF model.
The RMSE values of the WRF model showed a strong dependence on cloud cover over the study region. The RMSE values decreased from May under partly-cloudy to August under clear-sky conditions at all stations in the region. For example, the spatial average of RMSE for each WRF configuration in the region showed about 50% reduction from partly-cloudy (in May) to clear-sky conditions (in August). MBE also decreased from partly-cloudy to clear-sky conditions by about 30% in the region. Lara-Fanego et al. (2012) showed a similar result for Andalusia, Spain for a 24-h lead time.
The lowest rRMSE values in hourly GHI forecasts were found for the clear-sky case (August), in the range of 20-39% depending on the WRF configuration, while the maximum error occurred in WRF3. The rRMSE was 41-70% on partly-cloudy days. The results of the WRF0 configuration were slightly better than the others for the partly-cloudy case while the results for WRF1 configuration were much better than the others for the clear sky case. The only difference between WRF0 and WRF1 is the microphysics scheme (WRF0 uses WSM6, while WRF1 uses the Thompson scheme). In the partly-cloudy case, the rRMSE values were alike, except for the WRF3 configuration, which includes the Goddard scheme. The WRF model with the Goddard scheme, when compared to others, yields poor results for both partly cloudy and clear-sky cases.
Clouds are the biggest source of uncertainty for GHI forecasting. Extensive satellite and ground-based cloud measurements are needed for the evaluation.
In summary, the presented results show that the performance of WRF in predicting solar irradiance in the SEA region depends highly on the combination of the physical parameterization schemes in both clear-sky and partly-cloudy conditions. Additionally, all physical parameterization schemes, except the Goddard, present reasonably accurate values for both clear-sky and partly-cloudy conditions.
Clear-sky conditions are relatively easier to predict resulting in lower errors than partly-cloudy conditions.
This study shows that the state-of-the-art parameterization schemes in the WRF model still have scope for improvement.
As a future study, the cloud and dust sensitivity analyses may be performed.