Drainage Events Research Articles

Patterns and mechanisms of repeat drainages of glacier-dammed Dań Zhùr (Donjek) Lake, Yukon

Dań Zhùr (Donjek) Glacier is a surge-type glacier that undergoes cyclical periods of rapid advance over ∼1–2 years, followed by retreat for ∼10 years. Since the 1990s, the advances have caused the formation of ice-dammed Dań Zhùr Lake, which has filled and drained in summers following a surge event. Here, we report how these drainages initially occur through a subglacial channel under the glacier terminus, which then typically closes the following winter, enabling another lake to form and drain the next summer. However, our remote sensing and field observations indicate that after several drainage events, a subaerial ice canyon is formed through the glacier terminus, which prevents another lake from forming until after the glacier surges again. We predict that the next surge of Dań Zhùr Glacier will occur around the mid-2020s, causing the formation of a larger Dań Zhùr Lake during the following quiescent phase because, despite periodic advances, a long-term trend of glacier recession is exposing a larger basin for the lake to form in. However, each subsequent surge is causing less terminus advance than the previous one, until ultimately the surges will be insufficient to dam Dań Zhùr Chù’ (Donjek River), and lakes will cease to form.

Integrating local environmental observations and remote sensing to better understand the life cycle of a thermokarst lake in Arctic Alaska

ABSTRACT On 29 June 2022, local observers reported the drainage of a 0.5 ha lake near Qikiqtaġruk (Kotzebue), Alaska, that prompted this collaborative study on the life cycle of a thermokarst lake in the Arctic. Prior to its drainage, the lake expanded from 0.13 ha in 1951 to 0.54 ha in 2021 at lateral rates that ranged from 0.25 to 0.35 m/year. During the drainage event, we estimate that 18,500 m3 of water drained from the lake into Kotzebue Sound, forming a 125-m-long thermo-erosional gully that incised 2 to 14 m in ice-rich permafrost. Between 29 June and 18 August 2022, the drainage gully expanded from 1 m to >10 m wide, mobilizing ~8,500 m3 of material through erosion and thaw. By reconstructing a pre-lake disturbance terrain model, we show that thaw subsidence occurs rapidly (0.78 m/year) upon transition from tundra to lake but that over a seventy-year period it slows to 0.12 m/year. The combination of multiple remote sensing tools and local environmental observations provided a rich data set that allowed us to assess rates of lake expansion relative to rates of sub-lake permafrost thaw subsidence as well as hypothesizing about the potential role of beavers in arctic lake drainage.

Bedfast and floating-ice dynamics of thermokarst lakes using a temporal deep-learning mapping approach: case study of the Old Crow Flats, Yukon, Canada

Abstract. In light of the recent climate warming, monitoring of lake ice in Arctic and subarctic regions is becoming increasingly important. Many shallow Arctic lakes and ponds of thermokarst origin freeze to the bed in the winter months, maintaining the underlying permafrost in its frozen state. However, as air temperatures rise and precipitation increases, fewer lakes are expected to develop bedfast ice. In this work, we propose a novel temporal deep-learning approach to lake ice regime mapping from synthetic aperture radar (SAR) and employ it to study lake ice dynamics in the Old Crow Flats (OCF), Yukon, Canada, over the 1992/1993 to 2020/2021 period. We utilized a combination of Sentinel-1, ERS-1 and ERS-2, and RADARSAT-1 to create an extensive annotated dataset of SAR time series labeled as either bedfast ice, floating ice, or land, which was used to train a temporal convolutional neural network (TempCNN). The trained TempCNN, in turn, allowed us to automatically map lake ice regimes. The classified maps aligned well with the available field measurements and ice thickness simulations obtained with a thermodynamic lake ice model. Reaching a mean overall classification accuracy of 95 %, the TempCNN was determined to be suitable for automated lake ice regime classification. The fraction of bedfast ice in the OCF increased by 11 % over the 29-year period of analysis. Findings suggest that the OCF lake ice dynamics are dominated by lake drainage events, brought on by thermokarst processes accelerated by climate warming, and fluctuations in water level and winter snowfall. Catastrophic drainage and lowered water levels cause surface water area and lake depth to decrease and lake ice to often transition from floating to bedfast ice, while a reduction in snowfall allows for the growth of thicker ice. The proposed lake ice regime mapping approach allowed us to assess the combined impacts of warming, drainage, and changing precipitation patterns on transitions between bedfast and floating-ice regimes, which is crucial to understanding evolving permafrost dynamics beneath shallow lakes and drained basins in thermokarst lowlands such as the OCF.

