Spatial and Temporal Variations in Active Layer Thawing and Their Implication on Runoff Generation in Peat-Covered Permafrost Terrain

The distribution of frost table depths on a peat-covered permafrost slope was examined in a discontinuous permafrost region in northern Canada over 4 consecutive years at a variety of spatial scales, to elucidate the role of active layer development on runoff generation. Frost table depths were highly variable over relatively short distances (0.25–1 m), and the spatial variability was strongly correlated to soil moisture distribution, which was partly influenced by lateral flow converging to frost table depressions. On an interannual basis, thaw rates were temporally correlated to air temperature and the amount of precipitation input. Simple simulations show that lateral subsurface flow is governed by the frost table topography having spatially variable storage that has to be filled before water can spill over to generate flow downslope, in a similar manner that bedrock topography controls subsurface flow. However, unlike the bedrock surface, the frost table is variable with time and strongly influenced by the heat transfer involving water. Therefore, it is important to understand the feedback between thawing and subsurface water flow and to properly represent the feedback in hydrological models of permafrost regions

Effects of Freezing on Soil Temperature, Freezing Front Propagation and Moisture Redistribution in Peat: Laboratory Investigations

There are not many studies that report water movement in freezing peat. Soil column studies under controlled laboratory settings can help isolate and understand the effects of different factors controlling freezing of the active layer in organic covered permafrost terrain. In this study, four peat Mesocosms were subjected to temperature gradients by bringing the Mesocosm tops in contact with subzero air temperature while maintaining a continuously frozen layer at the bottom (proxy permafrost). Soil water movement towards the freezing front (from warmer to colder regions) was inferred from soil freezing curves, liquid water content time series and from the total water content of frozen core samples collected at the end of freezing cycle. A substantial amount of water, enough to raise the upper surface of frozen saturated soil within 15 cm of the soil surface at the end of freezing period appeared to have moved upwards during freezing. Diffusion under moisture gradients and effects of temperature on soil matric potential, at least in the initial period, appear to drive such movement as seen from analysis of freezing curves. Freezing front (separation front between soil zones containing and free of ice) propagation is controlled by latent heat for a long time during freezing. A simple conceptual model describing freezing of an organic active layer initially resembling a variable moisture landscape is proposed based upon the results of this study. The results of this study will help in understanding, and ultimately forecasting, the hydrologic response of wetland-dominated terrain underlain by discontinuous permafrost

Examining the Effect of Pore Size Distribution and Shape on Flow through Unsaturated Peat using Computer Tomography

The hydraulic conductivity of unsaturated peat soil is controlled by the air-filled porosity, pore size and geometric distribution as well as other physical properties of peat materials. This study investigates how the size and shape of pores affects the flow of water through peat soils. In this study we used X-ray Computed Tomography (CT), at 45μm resolution under 5 specific soil-water pressure head levels to provide 3-D, high-resolution images that were used to detect the inner pore structure of peat samples under a changing water regime. Pore structure and configuration were found to be irregular, which affected the rate of water transmission through peat soils. The 3-D analysis suggested that pore distribution is dominated by a single large pore-space. At low pressure head, this single large air-filled pore imparted a more effective flowpath compared to smaller pores. Smaller pores were disconnected and the flowpath was more tortuous than in the single large air-filled pore, and their contribution to flow was negligible when the single large pore was active. We quantify the pore structure of peat soil that affects the hydraulic conductivity in the unsaturated condition, and demonstrate the validity of our estimation of peat unsaturated hydraulic conductivity by making a comparison with a standard permeameter-based method. Estimates of unsaturated hydraulic conductivities were made for the purpose of testing the sensitivity of pore shape and geometry parameters on the hydraulic properties of peats and how to evaluate the structure of the peat and its affects on parameterization. We also studied the ability to quantify these factors for different soil moisture contents in order to define how the factors controlling the shape coefficient vary with changes in soil water pressure head. The relation between measured and estimated unsaturated hydraulic conductivity at various heads shows that rapid initial drainage, that changes the air-filled pore properties, creates a sharp decline in hydraulic conductivity. This is because the large pores readily lose water, the peat rapidly becomes less conductive and the flow path among pores, more tortuous

Comparison of Algorithms and Parameterisations for Infiltration into Organic-Covered Permafrost Soils

