Oceanic Melting Research Articles

Bathymetry of Northwest Greenland Using “Ocean Melting Greenland” (OMG) High-Resolution Airborne Gravity and Other Data

Marine-terminating glaciers dominate the evolution of the Greenland Ice Sheet (GrIS) and its contribution to sea-level rise. Widespread glacier acceleration has been linked to the warming of ocean waters around the periphery of Greenland but a lack of information on the bathymetry of the continental shelf and glacial fjords has limited our ability to understand how subsurface, warm, salty ocean waters of Atlantic origin (AW) reach the glaciers and melt them from below. Here, we employ high-resolution, airborne gravity data (AIRGrav) in combination with multibeam echo sounding (MBES) data, to infer the bathymetry of the coastal areas of Northwest Greenland for NASA’s Ocean Melting Greenland (OMG) mission. High-resolution, AIRGrav data acquired on a 2 km spacing, 150 m ground clearance, with 1.5 mGal crossover error, is inverted in three dimensions to map the bathymetry. To constrain the inversion away from MBES data, we compare two methods: one based on the Direct Current (DC) shift of the gravity field (absolute minus observed gravity) and another based on the density of the bedrock. We evaluate and compare the two methods in areas with complete MBES coverage. We find the lowest standard error in bed elevation (±60 m) using the DC shift method. When applied to the entire coast of Northwest Greenland, the three-dimensional inversion reveals a complex network of connected sea bed channels, not known previously, that provide natural and varied pathways for AW to reach the glaciers across the continental shelf. The study demonstrates that the gravity approach offers an efficient and practical alternative to extensive ship mapping in ice-filled waters to obtain information critical to understanding and modeling ice-ocean interaction along ice sheet margins.

Localized Plumes Drive Front‐Wide Ocean Melting of A Greenlandic Tidewater Glacier

AbstractRecent acceleration of Greenland's ocean‐terminating glaciers has substantially amplified the ice sheet's contribution to global sea level. Increased oceanic melting of these tidewater glaciers is widely cited as the likely trigger, and is thought to be highest within vigorous plumes driven by freshwater drainage from beneath glaciers. Yet melting of the larger part of calving fronts outside of plumes remains largely unstudied. Here we combine ocean observations collected within 100 m of a tidewater glacier with a numerical model to show that unlike previously assumed, plumes drive an energetic fjord‐wide circulation which enhances melting along the entire calving front. Compared to estimates of melting within plumes alone, this fjord‐wide circulation effectively doubles the glacier‐wide melt rate, and through shaping the calving front has a potential dynamic impact on calving. Our results suggest that melting driven by fjord‐scale circulation should be considered in process‐based projections of Greenland's sea level contribution.

Ocean‐Induced Melt Triggers Glacier Retreat in Northwest Greenland

AbstractIn recent decades, tidewater glaciers in Northwest Greenland contributed significantly to sea level rise but exhibited a complex spatial pattern of retreat. Here we use novel observations of bathymetry and water temperature from NASA's Ocean Melting Greenland mission to quantify the role of warm, salty Atlantic Water in controlling the evolution of 37 glaciers. Modeled ocean‐induced undercutting of calving margins compared with ice advection and ice front retreat observed by satellites from 1985 to 2015 indicate that 35 glaciers retreated when cumulative anomalies in ocean‐induced undercutting rose above the range of seasonal variability of calving‐front positions, while two glaciers standing on shallow sills and colder water did not retreat. Deviations in the observed timing of retreat are explained by residual uncertainties in bathymetry, inefficient mixing of waters in shallow fjords, and the presence of small floating sections. Overall, warmer ocean temperature triggered the retreat, but calving processes dominate ablation (71%).

Modelled fracture and calving on the Totten Ice Shelf

