Meltwater Percolation Research Articles

Direct measurements of firn-density evolution from 2016 to 2022 at Wolverine Glacier, Alaska

Abstract Knowledge of snow and firn-density change is needed to use elevation-change measurements to estimate glacier mass change. Additionally, firn-density evolution on glaciers is closely connected to meltwater percolation, refreezing and runoff, which are key processes for glacier mass balance and hydrology. Since 2016, the U.S. Geological Survey Benchmark Glacier Project has recovered firn cores from a site on Wolverine Glacier in Alaska's Kenai Mountains. We use annual horizons in repeat cores to track firn densification and meltwater retention over seasonal and interannual timescales, and we use density measurements to quantify how the firn air content (FAC) changes through time. The results suggest the firn is densifying due primarily to compaction rather than refreezing. Liquid-water retention in the firn is transient, likely due to gravity-fed drainage and irreducible-water-content decreases that accompany decreasing porosity. We show that the uncertainty (±60 kg m−3) in the commonly used volume-to-mass conversion factor of 850 kg m−3 is an underestimation when glacier-wide FAC variability exceeds 12% of the glacier-averaged height change. Our results demonstrate how direct measurements of firn properties on mountain glaciers can be used to better quantify the uncertainty in geodetic volume-to-mass conversions.

How water, temperature, and seismicity control the preconditioning of massive rock slope failure (Hochvogel)

Abstract. The anticipation of massive rock slope failures is a key mitigation strategy in a changing climate and environment requiring a precise understanding of pre-failure process dynamics. Here we exploit >4 years of multi-method high-resolution monitoring data from a large rock slope instability close to failure. To quantify and understand the effect of possible drivers (water from rain and snowmelt, internal rock fracturing, and earthquakes), we correlate slope displacements with environmental data, local seismic recordings, and earthquake catalogues. During the snowmelt phase, displacements are controlled by meltwater infiltration with high correlation and a time lag of 4–9 d. During the snow-free summer, rainfall induces accelerations with a time lag of 1–16 h for up to several days without a minimum activation rain sum threshold. Rock fracturing, linked to temperature and freeze–thaw cycles, is predominantly near the surface and unrelated to displacement rates. A classic Newmark analysis of recent and historic earthquakes indicates a low potential for immediate triggering of a major failure at the case site, unless it is already very close to failure. Seismic topographic amplification of the peak ground velocity (PGV) at the summit ranges from a factor of 2–11 and is spatially heterogeneous, indicating a high criticality of the slope. The presented in-depth monitoring data analysis enables a comprehensive rockfall driver evaluation and indicates where future climatic changes, e.g. in precipitation intensity and frequency, may alter the preconditioning of major rock slope failures.

Mapping the vertical heterogeneity of Greenland's firn from 2011–2019 using airborne radar and laser altimetry

Abstract. The firn layer on the Greenland Ice Sheet (GrIS) plays a crucial role in buffering surface meltwater runoff, which is constrained by the available firn pore space and impermeable ice layers that limit deeper meltwater percolation. Understanding these firn properties is essential for predicting current and future meltwater runoff and its contribution to global sea-level rise. While very-high-frequency (VHF) radars have been extensively used for surveying the GrIS, their lower bandwidth restricts direct firn stratigraphy extraction. In this study, we use concurrent VHF airborne radar and laser altimetry data collected as part of Operation IceBridge over the 2011–2019 period to investigate our hypothesis that vertical heterogeneities in firn (i.e. ice layers) cause vertical offsets in the radar surface reflection (dz). Our results, corroborated by modelling and firn core analyses, show that a dz larger than 1 m is strongly related to the vertical heterogeneity of a firn profile and effectively delineates between vertically homogeneous and vertically heterogeneous firn profiles over a depth range of ∼ 4 m. Temporal variations in dz align with climatic events and reveal an expansion of heterogeneous firn between 2011–2013 covering an area of ∼ 350 815 km2, followed by firn replenishment over the years 2014–2019 spanning an area of ∼ 667 725 km2. Our approach reveals the firn evolution of key regions on the Greenland Ice Sheet, providing valuable insights for detecting potential alterations in meltwater runoff patterns.

A cold laboratory hyperspectral imaging system to map grain size and ice layer distributions in firn cores

