Soil freeze–thaw processes are remarkable features of the land surface across the Tibetan Plateau (TP). Soil moisture and temperature fluctuate during the freeze–thaw cycle, affecting the soil water and energy exchange between the land and atmosphere. This study investigates variations in the soil temperature, humidity, and freeze–thaw state and their responses to air temperature and precipitation on the TP from 1981 to 2016. Regional simulations of the TP using Community Land Model version 4.5 demonstrate that the climate of the TP has become warmer and wetter over the past 37 years, with increases in both regional average temperature and precipitation. Using empirical orthogonal function analysis and the Mann–Kendall trend test of air temperature, we show that 1980–1998 was relatively cold, and 1999–2016 was relatively warm. Soil temperature and moisture in most areas of the TP were affected by air temperature and precipitation, and both showed an upward trend during the past 37 years. Overall, from 1981 to 2016, the freezing date of the TP has become delayed, the thawing date has been hastened, and the duration of the freeze–thaw state has shortened. The surface soil freezes and thaws first, and these processes pervade deeper soil with the passage of time; freeze–thaw processes have an obvious hysteresis. Precipitation and air temperature had marked effects on the freeze–thaw processes. Higher air temperatures delay the freezing date, hasten the thawing date, and shorten the freeze–thaw period. Areas with the highest precipitation saw late soil freeze, early thaw, and the shortest freeze–thaw duration. Areas with less vegetation froze earlier and thawed later. The freeze–thaw duration increased in the northwest of the plateau and decreased on the rest of the plateau. This article informs research on frozen soil change in the context of global warming.

Through the observation and study of the raindrop spectrum, we can not only explore the evolution law of precipitation but also understand the microphysical characteristics of different types of precipitation clouds. This paper uses the raindrop spectrum data observed at Naqu Station, Yushu Station, Linzhi Station, and Emei Mount Station, as well as cloud radar data in the Yushu region in the Tibetan Plateau. Raindrop spectral characteristics are studied, and the raindrop size distribution (DSD) characteristics of the four stations are analyzed. The results are as follows: (1) The overall raindrop spectral concentration of the four stations decreases with the increase in particle size after the peak value. The downtrend is most gentle at the peak. (2) All rain intensity levels show the characteristics of multiple vertices. At Linzhi Station, Naqu Station and Yushu Station, the particles with a diameter of <1 mm contributed the most to the precipitation rate, while the particles with a diameter of <1 mm at Emei Mount Station contributed the least to the precipitation rate. (3) The precipitation in the central and eastern margins of the Tibetan Plateau is dominated by small and medium-sized particles, accounting for 95–99% of the precipitation particle number concentration. (4) The raindrop spectra of the three types of precipitation clouds at the four sites are all cumulonimbus precipitation Cumulonimbus in the uppermost layer, stratiform cloud precipitation stratiform cloud in the lowermost layer, and mixed cloud precipitation stratiform cloud in between. (5) Cloud radar analysis of precipitation structure shows that cumulonimbus and mixed clouds develop vigorously, and the highest cloud height can reach 13 km, but the average precipitation duration of mixed clouds is shorter than that of cumulonimbus; stratiform clouds develop relatively smoothly, with cloud heights of 6–7 km, the average precipitation duration is the shortest.

Based on 35-year (1981–2015) ensemble (24 members) hindcasts of the IAP/LASG global seasonal prediction system named FGOALS-f2 V1.0 (FGOALS-f2), the tropical cyclone (TC) seasonal prediction skills over the South China Sea (SCS) during the TC peak season (July–November) are evaluated. Starting the prediction from June 20th, FGOALS-f2 can well capture the seasonal mean characteristics for both the genesis location and track of TCs over the SCS. For seasonal anomalous TC numbers, FGOALS-f2 underestimates the maximum and minimum of the TC number compared to the observation. The temporal correlation coefficients (TCCs) between FGOALS-f2 and the observation are 0.39 for the TC number and 0.51 for accumulated cyclone energy (ACE) over the SCS, respectively, which are both above the 95% significant level. Additionally, FGOALS-f2 has acceptable prediction skill for the seasonal mean number of TCs landing on three areas (coastal southeastern China, Indochina Peninsula, and Philippines) surrounding the SCS. The skillful prediction of SCS TCs could be ascribed to the well-predicted tropical anomaly of sea surface temperature (SSTA), TC and El Niño-Southern Oscillation (TC-ENSO) relations, and Genesis potential index (GPI).