Abstract. The Totten Ice Shelf (IS) has a large drainage basin, much of which is grounded below sea level, leaving the glacier vulnerable to retreat through the marine ice sheet instability mechanism. The ice shelf has also been shown to be sensitive to changes in calving rate, as a very small retreat of the calving front from its current position is predicted to cause a change in flow at the grounding line. Therefore understanding the processes behind calving on the Totten IS is key to predicting its future sea level rise contribution. Here we use the Helsinki Discrete Element Model (HiDEM) to show that not all of the fractures visible at the front of the Totten IS are produced locally, but that the across-flow basal crevasses, which are part of the distinctive cross-cutting fracture pattern, are advected into the calving front area from upstream. A separate simulation of the grounding line shows that re-grounding points may be key areas of basal crevasse production, and can produce basal crevasses in both an along- and across-flow orientation. The along-flow basal crevasses at the grounding line may be a possible precursor to basal channels, while we suggest the across-flow grounding-line fractures are the source of the across-flow features observed at the calving front. We use two additional models to simulate the evolution of basal fractures as they advect downstream, demonstrating that both strain and ocean melt have the potential to deform narrow fractures into the broad basal features observed near the calving front. The wide range of factors which influence fracture patterns and calving on this glacier will be a challenge for predicting its future mass loss.

Evidence for a sulfur-undersaturated lunar interior from the solubility of sulfur in lunar melts and sulfide-silicate partitioning of siderophile elements

Sulfur concentrations at sulfide saturation (SCSS) were determined for a range of low- to high-Ti lunar melt compositions (synthetic equivalents of Apollo 14 black and yellow glass, Apollo 15 green glass, Apollo 17 orange glass and a late-stage lunar magma ocean melt, containing between 0.2 and 25 wt.% TiO2) as a function of pressure (1–2.5 GPa) and temperature (1683–1883 K). For the same experiments, sulfide-silicate partition coefficients were derived for elements V, Cr, Mn, Co, Cu, Zn, Ga, Ge, As, Se, Mo, Sn, Sb, Te, W and Pb. The SCSS is a strong function of silicate melt composition, most notably FeO content. An increase in temperature increases the SCSS and an increase in pressure decreases the SCSS, both in agreement with previous work on terrestrial, lunar and martian compositions. Previously reported SCSS values for high-FeO melts were combined with the experimental data reported here to obtain a new predictive equation to calculate the SCSS for high-FeO lunar melt compositions. Calculated SCSS values, combined with previously estimated S contents of lunar low-Ti basalts and primitive pyroclastic glasses, suggest their source regions were not sulfide saturated. Even when correcting for the currently inferred maximum extent of S degassing during or after eruption, sample S abundances are still > 700 ppm lower than the calculated SCSS values for these compositions. To achieve sulfide saturation in the source regions of low-Ti basalts and lunar pyroclastic glasses, the extent of degassing of S in lunar magma would have to be orders of magnitude higher than currently thought, inconsistent with S isotopic and core-to-rim S diffusion profile data. The only lunar samples that could have experienced sulfide saturation are some of the more evolved A17 high-Ti basalts, if sulfides are Ni- and/or Cu rich.Sulfide saturation in the source regions of lunar melts is also inconsistent with the sulfide-silicate partitioning systematics of Ni, Co and Cu. Segregation of significant quantities of (non)-stoichiometric sulfides during fractional crystallization would result in far larger depletions of Ni, Co and Cu than observed, whereas trends in their abundances are more likely explained by olivine fractionation. The sulfide exhaustion of the lunar magma source regions agrees with previously proposed low S abundances in the lunar core and mantle, and by extension with relatively minor degassing of S during the Moon-forming event. Our results support the hypothesis that refractory chalcophile and highly siderophile element systematics of low-Ti basalts and pyroclastic glasses reflect the geochemical characteristics of their source regions, instead of indicating the presence of residual sulfides in the lunar interior.

Vulnerability of Southeast Greenland Glaciers to Warm Atlantic Water From Operation IceBridge and Ocean Melting Greenland Data.

We employ National Aeronautics and Space Administration (NASA)'s Operation IceBridge high‐resolution airborne gravity from 2016, NASA's Ocean Melting Greenland bathymetry from 2015, ice thickness from Operation IceBridge from 2010 to 2015, and BedMachine v3 to analyze 20 major southeast Greenland glaciers. The results reveal glacial fjords several hundreds of meters deeper than previously thought; the full extent of the marine‐based portions of the glaciers; deep troughs enabling warm, salty Atlantic Water (AW) to reach the glacier fronts and melt them from below; and few shallow sills that limit the access of AW. The new oceanographic and topographic data help to fully resolve the complex pattern of historical ice front positions from the 1930s to 2017: glaciers exposed to AW and resting on retrograde beds have retreated rapidly, while glaciers perched on shallow sills or standing in colder waters or with major sills in the fjords have remained stable.