Abstract. The Greenland and Antarctic ice sheets are covered in a layer of porous firn. Knowledge of firn structure improves our understanding of ice sheet mass balance, supra- and englacial hydrology, and ice core paleoclimate records. While macroscale firn properties, such as firn density, are relatively easy to measure in the field or lab, more intensive measurements of microstructural properties are necessary to reduce uncertainty in remote sensing observations of mass balance, model meltwater infiltration, and constrain ice age – gas age differences in ice cores. Additionally, as the duration and extent of surface melting increases, refreezing meltwater will greatly alter firn structure. Field observations of firn grain size and ice layer stratigraphy are required to test and validate physical models that simulate the ice-sheet-wide evolution of the firn layer. However, visually measuring grain size and ice layer distributions is tedious, is time-consuming, and can be subjective depending on the method. Here we demonstrate a method to systematically map firn core grain size and ice layer stratigraphy using a near-infrared hyperspectral imager (NIR-HSI; 900–1700 nm). We scanned 14 firn cores spanning ∼ 1000 km across western Greenland’s percolation zone with the NIR-HSI mounted on a linear translation stage in a cold laboratory. We leverage the relationship between effective grain size, a measure of NIR light absorption by firn grains, and NIR reflectance to produce high-resolution (0.4 mm) maps of effective grain size and ice layer stratigraphy. We show the NIR-HSI reproduces visually identified ice layer stratigraphy and infiltration ice content across all cores. Effective grain sizes change synchronously with traditionally measured grain radii with depth, although effective grains in each core are 1.5× larger on average, which is largely related to the differences in measurement techniques. To demonstrate the utility of the firn stratigraphic maps produced by the NIR-HSI, we track the 2012 melt event across the transect and assess its impact on deep firn structure by quantifying changes to infiltration ice content and grain size. These results indicate that NIR-HSI firn core analysis is a robust technique that can document deep and long-lasting changes to the firn column from meltwater percolation while quickly and accurately providing detailed firn stratigraphy datasets necessary for firn research applications.

The application and modification of WRF-Hydro/Glacier to a cold-based Antarctic glacier

Abstract. The McMurdo Dry Valleys (MDV) are home to a unique microbial ecosystem that is dependent on the availability of freshwater. This is a polar desert and freshwater originates almost entirely from surface and near-surface melt of the cold-based glaciers. Understanding the future evolution of these environments requires the simulation of the full chain of physical processes from net radiative forcing, surface energy balance, melt, runoff and transport of meltwater in stream channels from the glaciers to the terminal lakes where the microbial community resides. To establish a new framework to do this, we present the first application of WRF-Hydro/Glacier in the MDV, which as a fully distributed hydrological model has the capability to resolve the streams from the glaciers to the bare land that surround them. Given that meltwater generation in the MDV is almost entirely dependent on small changes in the energy balance of the glaciers, the aim of this study is to optimize the multi-layer snowpack scheme that is embedded in WRF-Hydro/Glacier to ensure that the feedbacks between albedo, snowfall and melt are fully resolved. To achieve this, WRF-Hydro/Glacier is implemented at a point scale using automatic weather station data on Commonwealth Glacier to physically model the onset, duration and end of melt over a 7-month period (1 August 2021 to 28 February 2022). To resolve the limited energetics controlling melt, it was necessary to (1) limit the percolation of meltwater through the ice layers in the multi-layer snowpack scheme and (2) optimize the parameters controlling the albedo of both snow and ice over the melt season based on observed spectral signatures of albedo. These modifications enabled the variability of broadband albedo over the melt season to be accurately simulated and ensured that modelled surface and near-surface temperatures, surface height change and runoff were fully resolved. By establishing a new framework that couples a detailed snowpack model to a fully distributed hydrological model, this work provides a stepping stone to model the spatial and temporal variability of melt and streamflow in the future, which will enable some of the unknown questions about the hydrological connectivity of the MDV to be answered.

Monitoring hydrological changes with satellite data: Spring thaw's effect on soil moisture and groundwater in seasonal Freezing-Thawing zones

