全天空輻射量(global horizontal irradiance, GHI)為光電能源業者用來監測太陽能系統發電效率的主要參考依據，然而GHI容易會受到天氣條件的影響而有明顯的變化，例如：太陽位置、雲和氣膠…等。目前大部分GHI產品的時間解析度為一小時，這樣的時間解析度不足以反應出GHI短時間內快速變化的性質。在本研究中，我們使用向日葵8號衛星(Himawari-8, H8)十分鐘解析度的資料以及應用單層輻射傳送模式反演高時間解析度的GHI資料。本研究的模式中考慮了氣膠、臭氧、水氣、氣體及雲對輻射量的吸收及散射效應，並利用了地面觀測GHI資料及衛星資料透過經驗公式來建立太陽輻射波段及衛星波段之間的關係。同時也利用晴空判斷法及晴空指數將天空條件分為晴空條件、部分晴空條件及多雲條件，並以此改進模式結果之精確度。
本研究選擇中央氣象局嘉義氣象站2018年1月1日至2018年4月12日之資料作為研究站點及研究時間，利用1月1日至2月28日建立太陽輻射波段及衛星波段之間的線性關係，而後反演3月1日至4月12日之地面輻射量，最後與氣象局地面氣候自動觀測系統(Automatic Climate Observation System, ACOS)之地面輻射量觀測資料進行誤差分析並探討模式之反演表現。結果顯示，晴空下的rMBE(相對平均偏差)、rRMSE(相對均方根誤差)及r2 (判定係數)為 -1.7%、6.8%及0.98，非晴空下的rMBE、rRMSE及r2為 -5.9%、22.7%及0.84，模式在晴空下有極佳的表現。在反演策略中，吾人採用將天空條件分為晴空、部分晴空及多雲分別進行模式反演，整體表現的rMBE、rRMSE及r2為 -5.1%、20.5%及0.86，相較於沒有使用天空條件判斷下的結果(rMBE、rRMSE及r2為 -7.9%、22.7%及0.84)，有較佳的表現。透過不同天空條件下的誤差分析，吾人發現衛星反演地面輻射在計算中最大的誤差由雲所造成，可使反演之輻射量的r2值由0.98下降至0.47。
本研究反演之GHI的時間解析度為十分鐘，相較氣象局所提供的一小時時間解析度GHI反演資料，十分鐘時間解析度較能夠模擬出短時間內雲經過時地面輻射量的變化。為了能夠與氣象局衛星反演產品比對，進一步將ACOS地面觀測GHI與本研究反演之GHI的時間解析度(分別為一分鐘及十分鐘)皆平均至一小時，與平均後的ACOS地面觀測進行誤差分析後，結果顯示本研究反演之GHI平均後的rMBE、rRMSE及r2為 -6.2%、14.3%及0.92，而氣象局所衛星反演之GHI產品的rMBE、rRMSE及r2為-1.5%、11.2%及0.94，可以發現模式結果能與氣象局衛星產品相近。進一步評估本反演法的時空適用性，吾人發現距離建立線性關係式的氣象觀測站越遠，模式所反演之輻射量的誤差也會隨之增加，因此若能定義出台灣各個氣象觀測站代表的影響範圍，便能針對不同區域建立不同的線性回歸式以完成更大範圍的地面輻射量反演。

