In order to alleviate global warming and the energy crisis, it is of great significance to develop and utilize solar energy resources. Sunshine duration (SD) is considered to be the best substitute for solar radiation and a key factor in evaluating solar energy resources. Therefore, the spatial and temporal characteristics of SD and the reasons for its changes have received extensive attention and discussion. Based on the data of 415 meteorological stations from 1970 to 2019, this paper uses linear trend analysis, Mann–Kendall mutation analysis, the Hurst index, empirical orthogonal decomposition, correlation analysis and partial correlation analysis to analyze the spatiotemporal characteristics of SD and its relationship with influencing factors. The results show that the annual SD in China shows a downward trend, with a climate trend rate of −37.93 h/10a, and a significant decline from 1982 to 2019. The seasonal SD shows a downward trend, and the downward trend is most obvious in summer. The annual and seasonal SD will still show a downward trend in the future. The spatial distribution of SD not only has an overall consistent distribution but also takes the Yellow River from Ningxia to Shandong as the boundary, showing a north–south opposite distribution. Annual SD has a significant positive correlation, a significant negative correlation, a positive correlation and a negative correlation with wind speed, precipitation, temperature and relative humidity, respectively, and it is most closely related to wind speed and precipitation. In addition, the change in SD may also be related to human activities.

Through the observation of ground-based LIDAR and satellite sensors, the weather conditions of continuous foggy-hazy alternations in the Huainan region from 26 December 2016 to 5 January 2017 were analyzed and observed. In this study, the formation and influence of this event were discussed by analyzing pollutant concentrations, meteorological factors and aerosol optical characteristics. The concentrations of PM₁₀ and PM_2.5 increased significantly. The maximum value of PM₁₀ was 412 μg/m³, and the maximum value of PM_2.5 was 258 μg/m³. The transportation of pollutants and the production of man-made pollutants promote the accumulation of pollutants. In this weather process, meteorological factors such as the surface wind speed, humidity, surface temperature, and inversion also promote the accumulation of pollutants, which is the main reason for the formation of this weather process. Furthermore, the near surface air mass mainly came from the cities near the Huainan region and the heavily polluted areas in the north, while the upper air mass came from Inner Mongolia. In this paper, piecewise inversion was adopted to achieve accurate all-weather extinction coefficient profile inversion by reasonably selecting a cloud LIDAR ratio through a backscatter ratio, and the LIDAR ratio of cloud in this period was 22.57–34.14 Sr. By means of extinction coefficient inversion and correlation analysis, the correlation index of PM_2.5 and the aerosol optical depth (AOD) was 0.7368, indicating that there was a positive correlation between PM_2.5 and AOD, and AOD can also reflect the pollution condition of this region. The formation process of foggy-hazy weather in the Huainan region studied in this paper can provide a research basis for foggy-hazy pollution in this region.

Given the increasing impact of extreme rainfall and flooding on human life, studying and predicting changes in atmospheric water vapor (AWV) becomes particularly important. This paper analyzes the moderate-resolution imaging spectroradiometer (MODIS) data of the East Asian region from January 2003 to February 2023. The AWV data are examined in the time and frequency domain using methods such as empirical orthogonal function (EOF), Mann–Kendall (MK) analysis, and others. Additionally, four prediction models are applied to forecast the monthly average AWV data for the next year. The accuracy of these models is evaluated using metrics such as mean square error (MSE), mean absolute error (MAE), root mean square error (RMSE), and coefficient of determination (R²). The findings reveal several key insights: (1) The East Asian region exhibits highly variable seasonal variability in AWV, with identified mutation points after the MK test. (2) Spatial analysis shows high AWV data in the southern coastal areas of China, Thailand, Myanmar, Nansha Islands, and other regions during winter, while the Qinghai-Tibet Plateau region experiences low AWV during summer. (3) The first mode obtained through EOF decomposition contributes over 60% of the variance. Analysis of this mode reveals an increasing trend in AWV data for regions such as the Indian peninsula, Mongolia, and central and northeastern China over the past nine years. Conversely, the Bay of Bengal, Spratly Islands, eastern coast, and certain areas display a decreasing trend. (4) Employing the ensemble empirical mode decomposition (EEMD), the study identifies AWV data as a non-stationary series with an overall decreasing trend from 2003 to 2022. The filtered AWV series undergoes fast Fourier transform (FFT), uncovering periodicities of 2.6 years, 5 years, and 19 years. (5) Among the four forecasting models compared, the seasonal autoregressive integrated moving average model (SARIMA) demonstrates superior performance with the smallest MSE of 0.00782, MAE of 0.06977, RMSE of 0.08843, and the largest R² value of 0.98454. These results clearly indicate that the SARIMA model provides the best fit. Therefore, the SARIMA forecasting model can be effectively utilized for forecasting AWV data, offering valuable insights for studying weather variability.

