4.1. Overall Accuracy
Take GNSS PWV as the reference, the Bias, MAE, and RMSE [
70] between ERA5 PWV adjusted by models GLPW, GFPW, and GMPW and the unadjusted ERA5 PWV are compiled in
Table 6. The variations of these metrics across different subregions and four seasons are depicted in
Figure 5 and
Figure 6.
As shown in
Table 6, the Bias adjusted by the models GLPW, GFPW, and GMPW are basically close to 0 mm, which shows that the systematic differences between ERA5 PWV and GNSS PWV are basically eliminated. For MAE, the ranges of the GLPW model in the NWC, TP, NC, and SC subregions are 0.68–1.71 mm, 1.07–2.02 mm, 0.74–2.00 mm, and 1.35–2.28 mm, respectively. For MAE, the ranges of the GFPW model in the NWC, TP, NC, and SC subregions are 0.41–0.72 mm, 0.60–0.86 mm, 0.55–1.70 mm, and 0.75–1.03 mm, respectively, while the value range of the GMPW model is 0.02–0.49 mm, 0.03–0.64 mm, 0.18–1.33 mm and 0.10–0.55 mm. Compared with unadjusted ERA5 PWV, the adjusted MAE of GLPW, GFPW, and GMPW models are reduced, and the GMPW model performs better.
As shown in
Figure 5, the adjusted Bias values of the GLPW, GFPW, and GMPW models decrease significantly in each season, indicating that the systematic differences between ERA5 PWV and GNSS PWV have been effectively eliminated. Therefore, the unadjusted ERA5 PWV exhibits significant negative Biases in spring, summer, and fall in the TP and SC subregions. The Biases of the adjusted GFPW and GMPW models generally show positive Biases in the four subregions and four seasons. The Bias of the GMPW model is lower than that of the GFPW model, indicating that the stability of the GMPW model is better. This also shows that the accuracy of ERA5 PWV before adjustment varies in different subregions in each season and is unevenly distributed.
As shown in
Figure 6, the RMSE of unadjusted ERA5 PWV in the NWC subregion compared with GNSS PWV and the RMSE after adjustment by the GLPW model across four seasons is 1.43/0.83 mm (spring), 2.09/1.07 mm (summer), 1.40/0.84 mm (fall) and 0.92/0.53 mm (winter), respectively. Compared to the unadjusted RMSE, the RMSE of GFPW and GMPW models in the corresponding seasons decreased by 0.72/0.86 mm (50.34%/60.13%), 1.11/1.26 mm (53.11%/60.28%), 0.75/0.76 mm (53.57%/54.28%) and 0.44/0.46 mm (47.82%/50.00%), respectively. The optimization performance of the GMPW model is commensurate in four seasons, and the accuracy optimization of the GFPW model is more pronounced in summer and winter. The RMSE improvement of the GMPW model in the four seasons is superior to that of the GFPW model by 0.14 mm (9.79%), 0.15 mm (7.17%), 0.01 mm (0.71%) and 0.02 mm (2.17%), correspondingly. This indicates that the accuracy of the GMPW model is better in spring and summer, and the accuracy of the two models is equivalent in fall and winter.
Compared with the GNSS PWV of the TP subregion, the unadjusted RMSE and the RMSE after adjustment by the GLPW model for the four seasons are 1.73/1.04 mm, 2.42/1.34 mm, 1.82/1.16 mm, and 1.33/1.08 mm, respectively. The RMSE adjusted by GFPW and GMPW models reduced by 0.8/0.87 mm (46.24%/50.28%), 1.50/1.58 mm (61.98%/65.28%), 0.67/0.94 mm (36.81%/51.64%) and 0.53/0.82 mm (39.84%/61.65%) compared to the unadjusted RMSE, respectively in the four seasons. This implies that GFPW and GMPW models show significant optimization performance in spring and summer. The RMSE improvement values between the GMPW and GFPW models for each season differ by 0.07 mm (4.04%), 0.08 mm (3.30%), 0.27 mm (14.83%) and 0.29 mm (21.80%), respectively. Therefore, the accuracy of the GMPW model in the TP subregion is equivalent in four seasons.
Compared with the GNSS PWV of the NC subregion, the unadjusted RMSE and the RMSE after adjustment by the GLPW model for the four seasons are 1.85/1.62 mm, 2.73/2.15 mm, 1.81/1.77 mm, and 1.17/1.13 mm, respectively. The RMSE adjusted by GFPW and GMPW models reduced by 0.56/0.68 mm (30.27%/36.75%), 0.85/0.99 mm (31.13%/36.26%), 0.45/0.99 mm (24.86%/54.69%) and 0.64/0.69 mm (54.70%/58.97%), respectively in each season. This implies that GFPW and GMPW models show significant optimization performance in four seasons, with improvements of 24–54% and 36–58%, respectively. The RMSE improvement values between the GFPW and GMPW models for each season differ by 0.12 mm (6.48%), 0.14 mm (5.12%), 0.54 mm (29.83%) and 0.05 mm (4.27%), respectively. This indicates that the accuracy of the GMPW model in the NC subregion is better in fall, while it is equivalent in spring, summer, and winter.
Compared with the GNSS PWV of the SC subregion, the unadjusted RMSE and the RMSE after adjustment by the GLPW model for the four seasons are 2.36/2.14 mm, 2.86/2.70 mm, 2.59/2.03 mm, and 1.99/1.91 mm, respectively. The RMSE adjusted by GFPW and GMPW models reduced by 1.02/1.44 mm (43,22%/58.47%), 1.44/1.83 mm (50.34%/63.98%), 1.11/1.52 mm (42.85%/58.65%) and 0.80/1.23 mm (40.20%/61.80%) compared to the unadjusted RMSE, respectively. This implies that GFPW and GMPW models show significant optimization performance in summer and fall. The RMSE improvement values between the GFPW and GMPW models for each season differ by 0.36 mm (15.25%), 0.39 mm (13.63%), 0.41 mm (15.83%), and 0.43 mm (21.60%). Therefore, the accuracy of the GMPW model in the SC subregion is equivalent in four seasons.
The results indicate that the accuracy of ERA5 PWV has been significantly improved after model adjustment, with the systematic differences between ERA5 PWV and GNSS PWV almost disappearing. Moreover, the GMPW model outperforms the GFPW model overall, demonstrating a significant improvement in RMSE accuracy during the summer and a slight enhancement in the fall and winter. The larger PWV errors in summer might result from higher water vapor content, while the reduced errors in fall and winter could stem from less variability in meteorological parameters and water vapor. Thus, PWV adjustment models that account for seasonal differences exhibit better performance and contribute to the production of high-quality PWV products in China.
4.2. Spatiotemporal Properties Analysis
To evaluate the temporal and spatial characteristics of the optimization performance of the PWV adjustment model, we applied seasonal PWV subregion adjustment models to adjust the corresponding ERA5 PWV.
Figure 7 and
Figure 8 show the distributions of Bias and RMSE between ERA5 PWV and GNSS PWV before and after optimization. The Bias and RMSE between unadjusted ERA5 PWV and GNSS PWV are notably high, especially in the southeastern coastal subregions. The adjusted Bias is nearly zero, and the overall RMSE has significantly decreased, suggesting that the PWV adjustment model effectively reduces systematic errors.
As shown in
Figure 7, the unadjusted ERA5 PWV in the NC subregion shows a pronounced positive Bias in spring and winter while exhibiting a significant negative Bias in the other subregions. This indicates that the accuracy of unadjusted ERA5 PWV at different subregional sites varies and exhibits notable regional distribution characteristics. Moreover, the Bias in the adjusted ERA5 PWV has significantly decreased across all seasons, indicating that the spatial variability and land–sea variability of PWV have improved from high latitudes to low latitudes and from coastal to inland areas. Therefore, it can be concluded that the GFPW and GMPW models can effectively adjust the Bias between ERA5 PWV and GNSS PWV.
As shown in
Figure 8, the RMSE after adjustment by the GLPW model is basically consistent with that before adjustment. This shows that the linear model is not very suitable for ERA5PWV adjustment in China. After adjustment by GFPW and GMPW, the seasonal RMSE of ERA5 PWV in each subregion is significantly reduced to a range of 0–2 mm, with an especially notable improvement in the SC subregion. It can be seen from
Figure 7 that the adjustment effect of the GMPW model on the NC subregion in summer is lower than the adjustment accuracy of the GFPW model. This may be because the summer climate changes are complex, and the GMPW model does not apply the ERA5 PWV adjustment for this site. The variations in PWV across different subregions are associated with geographical conditions, typically exhibiting a significant decrease from the southeastern coastal areas to the northwestern inland subregions, with the highest values observed in the southeastern coastal areas.
Furthermore, the GMPW model effectively adjusted the PWV differences between coastal and inland areas in the SC subregion, resulting in a more consistent RMSE distribution of ERA5 PWV across these areas compared to the unadjusted ERA5 PWV. The RMSE of ERA5 PWV adjusted by the GFPW model in the southwestern coastal subregion shows a significant difference compared to the inland areas, especially in three seasons other than winter.
The Bias and RMSE of unadjusted ERA5 PWV across the entire area range from −1.42 to 0.31 mm and 0.92 to 2.86 mm, respectively. After adjustment by the GLPW model, the RMSE distribution of ERA5 PWV in each subregion is generally consistent, ranging from 0.83 to 2.15 mm. For the GFPW model, the RMSE of adjusted ERA5 PWV in the NC, TP, and NWC subregions is concentrated within the range of 0.41 mm to 1.74 mm. After adjustment by the GMPW model, the RMSE distribution of ERA5 PWV in each subregion is generally consistent, ranging from 0.14 to 1.17 mm. In the SC subregion, following adjustment by the GFPW model, the RMSE values differ significantly between the southwest coastal area and the inland sites, with differences ranging from 1 to 2 mm. Compared to the unadjusted ERA5 PWV in the SC subregion, the RMSE differences for GFPW and GMPW models in spring, summer, fall, and winter are 1.34/0.78, 1.42/0.83, 1.48/0.85, and 1.19/0.06 mm, respectively. This may be related to the obvious climate change in the SC subregion and the more obvious water vapor fluctuations in the three seasons except winter [
71].
In conclusion, the GLPW model has the worst accuracy in adjusting ERA5 PWV. Both the GFPW and GMPW adjustment models effectively address the accuracy differences between ERA5 PWV and GNSS PWV, demonstrating robust accuracy and applicability. The GMPW model notably enhances the accuracy disparities between land and sea in ERA5 PWV and exhibits superior applicability in the SC subregion.
To evaluate the PWV adjustment model more comprehensively, we randomly selected 70 sites that did not participate in the modeling and applied seasonal PWV subregion adjustment models to adjust the corresponding ERA5 PWV. Taking the GNSS PWV as the reference, the Bias, MAE, and RMSE between the adjusted ERA5 PWV and the unadjusted ERA5 PWV of the GLPW, GFPW, and GMPW models are summarized in
Table 7.
As shown in
Table 7, the Bias adjusted by the models GFPW and GMPW are basically close to 0 mm, which shows that the systematic differences between ERA5 PWV and GNSS PWV are basically eliminated. The Bias adjusted by model GLPW is slightly lower than before unadjusted, but there is still a certain Bias, which shows that the GLPW model is less effective in adjusting ERA5 PWV. For MAE, the ranges of the GLPW model in the NWC, TP, NC, and SC subregions are 0.71–2.54 mm, 1.29–2.39 mm, 1.02–3.36 mm, and 1.76–3.95 mm, respectively. For MAE, the ranges of the GFPW model in the NWC, TP, NC, and SC subregions are 0.68–1.82 mm, 1.24–2.23 mm, 0.76–2.05 mm, and 1.47–2.60 mm, respectively, while the value range of the GMPW model is 0.65–1.75 mm, 0.96–1.89 mm, 0.74–2.00 mm, 1.43–2.28 mm. Compared with unadjusted PWV, the adjusted MAE of GLPW, GFPW, and GMPW models are reduced, and the GMPW model performs better.
As shown in
Figure 9 and
Figure 10, the Bias and RMSE spatial distribution between ERA5 PWV and GNSS PWV before and after adjustment. In
Figure 8, the Bias after adjustment of the GLPW model is basically the same as before the adjustment, which shows that the adjustment accuracy of the linear model is poor and is not suitable for ERA5 PWV water vapor adjustment. In addition, the Bias of the ERA5 PWV corrected by the GFPW and GMPW models is significantly reduced in each subregion, and the GMPW model adjustment results are better. After adjustment by GFPW and GMPW, the seasonal Bias of each partition of ERA5 PWV was significantly reduced to approximately 0, with the NC subregion improving particularly significantly. This shows that both the GFPW and GMPW models can adjust the ERA5 PWV of external sites well.
As shown in
Figure 10, the RMSE of ERA5 PWV at the external site of the GLPW model is basically the same after adjustment as before adjustment, which indicates that the linear model is not suitable for PWV adjustment in China. After adjustment by the GFPW and GMPW models, the seasonal RMSE in each subregion of the ERA5 PWV was significantly reduced to the range of 1–2 mm, with the NC subregion improving particularly significantly in fall and winter. In spring and summer, the adjusted RMSE of the GLPW model is around 2.5 mm, and the adjusted RMSE of the GFPW and GMPW models is around 2.0 mm. In addition, in the summer of the GFPW model, the RMSE value of the QHTR site on the west side of the TP subregion is very large, about 5 mm. This may be because the site is surrounded by high mountains near its location, resulting in complex precipitation there [
72]. Therefore, the adjusted RMSE of ERA5 PWV is larger. In conclusion, the GFPW and GMPW models are more suitable for China ERA5 PWV adjustment and have higher accuracy than the GLPW model.