Abstract

The amount of solar insolation that reaches the Earth in one hour is sufficient to fulfill its annual energy budget. One of the challenges for harvesting this energy is due to a lack of relevant data. In the least developed countries like Nepal, the number of observation stations is insufficient. This data gap can be filled by employing credible empirical models to estimate solar insolation in regions where insolation measurements are not available. In this paper, Angstrom–Prescott model parameters are estimated for fifteen different locations of Nepal. Then, correlation is developed for the prediction of solar insolation using only sunshine hour data. The different statistical parameters such as root mean square error (RMSE = 1.958), mean bias error (MBE = −0.018), mean percentage error (MPE = 2.973), coefficient of residual mass (CRM = 0.001), and correlation coefficient (r = 0.909) were used to validate the developed coefficients. The resulting Angstrom–Prescott coefficients are a = 0.239 and b = 0.508. These coefficients can be utilized for the prediction of solar energy at different parts of the country in similar weather conditions.

1. Introduction

The knowledge of solar radiation is important for many applications such as solar power plants, engineering designs, building energy systems, irrigation system development, climatological studies, solar energy systems, evapotranspiration estimation, and regional crop growth modeling [1, 2]. In all these field studies, reliable solar radiation data is vital for a good result.

Nepal earns ample solar radiation from the sun throughout the country as it lies in the most favorable latitude (15°–35°) on the global map. The average global solar radiation in Nepal varies from 3.6 to 6.2 kWhr/m2/day. Likewise, in 2003, there are about 300 sunny days [3] and the annual average sunshine hour and solar energy are about 6.8 hours per day and 4.7 kWh/m2/day in Nepal, respectively [4]. Later in 2010, it was found that the annual average solar insolation is 4.23 kWhr/m2 in Nepal [5] which is very high in comparison with many European countries [6]. This clean energy not only is used for rural electrification but also helps to improve the quality of education, public health, and small-scale cottage industries. That helps to reduce import in LPG gases and petroleum products. In the end, it supports the growth of the economy of the nation. In addition, solar energy is clean renewable energy which helps to decrease the effect of pollution due to the use of clean energy in place of a large amount of traditional energy. In recent years, more than 70 MW of solar PV electricity has been generated in different parts of the country. This solar energy has been utilized for rural electrification and also supply to national grid (30.14 MW) which helps to solve the energy crisis in the nation in a small step. The above mentioned data shows that there is a large amount of solar energy potential available to be produced at different parts of the country to solve the basic need for energy for the holistic development of the country [7, 8].

It is crucial to estimate global solar radiation (GSR) for development of the solar power projects. The information about the GSR can be acquired by installing the pyranometer at different locations. However, it is quite expensive and lengthy process. In Nepal, there are very few stations where a pyranometer is installed by the Department of Hydrology and Meteorology, Government of Nepal. Due to financial and technical constraints, there is an alternative way available to us, which is to use empirical models, RadEst-based model [9], ANN-based model [10, 11], machine learning [12, 13], etc. The best way is to use mathematical models to estimate solar insolation using sunshine hours [14]. Due to the lack of continuous data for a long time and the few meteorological stations that are available in Nepal, there is a lack of research interest. Even if few research activities are found, there is no continuity. This finding could fill a little hole in the realm of solar energy research.

The bright sunshine hour-based model provides better GSR estimation when compared to temperature and cloud-based model. This is because the amount of GSR reaching the Earth’s surface is closely related to a bright sunshine hour [15]. The virtue of bright sunshine hour is that it is a function of latitude, solar declination, and cloudiness. Using the ratio of the bright sunshine hour to sunshine duration called relative sunshine hour data (n/N) has the advantages of eradicating the first two factors leaving the cloudiness as the only one to be considered. The relative sunshine hour follows the necessity of relative radiation, i.e., the ratio of GSR to extraterrestrial solar radiation [16]. This type of model is found in the Angstrom–Prescott model [17, 18]. Most of the long wavelength solar radiation reaching the Earth surface is absorbed and transmitted by the atmosphere and certain part is reflected to space. The GSR is the total amount of solar insolation reaching the Earth’s surface and has less dependence on temperature [19].

