Humidity is a basic and crucial meteorological indicator commonly measured in several forms, including specific humidity, relative humidity, and absolute humidity. These different forms can be inter-derived based on the saturation vapor pressure (SVP). In past decades, dozens of formulae have been developed to calculate the SVP with respect to, and in equilibrium with, liquid water and solid ice surfaces, but many prior studies use a single function for all temperature ranges, without considering the distinction between over the liquid water and ice surfaces. These different approaches can result in humidity estimates that may impact our understanding of surface-subsurface thermal-hydrological dynamics in cold regions. In this study, we compared the relative humidity (RH) downloaded and calculated from four data sources in Alaska based on five commonly used SVP formulas. These RHs, along with other meteorological indicators, were then used to drive physics-rich land surface models at a permafrost-affected site. We found that higher values of RH (up to 40 %) were obtained if the SVP was calculated with the over-ice formulation when air temperatures were below freezing, which could lead to a 30 % maximum difference in snow depths. The choice of whether to separately calculate the SVP over an ice surface in winter also produced a significant range (up to 0.2 m) in simulated annual maximum thaw depths. The sensitivity of seasonal thaw depth to the formulation of SVP increases with the rainfall rate and the height of above-ground ponded water, while it diminishes with warmer air temperatures. These results show that RH variations based on the calculation of SVP with or without over-ice calculation meaningfully impact physicallybased predictions of snow depth, sublimation, soil temperature, and active layer thickness. Under particular conditions, when severe flooding (inundation) and cool air temperatures are present, care should be taken to evaluate how humidity data is estimated for land surface and earth system modeling

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

Arctic and boreal permafrost soil organic carbon (SOC) decomposition has been slower than carbon inputs from plant growth since the last glaciation. Anthropogenic climate warming has threatened this historical trend by accelerating SOC decomposition and altering wildfire regimes. We accurately modeled observed plant biomass and carbon emissions from wildfires in Alaskan ecosystems under current climate conditions. In projections to 2300 under the RCP8.5 climate scenario, we found that warming and increased atmospheric CO2 will result in plant biomass gains and higher litterfall. However, increased carbon losses from (a) wildfire combustion and (b) rapid SOC decomposition driven by increased deciduous litter production, root exudation, and active layer depth will lead to about 4.4 PgC of soil carbon losses from Alaska by 2300 and most (88%) of these loses will be from the top 1 m of soil. These SOC losses offset plant carbon gains, causing the ecosystem to transition to a net carbon source after 2200. Simulations excluding wildfire increases yielded about a factor of four lower SOC losses by 2300. Our results show that projected wildfire and its direct and indirect effects on plant and soil carbon may accelerate high-latitude soil carbon losses, resulting in a positive feedback to climate change.

Northern peatlands are a major component of the global carbon (C) cycle. Widespread climate-driven ecohydrological changes in these ecosystems can have major consequences on their C sequestration function. Here, we synthesize plant macrofossil data from 33 surficial peat cores from different ecoclimatic regions, with high-resolution chronologies. The main objectives were to document recent ecosystem state shifts and explore their impact on C sequestration in high-latitude undisturbed peatlands of northeastern Canada. Our synthesis shows widespread recent ecosystem shifts in peatlands, such as transitions from oligotrophic fens to bogs and Sphagnum expansion, coinciding with climate warming which has also influenced C accumulation during the last similar to 100 years. The rapid shifts towards drier bog communities and an expansion of Sphagnum sect. Acutifolia after 1980 CE were most pronounced in the northern subarctic sites and are concurrent with summer warming in northeastern Canada. These results provide further evidence of a northward migration of Sphagnum-dominated peatlands in North America in response to climate change. The results also highlight differences in the timing of ecosystem shifts among peatlands and regions, reflecting internal peatland dynamics and varying responses of vegetation communities. Our study suggests that the recent rapid climate-driven shifts from oligotrophic fen to drier bog communities have promoted plant productivity and thus peat C accumulation. We highlight the importance of considering recent ecohydrological trajectories when modelling the potential contribution of peatlands to climate change. Our study suggests that, contrary to expectations, peat C sequestration could be promoted in high-latitude non-permafrost peatlands where wet sedge fens may transition to drier Sphagnum bog communities due to warmer and longer growing seasons.