Processes In Cold Regions Research Articles

Spatially distributed tracer-aided runoff modelling and dynamics of storage and water ages in a permafrost-influenced catchment

Abstract. Permafrost strongly controls hydrological processes in cold regions. Our understanding of how changes in seasonal and perennial frozen ground disposition and linked storage dynamics affect runoff generation processes remains limited. Storage dynamics and water redistribution are influenced by the seasonal variability and spatial heterogeneity of frozen ground, snow accumulation and melt. Stable isotopes are potentially useful for quantifying the dynamics of water sources, flow paths and ages, yet few studies have employed isotope data in permafrost-influenced catchments. Here, we applied the conceptual model STARR (the Spatially distributed Tracer-Aided Rainfall–Runoff model), which facilitates fully distributed simulations of hydrological storage dynamics and runoff processes, isotopic composition and water ages. We adapted this model for a subarctic catchment in Yukon Territory, Canada, with a time-variable implementation of field capacity to include the influence of thaw dynamics. A multi-criteria calibration based on stream flow, snow water equivalent and isotopes was applied to 3 years of data. The integration of isotope data in the spatially distributed model provided the basis for quantifying spatio-temporal dynamics of water storage and ages, emphasizing the importance of thaw layer dynamics in mixing and damping the melt signal. By using the model conceptualization of spatially and temporally variable storage, this study demonstrates the ability of tracer-aided modelling to capture thaw layer dynamics that cause mixing and damping of the isotopic melt signal.

Evaluation of calculation models for the unfrozen water content of freezing soils

Evaluation of calculation models for the unfrozen water content of freezing soil is very important for numerical simulations. Based on the published 250 laboratory data and around 5 years field observed data, this paper evaluated the five widely used unfrozen water content models, i.e., the Zhang et al. (2017a) model, the Anderson and Tice (1972) model, the Michalowski (1993) model, the Mckenzie et al. (2007) model, and the Kozlowski (2007) model. The results show that the five models can all be used to calculate the unfrozen water content of laboratory soil samples. The Anderson and Tice (1972) model, the Michalowski (1993) model, and the Kozlowski (2007) model are descripted by gravimetric unfrozen water content. However, it is difficult to directly measure the gravimetric unfrozen water content during the freezing process. The Mckenzie et al. (2007) model and the Zhang et al. (2017a) model are descripted by volumetric unfrozen water content which can be directly measured during the freezing process. Besides, the parameters of the Zhang et al. (2017a) model are relatively simple, although the model slightly overestimates the volumetric unfrozen water content for some laboratory soils. Therefore, the Zhang et al. (2017a) model is more convenient for numerical simulations. Furthermore, all the models have great potential to calculate unfrozen water content of freezing soils in field observations. This paper would provide significant reference for the selection of unfrozen water content models in the description of hydro-thermo coupling process in cold regions.

IRASSI infrared space interferometer: Formation geometry and relative dynamics analysis

Space-based interferometry has gained prominence in recent years, largely because higher spatial resolutions of celestial observations can be achieved with multi-telescope formations compared to those achieved with a fixed-aperture, single telescope. IRASSI is a space interferometer composed of five spacecraft, whose aim is to observe particular chemical and physical processes in cold regions of space, such as dust clouds and stellar disks, in the far-infrared frequencies. Ultimately, the goal is to study the genesis of planets, star formation and evolution processes in these cold regions and to understand how prebiotic conditions in Earth-like planets are created. IRASSI will orbit the second Lagrange point, L2, of Sun-Earth/Moon system. The operating principle of IRASSI is based on free-drifting baselines, which dynamically change during the observations and measure therefore the incoming wavefront of a celestial target at different locations in space. This process relies on very accurate measurements of the baselines - at micrometer level - rather than on precise control of the formation. Naturally, a free-flying formation comes with a set of challenges, namely identifying a nominal formation geometry, that is, a suitable dispersion of the telescopes in three-dimensional space. In addition, understanding how this free-drifting geometry is expected to change is crucial, particularly if this may affect the operation of the telescope instruments and thus the quality of the final synthesized images. The present paper describes therefore the major requirements for establishing a desired formation geometry and proposes a preliminary nominal formation for the observations. The relative dynamics of the free-drifting spacecraft are modeled and evaluated. Final considerations regarding formation control are presented and the paper concludes with a summary of the work and outlook for the future.

