Guojie Hu

Modeling permafrost properties in the Qinghai-Xizang (Tibet) Plateau

Water and heat dynamics in the active layer at a monitoring site in the Tanggula Mountains, located in the permafrost region of the Qinghai-Xizang (Tibet) Plateau (QXP), were studied using the physical-process-based COUPMODEL model, including the interaction between soil temperature and moisture under freeze-thaw cycles. Meteorological, ground temperature and moisture data from different depths within the active layer were used to calibrate and validate the model. The results indicate that the calibrated model satisfactorily simulates the soil temperatures from the top to the bottom of the soil layers as well as the moisture content of the active layer in permafrost regions. The simulated soil heat flux at depths of 0 to 20 cm was consistent with the monitoring data, and the simulations of the radiation balance components were reasonable. Energy consumed for phase change was estimated from the simulated ice content during the freeze/thaw processes from 2007 to 2008. Using this model, the active layer thickness and the energy consumed for phase change were predicted for future climate warming scenarios. The model predicts an increase of the active layer thickness from the current 330 cm to approximately 350–390 cm as a result of a 1–2°C warming. However, the effect active layer thickness of more precipitation is limited when the precipitation is increased by 20%–50%. The COUPMODEL provides a useful tool for predicting and understanding the fate of permafrost in the QXP under a warming climate.

Numerical Modeling of the Active Layer Thickness and Permafrost Thermal State Across Qinghai‐Tibetan Plateau

The dynamics of permafrost (including the permafrost thermal state and active layer thicknesses (ALT)) across the Qinghai-Tibetan Plateau (QTP) have not been well understood on a large scale. Here, we simulate the ALT and permafrost thermal state using the Geophysical Institute Permafrost Lab Version 2 (GIPL2) model across the QTP. Based on the single-point simulations, the model is upscaled to the entire QTP. The upscaled model is validated with five investigated regions (IRs), including Wenquan (WQIR), Gaize (GZIR), Aerjin (AEJIR), Xikunlun (XKLIR) and Qinghai-Tibetan highway (G109IR). The results show that the modified GIPL2 model improves the accuracy of the permafrost thermal state simulations. Due to our simulated results on the QTP, the average ALT is of 2.30 m (2.21 - 2.40 m). The ALT decreases with an increase in the altitude and decreases from the southeast to the northwest. The ALT is thin in the central QTP, but it is thick in the high-elevation mountain areas and some areas surrounding glaciers and lakes. The largest ALT is found in the border areas between permafrost and seasonally frozen ground regions. The simulated results of the MAGT (the mean annual ground temperature) indicate that most of the permafrost is sub-stable, which is sensitive to climate warming. The simulated results would be of great significance on assessing the impacts of permafrost dynamics on local hydrology, ecology, and engineering construction.

Modeling Hydrothermal Transfer Processes in Permafrost Regions of Qinghai-Tibet Plateau in China

Hydrothermal processes are key components in permafrost dynamics; these processes are integral to global warming. In this study the coupled heat and mass transfer model for (CoupModel) the soil-plant-atmosphere-system is applied in high-altitude permafrost regions and to model hydrothermal transfer processes in freeze-thaw cycles. Measured meteorological forcing and soil and vegetation properties are used in the CoupModel for the period from January 1, 2009 to December 31, 2012 at the Tanggula observation site in the Qinghai-Tibet Plateau. A 24-h time step is used in the model simulation. The results show that the simulated soil temperature and water content, as well as the frozen depth compare well with the measured data. The coefficient of determination (R2) is 0.97 for the mean soil temperature and 0.73 for the mean soil water content, respectively. The simulated soil heat flux at a depth of 0–20 cm is also consistent with the monitored data. An analysis is performed on the simulated hydrothermal transfer processes from the deep soil layer to the upper one during the freezing and thawing period. At the beginning of the freezing period, the water in the deep soil layer moves upward to the freezing front and releases heat during the freezing process. When the soil layer is completely frozen, there are no vertical water exchanges between the soil layers, and the heat exchange process is controlled by the vertical soil temperature gradient. During the thawing period, the downward heat process becomes more active due to increased incoming shortwave radiation at the ground surface. The melt water is quickly dissolved in the soil, and the soil water movement only changes in the shallow soil layer. Subsequently, the model was used to provide an evaluation of the potential response of the active layer to different scenarios of initial water content and climate warming at the Tanggula site. The results reveal that the soil water content and the organic layer provide protection against active layer deepening in summer, so climate warming will cause the permafrost active layer to become deeper and permafrost degradation.