Many authors estimated the value of “a” and “b” for different sites. Black et al. collected records of solar radiation and sunshine duration have been collected for 32 stations and obtained regression constants a = 0.23 and b = 0.48 [20]. Mabasa et al. estimated Angstrom–Prescott regression constants for six locations in South Africa [19]. Srivastava and Panday estimated Angstrom–Prescott regression constants for India using records of seven different locations [21]. Muzathik et al. estimated empirical constants based on the monthly average record for Terengganu state, Malaysia [22]. Similarly, Amorox et al. [23], Janjai and Tohong [24], Chegaar et al. [25], Junliang [26, 27], Ogleman et al. [28], Samuel [29], and Rivington et al. [30] estimated the values of “a” and “b” for different sites at different parts of the world. In the case of Nepal, Poudyal et al. [5] estimated Angstrom–Prescott coefficients (a = 0.21, b = 0.25) for Kathmandu using monthly mean daily GSR data for year 2009. Similarly, Adhikari et al. [31] did the same for four locations, Kathmandu, Biratnagar, Pokhara, and Jumla, from 2011 to 2012, and compared it with temperature and humidity dependent model. Joshi et al. [32] estimated solar insolation (4.16 kWh/m2/day) and empirical constants (a = 0.325, and b = 0.272) for Kathmandu using the A-P model and other meteorological parameters dependent model on the basis of daily mean hourly data of GSR from 2010 through 2015.

The aim of this work is to develop empirical constants for the Angstrom–Prescott model for different locations of Nepal. Then, correlations are developed for the prediction of solar insolation using only sunshine hour. The obtained empirical constants can be used to predict GSR with radially available sunshine hour in coming years for given locations and similar geographical locations. The value of empirical constants obtained from plot can be used to predict GSR for all over the locations of Nepal. The estimated GSR for any location can promote clean and renewable energy technology consequently reducing pollution due to use of large amount of traditional energy resources. Again, it supports the growth of economy of nation by lowering import of LPG gas and petroleum products and reinforces rural electrification. However, this research does not cover other meteorological parameters such as temperature, relative humidity, rainfall, altitude, and dewpoints.

2. Materials and Methods

2.1. Site Selection

Nepal (Lat. 26°22′N–30°27′N and longitude of 80°04′E to 88°12′E) is situated at the complex terrain of the Trans-Himalaya region and is landlocked between India and China. It is about 800 km long and 200 km wide with an area of 147,516 km2. Ecologically, Nepal is divided into three regions: Low-land, Mid-land, and High-land. The geographical locations of the fifteen measuring sites are presented in Figure 1 and Table 1.

The solar insolation and meteorological parameter data were obtained from local government authorities such as Alternate Energy Promotion Centre (AEPC), Government of Nepal, World Bank, and Department of Hydrology and Meteorology (DHM), Government of Nepal, for the year 2018–2020. For five stations (Dharan, Pulchowk, Lumle, Nepalgunj, and Jumla), the data of solar radiation were obtained from stations installed by AEPC in collaboration with World Bank. The data in these stations were logged every minute. For the remaining ten locations, the data of solar radiation were obtained for one hour-interval from Department of Hydrology and Meteorology (DHM), Government of Nepal. These data were converted into daily average global solar radiation using Simpson’s one-third rule implemented using Python programming language. There were negligible missing or unknown data. Such missing and unknown data points were filled with the average of the adjacent values of the same variable. Same was done for the outliers [13]. Kipp and Zonen pyranometer is used to measure solar insolation. Campbell-Stock sunshine recorder is used to measure bright sunshine duration. Sunshine duration is the period during which direct solar irradiance exceeds a threshold value of 120 W/sq m [34].

2.2. Model

There exists linear relationship between GSR and sunshine duration was presented by Angstrom in 1924 [17, 35]. The equation iswhere is monthly average daily global solar radiation measured on horizontal surface, is monthly average clear sky daily global solar radiation measured on horizontal surface, is monthly average daily bright sunshine hours, is monthly average maximum possible sunshine hours, and and are empirical constants.