Infiltration into frozen soil: From core‐scale dynamics to hillslope‐scale connectivity

AbstractInfiltration into frozen soil is a key hydrological process in cold regions. Although the mechanisms behind point‐scale infiltration into frozen soil are relatively well understood, questions remain about upscaling point‐scale results to estimate hillslope‐scale run‐off generation. Here, we tackle this question by combining laboratory, field, and modelling experiments. Six large (0.30‐m diameter by 0.35‐m deep) soil cores were extracted from an experimental hillslope on the Canadian Prairies. In the laboratory, we measured run‐off and infiltration rates of the cores for two antecedent moisture conditions under snowmelt rates and diurnal freeze–thaw conditions observed on the same hillslope. We combined the infiltration data with spatially variable data from the hillslope, to parameterise a surface run‐off redistribution model. We used the model to determine how spatial patterns of soil water content, snowpack water equivalent (SWE), and snowmelt rates affect the spatial variability of infiltration and hydrological connectivity over frozen soil. Our experiments showed that antecedent moisture conditions of the frozen soil affected infiltration rates by limiting the initial soil storage capacity and infiltration front penetration depth. However, shallow depths of infiltration and refreezing created saturated conditions at the surface for dry and wet antecedent conditions, resulting in similar final infiltration rates (0.3 mm hr−1). On the hillslope‐scale, the spatial variability of snowmelt rates controlled the development of hydrological connectivity during the 2014 spring melt, whereas SWE and antecedent soil moisture were unimportant. Geostatistical analysis showed that this was because SWE variability and antecedent moisture variability occurred at distances shorter than that of topographic variability, whereas melt variability occurred at distances longer than that of topographic variability. The importance of spatial controls will shift for differing locations and winter conditions. Overall, our results suggest that run‐off connectivity is determined by (a) a pre‐fill phase, during which a thin surface soil layer wets up, refreezes, and saturates, before infiltration excess run‐off is generated and (b) a subsequent fill‐and‐spill phase on the surface that drives hillslope‐scale run‐off.

Spatiotemporal Changes in Active Layer Thickness under Contemporary and Projected Climate in the Northern Hemisphere

Abstract Variability of active layer thickness (ALT) in permafrost regions is critical for assessments of climate change, water resources, and engineering applications. Detailed knowledge of ALT variations is also important for studies on ecosystem, hydrological, and geomorphological processes in cold regions. The primary objective of this study is therefore to provide a comprehensive 1971–2000 climatology of ALT and its changes across the entire Northern Hemisphere from 1850 through 2100. To accomplish this, in situ observations, the Stefan solution based on a thawing index, and the edaphic factor (E factor) are employed to calculate ALT. The thawing index is derived from (i) the multimodel ensemble mean of 16 models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) over 1850–2005, (ii) three representative concentration pathways (RCP2.6, RCP4.5, and RCP8.5) for 2006–2100, and (iii) Climatic Research Unit (CRU) gridded observations for 1901–2014. The results show significant spatial variability in in situ ALT that generally ranges from 40 to 320 cm, with some extreme values of 900 cm in the Alps. The differences in the ALT climatology between the three RCPs and the historical experiments ranged from 0 to 200 cm. The biggest increases, of 120–200 cm, are on the Qinghai–Tibetan Plateau, while the smallest increases of less than 20 cm are in Alaska. Averaged over all permafrost regions, mean ALT from CMIP5 increased significantly at 0.57 ± 0.04 cm decade−1 during 1850–2005, while 2006–2100 projections show ALT increases of 0.77 ± 0.08 cm decade−1 for RCP2.6, 2.56 ± 0.07 cm decade−1 for RCP4.5, and 6.51 ± 0.07 cm decade−1 for RCP8.5.