Infiltration into frozen and unfrozen soils is critical in hydrology, controlling active layer soil water dynamics and influencing runoff. Few Land Surface Models (LSMs) and Hydrological Models (HMs) have been developed, adapted or tested for frozen conditions and permafrost soils. Considering the vast geographical area influenced by freeze/thaw processes and permafrost, and the rapid environmental change observed worldwide in these regions, a need exists to improve models to better represent their hydrology. In this study, various infiltration algorithms and parameterisation methods, which are commonly employed in current LSMs and HMs were tested against detailed measurements at three sites in Canada’s discontinuous permafrost region with organic soil depths ranging from 0.02 to 3 m. Field data from two consecutive years were used to calibrate and evaluate the infiltration algorithms and parameterisations. Important conclusions include: (1) the single most important factor that controls the infiltration at permafrost sites is ground thaw depth, (2) differences among the simulated infiltration by different algorithms and parameterisations were only found when the ground was frozen or during the initial fast thawing stages, but not after ground thaw reaches a critical depth of 15 to 30 cm, (3) despite similarities in simulated total infiltration after ground thaw reaches the critical depth, the choice of algorithm influenced the distribution of water among the soil layers, and (4) the ice impedance factor for hydraulic conductivity, which is commonly used in LSMs and HMs, may not be necessary once the water potential driven frozen soil parameterisation is employed. Results from this work provide guidelines that can be directly implemented in LSMs and HMs to improve their application in organic covered permafrost soils

The positive net radiative greenhouse gas forcing of increasing methane emissions from a thawing boreal forest-wetland landscape

At the southern margin of permafrost in North America, climate change causes widespread permafrost thaw. In boreal lowlands, thawing forested permafrost peat plateaus (‘forest’) lead to expansion of permafrost‐free wetlands (‘wetland’). Expanding wetland area with saturated and warmer organic soils is expected to increase landscape methane (CH4) emissions. Here, we quantify the thaw‐induced increase in CH4 emissions for a boreal forest‐wetland landscape in the southern Taiga Plains, Canada, and evaluate its impact on net radiative forcing relative to potential long‐term net carbon dioxide (CO2) exchange. Using nested wetland and landscape eddy covariance net CH4 flux measurements in combination with flux footprint modeling, we find that landscape CH4 emissions increase with increasing wetland‐to‐forest ratio. Landscape CH4 emissions are most sensitive to this ratio during peak emission periods, when wetland soils are up to 10 °C warmer than forest soils. The cumulative growing season (May–October) wetland CH4 emission of ~13 g CH4 m−2 is the dominating contribution to the landscape CH4 emission of ~7 g CH4 m−2. In contrast, forest contributions to landscape CH4 emissions appear to be negligible. The rapid wetland expansion of 0.26 ± 0.05% yr−1 in this region causes an estimated growing season increase of 0.034 ± 0.007 g CH4 m−2 yr−1 in landscape CH4 emissions. A long‐term net CO2 uptake of >200 g CO2 m−2 yr−1 is required to offset the positive radiative forcing of increasing CH4 emissions until the end of the 21st century as indicated by an atmospheric CH4 and CO2 concentration model. However, long‐term apparent carbon accumulation rates in similar boreal forest‐wetland landscapes and eddy covariance landscape net CO2 flux measurements suggest a long‐term net CO2 uptake between 49 and 157 g CO2 m−2 yr−1. Thus, thaw‐induced CH4 emission increases likely exert a positive net radiative greenhouse gas forcing through the 21st century

Evaluation of the Algorithms and Parameterizations for Ground Thawing and Freezing Simulation in Permafrost Regions

Ground thawing and freezing depths (GTFDs) strongly influence the hydrology and energy balances of permafrost regions. Current methods to simulate GTFD differ in algorithm type, soil parameterization, representation of latent heat, and unfrozen water content. In this study, five algorithms (one semiempirical, two analytical, and two numerical), three soil thermal conductivity parameterizations, and three unfrozen water parameterizations were evaluated against detailed field measurements at four field sites in Canada’s discontinuous permafrost region. Key findings include: (1) de Vries’ parameterization is recommended to determine the thermal conductivity in permafrost soils; (2) the three unfrozen water parameterization methods exhibited little difference in terms of GTFD simulations, yet the segmented linear function is the simplest to be implemented; (3) the semiempirical algorithm reasonably simulates thawing at permafrost sites and freezing at seasonal frost sites with site-specific calibration. However, large interannual and intersite variations in calibration coefficients limit its applicability for dynamic analysis; (4) when driven by surface forcing, analytical algorithms performed marginally better than the semiempirical algorithm. The inclusion of bottom forcing improved analytical algorithm performance, yet their results were still poor compared with those achieved by numerical algorithms; (5) when supplied with the optimal inputs, soil parameterizations, and model configurations, the numerical algorithm with latent heat treated as an apparent heat capacity achieved the best GTFD simulations among all algorithms at all sites. Replacing the observed bottom temperature with a zero heat flux boundary condition did not significantly reduce simulation accuracy, while assuming a saturated profile caused large errors at several sites

Structure of peat soils and implications for water storage, flow and solute transport: A review update for geochemists