There is ambiguity in (1) as there is uncertainty in the definition of and because there may be a problem in calculating clear sky radiation accurately. This model was modified by Prescott in 1940 [18] by replacing by , extraterrestrial solar radiation. The new model is known as Angstrom–Prescott model which is the most popular and commonly used model, given bywhere is monthly average daily extraterrestrial solar insolation for the location and and are empirical constants.

The empirical constants and depend upon location. The coefficients and represent the fraction of extraterrestrial radiation on overcast days and average days, respectively. The ratio Hg/Ho is the clearness index and n/N is the cloudless index. It gives information about the atmospheric characteristics and conditions of the study area. In this way, the empirical relations can be used to generate solar radiation for the implementation of solar energy and solar thermal technologies where there are no alternative means of energy in all parts of the world [36].

The daily extraterrestrial solar radiation on a horizontal surface (Ho) in MJ/m2/day is computed from the following equations [35, 36]:where Isc is the solar constant (=1367 W m−2), is the latitude of the site (rad), δ is the solar declination (rad), ω is the mean sunrise hour angle for the given month, and nd is the number of days of the year starting from the 1st of January (nd = 1) to 31st December (nd = 365).

The solar declination (δ) and the mean sunrise hour angle (ω) can be computed by the following equations [35, 36]:

The maximum possible sunshine duration (day length) in hours can be computed by [35, 36]

The clearness index (KT) is the ratio of the measured horizontal solar insolation () to the extraterrestrial solar radiation.

The daily extraterrestrial solar insolation and day length N for 15 meteorological stations each are mentioned above using (3) and (6), respectively. The data of solar insolation and bright sunshine hour of each station over the period 2018–2020 are analyzed and prepared in the form of hourly averaged daily solar insolation. These data were used in (7) to estimate daily solar insolation () on the surface, using regression techniques.

is the hourly average daily measured solar radiation; is daily extraterrestrial solar insolation for the location; n is the daily bright sunshine hour. The regression coefficients a and b of the Angstrom–Prescott model are the intercept on axes and slope of regression line, respectively. The validation of estimated solar insolation is done by comparing estimated annual solar insolation with measured solar radiation.

2.3. Statistical Analysis

The statistical tools used to validate estimated data are room mean square error (RMSE), mean bias error (MBE), mean percentage error (MPE), coefficient of residual mass (CRM), and correlation coefficient (r) [37]:where is measured GSR, is estimated GSR, and is number of data points.

Then the yearly average of correlation between estimated and measured values of solar insolation is determined. The correlation between Angstrom–Prescott coefficients and sunshine duration ratio for all given places are determined using trend lines in the plot between them. With the help of a correlation equation, which uses only sunshine hour data, it will be possible to estimate solar insolation for the places where solar radiation data is not available. The proposed methodology is expressed as flow chart in Figure 2.

3. Results and Discussion

The extraterrestrial daily solar radiation and maximum day length are calculated by using (3) and (6), respectively, for all locations separately. The empirical constants also called model coefficients a and b are determined by linear regression analysis technique for different locations for the years 2018 and 2020 on average. The regression coefficient ‘a’ is the intercept and b is the slope of the regression line. Regression equations obtained for different meteorological stations are listed in Table 1 for daily solar insolation and sunshine data. These equations are used for the estimation of daily GSR for different years at different geographical locations separately. The estimated GSR was compared with corresponding values of measured GSR. The validation of a model for different locations is performed by different statistical tests such as correlation coefficient, root mean square error (RMSE), mean bias error (MBE), and mean percentage error (MPE).

