Subsurface Meltwater Research Articles

Numerical Modelling of the Influence of Tidal Stresses on Fracture Patterns on the Surface of Europa

AbstractJupiter’s satellite Europa is covered by an ice shell exhibiting many surface features, including linear structures called lineae, which in this work are treated as fractures. A three-dimensional finite-element simulator is used to numerically model fracture nucleation, growth, and interaction, assuming the ice is an isotropic, linear elastic medium. Tidal stresses are exerted upon the ice through the Jovian orbital relationship. These stresses are calculated using a closed-form model derived from first principles. The fractures grow in response to stress concentration around their tips, and a damage criterion models the weakening of the ice matrix. Three-dimensional non-planar multiple fracture growth is modelled as a function of geometric multi-modal stress intensity factors computed at the fracture tips. Fracturing is evaluated over multi-scale periods, from days to millions of years, thus capturing multiple tidal effects. Fracture behaviour is modelled across the Europan surface in one domain. The patterns are dense clusters of lineae about the stress maxima with diffuse fracturing in outlying regions. Fractures are also modelled in the vicinity of subsurface meltwater lenses, where fractures form parallel to the surface in contrast to the usual perpendicular orientation. The resultant fracture patterns are qualitatively compared against images from NASA’s Galileo mission. This work contributes to the understanding of Europan lineae by illustrating how they behave in a fracture mechanics framework, and suggests interesting results regarding lineae interaction with meltwater lenses. This work is also a proof of concept for this modelling approach, and will serve as the framework for future work.

Pathways of Inter‐Basin Exchange From the Bellingshausen Sea to the Amundsen Sea

AbstractThe West Antarctic Ice Sheet is experiencing rapid thinning of its floating ice shelves, largely attributed to oceanic basal melt. Numerical models suggest that the Bellingshausen Sea has a key role in setting water properties in the Amundsen Sea and further downstream. Yet, observations confirming these pathways of volume and tracer exchange between coast and shelf break and their impact on inter‐sea exchange remain sparse. Here we analyze the circulation and distribution of glacial meltwater at the boundary between the Bellingshausen Sea and the Amundsen Sea using a combination of glider observations from January 2020 and hydrographic data from instrumented seals. Meltwater distributions over previously unmapped western regions of the continental shelf and slope reveal two distinct meltwater cores with different optical backscatter properties. At Belgica Trough, a subsurface meltwater peak is linked with hydrographic properties from Venable Ice Shelf. West of Belgica Trough, the vertical structure of meltwater concentration changes, with peak values occurring at greater depths and denser isopycnals. Hydrographic analysis suggests that the western (deep) meltwater core is supplied from the eastern part of Abbot Ice Shelf, and is exported to the shelf break via a previously‐overlooked bathymetric trough (here named Seal Trough). Hydrographic sections constructed from seal data reveal that the Antarctic Coastal Current extends west past Belgica Trough, delivering meltwater to the Amundsen Sea. Each of these circulation elements has distinct dynamical implications for the evolution of ice shelves and water masses both locally and downstream, in the Amundsen Sea and beyond.

Spatial and temporal differences in surface and subsurface meltwater distribution over Greenland ice sheet using multi-frequency passive microwave observations

Increasingly larger portions of the Greenland ice sheet are undergoing seasonal melting-refreeze cycles due to global climate warming. The cycle begins with the arrival of high temperatures and increased solar radiation in the spring and summer seasons generating meltwater on the ice sheet's surface. Meltwater percolates to deeper ice layers, either refreezing within the firn, creating longer-term meltwater pockets (firn aquifers), or generating peripheral runoff. Depending on the location and climate, the refreeze duration, the depth of infiltration, and meltwater persistence are temporally and spatially complex. Our recent study showed that multi-frequency passive microwave measurements in the 1.4 GHz to 36.5 GHz range effectively distinguished seasonal meltwater between the immediate surface and deeper firn layers at an experiment site in the accumulation zone of the southwestern Greenland ice sheet. Here, we further explored the vertically and horizontally polarized multi-frequency melt response at the pan-Greenland scale. We employed 1.4 GHz brightness temperature (TB) measurements from the NASA Soil Moisture Active Passive (SMAP) satellite and 6.9, 10.7, 18.9, and 36.5 GHz TB measurements from the JAXA Global Change Observation Mission-Water Shizuku (GCOM-W) satellite. The results show that the frequency-dependent response was consistent across the ice sheet. The multi-frequency melt indications match with lasting seasonal subsurface meltwater with delayed refreezing compared to the surface. These results suggest persistent seasonal subsurface meltwater occurrences that are spatially and temporally significant but concealed from the high-frequency observations. Retrieving the meltwater evolution in snow and firn presents a complex problem; this work represents an initial step toward developing an ice-sheet-wide algorithm for more comprehensive retrieval of the meltwater profile.

