In permafrost regions, vegetation growth is influenced by both climate conditions and the effects of permafrost degradation. Climate factors affect multiple aspects of the environment, while permafrost degradation has a significant impact on soil moisture and nutrient availability, both of which are crucial for ecosystem health and vegetation growth. However, the quantitative analysis of climate and permafrost remains largely unknown, hindering our ability to predict future vegetation changes in permafrost regions. Here, we used statistical methods to analyze the NDVI change in the permafrost region from 1982 to 2022. We employed correlation analysis, multiple regression residual analysis and partial least squares structural equation modeling (PLS-SEM) methods to examine the impacts of different environmental factors on NDVI changes. The results show that the average NDVI in the study area from 1982 to 2022 is 0.39, with NDVI values in 80% of the area remaining stable or exhibiting an increasing trend. NDVI had the highest correlation with air temperature, averaging 0.32, with active layer thickness coming in second at 0.25. Climate change plays a dominant role in NDVI variations, with a relative contribution rate of 89.6%. The changes in NDVI are positively influenced by air temperature, with correlation coefficients of 0.92. Although the active layer thickness accounted for only 7% of the NDVI changes, its influence demonstrated an increasing trend from 1982 to 2022. Overall, our results suggest that temperature is the primary factor influencing NDVI variations in this region.

Accurately quantifying the impact of permafrost degradation and soil freeze-thaw cycles on hydrological processes while minimizing the reliance on observational data are challenging issues in hydrological modeling in cold regions. In this study, we developed a modular distributed hydro-thermal coupled hydrological model for cold regions (DHTC) that features a flexible structure. The DHTC model couples heat-water transport processes by employing the conduction-advection heat transport equation and Richard equation considering ice-water phase change. Additionally, the DHTC model integrates the influence of organic matter into the hydrothermal parameterization scheme and includes a subpermafrost module based on the flow duration curve analysis to estimate cold-season streamflow sustained by subpermafrost groundwater. Moreover, we incorporated energy consumption due to ice phase changes to the available energy, enhancing the accuracy of evaporation estimation in cold regions. A comprehensive evaluation of the DHTC model was conducted. At the point scale, the DHTC model accurately replicates daily soil temperature and moisture dynamics at various depths, achieving average R-2 of 0.98 and 0.87, and average RMSE of 0.61degree celsius and 0.03 m(3)m(-3), respectively. At the basin scale, DHTC outperformed (Daily: R-2 = 0.66, RMSE = 0.75 mm; Monthly: R-2 = 0.90, RMSE = 15.7 mm) the GLDAS/FLDAS Noah, GLDAS/VIC, and PML-V2 models in evapotranspiration simulation. The DHTC model also demonstrated reasonable performance in simulating daily (NSE = 0.70, KGE = 0.84), monthly (NSE = 0.86, KGE = 0.90), and multi-year monthly (NSE = 0.97, KGE = 0.93) streamflow in the Source Regions of Yangtze River. DHTC also successfully reproduced the snow depth in basin-averaged time series and spatial distributions (RMSE = 0.86 cm). The DHTC model provides a robust tool for exploring the interactions between permafrost and hydrological processes, and their responses to climate change.

Soil freeze-thaw cycles play a critical role in ecosystem, hydrological and biogeochemical processes, and climate. The Tibetan Plateau (TP) has the largest area of frozen soil that undergoes freeze-thaw cycles in the low-mid latitudes. Evidence suggests ongoing changes in seasonal freeze-thaw cycles during the past several decades on the TP. However, the status of diurnal freeze-thaw cycles (DFTC) of shallow soil and their response to climate change largely remain unknown. In this study, using in-situ observations, the latest reanalysis, machine learning, and physics-based modeling, we conducted a comprehensive assessment of the spatiotemporal variations of DFTC and their response to climate change in the upper Brahmaputra (UB) basin. About 24 +/- 8% of the basin is subjected to DFTC with a mean frequency of 87 +/- 55 days during 1980-2018. The area and frequency of DFTC show small long-term changes during 1980-2018. Air temperature impacts on the frequency of DFTC changes center mainly around the freezing point (0 degrees C). The spatial variations in the response of DFTC to air temperature can primarily be explained by three factors: precipitation (30.4%), snow depth (22.6%) and seasonal warming/cooling rates (14.9%). Both rainfall and snow events reduce diurnal fluctuations of soil temperature, subsequently reducing DFTC frequency, primarily by decreasing daytime temperature through evaporation-cooling and albedo-cooling effects, respectively. These results provide an in-depth understanding of diurnal soil freeze-thaw status and its response to climate change. Freeze-thaw transitions of terrestrial landscapes are a common phenomenon in cold regions. The seasonal and diurnal freeze-thaw cycles (DFTC) of shallow soil exhibit substantial differences in response to climate. Understanding of the spatiotemporal patterns of DFTC and their response to climate change remains limited over the Tibetan Plateau (TP), which is characterized by the largest areas of freeze-thaw terrain in the mid- and low-latitudes of the world. We found the frequency and area of DFTC show a slight increase trend in a significantly warming climate in upper Brahmaputra (UB) basin, the largest river basin of the TP. The variation of DFTC depends on climatic conditions, with soils near the freezing point (0 degrees C) being more susceptible to changes in DFTC. Precipitation, snow depth and seasonal warming/cooling rates are the top three factors influencing the response of DFTC to air temperature changes. Snowfall plays a more important role in the temporal variability of DFTC frequency than rainfall. The number of diurnal freeze-thaw cycles (DFTC) in shallow soil increase slightly during the period 1980-2018 in the upper Brahmaputra (UB) basin Air temperature effects on the changes in DFTC frequency center on the freezing point Snowfall plays a more important role in the temporal variability of DFTC than rainfall

