Explaining the presence of perennial liquid water bodies in the firn of the Greenland Ice Sheet

pre-printRecent observations have shown that the firn layer on the Greenland Ice Sheet features subsurface bodies of liquid water at the end of the winter season. Using a model with basic firn hydrology, thermodynamics, and compaction in one dimension, we find that a combination of moderate to strong surface melt and a high annual accumulation rate is required to form such a perennial firn aquifer. The high accumulation rate ensures that there is pore space available to store water at a depth where it is protected from the winter cold. Low-accumulation sites cannot provide sufficiently deep pore space to store liquid water. However, for even higher accumulation rates, the total cold content of the winter accumulation becomes sufficient to refreeze the total amount of liquid water. As a consequence, wintertime or springtime observations of subsurface liquid water in these specific accumulation conditions cannot distinguish between a truly perennial firn aquifer and water layers that will ultimately refreeze completely

Systems Analysis of complex glaciological processes and application to calving of Amery Ice Shelf, East Antarctica

Calving is a complex process subject to several cooperating atmospheric, oceanographic and glaciological forcings that vary in space and time, and whose relative effects are challenging to separate. Statistical ‘Systems Analysis’ is commonly used in engineering and economics to extricate complex ‘force–response’ relationships. Here we apply Systems Analysis to the Amery rift system, East Antarctica. We develop a scalable ‘System Model’ driven by a coarsely-sampled dataset characteristic of glaciological observations in remote locations, and validate it using rift lengths observed in 2000–06 and 2012. In this initial demonstration, we forecast a detachment date of ∼2019 ± 5 years for the large tabular iceberg colloquially known as the ‘Loose Tooth’, for which relative humidity surprisingly emerges as the best statistical predictor. RACMO2 climate modelling reveals that relative humidity correlates best with surface albedo and snowmelt, both of which are intimately linked to firn compaction and ice shelf temperature and flow. We postulate that relative humidity can therefore serve as a proxy for internal stress, a known key control of ‘Loose Tooth’ calving. Although no physical causality is implied in Systems Analysis, postulates such as this can aid in setting priorities in studies of complex glaciological processes

Spatial Response of Greenland's Firn Layer to NAO Variability

Firn on the Greenland Ice Sheet (GrIS) buffers meltwater, and has a variable thickness, complicating observations of volume change to mass change. In this study, we use a firn model (IMAU-FDM v1.2G) forced by a regional climate model (RACMO2.3p2) to investigate how the GrIS firn layer thickness and pore space have evolved since 1958 in response to variability in the large-scale atmospheric circulation. On interannual timescales, the firn layer thickness and pore space show a spatially heterogeneous response to variability in the North Atlantic Oscillation (NAO). Notably, a stronger NAO following the record warm summer of 2012 led the firn layer in the south and east of the ice sheet to regain thickness and pore space after a period of thinning and reduced pore space. In the southwest, a decrease in melt dominated after 2012, whereas in the east, the main driver was an increase in snow accumulation. At the same time, the firn in the northwestern ice sheet continued to lose pore space. The NAO also varies on intra-annual timescales, being typically stronger in winter than in summer. This impacts the amplitude of the seasonal cycle in GrIS firn thickness and pore space. In the wet southeastern GrIS, most of the snow accumulates during the winter, when melting and densification are relatively weak, leading to a large seasonal cycle in thickness and pore space. The opposite occurs in other regions, where snowfall peaks in summer or autumn. This dampens the seasonal amplitude of firn thickness and pore space

The influence of föhn winds on annual and seasonal surface melt on the Larsen C Ice Shelf, Antarctica