Ocean-ice sheet interaction along the SE Nordic Seas margin from 35 to 15 ka BP

Sediment cores from the south-eastern Nordic Seas simultaneously archive the variability of the Fennoscandian Ice Sheet (FIS), the British-Irish Ice Sheet (BIIS) and the regional oceanic conditions. This study aims to contribute to our understanding of the marine-based section of the FIS and the BIIS between 35,000 and 15,000years BP, by using cores MD99-2283, MD99-2284 and MD99-2289, retrieved along the upper continental slope between the Faroe-Shetland Channel and the Vøring Plateau. For this, we present a revised, radiocarbon based, Bayesian modelled chronological framework and a compilation of new and published sedimentological, geochemical and micro-paleontological datasets. Our results show a possibly first Weichselian FIS/BIIS confluence at ca. 25,500years BP in the central North Sea, which buttressed the BIIS to the East, potentially leading to a northwards BIIS deflection via the Shetlands. The Norwegian Channel Ice Stream (NCIS) most likely only reached the shelf edge after 23,300±500years BP, possibly for the first time during the Weichselian. The NCIS onset directly preceded a pronounced influx of warm Atlantic water to the northern North Sea margin possibly implying forcing through ocean melt. We find a highly variable NCIS, with three ~1400 yearlong episodes of increased ice rafted debris flux interrupted by ~600 yearlong minima. When compared to other sides of the European Ice Sheet, these episodes appear to correlate well, suggesting a common forcing mechanism. In conclusion, our data supports recent suggestions that the last glacial stage of the BIIS was more extensive in the central North Sea and the confluence later than previously thought.

Mathematical Models for Global Mean Sea Level Rise

SummaryGlobal mean sea level rise (GMSLR) is one consequence of global warming. We develop simple, deterministic, and concise models for GMSLR. We consider just two factors, ocean thermal expansion and glacier melt. The resulting calculus-based models give values within the margin of error of official estimates of these two major components of GMSLR.

Seasonal and interannual changes in the sound scattering layer at deep-sea floor in the Amundsen Sea, Antarctica

Vertical migration of zooplankton is ubiquitous behavior in marine plankton community; however, seasonal and interannual behavior are little observed in the deep sea under seasonal varying sea ice. Here, the first evidence that sound scattering layers of zooplankton can support the knowledge for understanding the effect of climate change is presented, based on four-years acoustic backscattering strengths in the Amundsen Sea, Antarctica. Amundsen Sea is a biological hotspot region with the rapid oceanic melting of the ice shelf as well as the most productive (per unit area) in the Antarctic. High-temporal resolution profiles of acoustic backscattering strength collected from a bottom-moored, upward looking Acoustic Doppler Current Profiler were examined to describe the temporal variation of sound scattering layers. Our observations show that sound scattering layers exhibited clear diel, seasonal, and interannual pattern associated with solar radiation, sea ice concentration, and phytoplankton biomass. The timing of seasonal variation is also closely related with climate index for Southern Annular Mode and El Niño Southern Oscillation. This observation points out that acoustic signals from sound scattering layers could be a proxy to understand zooplankton ecosystem response to climate shifts.

Ice loss processes in the Seal Nunataks ice shelf region from satellite altimetry and imagery

ABSTRACTThe Seal Nunataks ice shelf (SNIS, ~743 km2 in 2013) is an unofficial name for a remnant area between the former Larsen A and Larsen B ice shelves off the northeastern Antarctic Peninsula. Analyses using Landsat 7 ETM+ and Terra ASTER images from 2001 to 13 and ICESat altimetry from 2003 to 09 show it has retreated and thinned following the Larsen A (1995) and Larsen B (2002) disintegrations. Despite some regional cooling and more fast ice since 2008, SNIS continues to lose ice along its margins and may be losing contact with some nunataks. Detailed analysis of data from four ICESat tracks indicates that ice shelf thinning rates range between 1.9 and 2.7 m a−1, and generally increase from west to east. An ICESat repeat track crossing the adjacent Robertson Island shows a mean elevation loss of 1.8 m a−1. Two tracks crossing the SNIS's remaining tributary, Rogosh Glacier, show sub-meter elevation losses. Comparing shelf remnant and grounded ice thinning rates implies that basal ocean melting augments SNIS thinning by ~1 m a−1, a rate that is consistent with other estimates of ocean-driven shelf thinning in the region.

