Abstract

This study mapped the solar radiation received by slopes for all of Korea, including areas that are not measured by ground station measurements, through using satellites and topographical data. When estimating insolation with satellite, we used a physical model to measure the amount of hourly based solar surface insolation. Furthermore, we also considered the effects of topography using the Global Land One-Kilometer Base Elevation (GLOBE) digital elevation model (DEM) for the actual amount of incident solar radiation according to solar geometry. The surface insolation mapping, by integrating a physical model with the Communication, Ocean, and Meteorological Satellite (COMS) Meteorological Imager (MI) image, was performed through a comparative analysis with ground-based observation data (pyranometer). Original and topographically corrected solar radiation maps were created and their characteristics analyzed. Both the original and the topographically corrected solar energy resource maps captured the temporal variations in atmospheric conditions, such as the movement of seasonal rain fronts during summer. In contrast, although the original solar radiation map had a low insolation value over mountain areas with a high rate of cloudiness, the topographically corrected solar radiation map provided a better description of the actual surface geometric characteristics.

1. Introduction

Renewable energy has attracted considerable attention as we attempt to meet the dual goals of reducing greenhouse gas emissions and securing energy resources, while limiting future extreme weather and climate impacts, and ensuring reliable, timely, and cost-efficient delivery of energy [1]. Among renewable energy sources such as wind and solar energies, solar energy is one of the oldest, cleanest, and most reliable renewable energy source. In the case of solar radiation, there are many different applications: analysis of the thermal load on buildings, solar energy collecting systems, crop growth models, atmospheric energy balance studies, and so forth [2], with the result that the solar radiation levels form the most important parameter for solar energy applications [3]. Therefore, an accurate model predicting the levels of solar irradiation is required [4].

The incident of solar radiation on the Earth’s surface has high temporal and spatial variations, affected by the environmental system, indicating that the areas in which solar energy capture is efficient must be determined. Although pyranometer measurements could accurately determine the incident solar radiation in the desired locations, the number of such stations is usually insufficient to provide solar radiation data for the desired areas; the spatial representativeness is low [5]. An alternative solution to this problem is to use a physical model to derive direct normal irradiation, by integrating satellite data [6]. Use of satellite data has an advantage over radiation data from ground-based measurements in terms of better temporal and spatial coverage [7]. Zhang et al. [8] presented a practical scheme to create the Global LA Surface Satellite (GLASS) incident shortwave radiation (ISR) and photosynthetically active radiation (PAR) products. These were developed using a look-up table (LUT) containing data from polar and geostationary satellites for global coverage. He et al. [9] and Bisht and Bras [10] developed a physics-based algorithm for estimating incident shortwave radiation using satellite data, demonstrating its usefulness in a separate but related context. The use of these data, in conjunction with quality solar ground datasets and other meteorological data, has become an effective way of developing site-time specific solar resource assessment over large areas. In particular, deriving solar radiation data from geostationary satellite images is advantageous for solar resource characterization over large areas [11]. Satellite imagery is useful for monitoring the complex variations of cloud because clouds are the most important atmospheric phenomena influencing the spatial and temporal distribution of solar radiation [12]. In particular, geostationary satellites provide information with high spatial and temporal resolution about the Earth’s atmosphere and cloud cover over a large area [11].

However, the amount of solar energy incident on the ground strongly depends on the sky condition. Larger amounts of radiation are received at ground level when the sky is free of clouds. Moreover, when clouds are present, incident radiation depends on cloud type. This being so, most of the solar insolation algorithm models (e.g., Festa and Ratto [13] and Badescu [14]) consider the state of the sky [15]. The amount of cloud in the same time zone differs regionally, and due to the recent climate change, the incidence of clouds is changing regionally. Therefore, the importance of satellites, which can observe the characteristics of these complex clouds in real time, has been emphasized. Recent studies integrating a physical model and geostationary data showed reasonable results when compared with ground station measurements and revealed high temporally and spatially insolation [16, 17]. However, the problems of existing studies are that they did not take topography elements, such as slope and aspect, into consideration. It is important to consider topographic effects when calculating the amount of solar radiation because the amount of solar radiation actually arriving at a particular point is called global (or total) solar radiation and depends mainly on cloudiness, time of year, latitude, and surface geometry [18].