Changes of soil thermal and hydraulic regimes in the Heihe River Basin.

Soil thermal and hydraulic regimes are critical factors influencing terrestrial processes in cold regions. Collection of field data from frozen ground has occurred at point scales, but limited data exist that characterize changes of soil thermal and hydraulic regimes at the scale of the whole Heihe River Basin. This study uses a long-term regional climate model coupled with land surface model to investigate the soil thermal and hydraulic regime changes at a large spatial scale. It also explores potential factors, including the climate and non-climate factors. Results show that there is significant variability in mean annual air temperature (MAAT) of about 0.47°C/decade during 1980-2013. A time series of area-averaged mean annual soil temperature (MAST) over the whole Heihe River Basin shows a significant increase between 0.25 and 0.36°C/decade during 1984-2013, with a net change of 0.9°C. A trend of increasing wetness is found in soil moisture. Frozen days (FD) decreased significantly both in seasonally frozen ground (SFG) regions and permafrost regions, with a net change between 7 and 13days during 1984-2013. Freezing index (FI) had a positive effect on FD, while thawing index (TI), MAAT, precipitation, and normalized difference vegetation index (NDVI) had a negative effect. These results are important to understand dynamic mechanisms of soil freeze/thaw cycles.

Development of WEP-COR model to simulate land surface water and energy budgets in a cold region

Abstract The Water and Energy transfer Processes in Cold Regions (WEP-COR) model is an improved version of the Water and Energy transfer Processes in Large basins (WEP-L) model that integrates a multi-layer frozen soil model to simulate the hydrological processes in cold regions and the heat fluxes at different depths of frozen soil. The temperature, water content, freezing depth of the soil, and daily discharge were simulated and compared with observations. The simulated and observed data were used to analyze the runoff flow components. The results showed that the WEP-COR model can effectively simulate the distributions of the soil temperature and water content. The average root mean squared errors of the temperature, unfrozen water content, total water content, and freezing depth of the soil were 1.21 °C, 0.035 cm3/cm3, 0.034 cm3/cm3, and 17.6 cm, respectively. The mean Nash–Sutcliffe efficiency and relative error of the daily discharge were 0.64 and 6.58%, respectively. Compared with the WEP-L model, the WEP-COR model simulated the discharge with higher accuracy, especially during the soil thawing period. This improvement was mainly due to the addition of the frozen soil mechanism. The WEP-COR model can provide support for agricultural and water resources management in cold regions.

Chapter 2 The Quaternary

Abstract The Quaternary is the youngest geological period, beginning 2.58 Ma ago and including the present day; it is therefore the only geological period that is continuously growing. During the first epoch of the Quaternary, the Pleistocene, extremely cold and warm conditions alternated, frequently over short periods of time. This resulted in processes currently only operating in cold (polar and high-mountain) environments extending to and affecting the mid-latitudes, including the currently densely populated areas of North America and Europe. In Britain every region has been affected by cold-region processes, which have produced unique sedimentary and geomorphological signatures. Hence, an intimate knowledge of these processes is of direct relevance to engineering geologists and anyone working with natural materials. This chapter reviews the state of the art of (a) the stratigraphic (nomenclatorial) framework of the Quaternary, (b) prominent concepts that are of direct relevance to understanding the detailed overviews in Chapters 3–5; and (c) key findings on the dynamics of these processes and their implications for engineering-geological questions and problems.

An integrated simulation and optimization approach for reducing CO2 emissions from on-site construction process in cold regions

