Several energy management applications rely on accurate photovoltaic generation forecasts. Common metrics like mean absolute error or root-mean-square error, omit error-distribution details needed for stochastic optimization. In addition, several approaches use weather forecasts as inputs without analyzing the source of the prediction error. To overcome this gap, we decompose forecasting into a weather forecast model for environmental parameters such as solar irradiance and temperature and a plant characteristic model that captures site-specific parameters like panel orientation, temperature influence, or regular shading. Satellite-based weather observation serves as an intermediate layer. We analyze the error distribution of the high-resolution rapid-refresh numerical weather prediction model that covers the United States as a black-box model for weather forecasting and train an ensemble of neural networks on historical power output data for the plant characteristic model. Results show mean absolute error increases by 11% and 68% for two selected photovoltaic systems when using weather forecasts instead of satellite-based ground-truth weather observations as a perfect forecast. The generalized hyperbolic and Student's t distributions adequately fit the forecast errors across lead times.
Figures & tables
Figure 1: Inputs/outputs of the plant characteristic model on the time axis.
Figure 2: Plant characteristic model: ensemble of MLP .
Figure 3: Location of selected PV systems and the SOLRAD irradiance observation network.
Figure 4: Pre-processed AC power data from the selected PV systems with training and test data split.
System
33
1423
both
Corr.
P
S
P
S
P
S
GHI
0.91
0.88
0.93
0.97
0.91
0.89
s
0.79
0.86
0.87
0.96
0.81
0.87
Temp.
0.40
0.54
0.23
0.29
0.34
0.34
pt−1
0.90
0.90
0.93
0.93
0.91
0.91
pt−24
0.80
0.89
0.93
0.98
0.82
0.92
Table 1: Pearson (P) and Spearman (S) correlation of input features to PV power output.
Figure 5: Spearman correlation of input features to the AC power output of both PV systems.
Figure 6: Impact of hidden layer width on the R2 score for system 33. The red indicator shows the selected width, highlighting the start of the yellowish plateau.
Figure 7: Errors of the plant characteristic model with a forecast at 23:00 of system 33. Orange curves show the fitted normal distribution in the probability density plots.
Plant
Weather
Combined
System
Power
GHI
Temp
Power
W/kWp
W/m2
\SIUnitSymbolCelsius
W/kWp
MBE
33
-0.12
20.68
3.88
18.34
1423
-5.53
12.66
-1.02
6.57
MAE
33
28.12
39.94
3.97
47.28
1423
28.04
22.74
1.61
31.16
Table 2: Error metrics for 24 leads starting at 6:00 UTC.
Figure 8: Errors of the plant characteristic model with a forecast at 22:00 of system 1423. Orange curves show the fitted normal distribution in the probability density plots.
Figure 9: Errors of the weather forecast model with a forecast at 23:00 of system 33. Orange curves show the fitted normal distribution in the probability density plots.
GHI in W/m2
Temperature in \SIUnitSymbolCelsius
Station
MBE
MAE
MBE
MAE
System 33
20.68
39.94
3.88
3.97
System 1423
12.66
22.74
-1.02
1.61
Albuquerque
11.65
27.89
0.77
1.41
Bismarck
14.00
31.05
0.26
1.34
Hanford
6.04
20.17
-1.88
2.10
Table 3: HRRRv4 forecast accuracy at SOLRAD network stations and selected PV systems with 24 leads at 6:00 UTC.
Figure 10: Errors of the weather forecast model with a forecast at 23:00 of system 1423. Orange curves show the fitted normal distribution in the probability density plots.
Figure 11: Errors of the combined model with a forecast at 23:00 of system 33. Orange curves show the fitted normal distribution in the probability density plots.
Figure 12: Temporal correlations of lead distance using SOLRAD stations locations forecasts for every 6 hours
Figure 13: Errors of the combined model with a forecast at 22:00 of system 1423. Orange curves show the fitted normal distribution in the probability density plots.
Normal
Student’s t
Generalized Hyperbolic
Lead
K-S
C-M
K-S
C-M
K-S
C-M
6
0.615
0.445
0.591
0.419
0.576
0.474
7
0.160
0.092
0.090
0.158
-
0.009
8
0.735
0.866
0.735
0.866
0.597
0.811
9
0.634
0.612
0.590
0.596
0.864
0.845
10
0.000
0.000
0.049
0.104
0.534
0.511
Table 4: Goodness-of-Fit of different distributions for system 1423 showing the p-values of the Kolmogorov-Smirnov (K-S) and Cramér-von Mises (C-M) test at 22:00 local time.
Figure 14: Errors of the weather forecast model for air temperature with a forecast at 6:00 UTC.
Figure 15: Q-Q plots on empirical data against fitted distributions for the forecast at 22:00 for system 1423.
Solar power forecasting predicts photovoltaic generation and irradiance from weather, imagery, plant data, and recent output to support reliable grid and energy decisions.
Short-term photovoltaic power forecasting requires models to represent regular solar cycles and weather-driven fluctuations whose importance changes with the forecast horizon. This study develops a hierarchical ensemble that combines temporal neural models, historical analogs, state climatology, and gradient-boosted trees. Solar geometry and numerical weather forecasts describe the expected generation conditions, while horizon-specific convex weights combine complementary predictions. A separate calibration step uses available historical forecast errors to account for recent bias. The framework is evaluated on public PVDAQ data at 15--240-minute horizons and on three GEFCom2014 solar zones at hourly horizons up to four hours. On PVDAQ, the ensemble achieves a daylight capacity-normalized mean absolute error of 4.315%, reducing error by 4.11% relative to full-feature LightGBM and by 6.03% relative to fine-tuned Chronos-2 under identical calibration. Expert-removal experiments identify redundancy within the ensemble. Across three training seeds on GEFCom2014, learned fusion improves upon equal weighting but performs comparably to LightGBM. The results support horizon-specific combination as a useful forecasting strategy while showing that its advantage over strong individual models depends on the dataset and evaluation period.
Xu Yuqing, Zhou Liguo, Sun Ze +2
School of Computer Science and Technology, Huaibei Normal University, China
An ensemble can improve photovoltaic forecasts while adding components that contribute little or increase computation. We assess these effects through matched comparisons and ablations of a fixed heterogeneous predictor bank. Hourly experiments use GEFCom2014 and three additional public datasets, with chronological partitions and three seeds. Under retrospective ERA5 assistance, static fusion reduces scaled mean absolute error against matched boosting by 1.11%, 4.41%, and 1.63% on PVDAQ, OPSD, and Ausgrid; only OPSD remains supported after multiple-comparison correction. Weather gating offers no consistent incremental benefit. Exploratory member removals show group-level dependence alongside individual redundancy. A separate, previously inspected fifteen-minute case replaces one neural member with a tree predictor: normalized error falls by 1.72%, but measured inference is slower. These findings support component-wise evaluation with explicit limits on weather availability and test-set reuse.
Sun Ze, Zhou Liguo, Xu Yuqing +2
aSchool of Computer Science and Technology, Huaibei Normal University, China
Reliable photovoltaic (PV) forecasts can support low-carbon energy systems, but deployed sites may have only short and incomplete records. Physical and hybrid methods can be sensitive to weather inputs, calibration, and timestamp-alignment, while individual machine learning models may capture different parts of the forecasting problem. We study hourly day-ahead PV forecasting at a United Kingdom charging station using one year of inverter measurements, with 9.25% of hours missing. The pipeline checks timestamp-alignment, derives solar and clearness features, adds short-term weather context, and combines five complementary models using non-negative least squares stacking, with the combination fitted only on validation observations. We compare against smart persistence, a weather-scaled baseline that carries the previous day's PV behaviour forward using target-day irradiance. With retrospective weather, the combined model reduces daylight normalised root mean square error (RMSE) by 31.2% under random day-fold evaluation and by 2.9% under rolling-origin evaluation, although the latter improvement is not robust across days. It also improves by 3.0% over the single model selected from validation performance. Replacing retrospective weather with a public product sampled at a constant 24-hour lead increases daylight RMSE by 13.1% and 4.2% under the two protocols, while retaining positive skill over smart persistence.
Fariba Dehghan, Sebastian Stein, Vahid Yazdanpanah +2