The linear variations of daily average hourly measured and estimated GSR for the model for all locations are shown in Figure 3. Similarly, daily variations of measured and estimated GSR for given locations for the models are shown in Figure 4. Figure 3 shows that there is remarkable agreement between daily average measured and estimated solar insolation with a highly acceptable coefficient of determination (R2) greater than 0.705 except Okhaldhunga. Jumla has the highest R2 equal to 0.906. Likewise, Figure 4 shows that there is good agreement between daily average measured and estimated GSR. The annual average highest value of GSR is 19.11 MJ/m2/day for Jumla and that lowest value is 13.06 MJ/m2/day for Ilam. This can also be observed in Figure 5 that gives the variation of GSR, clearness index, and relative sunshine hour. Generally, with the increase of clearness index and relative sunshine hour, GSR is increased following all stations except for Kankai. Here, even at large relative sunshine hour, clearness index is low since the sky is cloudy and dusty. In Figure 4, it is shown that at the time of monsoon season (June, July, August, and September) GSR decreases, due to increase in cloudy and rainy days. This effect occurs highest in Lumle, Pokhara, Ilam, and Jiri. However, this effect is lowest in Nepalgunj and Dharchula. From Figure 6, it is observed that Jumla has maximum clear days (189 days) (clearness index >0.65) and 25 cloudy (clearness index <0.34) days as there is less pollution and it is far from the urban areas [38]. Similarly, in eastern part of Nepal, Ilam has 141 (maximum) days, where 60 days are cloudy and clear days, respectively.

Furthermore, it is validated by comparing estimated annual GSR for the same meteorological locations and measure annual GSR for that location. Table 2 gives the estimates of annual solar insolation for selected meteorological stations in Nepal with measured annual solar insolation. It shows that there is good agreement between measured and estimated values. The relative percentage change between measured and estimated GSR varies numerically from 0.117% (Lumle) to 1.136% (Dhunibesi). This validation can be seen in Figure 7(a) and Figure 7(b). There is a slight increase in GSR with altitude except for Dhunibesi and Ilam due to local weather conditions.

The values of statistical errors, RMSE, MBE, MPE, CRM, and r, for fifteen meteorological stations are listed in Table 3. These errors give performance of the Angstrom–Prescott model in different locations [12]. The correlation coefficient lies between 0.755 and 0.972 and highest values occurred for station at Jumla and lowest at Okhaldhunga. Similarly, the coefficients of determination (R2) are also highest at Jumla as shown in Figure 3. The lower values of RMSE, MBE, MPE, and CRM are preferred as they indicated differences between estimated and measured value of solar radiation. The lowest value of RMSE is 1.344 MJ/m2 /day for the station at Nepalgunj. Similarly, the lowest values of MBE and MPE are 0.005 MJ/m2/day and 0.645% for the stations at Kankai and Musikot, respectively. Lower values of MBE and MPE indicate goodness of fit between clearness index and relative sunshine hour. Also, the value of MPE is less than 7.746% at any stations considered. The CRM ranges from 0 to 0.01 which indicates perfect estimation. This suggested that the Angstrom–Prescott correlation model is good model to estimate the solar radiation in Nepal.

Now from table, the sum of regression coefficients (a + b) is the transmissivity of the atmosphere for solar insolation under perfectly clear sky condition. The clear sky (day) means n/N = 1; then equation (7) becomes

For complete overcast day, n/N = 0, then (7) is

Thus, the empirical constant “a” can be interpreted as transmissivity of an overcast (day) atmosphere [39]. At this time insolation is due to different components. The values of the sum of the empirical constants (a + b) representing the max clearness index (for n/N = 1) are found to be almost equal for all stations. The average value of (a + b) is 0.747. The highest value of (a + b) is 0.833 found at Jumla. For Darchula, Okaldhunga, Lumle, Jiri, and Jumla, the value of (a + b) is greater than 0.802 indicating clear sky due to less pollution at those places. In the paper of Martinez-Lozano, they found parameters “a” and “b” ranging from (0.016 to 0.44) and (0.19 to 0.87), respectively, for 101 locations using monthly average data. Similarly, they found parameters “a” and “b” ranging from (0.19 to 0.36) and (0.43 to 0.62), respectively, for 57 stations using daily data. Thus, the obtained values of regression coefficients are in close agreement with this paper [16]. The observed empirical constant a is in close agreement with value observed in paper of Poudyal (a = 0.21, b = 0.26) for Kathmandu on the basis of data of 2010 [6] but second constant is different. Table 4 lists the empirical constants of different locations to compare with literature obtained values. It is observed that the empirical constant a is almost same for all mentioned and literature locations. There are fluctuations of value of b for Indian cities.