Carbon dioxide (CO2) reduction has obtained worldwide attention due to the global warming effects and the mitigating means such as carbon tax. Among the CO2 emission sources, construction and building related activities consume a considerable amount of energy and produce a significant volume of greenhouse gases (GHGs). This paper proposes an integrated simulation and optimization approach to reducing CO2 emissions from the on-site construction process, which consists of three main steps: (i) construction process simulation, (ii) CO2 emission quantification, and (iii) labour allocation optimization. The construction process is simulated utilizing a discrete-event simulation (DES) platform, Simphony.NET, in order to quantify and analyze the amount of CO2 emissions produced by the construction process. Based on the simulation model, the crew sizes for labour-intensive activities are analyzed and optimized with respect to CO2 emissions, with a genetic algorithm (GA) employed as the method of optimization. The integrated simulation and optimization approach is capable of reducing CO2 emission through (i) the choice of the start-date of the construction, and (ii) labour allocation optimization. A case study in Edmonton, Canada, is used to demonstrate the effectiveness and applicability of the developed approach. Results indicate that the total measured CO2 emissions from on-site construction can be reduced by 21.7% in winter by optimizing labour allocation. This research thus proposes a generic methodology by which to measure and reduce CO2 emissions produced from construction process through labour allocation.

Hydrological responses to climatic variability in a cold agricultural region

Extended severe dry and wet periods are frequently observed in the northern continental climate of the Canadian Prairies. Prairie streamflow is mainly driven by spring snowmelt of the winter snowpack, whilst summer rainfall is an important control on evapotranspiration and thus seasonality affects the hydrological response to drought and wet periods in complex ways. A field-tested physically based model was used to investigate the influences of climatic variability on hydrological processes in this region. The model was set up to resolve agricultural fields and to include key cold regions processes. It was parameterized from local and regional measurements without calibration and run for the South Tobacco Creek basin in southern Manitoba, Canada. The model was tested against snow depth and streamflow observations at multiple scales and performed well enough to explore the impacts of wet and dry periods on hydrological processes governing the basin scale hydrological response. Four hydro-climatic patterns with distinctive climatic seasonality and runoff responses were identified from differing combinations of wet/dry winter and summer seasons. Water balance analyses of these patterns identified substantive multiyear subsurface soil moisture storage depletion during drought (2001–2005) and recharge during a subsequent wet period (2009–2011). The fractional percentage of heavy rainfall days was a useful metric to explain the contrasting runoff volumes between dry and wet summers. Finally, a comparison of modeling approaches highlights the importance of antecedent fall soil moisture, ice lens formation during the snowmelt period, and peak snow water equivalent in simulating snowmelt runoff.

Improving snow process modeling with satellite‐based estimation of near‐surface‐air‐temperature lapse rate

AbstractIn distributed hydrological modeling, surface air temperature (Tair) is of great importance in simulating cold region processes, while the near‐surface‐air‐temperature lapse rate (NLR) is crucial to prepare Tair (when interpolating Tair from site observations to model grids). In this study, a distributed biosphere hydrological model with improved snow physics (WEB‐DHM‐S) was rigorously evaluated in a typical cold, large river basin (e.g., the upper Yellow River basin), given a mean monthly NLRs. Based on the validated model, we have examined the influence of the NLR on the simulated snow processes and streamflows. We found that the NLR has a large effect on the simulated streamflows, with a maximum difference of greater than 24% among the various scenarios for NLRs considered. To supplement the insufficient number of monitoring sites for near‐surface‐air‐temperature at developing/undeveloped mountain regions, the nighttime Moderate Resolution Imaging Spectroradiometer land surface temperature is used as an alternative to derive the approximate NLR at a finer spatial scale (e.g., at different elevation bands, different land covers, different aspects, and different snow conditions). Using satellite‐based estimation of NLR, the modeling of snow processes has been greatly refined. Results show that both the determination of rainfall/snowfall and the snowpack process were significantly improved, contributing to a reduced summer evapotranspiration and thus an improved streamflow simulation.

Effect of macropores on soil freezing and thawing with infiltration

AbstractAn understanding of heat transport and water flow in unsaturated soils experiencing freezing and thawing is important when considering hydrological and thermal processes in cold regions. Macropores, such as cracks, roots, and animal holes, provide efficient conduits for enhanced infiltration, resulting in a unique distribution of water content. However, the effects of macropores on soil freezing and thawing with infiltration have not been well studied. A one‐directional soil‐column freezing and thawing experiment was conducted using unsaturated sandy and silt loams with different sizes and numbers of macropores. During freezing, macropores were found to retard the formation of the frozen layer, depending on their size and number. During thawing, water flowed through macropores in the frozen layer and reached the underlying unfrozen soil. However, infiltrated water sometimes refroze in a macropore. The ice started to form at near inner wall of the macropore, grew to the centre, and blocked flow through the macropore. The blockage ice in the macropore could not melt until the frozen layer disappeared. Improving a soil freezing model to consider these macropore effects is required. Copyright © 2016 John Wiley & Sons, Ltd.

