Increasing greenhouse gas levels drive extensive changes in Arctic and cold-dominated environments, leading to a warmer, more humid, and variable climate. Associated permafrost thaw creates new groundwater flow paths in cold regions that are causing unprecedented environmental changes. This review of recent advances in groundwater research in cold environments has revealed that a new paradigm is emerging where groundwater is at the center of these changes. Groundwater flow and associated heat and solute transport are now used as a basis to understand hydrological changes, permafrost dynamics, water quality, integrity of infrastructure along with ecological impacts. Although major advances have been achieved in cold regions' cryohydrogeological research, the remaining knowledge gaps are numerous. For example, groundwater as a drinking water source is poorly documented despite its social importance. Lateral transport processes for carbon and contaminants are still inadequately understood. Numerical models are improving, but the highly complex physical-ecological changes occurring in the arctic involve coupled thermal, hydrological, hydrogeological, mechanical, and geochemical processes that are difficult to represent and hamper quantitative analysis and limit predictive capacity. Systematic long-term observatories where measurements involving groundwater are considered central are needed to help resolve these research gaps. Innovative transdisciplinary research will be critical to comprehend and predict these complex transformations.

The warming trend presents a significant threat to the underlying permafrost. Talik formation is widely recognized as a significant mechanism of permafrost degradation. Our research indicates that the term talik has undergone a long period of development and gradually formed, referring to unfrozen layers in permafrost. The talik has already resulted in extensive damage to the infrastructure built in permafrost areas. Here, we provide a brief overview of the current research status of talik. Accurately identifying talik presents a significant challenge. However, by integrating multiple identification tools with technology, the precision of talik detection can be enhanced, resulting in more accurate results. This paper discusses the strengths and weaknesses of each approach. While numerical simulations can enhance our understanding of the development mechanism and evolution process of taliks, most simulations focus on the evolution of taliks beneath lakes. These simulations emphasize the impact of subpermafrost groundwater flow on the development of lake taliks and the surrounding permafrost thickness. Today, there is a scarcity of relevant studies about taliks in cold zone engineering. The presence of talik exacerbates the occurrence of permafrost-related subgrade diseases, which are chronic and irreversible. Additionally, it poses a threat to the stability of the subgrades and worsens settlement issues. Therefore, we have analyzed the causes and distribution characteristics of talik beneath the subgrade and proposed a novel measure for preventing and controlling it. This measure aims to enhance the long-term service performance of subgrade in permafrost regions. The modified polyurethane material is injected into the talik through grouting technology as a replacement. This material has low thermal conductivity, strong water resistance, and certain strength. It effectively improves the hydrothermal environment conditions necessary for talik formation, preventing the formation of new taliks or impeding their development. As a result, the subgrade performance is enhanced.

The presence of taliks (perennially unfrozen zones in permafrost areas) adversely affects the thermal stability of infrastructure in cold regions, including roads. The role of heat advection on talik development and feedback on permafrost degradation has not been quantified methodically in this context. We incorporate a surface energy balance model into a coupled groundwater flow and energy transport numerical model (SUTRA-ice). The model, calibrated with long-term observations (1997-2018 on the Alaska Highway), is used to investigate and quantify the role of heat advection on talik initiation and development under a road embankment. Over the 25-year simulation period, the new model is driven by reconstructed meteorological data and has a good agreement with near surface soil temperatures. The model successfully reproduces the increasing depth to the permafrost table (mean absolute error <0.2 m), and talik development. The results demonstrate that heat advection provides an additional energy source that expedites the rate of permafrost thaw and roughly doubles the rate of permafrost table deepening, compared to purely conductive thawing. Talik initially formed and grew over time under the combined effect of water flow, snow insulation, road construction and climate warming. Talik formation creates a new thermal state under the road embankment, resulting in acceleration of underlying permafrost degradation, due to the positive feedback of heat accumulation created by trapped unfrozen water. In a changing climate, mobile water flow will play a more important role in permafrost thaw and talik development under road embankments, and is likely to significantly increase maintenance costs and reduce the long-term stability of the infrastructure.

Global warming has aggravated the problem of permafrost degradation, the long-term variation of which is not well estimated. In this study, a fully coupled hydrothermal dynamics cold region hydrologic model was used to quantify the historical spatial-temporal characteristics of permafrost degradation by estimating the depth to permafrost table (DPT), taliks in different subbasins, maximum frozen depth (MFD) for seasonally frozen ground (SFG), thaw date (TD), and annual number of frozen days (NFD) in the headwaters of the Yellow River (HWYR) in eastern High Mountain Asia. The model considered the taliks and successfully estimated four historical finescale permafrost maps. Permafrost degradation began with the enlargement of active layer thickness (ALT), formation of taliks, rise in the permafrost bottom, and merging of the permafrost table and bottom, and ended with a decrease in MFD. From 1980 to 2014, the permafrost area (PA) in the study area decreased by a total of 17,849 km2 at a rate of more than -7100 km2/decade (-5.8 %/decade). The mean DPT increased by 0.8 m at 31.76 cm/decade and more than 100 cm/decade at the grid scale. The absolute MFD (AMFD) decreased by 0.11 m at -14.47 cm/decade. The mean TD decreased by 14.28 days at -5.56 days/decade and more than -20 days/ decade at the grid scale. The mean NFD decreased by 25.51 days at -9.97 days/decade. The lowland area had a smaller NFD and an earlier TD than the mountainous area. Basins with a large air temperature (AT) and small elevation had a small DPT, MFD, TD, NFD, a large rate of decrease in TD and NFD, and increase in absolute DPT (ADPT), and a small rate of decrease in AMFD. The rates of increase in ADPT and decrease in AMFD were close to 0 when the change in AT was lower than 0.55 degrees C/decade.