Soil moisture and texture primarily control the soil nutrient stoichiometry across the Tibetan grassland

Soil nutrient stoichiometry and its environmental controllers play vital roles in understanding soil-plant interaction and nutrient cycling under a changing environment, while they remain poorly understood in alpine grassland due to lack of systematic field investigations. We examined the patterns and controls of soil nutrients stoichiometry for the top 10cm soils across the Tibetan ecosystems. Soil nutrient stoichiometry varied substantially among vegetation types. Alpine swamp meadow had larger topsoil C:N, C:P, N:P, and C:K ratios compared to the alpine meadow, alpine steppe, and alpine desert. In addition, the presence or absence of permafrost did not significantly impact soil nutrient stoichiometry in Tibetan grassland. Moreover, clay and silt contents explained approximately 32.5% of the total variation in soil C:N ratio. Climate, topography, soil properties, and vegetation combined to explain 10.3-13.2% for the stoichiometry of soil C:P, N:P, and C:K. Furthermore, soil C and N were weakly related to P and K in alpine grassland. These results indicated that the nutrient limitation in alpine ecosystem might shifts from N-limited to P-limited or K-limited due to the increase of N deposition and decrease of soil P and K contents under the changing climate conditions and weathering stages. Finally, we suggested that soil moisture and mud content could be good predictors of topsoil nutrient stoichiometry in Tibetan grassland.

Permafrost and land cover as controlling factors for light fraction organic matter on the southern Qinghai-Tibetan plateau

Permafrost degradation can stimulate the decomposition of organic soil matter and cause a large amount of greenhouse gas emissions into the atmosphere. The light fraction organic matter (LFOM) is a labile substrate for microbial decomposition and probably plays an important role in future permafrost carbon cycles. However, little is known about the distribution of LFOM and its relationship with permafrost and environmental factors. Here, we investigated the light fraction carbon (LFC) and nitrogen (LFN) contents and stocks under meadows and wet meadows with different permafrost conditions on the southern Qinghai-Tibetan Plateau. Our results showed that LFC and LFN were mainly distributed in the upper 30cm of soils, and the sites with permafrost had significantly higher contents of LFC and LFN than those from the sites without existing permafrost. The LFC and LFN decreased sharply with depth, suggesting that the soil organic matter (SOM) in this area was highly decomposed in deep soils. Soil moisture and bulk density explained approximately 50% of the variances in LFC and LFN for all the sampling sites, while soil moisture explained approximately 30% of the variance in permafrost sites. Both the C:N ratios and LFC:LFN ratios in the sites with permafrost were higher than those in the sites without permafrost. The results suggested that the permafrost and land cover types are the main factors controlling LFOM content and stock, and that permafrost degradation would lead to a decrease of LFOM and soil C:N ratios, thus accelerating the decomposition of SOM.

An analytical model for estimating soil temperature profiles on the Qinghai-Tibet Plateau of China

Soil temperature is a key variable in the control of underground hydro-thermal processes. To estimate soil temperature more accurately, this study proposed a solution method of the heat conduction equation of soil temperature (improved heat conduction model) by applying boundary conditions that incorporate the annual and diurnal variations of soil surface temperature and the temporal variation of daily temperature amplitude, as well as the temperature difference between two soil layers in the Tanggula observation site of the Qinghai-Tibet Plateau of China. We employed both the improved heat conduction model and the classical heat conduction model to fit soil temperature by using the 5 cm soil layer as the upper boundary for soil depth. The results indicated that the daily soil temperature amplitude can be better described by the sinusoidal function in the improved model, which then yielded more accurate soil temperature simulating effect at the depth of 5 cm. The simulated soil temperature values generated by the improved model and classical heat conduction model were then compared to the observed soil temperature values at different soil depths. Statistical analyses of the root mean square error (RMSE), the normalized standard error (NSEE) and the bias demonstrated that the improved model showed higher accuracy, and the average values of RMSE, bias and NSEE at the soil depth of 10–105 cm were 1.41°C, 1.15°C and 22.40%, respectively. These results indicated that the improved heat conduction model can better estimate soil temperature profiles compared to the traditional model.