Observation of a rapid lake‐drainage event in the Arctic: Set‐up and trigger mechanisms, outburst flood behaviour, and broader fluvial impacts

AbstractLakes set in arctic permafrost landscapes can be susceptible to rapid drainage and downstream flood generation. Of many thousands of lakes in northern Alaska, hundreds have been identified as having high drainage potential directly to river systems and 18 such drainage events have been documented since 1955. In 2018 we began monitoring a large lake with high drainage potential as part of a long‐term hydrological observation network designed to evaluate impacts of land use and climate change. In early June 2022, surface water was observed flowing over a 30‐m wide bluff, with active headward erosion of ice‐rich permafrost soils apparent by late June. This overflow point breached rapidly in early July, draining almost the entire lake within 12 h and generating a 191 m3/s flood to a downstream creek. Water level and turbidity sensors and time‐lapse cameras captured this rapid lake‐drainage event at high resolution. A wind‐driven surface seiche and warming waters following ice‐out helped trigger the initial thermomechanical breach. We estimate at least 600 MT of lake sediment was eroded, mobilized, and transported downstream. A flood wave peaking at 42 m3/s arrived 14 h after the initial breach at a river gauge 9‐km downstream. Comparing this event with three other quantified arctic lake‐drainage floods suggests that lake surface area coupled with drainage gradient height can predict outburst flood magnitude. Using this relationship we estimated future flood hazards from the 146 lakes in the Arctic Coastal Plain of northern Alaska (ACP) with high drainage potential, of which 20% are expected to generate outburst floods exceeding 100 m3/s to downstream rivers. This fortunate and detailed drainage‐event observation adds to a growing body of research on the impact of lakes on arctic hydrology, hazard forecasting in a region with an increasing human footprint, and broader processes of landscape evolution in arctic lowlands.

Sedimentary architecture and glacial hydrodynamic significance of the stratified Oak Ridges Moraine, southern Ontario, Canada

High‐quality subsurface data provide new insights into the formation of Oak Ridges Moraine (ORM), an ~80 km3 sequence of stratified meltwater deposits resting >200 m above adjacent Lake Ontario. The ORM sedimentary succession comprises a three‐part regional architecture: (i) ~north–south channel sand–gravel; (ii) channel‐capping rhythmites; and (iii) east–west ridge sediments. The ORM depositional sequence overlies a regional unconformity with a cross‐cutting channel network resulting from ~north–south meltwater floods that transitioned progressively (falling stage) from a ~NNE to ENE flow direction (parallels Lake Ontario depression). Seismic profiles delineate the channels and channel fill characteristics of bank‐to‐bank channel sedimentation of thick gradational gravel–sand–mud sequences. Channel‐capping mud (~100–236 rhythmites) within multiple channels beneath the ORM landform mark a widespread interval of low‐energy, seasonally controlled subglacial pond deposition. During this quiescent period ice‐sheet thickness adjusted to flood‐induced stretching/thinning and re‐profiled slopes. New ice gradients led to east–west flow and deposition of the overlying third element, a sequence of high‐energy confined esker–fan sediments along ORM ridge. Close, sequential timing (~329 varve years) of channel, basin and ridge‐forming architectural elements supports naming this assemblage the ORM formation. Proposed ORM floods are analogous to Icelandic jökulhlaups based on the size, geometry and sedimentology. The observed rhythmite interval between flood events represents a short period (~236 years) of regional meltwater storage prior to east–west ORM flooding. The ORM channel and overlying esker‐fan sediment ridge represent two closely timed meltwater drainage events rather than formation by coalescing ice streams. The scale and timing of the ORM flood events are linked to rapid sea‐level rise, ~13.5 ka BP. This high‐resolution ORM sedimentological record may provide insights into depositional and glaciogenic controls of other large, stratified moraines. The ORM data indicate deposition in response to hydrodynamic events (outbreak floods, re‐profiled ice) rather than direct climate forcing.

Sedimentary facies and clay mineralogy of the late Pleistocene Landsort Deep sediments, Baltic Sea — Implications for the Baltic Ice Lake development