Rising temperatures in the Arctic and subarctic are driving the rapid thaw of permafrost by reducing permafrost cooling, increasing active layer thickness, and promoting talik formation. In this study, the cyrohydrogeology of a permafrost mound located within the discontinuous permafrost zone near Umiujaq (Nunavik, Quebec, Canada) is characterized through the analysis of a dataset covering more than two decades of monitoring. This dataset captures a high degree of interannual variability in air temperature and ground thermal conditions, as well as the formation and closure of a supra-permafrost talik. Data indicate that variable saturation and advective heat transport directly contribute to the expansion and contraction of the talik. Data further indicate the presence of two distinct thermo-hydrologic settings resulting from differences in surface conditions, as well as subsurface thermal and flow regimes. The first, found at the top of the mound feature, is characterized by very low moisture contents (& lt;0.05 m(3)/m(3)), while the second, found at the side of the mound feature, shows higher annual moisture contents that strongly influence the dynamics of heat and groundwater flow. The data were synthesized into a detailed conceptual model of the cyrohydrogeological dynamics that highlights the important role of hydrogeological characterization and long-term data sets in understanding the effects of groundwater flow on seasonal frost and permafrost dynamics. Specifically, the results presented here show that in the absence of long-term data sets, longer-period transient phenomena such as talik opening and closure may be misrepresented as uni-directional feedback loops, as opposed to highly dynamic temporary phenomena.

Lakes represent as much as similar to 25% of the total land surface area in lowland permafrost regions. Though decreasing lake area has become a widespread phenomenon in permafrost regions, our ability to forecast future patterns of lake drainage spanning gradients of space and time remain limited. Here, we modeled the drivers of gradual (steady declining lake area) and catastrophic (temporally abrupt decrease in lake area) lake drainage using 45 years of Landsat observations (i.e. 1975-2019) across 32 690 lakes spanning climate and environmental gradients across northern Alaska. We mapped lake area using supervised support vector machine classifiers and object based image analyses using five-year Landsat image composites spanning 388 968 km(2). Drivers of lake drainage were determined with boosted regression tree models, using both static (e.g. lake morphology, proximity to drainage gradient) and dynamic predictor variables (e.g. temperature, precipitation, wildfire). Over the past 45 years, gradual drainage decreased lake area between 10% and 16%, but rates varied over time as the 1990s recorded the highest rates of gradual lake area losses associated with warm periods. Interestingly, the number of catastrophically drained lakes progressively decreased at a rate of similar to 37% decade(-1) from 1975-1979 (102-273 lakes draining year(-1)) to 2010-2014 (3-8 lakes draining year(-1)). However this 40 year negative trend was reversed during the most recent time-period (2015-2019), with observations of catastrophic drainage among the highest on record (i.e. 100-250 lakes draining year(-1)), the majority of which occurred in northwestern Alaska. Gradual drainage processes were driven by lake morphology, summer air and lake temperature, snow cover, active layer depth, and the thermokarst lake settlement index (R (2) (adj) = 0.42, CV = 0.35, p < 0.0001), whereas, catastrophic drainage was driven by the thawing season length, total precipitation, permafrost thickness, and lake temperature (R (2) (adj) = 0.75, CV = 0.67, p < 0.0001). Models forecast a continued decline in lake area across northern Alaska by 15%-21% by 2050. However these estimates are conservative, as the anticipated amplitude of future climate change were well-beyond historical variability and thus insufficient to forecast abrupt 'catastrophic' drainage processes. Results highlight the urgency to understand the potential ecological responses and feedbacks linked with ongoing Arctic landscape reorganization.