Ice Sheet Surface and Subsurface Melt Water Discrimination Using Multi‐Frequency Microwave Radiometry

AbstractFor understanding englacial hydrology and its impact on ice sheet mass balance, observations of the liquid water content (LWC) within the ice sheets are needed. We combined 1.4–10.7 GHz passive microwave measurements with traditional 18.7–36.5 GHz measurements to detect subsurface LWC. In situ measurements from the DYE‐2 experiment site in Greenland and a modeled LWC at this site were used to calibrate and validate the method. Our analysis showed sensitivity of the lower microwave frequencies to LWC in surface and subsurface layers down to at least 2 m, enabling detection of seasonal subsurface LWC and its refreezing. A simplified retrieval detected a delayed refreezing of subsurface LWC following surface freezing, while also capturing total seasonal meltwater production. These advancements open the door to detection of subsurface meltwater and refreezing twice a day at pan‐Greenland scale, thereby enabling improved estimates of ice sheet contributions to global sea level rise.

An empirical algorithm to map perennial firn aquifers and ice slabs within the Greenland Ice Sheet using satellite L-band microwave radiometry

Abstract. Perennial firn aquifers are subsurface meltwater reservoirs consisting of a meters-thick water-saturated firn layer that can form on spatial scales as large as tens of kilometers. They have been observed within the percolation facies of glaciated regions experiencing intense seasonal surface melting and high snow accumulation. Widespread perennial firn aquifers have been identified within the Greenland Ice Sheet (GrIS) via field expeditions, airborne ice-penetrating radar surveys, and satellite microwave sensors. In contrast, ice slabs are nearly continuous ice layers that can also form on spatial scales as large as tens of kilometers as a result of surface and subsurface water-saturated snow and firn layers sequentially refreezing following multiple melting seasons. They have been observed within the percolation facies of glaciated regions experiencing intense seasonal surface melting but in areas where snow accumulation is at least 25 % lower as compared to perennial firn aquifer areas. Widespread ice slabs have recently been identified within the GrIS via field expeditions and airborne ice-penetrating radar surveys, specifically in areas where perennial firn aquifers typically do not form. However, ice slabs have yet to be identified from space. Together, these two ice sheet features represent distinct, but related, sub-facies within the broader percolation facies of the GrIS that can be defined primarily by differences in snow accumulation, which influences the englacial hydrology and thermal characteristics of firn layers at depth. Here, for the first time, we use enhanced-resolution vertically polarized L-band brightness temperature (TVB) imagery (2015–2019) generated using observations collected over the GrIS by NASA's Soil Moisture Active Passive (SMAP) satellite to map perennial firn aquifer and ice slab areas together as a continuous englacial hydrological system. We use an empirical algorithm previously developed to map the extent of Greenland's perennial firn aquifers via fitting exponentially decreasing temporal L-band signatures to a set of sigmoidal curves. This algorithm is recalibrated to also map the extent of ice slab areas using airborne ice-penetrating radar surveys collected by NASA's Operation IceBridge (OIB) campaigns (2010–2017). Our SMAP-derived maps show that between 2015 and 2019, perennial firn aquifer areas extended over 64 000 km2, and ice slab areas extended over 76 000 km2. Combined together, these sub-facies are the equivalent of 24 % of the percolation facies of the GrIS. As Greenland's climate continues to warm, seasonal surface melting will increase in extent, intensity, and duration. Quantifying the possible rapid expansion of these sub-facies using satellite L-band microwave radiometry has significant implications for understanding ice-sheet-wide variability in englacial hydrology that may drive meltwater-induced hydrofracturing and accelerated ice flow as well as high-elevation meltwater runoff that can impact the mass balance and stability of the GrIS.

The potential of synthetic aperture radar interferometry for assessing meltwater lake dynamics on Antarctic ice shelves