Surface solar radiation (or global horizontal irradiance, GHI) data is critical for photovoltaic and electric power companies to monitor power generation efficiency of photovoltaic system. However, GHI highly varied with sky conditions which mainly due to solar position and clouds variations. The typically temporal resolution of GHI product is 1-hour resolution, which might not be enough to represent the GHI variability. In this study, we used Himawari-8 (H8) satellite 10-min resolution data and applied a one-layer radiation transfer model to derived high temporal resolution GHI data. Our model includes the calculations for the scattering and absorption from aerosol, ozone, water vapor, gases and clouds. The observational GHI data has been used to construct relationships between solar spectrum and satellite bands, and further create an empirical function. Furthermore, a clear-sky identification method with sky index (i.e., clear sky, partly clear sky and cloudy sky) is proposed to evaluate the model derived-GHI performance in different sky conditions.
Data from the CWB Chiayi station was used to evaluate the performance of our model deriving results. Statistical results show that the rMBE, rRMSE, and r2 are -1.7% vs -5.9%, 6.8% vs 22.7% and 0.98 vs 0.84 for under clear sky and unclear sky, respectively. The performance of model under clear sky is outstanding, and the performance of model under unclear sky is moderately, and according to the analysis of the results under different sky conditions, we found that the largest error in the model is caused by the clouds, which could reduce the r2 value of derived GHI from 0.98 to 0.47.
We also found that when we consider sky-index-dependent empirical functions using in the model, the r2 value improved from 0.84 to 0.86. Compared to the CWB 1-hr resolution GHI product (rMBE: -1.5%, rRMSE: 11.2% and r2: 0.94), our high-resolution (10-min computation and average to 1-hr for comparison) show the rMBE, rRMSE and r2 value are -6.2%, 14.3% and 0.92, respectively. In order to evaluate the spatial and temporal applicability of our model, we found that the farther away from the station which used to establish the linear relationships, the error of the derived GHI would increase. With the basement of our model, if we could define the effective radius of weather stations in Taiwan, we can establish different linear regression equations at different places, so that GHI can be derived in more region.

王聖翔 (2020)，台灣地區標準地面輻射觀測網絡系統發展案，交通部中央氣象局一百零九年資訊委託服務計畫期末報告。

ASTM E490-00a (2014), Standard Solar Constant and Zero Air Mass Solar Spectral Irradiance Tables, ASTM International, West Conshohocken, PA, USA.

ASTM G173-03 (2012), Standard Tables for Reference Solar Spectral Irradiances: Direct Normal and Hemispherical on 37° Tilted Surface, ASTM G173-03, ASTM International, West Conshohocken, PA, USA.

Bessho, K., et al. (2016). "An introduction to Himawari-8/9—Japan’s new-generation geostationary meteorological satellites." Journal of the Meteorological Society of Japan. Ser. II 94(2): 151-183.

Bhartia, P. K. (2012). OMI/Aura Ozone (O3) Total Column Daily L2 Global Gridded 0.25 degree x 0.25 degree V3, Goddard Earth Sciences Data and Information Services Center (GES DISC).

Driemel, A., et al. (2018). "Baseline Surface Radiation Network (BSRN): structure and data description (1992–2017)." Earth System Science Data 10(3): 1491-1501.

Dubovik, O., et al. (2000). "Accuracy assessments of aerosol optical properties retrieved from Aerosol Robotic Network (AERONET) Sun and sky radiance measurements." Journal of Geophysical Research: Atmospheres 105(D8): 9791-9806.

Duffie, J. and W. Beckman (1991). "Solar engineering of thermal process." BW Duffie John, Solar Engineering of Thermal Process. 2a Ed., Wiley-Interscience.

Eck, T. F., et al. (1999). "Wavelength dependence of the optical depth of biomass burning, urban, and desert dust aerosols." Journal of Geophysical Research: Atmospheres 104(D24): 31333-31349.

Engerer, N. A., et al. (2017). Himawari-8 enabled real-time distributed PV simulations for distribution networks. 2017 IEEE 44th Photovoltaic Specialist Conference (PVSC), IEEE.

Gueymard, C. (1989). "A two-band model for the calculation of clear sky solar irradiance, illuminance, and photosynthetically active radiation at the earth's surface." Solar energy 43(5): 253-265.

Guyonnet, A., et al. (2019). "Local monitoring of atmospheric transparency from the NASA MERRA-2 global assimilation system." Journal of Astronomical Instrumentation 8(04): 1950013.

He, Q., et al. (2010). "Validation of MODIS derived aerosol optical depth over the Yangtze River Delta in China." Remote Sensing of Environment 114(8): 1649-1661.

Holben, B. N., et al. (1998). "AERONET - A federated instrument network and data archive for aerosol characterization." Remote Sensing of Environment 66(1): 1-16.

Ineichen, P. and R. Perez (2002). "A new airmass independent formulation for the Linke turbidity coefficient." Solar energy 73(3): 151-157.

IPCC (2013). Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, United Kingdom and New York, NY, USA, Cambridge University Press.

Iqbal, M. (1983). "An Introduction to Solar Radiation Academic." New York 19832.