Past and present dynamics of Skelton Glacier, Transantarctic Mountains

AbstractAny future changes in the volume of Antarctica’s ice sheets will depend on the dynamic response of outlet glaciers to shifts in environmental conditions. In the Transantarctic Mountains, this response is probably heavily dependent on the geometry of the system, but few studies have quantified the sensitivity of these glaciers to environmental forcings. Here we investigated the controls, along-flow sensitivity and time-dependent dynamics of Skelton Glacier. Three key outcomes were: i) present-day flow is governed primarily by surface slope, which responds to reduced valley width and large bed undulations, ii) Skelton Glacier is more susceptible to changes in atmospheric temperature than precipitation through its effect on basal sliding near the grounding line, and iii) under conditions representative of Pliocene and Quaternary climates large changes in ice thickness and velocity would have occurred in the lower reaches of the glacier. Based on these new quantitative predictions of the past and present dynamics of Skelton Glacier, we suggest that similar Transantarctic Mountain outlet glaciers could experience greater ice loss in their confined, lower reaches through increased basal sliding and ocean melt under warmer-than-present conditions. These effects are greatest where overdeepenings exist near the grounding line.

Oceanic and atmospheric forcing of Larsen C Ice-Shelf thinning

Abstract. The catastrophic collapses of Larsen A and B ice shelves on the eastern Antarctic Peninsula have caused their tributary glaciers to accelerate, contributing to sea-level rise and freshening the Antarctic Bottom Water formed nearby. The surface of Larsen C Ice Shelf (LCIS), the largest ice shelf on the peninsula, is lowering. This could be caused by unbalanced ocean melting (ice loss) or enhanced firn melting and compaction (englacial air loss). Using a novel method to analyse eight radar surveys, this study derives separate estimates of ice and air thickness changes during a 15-year period. The uncertainties are considerable, but the primary estimate is that the surveyed lowering (0.066 ± 0.017 m yr−1) is caused by both ice loss (0.28 ± 0.18 m yr−1) and firn-air loss (0.037 ± 0.026 m yr−1). The ice loss is much larger than the air loss, but both contribute approximately equally to the lowering because the ice is floating. The ice loss could be explained by high basal melting and/or ice divergence, and the air loss by low surface accumulation or high surface melting and/or compaction. The primary estimate therefore requires that at least two forcings caused the surveyed lowering. Mechanisms are discussed by which LCIS stability could be compromised in the future. The most rapid pathways to collapse are offered by the ungrounding of LCIS from Bawden Ice Rise or ice-front retreat past a "compressive arch" in strain rates. Recent evidence suggests that either mechanism could pose an imminent risk.

Potential Antarctic Ice Sheet retreat driven by hydrofracturing and ice cliff failure

Geological data indicate that global mean sea level has fluctuated on 103 to 106 yr time scales during the last ∼25 million years, at times reaching 20 m or more above modern. If correct, this implies substantial variations in the size of the East Antarctic Ice Sheet (EAIS). However, most climate and ice sheet models have not been able to simulate significant EAIS retreat from continental size, given that atmospheric CO2 levels were relatively low throughout this period. Here, we use a continental ice sheet model to show that mechanisms based on recent observations and analysis have the potential to resolve this model–data conflict. In response to atmospheric and ocean temperatures typical of past warm periods, floating ice shelves may be drastically reduced or removed completely by increased oceanic melting, and by hydrofracturing due to surface melt draining into crevasses. Ice at deep grounding lines may be weakened by hydrofracturing and reduced buttressing, and may fail structurally if stresses exceed the ice yield strength, producing rapid retreat. Incorporating these mechanisms in our ice-sheet model accelerates the expected collapse of the West Antarctic Ice Sheet to decadal time scales, and also causes retreat into major East Antarctic subglacial basins, producing ∼17 m global sea-level rise within a few thousand years. The mechanisms are highly parameterized and should be tested by further process studies. But if accurate, they offer one explanation for past sea-level high stands, and suggest that Antarctica may be more vulnerable to warm climates than in most previous studies.

Venice: Rising Water or Sinking Land?

