Yonghua Zhao

Surface heat sources play a critical role in shaping meteorological conditions and permafrost dynamics on the Tibetan Plateau. To better understand surface heat source variability in the permafrost region, multiyear observational data from an isolated permafrost site and a continuous permafrost site were used to analyse the variability. The results indicated that the surface heat source at the isolated permafrost site remained relatively stable, whereas that at the continuous permafrost site increased significantly, at a rate of approximately 2.2 Wm⁻²yr⁻¹. The proportions of sensible and latent heat fluxes in the surface heat source budget varied with the season. Overall, the proportion of latent heat flux was greater than that of sensible heat flux in summer and autumn, whereas the proportion of sensible heat flux was greater than that of latent heat flux in winter and spring. The peak proportion for both fluxes exceeded 80%. A random forest model effectively captured the variations in the surface heat source. And the machine learning simulation indicated that soil temperature, downward shortwave radiation and albedo were identified as major contributors to the surface heat source, collectively accounting for more than 88% of the overall variability at both sites. The impacts of snow cover on the surface heat source varied with intensity. The increasing trend of the surface heat source was closely related to climate warming in autumn and winter. A significant but weak positive correlation was observed between vegetation and the surface heat source.

The Complementary Relationship (CR) principle of evapotranspiration provides an efficient approach for estimating actual evapotranspiration (ET_a), owing to its simplified computation and effectiveness in utilizing meteorological factors. Accurate estimation of actual evapotranspiration (ET_a) is crucial for understanding surface energy and water cycles, especially in permafrost regions. This study aims to evaluate the applicability of two Complementary Relationship (CR)-based methods—Bouchet’s in 1963 and Brutsaert’s in 2015—for estimating ETa on the Qinghai–Tibetan Plateau (QTP), using observations from Eddy Covariance (EC) systems. The potential evapotranspiration (ET_p) was calculated using the Penman equation with two wind functions: the Rome wind function and the Monin–Obukhov Similarity Theory (MOST). The comparison revealed that Bouchet’s method underestimated ET_a during frozen soil periods and overestimated it during thawed periods. In contrast, Brutsaert’s method combined with the MOST yielded the lowest RMSE values (0.67–0.70 mm/day) and the highest correlation coefficients (r > 0.85), indicating superior performance. Sensitivity analysis showed that net radiation (Rn) had the strongest influence on ET_a, with a daily sensitivity coefficient of up to 1.35. This study highlights the improved accuracy and reliability of Brutsaert’s CR method in cold alpine environments, underscoring the importance of considering freeze–thaw dynamics in ET modeling. Future research should incorporate seasonal calibration of key parameters (e.g., ε) to further reduce uncertainty.

Accurate permafrost mapping in mountainous regions is hindered by sparse in situ observations and heterogeneous terrain. This study develops a GPR-augmented machine learning framework to map mountain permafrost in the northeastern Qinghai–Tibet Plateau. A total of 1037 presence–absence samples were compiled from boreholes, soil pits, 128 GPR transects collected in 2009, and 22 additional empirical points above 4700 m, covering diverse topographic and thermal conditions. Thirteen classification algorithms were evaluated using 5-fold cross-validation repeated 40 times, with LightGBM, CatBoost, XGBoost, and RF achieving top performance (F1 > 0.98). Elevation-based spatial comparisons revealed that LightGBM and CatBoost produced more terrain-adaptive predictions at high altitudes and slope transitions. Aspect-controlled permafrost boundaries were captured, with modeled lower elevation limits varying by >200 m across slope directions. SHAP analysis showed that climate and soil variables contributed nearly 80% to model outputs, with LST, FDD, BD, and TDD being dominant. Several predictors exhibited threshold or nonlinear responses, reinforcing their physical relevance. Additional experiments confirmed that integration of GPR and high-elevation constraint samples significantly improved model generalization, especially in underrepresented terrain zones. This study demonstrates that a GPR-augmented machine learning framework can support cost-effective, physically informed mapping of frozen ground in complex alpine environments.

The functional relationship between soil permittivity and soil water content serves as the theoretical foundation for electromagnetic wave-based techniques used to determine soil moisture levels. However, the response of permittivity to changes in soil water content varies significantly across different soil types. Current models that utilize soil permittivity to estimate soil water content are often based on empirical statistical relationships specific to particular soil types. Moreover, existing physical models are hindered by an excessive number of parameters, which can be difficult to measure or calculate. This study introduces a universal model, termed the Soil Refractive Index (SRI) model, to describe the relationship between soil permittivity and soil water content. The SRI model is derived from the propagation velocity of electromagnetic waves in various soil components and the functional relationship between electromagnetic wave velocity and relative permittivity. The SRI model expresses soil water content as a linear function of the square root of the relative permittivity for any soil type with the slope and intercept as the two undetermined parameters. The slope is primarily influenced by the relative permittivity of soil water, while the intercept is mainly affected by both the slope and the soil porosity. The applicability of the SRI model is validated through tested soil samples and comparison with previously published empirical statistical models. For dielectric lossless soil, the theoretical value of the slope is calculated to be 0.126. The intercept varies across different soil types and increases linearly with soil porosity. The SRI model provides a theoretical basis for calculating soil water content using permittivity across various soil types.