Permafrost thaw leads to thermokarst lake formation and talik growth tens of meters deep, enabling microbial decomposition of formerly frozen organic matter (OM). We analyzed two 17-m-long thermokarst lake sediment cores taken in Central Yakutia, Russia. One core was from an Alas lake in a Holocene thermokarst basin that underwent multiple lake generations, and the second core from a young Yedoma upland lake (formed similar to 70 years ago) whose sediments have thawed for the first time since deposition. This comparison provides a glance into OM fate in thawing Yedoma deposits. We analyzed total organic carbon (TOC) and dissolved organic carbon (DOC) content, n-alkane concentrations, and bacterial and archaeal membrane markers. Furthermore, we conducted 1-year-long incubations (4 degrees C, dark) and measured anaerobic carbon dioxide (CO2) and methane (CH4) production. The sediments from both cores contained little TOC (0.7 +/- 0.4 wt%), but DOC values were relatively high, with the highest values in the frozen Yedoma lake sediments (1620 mg L-1). Cumulative greenhouse gas (GHG) production after 1 year was highest in the Yedoma lake sediments (226 +/- 212 mu g CO2-C g(-1) dw, 28 +/- 36 mu g CH4-C g(-1) dw) and 3 and 1.5 times lower in the Alas lake sediments, respectively (75 +/- 76 mu g CO2-C g(-1) dw, 19 +/- 29 mu g CH4-C g(-1) dw). The highest CO2 production in the frozen Yedoma lake sediments likely results from decomposition of readily bioavailable OM, while highest CH4 production in the non-frozen top sediments of this core suggests that methanogenic communities established upon thaw. The lower GHG production in the non-frozen Alas lake sediments resulted from advanced OM decomposition during Holocene talik development. Furthermore, we found that drivers of CO2 and CH4 production differ following thaw. Our results suggest that GHG production from TOC-poor mineral deposits, which are widespread throughout the Arctic, can be substantial. Therefore, our novel data are relevant for vast ice-rich permafrost deposits vulnerable to thermokarst formation.

The impact of permafrost thaw on hydrologic, thermal, and biotic processes remains uncertain, in part due to limitations in subsurface measurement capabilities. To better understand subsurface processes in thermokarst environments, we collocated geophysical and biogeochemical instruments along a thaw gradient between forested permafrost and collapse-scar bogs at the Alaska Peatland Experiment site near Fairbanks, Alaska. Ambient seismic noise monitoring provided continuous high-temporal resolution measurements of water and ice saturation changes. Maps of seismic velocity change identified areas of large summertime velocity reductions nearest the youngest bog, indicating potential thaw and expansion at the bog margin. These results corresponded well with complementary borehole nuclear magnetic resonance measurements of unfrozen water content with depth, which showed permafrost soils nearest the bog edges contained the largest amount of unfrozen water along the study transect, up to 25% by volume. In situ measurements of methane within permafrost soils revealed high concentrations at these bog-edge locations, up to 30% soil gas. Supra-permafrost talik zones were observed at the bog margins, indicating talik formation and perennial liquid water may drive lateral bog expansion and enhanced permafrost carbon losses preceding thaw. Comparison of seismic monitoring with wintertime surface carbon dioxide fluxes revealed differential responses depending on time and proximity to the bogs, capturing the controlling influence of subsurface water and ice on microbial activity and surficial emissions. This study demonstrates a multidisciplinary approach for gaining new understanding of how subsurface physical properties influence greenhouse gas production, emissions, and thermokarst development.

Climate warming and anthropogenic impact causes transformation of geocryological conditions in the river basins of the North-East of Russia. Changes in the thickness of the active layer, configuration of taliks, types of landscapes and other factors lead to transformation of water exchange processes between surface and groundwater runoff. This is manifested in the seasonal redistribution of the components of the water balance, accelerated melting of aufeis, change in the ratio of waters of different genesis in the structure of river runoff. As a result, natural and anthropogenic risks that affect the safe and efficient development of infrastructure and socio-economic processes are increasing. At the same time the system of observations developed in the Soviet period has been practically destroyed in the region. This paper offers a vision of organizing complex multidisciplinary research to assess and project the changes in the conditions of underground and surface water interaction in natural and disturbed river basins of the cryolithozone of the North-East of Russia, including for solving applied problems, based on permafrost, hydrology. hydrogeology, landscape science and geophysics with applications of remote sensing and field research integrated through mathematical modeling methods. To achieve the goal, the identification of natural and disturbed landscapes using remote sensing data. and key areas for detailed research will be selected. Geophysical and drilling works will be carried out within the sites to establish permafrost-hydrogeological conditions, monitoring stations will be equipped to determine hydrogeological, hydrometeorological and geocryological characteristics, including sampling for isotopic and hydrogeochemical studies. As the main key sites, it is proposed to use the area of the Kolyma water-balance station and the site on Anmangynda aufeis, for which there are long-term observation series in the 20th century. Field data will become the basis for improving the mathematical model of runoff formation, considering the relationship between groundwater and river runoff in the conditions of permafrost. Mathematical modeling will make it possible to quantitatively analyze the water balance of rivers considering various factors and project water availability both for specific industrial facilities and for the region as a whole.

As the Arctic warms and wildfire occurrence increases, talik formation in permafrost regions is projected to expand and affect the cycling of water and carbon. Yet, few unified field and modeling studies have examined this process in detail, particularly in areas of continuous permafrost. We address this gap by presenting multimethod, multiseasonal geophysical measurements of permafrost and liquid-water content that reveal substantial talik development in response to recent wildfire in continuous permafrost of boreal Alaska. Results from observation-based cryohydrogeologic model simulations suggest that predisturbance subsurface conditions are key factors influencing thaw response to fire disturbance and air temperature warming. Our high-resolution integrated study illustrates enhanced vulnerability of boreal continuous permafrost, with observed talik formation that exceeds coarse-scale model projections by similar to 100 years even under the most extreme future emissions scenario. Results raise important scaling questions for representing extreme permafrost thaw phenomena of growing widespread importance in large-scale predictive models.