AbstractThe recent paper by Bellfiore, Ghezzo, Tagliapietra and Umgiesser [1] gives us the opportunity to discuss the relative sea level rise in Venice (Venezia). They accept that the Venetian sea level rise will be part of a global rise caused by ocean thermal expansion and melting of icecaps caused in turn by global warming induced by anthropogenic carbon dioxide production. In reality the sea levels have been only oscillating during the last few decades. Theworldwide average sea level rise is free of any acceleration. This is also true for theVenice and Trieste tide gauges. The sea level records of Venice and Trieste are very different, and it is impossible for two places so close together to have such differences from a common rise in sea level. Venice is sinking relative to Trieste.

Basal terraces on melting ice shelves

AbstractOcean waters melt the margins of Antarctic and Greenland glaciers, and individual glaciers' responses and the integrity of their ice shelves are expected to depend on the spatial distribution of melt. The bases of the ice shelves associated with Pine Island Glacier (West Antarctica) and Petermann Glacier (Greenland) have similar geometries, including kilometer‐wide, hundreds‐of‐meter high channels oriented along and across the direction of ice flow. The channels are enhanced by, and constrain, oceanic melt. New meter‐scale observations of basal topography reveal peculiar glaciated landscapes. Channel flanks are not smooth, but are instead stepped, with hundreds‐of‐meters‐wide flat terraces separated by 5–50 m high walls. Melting is shown to be modulated by the geometry: constant across each terrace, changing from one terrace to the next, and greatly enhanced on the ~45° inclined walls. Melting is therefore fundamentally heterogeneous and likely associated with stratification in the ice‐ocean boundary layer, challenging current models of ice shelf‐ocean interactions.

Geochemical constraints on the contribution of Louisville seamount materials to magmagenesis in the Lau back-arc basin, SW Pacific

Seamounts are an integral part of element recycling in global subduction zones. The published trace element and Pb-Sr-Nd isotope data for basaltic lavas from three key segments (Central Lau Spreading Ridge (CLSR), Eastern Lau Spreading Ridge (ELSR), and Valu Fa Ridge (VFR)) of the Lau back-arc basin were compiled to evaluate the contribution of Louisville seamount materials to their magma genesis. Two geochemical transitions, separating three provinces with distinct geochemical characteristics independent of ridge segmentation, were identified based on abrupt geochemical shifts. The origin of the geochemical transitions was determined to be the result of drastic compositional changes of subduction components added into the mantle source, rather than the transition from Indian to Pacific mid-ocean ridge basalt (MORB) mantle, or due to variable mantle fertilities. The most likely explanation for the drastic shifts in subduction input is the superimposition of Louisville materials on ‘normal’ subduction components consisting predominantly of aqueous fluids liberated from the down-going altered oceanic crust and minor pelagic sediment melts. Quantitative estimation reveals that Louisville materials contributed 0–74% and 21–83% of the Th budget, respectively, to CLSR and VFR lavas, but had no definite contribution to the lavas from the ELSR, which lies farthest away from the subducted Louisville seamount chain (LSC). The spatial association of the subducted LSC with the Louisville-affected segments suggests that the Louisville signature is regionally but not locally available in the Tonga subduction zone. Besides, the preferential melting of subducted old Cretaceous LSC crust instead of the old normal Pacific oceanic crust at similar depths implies that elevated temperature across the subduction interface or seamount erosion and rupture were required to trigger melting. A wider implication of this study, thus, is that seamount subduction may promote efficiency of element recycling in subduction zones.

Marine ice sheet collapse potentially under way for the Thwaites Glacier Basin, West Antarctica.

Resting atop a deep marine basin, the West Antarctic Ice Sheet has long been considered prone to instability. Using a numerical model, we investigated the sensitivity of Thwaites Glacier to ocean melt and whether its unstable retreat is already under way. Our model reproduces observed losses when forced with ocean melt comparable to estimates. Simulated losses are moderate (<0.25 mm per year at sea level) over the 21st century but generally increase thereafter. Except possibly for the lowest-melt scenario, the simulations indicate that early-stage collapse has begun. Less certain is the time scale, with the onset of rapid (>1 mm per year of sea-level rise) collapse in the different simulations within the range of 200 to 900 years.

The effect of basal channels on oceanic ice‐shelf melting

