In this study, we used hourly observations to investigate the cooling effect of summer rainfall on surface air temperature (Ta) in a subtropical area, Guangdong province, South China. Data were categorized step-by-step by rainfall system (convection, monsoon, and typhoon), daily rainfall amount, and relative humidity (RH) level. Moreover, the average hourly Ta variation due to solar radiation was removed from all observations before statistical analysis. The results showed that the linear relationship between hourly Ta variation and rainfall intensity did not exist. However, the cooling effect of rainfall on Ta variation was dominant. In addition, convective rainfall does cause a greater temperature drop than the other two rainfall systems. After further partitioning all samples by RH level preceding the rainfall, the relationship between hourly Ta variation and rainfall intensity became distinctive. When RH was below 70%, rainfall-induced cooling became more substantial and scaled linearly with event intensity, but when RH exceeded 70%, the rainfall cooling effect was generally restrained by the RH increase. A strong correlation between hourly Ta variation and RH level preceding the rainfall suggests the importance of RH on the rainfall cooling effect.

In the context of global warming, compound dry-hot events (CDHEs) are intensifying in Guangdong, yet CMIP6 projections remain uncertain. This study employs CMIP6 data and the Standardized Compound Event Indicator (SCEI) to quantify CDHEs severity, applying an observational constraint approach to reduce inter-model uncertainty. The results show that, after observational constraint, uncertainties decrease by about 63% and 77% in Period I and II under SSP126 and by about 57% and 59% under SSP585, greatly improving projection robustness. CDHE risk is highest in SSP585-Period II. Future dry-hot intensification in Guangdong generally increases from north to south, with western Guangdong most strongly affected. Although CDHEs weaken in other periods, western Guangdong shows persistent aggravation. Mechanism analyses indicate that SSP585-Period I is mainly linked to cold sea surface temperature (SST) anomalies in the South Atlantic and waters near Australia. After correction, dry-hot conditions show a marked weakening across Guangdong, although slight intensification persists over the Leizhou Peninsula. SSP585-Period II is primarily influenced by warm SST anomalies in the eastern Pacific and South Atlantic and cold anomalies in the North Pacific. The two SSP126 periods are associated with warm SST anomalies in the South Atlantic and waters near Australia and with cold anomalies in the South Atlantic, North Pacific, and North Atlantic, respectively. After correction, CDHEs generally weaken across Guangdong, although southern and south-central areas remain relatively severe. These findings indicate that historical key SST biases can strongly influence future CDHEs projections in Guangdong by modulating large-scale atmospheric circulation, including the Pacific-South American wave train, Indian Ocean SST anomalies, and the Western North Pacific Subtropical Anticyclone.

Shelterbelts are increasingly used to mitigate strong wind damage, but the complex canopy structures create challenges for numerical studies of windbreak effectiveness, such as the trade-off between computational cost and accuracy of results. To address these challenges and accurately investigate the downstream wind fields, most conventional studies represent shelterbelts as rectangular porous media with a uniformly distributed aerodynamic resistance coefficient. However, due to the vertical variation in canopy diameter and the irregular distribution of leaf density, the aerodynamic resistance of natural shelterbelts becomes nonuniform accordingly. To quantify the discrepancies arising from this simplification, this study first proposes a non-destructive approach to calculate canopy porosity profiles, which are further used to derive aerodynamic resistance at different heights. Then, by comparing the results obtained from the conventional and proposed approaches in Large-Eddy Simulations, the discrepancies caused by ignoring the vertical variation in canopy structures are analyzed. Finally, these discrepancies are further investigated for double-row shelterbelts. The results show that ignoring the vertical variation in canopy diameter leads to significant differences in windbreak effectiveness, especially for the downstream velocity and pressure fields at the top and middle heights of the canopy. The proposed approach provides a computationally efficient and more accurate representation of near-surface wind fields downstream of shelterbelts, thereby contributing to the accurate prediction of local wind fields for meteorological services.

Accurate prediction of electricity demand during hot seasons is essential for maintaining power system reliability, particularly in humid subtropical regions such as Guangdong, China, where high temperatures strongly influence consumption. However, many models rely primarily on air temperature and may not fully capture combined atmospheric effects. This study evaluates the potential of biometeorological thermal indices for improving summer electricity load forecasting. Daily maximum load and meteorological data during May–September 2019–2021 were analyzed using Back-Propagation Neural Network (BP), Random Forest (RF), and a Stacking ensemble model. Three indices—Effective Temperature (ET), Physiological Equivalent Temperature (PET), and the Universal Thermal Climate Index (UTCI)—were introduced as predictors. The ensemble model achieved the best performance, with Ensemble–UTCI yielding the highest accuracy (R² = 0.559, RMSE = 60.96 × 10⁴ kW, MAE = 45.10 × 10⁴ kW). Compared with temperature-based models, biometeorological indices consistently improved predictions, with UTCI performing best (average RMSE = 62.81 × 10⁴ kW). Bayesian analysis shows strong evidence of improvement in RF and ensemble models, but not in BP or linear models, indicating model dependence. During the July 2021 heat event, RF showed greater robustness, with PET–RF achieving the lowest error (MAPE = 3.03%). These results demonstrate the value of biometeorological indices for load forecasting in humid subtropical regions.