ABSTRACT The Landsort Deep is the deepest part of the Baltic Sea and contains an excellent high-resolution late Pleistocene sediment record suitable to study the retreat history of the southern margin of the Fennoscandian Ice Sheet and the development of the ice-marginal Baltic Ice Lake (BIL) from ~13.5 to 10.5 ka BP. The studied cores are from the lithostratigraphic Units V and VI of Hole M0063C that were recovered during the Integrated Ocean Drilling Program (IODP) Expedition 347. The subdivision and interpretation of the studied BIL record are based on the identification of seven distinct sedimentary facies, grain size characteristics, detrital clay minerals, water and carbon content, and physical properties. The ice-rafted debris (IRD) in the lowermost part of Unit VI indicates a proximal glaciolacustrine environment. The onset of the Fennoscandian Ice Sheet deglaciation during the warm Bølling/Allerød interstadial caused a rise of the proglacial lake water level and the transition from ice-proximal (IRD-containing sediments) to ice-distal varved sediments (rhythmically laminated silty clays) detected within the middle part of Unit VI. This depositional unit of constant sedimentation was followed by the first drainage event of the BIL, which led to the deposition of a distinct clay-rich interval representing erosion and redeposition of the emerged lake bottom sediments. The reduced grain size of the overlying sediments is associated to the onset of cooler climate conditions and a possible ice sheet readvance during the Younger Dryas that resulted in a decrease of meltwater release and sediment availability. The higher sand content and kaolinite peak are associated to a rapid ice retreat and the release of abundant sediment laden meltwater plumes at the end of the Younger Dryas. The final drainage of the BIL enabled the erosion of coastal clay sediments, and the termination of the ice lake stage led to enhanced clay sedimentation to the Landsort Deep.

A greater source of methane from drainage rivers than from rice paddies with drainage practices in southeast China

Rice paddies and rivers draining agricultural watersheds have been documented as two major sources of atmospheric methane (CH4), while parallel flux measurements have been rarely taken on CH4 from rice paddies against the drainage rivers in agricultural watersheds. Moreover, the drainage events during the rice growing season deliver sufficient organic and inorganic nutrients to rivers, which makes there an increasing concerned source of atmospheric CH4. Here, we compared CH4 emissions from rice paddies and the associated rivers draining agricultural watersheds by a parallel field experiment in southeast China. Over the whole rice–growing season (June–November), CH4 fluxes from drainage rivers averaged 25 mg m–2 h–1, which were 82 % higher than those from rice paddies (13.7 mg m–2 h–1). Besides the dependence of riverine CH4 fluxes on water dissolved oxygen (DO), both ecosystem CH4 fluxes were significantly correlated with water/soil dissolved organic carbon (DOC), oxidation-reduction potential (Eh) and the microbial functional genes that drive CH4 production or uptake [methanogens (mcrA and methanogenic archaeal 16S rRNA) and methanotrophs (pmoA)]. Our results highlighted that the CH4 mitigation potential benefited from the drainage practice in rice paddies might have been totally offset by its induced increases in CH4 emissions from rivers draining agricultural watersheds at the regional scale.

Mapping Chini Lake (Pahang, Malaysia) using Sentinel-2 images to determine the effect of acid mine drainage in the pre- to post-COVID-19 restriction period.

Mining activities in the Chini Lake catchment area have been extensive for several years, contributing to acid mine drainage (AMD) events with high concentrations of iron (Fe) and other heavy metals impacting the surface water. However, during the restriction period due to the COVID-19 outbreak, anthropogenic activities have been suspended, which clearly shows a good opportunity for a better environment. Therefore, we aimed to analyze the variation of AMD-associated water pollution in three main zones of the Chini Lake catchment area using Sentinel-2 data for the periods pre-movement control order (MCO), during MCO, and post-MCO from 2019 to 2021. These three zones were chosen due to their proximity to mining areas: zone 1 in the northeastern part, zone 2 in the southeastern part, and zone 3 in the southern part of the Chini Lake area. The acid mine water index (AMWI) was a specific index used to estimate acid mine water. The AMWI values from Sentinel-2 images exhibited that the mean AMWI values in all zones during the MCO period decreased by 14% compared with the pre-MCO period. The spatiotemporal analysis found that the highest polluted zones were recorded in zone 1, followed by zone 3 and zone 2. As compared with during the MCO period, the maximum percentage of increment during post-MCO in all zones was up to 25%. The loosened restriction policy has resulted in more AMD flowing into surface water and increased pollution in Chini Lake. As a whole, our outputs revealed that Sentinel-2 data had a major potential for assessing the AMD-associated pollution of water.

