A firn aquifer in the Helheim Glacier catchment of Southeast Greenland lies directly upstream of a crevasse field. Previous measurements show that a 3.5-km long segment of the aquifer lost a large volume of water (26,000-65,000 m 2 in cross section) between spring 2012 and spring 2013, compared to annual meltwater accumulation of 6000-15,000 m 2 . The water is thought to have entered the crevasses, but whether the water reached the bed or refroze within the ice sheet is unknown. We used a thermo-visco-elastic model for crevasse propagation to calculate the depths and volumes of these water-filled crevasses. We compared our model output to data from the Airborne Topographic Mapper (ATM), which reveals the near-surface geometry of specific crevasses, and WorldView images, which capture the surface expressions of crevasses across our 1.5-km study area. We found a best fit with a shear modulus between 0.2 and 1.5 GPa within our study area. We show that surface meltwater can drive crevasses to the top surface of the firn aquifer (∼20 m depth), whereupon it receives water at rates corresponding to the water flux through the aquifer. Our model shows that crevasses receiving firn-aquifer water hydrofracture through to the bed, ∼1000 m below, in 10-40 days. Englacial refreezing of firn-aquifer water raises the average local ice temperature by ∼4 • C over a ten-year period, which enhances deformational ice motion by ∼50 m year −1 , compared to the observed surface velocity of ∼200 m year −1 . The effect of the basal water on the sliding velocity remains unknown. Were the firn aquifer not present to concentrate surface meltwater into crevasses, we find that no surface melt would reach the bed; instead, it would refreeze annually in crevasses at depths <500 m. The crevasse field downstream of the firn aquifer likely allows a large fraction of the aquifer water in our study area to reach the bed. Thus, future studies should consider the aquifer and crevasses as part of a common system. This system may uniquely affect ice-sheet dynamics by routing a large volume of water to the bed outside of the typical runoff period.

The South Pole Ice Core (SPICEcore) was drilled at least 180 km from an ice‐flow divide. Thus, the annual‐equivalent layer thicknesses in the core are affected by spatial variations in accumulation upstream in addition to temporal variations in regional accumulation. We use a new method to compare the SPICEcore accumulation record, derived by correcting measured layer thicknesses for thinning, with an accumulation record derived from new GPS and radar measurements upstream. When ice speeds are modeled as increasing by 15% since 10 ka, the upstream accumulation explains 77% of the variance in the SPICEcore‐derived accumulation (versus 22% without speedup). This result demonstrates that the ice‐flow direction and spatial pattern of accumulation were stable throughout the Holocene. The 15% speedup in turn suggests a slight (3–4%) steepening or thickening of the ice‐sheet interior and provides a new constraint on the evolution of the East Antarctic Ice Sheet following the glacial termination.