Rapid surface and subsurface changes in the Arctic polygonal tundra landscapes due to the melting of ice wedges, known as thermokarst processes, have significant implications for Arctic ecosystems. However, the integration of thermokarst processes into widely used global climate models for projections poses an important question. Here we use an integrated permafrost thermal hydrology model to explore the decoupled nature of two thermokarst processes - microtopography evolution and ground subsidence - in six Arctic locations. Our study specifically investigates this decoupled nature during the transformation of poorly drained low-centered polygons to welldrained high-centered polygons. Spanning diverse climates in polygonal tundra landscapes under the RCP8.5 climate scenario, our findings reveal small variations in permafrost thaw and ground subsidence rates - 2-10 % and 2-4 %, respectively - with and without the representation of microtopography evolution. This suggests that neglecting surface microtopography and its evolution is unlikely to have significant impacts on permafrost projections, regardless of the climate and location. As a result, we suggest the representation of microtopography in Earth System Models may not be imperative. Disclaimer: Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the author(s) and do not necessarily reflect the view of the U.S. Department of Commerce, National Oceanic and Atmospheric Administration.

Arctic land is characterized by a high surface and subsurface heterogeneity on different scales. However, the effects of land surface model resolution on fluxes and soil state variables in the Arctic have never been systematically studied, even though smaller scale heterogeneities are resolved in high-resolution land boundary condition datasets. Here, we compare 210 km and 5 km setups of the land surface model JSBACH3 for an idealized case study in eastern Siberia to investigate the effects of high versus low-resolution land boundary conditions on simulating the interactions of soil physics, hydrology and vegetation. We show for the first time that there are differences in the spatial averages of the simulated fluxes and soil state variables between resolution setups. Most differences are small in the summer mean, but larger within individual months. Heterogeneous soil properties induce large parts of the differences while vegetation characteristics play a minor role. Active layer depth shows a statistically significant increase of +20% in the 5 km setup relative to the 210 km setup for the summer mean and +43% for August. The differences are due to the nonlinear vertical discretization of the soil column amplifying the impact of the heterogeneous distributions of soil organic matter content and supercooled water. Resolution-induced differences in evaporation fluxes amount to +43% in July and are statistically significant. Our results show that spatial resolution significantly affects model outcomes due to nonlinear processes in heterogenous land surfaces. This suggests that resolution needs to be accounted in simulations of land surface models in the Arctic.

This study uses a new dataset on gauge locations and catchments to assess the impact of 21st-century climate change on the hydrology of 221 high-mountain catchments in Central Asia. A steady-state stochastic soil moisture water balance model was employed to project changes in runoff and evaporation for 2011-2040, 2041-2070, and 2071-2100, compared to the baseline period of 1979-2011. Baseline climate data were sourced from CHELSA V21 climatology, providing daily temperature and precipitation for each subcatchment. Future projections used bias-corrected outputs from four General Circulation Models under four pathways/scenarios (SSP1 RCP 2.6, SSP2 RCP 4.5, SSP3 RCP 7.0, SSP5 RCP 8.5). Global datasets informed soil parameter distribution, and glacier ablation data were integrated to refine discharge modeling and validated against long-term catchment discharge data. The atmospheric models predict an increase in median precipitation between 5.5% to 10.1% and a rise in median temperatures by 1.9 degrees C to 5.6 degrees C by the end of the 21st century, depending on the scenario and relative to the baseline. Hydrological model projections for this period indicate increases in actual evaporation between 7.3% to 17.4% and changes in discharge between + 1.1% to -2.7% for the SSP1 RCP 2.6 and SSP5 RCP 8.5 scenarios, respectively. Under the most extreme climate scenario (SSP5-8.5), discharge increases of 3.8% and 5.0% are anticipated during the first and second future periods, followed by a decrease of -2.7% in the third period. Significant glacier wastage is expected in lower-lying runoff zones, with overall discharge reductions in parts of the Tien Shan, including the Naryn catchment. Conversely, high-elevation areas in the Gissar-Alay and Pamir mountains are projected to experience discharge increases, driven by enhanced glacier ablation and delayed peak water, among other things. Shifts in precipitation patterns suggest more extreme but less frequent events, potentially altering the hydroclimate risk landscape in the region. Our findings highlight varied hydrological responses to climate change throughout high-mountain Central Asia. These insights inform strategies for effective and sustainable water management at the national and transboundary levels and help guide local stakeholders.