Terminus change of Kaskawulsh Glacier, Yukon, under a warming climate: retreat, thinning, slowdown and modified proglacial lake geometry

AbstractLinks between proglacial lakes and glacier dynamics are poorly understood but are necessary to predict how mountain glaciers will react to a warmer, wetter climate, where such lakes are expected to increase both in number and volume. Here, we examine a long-term (~120 year) record of terminus retreat, thinning and surface velocities from in-situ and remote sensing observations at the terminus of Kaskawulsh Glacier, Yukon, Canada, and determine the impact of a local proglacial hydrological reorganisation on glacier dynamics. After an initial deceleration during the late 1990s, terminus velocities increased at a rate of 3 m a−2from 2000–12, while proglacial Slims Lake area increased simultaneously. The rapid drainage of the lake in May 2016 substantially altered the velocity profile, decreasing annual velocities by 48% within 3 km of the terminus between 2015 and 2021, at an average rate of ~ 12.5 m a−2. A key cause of the rapid drop in glacier motion was a reduction in flotation of the lower part of the glacier terminus after lake drainage. This has important implications for glacier dynamics and provides one of the first assessments of the impacts of a rapid proglacial lake drainage event on local terminus velocities.

Three-Dimensional Analysis of Air-Admission Orifices in Pipelines during Hydraulic Drainage Events

Air valves operate as protection devices in pipelines during drainage processes in order to mitigate vacuum pressures and control the transient flows. Currently, different authors have proposed one-dimensional models to predict the behaviour of orifices during filling and draining events, which offer good numerical results. However, the three-dimensional dynamic behaviour of air-admission orifices during drainage processes has not been studied in depth in the literature. In this research, the effects of air inflow on an orifice installed in a single pipe during drainage events are analysed using a three-dimensional computational fluid dynamics model by testing orifices with diameters of 1.5 and 3.0 mm. This model was validated with different experimental measurements associated to the vacuum pressure, obtaining good fits. The three-dimensional model predicts additional information associated to the aerodynamic effects that occur during the air-admission processes, which is studied. Subsonic flows are observed in different orifices with Mach numbers between 0.18 and 0.30. In addition, it is shown that the larger-diameter orifice ensures a more effective airflow control compared to the smaller-diameter orifice.

Supraglacial lake expansion, intensified lake drainage frequency, and first observation of coupled lake drainage, during 1985–2020 at Ryder Glacier, Northern Greenland

Along the Greenland Ice Sheet margin, supraglacial lakes store and redistribute ice sheet surface run off, and comprise an important potential hydrological link between the ice surface and the base, with ramifications for subglacial drainage systems and ice flow. As a consequence of increasing global mean surface air temperatures, these lakes have been predicted to expand further inland and to affect larger areas of the ice sheet. However, as contemporary dynamics of such supraglacial lake expansion are not well studied, any assessment of their future implications remains afflicted with uncertainty. Here, recent changes in supraglacial lake distribution and expansion, and in their drainage behavior and frequency, are presented for Ryder Glacier, Northern Greenland, as concluded from a remote sensing based analysis. The 35-year time span covered in the analysis allows for the detection of trends in lake processes and ice velocity, which otherwise were found to exhibit large inter-annual variability. It also reveals the first occurrence of a coupled lake drainage event in 2002. By linking supraglacial lake expansion, drainage modes, and drainage frequency to the efficiency of the subglacial drainage system and ice flow on seasonal and decadal timescales, a contribution is made to better understand the complexity of coupled glacio-hydrological processes, and to help reduce uncertainties in predictions of future mass loss from the Greenland Ice Sheet.

The effect of original and initial saturation on residual nonwetting phase capillary trapping efficiency