Brilliant animal colors often are produced from light interacting with intricate nano-morphologies present in biological materials such as butterfly wing scales. Surveys across widely divergent butterfly species have identified multiple mechanisms of structural color production; however, little is known about how these colors evolved. Here, we examine how closely related species and populations of Bicyclus butterflies have evolved violet structural color from brown-pigmented ancestors with UV structural color. We used artificial selection on a laboratory model butterfly, B. anynana, to evolve violet scales from UV brown scales and compared the mechanism of violet color production with that of two other Bicyclus species, Bicyclus sambulos and Bicyclus medontias, which have evolved violet/blue scales independently via natural selection. The UV reflectance peak of B. anynana brown scales shifted to violet over six generations of artificial selection (i.e., in less than 1 y) as the result of an increase in the thickness of the lower lamina in ground scales. Similar scale structures and the same mechanism for producing violet/blue structural colors were found in the other Bicyclus species. This work shows that populations harbor large amounts of standing genetic variation that can lead to rapid evolution of scales' structural color via slight modifications to the scales' physical dimensions.thin film | constructive interference | parallel evolution | photonics O rganisms produce colors in two basic ways: by synthesizing pigments that selectively absorb light of certain spectral bands so that only light outside the absorption bands is backscattered (chemical color) or by developing nanomorphologies that enhance the reflection of light of certain wavelengths by interference (physical color or structural color). Structural colors play major roles in natural and sexual selection in many species (1) and have a broad range of applications in color display, paint, cosmetics, and textile industries (2). Structural color surveys across widely divergent species have revealed a large diversity of color-producing mechanisms (3-9). However, there has been a lack of systematic study and comparison of how different colors from closely related species or within populations of a single species evolve, even though these colors can vary dramatically. By examining how these species/populations evolve different colors, it is possible to identify the minimal amount of morphological change that results in significant color variation. Furthermore, this research may serve as an inspiration for future application of similar evolutionary principles to the design of photonic devices for color tuning, light trapping, or beam steering (2,(10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20). From an evolutionary biology point of view, we are curious to examine how structural colors respond to selection pressure and whether there is sufficient standing genetic variation in natural populations to allow the rapid evolution of novel colors. Here we focus on d...

Abstract. Crosson and Dotson ice shelves are two of the most rapidly changing outlets in West Antarctica, displaying both significant thinning and grounding-line retreat in recent decades. We used remotely sensed measurements of velocity and ice geometry to investigate the processes controlling their changes in speed and grounding-line position over the past 20 years. We combined these observations with inverse modeling of the viscosity of the ice shelves to understand how weakening of the shelves affected this speedup. These ice shelves have lost mass continuously since the 1990s, and we find that this loss results from increasing melt beneath both shelves and the increasing speed of Crosson. High melt rates persisted over the period covered by our observations (1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014), with the highest rates beneath areas that ungrounded during this time. Grounding-line flux exceeded basin-wide accumulation by about a factor of 2 throughout the study period, consistent with earlier studies, resulting in significant loss of grounded as well as floating ice. The near doubling of Crosson's speed in some areas during this time is likely the result of weakening of its margins and retreat of its grounding line. This speedup contrasts with Dotson, which has maintained its speed despite increasingly high melt rates near its grounding line, likely a result of the sustained competency of the shelf. Our results indicate that changes to melt rates began before 1996 and suggest that observed increases in melt in the 2000s compounded an ongoing retreat of this system. Advection of a channel along Dotson, as well as the grounding-line position of Kohler Glacier, suggests that Dotson experienced a change in flow around the 1970s, which may be the initial cause of its continuing retreat.

Ice streams are bounded by abrupt transitions in speed called shear margins. Some shear margins are fixed by subglacial topography, but others are thought to be self‐organizing, evolving by thermal feedback to ice viscosity and basal drag which govern the stress balance of ice sheets. Resistive stresses (and properties governing shear‐margin formation) manifest nonuniquely at the surface, motivating the use of subsurface observations to constrain ice sheet models. In this study, we use ice‐penetrating radar data to evaluate three 3‐D thermomechanical models of the Northeast Greenland Ice Stream, focusing on model reproductions of ice temperature (a primary control on viscosity) and subsurface velocity. Data/model agreement indicates elevated temperatures in the Northeast Greenland Ice Stream margins, with depth‐averaged temperatures between 2 °C and 6 °C warmer in the southeast margin compared to ice in streaming flow, driven by vertical heat transport rather than shear heating. This work highlights complexity in ice divergence across stagnant/streaming transitions.

Abstract. The area near Dome C, East Antarctica, is thought to be one of the most promising targets for recovering a continuous ice-core record spanning more than a million years. The European Beyond EPICA consortium has selected Little Dome C (LDC), an area ∼ 35 km southeast of Concordia Station, to attempt to recover such a record. Here, we present the results of the final ice-penetrating radar survey used to refine the exact drill site. These data were acquired during the 2019–2020 austral summer using a new, multi-channel high-resolution very high frequency (VHF) radar operating in the frequency range of 170–230 MHz. This new instrument is able to detect reflectors in the near-basal region, where previous surveys were largely unable to detect horizons. The radar stratigraphy is used to transfer the timescale of the EPICA Dome C ice core (EDC) to the area of Little Dome C, using radar isochrones dating back past 600 ka. We use these data to derive the expected depth–age relationship through the ice column at the now-chosen drill site, termed BELDC (Beyond EPICA LDC). These new data indicate that the ice at BELDC is considerably older than that at EDC at the same depth and that there is about 375 m of ice older than 600 kyr at BELDC. Stratigraphy is well preserved to 2565 m, ∼ 93 % of the ice thickness, below which there is a basal unit with unknown properties. An ice-flow model tuned to the isochrones suggests ages likely reach 1.5 Myr near 2500 m, ∼ 65 m above the basal unit and ∼ 265 m above the bed, with sufficient resolution (19 ± 2 kyr m−1) to resolve 41 kyr glacial cycles.

Bulk directional enhancement factors are determined for axisymmetric (girdle and single-maximum) orientation fabrics using a transversely isotropic grain rheology with an orientation-dependent non-linear grain fluidity. Compared to grain fluidities that are simplified as orientation independent, we find that bulk strain-rate enhancements for intermediate-to-strong axisymmetric fabrics can be up to a factor of ten larger, assuming stress homogenization over the polycrystal scale. Our work thus extends previous results based on simple basal slip (Schmid) grain rheologies to the transversely isotropic rheology, which has implications for large-scale anisotropic ice-flow modelling that relies on a transversely isotropic grain rheology. In order to derive bulk enhancement factors for arbitrary evolving fabrics, we expand the c-axis distribution in terms of a spherical harmonic series, which allows the rheology-required structure tensors through order eight to easily be calculated and provides an alternative to current structure-tensor-based modelling.

The crystal structure of glaciers and ice sheets records past flow and affects present-day ice deformation. Much work on ice-crystal orientation fabric has focused on the fabric's effect upon deformation (e.g.