Abstract. Warm, dry föhn winds are observed over the Larsen C Ice Shelf year-round and are thought to contribute to the continuing weakening and collapse of ice shelves on the eastern Antarctic Peninsula (AP). We use a surface energy balance (SEB) model, driven by observations from two locations on the Larsen C Ice Shelf and one on the remnants of Larsen B, in combination with output from the Antarctic Mesoscale Prediction System (AMPS), to investigate the year-round impact of föhn winds on the SEB and melt from 2009 to 2012. Föhn winds have an impact on the individual components of the surface energy balance in all seasons and lead to an increase in surface melt in spring, summer and autumn up to 100 km away from the foot of the AP. When föhn winds occur in spring they increase surface melt, extend the melt season and increase the number of melt days within a year. Whilst AMPS is able to simulate the percentage of melt days associated with föhn with high skill, it overestimates the total amount of melting during föhn events and non-föhn events. This study extends previous attempts to quantify the impact of föhn on the Larsen C Ice Shelf by including a 4-year study period and a wider area of interest and provides evidence for föhn-related melting on both the Larsen C and Larsen B ice shelves. </jats:p

Trends in Antarctic Peninsula surface melting conditions from observations and regional climate modeling

Multidecadal meteorological station records and microwave backscatter time-series from the SeaWinds scatterometer onboard QuikSCAT (QSCAT) were used to calculate temporal and spatial trends in surface melting conditions on the Antarctic Peninsula (AP). Four of six long-term station records showed strongly positive and statistically significant trends in duration of melting conditions, including a 95% increase in the average annual positive degree day sum (PDD) at Faraday/Vernadsky, since 1948. A validated, threshold-based melt detection method was employed to derive detailed melt season onset, extent, and duration climatologies on the AP from enhanced resolution QSCAT data during 1999–2009. Austral summer melt on the AP was linked to regional- and synoptic-scale atmospheric variability by respectively correlating melt season onset and extent with November near-surface air temperatures and the October–January averaged index of the Southern Hemisphere Annular Mode (SAM). The spatial pattern, magnitude, and interannual variability of AP melt from observations was closely reproduced by simulations of the regional model RACMO2. Local discrepancies between observations and model simulations were likely a result of the QSCAT response to, and RACMO2 treatment of, ponded surface water, and the relatively crude representation of coastal climate in the 27 km RACMO2 grid

Centuries of intense surface melt on Larsen C Ice Shelf

Following a southward progression of ice-shelf disintegration along the Antarctic Peninsula, Larsen C Ice Shelf is the focus of ongoing investigation regarding its future stability. The ice shelf is known to be experience surface melt, and commonly features surface meltwater ponds. Here, we use a flowline model and a firn density model to date and interpret observations of melt-affected ice layers found within five 90?m boreholes distributed across the ice shelf. We find that units of ice within the boreholes, which have densities exceeding those expected under normal compaction metamorphism, correspond to two climatic warm periods within the last 300 years on the Antarctic Peninsula. The more recent warm period, from the 1960s onwards, has generated distinct sections of dense ice in two boreholes in Cabinet Inlet, close to the Antarctic Peninsula mountains ? a region currently affected by f?hn winds. Previous work has classified these layers as refrozen pond ice, requiring large quantities of mobile liquid water to form. Our flowline model shows that, whilst preconditioning of the ice began in the late 1960s, it was probably not until the early 1990s that twentieth-century ponding began. The earlier warm period occurred during the 18th century and resulted in two additional sections of anomalously dense ice deep within the boreholes. The first, in one of the Cabinet Inlet boreholes, consists of ice characteristic of refrozen ponds and must have formed in an area currently featuring ponding. The second, in a mid-shelf borehole, formed at the same time in an area which now experiences significant annual melt. Further south on the shelf, the boreholes sample ice that is of an equivalent age but which does not exhibit the same degree of melt influence. This west?east and north?south gradient in past melt distribution resembles current spatial patterns of surface melt intensity. Using flowlines to trace the advection and submergence of continental ice identified in boreholes, we demonstrate that, even by the time the ice reaches the calving front, only the upper 40 to 50?% of the shelf is composed of meteoric ice accumulated on the shelf. This vertical composition implies that basal crevasses must be confined within continental and/or basally accreted ice, and therefore will be unaffected by current climate-induced firn compactionauthorsversio