Based on empirical orthogonal function (EOF) analysis, the dominant modes of variations in summer surface sensible heating (SH) over the Tibetan Plateau (TP), as well as the associated atmospheric circulation anomalies, were investigated in this study. The results show that the first dominant mode of summer SH presents a feature of decadal reduction over the whole TP on an interdecadal time scale, and the second dominant mode is characterized by a zonally asymmetric pattern with positive (negative) SH anomalies in the western (eastern) TP on an interannual time scale. The variations of summer SH are dominated by anomalies in downwelling surface shortwave radiation (DSWR), which are associated with atmospheric circulation changes. The first dominant mode of variation in SH is connected to the interdecadal variation of the Silk Road Pattern (SRP). Further analysis reveals that the interdecadal phase shift of the SRP induces anticyclone circulation to the northeast of the TP, leading to enhanced water vapor supply and convergence over the TP. This can lead to an increase in the total cloud cover, and a reduction in DSWR, contributing to the decadal reduction in SH over the TP. The second dominant mode of variation in SH is related to a stationary teleconnection pattern over the Eurasian continent named the North Atlantic-East and North Asia pattern (NAENA). Corresponding to the positive phase of the NAENA, there is a cyclone anomaly to the west TP, leading to anomalous water vapor convergence (divergence) over the eastern (western) TP. This can result in enhanced (decreased) cloud cover, reduced (increased) DSWR, and therefore, an anomalous decrease (enhancement) in SH over the east (west) of the TP. Furthermore, the southwesterly wind anomaly, which is accompanied by the anomalous cyclone to the west TP, leads to positive SH in the western TP.

Changes in the surface fluxes cause changes in the annular flow field over a region, and they affect the transport of water vapor. To study the influence of the changes in the surface flux on the water vapor transport in the upper layer in the canyon area of southeastern Tibet, in this study, the water vapor transport characteristics were analyzed using the HYSPLIT_v4 backward trajectory model at Danka and Motuo stations in the canyons in the southeastern Tibetan Plateau from November 2018 to October 2019. Then, using ERA-5 reanalysis data from 1989 to 2019 and the characteristics of the high-altitude water vapor transportation, the impact of the surface flux changes on the water vapor transportation was analyzed using singular value decomposition (SVD). The results show that the main sources of the water vapor in the study area were from the west and southwest during the non-Asian monsoon (non-AMS), while there was mainly southwest air flow and a small amount of southeast air flow in the lower layer during the Asian monsoon (AMS) at the stations in southeastern Tibet. The water vapor transmission channel of the westward airflow is higher than 3000 m, and the water vapor transmission channel of the southwestward and southeastward airflow is about 2000 m. The sensible heat and latent heat are negatively correlated with water vapor flux divergence. The southwest boundary of southeastern Tibet is a key area affecting water vapor flux divergence. When the sensible heat and latent heat exhibit downward trends during the non-Asian monsoon season, the eastward water vapor flux exhibits an upward trend. During the Asian monsoon season, when the sensible heat and latent heat in southeastern Tibet increase as a whole, the eastward water vapor flux in the total-column of southeastern Tibet increases.

Mount Emei is located on the eastern edge of the Tibetan Plateau, on the transition zone between the main body of the Tibetan Plateau and the Sichuan Basin in China. It is not only the necessary place for the eastward movement of the plateau system but also the place where the southwest vortex begins to develop. Its special geographical location makes it particularly important to understand the turbulence characteristics and surface energy balance of this place. Based on the Atmospheric Boundary Layer (ABL) tower data, radiation observation data and surface flux data of Mount Emei station on the eastern edge of the Tibetan Plateau from December 2019 to February 2022, the components of surface equilibrium are estimated by the eddy correlation method and Thermal Diffusion Equation and Correction (TDEC) method, the characteristics of surface energy exchange in the Mount Emei area are analyzed, and the aerodynamic and thermodynamic parameters are estimated. The results show that the annual average value of zero-plane displacement $\mathrm{d}$ is 10.45 m, the annual average values of aerodynamic roughness $Z_{0m}$ and aerothermal roughness $Z_{0h}$ are 1.61 and 1.67 m, respectively, and the annual average values of momentum flux transport coefficient $C_D$ and sensible heat flux transport coefficient $C_H$ are $1.58\times {10}^{-2}$ and $3.79\times {10}^{-3}$, respectively. The dimensionless vertical wind fluctuation variance in the Mount Emei area under unstable conditions can better conform to the 1/3rd power law of the Monin–Obukhov similarity theory, while the dimensionless horizontal wind fluctuation variance under unstable lamination and the dimensionless 3D wind fluctuation variance under stable condition does not conform to this law. In the near-neutral case, the dimensionless velocity variance in the vertical direction in this area is 1.314. The daytime dominance of sensible and latent heat fluxes varied seasonally, with latent heat fluxes dominating in summer and sensible heat transport dominating in winter. he surface albedo of Mount Emei in four seasons is between 0.04 and 0.08. The surface albedo in summer and autumn is higher than that in Mount Emei. The influence of the underlying surface on surface reflectance is much greater than other factors, such as altitude, longitude and latitude. The non-closure phenomenon is significant in the Mount Emei area. The energy closure rates before and after considering canopy thermal storage are 46% and 48%, respectively. The possible reason for the energy non-closure in this area is that the influence of horizontal advection and vertical advection on the energy closure is not considered.