Subsurface flow measurements using passive flux meters in variably-saturated cold-regions landscapes

To date, passive flux meters have predominantly been applied in temperate environments for tracking the movement of contaminants in groundwater. This study applies these instruments to reduce uncertainty in (typically instantaneous) flux measurements made in a low-gradient, wetland dominated, discontinuous permafrost environment. This method supports improved estimation of unsaturated and over-winter subsurface flows which are very difficult to quantify using hydraulic gradient-based approaches. Improved subsurface flow estimates can play a key role in understanding the water budget of this landscape.     

Mercury and methylmercury biogeochemistry in a thawing permafrost wetland complex,Northwest Territories,Canada

In arctic and sub‐arctic environments, mercury (Hg), more specifically toxic methylmercury (MeHg), is of growing concern to local communities because of its accumulation in fish. In these regions, there is particular interest in the potential mobilization of atmospherically deposited Hg sequestered in permafrost that is thawing at unprecedented rates. Permafrost thaw and the resulting ground surface subsidence transforms forested peat plateaus into treeless and permafrost‐free thermokarst wetlands where inorganic Hg released from the thawed permafrost and draining from the surrounding peat plateaus may be transformed to MeHg. This study begins to characterize the spatial distribution of MeHg in a peat plateau–thermokarst wetland complex, a feature that prevails throughout the wetland‐dominated southern margin of thawing discontinuous permafrost in Canada's Northwest Territories. We measured pore water total Hg, MeHg, dissolved organic matter characteristics and general water chemistry parameters to evaluate the role of permafrost thaw on the pattern of water chemistry. A gradient in vegetation composition, water chemistry and dissolved organic matter characteristics followed a toposequence from the ombrotrophic bogs near the crest of the complex to poor fens at its downslope margins. We found that pore waters in poor fens contained elevated levels of MeHg, and the water draining from these features had dissolved MeHg concentrations 4.5 to 14.5 times higher than the water draining from the bogs. It was determined through analysis of historical aerial images that the poor fens in the toposequence had formed relatively recently (early 1970s) as a result of permafrost thaw. Differences between the fens and bogs are likely to be a result of their differences in groundwater function, and this suggests that permafrost thaw in this landscape can result in hotspots for Hg methylation that are hydrologically connected to downstream ecosystems. Copyright © 2016 John Wiley & Sons, Ltd.     

Influence of pore size and geometry on peat unsaturated hydraulic conductivity computed from 3D computed tomography image analysis

In organic soils, hydraulic conductivity is related to the degree of decomposition and soil compression, which reduce the effective pore diameter and consequently restrict water flow. This study investigates how the size distribution and geometry of air‐filled pores control the unsaturated hydraulic conductivity of peat soils using high‐resolution (45 µm) three‐dimensional (3D) X‐ray computed tomography (CT) and digital image processing of four peat sub‐samples from varying depths under a constant soil water pressure head. Pore structure and configuration in peat were found to be irregular, with volume and cross‐sectional area showing fractal behaviour that suggests pores having smaller values of the fractal dimension in deeper, more decomposed peat, have higher tortuosity and lower connectivity, which influences hydraulic conductivity. The image analysis showed that the large reduction of unsaturated hydraulic conductivity with depth is essentially controlled by air‐filled pore hydraulic radius, tortuosity, air‐filled pore density and the fractal dimension due to degree of decomposition and compression of the organic matter. The comparisons between unsaturated hydraulic conductivity computed from the air‐filled pore size and geometric distribution showed satisfactory agreement with direct measurements using the permeameter method. This understanding is important in characterizing peat properties and its heterogeneity for monitoring the progress of complex flow processes at the field scale in peatlands. Copyright © 2010 John Wiley & Sons, Ltd.     

Permafrost‐thaw‐induced land‐cover change in the Canadian subarctic: implications for water resources

