With the increasing problem of global warming caused by the massive use of fossil fuels, biomass energy as a renewable energy source has attracted widespread attention throughout the globe. In this paper, we analyzed the spatial and temporal variation in wind energy in the East Asia and Western Pacific areas using IGRA site data, ERA5, and NCEP/NCAR reanalysis data from 2000 to 2021, and multi-variate empirical orthogonal function (MV-EOF) decomposition with the Pettitt mutation test, and the seasonal autoregression integrated moving average (SARIMA) model was used to predict the trend of wind speed. The spatial and temporal variations in wind energy in East Asia and Western Pacific areas were analyzed, and it was found that the richer wind-energy resources were mainly concentrated in the “Three Norths” (North China, Northwest China, and Northeast China) and Mongolia, followed by the Western Pacific areas. In addition, the T’ai-hang Mountains and the Qinghai-Tibet Plateau in China block the wind resources in the eastern and southern regions of East Asia, resulting in a shortage of wind resources in this region. In addition, the summer wind speed is significantly lower than in the other three seasons. The first-mode contributions of the MV-EOF wind field and geopotential heights, respectively, are 29.47% and 37.75%. The results show that: (1) There are significant seasonal differences in wind-energy resources in the study area, with the lowest wind speed in summer and the highest wind speed in winter. (2) The wind energy in the study area has significant regional characteristics. For example, China’s Qinghai-Tibet Plateau, Inner Mongolia, Xinjiang region, and Mongolia are rich in wind-energy resources. (3) Wind-energy resources in the study area have gradually increased since 2010, mainly due to changes in large-scale oceanic and atmospheric circulation patterns caused by global warming.

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.