Permafrost distribution is of great significance for the study of climate, ecology, hydrology, and infrastructure construction in high-cold mountain regions with complex topography. Therefore, updated high-resolution permafrost distribution mapping is necessary and highly demanded in related fields. This case study conducted in a small catchment in the northeast of the Qinghai Tibet Plateau proposes a new method of using ground-penetrating radar (GPR) to detect the stratigraphic structure, interpret the characteristics of frozen ground, and extract the boundaries of permafrost patches in mountain areas. Thus, an empirical–statistical model of mountain frozen ground zonation, along with aspect (ASP) adjustment, is established based on the results of the GPR data interpretation. The spatial mapping of the frozen ground based on this model is compared with a field survey dataset and two existing permafrost distribution maps, and their consistencies are all higher than 80. In addition, the new map provides more details on the distribution of frozen ground. In this case, the influence of ASP on the distribution of permafrost in mountain areas is revealed: the adjustment of ASP on the lower limit of continuous and discontinuous permafrost is 180–200 m, the difference in the annual mean ground temperature between sunny and shady slopes is up to 1.4–1.6 °C, and the altitude-related temperature variation and uneven distribution of solar radiation in different ASPs comprehensively affect the zonation of mountain frozen ground. This work supplements the traditional theory of mountain permafrost zonation, the results of which are of value to relevant scientific studies and instructive to engineering construction in this region.

In order to improve the utilization efficiency of cultural buildings, architects often design them as multi-purpose venues. This paper introduces the design of a multi-purpose concert hall built in Fuzhou New District. In order to meet multipurpose performances, including large-scale symphony and organ music, small-scale orchestral and solo performance, adjustable acoustic curtains have been implemented into this hall for tuning its acoustical response. The results of simulation and measurement show that this hall has adjustable reverberation in suitable range and is respond to the set design demands.

Spatial information of particle size fractions (PSFs) is primary for understanding the thermal state of permafrost in the Qinghai-Tibet Plateau (QTP) in response to climate change. However, the limitation of field observations and the tremendous spatial heterogeneity hamper the digital mapping of PSF. This study integrated log-ratio transformation approaches, variable searching methods, and machine learning techniques to map the surficial soil PSF distribution of two typical permafrost regions. Results showed that the Boruta technique identified different covariates but retained those covariates of vegetation and land surface temperature in both regions. Variable selection techniques effectively decreased the data redundancy and improved model performance. In addition, the spatial distribution of soil PSFs generated by four log-ratio models presented similar patterns. Isometric log-ratio random forest (ILR-RF) outperformed the other models in both regions (i.e., R² ranged between 0.36 to 0.56, RMSE ranged between 0.02 and 0.10). Compared with three legacy datasets, our prediction better captured the spatial pattern of PSFs with higher accuracy. Although this study largely improved the accuracy of spatial distribution of soil PSFs, further endeavors should also be made to improve model accuracy and interpretability for a better understanding of the interaction and processes between environmental predictors and soil PSFs at permafrost regions.

Ground-penetrating radar (GPR) is a convenient geophysical technique for active-layer soil moisture detection in permafrost regions, which is theoretically based on the petrophysical relationship between soil moisture (θ) and the soil dielectric constant (ε). The θ–ε relationship varies with soil type and thus must be calibrated for a specific region or soil type. At present, there is lack of such a relationship for active-layer soil moisture estimation for the Qinghai–Tibet plateau permafrost regions. In this paper, we utilize the Complex Refractive Index Model to establish such a calibration equation that is suitable for active-layer soil moisture estimation with GPR velocity. Based on the relationship between liquid water, temperature, and salinity, the soil water dielectric constant was determined, which varied from 84 to 88, with an average value of 86 within the active layer for our research regions. Based on the calculated soil-water dielectric constant variation range, and the exponent value range within the Complex Refractive Index Model, the exponent value was determined as 0.26 with our field-investigated active-layer soil moisture and dielectric data set. By neglecting the influence of the soil matrix dielectric constant and soil porosity variations on soil moisture estimation at the regional scale, a simple active-layer soil moisture calibration curve, named CRIM, which is suitable for the Qinghai–Tibet plateau permafrost regions, was established. The main shortage of the CRIM calibration equation is that its calculated soil-moisture error will gradually increase with a decreasing GPR velocity and an increasing GPR velocity interpretation error. To avoid this shortage, a direct linear fitting calibration equation, named as υ-fitting, was acquired based on the statistical relationship between the active-layer soil moisture and GPR velocity with our field-investigated data set. When the GPR velocity interpretation error is within ±0.004 m/ns, the maximum moisture error calculated by CRIM is within 0.08 m³/m³. While when the GPR velocity interpretation error is larger than ±0.004 m/ns, a piecewise formula calculation method, combined with the υ-fitting equation when the GPR velocity is lower than 0.07 m/ns and the CRIM equation when the GPR velocity is larger than 0.07 m/ns, was recommended for the active-layer moisture estimation with GPR detection in the Qinghai–Tibet plateau permafrost regions.