Climate warming and human disturbance in north‐western Canada have been accompanied by degradation of permafrost, which introduces considerable uncertainty to the future availability of northern freshwater resources. This study demonstrates the rate and spatial pattern of permafrost loss in a region that typifies the southern boundary of permafrost. Remote‐sensing analysis of a 1·0 km² area indicates that permafrost occupied 0·70 km² in 1947 and decreased with time to 0·43 km² by 2008. Ground‐based measurements demonstrate the importance of horizontal heat flows in thawing discontinuous permafrost, and show that such thaw produces dramatic land‐cover changes that can alter basin runoff production in this region. A major challenge to northern water resources management in the twenty‐first century therefore lies in predicting stream flows dynamically in the context of widely occurring permafrost thaw. The need for appropriate water resource planning, mitigation, and adaptation strategies is explained. Copyright © 2010 John Wiley & Sons, Ltd.     

Sensitivity analysis of EUROSEM using Monte Carlo simulation I: hydrological,soil and vegetation parameters

Knowledge about model uncertainty is essential for erosion modelling and provides important information when it comes to parameterizing models. In this paper a sensitivity analysis of the European soil erosion model (EUROSEM) is carried out using Monte Carlo simulation, suitable for complex non‐linear models, using time‐dependent driving variables. The analysis revealed some important characteristics of the model. The variability of the static output parameters was generally high, with the hydrologic parameters being the most important ones, especially saturated hydraulic conductivity and net capillary drive followed by the percentage basal area for the hydrological and vegetation parameters and detachability and cohesion for the soil erosion parameters. Overall, sensitivity to vegetation parameters was insignificant. The coefficient of variation for the sedigraph was higher than for the hydrograph, especially from the beginning of the rainstorm and up to the peak, and may explain difficulties encountered when trying to match simulated hydrographs and sedigraphs with observed ones. The findings from this Monte Carlo simulation calls for improved within‐storm modelling of erosion processes in EUROSEM. Information about model uncertainty will be incorporated in a new EUROSEM user interface. Copyright © 2000 John Wiley & Sons, Ltd.     

Transient and Transition Factors in Modeling Permafrost Thaw and Groundwater Flow

Permafrost covers approximately 24% of the Northern Hemisphere, and much of it is degrading, which causes infrastructure failures and ecosystem transitions. Understanding groundwater and heat flow processes in permafrost environments is challenging due to spatially and temporarily varying hydraulic connections between water above and below the near-surface discontinuous frozen zone. To characterize the transitional period of permafrost degradation, a three-dimensional model of a permafrost plateau that includes the supra-permafrost zone and surrounding wetlands was developed. The model is based on the Scotty Creek basin in the Northwest Territories, Canada. FEFLOW groundwater flow and heat transport modeling software is used in conjunction with the piFreeze plug-in, to account for phase changes between ice and water. The Simultaneous Heat and Water (SHAW) flow model is used to calculate ground temperatures and surface water balance, which are then used as FEFLOW boundary conditions. As simulating actual permafrost evolution would require hundreds of years of climate variations over an evolving landscape, whose geomorphic features are unknown, methodologies for developing permafrost initial conditions for transient simulations were investigated. It was found that a model initialized with a transient spin-up methodology, that includes an unfrozen layer between the permafrost table and ground surface, yields better results than with steady-state permafrost initial conditions. This study also demonstrates the critical role that variations in land surface and permafrost table microtopography, along with talik development, play in permafrost degradation. Modeling permafrost dynamics will allow for the testing of remedial measures to stabilize permafrost in high value infrastructure environments.     

Ultra‐rapid topographic surveying for complex environments: the hand‐held mobile laser scanner (HMLS)

Terrestrial laser scanning is the current technique of choice for acquiring high resolution topographic data at the site scale (i.e. over tens to hundreds of metres), for accurate volume measurements or process modelling. However, in regions of complex topography with multiple local horizons, restricted lines of sight significantly hinder use of such tripod‐based instruments by requiring multiple setups to achieve full coverage of the area. We demonstrate a novel hand‐held mobile laser scanning technique that offers particular promise for site‐scale topographic surveys of complex environments. To carry out a survey, the hand‐held mobile laser scanner (HMLS) is walked across a site, mapping around the surveyor continuously en route. We assess the accuracy of HMLS data by comparing survey results from an eroding coastal cliff site with those acquired by a state‐of‐the‐art terrestrial laser scanner (TLS) and also with the results of a photo‐survey, processed by structure from motion and multi‐view stereo (SfM‐MVS) algorithms. HMLS data are shown to have a root mean square (RMS) difference to the benchmark TLS data of 20 mm, not dissimilar to that of the SfM‐MVS survey (18 mm). The efficiency of the HMLS system in complex terrain is demonstrated by acquiring topographic data covering ~780 m² of salt‐marsh gullies, with a mean point spacing of 4.4 cm, in approximately six minutes. We estimate that HMLS surveying of gullies is approximately 40 times faster than using a TLS and six times faster than using SfM‐MVS. © 2013 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Effects of soil compaction,rain exposure and their interaction on soil carbon dioxide emission

