Abstract. Climate change is profoundly transforming the carbon-rich Arctic tundra landscape, potentially moving it from a carbon sink to a carbon source by increasing the thickness of soil that thaws on a seasonal basis. However, the modeling capability and precise parameterizations of the physical characteristics needed to estimate projected active layer thickness (ALT) are limited in Earth system models (ESMs). In particular, discrepancies in spatial scale between field measurements and Earth system models challenge validation and parameterization of hydrothermal models. A recently developed surface-subsurface model for permafrost thermal hydrology, the Advanced Terrestrial Simulator (ATS), is used in combination with field measurements to achieve the goals of constructing a process-rich model based on plausible parameters and to identify fine-scale controls of ALT in ice-wedge polygon tundra in Barrow, Alaska. An iterative model refinement procedure that cycles between borehole temperature and snow cover measurements and simulations functions to evaluate and parameterize different model processes necessary to simulate freeze-thaw processes and ALT formation. After model refinement and calibration, reasonable matches between simulated and measured soil temperatures are obtained, with the largest errors occurring during early summer above ice wedges (e.g., troughs). The results suggest that properly constructed and calibrated onedimensional thermal hydrology models have the potential to provide reasonable representation of the subsurface thermal response and can be used to infer model input parameters and process representations. The models for soil thermal conductivity and snow distribution were found to be the most sensitive process representations. However, information on lateral flow and snowpack evolution might be needed to constrain model representations of surface hydrology and snow depth.

Through taliks-thawed zones extending through the entire permafrost layer-represent a critical type of heterogeneity that affects water redistribution and heat transport, especially in sloping landscapes. The formation of through taliks as part of the transition from continuous to discontinuous permafrost creates new hydrologic pathways connecting the active layer to sub-permafrost regions, with significant hydrological and biogeochemical consequences. At hilly field sites in the southern Seward Peninsula, AK, patches of deep snow in tall shrubs are associated with higher winter ground temperatures and an anomalously deep active layer. To better understand the thermal-hydrologic controls and consequences of through taliks, we used the coupled surface/subsurface permafrost hydrology model ATS (Advanced Terrestrial Simulator) to simulate through taliks associated with preferentially distributing snow. Scenarios were developed based on an intensively studied hillslope transect on the southern Seward Peninsula, which predominately has taller shrubs midslope and tundra in upslope and downslope areas. The model was forced with detrended meteorological data with snow preferentially distributed at the midslope of the domain to investigate the potential role of vegetation-induced snow trapping in controlling through talik development under conditions typical of the current-day Seward Peninsula. We simulated thermal hydrology and talik development for five permafrost conditions ranging in thickness from 17-45 m. For the three thinnest permafrost configurations, a through talik developed, which allowed water from the seasonally thawed layer into sub-permafrost waters, increasing sub-permafrost groundwater flow. These numerical experiments suggest that in the transition from continuous to discontinuous permafrost, through taliks may appear at locations that preferential trap snow and that the appearance of those through taliks may drive significant changes in permafrost hydrology.

Modeling and observation of ground temperature dynamics are the main tools for understanding current permafrost thermal regimes and projecting future thaw. Until recently, most studies on permafrost have focused on vertical ground heat fluxes. Groundwater can transport heat in both lateral and vertical directions but its influence on ground temperatures at local scales in permafrost environments is not well understood. In this study we combine field observations from a subarctic fen in the sporadic permafrost zone with numerical simulations of coupled water and thermal fluxes. At the Tavvavuoma study site in northern Sweden, ground temperature profiles and groundwater levels were observed in boreholes. These observations were used to set up one‐ and two‐dimensional simulations down to 2 m depth across a gradient of permafrost conditions within and surrounding the fen. Two‐dimensional scenarios representing the fen under various hydraulic gradients were developed to quantify the influence of groundwater flow on ground temperature. Our observations suggest that lateral groundwater flow significantly affects ground temperatures. This is corroborated by modeling results that show seasonal ground ice melts 1 month earlier when a lateral groundwater flux is present. Further, although the thermal regime may be dominated by vertically conducted heat fluxes during most of the year, isolated high groundwater flow rate events such as the spring freshet are potentially important for ground temperatures. As sporadic permafrost environments often contain substantial portions of unfrozen ground with active groundwater flow paths, knowledge of this heat transport mechanism is important for understanding permafrost dynamics in these environments.