Injection of supercritical carbon dioxide (CO2) into geological formations is a strategy for both atmospheric greenhouse gas reduction (climate change mitigation) and enhanced oil recovery. To understand CO2 trapping efficiency, the capillary trapping behaviors that immobilize subsurface fluids were analyzed at the pore-scale using pairs of proxy fluids representing the range of in situ nonwetting and wetting fluid properties encountered in geologic storage reservoirs. The pairs of fluids were cycled through imbibition and drainage processes using a flow cell apparatus containing a sintered glass bead column. Computed x-ray microtomography (microCT) was used to identify immobilized nonwetting fluid after imbibition and drainage events. From microCT images, the trapped residual (post-secondary imbibition) nonwetting phase was spatially correlated to both the original (post-primary imbibition) and the initial (post-primary drainage) nonwetting phase; this relationship is referred to here as the original saturation dependence (SO-dependence) and initial saturation dependence (SI-dependence), respectively. Significant trends of decreasing SO- and SI-dependence with increasing wetting and nonwetting fluid phase viscosities were observed. This finding implies that the amount of CO2 injected and ultimately trapped is dependent on the nonwetting phase (e.g. oil or gas) already present in the formation, as well as on the manner in which supercritical CO2 is initially injected, which are factors not considered in current trapping models. To explore the potential effect that varying viscosity and interfacial tension (IFT) might have on this process, we report on a variety of fluid pairs with different viscosities and IFT.

Experimental approach of free flow to pressurized flow in sewer pipes: upstream to downstream pressurization

ABSTRACT Sewer pipe pressurization is a phenomenon that can cause structural damage in urban drainage systems and can result in a public health issue, especially in cities that have combined sewer systems. Studies have classified the phenomena of downstream to upstream pressurization in stages, proposing several mathematical models to predict pressure variations and the effect generated by air entrapment. This paper considers the opposite pressurization process, upstream to downstream, caused by a sudden increase in discharge. Using a physical model in which drainage time and flow rate magnitude are controlled, severe flood drainage events are simulated. Results show that the increases of water level are gradual, and when the depth reaches full pipe flow, pressure increases are also gradual without any transients. The contribution of this work relies on concluding that it is highly improbable that the pressurization process in this case could cause damage problems in the drainage infrastructure. The same results were obtained regardless of the rate of increase in flow and the maximum final pressure in the pipe, which differs from the downstream-upstream case. There are no pressure transients.

Large‐scale lateral saturated soil hydraulic conductivity as a metric for the connectivity of subsurface flow paths at hillslope scale

AbstractLateral saturated soil hydraulic conductivity, Ks,l, is the soil property governing subsurface water transfer in hillslopes, and the key parameter in many numerical models simulating hydrological processes at the hillslope and catchment scales. Likewise, the hydrological connectivity of the lateral flow paths plays a significant role in determining the rate of the subsurface flow at various spatial scales. This study investigates the relationship between Ks,l and hydrological connectivity at the hillslope spatial scale. Ks,l was determined by the subsurface flow rates intercepted by drains and water table depths observed in a well network. The hydrological connectivity was evaluated by the synchronicity among water table peaks, and between these and the peaks of the drained flow. Rainfall and soil moisture were used to investigate the influence of the transient hydrological soil condition on connectivity and Ks,l. As the synchronicity of the water table response between wells increased, the lag times between the peaks of water levels and those of the drained subsurface flow decreased. Moreover, the most synchronic water table rises determined the highest drainage rates. The relationships between Ks,l and water table depths were highly non‐linear, with a sharp increase in the values for water table levels close to the soil surface. Estimated Ks,l values for the full saturated soil were in the order of thousands of mm h−1, suggesting the activation of macropores in the root zone. The Ks,l values determined at the peak of the drainage events were correlated with the indicators of synchronicity. The sum of cumulative rainfall and antecedent soil moisture was correlated with the connectivity indicators and Ks,l. We suggest that, for simulating realistic processes at the hillslope scale, the hydrological connectivity could be implicitly considered in hydrological modelling through an evaluation of Ks,l at the same spatial scale.

The Last Glacial Maximum and Deglacial History of the Seno Skyring Ice Lobe (52°S), Southern Patagonia