Soil organic carbon and total nitrogen pools in permafrost zones of the Qinghai-Tibetan Plateau

There are several publications related to the soil organic carbon (SOC) on the Qinghai-Tibetan Plateau (QTP). However, most of these reports were from different parts of the plateau with various sampling depth. Here, we present the results from a systematic sampling and analysis of 200 soil pits. Most of the pits were deeper than 2 m from an east-west transect across the plateau. The SOC and total nitrogen (TN) pools of the 148 × 104 km2, the area of the permafrost zone, for the upper 2 m soils calculated from the vegetation map were estimated to be 17.07 Pg (interquartile range: 11.34–25.33 Pg) and 1.72 Pg (interquartile range: 1.08–2.06 Pg), respectively. We also predicted the distribution of land cover types in 2050 and 2070 using decision tree rules and climate scenarios, and then predicted SOC and TN pools of this region. The results suggested that the SOC and TN pools will decrease in the future. The results not only contribute to the carbon and nitrogen storage and stocks in the permafrost regions as a whole but most importantly, to our knowledge of the possible changes of C and N storage on the QTP in the future.

Soil infiltration processes of different underlying surfaces in the permafrost region on the Tibetan Plateau

ABSTRACT Soil infiltration processes were evaluated under field conditions by double-ring infiltrometers with different underlying surfaces in permafrost regions of the Tibetan Plateau. The results show that initial infiltration rates, stable soil infiltration rates and cumulative soil infiltration are strongly dependent on the underlying surface types, with the highest initial and stable soil infiltration rates in the alpine desert steppe, and the lowest in alpine meadow. The effects of soil moisture and texture on infiltration processes were also assessed. Within the same underlying surfaces, the values of infiltration parameters increased with the amount of vegetation cover, while soil moisture and soil infiltration rates displayed opposing trends, with fitting slopes of −0.03 and −0.01 for the initial and stable soil infiltration rates, respectively. The accuracies of the five models in simulating soil infiltration rates and seven models in predicting cumulative infiltration rates were evaluated against data generated from field experiments at four sites. Based on a comparative analysis, the Horton model provided the most complete understanding of the underlying surface effects on soil infiltration processes. Altogether, these findings show that different underlying surfaces can alter soil infiltration processes. This study provides a useful reference for understanding the parameterization of land surface processes for simulating changes in hydrological processes under global warming conditions in the permafrost region on the Tibetan Plateau.

Assessment of reanalysis soil moisture products in the permafrost regions of the central of the Qinghai–Tibet Plateau

The long-term and large-scale soil moisture (SM) record is important for understanding land atmosphere interactions and their impacts on the weather, climate, and regional ecosystem. SM products are one of the parameters used in some earth system models, but these records require evaluation before use. The water resources on the Qinghai-Tibet Plateau (QTP) are important to the water security of billions of people in Asia. Therefore, it is necessary to know the SM conditions on the QTP. In this study, the evaluation metrics of multilayer (0 - 10, 10 - 40 and 40 - 100 cm) SM in different reanalysis datasets of the European Centre for Medium-Range Weather Forecasts interim reanalysis (ERA), National Centers for Environmental Prediction Climate Forecast System and the Climate Forecast System version 2 (CFSv2), and China Meteorological Administration Land Data Assimilation System (CLDAS) are compared with in situ observations at five observation sites, which represent alpine meadow, alpine swamp meadow, alpine grassy meadow, alpine desert steppe, and alpine steppe environments during the thawing season from 1 Jan 2011 to 31 Dec 2013 on the QTP. The ERA SM remains constant at approximately 0.2 m3⋅m-3 at all observation sites during the entire thawing season. The CLDAS and CFSv2 SM products show similar patterns with those of the in situ SM observations during the thawing season. The CLDAS SM product performs better than the CFSv2 and ERA for all vegetation types except the alpine swamp meadow. The results indicate that the soil texture and land cover types play a more important role than the precipitation to increase the biases of the CLDAS SM product on the QTP.