Soils release more carbon, primarily as carbon dioxide (CO₂), per annum than current global anthropogenic emissions. Soils emit CO₂ through mineralization and decomposition of organic matter and respiration of roots and soil organisms. Given this, the evaluation of the effects of abiotic factors on microbial activity is of major importance when considering the mitigation of greenhouse gases emissions. Previous studies demonstrate that soil CO₂ emission is significantly affected by temperature and soil water content. A limited number of studies have illustrated the importance of bulk density and soil surface characteristics as a result of exposure to rain on CO₂ emission, however, none examine their relative importance. Therefore, this study investigated the effects of soil compaction and exposure of the soil surface to rainfall and their interaction on CO₂ release. We conducted a factorial laboratory experiment with three soil types after sieving (clay, silt and sand soil), three different bulk densities (1·1 g cm^–3, 1·3 g cm^–3, 1·5 g cm^–3) and three different exposures to rainfall (no rain, 30 minutes and 90 minutes of rainfall). The results demonstrated CO₂ release varied significantly with bulk density, exposure to rain and time. The relationship between rain exposure and CO₂ is positive: CO₂ emission was 53% and 42% greater for the 90 minutes and 30 minutes rainfall exposure, respectively, compared to those not exposed to rain. Bulk density exhibited a negative relationship with CO₂ emission: soil compacted to a bulk density of 1·1 g cm^–3 emitted 32% more CO₂ than soil compacted to 1·5 g cm^–3. Furthermore we found that the magnitude of CO₂ effluxes depended on the interaction of these two abiotic factors. Given these results, understanding the influence of soil compaction and raindrop impact on CO₂ emission could lead to modified soil management practices which promote carbon sequestration. Copyright © 2012 John Wiley & Sons, Ltd.     

Coherent structures and flux contribution over an inhomogeneously irrigated cotton field

The turbulence data measured at two levels (i.e., 8.7 and 2.7?m) in the Energy Balance Experiment (EBEX), which was conducted in San Joaquin Valley in California during the period from July 20 to August 24, 2000, are used to study the characteristics of coherent structures over an irrigated cotton field. Patch-to-patch irrigation in the field generated the dry-to-wet horizontal advection and the oasis effects, leading to the development of a stably internal boundary layer (SIBL) in the late mornings or the early afternoons. The SIBL persisted in the rest of the afternoons. Under this circumstance, a near-neutral atmospheric surface layer (ASL) developed during the period with a stratification transition from the unstable to stable conditions during the daytime. Therefore, EBEX provides us with unique datasets to investigate the features of coherent structures that were generated over the patches upstream and passed by our site in the unstable ASL, the near-neutral ASL, and the SIBL. We use an objective detection technique and the conditional average method that is developed based on the wavelet analysis. Our data reveal some consistencies and inconsistencies in the characteristics of coherent structures as compared with previous studies. Ramp-like structures and sweep?Cejection cycles under the daytime SIBL have similar patterns to those under the nocturnal stable ASL. However, some features (i.e., intermittence) are different from those under the nocturnal stable ASL. Under the three stratifications, thermal and mechanical factors in the ASL perform differently in affecting the ramp intensity for different quantities (i.e., velocity components, temperature, and specific humidity), leading to coherent structures that modulate turbulence flow and alter turbulent transfer differently. It is also found that coherent structures contribute about 10?C20% to the total fluxes in our case with different flux contributions under three stratifications and with higher transporting efficiency in sensible heat flux than latent heat and momentum fluxes.