In this study, solar radiation received in slopes maps was created using the Communication, Ocean, and Meteorological Satellite (COMS) data and calculated based on a physical model. Additionally, to reflect actual surface geometry topographic information was also considered for the efficiency of solar insolation on the Korean peninsula. An accurate solar radiation from satellite would be helpful for estimating solar energy resource maps and deciding the location of solar energy facilities and would enable analysis of suitable land to supply energy to poor areas where electrical equipment is insufficient.

2. Materials and Methods

2.1. Data Used

COMS is the Republic of Korea’s first geostationary orbit composite satellite to perform the duties of ocean and meteorological observation and communication services. This satellite provides observation of the Korean Peninsula and the larger Asia Pacific region at up to 8 min intervals for 24 h and is a significant contributor to weather forecasts and early warnings of high-impact weather events, such as typhoons and torrential rains. The COMS Meteorological Imager (MI) observes one visible band and four infrared bands. Its spatial resolution is 1 × 1 km for the visible area, and it obtains infrared information from a 4 × 4 km area per pixel (Table 1). The data used were from 1 January 2012 to 31 December 2012 over the Korean peninsula. The study area was defined as the southern parts of the Korean peninsula from 34.00°N to 38.41°N and 125.36°E to 129.77°E with a 1 km spatial resolution using the UTM (Universal Transverse Mercator) World Geodetic System (WGS) 84. For physical observations made from the COMS satellite, we used the digital number (DN) to radiance and DN to brightness temperature conversion tables from the National Meteorological Satellite Center, Korea Meteorological Agency (http://nmsc.kma.go.kr/).

In this study, 37 ground station pyranometers of the CM21 model produced by Kipp and Zonen, established in the Korean peninsula, were used as validation data for accuracy assessment of satellite-based insolation (Figure 1). A detailed specification of the CM21 pyranometer, including its associated uncertainties, is presented in Table 2 [19].

The pyranometer data distributed by the Korea Meteorological Agency (KMA) provided hourly measurements based on an upper-hemisphere observation system.

In this study, to consider topography effects, the Global Land One-Kilometer Base Elevation (GLOBE) was used, which is an internationally designed, developed, and independently peer-reviewed global elevation model (DEM), with a latitude-longitude grid spacing of 30 arc-seconds [20]. Figure 2(a) shows GLOBE images over the study area, geometrically resampled to match the defined projection. It shows that the mountain range stretches from the northeast to the south of the Korean peninsula. Based on the DEM of GLOBE, the corresponding slope and aspect were retrieved by using the Environment for Visualizing Images (ENVI) version 5.0 topographic modeling function (Figures 2(b)-2(c)).

2.2. Insolation from a Physical Model with Topography

In contrast to previous empirical models, satellite-based physical insolation models utilize instantaneous satellite observation and atmospheric information [21]. To simplify the estimation of solar radiation using satellite data, we assumed that the hemispherical angular integral of the incident solar radiation consisted mainly of the direct solar radiation and diffuse radiation attenuated by the atmosphere. In other words, we did not consider the radiation effects of atmospheric constituents and cloud from hemispherical sphere due to the complicated physical characteristics and time-consuming calculations involved. Instead, for simplification, we adopted a pixel-based physical model to estimate the solar radiation. Considering the topographical effects of solar radiation, satellite-based incident radiation on a topographically effective pixel was also used to estimate the solar radiation received by slopes corresponding to the slope and aspect pixel values.

Atmospheric constituents such as Rayleigh scattering, aerosols, water vapor, and total ozone are used as input parameters to simulate the atmospheric effects due to absorption, scattering, and reflection. When simulating the satellite-based physical insolation with COMS MI data, we separated the atmospheric constituent products as input for the physical model since the improved atmospheric parameter could easily be replaced to obtain better results with the development of satellite sensors and ancillary data. COMS MI observation is mainly used for interpreting the effects of attenuation by cloud. Cloud is the main obstacle to acquiring reliable incident solar radiation on the surface, due to the high sensitivity of its temporal and spatial variations.

The suggested satellite-based insolation is based on the Kawamura physical model [21] and we improved its cloud factor by considering the satellite visible reflectance and solar zenith angle instead of using brightness temperature, because the pass depth of cloud is more sensitive to the amount of irradiance attenuation [22]. When determining the satellite-based insolation physical model over the study area, we placed a high priority on product consecution by using different satellites. The algorithms for solar insolation obtained with the geostationary satellites GMS-5 [21], MTSAT-1R [16], and COMS MI [23] were based on the Kawamura physical model and improved according to the characteristics of the onboard sensors. In this study, for consecution of the insolation over the study area, the Kawamura physical model was used. The details of the physical model are as follows [16, 17, 2125]: where ,  ,  , and are the total insolation, direct irradiance, diffuse irradiance due to Rayleigh scattering, and diffuse irradiance due to scattering by aerosols, respectively. Detailed descriptions of the parameters used for the physical model of satellite-based insolation are presented in Table 3 [22].

When estimating insolation with satellite data, atmospheric variables are crucial parameters, the validity of which will determine the accuracy of the final products. We studied the variations in attenuation by atmospheric constituents such as ozone, aerosol, water vapor, and cloud. Before considering atmospheric attenuation along each of these parameters, the relative air mass was estimated (Table 3). Relative air mass is used in most atmospheric parameterization equations, because the level of attenuation due to existing constituents (ozone, aerosols, Rayleigh scattering, and water vapor) strongly depends on the relative optical depth determined by the sensor-target-solar geometry [22]. In the case of transmittance due to absorption by ozone, the following equation is used:where is the transmittance due to absorption by ozone, is the amount of ozone, and is relative air mass. The total ozone amount is derived from Ozone Monitoring Instrument (OMI)/Total Ozone Measurement Satellite (TOMS) data from the Earth Observing System Data and Information System, in the Distributed Active Archive Center at the Goddard Space Flight Center. In this study, the daily total ozone amount was used because the sensitivity of SSI to the total quantity of ozone is low, especially when the solar zenith angle is high [16, 26].

For the transmittance due to attenuation by aerosols, the method developed in Frouin et al. [27] was used:where is visibility. In this study, we used Moderate Resolution Imaging Spectroradiometer (MODIS) aerosol optical depth (AOD) products for inputting the visibility by resampling its coarse 10 km spatial resolution for the corresponding projection. When no MODIS products were available, visibility was kept at a fixed value of 20 km indicating nominal weather conditions in South Korea.

In the case of water vapor absorption, the absorptance of water vapor, , was estimated using precipitable water, , which is derived using COMS MI infrared bands. To calculate absorptance, we use the equations from Lacis and Hansen [28] and Chesters et al. [29]: where and are the absorptances of water vapor and precipitable water, respectively. is brightness temperature of COMS MI at 10.2 μm, and is brightness temperature of COMS MI at 12.5 μm. is the air temperature, which is estimated as from Chesters et al. [29].

In this study, for actual solar energy resource mapping of the study area, topographic information was considered using GLOBE topography data. To apply topographic effects of the study areas to COMS-based insolation, the following equations were used [3032]:where and are the slope angle and azimuth of slope, respectively. For direct irradiance, (5) was applied to consider the effect of slope and aspect. For diffuse irradiance on slope, assuming isotropic sky radiation, diffuse insolation was retrieved using (6). The process of physical modeling involved mainly separating images according to cloudiness-corresponding pixels. In the case of clear sky conditions in (7), total insolation, which consists of and , was used as the final insolation product. In cloudy conditions, cloud factor was applied to total insolation using (8):where cloud factor is the ratio of attenuation by cloud for incident solar radiation using COMS visible band and solar zenith angle [22]. The COMS visible band was used to indirectly determine the thickness of the target clouds, where a higher cloud reflectance indicates a lower cloud penetration value. Additionally, the solar zenith angle over the cloudy area is used to determine cloud effects, because the attenuation of incident solar radiation by clouds depends on the pass depth of the clouds [16]. To determine the cloud factor, reference match-up datasets between satellite observations (top of the atmosphere reflectance and solar zenith angle) and cloudy pyranometer measurements were first analyzed and the cloud attenuation was estimated as a function of cloudy reflectance and solar zenith angle by using a simple look-up table [22]. Table 4 shows the simplified cloud LUT for solar radiation with COMS MI.

3. Results

To compare satellite-based insolation with ground station-based pyranometer, temporally and spatially systematic differences were considered to match up datasets. From the point of view of temporal discrepancy, satellite-based insolation was estimated by using instantaneous observed images, while incident irradiances observed by pyranometer measurements every 2 min were integrated to estimate an hourly based ground measurement. By adjusting the hour factor in the physical model, the temporal discrepancy was reduced when establishing matched-up data [33]. For spatial discrepancies induced by systematic differences between pixel-based satellite insolation and hemisphere upward-looking based pyranometer measurement, a representative test on the pyranometer measurement of insolation through satellite observation was performed to determine the spatial representativeness of the ground stations. In this study, we averaged the satellite-based insolation using windows of sizes 1 × 1, 3 × 3, and 5 × 5 by comparing the calculated insolation values with pyranometer measurements for both clear and cloudy sky conditions. The results of the comparison are shown in Figure 3. By using averaged values for each window size, we mitigated the effects of possible cloud drift between adjacent pixels during the hour-long averaging interval [34]. In Figure 3, scatterplots show comparisons of the satellite-based measurements of insolation with ground measurements for window sizes from 1 × 1 to 5 × 5 for clear (a, c, e) and cloudy (b, d, f) conditions. The relationship between satellite-based measurements of insolation and pyranometer measurement appears to be linear, with most of the data points located close to the reference line , though many data points deviate significantly from this line. In the case of the clear sky conditions shown in the (a, c, e) plots of Figure 3, the estimated satellite-based insolation was more accurate than it was for cloudy conditions because we did not have to take the effects of clouds into account. The statistical values were also more accurate in this case, as the window size increased; we obtained a root mean square error (RMSE) value of 72.49 W/m2, a mean bias error (MBE) of 6.64 W/m2, and a value for (determination coefficients) of 0.903. Since the effect of cloud is to attenuate incident solar radiation, black dots show areas of low insolation under cloudy conditions (b, d, f). The statistical values for the RMSE, MBE, and decreased to 106.94 W/m2, −4.99 W/m2, and 0.712, respectively. Generally, as collocated window sizes increased, RMSE and MBE decreased for both sky conditions. As window sizes increased beyond 3 × 3, the accuracy of satellite-based insolation values stabilized. Therefore, in this study, 5 × 5 widow sizes were determined to be adequately representative of the ground stations. Based on the comparisons with the pyranometer, an accuracy assessment for the physical model with COMS MI images was performed and showed reliable accuracy for mapping of solar energy resources over the study area. In this study, we did not compare topographically corrected satellite-based insolation with pyranometer data for accuracy assessment, because all of the pyranometer instruments were carefully installed on a horizontal plane, indicating that they did not reflect the topography effect. That is to say, most of the measurement sites were established based on World Meteorological Organization (WMO) criteria (Guide to Meteorological Instruments and Methods of Observation WMO-No. 8), which meets requirements for reducing adjacent effects.

Before analyzing the spatial characteristics of the estimated insolation, we estimated the sensitivity of the results to variation in the values for water vapor, visibility, and total ozone input as atmospheric variables into the physical model. Figure 4 shows contour diagrams of the satellite-based insolation as a function of water vapor, visibility, and total ozone. As shown in Figure 4(a), when comparing water vapor with visibility, insolation was more sensitive to water vapor levels. In the bottom contour diagram in Figure 4(a), the physical model was more sensitive to the water vapor content than visibility when the visibilities were greater than 10 km. Ozone levels vary more than water vapor, but the model is less sensitive to ozone variation. When comparing visibility to ozone, visibility was the dominant factor, with the model indicating the highest sensitivity to visibilities of less than 10 km. We concluded that the simulation model was most sensitive to water vapor, but visibility was also a significant consideration in the case of thick aerosol optical depth.

Lastly, we produced an integrated monthly and yearly solar radiation map from COMS hourly based insolation data for 2012 over the Korean peninsula. The original solar radiation monthly maps, without considering topography, are presented in Figure 5. As seen in the figures, temporal variations are well described, showing lower values (blue color) close to winter and higher insolation (red color) close to summer. Because integrating the physical model with the COMS MI image is designed to reflect the temporal variation of atmospheric conditions, especially for cloud, the detailed patterns of solar radiation map are well represented. For example, the solar radiance map for the south and east parts of the peninsula had lower values (yellow color) in June 2012 despite it being summer, which has the highest extra solar irradiance. In July, the solar radiance map shows a low value in the north (green color) and central (yellow color) areas. Those lower areas during the June and July period, compared with May, mainly originate from the rainy monsoon season, which is called Jang-Ma in the Republic of Korea. According to KMA weather reports (http://www.climate.go.kr/), Jang-Ma started on 18 June 2012 and lasted until 17 July 2012 and moved from south to north. This means that the regional atmospheric conditions, such as the movement of the seasonal rain front, are well captured in the solar energy resource map. Furthermore, spatial variations are also well visualized in the original solar radiance map. For the whole maps, because the high mountain range area has severe weather conditions, meaning a high rate of cloudiness, constantly low value areas are well matched to the high elevation of the mountain ranges. The northeast coast area, close to the highest mountain range stretching from north to south, has reverse solar radiance patterns between February and August. That abnormal pattern is inferred by the seasonal variation of foehn wind, a dry and hot downward wind that occurs on the leeside of a mountain range. The area opposite (upwind) to the mountain site of the foehn wind appears to be cloudy, which reduces insolation. Generally, since periodic winds blowing in February and August occur from opposite sides, the foehn wind shows similar seasonal variation. Therefore, in Figure 5, the solar radiance map of the northeast shows the opposite pattern for the mountain range between February and August.

The topographically corrected solar radiance maps are presented in Figure 6. The most effective influence on application of topography information on solar radiance maps is the detailed surface geometry of the terrain. In the case of the original solar energy resource maps in Figure 5, because their spatial resolution was 4 km from use of the COMS MI image a coarse spatial pattern is shown, while the topographically corrected solar radiance map with a 1 km resolution from the GLOBE DEM captures the detailed terrain effects and resembles a high resolution map. The overall temporal characteristics of topographically corrected solar energy resource maps tend to not only show the similar temporal pattern of the original solar radiance maps but also well describe the topography of the study area under benign weather conditions. For example, the topographically corrected solar energy resource map shows a better result with regard to the spatial characteristics of foehn wind, showing more detailed reverse patterns in the northeast part of the mountain range. However, the topographically corrected solar radiance map taken in the summer season indicates reduced insolation induced by the rainy monsoon, which can be found as coarse spatial resolution similar to the original map. In other words, the detailed surface geometry of the topographically corrected solar radiance is faint during the summer monsoon season. This happens because cloud diffuses incident solar radiation, which contributes overall to the reduction of the topographical effect in our physical model.

Table 5 compares the monthly ground measurements with the original and topographically corrected COMS MI-based solar radiation for clear and cloudy conditions. According to the results from Table 5, the monthly solar radiation for the clear sky condition was more accurate than for cloudy conditions, similar to previous scatterplot results. In terms of the seasonal variation in the model accuracy, summer, which is the season with the highest rate of cloudy conditions, has the lowest accuracy with a 0.340 MJ/m2 RMSE in July. The seasonal accuracy of solar radiation under cloudy conditions was similar, with a high RMSE during summer. The original solar radiation from COMS MI seemed to be more accurate than the topographically corrected solar radiation since most of the selected manned ground measurement sites were established based on WMO criteria. That is to say, pyranometer measurements did not reflect the topographical effects for each selected ground site.

In this study, we also produced yearly solar radiation maps for the study area. As seen in Figure 7, the most outstanding difference between the original and topographically corrected solar radiation maps was the spatial pattern we expected. Similar to the monthly map, DEM was an effective parameter in terms of determining the value of insolation. Although the original solar radiation map did not consider topography, the location of mountains having severe weather conditions could have been included (Figure 7(a)). Therefore, we can conclude that the topographically corrected solar radiation map not only reflects the temporal variations of regional atmospheric conditions but also considers the topographical effect. Based on the suggested solar maps, we can determine the most efficient areas for solar resources in the Korean peninsula; some parts of the western and southeastern regions (shown as blue rectangles in Figure 7(b)) are considered the most efficient solar energy resource area.

4. Summary and Conclusions

In this study, solar radiation received in slopes maps, based on integrating a physical model with COMS MI images, were estimated by additionally considering topography from GLOBE DEM during 2012 for the Republic of Korea. Before presenting a solar radiation map, we assessed the accuracy of the results compared with measurements from 37 ground stations. The temporal and spatial system differences between satellite sensor and pyranometers were carefully considered and showed reasonable accuracy for clear and cloudy sky conditions. In the case of clear sky conditions, a benign weather condition, the model was more accurate than it was in the case of cloudy conditions, yielding an RMSE of 106.94 W/m2 and an MBE of −4.99 W/m2 for 3 × 3 window sizes. The RMSE and MBE in overcast conditions over 3 × 3 windows were 72.49 W/m2 and 6.64 W/m2, respectively, indicating a greater reduction in accuracy than for clear conditions, but reasonable in comparison with previous studies [27, 31, 32]. Based on scatterplots for both sky conditions, we conclude that integrating a physical model with COMS MI image is a reliable method of mapping solar radiation without considering topography. Before analyzing the spatial characteristics of the estimated insolation, we analyzed the sensitivity of atmospheric constituents and showed that the simulation model was most sensitive to water vapor. In this study, topography from the 1 km GLOBE DEM was applied to the original solar radiation map, since the amount of solar radiation actually depends mainly on surface geometry [29]. Monthly and yearly original and topographically corrected solar radiation maps were estimated for 2012. In the case of the monthly maps, both original and corrected maps captured the temporal variations of regional atmospheric conditions such as the movement of rainy monsoon fronts. While, as we expected, the topographically corrected solar radiation map showed more detailed and sophisticated effects of the Korean surface geometry, indicating they could reflect more accurately the real surface terrain. However, during the summer season, which has a high rate of cloudiness, the topographically corrected solar radiation maps did not capture the surface geometry in the same way as the original maps, for two reasons. First, our suggested physical model was not designed to consider downward diffused radiance by cloud on surface geometry. Second, the physical characteristics of complex cloud diffusion are very difficult to parameterize.

In addition, the monthly RMSEs of the original solar radiation seemed to have better accuracy when compared with the topographically corrected values for selected manned ground pyranometers. We inferred that horizontally deployed pyranometers that do not consider environmental topography do not provide good reference data for the data to not reflect the actual topographical effects of solar radiation received by slopes. Finally, the temporally and spatially different patterns found in the yearly solar radiation map are also present in the monthly map. We conclude that our final solar radiation map is useful for determining the areas of greatest solar energy resources, since it not only captures atmospheric variations but also reflects the surface geometry effects over the study domain.

In future works, for a more detailed solar radiation received in slopes map, we plan to perform further analysis by using GIS information based on our final solar energy resource map. To validate the topographically corrected insolation derived from satellite images, we plan to establish pyranometers which would enable us to determine the effects of regional variations in terrain. We plan to install pyranometers in mountainous areas and valleys, anticipating complicated topographical effects in such terrain. We also recognize the importance of solar radiation climatology and plan to study solar radiation climatology when sufficient COMS MI data are secured.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was funded by the Korea Meteorological Administration Research and Development Program under Grant KMIPA 2015-5040.