Ozone is a very important trace gas in the atmosphere, it is like a “double-edged sword”. Because the ozone in the stratosphere can effectively help the earth’s organisms to avoid the sun’s ultraviolet radiation damage, the ozone near the ground causes pollution. Therefore, it is essential to explore the time-frequency domain variation characteristics of total column ozone and have a better understanding of its cyclic variation. In this paper, based on the monthly scale dataset of total column ozone (TCO) (September 2002 to February 2023) from Atmospheric Infrared Sounder (AIRS) carried by NASA’s Aqua satellite, linear regression, coefficient of variation, Mann-Kendall (M-K) mutation tests, wavelet analysis, and empirical orthogonal function decomposition (EOF) analysis were used to analyze the variation characteristics of the TCO in China from the perspectives of time domain, frequency domain, and spatial characteristics. Finally, this study predicted the future of TCO data based on the seasonal autoregressive integrated moving average (SARIMA) model in the time series algorithm. The results showed the following: (1) From 2003 to 2022, the TCO in China showed a slight downward trend, with an average annual change rate of −0.29 DU/a; the coefficient of variation analysis found that TCO had the smallest intra-year fluctuations in 2008 and the largest intra-year fluctuations in 2005. (2) Using the M-K mutation test, it was found that there was a mutation point in the total amount of column ozone in 2016. (3) Using wavelet analysis to analyze the frequency domain characteristics of the TCO, it was observed that TCO variation in China had a combination of 14-year, 6-year, and 4-year main cycles, where 14 years is the first main cycle with a 10-year cycle and 6 years is the second main cycle with a 4-year cycle. (4) The spatial distribution characteristics of the TCO in China were significantly different in each region, showing a distribution characteristic of being high in the northeast and low in the southwest. (5) Based on the EOF analysis of the TCO in China, it was found that the variance contribution rate of the first mode was as high as 52.85%, and its spatial distribution of eigenvectors showed a “-” distribution. Combined with the trend analysis of the time coefficient, this showed that the TCO in China has declined in the past 20 years. (6) The SARIMA model with the best parameters of (1, 1, 2) × (0, 1, 2, 12) based on the training on the TCO data was used for prediction, and the final model error rate was calculated as 1.34% using the mean absolute percentage error (MAPE) index, indicating a good model fit.

To comprehensively explore the characteristics of global SST anomalies, a novel time–frequency combination method, based on the COBE data and NCEP/NCAR reanalysis products in the past 100 years, was developed. From the view of the time domain, the global SST generally showed an upward trend from 1920 to 2019, the upward trend was significant after 1988, and the growth mutation occurred in 1930, according to the Mann–Kendall (MK) mutation test. Moreover, we extracted spatiotemporal modes of SST anomalies’ variability by empirical orthogonal function (EOF) analysis and obtained global spatial EOFs that closely correspond to regionally defined climate modes. Our results demonstrated that El Niño–Southern Oscillation (ENSO) is the typical character for the first mode of SST anomaly EOF, and Atlantic multidecadal oscillation (AMO) for the second. From the view of the frequency domain, our data suggested that there is a multi-period nesting phenomenon in global SST variations, in which the first main cycle with the most obvious oscillation was a 30-year cycle and changed in 20-year cycles, and the second cycle was a 15-year cycle and changed in 10-year cycles through wavelet analysis. As for the perspective of time–frequency characteristics, the dominant period of ENSO in the first mode of EOF is 4 years, obtained through filtering and cross wavelet transform. In addition, SST anomalies will maintain an upward trend for the next 60 months, according to the seasonal autoregressive integrated moving average (SARIMA) model, which has the potential value for predicting ENSO.

Outgoing longwave radiation (OLR) data are one of the key factors in studying the radiation balance of the earth–atmosphere system in East Asia. It is of great significance to explore the influence factors on OLR. This paper processes the data of nearly 19 years, from September 2002 to February 2022, and conducts in-depth research using the exponential smoothing method, empirical orthogonal decomposition (EOF), correlation analysis, and other methods. We found that the spatial distribution of OLR is zonal symmetry and gradually decreases with the increase of latitude. Using EOF analysis, it is found that the total variance contribution of the first four decomposed spatial features exceeds 70%, and the overall change trend of the four-time coefficients in the past 19 years all show a downward trend. OLR is positively correlated with total column water vapor (TCWV), air temperature (AT), and cloud top temperature (CTT), but negatively correlated with cloud top pressure (CTP). OLR has a similar spatial correlation distribution with TCWV and AT, while the spatial correlation between OLR and CTP is opposite to the first two parameters. In most parts of East Asia, the spatial correlation with CTT exceeds 0.8. The change in OLR value is affected by various meteorological parameters. In East Asia, the positive correlation between 30° N and 60° N is significantly affected by TCWV, AT, and CTT; and the negative correlation is more significantly affected by CTP. At 0–25° N, the positive correlation is significantly affected by CTP and CTT, while the negative correlation is significantly affected by TCWV and AT.