Global warming has shown an Arctic amplification effect in recent decades, leading to pronounced changes in pan-Arctic soil surface temperature (SST). SST plays a direct role in energy exchange between soil and atmosphere and serves as an indicator of the land-atmosphere energy balance. Remote sensing land surface temperature (LST) data is able to indicate near-surface temperature, but influences from environment factors, such as vegetation and snow, can introduce biases between LST and SST. In this study, the importances of five environment factors (vegetation, snow, surface soil composition, topography, and solar radiation) to monthly mean SST estimation from MODIS LST in pan-Arctic were analyzed. Then a method for pan-Arctic monthly mean SST estimation from MODIS LST by incorporating these environment factors and monthly-based modeling based on random forest (RF) algorithm was proposed. The results reveal that all the selected environment factors contribute to monthly-based modeling, with vegetation exerting the greatest importance from May to October and snow in March and April. The root mean square error (RMSE) of pan-Arctic monthly SST estimated by the proposed method from 2003 to 2022 ranges from 0.89 to 1.88 degrees C, which is a 42.95---53.35 % reduction compared to the widely used season-based multivariate linear regression (MLR) models based solely on LST (RMSE between 1.56 and 4.03 degrees C). The accuracy is notably improved in areas with lower and no vegetation (grassy woodlands, grasslands, permanent wetlands, and barrens) in the cold season (September to the following April), and in higher vegetation (forests) areas in the warm season (May to August). The proposed method can contribute to producing high-precision monthly mean SST data from LST, estimating permafrost extent and active layer thickness, and understanding the land-atmosphere energy balance in pan-Arctic.

Modeling Arctic-Boreal vegetation is a challenging but important task, since this highly dynamic ecosystem is undergoing rapid and substantial environmental change. In this work, we synthesized information on 18 dynamic vegetation models (DVMs) that can be used to project vegetation structure, composition, and function in North American Arctic-Boreal ecosystems. We reviewed the ecosystem properties and scaling assumptions these models make, reviewed their applications from the scholarly literature, and conducted a survey of expert opinion to determine which processes are important but lacking in DVMs. We then grouped the models into four categories (specific intention models, forest species models, cohort models, and carbon tracking models) using cluster analysis to highlight similarities among the models. Our application review identified 48 papers that addressed vegetation dynamics either directly (22) or indirectly (26). The expert survey results indicated a large desire for increased representation of active layer depth and permafrost in future model development. Ultimately, this paper serves as a summary of DVM development and application in Arctic-Boreal environments and can be used as a guide for potential model users, thereby prioritizing options for model development.

Fluoroquinolones, a class of antibiotics, have been detected in various aquatic environments, including those experiencing freeze-thaw cycles. This study investigated the adsorption of ciprofloxacin (CIP) in frozen (-21 degrees C) and aqueous (25 degrees C) solutions under varying pH levels, electrolyte types, and ionic strengths. CIP sorption on goethite was found to be transient, as freezing re-establishes equilibrium, nearly doubling CIP loadings at acidic to circumneutral pH values. The original equilibrium was restored by thawing. Our investigation reveals that ion pairs, formed between the positively charged piperazine group of CIP and anions (Cl-, Br-, and NO3-), create a charge-shielding effect, explaining the transient nature of CIP sorption equilibrium at goethite-water interfaces. In situ ATR-FTIR observations and model predictions further confirm the significant role of ion-paired surface complexes in transient CIP sorption. The transience of CIP sorption equilibrium in frozen and aqueous solutions is attributed to the local concentrations of anions, which undergo freeze-concentration into liquid intergrain boundaries and dilution by reversible ice nucleation and thawing. As the interaction between the hydrosphere and cryosphere intensifies with climate change, these findings have significant implications for evaluating the fate of contaminants in both terrestrial and aquatic environments.

Recent research on the Himalayan cryosphere has increasingly been focused on climate uncertainty and regional variations, considering features such as glacier recession, lake expansion, outburst floods, and regional hazards. The Bhilangana river basin, located in the central Himalayas, is predominantly characterized by increased elevation-dependent warming and declining seasonal precipitation. Our study shows that high-elevation temperature increased from 2000 to 2022 (0.05(degrees)C/year, p = 20 m/sec). Quantification of the regional hazard reveals potentially severe downstream challenges for low-to-medium-scale hydropower stations, local settlements, and road and railway bridges near Devling and Ghuttu villages.