Variations of total dissolved iron and its impacts during an extreme spring flooding event in the Songhua River

Total dissolved iron (TDFe) is a controlling factor in primary productivity in marine ecosystems, hence the variation of riverine TDFe output is attracting an increasing attention. Spring flooding is a key hydrological process in cold regions and the extreme spring flooding event (ESE) is thought to influence the output of TDFe considerably. In 2013, an ESE arose in the Songhua River, to reveal its impact on TDFe output and its species, water samples were collected in the river during this ESE as well as normal spring flooding period (NSP) in 2014 and 2015 to analyze the concentration and species of TDFe as well as other basic parameters. The speciation analysis was conducted by filtration and ultrafiltration methods. The results indicated the concentration of TDFe did not represent a significant difference during ESE, with an average concentration of 0.28mg/L in ESE and 0.30mg/L in NSP. The steady trend of TDFe concentration can be contributed to the limited intensity and duration of erosion by the snowmelt runoff and the construction of hydraulic works. Whereas, ESE intensified TDFe output increasing from 27.93ton/day during NSP to 48.56ton/day during ESE due to its high discharge. The species of TDFe was dominated by Fe(III) with the molecular weight lower than 10kDa. The correlation analysis indicated that content and property of DOM were the main controlling factors of TDFe species. Furthermore, the dynamics and migration of TDFe were accompanied by the nutrient and heavy metal transport which may potentially influence the water quality in the river basin.

The effects of freeze–thaw (F–T) cycles on the Gödene travertine used in historical structures in Konya (Turkey)

Travertine, a building material used around the globe for its decorative and easy-to-be-worked features from past to present, is naturally exposed to the freeze–thaw (F–T) process in cold regions in winters. As a result of this process, building stones deteriorate partially or completely. The greatest factor during the F–T process is the existence of water. The water contained in the body (in discontinuities and pores) of a building stone freezes along with the fall of temperature below 0°C. Following freezing, the volume of water contained in discontinuities and pores increases. New fractures develop on the surfaces of pores as a result of outward pressure due to this increase in volume. As the temperature drops below 0°C, the volume of water increases again, due to freezing. This process repeats itself with the decrease and increase in temperature and causes undesired deteriorations in the bodies of stones.This issue is quite important from the point of view of preservation and restoration of historical buildings, and foreseeing the effects of the F–T process on new buildings to be designed.In this study, the effects of the F–T process on Gödene travertine were investigated. Gödene travertine is a stone that has been widely used in historical buildings in Konya, a city that also covers Çatal Höyük, which is known as one of the oldest human settlements in Anatolia. The stone is still being used in buildings constructed in the region. Fresh travertine stones obtained from the quarry were exposed to freeze and thaw cycles in different numbers. Values of porosity (n), uniaxial compressive strength (σu), point load strength (Is(50)), Brazilian tensile strength (σt), Böhme abrasion loss (BA) and P wave velocity (Vp) were statistically evaluated and thus, the effects of the number of cycles on basic physical and mechanical characteristics of the Gödene stone were investigated. Furthermore, deteriorations that took place in historical structures in the region were case-studied.

Constraining Parameter Uncertainty in Simulations of Water and Heat Dynamics in Seasonally Frozen Soil Using Limited Observed Data