Due to the dense vegetation coverage in mountainous areas and difficulties in field monitoring in mid-high-latitude freezing-thawing zones, satellite and remote sensing technologies can be used to monitor long-term global hydrological changes. This study uses multiple satellite data of snow water equivalent (SWE), soil moisture (SM), and groundwater storage anomaly (GWSA) to analyse their spatial–temporal changes in the seasonal freezing-thawing Tumen River Basin. Water transformation during the thawing period is determined through groundwater level fluctuations and correlation analysis methods; thus, the recharge effects of snowmelt and seasonally frozen soil thawing on soil moisture and groundwater are identified. The snow begins to melt in February, and snow water equivalent decreases gradually from upstream to downstream after reaching the peak. Furthermore, seasonally frozen soil melts gradually from March. Soil moisture downstream is higher than that in the upper and middle reaches. Additionally, the GWSA increases gradually in the thawing period, which is high in mountainous areas and low in plains. The Spearman correlation coefficient among ΔSWE, ΔSM, and ΔGWSA is more significant than 0.5 and dependent on altitude. With snow and frozen soil melting, vertical infiltration and meltwater percolation increase soil moisture and GWSA. This phenomenon is predominant at low altitudes (<600 m) and pore phreatic groundwater in the middle and lower reaches. However, snowmelt is the dominant factor in groundwater recharge at high-altitude upper reaches. For seasonal freezing-thawing zones, a close water transformation was observed between snow melting, soil water thawing, and groundwater in spring. Thus, the environmental and ecological impact of spring snowmelt and floods in freezing-thawing zones can be monitored using multi-source data.

Modern air, englacial and permafrost temperatures at high altitude on Mt Ortles (3905 m a.s.l.), in the eastern European Alps