Abstract. The effects of soil property uncertainties on permafrost thaw projections are studied using a three-phase subsurface thermal hydrology model and calibration-constrained uncertainty analysis. The null-space Monte Carlo method is used to identify soil hydrothermal parameter combinations that are consistent with borehole temperature measurements at the study site, the Barrow Environmental Observatory. Each parameter combination is then used in a forward projection of permafrost conditions for the 21st century (from calendar year 2006 to 2100) using atmospheric forcings from the Community Earth System Model (CESM) in the Representative Concentration Pathway (RCP) 8.5 greenhouse gas concentration trajectory. A 100-year projection allows for the evaluation of predictive uncertainty (due to soil property (parametric) uncertainty) and the inter-annual climate variability due to year to year differences in CESM climate forcings. After calibrating to measured borehole temperature data at this well-characterized site, soil property uncertainties are still significant and result in significant predictive uncertainties in projected active layer thickness and annual thaw depth-duration even with a specified future climate. Inter-annual climate variability in projected soil moisture content and Stefan number are small. A volume- and time-integrated Stefan number decreases significantly, indicating a shift in subsurface energy utilization in the future climate (latent heat of phase change becomes more important than heat conduction). Out of 10 soil parameters, ALT, annual thaw depth-duration, and Stefan number are highly dependent on mineral soil porosity, while annual mean liquid saturation of the active layer is highly dependent on the mineral soil residual saturation and moderately dependent on peat residual saturation. By comparing the ensemble statistics to the spread of projected permafrost metrics using different climate models, we quantify the relative magnitude of soil property uncertainty to another source of permafrost uncertainty, structural climate model uncertainty. We show that the effect of calibration-constrained uncertainty in soil properties, although significant, is less than that produced by structural climate model uncertainty for this location.

[1] This study presents a stochastic inverse method for aquifer structure identification using sparse geophysical and hydraulic response data. The method is based on updating structure parameters from a transition probability model to iteratively modify the aquifer structure and parameter zonation. The method is extended to the adaptive parameterization of facies hydraulic parameters by including these parameters as optimization variables. The stochastic nature of the statistical structure parameters leads to nonconvex objective functions. A multi-method genetically adaptive evolutionary approach (AMALGAM-SO) was selected to perform the inversion given its search capabilities. Results are obtained as a probabilistic assessment of facies distribution based on indicator cokriging simulation of the optimized structural parameters. The method is illustrated by estimating the structure and facies hydraulic parameters of a synthetic example with a transient hydraulic response.

Abstract. The goal of this research is to constrain the influence of ice wedge polygon microtopography on near-surface ground temperatures. Ice wedge polygon microtopography is prone to rapid deformation in a changing climate, and cracking in the ice wedge depends on thermal conditions at the top of the permafrost; therefore, feedbacks between microtopography and ground temperature can shed light on the potential for future ice wedge cracking in the Arctic. We first report on a year of sub-daily ground temperature observations at 5 depths and 9 locations throughout a cluster of lowcentered polygons near Prudhoe Bay, Alaska, and demonstrate that the rims become the coldest zone of the polygon during winter, due to thinner snowpack. We then calibrate a polygon-scale numerical model of coupled thermal and hydrologic processes against this dataset, achieving an RMSE of less than 1.1 • C between observed and simulated ground temperature. Finally, we conduct a sensitivity analysis of the model by systematically manipulating the height of the rims and the depth of the troughs and tracking the effects on ice wedge temperature. The results indicate that winter temperatures in the ice wedge are sensitive to both rim height and trough depth, but more sensitive to rim height. Rims act as preferential outlets of subsurface heat; increasing rim size decreases winter temperatures in the ice wedge. Deeper troughs lead to increased snow entrapment, promoting insulation of the ice wedge. The potential for ice wedge cracking is therefore reduced if rims are destroyed or if troughs subside, due to warmer conditions in the ice wedge. These findings can help explain the origins of secondary ice wedges in modern and ancient polygons. The findings also imply that the potential for re-establishing rims in modern thermokarst-affected terrain will be limited by reduced cracking activity in the ice wedges, even if regional air temperatures stabilize.

Active layer thickness (ALT), the uppermost layer of soil that thaws on an annual basis, is a direct control on the amount of organic carbon potentially available for decomposition and release to the atmosphere as carbon‐rich Arctic permafrost soils thaw in a warming climate. We investigate how key site characteristics affect ALT using an integrated surface/subsurface permafrost thermal hydrology model. ALT is most sensitive to organic layer thickness followed by snow depth but is relatively insensitive to the amount of water on the landscape with other conditions held fixed. The weak ALT sensitivity to subsurface saturation suggests that changes in Arctic landscape hydrology may only have a minor effect on future ALT. However, surface inundation amplifies the sensitivities to the other parameters and under large snowpacks can trigger the formation of near‐surface taliks.

Identification of the pumping influences at monitoring wells caused by spatially and temporally variable water supply pumping can be a challenging, yet an important hydrogeological task. The information that can be obtained can be critical for conceptualization of the hydrogeological conditions and indications of the zone of influence of the individual pumping wells. However, the pumping influences are often intermittent and small in magnitude with variable production rates from multiple pumping wells. While these difficulties may support an inclination to abandon the existing dataset and conduct a dedicated cross-hole pumping test, that option can be challenging and expensive to coordinate and execute. This paper presents a method that utilizes a simple analytical modeling approach for analysis of a long-term water level record utilizing an inverse modeling approach. The methodology allows the identification of pumping wells influencing the water level fluctuations. Thus, the analysis provides an efficient and cost-effective alternative to designed and coordinated cross-hole pumping tests. We apply this method on a dataset from the Los Alamos National Laboratory site. Our analysis also provides (1) an evaluation of the information content of the transient water level data; (2) indications of potential structures of the aquifer heterogeneity inhibiting or promoting pressure propagation; and (3) guidance for the development of more complicated models requiring detailed specification of the aquifer heterogeneity.