Janjai, S., et al. (2013). "Mapping global solar radiation from long-term satellite data in the tropics using an improved model." International Journal of Photoenergy 2013.

Kasten (1980). "A simple parameterization of the pyrheliometric formula for determining the Linke turbidity factor."

Kasten and A. T. Young (1989). "Revised optical air mass tables and approximation formula." Applied optics 28(22): 4735-4738.

Kuo, P.-H., et al. (2018). "Solar radiation estimation algorithm and field verification in Taiwan." Energies 11(6): 1374.

Lave, M. and J. Kleissl (2013). "Cloud speed impact on solar variability scaling–Application to the wavelet variability model." Solar energy 91: 11-21.

Leckner, B. (1978). "The spectral distribution of solar radiation at the earth's surface—elements of a model." Solar energy 20(2): 143-150.

Letu, H., et al. (2020). "High-resolution retrieval of cloud microphysical properties and surface solar radiation using Himawari-8/AHI next-generation geostationary satellite." Remote Sensing of Environment 239: 111583.

Linke, F. (1992). "Transmission-koeffizient und trubungsfaktor." Beitr. Phys. Atomos. 10: 91-103.

Long, C. N. and T. P. Ackerman (2000). "Identification of clear skies from broadband pyranometer measurements and calculation of downwelling shortwave cloud effects." Journal of Geophysical Research-Atmospheres 105(D12): 15609-15626.

McPeters, R., et al. (2008). "Validation of the Aura Ozone Monitoring Instrument total column ozone product." Journal of Geophysical Research: Atmospheres 113(D15).

Mermoud, A. and B. Wittmer (2014). "PVSYST user's manual." Switzerland, January.

Mira, M., et al. (2015). "The MODIS (collection V006) BRDF/albedo product MCD43D: Temporal course evaluated over agricultural landscape." Remote Sensing of Environment 170: 216-228.

Ohmura, A., et al. (1998). "Baseline Surface Radiation Network (BSRN/WCRP): New precision radiometry for climate research." Bulletin of the American Meteorological Society 79(10): 2115-2136.

Reno, M. J. and C. W. Hansen (2016). "Identification of periods of clear sky irradiance in time series of GHI measurements." Renewable Energy 90: 520-531.

Reno, M. J., et al. (2012). "Global horizontal irradiance clear sky models: Implementation and analysis." SANDIA report SAND2012-2389.

Rienecker, M. M., et al. (2011). "MERRA: NASA’s modern-era retrospective analysis for research and applications." Journal of climate 24(14): 3624-3648.

Rigollier, C., et al. (2004). "The method Heliosat-2 for deriving shortwave solar radiation from satellite images." Solar energy 77(2): 159-169.

Salazar, G. and C. Raichijk (2014). "Evaluation of clear-sky conditions in high altitude sites." Renewable energy 64: 197-202.

Sayer, A. M., et al. (2013). "Validation and uncertainty estimates for MODIS Collection 6 “Deep Blue” aerosol data." Journal of Geophysical Research: Atmospheres 118(14): 7864-7872.

Schaaf, C. and Z. Wang (2015). "MCD43C3 MODIS/Terra+ Aqua BRDF/Albedo Albedo Daily L3 Global 0.05 Deg CMG V006." NASA EOSDIS Land Processes DAAC. doi 10.

Schmid, B., et al. (1999). "Comparison of aerosol optical depth from four solar radiometers during the fall 1997 ARM intensive observation period." Geophysical Research Letters 26(17): 2725-2728.

Tang, W., et al. (2016). "Retrieving high-resolution surface solar radiation with cloud parameters derived by combining MODIS and MTSAT data." Atmospheric Chemistry and Physics 16(4): 2543-2557.

Waliser, D. E., et al. (1996). "An estimate of the surface shortwave cloud forcing over the western Pacific during TOGA COARE." Geophysical research letters 23(5): 519-522.

Wang, K., et al. (2004). "Validation of the MODIS global land surface albedo product using ground measurements in a semidesert region on the Tibetan Plateau." Journal of Geophysical Research: Atmospheres 109(D5).

Yao, C.-M. (2017). Mapping Surface Solar Radiation with Satellite Data over Taiwan, National Central University.