Global peatlands are a valuable but vulnerable resource. They represent a significant carbon and energy reservoir and play major roles in water and biogeochemical cycles. Peat soils are highly complex porous media with distinct characteristic physical and hydraulic properties. Pore sizes in undecomposed peat can exceed 5 mm, but significant shrinkage occurs during dewatering, compression and decomposition, reducing pore-sizes. The structure of peat soil consists of pores that are open and connected, dead-ended or isolated. The resulting dual-porosity nature of peat soils affects water flow and solute migration, which influence reactive transport processes and biogeochemical functions. Advective movement of aqueous and colloidal species is restricted to the hydrologically active (or mobile) fraction of the total porosity, i.e. the open and connected pores. Peat may attenuate solute migration through molecular diffusion into the closed and dead-end pores, and for reactive species, also through sorption and degradation reactions. Slow, diffusion-limited solute exchanges between the mobile and immobile regions may give rise to pore-scale chemical gradients and heterogeneous distributions of microbial habitats and activity in peat soils. While new information on the diversity and functionalities of peat microbial communities is rapidly accumulating, the significance of the geochemical and geomicrobial study on peat stands to benefit from a basic understanding of the physical structure of peat soils. In this paper, we review the current knowledge of key physical and hydraulic properties related to the structure of globally available peat soils and briefly discuss their implications for water storage, flow and the migration of solutes. This paper is intended to narrow the gap between the ecohydrological and biogeochemical research communities working on peat soils.Canada Excellence Research Chair (CERC) progra

Permafrost Landscape History Shapes Fluvial Chemistry, Ecosystem Carbon Balance, and Potential Trajectories of Future Change

Intensifying permafrost thaw alters carbon cycling by mobilizing large amounts of terrestrial substrate into aquatic ecosystems. Yet, few studies have measured aquatic carbon fluxes and constrained drivers of ecosystem carbon balance across heterogeneous Arctic landscapes. Here, we characterized hydrochemical and landscape controls on fluvial carbon cycling, quantified fluvial carbon fluxes, and estimated fluvial contributions to ecosystem carbon balance across 33 watersheds in four ecoregions in the continuous permafrost zone of the western Canadian Arctic: unglaciated uplands, ice-rich moraine, and organic-rich lowlands and till plains. Major ions, stable isotopes, and carbon speciation and fluxes revealed patterns in carbon cycling across ecoregions defined by terrain relief and accumulation of organics. In previously unglaciated mountainous watersheds, bicarbonate dominated carbon export (70% of total) due to chemical weathering of bedrock. In lowland watersheds, where soil organic carbon stores were largest, lateral transport of dissolved organic carbon (50%) and efflux of biotic CO2 (25%) dominated. In watersheds affected by thaw-induced mass wasting, erosion of ice-rich tills enhanced chemical weathering and increased particulate carbon fluxes by two orders of magnitude. From an ecosystem carbon balance perspective, fluvial carbon export in watersheds not affected by thaw-induced wasting was, on average, equivalent to 6%–16% of estimated net ecosystem exchange (NEE). In watersheds affected by thaw-induced wasting, fluvial carbon export approached 60% of NEE. Because future intensification of thermokarst activity will amplify fluvial carbon export, determining the fate of carbon across diverse northern landscapes is a priority for constraining trajectories of permafrost region ecosystem carbon balance

Functionally compromised CHD7 alleles in patients with isolated GnRH deficiency

Inactivating mutations in chromodomain helicase DNA binding protein 7 (CHD7) cause CHARGE syndrome, a severe multiorgan system disorder of which Isolated gonadotropin-releasing hormone (GnRH) deficiency (IGD) is a minor feature. Recent reports have described predominantly missense CHD7 alleles in IGD patients, but it is unclear if these alleles are relevant to causality or overall genetic burden of Kallmann syndrome (KS) and normosmic form of IGD. To address this question, we sequenced CHD7 in 783 well-phenotyped IGD patients lacking full CHARGE features; we identified nonsynonymous rare sequence variants in 5.2% of the IGD cohort (73% missense and 27% splice variants). Functional analyses in zebrafish using a surrogate otolith assay of a representative set of these CHD7 alleles showed that rare sequence variants observed in controls showed no altered function. In contrast, 75% of the IGD-associated alleles were deleterious and resulted in both KS and normosmic IGD. In two families, pathogenic mutations in CHD7 coexisted with mutations in other known IGD genes. Taken together, our data suggest that rare deleterious CHD7 alleles contribute to the mutational burden of patients with both KS and normosmic forms of IGD in the absence of full CHARGE syndrome. These findings (i) implicate a unique role or preferential sensitivity for CHD7 in the ontogeny of GnRH neurons, (ii) reiterate the emerging genetic complexity of this family of IGD disorders, and (iii) demonstrate how the coordinated use of well-phenotyped cohorts, families, and functional studies can inform genetic architecture and provide insights into the developmental biology of cellular systems