Abstract. Surface meltwater drains on several Antarctic ice shelves, resulting in surface and sub-surface lakes that are potentially critical for the ice shelf collapse. Despite these phenomena, our understanding and assessment of the drainage and refreezing of these lakes is limited, mainly due to lack of field observations and to the limitations of optical satellite imagery during polar night and in cloudy conditions. This paper explores the potential of backscatter intensity and of interferometric coherence and phase from synthetic aperture radar (SAR) imagery as an alternative to assess the dynamics of meltwater lakes. In four case study regions over Amery and Roi Baudouin ice shelves, East Antarctica, we examine spatial and temporal variations in SAR backscatter intensity and interferometric (InSAR) coherence and phase over several lakes derived from Sentinel-1A/B C-band SAR imagery. Throughout the year, the lakes are observed in a completely frozen state, in a partially frozen state with a floating ice lid and as open-water lakes. Our analysis reveals that the meltwater lake delineation is challenging during the melting period when the contrast between melting snow and lakes is indistinguishable. Despite this finding, we show using a combination of backscatter and InSAR observations that lake dynamics can be effectively captured during other non-summertime months. Moreover, our findings highlight the utility of InSAR-based observations for discriminating between refrozen ice and sub-surface meltwater and indicate the potential for phase-based detection and monitoring of rapid meltwater drainage events. The potential of this technique to monitor these meltwater change events is, however, strongly determined by the satellite revisit interval and potential changes in scattering properties due to snowfall or melt events.

Spatiotemporal Drivers of Hydrochemical Variability in a Tropical Glacierized Watershed in the Andes

AbstractThere is a critical knowledge gap about how glacier retreat in remote and rapidly warming tropical montane watersheds will impact solute export, which has implications for downstream geochemical cycling and ecological function. Because tropical glacierized watersheds are often uniquely characterized by year‐round ablation, upslope vegetation migration, and significant groundwater flow, baseline understanding is needed of how spatiotemporal variables within these watersheds control outlet hydrochemistry. We implemented a recently developed reactive transport watershed model, BioRT‐Flux‐PIHM, for a sub‐humid glacierized watershed in the Ecuadorian Andes with young volcanic soils and fractured bedrock. We found a unique simulated concentration and discharge (C‐Q) pattern that was mostly chemostatic but superimposed by dilution episodes. The chemostatic background was attributed to large simulated contributions of groundwater (subsurface lateral flow) to streamflow, of which a notable fraction (37%) comprised infiltrated ice‐melt. Relatively constant concentrations were further maintained in the model because times and locations of lower mineral surface wetting and dissolution were offset by concentrating effects of greater evapotranspiration. Ice‐melt did not all infiltrate in simulations, especially during large precipitation events, when high surface runoff contributions to discharge triggered dilution episodes. In a model scenario without ice‐melt, major ion concentrations, including Na+, Ca2+, and Mg2+, became more strongly chemostatic and higher, but weathering rates decreased, attenuating export by 23%. We expect this reduction to be exacerbated by higher evapotranspiration and drier conditions with expanded vegetation. This work brings to light the importance of subsurface meltwater flow, ecohydrological variability, and interactions between melt and precipitation for controlling hydrochemical processes in tropical watersheds with rapidly retreating glaciers.

Export of Strongly Diluted Greenland Meltwater From a Major Glacial Fjord

AbstractThe Greenland Ice Sheet has been, and will continue, losing mass at an accelerating rate. The influence of this anomalous meltwater discharge on the regional and large‐scale ocean could be considerable but remains poorly understood. This uncertainty is in part a consequence of challenges in observing water mass transformation and meltwater spreading in coastal Greenland. Here we use tracer observations that enable unprecedented quantification of the export, mixing, and vertical distribution of meltwaters leaving one of Greenland's major glacial fjords. We find that the primarily subsurface meltwater input results in the upwelling of the deep fjord waters and an export of a meltwater/deepwater mixture that is 30 times larger than the initial meltwater release. Using these tracer data, the vertical structure of Greenland's summer meltwater export is defined for the first time showing that half the meltwater export occurs below 65 m.

Evidence for an explosive origin of central pit craters on Mars