The presence of ice‐shelf basal channels has been noted in a number of Antarctic and Greenland ice shelves, but their impact on basal melting is not fully understood. Here we use the Massachusetts Institute of Technology general circulation model to investigate the effect of ice‐shelf basal channels on oceanic melt rate for an idealized ice shelf resembling the floating tongue of Petermann Glacier in Greenland. The introduction of basal channels prevents the formation of a single geostrophically balanced boundary current; instead the flow is diverted up the right‐hand (Coriolis‐favored) side of each channel, with a return flow in the opposite direction on the left‐hand side. As the prescribed number of basal channels is increased the mean basal melt rate decreases, in agreement with previous studies. For a small number of relatively wide channels the subice flow is found to be a largely geostrophic horizontal circulation. The reduction in melt rate is then caused by an increase in the relative contribution of weakly melting channel crests and keels. For a larger number of relatively narrow channels, the subice flow changes to a vertical overturning circulation. This change in circulation results in a weaker sensitivity of melt rates to channel size. The transition between the two regimes is governed by the Rossby radius of deformation. Our results explain why basal channels play an important role in regulating basal melting, increasing the stability of ice shelves.

Pine Island glacier ice shelf melt distributed at kilometre scales

Abstract. By thinning and accelerating, West Antarctic ice streams are contributing about 10% of the observed global sea level rise. Much of this ice loss is from Pine Island Glacier, which has thinned since at least 1992, driven by changes in ocean heat transport beneath its ice shelf and retreat of the grounding line. Details of the processes driving this change, however, remain largely elusive, hampering our ability to predict the future behaviour of this and similar systems. Here, a Lagrangian methodology is developed to measure oceanic melting of such rapidly advecting ice. High-resolution satellite and airborne observations of ice surface velocity and elevation are used to quantify patterns of basal melt under the Pine Island Glacier ice shelf and the associated adjustments to ice flow. At the broad scale, melt rates of up to 100 m yr−1 occur near the grounding line, reducing to 30 m yr−1 just 20 km downstream. Between 2008 and 2011, basal melting was largely compensated by ice advection, allowing us to estimate an average loss of ice to the ocean of 87 km3 yr−1, in close agreement with 2009 oceanographically constrained estimates. At smaller scales, a network of basal channels typically 500 m to 3 km wide is sculpted by concentrated melt, with kilometre-scale anomalies reaching 50% of the broad-scale basal melt. Basal melting enlarges the channels close to the grounding line, but farther downstream melting tends to diminish them. Kilometre-scale variations in melt are a key component of the complex ice–ocean interaction beneath the ice shelf, implying that greater understanding of their effect, or very high resolution models, are required to predict the sea-level contribution of the region.

Composition and inner structure of the third layer in the oceanic crust in the subequatorial segment of the Mid-Atlantic Ridge (5°–7° N)

The paper presents data on the major-component, trace-element, and mineralogical composition of plutonic rocks, and the composition of their minerals, from the Sierra Leone region in the crest zone of the Mid-Atlantic Ridge between the Strakhov and Bogdanov fracture zones. According to their relations with seafloor structures, the rock associations are subdivided into those of rift valleys and nontransform offset zones. The troctolites and olivine gabbro composing the rift association were produced early in the fractionation course of oceanic tholeiite melt in unstationary and relatively small magmatic chambers. Most rocks beneath the nontransform offset zones crystallized during the long-lasting fractionation of the melt in large chambers hosted in serpentinized peridotites. This part consists of various cumulates, ranging from troctolites to gabbroids. Where deep tectonic detachments entered partly consolidated portions of the chambers, the melt interacted with the wall rocks. Fluid that was generated via the dehydration of serpentine and concentrated hydrophile elements, locally modified the composition of the melt and resulted in amphibole-bearing rocks. Under stress, the intercumulus melts were squeezed into tectonically weakened zones, mixed there, and also interacted with the wall rocks. These mix melts produced (with the participation of fractional crystallization) mineralized Fe-Ti gabbroids. Residual portions of the melts generated most of the diorites and plagiogranites. The high-Na diorites likely crystallized from acid melts that were derived via the partial melting of older gabbroids where aqueous fluids circulated; these fluids were generated by the deserpenitization of the host rocks in tectonized zones cutting through the chambers.