Now the averages of the correlation obtained in Table 3 for all meteorological locations are a = 0.239 and b = 0.508. Then, we get Angstrom–Prescott model for Nepal:

From this relation, solar insolation can be found if sunshine duration is known at any part of Nepal. Also, the new correlation between “a” and “n/N” and “b” and “n/N” is shown in Figure 8 and Figure 9. On fitting curve, the empirical constants “a” and “b” fit into second-order polynomials with the values of R2 greater than 64%. Obtained values of coefficient of determination were considerably high and thus fitted quadratic equation for calculating the values of empirical constants for a given location can be justified. These values of empirical constants are extensively used for convenient calculation of solar insolation.

The solar insolation may thus be estimated for any site in Nepal using (16a), (16b), and (7), even when solar radiation measurements are not available. Sunshine hour data suffices for the calculation of coefficients of the Angstrom–Prescott model.

The solar insolation is highest in the months of May, June, July, and August (about 140 days). Jumla, Lumle, Nepalgunj, Kathmandu, and Ilam have a high effect of rainfall during monsoon causing a reduction in solar insolation. However, at Dharan and Tikapur, this effect is comparatively less. The average yearly solar insolation is maximum (19.113 MJ/m2/day) at Jumla and minimum (13.064 MJ/m2/day) at Ilam. At Jumla, the total solar radiation was 6976.32 MJ/m2/day and at Ilam was 4768 MJ/m2/day. Furthermore, solar insolation in western Nepal represented by Jumla was 5.309 kWh/m2, in mid-Nepal represented by Pulchowk was 4.499 kWh/m2, and in eastern Nepal represented by Ilam was 3.63 kWh/m2 [4]. It happened due to precipitation trends and local weather conditions. At the same time, it is noted that there is gradually lesser precipitation from the eastern part to the western part of Nepal, except in Pokhara. The annual average solar insolation of 4.31 kWh/m2/day is found which is slightly higher than Poudyal, 2015, due to the lockdown effect of COVID-19 [6, 42].

At the end, this type of research work of prediction of empirical constants to find the solar energy is novel work in our complex terrain of Himalaya. This location is not only vulnerable in terms of climate change, landslide, floods, fast rate of snow melt, and changes at biodiversity but also geographically very young mountain and still rising. So, this type of study is essential not only to promote carbon zero emission energy resources but also to solve the energy crisis at a local as well global scale [43].

4. Conclusion

Solar energy is one of the most effective and economical alternative energy sources. Estimation of solar insolation is essential for designing and sizing the solar energy system. In this study, regression technique was used to calculate annual average Angstrom–Prescott coefficients. The empirical constants were found to be a = 0.230 and b = 0.508, respectively, for Nepal and the annual solar insolation was 4.31 kWh/m2/day. The statistical analysis confirmed that there is a good harmony between measured and estimated solar insolation. In this result, second-order polynomial equations based on the relative sunshine hour had been obtained for each of the empirical constants. The empirical constants and equations developed in this study might be used to calculate solar insolation where sunshine duration values are readily available. The outcome of this research is supportive to make plans, policies, and programs to promote clean and renewable energy technology in Nepal. Lastly, the sunshine-based model is best for majority of the study sites and in order to account for complexity of terrain further meteorological parameters might need to be included in some cases.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to express their sincere thanks to the Department of Hydrology and Meteorology (DHM), Government of Nepal, and Alternate Energy Promotion Center (AEPC), Government of Nepal, for providing relevant meteorological data and solar insolation data. The authors would like to give special thanks to Rajesh K. Bachchan, Bishnu Maharjan, faculty members and staff of the Dept. of Physics, Patan Multiple Campus, TU, Nepal, and faculty members of Applied Sciences, IOE, Pulchowk Campus, TU. The author would like to thank NAST for providing partial financial support to forward Ph.D. research work.