Kilometer-scale pits are nested in the centers of many impact craters on Mars as well as on icy satellites. They have been inferred to form in the presence of a water–ice rich substrate; however, the process(es) responsible for their formation is still debated. Previous models invoke origins by either explosive excavation of potentially water-bearing crustal material, or by subsurface drainage of meltwater and/or collapse. If explosive excavation forms central pits, pit-derived ejecta should be draped around the pits, whereas internal collapse should not deposit significant material outside pit rims. Using visible wavelength images from the Mars Reconnaissance Orbiter (MRO) Context Camera (CTX) and High Resolution Imaging Science Experiment (HiRISE) instruments and thermal infrared images from the Odyssey Thermal Emission Imaging System (THEMIS) instrument, we conducted a survey to characterize, in detail, the global population of central pits in impact craters ⩾10km in diameter. We specifically examined the morphology and thermophysical characteristics of the pits for evidence of pit ejecta. Our analysis of thermal images suggests that coarse-grained materials are distributed proximally around many central pits on the floors of their host craters. The decrease in average grain size with distance from pit rims is consistent with pit-derived ejecta. These observations and interpretations better support an explosive origin for central pits on Mars than they do an origin of subsurface meltwater drainage and collapse of the overlying substrate. A major weakness to previous explosive central pit formation models is the inability for them to form pits late enough in the impact process to be preserved. To address this, we present an alternative “melt contact model” where a central uplift brings ice-bearing substrate into contact with impact melt to generate steam explosions and excavate central pits during the impact modification stage. Theoretical calculations show that more than enough thermal energy is available via impact melt from the host crater to form central pits by steam explosions, and such explosions would require only a modest amount (2–6% by volume) of uplifted water–ice. We therefore propose that central pits on Mars could have formed explosively by the interaction of impact melt and subsurface water–ice.

Three types of liquid water in icy surfaces of celestial bodies

It is shown that, at temperatures far below the triple point and under appropriate conditions, liquid water can stably or temporarily exist in upper ice-covered surfaces of planetary bodies (like Mars) in three different types: (i) undercooled interfacial water (due to freezing point depression by van der Waals forces and “premelting”), (ii) water in brines (due to freezing point depression in solutions), and (iii) sub-surface melt water (due to a solid-state greenhouse effect driven heating). The physics behind and the related conditions for these liquid waters to evolve and to exist, and possibly related consequences, are discussed. These calculations are mainly made in view of the possible presence of these sub-surface liquids in the upper surface of the present Mars.

Sensitivity of ice-cemented Antarctic soils to greenhouse-induced thawing: Are terrestrial archives at risk?

We modeled the sensitivity of six ice-cemented slope deposits from the western McMurdo Dry Valleys, Antarctica to failure by shallow, thaw-induced planar sliding. The deposits examined have purportedly remained physically stable, without morphologic evidence for downslope movement, for millions of years. Could they fail in the near future from greenhouse-induced warming? To address this question, we first prescribed various increases in mean summertime soil surface temperature (MSSST) and modeled numerically the resultant changes in soil thaw depths using a one-dimensional heat diffusion equation (including the effects of latent heat of fusion). Second, we calculated the minimum thaw depths required to facilitate failure by shallow planar sliding for each deposit; for all numerical simulations, we maintained present soil-moisture conditions and used a Mohr–Coulomb-based equation of safety factor. Third, we calculated the rate of subsurface meltwater flow assuming Darcy's Law. Our results show that although most deposits contain sufficient subsurface ice to induce sliding upon thawing, lateral rates of water flow of as much as ∼ 40 m/day for some colluvial deposits prohibit the build-up of requisite pore pressures for failure. On the other hand, silty deposits, that contain gravimetric water ≥ 15%, occur on slopes > 20°, and possess low hydraulic conductivities (∼ 30–60 cm/day), common in the Dry Valleys region, could fail if MSSST, and by inference mean summertime atmospheric temperatures, increase by 4 to 9 °C. This temperature increase is similar to that predicted to occur from greenhouse-induced warming in Antarctica over the next century.

Geogenic sources and sinks of trace metals in the Larsemann Hills, East Antarctica: Natural processes and human impact

To evaluate the extent of human impact on a pristine Antarctic environment, natural baseline levels of trace metals have been established in the basement rocks of the Larsemann Hills, East Antarctica. From a mineralogical and geochemical point of view the Larsemann Hills basement is relatively homogeneous, and contains high levels of Pb, Th and U. These may become soluble during the relatively mild Antarctic summer and be transported to lake waters by surface and subsurface melt water. Melt waters may also be locally enriched in V, Cr, Co, Ni, Zn and Sn derived from weathering of metabasite pods. With a few notable exceptions, the trace metal concentrations measured in the Larsemann Hills lake waters can be entirely accounted for by natural processes such as sea spray and surface melt water input. Thus, the amount of trace metals released by weathering of basement lithologies and dispersed into the Larsemann Hills environment, and presumably in similar Antarctic environments, is, in general, not negligible, and may locally be substantial. The Larsemann Hills sediments are coarse-grained and contain minute amounts of clay-size particles, although human activities have contributed to the generation of fine-grained material at the most impacted sites. Irrespective of their origin, these small amounts of fine-grained clastic sediments have a relatively small surface area and charge, and are not as effective metal sinks as the abundant, thick cyanobacterial algal mats that cover the lake floors. Thus, the concentration of trace metals in the Larsemann Hills lake waters is regulated by biological activity and thawing–freezing cycles, rather than by the type and amount of clastic sediment supply.