There are still many uncertainties about the climatic forcing that drove the glacier fluctuations of the Patagonian Ice Sheet (PIS, 38–55°S) during the last glacial period. A key source of uncertainty is the asynchrony of ice lobe fluctuations between the northern, central, and southern PIS. To fully understand the regional trends requires careful mapping and extensive geochronological studies. This paper presents geomorphological and geochronological reconstructions of the glacial and deglacial landforms formed during the last glacial period at the Seno Skyring lobe, southernmost Patagonia (52°S, 71°W). We present a detailed geomorphological map, where we identify two moraine systems. The outer and older is named Laguna Blanca (LB) and the inner Río Verde (RV). The LB moraines were built subaerially, whereas parts of the RV were deposited subaqueously under the palaeo lake Laguna Blanca, which developed during deglaciation. We conducted surface exposure 10Be dating methods on boulder samples collected from LB and RV glacial margins. The moraine LB III and LB IV formed at 26.3 ± 2.3 ka (n = 5) and 24.3 ± 0.9 ka (n = 3), respectively. For the inner RV moraine, we obtained an age of 18.7 ± 1.5 ka (n = 6). For the palaeo Laguna Blanca evolution, we performed 10Be exposure ages on shoreline berms and optically stimulated luminesce dating to constrain the lake levels, and 10Be depth profile dating on an outwash deposit formed by a partial lake drainage event, which occurred at 22 ± 3 ka. For the RV moraine deglaciation, we performed radiocarbon dating of basal sediments in a peat bog, which indicates that the glacier retreated from the terminal RV moraine by at least c. 16.4 cal kyr BP. Our moraine geochronology shows an asynchrony in the maximum extents and a different pattern of ice advances between neighbouring lobes in southern Patagonia. We speculate that this may be due, at least in part, to the interaction between topography and the precipitation carried by the southern westerly wind belt. However, we found broad synchrony of glacial readvances contemporaneous with the RV moraine.

Remotely sensed time series of rapid terrace formation: Laguna del Viedma valley (Patagonia)

The Patagonian Andes were affected by a range of geophysical drivers of landscape incision during the Last Glacial Interglacial Transition and Early Holocene. Deciphering drivers of river system response during this period is complex, and magnitudes and timescales of landscape change are poorly constrained. Herein, a remotely sensed time series of modern terrace formation is investigated from the Laguna del Viedma valley as a modern analogue of Late Quaternary landscape evolution in Patagonia. The aim of the research was to constrain the timing of terrace formation following lake level fall of the Laguna del Viedma over a 35 year period from 1985-2019. The objectives were to: 1) use satellite imagery from 1985-2019 to document glacier and lake changes in the study area; 2) map landforms of the Laguna del Viedma valley; and 3) analyse terrace elevations. In total seven terrace surfaces were distinguished, with the oldest four pre-dating the ALOS PALSAR DEM (February 2000) used. Landform evidence shows the highest, and vegetated, T1 terrace surface (+40-75 m) grades to the highest lake level and was likely the elevation of the valley floor during Holocene neoglacials. Viedma glacier recession then led to a phase of lake regressions/transgressions with an overall trend of lake level fall. The DEM shows ~20 m incision from the 1985 floodplain level (T3) to the T4 level floodplain by 2000. This constrains a minimum rate of incision of 1.33 m/yr, however, the satellite time series demonstrates rapid T3 terrace formation, with evidence for mass movements contributing to lateral terrace erosion by 1986. The implications of the data are discussed within the context of the Late Quaternary palaeohydrology of Patagonia where lake level falls of 10s to 100s of metres occurred within many large river systems of the Patagonian Andes from 42-52⁰ S. The data herein demonstrate that base level falls from sudden lake drainage events were likely a major driver of rapid landscape change in Patagonia during deglaciation.

Antarctic Basal Water Storage Variation Inferred from Multi-Source Satellite Observation and Relevant Models

Antarctic basal water storage variation (BWSV) refers to mass changes of basal water beneath the Antarctic ice sheet (AIS). Identifying these variations is critical for understanding Antarctic basal hydrology variations and basal heat conduction, yet they are rarely accessible due to a lack of direct observation. This paper proposes a layered gravity density forward/inversion iteration method to investigate Antarctic BWSV based on multi-source satellite observations and relevant models. During 2003–2009, BWSV increased at an average rate of 43 ± 23 Gt/yr, which accounts for 29% of the previously documented total mass loss rate (−76 ± 20 Gt/yr) of AIS. Major uncertainty arises from satellite gravimetry, satellite altimetry, the glacial isostatic adjustment (GIA) model, and the modelled basal melting rate. We find that increases in basal water mainly occurred in regions with widespread active subglacial lakes, such as the Rockefeller Plateau, Siple Coast, Institute Ice Stream regions, and marginal regions of East Antarctic Ice Sheet (EAIS), which indicates the increased water storage in these active subglacial lakes, despite the frequent water drainage events. The Amundsen Sea coast experienced a significant loss during the same period, which is attributed to the basal meltwater discharging into the Amundsen Sea through basal channels.