Water and energy processes in frozen soils are important for better understanding hydrologic processes and water resources management in cold regions. To investigate the water and energy balance in seasonally frozen soils, CoupModel combined with the generalized likelihood uncertainty estimation (GLUE) method was used. Simulation work on water and heat processes in frozen soil in northern China during the 2012/2013 winter was conducted. Ensemble simulations through the Monte Carlo sampling method were generated for uncertainty analysis. Behavioral simulations were selected based on combinations of multiple model performance index criteria with respect to simulated soil water and temperature at four depths (5 cm, 15 cm, 25 cm, and 35 cm). Posterior distributions for parameters related to soil hydraulic, radiation processes, and heat transport indicated that uncertainties in both input and model structures could influence model performance in modeling water and heat processes in seasonally frozen soils. Seasonal courses in water and energy partitioning were obvious during the winter. Within the day-cycle, soil evaporation/condensation and energy distributions were well captured and clarified as an important phenomenon in the dynamics of the energy balance system. The combination of the CoupModel simulations with the uncertainty-based calibration method provides a way of understanding the seasonal courses of hydrology and energy processes in cold regions with limited data. Additional measurements may be used to further reduce the uncertainty of regulating factors during the different stages of freezing–thawing.

Ground-penetrating radar-derived measurements of active-layer thickness on the landscape scale with sparse calibration at Toolik and Happy Valley, Alaska

Active-layer thickness (ALT) is an important parameter for studying surface energy balance, ecosystems, and hydrologic processes in cold regions. We measured ALT along 10 routes with lengths ranging from 0.7 to 6.9 km located on the Alaska North Slope near Toolik Lake and the Happy Valley airstrip (between 68.475° and 69.150°N, and [Formula: see text] and [Formula: see text]). Using a ground-penetrating radar (GPR) system in a common-offset configuration, we measured the two-way traveltimes from the surface to the bottom of the active layer at the end of summer, when the thaw depth was greatest. We used 500 and 800 MHz antennas; the 500 MHz antenna provided suitable vertical resolution, while producing more unambiguous active-layer reflections in the presence of nonideal antenna coupling and active layer inhomogeneity. We derived ALT measurements and their uncertainties from GPR two-way traveltimes, with mechanical probing for velocity calibration. Using an empirical relationship between the wave velocity and soil volumetric water content (VWC), we found that the velocities were consistent with soil VWCs ranging from 0.46 to 0.63. In 31% of traces, the permafrost table horizon was identifiable, resulting in ALT measurements with uncertainties of generally less than 25%. The average ALT was 48.1 cm, with a standard deviation of 16.1 cm. We found distinct patterns of ALT spatial variability at different sites and different length scales. At some sites, the ALT at one point was effectively uncorrelated with ALT at other points separated by lag distances as small as tens of meters; for other sites, there was correlation at lag distances up to approximately 400 m. The ALT statistics were similar to nearby long-term in situ ALT measurements from the Circumpolar Active Layer Monitoring Network, through which yearly ALT measurements have been made since 1990.

Gravel transport by ice in a subarctic river from accurate laser scanning

For decades the importance of ice and the effects of cold-region processes on river channel morphology have been discussed, with a general consensus as to their importance emerging only recently. River ice cover, anchor ice, frazil ice, and ice jams may not only scour the channel bed and banks but also pick up, transport, and deposit fine sediments and gravels during winter, especially during the spring ice breakup period. However, knowledge of the interactions between coarse sediment transport and ice processes remains insufficient, particularly in rockier river reaches, with a lack of accurate and sufficiently extensive data hindering their quantification.The aim of this study was to quantify and analyse the impact of river ice on gravel transport in a subarctic river during one winter via the acquisition of laser scanning data for the river channel and ice surface. Terrestrial and mobile laser scanning were performed in 2012–2013 on the Tana River in northern Finland. Both of these techniques are considered accurate and applicable for detecting elevation and volumetric changes in river bed, defining gravel clast sizes, and detecting the movement of individual clasts. More importantly, ice surface, thickness, and decay during spring were also captured via laser scanning.In the winter of 2012–2013, a period characterised by an absence of ice jams and mid-winter ice-decay periods, with spring ice breakup discharges close to average yearly conditions, ice had the most significant role, greater than that of flowing water, in erosion and transport of coarse sediment from the channel bed and gently sloping banks. Changes in river bed elevation and volume were recorded throughout the study site, and erosion predominated. In addition to broader scale erosion, the movement of single clasts up to 2m in size occurred. However, the observed overall channel change patterns did not coincide with the areas of fastest ice decay. The obtained results could also be applied to the enhancement of riverine planning in subarctic environments.