Antarctic Surface and Subsurface Snow and Ice Melt Fluxes

Abstract This paper presents modeled surface and subsurface melt fluxes across near-coastal Antarctica. Simulations were performed using a physical-based energy balance model developed in conjunction with detailed field measurements in a mixed snow and blue-ice area of Dronning Maud Land, Antarctica. The model was combined with a satellite-derived map of Antarctic snow and blue-ice areas, 10 yr (1991–2000) of Antarctic meteorological station data, and a high-resolution meteorological distribution model, to provide daily simulated melt values on a 1-km grid covering Antarctica. Model simulations showed that 11.8% and 21.6% of the Antarctic continent experienced surface and subsurface melt, respectively. In addition, the simulations produced 10-yr averaged subsurface meltwater production fluxes of 316.5 and 57.4 km3 yr−1 for snow-covered and blue-ice areas, respectively. The corresponding figures for surface melt were 46.0 and 2.0 km3 yr−1, respectively, thus demonstrating the dominant role of subsurface over surface meltwater production. In total, computed surface and subsurface meltwater production values equal 31 mm yr−1 if evenly distributed over all of Antarctica. While, at any given location, meltwater production rates were highest in blue-ice areas, total annual Antarctic meltwater production was highest for snow-covered areas due to its larger spatial extent. The simulations also showed higher interannual meltwater variations for surface melt than subsurface melt. Since most of the produced meltwater refreezes near where it was produced, the simulated melt has little effect on the Antarctic mass balance. However, the melt contribution is important for the surface energy balance and in modifying surface and near-surface snow and ice properties such as density and grain size.

Meltwater production in Antarctic blue-ice areas: sensitivity to changes in atmospheric forcing

In the near coastal regions of Dronning Maud Land, Antarctica, below-surface ice-melt in blue-ice areas has been observed. The low scattering coefficients of the large-grained blue-ice allow penetration of solar radiation, thus providing an energy source below the ice surface. The sub-surface meltwater is significant enough to show up on remote-sensing imagery in the form of ice-covered lakes. Adjacent snow-accumulation areas have much higher scattering coefficients and consequently limit solar radiation penetration in these regions. These snow and ice surfaces are generally below freezing, and little surface melting occurs. To assess the response of these melt features to changes in atmospheric forcings such as cloudiness, air temperature, and snow accumulation, a physically-based model of the coupled atmosphere, radiation, snow, and blue-ice system has been developed. The model consists of a heat transfer equation with a spectrally-dependent solar-radiation source term. The penetration of radiation into the snow and blue-ice depends on the surface albedo, and the snow and blue-ice grain size and density. Model simulations show that ice melt occurring in this area is sensitive to potential variations in atmospheric forcing. Under certain conditions more traditional surface melting occurs and, under other conditions, the existing melt processes can be shut down completely. In light of the sensitivity of this system to variations in atmospheric forcing, and the ability to view melt-related features using remote sensing, a tool exists to efficiently monitor variations in Antarctic coastal climate.

A premilinary assessment of the chemical and hydrological interaction of acidic snowmelt water with the terrestrial portion of a Canadian Shield Catchment

This study examines certain aspects of the hydrological and chemical interaction of springmelt runoff with the terrestrial portion of a Canadian Shield catchment. Though NO3 − is the dominant anion in the snowpack for most of the winter, SO4 −2 deposition increases at the end of the winter season (March) such that it is the dominant anion in snowpack meltwater runoff. Pronounced formation office lenses in the snowpack are responsible for notable movement of H+, NO3 − and SO4 ions downslope prior to contact with the ground. Subsurface meltwater flow through the lower portion of the slope is much more acidic (pH 5.18 ) than that flowing within the subsurface in the upper slope ( pH 6.44 ).