Repeated Tidally Induced Hydrofracture of a Supraglacial Lake at the Amery Ice Shelf Grounding Zone

AbstractSurface melting and lakes are common to Antarctic ice shelves, and their existence and drainages have been invoked as a precursor for ice shelf collapse. Here, we present satellite observations over 2014–2020 of repeated, rapid drainages of a supraglacial lake at the grounding zone of Amery Ice Shelf, East Antarctica. Post‐drainage imagery in 2018 reveals lake bottom features characteristic of rapid, vertical lake drainage. Observed lake volumes indicate drainages are not associated with a threshold meltwater volume. Instead, drainages typically coincide with periods of high daily tidal amplitude, suggesting hydrofracture is assisted by tidally forced ice flexure inherent to the ice shelf grounding zone. Combined with observations of widespread grounding zone lake drainages on Amery, these findings indicate ice shelf meltwater accumulation may be inhibited by grounding zone drainage events, thus representing a potential stabilizing mechanism despite enhanced melting common to these regions.

Seismic Mapping of Subglacial Hydrology Reveals Previously Undetected Pressurization Event

AbstractUnderstanding the dynamic response of glaciers to climate change is vital for assessing water resources and hazards, and subglacial hydrology is a key player in glacier systems. Traditional observations of subglacial hydrology are spatially and temporally limited, but recent seismic deployments on and around glaciers show the potential for comprehensive observation of glacial hydrologic systems. We present results from a high‐density seismic deployment spanning the surface of Lemon Creek Glacier, Alaska. Our study coincided with a marginal lake drainage event, which served as a natural experiment for seismic detection of changes in subglacial hydrology. We observed glaciohydraulic tremor across the surface of the glacier that was generated by the subglacial hydrologic system. During the lake drainage, the relative changes in seismic tremor power and water flux are consistent with pressurization of the subglacial system of only the upper part of the glacier. This event was not accompanied by a significant increase in glacier velocity; either some threshold necessary for rapid basal motion was not attained, or, plausibly, the geometry of Lemon Creek Glacier inhibited speedup. This pressurization event would have likely gone undetected without seismic observations, demonstrating the power of cryoseismology in testing assumptions about and mapping the spatial extent of subglacial pressurization.

Single midseason drainage events decrease global warming potential without sacrificing grain yield in flooded rice systems

Rice (Oryza sativa L.) cultivation is an important part of global food security, yet it is also responsible for a significant portion of agricultural greenhouse gas (GHG) emissions, particularly methane (CH4). Midseason drainage of flooded rice fields can decrease CH4 emissions, but the magnitude of CH4 reduction and its effect on grain yield are variable due to variation in the timing and soil-drying severity of drainage across studies. Therefore, in this two-year study, we aimed to quantify the effect of timing and severity of a single midseason drainage event on seasonal GHG emissions and grain yields, compared to a continuously flooded (CF) control. Treatments varied in terms of soil-drying severity (low, medium, and high, corresponding to approximately 5, 8, and 12 days of drying, respectively) and the timing of when drainage events occurred (between 34–49 and 45–59 days after seeding, or roughly between tillering and panicle initiation). Soil moisture parameters (perched water table, volumetric water content, gravimetric water content (GWC), and soil water potential), soil mineral nitrogen, CH4 and nitrous oxide (N2O) emissions, grain yield, and yield components were all quantified. Midseason drainage reduced seasonal CH4 emissions by 38–66%, compared to the CF control. Seasonal CH4 emissions decreased with increasing drain severity, and for every 1% reduction in soil GWC during the drainage period, seasonal CH4 emissions were reduced by 2.5%. The timing of drainage had no significant impact on CH4 emissions. Emissions of N2O were low (average = 0.035 kg N2O-N ha−1) and accounted for only 0.5% of the seasonal global warming potential (GWP) across all drainage treatments. Within each year, drainage did not significantly affect grain yield compared to the CF control. Additionally, midseason drainage reduced both GWP and yield-scaled GWP by approximately the same amount as seasonal CH4 emissions, as N2O emissions were minimal and yields were similar across treatments. These results indicate that midseason drainage may be a viable GHG mitigation practice in flooded rice systems with limited risk for yield reduction, however, this practice should also be further tested under a broad range of soil types and different environments to determine its widespread adoptability.