Integrated simulation of snow and glacier melt in water and energy balance‐based, distributed hydrological modeling framework at Hunza River Basin of Pakistan Karakoram region

AbstractEnergy budget‐based distributed modeling of snow and glacier melt runoff is essential in a hydrologic model to accurately describe hydrologic processes in cold regions and high‐altitude catchments. We developed herein an integrated modeling system with an energy budget‐based multilayer scheme for clean glaciers, a single‐layer scheme for debris‐covered glaciers, and multilayer scheme for seasonal snow over glacier, soil, and forest within a distributed biosphere hydrological modeling framework. Model capability is demonstrated for Hunza River Basin (13,733 km2) in the Karakoram region of Pakistan on a 500 m grid for 3 hydrologic years (2002–2004). Discharge simulation results show good agreement with observations (Nash‐Sutcliffe efficiency = 0.93). Flow composition analysis reveals that the runoff regime is strongly controlled by the snow and glacier melt runoff (50% snowmelt and 33% glacier melt). Pixel‐by‐pixel evaluation of the simulated spatial distribution of snow‐covered area against Moderate Resolution Imaging Spectroradiometer‐derived 8 day maximum snow cover extent data indicates that the areal extent of snow cover is reproduced well, with average accuracy 84% and average absolute bias 7%. The 3 year mean value of net mass balance (NMB) was estimated at +0.04 myr−1. It is interesting that individual glaciers show similar characteristics of NMB over 3 years, suggesting that both topography and glacier hypsometry play key roles in glacier mass balance. This study provides a basis for potential application of such an integrated model to the entire Hindu‐Kush‐Karakoram‐Himalaya region toward simulating snow and glacier hydrologic processes within a water and energy balance‐based, distributed hydrological modeling framework.

Enhanced physicochemical-biological sewage treatment process in cold regions.

Biological treatment processes give relatively poor pollutant removal efficiencies in cold regions because microbial activity is inhibited at low temperatures. We developed an enhanced physicochemical-biological wastewater treatment process that involves micro-membrane filtration, anaerobic biofilter, and aerobic biofilter to improve the pollutant removal efficiencies that can be achieved under cold conditions. Full-scale experiments using the process were carried out in the northeast of China, at outdoor temperatures of around -30 °C. The average removal efficiencies achieved for chemical oxygen demand, total phosphorus, ammonia nitrogen, and suspended solids were 89.8, 92.9, 94.3, and 95.8%, respectively, using a polyaluminium chloride dosage of 50 mg L⁻¹. We concluded that the process is effective to treat sewage in cold regions.

Changes in microstructure and physical properties of rocks caused by artificial freeze–thaw action

The freeze–thaw action of water inside rocks’ pores, cracks, and joints plays a dominant role in physical rock weathering where the temperature fluctuates across 0 °C. To understand this process in cold regions, artificial weathering was simulated in the laboratory. Diorite, basalt, and tuff were used as specimen, and were maintained in saturated conditions to accelerate the effect of freeze–thaw weathering. X-ray computed tomography (CT) images and scanning electron microscope (SEM) photographs were obtained to investigate changes in the inner microstructure along with repeated freeze–thaw cycles. In addition, physical properties were measured during the cycles. Particle detachment, crack initiation, and propagation and porosity increase due to volumetric expansion of water inside rocks were detected using X-ray CT and SEM. The decrease in P-wave and S-wave velocities, and in weight loss were also observed. Tuff specimens, which have low tensile strength, high porosity, and loose structure, were deteriorated rapidly, while diorite specimens on the other hand, which have high tensile strength, low porosity, and compact structure, showed relatively slight changes. The frost susceptibility of basalt specimens was found to be intermediate between tuff and diorite. The developed process for cycling and freezing set up can be useful to test the various stones in terms of frost susceptibility.