Ice and firn heterogeneity within Larsen C Ice Shelf from borehole optical televiewing

Research was funded by the UK Natural Environmental Research Council grants NE/L006707/1 and NE/L005409/1 and a HEFCW/Aberystwyth University Capital Equipment Grant to B.H. Data will be available via the project website (www.projectmidas.org) and the UK Polar Data Centre (https://www.bas.ac.uk/data/uk-pdc/) from mid-2017.We use borehole optical televiewing (OPTV) to explore the internal structure of Larsen C Ice Shelf (LCIS). We report a suite of five ~90 m long OPTV logs, recording a light-emitting diode-illuminated, geometrically correct image of the borehole wall, from the northern and central sectors of LCIS collected during austral spring 2014 and 2015. We use a thresholding-based technique to estimate the refrozen ice content of the ice column and exploit a recently calibrated density-luminosity relationship to reveal its structure. All sites are dense and strongly influenced by surface melt, with frequent refrozen ice layers and mean densities, between the depths of 1.87 and 90 m, ranging from 862 to 894 kg m−3. We define four distinct units that comprise LCIS and relate these to ice provenance, dynamic history, and past melt events. These units are in situ meteoric ice with infiltration ice (U1), meteoric ice which has undergone enhanced densification (U2), thick refrozen ice (U3), and advected continental ice (U4). We show that the OPTV-derived pattern of firn air content is consistent with previous estimates, but that a significant proportion of firn air is contained within U4, which we interpret to have been deposited inland of the grounding line. The structure of LCIS is strongly influenced by the E-W gradient in föhn-driven melting, with sites close to the Antarctic Peninsula being predominantly composed of refrozen ice. Melting is also substantial toward the ice shelf center with >40% of the overall imaged ice column being composed of refrozen ice.Publisher PDFPeer reviewe

A new albedo parameterization for use in climate models over the Antarctic ice sheet

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/95255/1/jgrd16889.pd

Characteristics of the 1979–2020 Antarctic firn layer simulated with IMAU-FDM v1.2A

Firn simulations are essential for understanding Antarctic ice sheet mass change, as they enable us to convert satellite altimetry observed volume changes to mass changes and column thickness to ice thickness and to quantify the meltwater buffering capacity of firn. Here, we present and evaluate a simulation of the contemporary Antarctic firn layer using the updated semi-empirical IMAU Firn Densification Model (IMAU-FDM) for the period 1979–2020. We have improved previous fresh-snow density and firn compaction parameterizations and used updated atmospheric forcing. In addition, the model has been calibrated and evaluated using 112 firn core density observations across the ice sheet. We found that 62 % of the seasonal and 67 % of the decadal surface height variability are due to variations in firn air content rather than firn mass. Comparison of simulated surface elevation change with a previously published multi-mission altimetry product for the period 2003–2015 shows that performance of the updated model has improved, notably in Dronning Maud Land and Wilkes Land. However, a substantial trend difference (&gt;10 cm yr−1) remains in the Antarctic Peninsula and Ellsworth Land, mainly caused by uncertainties in the spin-up forcing. By estimating previous climatic conditions from ice core data, these trend differences can be reduced by 38 %.</p

Flow-line model code for accumulation of ice along velocity-based trajectories

The flow-line model was designed to enable estimation of the age and surface origin for various ice bodies identified within hot-water drilled boreholes on Larsen C Ice Shelf. Surface fluxes are accumulated, converted to thicknesses, and advected down flow from a fixed number of selected points. The model requires input datasets of surface mass balance, surface velocity, vertical strain rates, ice-shelf thickness, and a vertical density profile. This model is part of a larger project. Input datasets such as density profiles and trajectory vectors are available separately. Resolution is dependent on the input datasets. Funding was provided by the NERC grant NE/L005409/1