Abstract. The climatic response of mountain permafrost and glaciers located in high-elevation mountain areas has major implications for the stability of mountain slopes and related geomorphological hazards, water storage and supply, and preservation of palaeoclimatic archives. Despite a good knowledge of physical processes that govern the climatic response of mountain permafrost and glaciers, there is a lack of observational datasets from summit areas. This represents a crucial gap in knowledge and a serious limit for model-based projections of future behaviour of permafrost and glaciers. A new observational dataset is available for the summit area of Mt Ortles, which is the highest summit of South Tyrol, Italy. This paper presents a series of air, englacial, soil surface and rock wall temperatures collected between 2010 and 2016. Details are provided regarding instrument types and characteristics, field methods, and data quality control and assessment. The obtained data series are available through an open data repository (https://doi.org/10.5281/zenodo.8330289, Carturan et al., 2023). In the observed period, the mean annual air temperature at 3830 m a.s.l. was between −7.8 and −8.6 ∘C. The most shallow layers of snow and firn (down to a depth of about 10 m) froze during winter. However, melt water percolation restored isothermal conditions during the ablation season, and the entire firn layer was found at the melting pressure point. Glacier ice is cold, but only from about 30 m depth. Englacial temperature decreases with depth, reaching a minimum of almost −3 ∘C close to the bedrock, at 75 m depth. A small glacier located at 3470 m a.s.l., close to the summit of Mt Ortles, was also found in cold conditions down to a depth of 9.5 m. The mean annual ground surface temperature was negative for all but one monitored sites, indicating cold ground conditions and the existence of permafrost in nearly all debris-mantled slopes of the summit. Similarly, the mean annual rock wall temperature was negative at most monitored sites, except the lowest one at 3030 m a.s.l. This suggests that the rock faces of the summit are affected by permafrost at all exposures.

Integration of Sentinel–1 and Landsat–8 data for identifying deformation risk of hydraulic engineering in seasonally frozen salinization regions

Hydraulic engineering has positive effects on agricultural production and soil desalination. However, freeze–thaw cycles can cause instability in canal slopes, emphasizing the need for effective deformation monitoring. Using Sentinel-1B data acquired from November 2017 to August 2019 and applying the small baseline subset (SBAS) InSAR method, two unstable canal segments with vertical deformation rates exceeding –40 mm/year were identified in Songyuan, Northeast China. The canal slopes uplifted in winter and experienced accelerated subsidence in spring, with the shady slope experiencing a deformation rate approximately 10 mm/year faster than the sunny slope. Further investigation, integrating Landsat–8 data collected from November 2017 to August 2019, revealed that the shady slope received less thermal radiation than the sunny slope, with an average annual temperature difference on the land surface of approximately 1 ℃. In winter, the saline–alkali soil particles formed a thick unfrozen water film, which, driven by temperature and soil–water gradients, migrated towards the freezing front, froze, and caused soil frost heave. The shady slope experienced longer freezing periods and deeper frozen depths, up to 240 cm, causing more unfrozen water migration and greater uplift deformation than the sunny slope. In spring, as temperature rose, the soil thawed from the surface downward. Notably, the shady slope thawed at a slower rate than the sunny slope, forming a frozen lens near a depth of 100 cm. Frozen lens acted as a barrier for meltwater infiltration, impeding the dissipation of pore water pressure, making the shady slope prone to surface landslides. The integration of multiple data sources is beneficial for identifying deformation risk in hydraulic engineering.

Glide-snow avalanche characteristics at different timescales extracted from time-lapse photography

Abstract Glide-snow avalanches release due to a loss of friction at the snow–ground interface, which can result in large avalanches that endanger infrastructure in alpine regions. The underlying processes are still relatively poorly understood, in part due to the limited data available on glide processes. Here, we introduce a pixel-based algorithm to detect glide cracks in time-lapse photographs under changing illumination and shadow conditions. The algorithm was applied to 14 years of time-lapse images at Dorfberg (Davos, Switzerland). We analysed 947 glide-snow events at a high-spatial (0.5 m) and temporal (2–15 min) resolution. Avalanche activity and glide-crack opening dynamics were investigated across timescales ranging from seasonally to hourly. Events were separated into surface (meltwater percolation) and interface events (no meltwater percolation). The results show that glide activity is highly variable between and within seasons. Most avalanches released without crack formation or within 24 h after crack opening, and release was favoured in the afternoon hours. Glide rates often showed a stick–slip behaviour. The acceleration of glide rates and non-constant increases in glide crack aspect ratio were indicators for avalanche release. This comprehensive dataset provides the basis for further investigations into glide-snow avalanche drivers.

Spatial characterization of near-surface structure and meltwater runoff conditions across the Devon Ice Cap from dual-frequency radar reflectivity

Abstract. Melting and refreezing processes in the firn of the Devon Ice Cap control meltwater infiltration and runoff across the ice cap, but their full spatial extent and effect on near-surface structure is difficult to measure with surface-based traverses or existing satellite remote sensing. Here, we derive the coherent component of the near-surface return from airborne ice-penetrating radar surveys over the Devon Ice Cap, Canadian Arctic, to characterize firn containing centimeter- to meter-thick ice layers (i.e., ice slabs) formed from refrozen meltwater in firn. We assess the use of dual-frequency airborne ice-penetrating radar to characterize the spatial and vertical near-surface structure of the Devon Ice Cap by leveraging differences in range resolution of the radar systems. Comparison with reflectivities using a thin layer reflectivity model, informed by surface-based radar and firn core measurements, indicates that the coherent component is sensitive to the near-surface firn structure composed of quasi-specular ice and firn layers, limited by the bandwidth-constrained radar range resolution. Our results suggest that average ice slab thickness throughout the Devon Ice Cap percolation zone ranges from 4.2 to 5.6 m. This implies conditions that can enable lateral meltwater runoff and potentially contribute to the total surface runoff routed through supraglacial rivers down glacier. Together with the incoherent component of the surface return previously studied, our dual-frequency approach provides an alternative method for characterizing bulk firn properties, particularly where high-resolution radar data are not available.

Firn Core Evidence of Two‐Way Feedback Mechanisms Between Meltwater Infiltration and Firn Microstructure From the Western Percolation Zone of the Greenland Ice Sheet

AbstractThe relationship between firn microstructure and water movement is complex: firn microstructure controls the routing of meltwater through the firn while continuously being altered by liquid water flow processes. Importantly, microstructural transitions within the firn column can stall vertical meltwater percolation, which creates heterogeneities in liquid water content resulting in different rates of firn metamorphism. Physics‐based firn models aim to describe these processes to accurately predict ice layer or firn aquifer formation, but require detailed observations of firn structure to validate and inform percolation schemes. Here, we present grain size measurements and ice layer stratigraphy from seven firn cores collected in western Greenland's percolation zone during the 2016 Greenland Traverse for Accumulation and Climate Studies (GreenTrACS). Grain size transitions within the cores are negatively correlated with all temperature proxies for meltwater supply. Additionally, the number of grain size transitions are strongly anticorrelated with the number of ice layers within each core, despite these transitions, particularly fine‐over‐coarse transitions, promoting meltwater ponding and potential ice layer formation. To investigate if these negative correlations can be understood with firn model physics, we simulate water movement along stratigraphic transitions using the SNOWPACK model. We find that grain size transitions diminish from rapid grain growth in wet firn where ice layers can form, suggesting these microstructural transitions are unlikely to survive repeated meltwater infiltration. Incorporating these microstructure—meltwater feedbacks in firn models could improve their ability to model processes such as ice slab formation or firn aquifer recharge that require accurate predictions of meltwater infiltration depth.

Levels and distribution profiles of Per- and Polyfluoroalkyl Substances (PFAS) in a high Arctic Svalbard ice core

Per- and polyfluoroalkyl substances (PFAS) are a group of persistent organic contaminants of which some are toxic and bioaccumulative. Several PFAS can be formed from the atmospheric degradation of precursors such as fluorotelomer alcohols (FTOHs) as well as hydrochlorofluorocarbons (HFCs) and other ozone-depleting chlorofluorocarbon (CFC) replacement compounds. Svalbard ice cores have been shown to provide a valuable record of long-range atmospheric transport of contaminants to the Arctic. This study uses a 12.3 m ice core from the remote Lomonosovfonna ice cap on Svalbard to understand the atmospheric deposition of PFAS in the Arctic. A total of 45 PFAS were targeted, of which 26 were detected, using supercritical fluid chromatography (SFC) tandem mass spectrometry (MS/MS) and ultra-performance liquid chromatography (UPLC) MS/MS. C2 to C11 perfluoroalkyl carboxylic acids (PFCAs) were detected continuously in the ice core and their fluxes ranged from 2.5 to 8200 ng m−2 yr−1 (9.51–16,500 pg L−1). Trifluoroacetic acid (TFA) represented 71 % of the total mass of C2 – C11 PFCAs in the ice core and had increasing temporal trends in deposition. The distribution profile of PFCAs suggested that FTOHs were likely the atmospheric precursor to C8 – C11 PFCAs, whereas C2 – C6 PFCAs had alternative sources, such as HFCs and other CFC replacement compounds. Perfluorooctanesulfonic acid (PFOS) was also widely detected in 82 % of ice core subsections, and its isomer profile (81 % linear) indicated an electrochemical fluorination manufacturing source. Comparisons of PFAS concentrations with a marine aerosol proxy showed that marine aerosols were insignificant for the deposition of PFAS on Lomonosovfonna. Comparisons with a melt proxy showed that TFA and PFOS were mobile during meltwater percolation. This indicates that seasonal snowmelt and runoff from post-industrial accumulation on glaciers could be a significant seasonal source of PFAS to ecosystems in Arctic fjords.

Interactions between shrub encroachment and water infiltration on the hillslope of a typical steppe

AbstractShrub encroachment has been a worldwide phenomenon, but the interactions between the surface hydrological cycle and the shrub‐encroached landscape at different hillslope positions remain poorly investigated. The present study was undertaken to explore the interactions between the water infiltration patterns affected by shrub encroachment. At four slope positions of a hillslope encroached by Caragana microphylla Lam, soil water content and temperature were continuously measured at 10‐min intervals at four or five depths in both the presence (shrub patches) and absence (grass matrix) of shrubs. Then, the infiltration of rainfall and meltwater was estimated based on analyses of the above data. The rainfall infiltration ratios (IRs) of grass matrix were as high as 0.78 ± 0.08, except at the lower site, which had a value of only 0.47. The IRs of shrub patches increased from 0.38 at the top site to 0.77 at the lower site. In snowy years, snow and ice were blown upward by wind and accumulated in the shrub patches and their leeward areas, which resulted in higher water input for the shrub patches than for the grass matrix. Shrub encroachment changed the microtopography and soil properties of the hillslope and further affected the surface hydrological processes. The feedbacks between shrub encroachment and water infiltration varied among sites of the hillslope, which may in turn have affected the development of shrub patches.

Shallow Fracture Buffers High Elevation Runoff in Northwest Greenland

AbstractThe expansion of refrozen ice slabs in Greenland's firn may enhance meltwater runoff and increase surface mass loss. However, the impermeability of ice slabs and the pathways for meltwater export from these regions remain poorly characterized. Here, we present ice‐penetrating radar observations of extensive meltwater infiltration and refreezing beneath ice slabs in Northwest Greenland. We show that these buried ice complexes form where supraglacial streams or lakes drain through surface crevasses into relict firn beneath the ice slabs. This suggests that the firn can continue to buffer mass loss from surface meltwater runoff and limit meltwater delivery to the ice sheet bed even after ice slabs have formed. Therefore, a significant time lag may exist between the initial formation of ice slabs and the onset of complete surface runoff and seasonal meltwater drainage to the subglacial system in interior regions of the ice sheet.

Prefacing unexplored archives from Central Andean surface-to-bedrock ice cores through a multifaceted investigation of regional firn and ice core glaciochemistry

AbstractShallow firn cores, in addition to a near-basal ice core, were recovered in 2018 from the Quelccaya ice cap (5470 m a.s.l) in the Cordillera Vilcanota, Peru, and in 2017 from the Nevado Illimani glacier (6350 m a.s.l) in the Cordillera Real, Bolivia. The two sites are ~450 km apart. Despite meltwater percolation resulting from warming, particle-based trace element records (e.g. Fe, Mg, K) in the Quelccaya and Illimani shallow cores retain well-preserved signals. The firn core chronologies, established independently by annual layer counting, show a convincing overlap indicating the two records contain comparable signals and therefore capture similar regional scale climatology. Trace element records at a ~1–4 cm resolution provide past records of anthropogenic emissions, dust sources, volcanic emissions, evaporite salts and marine-sourced air masses. Using novel ultra-high-resolution (120 μm) laser technology, we identify annual layer thicknesses ranging from 0.3 to 0.8 cm in a section of 2000-year-old radiocarbon-dated near-basal ice which compared to the previous annual layer estimates suggests that Quelccaya ice cores drilled to bedrock may be older than previously suggested by depth-age models. With the information collected from this study in combination with past studies, we emphasize the importance of collecting new surface-to-bedrock ice cores from at least the Quelccaya ice cap, in particular, due to its projected disappearance as soon as the 2050s.

Dust dominates the summer melting of glacier ablation zones on the northeastern Tibetan Plateau

Dust and black carbon (BC) can darken snow and ice surface and play pivotal roles in glacier mass loss. Thus, a quantitative assessment of their contributions to glacier summer melting is critical. During the summer of 2018, surface snow and ice were sampled, and the albedo and mass balance were continuously measured in the ablation zone of Laohugou Glacier No. 12 in the western Qilian Mountains. The physical properties of dust and BC were measured in the laboratory, and their impacts on glacier surface albedo reduction and melting were simulated. The results indicate that the ice surface in the ablation zone was enriched with substantial amounts of particles, and the average particle concentrations of these samples were hundreds of times higher than those of fresh snow. The BC mass absorption cross-sections (MACs) ranged from 3.1 m2 g−1 at 550 nm for dirty ice to 4.6 m2 g−1 for fresh snow, largely owing to meltwater percolation and particle collapse. The spectral variations in dust MACs were significantly different in the visible light bands and near-infrared bands from those in the other areas. Moreover, the two-layer surface energy and mass balance model with the new albedo parameterization formula was validated and agreed well with the experimental measurements of spectral albedo, broadband albedo, and mass balance. BC and dust combined resulted in 26.7 % and 54.4 % of the total mass loss on the cleaner and dirtier (particle enriched) surfaces in the ablation zone, respectively, compared to particle-free surfaces, and although both impurities played vital roles, dust was the more prominent factor in accelerating glacier melting on the northeastern Tibetan Plateau. This study emphasizes the importance of dust in cryosphere changes where Tibetan glaciers are strongly affected by Asian dust deposition.

Roots and other residues from leys with or without red clover: Quality and effects on N2O emission factors in a partly frozen soil following autumn ploughing

Revised IPCC guidelines assume that a constant share of N in decomposing crop residues is directly emitted as N2O (emission factor: EFN2O), and calculate the amount of nitrogen (N) in non-removable residues of temporary grasslands proportionally to the average annual herbage yield. However, EFN2O depends on the intrinsic quality of the residues and their interactions with environmental conditions. Only a few field studies on N2O emissions from grassland renewal are available, and none have simultaneously quantified the N amount and quality of non-removable residues (roots and stubble).To gain insight into the effect of non-removable residue quality on EFN2O, we studied the amount and quality of roots and stubble and their effect on EFN2O following the ploughing of three-year-old swards.The measured amount of N in non-removable residues was approximately 20, 25, and 31 kg N per 1 Mg average annual dry matter yield in grass, red clover-grass, and red clover, and 70–83% of it was below ground. However, the EFN2O of non-removable residues measured over 252 days was lower (0.24%, SE = 14% for grass and red clover-grass) than the IPCC default value (0.6%, CV: 50%) for wet regions, although within the uncertainty margin, and was significantly lower than the EFN2O of incorporated herbage, which was related to differences in EFCO2. We advocate for more specific studies that separate the effects of belowground and aboveground residues (AGR), considering the possibility of simplifying the accounting of N2O emissions from belowground residues while improving that of non-removable AGR from temporary grasslands and other green crops.We observed the accumulation of N2O in the frozen soil under snow, which was released during diurnal percolation of meltwater. N2O emissions from frozen soil accounted for 30% or more of the total emissions.

Potential recharge sources and origin of solutes in groundwater in the central Qinghai–Tibet Plateau using hydrochemistry and isotopic data

Study regionThe Nagqu River Basin in the central part of the Qinghai–Tibet Plateau, China. Study focusGroundwater flow systems in permafrost regions are characterized by a unique circulation. Understanding the groundwater circulation mechanism is essential for studying the hydrological effect of permafrost degradation in alpine regions. In this study, potential water sources and origin of solutes in the groundwater of a typical alpine region were studied using hydrochemistry, 18O, 2H and 3H analyses. This research is essential for studying the hydrological effects of permafrost degradation. New hydrological insights for the regionUpon moving from the mountainous region to the plains, the renewal rate decreases from 10.5% to 0.6% and the mean transit time increases from 5–30 to 20–32 years. Groundwater is mainly recharged by ground ice meltwater and infiltration of glacier meltwater, snow meltwater and precipitation in mountainous areas, and by deep groundwater near faults. Discharge to surface water, evaporation, and overflow into springs are the main drainage modes. Water/ice–rock/soil interactions, such as dissolution of aluminosilicate minerals, are the main factors affecting the hydrochemical evolution. This study provides a data-driven approach for understanding groundwater recharge and evolution in similar systems.

Влияние полезащитных лесных полос ажурной конструкции на характер снегораспределения

Snow in agricultural landscapes is a reserve of moisture, as well as protection of soils from deep freezing and protection of wintering plants from destructive effects of low temperatures. In the zone of insufficient moisture to which the study area belongs, the snow cover is one of the main sources of water resources, as the moisture reserves in the soil, necessary for growth and development of agricultural crops, are formed to a greater extent by meltwater infiltration. Forest belts have a direct impact on snow accumulation and character of snow deposition. The purpose of this work is to study the effect of openwork forest shelter belts on snow accumulation and snow distribution in agricultural landscapes in Ilovlinsky District of Volgograd Oblast. The snow cover was studied in the top-down direction along the slope (perpendicular to the forest belt (FB)) by laying snow measuring profiles, along which the height of the snow cover was measured with a snow measuring rod M-46 in triplicate (in the field part every 4 m, in the forest belt and plume zones every meter). According to the results of the studies, the main zones of snow accumulation were identified: the first zone (I) – the forest belt, the second (II) – the zone of snow plume formation, the third (III) – the zone of FB influence (up to 30H, where H is the plantation height), the fourth (IV) – the entire field part of the agricultural landscape. It was found that the maximum values of snow cover were recorded in the forest belt, as well as in the plume zone (1-2H). The average snow height in the field part at the level is 7 cm, which is 14 cm lower than the snow height in the FB influence zone (1-30H). In some parts of the field, the snow cover height is reduced to 5 cm. The results obtained can be used to solve various agronomic as well as anti-erosion problems.

Shallow firn cores 1989–2019 in southwest Greenland's percolation zone reveal decreasing density and ice layer thickness after 2012

AbstractRefreezing of meltwater in firn is a major component of Greenland ice-sheet's mass budget, but in situ observations are rare. Here, we compare the firn density and total ice layer thickness in the upper 15 m of 19 new and 27 previously published firn cores drilled at 15 locations in southwest Greenland (1850–2360 m a.s.l.) between 1989 and 2019. At all sites, ice layer thickness covaries with density over time and space. At the two sites with the earliest observations (1989 and 1998), bulk density increased by 15–18%, in the top 15 m over 28 and 21 years, respectively. However, following the extreme melt in 2012, elevation-detrended density using 30 cores from all sites decreased by 15 kg m−3 a−1 in the top 3.75 m between 2013 and 2019. In contrast, the lowest elevation site's density shows no trend. Thus, temporary build-up in firn pore space and meltwater infiltration capacity is possible despite the long-term increase in Greenland ice-sheet melting.