Search

Ice core archives are well suited for reconstructing rapid past climate changes at high latitudes. Despite this, few records currently exist from coastal Greenlandic ice caps due to their remote nature, limiting our long-term understanding of past maritime and coastal climate variability across this rapidly changing Arctic region. Here, we reconstruct regionally representative glacier surface mass balance and climate variability over the last two thousand years (~169–2015 ce) using an ice core collected from the Nuussuaq Peninsula, west Greenland. We find indications of abrupt regional hydroclimate shifts, including an up to 20% decrease in average snow accumulation during the transition from the Medieval Warm Period (950–1250 ce) to Little Ice Age (1450–1850 ce), followed by a subsequent >40% accumulation increase from the early 18th to late 20th centuries ce. These coastal changes are substantially larger than those previously reported from interior Greenland records. Moreover, we show that the strong relationship observed today between Arctic temperature rise and coastal ice cap decay contrasts with that of the last millennium, during which periods of warming led to snowfall-driven glacial growth. Taken together with modern observations, the ice core evidence could indicate a recent reversal in the response of west Greenland ice caps to climate change. Coastal west Greenland ice caps fluctuated strongly compared to the interior in response to rapid Common Era changes in snow accumulation, according to modelling of proxy records developed from a Nuussuaq Peninsula ice core covering the last 2,000 years.

Greenland supraglacial lakes sometimes drain to the ice-sheet base only hours to days apart in time, leading to the hypothesis that stress transmission following one drainage may be sufficient to induce hydro-fracture-driven drainages of other lakes. However, available observations characterizing the time- and length-scales of drainage-induced stress perturbations are insufficient to evaluate this hypothesis, hindering our understanding of whether lower-elevation lake drainages could generate stress changes that initiate hydro-fracture beneath lakes in inland ice-sheet regions. Here, we invert ice-sheet surface displacement observations from a dense Global Positioning System (GPS) array deployed around three supraglacial lakes in the Greenland Ice Sheet ablation zone to estimate the space-time history of slip and opening along pre-defined planes within an elastic half-space that best reproduces GPS observations during two rapid drainage events. We then use these slip and opening estimates to forward model ice-sheet surface stresses across all three neighboring lake basins and investigate stress transmission between the basins. We find that drainage of a central lake places neighboring basins in either tensional or compressional surface stress relative to their individual hydro-fracture scarp orientations, which, respectively, promotes or inhibits hydro-fracture initiation beneath those lakes. Surface-stress-change direction is asymmetric across the lake-draining fracture of the central lake and highly dependent on where surface uplift occurs. On elastic timescales, lake-drainage-induced stresses within the extent of the GPS array in the inland direction remain low within ~3 km of the central lake. This short length-scale for stress coupling in the inland direction is consistent with idealized lake-drainage scenarios for different lake volumes and ice-sheet thicknesses. Our observations and idealized-model results support the stress-transmission hypothesis for inducing hydro-fracture-driven drainage of lakes located within the region of uplift produced by the initial drainage, but refute this hypothesis for distal lakes not located down the subglacial hydraulic catchment from the initial drainage.

Reconstruction of the North Atlantic jet stream (NAJ) presents a critical, albeit largely unconstrained, paleoclimatic target. Models suggest northward migration and changing variance of the NAJ under 21st-century warming scenarios, but assessing the significance of such projections is hindered by a lack of long-term observations. Here, we incorporate insights from an ensemble of last-millennium water isotope–enabled climate model simulations and a wide array of mean annual water isotope ([Formula: see text] O) and annually accumulated snowfall records from Greenland ice cores to reconstruct North Atlantic zonal-mean zonal winds back to the 8th century CE. Using this reconstruction we provide preobservational constraints on both annual mean NAJ position and intensity to show that late 20th- and early 21st-century NAJ variations were likely not unique relative to natural variability. Rather, insights from our 1,250 year reconstruction highlight the overwhelming role of natural variability in thus far masking the response of midlatitude atmospheric dynamics to anthropogenic forcing, consistent with recent large-ensemble transient modeling experiments. This masking is not projected to persist under high greenhouse gas emissions scenarios, however, with model projected annual mean NAJ position emerging as distinct from the range of reconstructed natural variability by as early as 2060 CE.

Surface meltwater reaching the base of the Greenland Ice Sheet transits through drainage networks, modulating the flow of the ice sheet. Dye and gas-tracing studies conducted in the western margin sector of the ice sheet have directly observed drainage efficiency to evolve seasonally along the drainage pathway. However, the local evolution of drainage systems further inland, where ice thicknesses exceed 1000 m, remains largely unknown. Here, we infer drainage system transmissivity based on surface uplift relaxation following rapid lake drainage events. Combining field observations of five lake drainage events with a mathematical model and laboratory experiments, we show that the surface uplift decreases exponentially with time, as the water in the blister formed beneath the drained lake permeates through the subglacial drainage system. This deflation obeys a universal relaxation law with a timescale that reveals hydraulic transmissivity and indicates a two-order-of-magnitude increase in subglacial transmissivity (from 0.8 ± 0.3 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rm{m}}{{\rm{m}}}^{3}$$\end{document}mm3 to 215 ± 90.2 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rm{m}}{{\rm{m}}}^{3}$$\end{document}mm3) as the melt season progresses, suggesting significant changes in basal hydrology beneath the lakes driven by seasonal meltwater input., Hydraulic transmissivity under the 1km-thick Greenland Ice Sheet was inferred by ice-sheet uplift relaxation after rapid lake drainage events. A two-order-of-magnitude increase in hydraulic transmissivity was found throughout the melt season.

We acknowledge support from the Woods Hole Oceanographic Institution Ocean and Climate Change Institute Arctic Research Initiative, and NSF OPP-1418256 and OPP-1743693, to Fiamma Straneo and Sarah B. Das. Till J. W. Wagner was further supported by NSF OPP award 1744835. Geospatial support for this work was provided by the Polar Geospatial Center under NSF OPP awards 1043681 and 1559691. DEMs provided by the Polar Geospatial Center under NSF OPP awards 1043681, 1559691, and 1542736. Donald A. Slater acknowledges the support of Scottish Alliance for Geoscience, Environment and Society early-career research exchange funding. The frontal flux balance of a medium-sized tidewater glacier in western Greenland in the summer is assessed by quantifying the individual components (ice flux, retreat, calving, and submarine melting) through a combination of data and models. Ice flux and retreat are obtained from satellite data. Submarine melting is derived using a high-resolution ocean model informed by near-ice observations, and calving is estimated using a record of calving events along the ice front. All terms exhibit large spatial variability along the ∼ 5 km wide ice front. It is found that submarine melting accounts for much of the frontal ablation in small regions where two subglacial discharge plumes emerge at the ice front. Away from the subglacial plumes, the estimated melting accounts for a small fraction of frontal ablation. Glacier-wide, these estimates suggest that mass loss is largely controlled by calving. This result, however, is at odds with the limited presence of icebergs at this calving front-suggesting that melt rates in regions outside of the subglacial plumes may be underestimated. Finally, we argue that localized melt incisions into the glacier front can be significant drivers of calving. Our results suggest a complex interplay of melting and calving marked by high spatial variability along the glacier front. Publisher PDF

A dense network of GPS observations shows that rapid lake drainage events on the western Greenland Ice Sheet are preceded by period of ice-sheet uplift and/or enhanced basal slip. Surface lakes are common during the warm season on the Greenland Ice Sheet, especially in middle to lower elevations. These lakes can drain rapidly, often preceded by a hydro-fracturing event. The mechanism behind the hydro-fracturing, and thus the connection between the lakes and the broader Greenland hydrological system, has been unclear. Now, Laura Stevens and colleagues use a dense network of GPS observations to show that rapid lake drainage events at their study site in western Greenland were preceded by uplift events, likely from water entering the basal system from neighboring moulins (vertical conduits connecting the surface and bedrock). The uplift created sufficient stress in the ice sheet to open the hydro-fractures, thus allowing the rapid lake drainage. Although speculative, the authors suggest that the comparative lack of crevassing in upper elevations might limit the geographical spread of the uplift-fracturing mechanism. Water-driven fracture propagation beneath supraglacial lakes rapidly transports large volumes of surface meltwater to the base of the Greenland Ice Sheet1. These drainage events drive transient ice-sheet acceleration1,2,3 and establish conduits for additional surface-to-bed meltwater transport for the remainder of the melt season1,4,5,6. Although it is well established that cracks must remain water-filled to propagate to the bed7,8,9, the precise mechanisms that initiate hydro-fracture events beneath lakes are unknown. Here we show that, for a lake on the western Greenland Ice Sheet, drainage events are preceded by a 6–12 hour period of ice-sheet uplift and/or enhanced basal slip. Our observations from a dense Global Positioning System (GPS) network allow us to determine the distribution of meltwater at the ice-sheet bed before, during, and after three rapid drainages in 2011–2013, each of which generates tensile stresses that promote hydro-fracture beneath the lake. We hypothesize that these precursors are associated with the introduction of meltwater to the bed through neighbouring moulin systems (vertical conduits connecting the surface and base of the ice sheet). Our results imply that as lakes form in less crevassed, interior regions of the ice sheet10,11,12,13,14, where water at the bed is currently less pervasive5,14,15,16, the creation of new surface-to-bed conduits caused by lake-draining hydro-fractures may be limited.

We present the first published uranium-series measurements from modern Greenland Ice Sheet (GrIS) runoff and proximal seawater, and investigate the influence of glacial melt on global seawater δ234U over glacial-interglacial (g-ig) timescales. Climate reconstructions based on closed-system uranium-thorium (U/Th) dating of fossil corals assume U chemistry of seawater has remained stable over time despite notable fluctuations in major elemental compositions, concentrations, and isotopic compositions of global seawater on g-ig timescales. Deglacial processes increase weathering, significantly increasing U-series concentrations and changing the δ234U of glacial meltwater. Analyses of glacial discharge from GrIS outlet glaciers indicate that meltwater runoff has elevated U concentrations and differing 222Rn concentrations and δ234U compositions, likely due to variations in subglacial residence time. Locations with high δ234U have the potential to increase proximal seawater δ234U. To better understand the impact of bulk glacial melt on global seawater δ234U over time, we use a simple box model to scale these processes to periods of extreme deglaciation. We account for U fluxes from the GrIS, Antarctica, and large Northern Hemisphere Continental Ice Sheets, and assess sensitivity by varying melt volumes, duration and U flux input rates based on modern subglacial water U concentrations and compositions. All scenarios support the hypothesis that global seawater δ234U has varied by more than 1‰ through time as a function of predictable perturbations in continental U fluxes during g-ig periods.

Meltwater and sediment-laden plumes at tidewater glaciers, resulting from the localized subglacial discharge of surface melt, influence submarine melting of the glacier and the delivery of nutrients to the fjord's surface waters. It is usually assumed that increased subglacial discharge will promote the surfacing of these plumes. Here, at a western Greenland tidewater glacier, we investigate the counterintuitive observation of a non-surfacing plume in July 2012 (a year of record surface melting) compared to the surfacing of the plume in July 2013 (an average melt year). We combine oceanographic observations, subglacial discharge estimates and an idealized plume model to explain the observed plumes' behavior and evaluate the relative impact of fjord stratification and subglacial discharge on plume dynamics. We find that increased fjord stratification prevented the plume from surfacing in 2012, show that the fjord was more stratified in 2012 due to increased freshwater content and speculate that this arose from an accumulation of ice sheet surface meltwater in the fjord in this record melt year. By developing theoretical scalings, we show that fjord stratification in general exerts a dominant control on plume vertical extent (and thus surface expression), so that studies using plume surface expression as a means of diagnosing variability in glacial processes should account for possible changes in stratification. We introduce the idea that, despite projections of increased surface melting over Greenland, the appearance of plumes at the fjord surface could in the future become less common if the increased freshwater acts to stratify fjords around the Greenland ice sheet. We discuss the implications of our findings for nutrient fluxes, trapping of atmospheric CO2 and the properties of water exported from Greenland's fjords.

Marine phytoplankton have a crucial role in the modulation of marine-based food webs1, fishery yields2 and the global drawdown of atmospheric carbon dioxide3. However, owing to sparse measurements before satellite monitoring in the twenty-first century, the long-term response of planktonic stocks to climate forcing is unknown. Here, using a continuous, multi-century record of subarctic Atlantic marine productivity, we show that a marked 10 ± 7% decline in net primary productivity has occurred across this highly productive ocean basin over the past two centuries. We support this conclusion by the application of a marine-productivity proxy, established using the signal of the planktonic-derived aerosol methanesulfonic acid, which is commonly identified across an array of Greenlandic ice cores. Using contemporaneous satellite-era observations, we demonstrate the use of this signal as a robust and high-resolution proxy for past variations in spatially integrated marine productivity. We show that the initiation of declining subarctic Atlantic productivity broadly coincides with the onset of Arctic surface warming4, and that productivity strongly covaries with regional sea-surface temperatures and basin-wide gyre circulation strength over recent decades. Taken together, our results suggest that the decline in industrial-era productivity may be evidence of the predicted5 collapse of northern Atlantic planktonic stocks in response to a weakened Atlantic Meridional Overturning Circulation6–8. Continued weakening of this Atlantic Meridional Overturning Circulation, as projected for the twenty-first century9,10, may therefore result in further productivity declines across this globally relevant region. A continuous, multi-century record of subarctic Atlantic marine productivity shows that a marked decline in net primary productivity has occurred across the subarctic Atlantic basin over the past two centuries.

Support was provided by the National Science Foundation (NSF) through PLR-1418256 and PLR-1744835, and through Woods Hole Oceanographic Institution (WHOI) Ocean and Climate Change Institute (OCCI) and the Clark Foundation. This work was also supported by a UK Natural Environmental Research Council (NERC) PhD studentship (NE/L501566/1) and Scottish Alliance for Geoscience, Environment & Society (SAGES) early career research exchange funding to D. A. S. Recent 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. Publisher PDF

Continuous seismic observations across the Ross Ice Shelf reveal ubiquitous ambient resonances at frequencies >5 Hz. These firn‐trapped surface wave signals arise through wind and snow bedform interactions coupled with very low velocity structures. Progressive and long‐term spectral changes are associated with surface snow redistribution by wind and with a January 2016 regional melt event. Modeling demonstrates high spectral sensitivity to near‐surface (top several meters) elastic parameters. We propose that spectral peak changes arise from surface snow redistribution in wind events and to velocity drops reflecting snow lattice weakening near 0°C for the melt event. Percolation‐related refrozen layers and layer thinning may also contribute to long‐term spectral changes after the melt event. Single‐station observations are inverted for elastic structure for multiple stations across the ice shelf. High‐frequency ambient noise seismology presents opportunities for continuous assessment of near‐surface ice shelf or other firn environments.

Marine phytoplankton have a crucial role in the modulation of marine-based food webs

We investigate the frontal mass budget of a medium-sized tidewater glacier in western Greenland. This is done by comparing the seasonal retreat of the glacier to ice advection and ablation along the front. Frontal ablation is partitioned into calving and submarine melting, both of which are estimated from in situ observations. We observe large spatial variability in all mass budget terms along the glacier front. In particular, we find that the ablation of the glacier front is characterized by two main regimes: melting dominated versus calving dominated. While melting-dominated segments appear to be associated with subglacial discharge plumes, calving-dominated regions occur outside such plumes. The melting-dominated segments are rather localized, and the majority of ablation is estimated to occur in the form of calving. However, we stress the large uncertainty in melt rate estimates and consider the possibility that current parameterizations substantially underestimate melting. Finally, we argue that localized melt incisions into the glacier front can be significant drivers of calving. Our results suggest a complex interplay of melting and calving marked by high spatial variability along the glacier front. Understanding the impact of such local variability on larger scale ice dynamics may help guide future mass balance projections for tidewater glaciers.

We use a numerical subglacial hydrology model and remotely sensed observations of Greenland Ice Sheet surface motion to test whether the inverse relationship between effective pressure and regional melt season surface speeds observed at individual sites holds on a regional scale. The model is forced with daily surface runoff estimates for 2009 and 2010 across an ~8,000‐km2 region on the western margin. The overall subglacial drainage system morphology develops similarly in both years, with subglacial channel networks growing inland from the ice sheet margin and robust subglacial pathways forming over bedrock ridges. Modeled effective pressures are compared to contemporaneous regional surface speeds derived from TerraSAR‐X imagery to investigate spatial relationships. Our results show an inverse spatial relationship between effective pressure and ice speed in the mid‐melt season, when surface speeds are elevated, indicating that effective pressure is the dominant control on surface velocities in the mid‐melt season. By contrast, in the early and late melt seasons, when surface speeds are slower, effective pressure and surface speed have a positive relationship. Our results suggest that outside of the mid‐melt season, the influence of effective pressures on sliding speeds may be secondary to the influence of driving stress and spatially variable bed roughness.

This study represents the first comprehensive noble gas study in glacial meltwater from the Greenland Ice Sheet. It shows that most samples are in disequilibrium with surface collection conditions. A preliminary Ne and Xe analysis suggests that about half of the samples equilibrated at a temperature of ~0°C and altitudes between 1000 m and 2000 m, with a few samples pointing to lower equilibration altitudes and temperatures between 2°C and 5°C. Two samples suggest an origin as melted ice and complete lack of equilibration with surface conditions. A helium component analysis suggests that this glacial meltwater was isolated from the atmosphere prior to the 1950s, with most samples yielding residence times ≤ 420 years. Most samples represent a mixture between a dominant atmospheric component originating as precipitation and basal meltwater or groundwater, which has accumulated crustal 4He over time.

Ice shelves modulate Antarctica’s contributions to sea-level rise. Regional-climate-model simulations and observations suggest historical ice melt intensification before collapse of Antarctic peninsula shelves, and project future melt evolution.

We analyzed geophone and GPS measurements collected within the ablation zone of the western Greenland Ice Sheet during a ~35 day period of the 2011 melt season to study changes in ice deformation before, during, and after a supraglacial lake drainage event. During rapid lake drainage, ice flow speeds increased to ~400% of winter values, and icequake activity peaked. At times >7 days after drainage, this seismicity developed variability over both diurnal and longer periods (~10 days), while coincident ice speeds fell to ~150% of winter values and showed nightly peaks in spatial variability. Approximately 95% of all detected seismicity in the lake basin and its immediate vicinity was triggered by fracture propagation within near-surface ice (

We produced four independent temperature time series derived from different sensors for the Siple Dome region of West Antarctica to investigate seasonal to interannual temperature variability over the last 20 years. We use data from automatic weather station air-temperature records (1997–99), European Centre for Medium-range Weather Forecasts surface temperature from the 15 year re-analyses (ERA-15, 1979–93), and emissivity-corrected brightness temperatures from the Special Sensor Microwave/Imager (1987–99) and the Scanning Multichannel Microwave Radiometer (1978–87). Each technique has limitations and errors, but all respond to temperature, and all agree in the large patterns of temperature variability over time. Our results show that there is high seasonal to interannual variability in both mean temperature and variance in the Siple Dome region during the study period. In particular, fluctuations in seasonal to interannual temperature variance occur on an approximately 5 year cycle and correlate with variations in the Southern Oscillation Index.

Insoluble aerosol particles trapped in glacial ice provide insight into past climates, but analysis requires information on climatically-relevant particle properties, such as size, abundance, and internal mixing. We present a new analytical method using a time-of-flight single particle mass spectrometer (SPMS) to determine the composition and size of insoluble particles in glacial ice over an aerodynamic size range of ~ 0.2–3.0 µm diameter. Using samples from two Greenland ice cores, we developed a procedure to nebulize insoluble particles suspended in melted ice, evaporate condensed liquid from those particles, and transport them to the SPMS for analysis. We further determined size-dependent extraction and instrument transmission efficiencies to investigate particle-class specific mass concentrations. We find SPMS can be used to provide constraints on the aerodynamic size, composition, and relative abundance of most insoluble particulate classes in ice core samples. We describe the importance of post-aqueous processing to particles, a process which occurs due to nebulization of aerosols from an aqueous suspension of originally soluble and insoluble aerosol components. This study represents an initial attempt to use SPMS as an emerging technique for the study of insoluble particulates in ice cores.

In this study, we present the first measurements of iron (Fe) stable isotopic composition (δ 56 Fe) of subglacial streams draining the Greenland Ice Sheet (GIS). We measure the δ 56 Fe values [(δ 56 Fe, ‰ = ( 56 Fe/ 54 Fe) sample /( 56 Fe/ 54 Fe) standard − 1) × 10 3 ] of both dissolved and suspended sediment Fe in subglacial outflows from five distinct land-terminating glaciers. Suspended sediments have δ 56 Fe values that lie within the crustal array (δ 56 Fe ∼ 0‰). In contrast, the δ 56 Fe values of dissolved Fe in subglacial outflows are consistently less than 0‰, reaching a minimum of −2.1‰ in the outflow from the Russell Glacier. The δ 56 Fe values of dissolved Fe vary geographically and on daily time scales. Major element chemistry and mineral saturation state modeling suggest that incongruent silicate weathering and sulfide oxidation are the likely drivers of subglacial stream Fe chemistry, and that the extent of chemical weathering influences the δ 56 Fe of dissolved Fe. The largest difference in δ 56 Fe between dissolved and suspended load is −2.1‰, and occurs in the subglacial system from the Russell glacier (southwest GIS). Major element chemistry indicates this outflow to be the least chemically weathered, while more mature subglacial systems (i.e., that exhibit greater extents of subglacial weathering) have dissolved loads with δ 56 Fe that are indistinguishable from suspended sediments (Δ 56 Fe suspended-dissolved ∼ 0‰). Ultimately, the dissolved Fe generated in some subglacial systems from the GIS is a previously unrecognized source of isotopically light Fe into the hydrosphere. The data illustrate that the dissolved Fe supplied by subglacial weathering can have variable δ 56 Fe values depending on the degree of chemical weathering. Thus, Fe isotopes have potential as a proxy for subglacial chemical weathering intensity or mode. Finally, based on our regional Fe concentration measurements from each glacial outflow, we estimate a flux weighted continental scale dissolved iron export of 2.1 Gg Fe yr −1 to the coastal ocean, which is within the range of previous estimates.

Methanesulfonic acid (MSA; CH3SO3H) in polar ice is a unique proxy of marine primary productivity, synoptic atmospheric transport, and regional sea-ice behavior. However, MSA can be mobile within the firn and ice matrix, a post-depositional process that is well known but poorly understood and documented, leading to uncertainties in the integrity of the MSA paleoclimatic signal. Here, we use a compilation of 22 ice core MSA records from Greenland and Antarctica and a model of soluble impurity transport in order to comprehensively investigate the vertical migration of MSA from summer layers, where MSA is originally deposited, to adjacent winter layers in polar ice. We find that the shallowest depth of MSA migration in our compilation varies over a wide range (∼ 2 to 400 m) and is positively correlated with snow accumulation rate and negatively correlated with ice concentration of Na+ (typically the most abundant marine cation). Although the considered soluble impurity transport model provides a useful mechanistic framework for studying MSA migration, it remains limited by inadequate constraints on key physico-chemical parameters – most notably, the diffusion coefficient of MSA in cold ice (DMS). We derive a simplified version of the model, which includes DMS as the sole parameter, in order to illuminate aspects of the migration process. Using this model, we show that the progressive phase alignment of MSA and Na+ concentration peaks observed along a high-resolution West Antarctic core is most consistent with 10−12 m2 s−1 DMS −11 m2 s−1, which is 1 order of magnitude greater than the DMS values previously estimated from laboratory studies. More generally, our data synthesis and model results suggest that (i) MSA migration may be fairly ubiquitous, particularly at coastal and (or) high-accumulation regions across Greenland and Antarctica; and (ii) can significantly change annual and multiyear MSA concentration averages. Thus, in most cases, caution should be exercised when interpreting polar ice core MSA records, although records that have undergone severe migration could still be useful for inferring decadal and lower-frequency climate variability.

The sources and transport pathways of aerosol species in Antarctica remain uncertain, partly due to limited seasonally resolved data from the harsh environment. Here, we examine the seasonal cycles of major ions in three high-accumulation West Antarctic ice cores for new information regarding the origin of aerosol species. A new method for continuous acidity measurement in ice cores is exploited to provide a comprehensive, charge-balance approach to assessing the major non-sea-salt (nss) species. The average nss-anion composition is 41% sulfate (SO42−), 36% nitrate (NO3−), 15% excess-chloride (ExCl−), and 8% methanesulfonic acid (MSA). Approximately 2% of the acid-anion content is neutralized by ammonium (NH4+), and the remainder is balanced by the acidity (Acy ≈ H+ − HCO3−). The annual cycle of NO3− shows a primary peak in summer and a secondary peak in late winter/spring that are consistent with previous air and snow studies in Antarctica. The origin of these peaks remains uncertain, however, and is an area of active research. A high correlation between NH4+ and black carbon (BC) suggests that a major source of NH4+ is midlatitude biomass burning rather than marine biomass decay, as previously assumed. The annual peak in excess chloride (ExCl−) coincides with the late-winter maximum in sea ice extent. Wintertime ExCl− is correlated with offshore sea ice concentrations and inversely correlated with temperature from nearby Byrd station. These observations suggest that the winter peak in ExCl− is an expression of fractionated sea-salt aerosol and that sea ice is therefore a major source of sea-salt aerosol in the region.

Assessments of climate sensitivity to projected greenhouse gas concentrations underpin environmental policy decisions, with such assessments often based on model simulations of climate during recent centuries and millennia1, 2, 3. These simulations depend critically on accurate records of past aerosol forcing from global-scale volcanic eruptions, reconstructed from measurements of sulphate deposition in ice cores4, 5, 6. Non-uniform transport and deposition of volcanic fallout mean that multiple records from a wide array of ice cores must be combined to create accurate reconstructions. Here we re-evaluated the record of volcanic sulphate deposition using a much more extensive array of Antarctic ice cores. In our new reconstruction, many additional records have been added and dating of previously published records corrected through precise synchronization to the annually dated West Antarctic Ice Sheet Divide ice core7, improving and extending the record throughout the Common Era. Whereas agreement with existing reconstructions is excellent after 1500, we found a substantially different history of volcanic aerosol deposition before 1500; for example, global aerosol forcing values from some of the largest eruptions (for example, 1257 and 1458) previously were overestimated by 20–30% and others underestimated by 20–50%.

The climate of West Antarctica is strongly influenced by remote forcing from the tropical Pacific. For example, recent surface warming over West Antarctica reflects atmospheric circulation changes over the Amundsen Sea, driven by an atmospheric Rossby wave response to tropical sea surface temperature (SST) anomalies. Here, it is demonstrated that tropical Pacific SST anomalies also influence the source and transport of marine-derived aerosols to the West Antarctic Ice Sheet. Using records from four firn cores collected along the Amundsen coast of West Antarctica, the relationship between sea ice–modulated chemical species and large-scale atmospheric variability in the tropical Pacific from 1979 to 2010 is investigated. Significant correlations are found between marine biogenic aerosols and sea salts, and SST and sea level pressure in the tropical Pacific. In particular, La Niña–like conditions generate an atmospheric Rossby wave response that influences atmospheric circulation over Pine Island Bay. Seasonal regression of atmospheric fields on methanesulfonic acid (MSA) reveals a reduction in onshore wind velocities in summer at Pine Island Bay, consistent with enhanced katabatic flow, polynya opening, and coastal dimethyl sulfide production. Seasonal regression of atmospheric fields on chloride (Cl−) reveals an intensification in onshore wind velocities in winter, consistent with sea salt transport from offshore source regions. Both the source and transport of marine aerosols to West Antarctica are found to be modulated by similar atmospheric dynamics in response to remote forcing. Finally, the regional ice-core array suggests that there is both a temporally and a spatially varying response to remote tropical forcing.

[1] We use an airborne-radar method, verified with ice-core accumulation records, to determine the spatiotemporal variations of snow accumulation over Thwaites Glacier, West Antarctica between 1980 and 2009. We also present a regional evaluation of modeled accumulation in Antarctica. Comparisons between radar-derived measurements and model outputs show that three global models capture the interannual variability well (r > 0.9), but a high-resolution regional model (RACMO2) has better absolute accuracy and captures the observed spatial variability (r = 0.86). Neither the measured nor modeled accumulation records over Thwaites Glacier show any trend since 1980. Although an increase in accumulation may potentially accompany the observed warming in the region, the projected trend is too small to detect over the 30 year record.

Glacial meltwater exports a unique type of organic carbon to marine systems, distinct from non-glacially derived riverine export, potentially capable of stimulating downstream marine primary productivity. Here, we describe for the first time the bulk-level dissolved organic carbon (DOC) and particulate organic carbon (POC) isotopic composition of glacial meltwater draining the Greenland ice sheet (GrIS). These data, in conjunction with an earlier study that investigated the molecular-level composition of GrIS dissolved organic matter, collectively describe the concentration, radiocarbon content, and lability of organic carbon in subglacial discharge from a land-terminating outlet glacier during a melt season. By scaling up our measurements across the ice sheet, we estimate that the annual DOC flux from the GrIS (0.08 Tg/y) is equivalent to that from a small Arctic river (discharge (Q) 3 /y), and that the annual POC flux from the GrIS (0.9 Tg/y) may be comparable to that of a large Arctic river (Q > 200 km 3 /y). The DOC flux is derived primarily from beneath the glacier (subglacial) (>75%) in the early season, and from surface ice-melt (up to 100%) transmitting through the subglacial environment at the peak of the meltseason. The POC flux is primarily derived from the subglacial environment throughout the meltseason. The early season (low flow) glacier discharge contains higher DOC concentrations (0.5–4.1 mg L −1 ), and exports more enriched carbon (Δ 14 C DOC ∼ −250‰) compared to the peak season (high flow) discharge, when the concentrations are lower (0.1–0.6 mg L −1 ) and the Δ 14 C is more depleted (Δ 14 C DOC ∼ −400‰). Conversely, the POC export (1.4–13.2 mg L −1 , Δ 14 C POC ∼ −250‰) shows no temporal variation in either concentration or radiocarbon content throughout the meltseason. Dissolved ion loads in concomitant samples reflected the seasonal evolution of the subglacial drainage system, confirming the influence of subglacial hydrology on the composition of the bulk carbon pools. Based on this work, we conclude that (1) different mechanisms control the DOC and POC flux from glacial systems; (2) chemically-distinct DOC pools are accessed by seasonally-evolving hydrological flow-paths; and (3) the GrIS can deliver labile carbon, which may also be 14 C-depleted, to downstream proglacial and marine environments.

The micronutrient iron is thought to limit primary production in large regions of the global ocean. Meltwater measurements suggest that the Greenland ice sheet serves as a significant source of potentially bioavailable iron to the surrounding coastal ocean

investigation of how regional SIC is recorded in the ice-sheet stratigraphy. Over the period 2002–2010 we find that the ice-sheet chemistry is significantly correlated with SIC variability within the AS and Pine Island Bay polynyas. Based on this result, we evaluate the use of icecore chemistry as a proxy for interannual polynya variability in this region, one of the largest and most persistent polynya areas in Antarctica. MSA concentrations correlate strongly with summer SIC within the polynya regions, consistent with MSA at this site being derived from marine biological productivity during the spring and summer. Cl – concentrations correlate strongly with winter SIC within the polynyas as well as some regions outside the polynyas, consistent with Cl – at this site originating primarily from winter sea-ice formation. Spatial correlations were generally insignificant outside of the polynya areas, with some notable exceptions. Ice-core glaciochemical records from this dynamic region thus may provide a proxy for reconstructing AS and Pine Island Bay polynya variability prior to the satellite era.

Cosmogenic 10Be in polar ice cores is a primary proxy for past solar activity. However, interpretation of the 10Be record is hindered by limited understanding of the physical processes governing its atmospheric transport and deposition to the ice sheets. This issue is addressed by evaluating two accurately dated, annually resolved ice core 10Be records against modern solar activity observations and instrumental and reanalysis climate data. The cores are sampled from the DSS site on Law Dome, East Antarctica (spanning 1936–2009) and the Das2 site, southeast Greenland (1936–2002), permitting inter-hemispheric comparisons. Concentrations at both DSS and Das2 are significantly correlated to the 11-yr solar cycle modulation of cosmic ray intensity, r xy = 0.54 with 95% CI [0.31; 0.70], and r xy = 0.45 with 95% CI [0.22; 0.62], respectively. For both sites, if fluxes are used instead of concentrations then correlations with solar activity decrease. The strength and spectral coherence of the solar activity signal in 10Be is enhanced when ice core records are combined from both Antarctica and Greenland. The amplitudes of the 11-yr solar cycles in the 10Be data appear inconsistent with the view that the ice sheets receive only 10Be produced at polar latitudes. Significant climate signals detected in the 10Be series include the zonal wave three pattern of atmospheric circulation at DSS, r xy = − 0.36 with 95% CI [−0.57; −0.10], and the North Atlantic Oscillation at Das2, r xy = − 0.42 with 95% CI [−0.64; −0.15]. The sensitivity of 10Be concentrations to modes of atmospheric circulation advises caution in the use of 10Be records from single sites in solar forcing reconstructions.

The Greenland ice sheet (GrIS) subglacial hydrological system may undergo a seasonal evolution, with significant geophysical and biogeochemical implications. We present results from a new isotope-mixing model to quantify the relative contributions of surface snow, glacial ice and delayed flow to the bulk meltwater discharge from a small (∼5 km2) land-terminating GrIS outlet glacier during melt onset (May) and at peak melt (July). We use radioactive (222Rn) and stable isotopes (18O, deuterium) to differentiate the water source contributions. Atmospherically derived 7Be further constrains meltwater transit time from the glacier surface to the ice margin. We show that (1) 222Rn is a promising tracer for glacial waters stored at the bed and (2) a quantitative chemical mixing model can be constructed by combining 222Rn and the stable water isotopes. Applying this model to the bulk subglacial outflow from our study area, we find a constant delayed-flow (stored) component from melt onset through peak melt. This component is diluted first by snowmelt and then by increasing glacial ice melt as the season progresses. Results from this pilot study are consistent with the hypothesis that subglacial drainage beneath land-terminating sections of the GrIS undergoes a seasonal evolution from a distributed to a channelized system.

Subsurface microbial oxidation of overridden soils and vegetation beneath glaciers and ice sheets may affect global carbon budgets on glacial–interglacial timescales. The likelihood and magnitude of this process depends on the chemical nature and reactivity of the subglacial organic carbon stores. We examined the composition of carbon pools associated with different regions of the Greenland ice sheet (subglacial, supraglacial, proglacial) in order to elucidate the type of dissolved organic matter (DOM) present in the subglacial discharge over a melt season. Electrospray ionization (ESI) Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry coupled to multivariate statistics permitted unprecedented molecular level characterization of this material and revealed that carbon pools associated with discrete glacial regions are comprised of different compound classes. Specifically, a larger proportion of protein-like compounds were observed in the supraglacial samples and in the early melt season (spring) subglacial discharge. In contrast, the late melt season (summer) subglacial discharge contained a greater fraction of lignin-like and other material presumably derived from underlying vegetation and soil. These results suggest (1) that the majority of supraglacial DOM originates from autochthonous microbial processes on the ice sheet surface, (2) that the subglacial DOM contains allochthonous carbon derived from overridden soils and vegetation as well as autochthonous carbon derived from in situ microbial metabolism, and (3) that the relative contribution of allochthonous and autochthonous material in subglacial discharge varies during the melt season. These conclusions are consistent with the hypothesis that, given sufficient time (e.g., overwinter storage), resident subglacial microbial communities may oxidize terrestrial material beneath the Greenland ice sheet.

Measurements of near-ice (

Measurements of near-ice (< 200 m) hydrography and near-terminus subglacial hydrology are lacking due in large part to the difficulty in working at the margin of calving glaciers. Here we pair detailed hydrographic and bathymetric measurements collected with an Autonomous Underwater Vehicle as close as 150 m from the ice/ocean interface of the Sarqardliup sermia/Sarqardleq Fjord system, West Greenland, with modeled and observed subglacial discharge locations and magnitudes. We find evidence of two main types of subsurface glacially modified water localized in space and with distinct properties that are consistent with runoff discharged at two locations along the grounded margin. These locations, in turn, correspond with two prominent subglacial subcatchments beneath Sarqardliup sermia. Thus, near-ice observations and subglacial discharge routing indicate that subglacial discharge from this glacier occurs at only two primary locations and gives rise to two distinct glacially modified waters. Furthermore, we show that the location with the largest discharge flux is associated with the lighter, fresher glacially modified watermass. This is qualitatively consistent with results from an idealized plume model.

Observed changes in the surface elevation of the Greenland ice sheet are caused by ice dynamics, basal elevation change, surface mass balance (SMB) variability, and by compaction of the overlying firn. The latter two contributions are quantified here using a firn model that includes compaction, meltwater percolation, and refreezing. The model is forced with surface mass fluxes and temperature from a regional climate model for the period 1960–2013. The model results agree with observations of surface density, density profiles from 62 firn cores, and altimetric observations from regions where ice-dynamical surface height changes are likely small. We find that the firn layer in the high interior is generally thickening slowly (1–5 cm yr−1). In the percolation and ablation areas, firn and SMB processes account for a surface elevation lowering of up to 20–50 cm yr−1. Most of this firn-induced marginal thinning is caused by an increase in melt since the mid-1990s, and partly compensated by an increase in the accumulation of fresh snow around most of the ice sheet. The total firn and ice volume change between 1980 and 2013 is estimated at −3900 ± 1030 km3 due to firn and SMB, corresponding to an ice-sheet average thinning of 2.32 ± 0.61 m. Most of this volume decrease occurred after 1995. The computed changes in surface elevation can be used to partition altimetrically observed volume change into surface mass balance and ice-dynamically related mass changes.

Annual and monthly snow accumulation for the Greenland Ice Sheet was derived from ECMWF forecasts [mainly 40-yr ECMWR Re-Analysis (ERA-40)] and further meteorological modeling. Modeled accumulation was validated using 58 ice core accumulation datasets across the ice sheet and was found to be 95% of the observed accumulation on average, with a mean correlation of 0.53 between modeled and observed. Many of the ice core datasets are new and are presented here for the first time. Central and northern interior parts of the ice sheet were found to be ∼10%–30% too dry in ERA-40, in line with earlier ECMWF analysis, although too much (>50% locally) snow accumulation was modeled for interior southern parts of Greenland. Nevertheless, 47 of 58 sites show significant correlation in temporal variability of modeled with observed accumulation. The model also captures the absolute amount of snow accumulation at several sites, most notably Das1 and Das2 in southeast Greenland. Mean modeled accumulation over the ice sheet was 0.279 (standard deviation 0.034) m yr−1 for 1958–2003 with no significant trend for either the ice sheet or any of the core sites. Unusually high accumulation in southeast Greenland in 2002/03 leads the authors to study meteorological synoptic forcing patterns and comment on the prospect of enhanced climate variability leading to more such events as a result of global warming. There is good agreement between precipitation measured at coastal meteorological stations in southern Greenland and accumulation modeled for adjacent regions of the ice sheet. There is no significant persistent relation between the North Atlantic Oscillation index and whole or southern Greenland accumulation.

Surface melting rarely occurs across most of the Antarctic ice sheet, away from the warmer coastal regions. Nonetheless, isolated melt features are preserved in the firn and ice in response to infrequent and short-lived melting events. An understanding of the formation and occurrence of these melt layers will help us to interpret records of past melt occurrences from polar ice cores such as the Siple Dome ice-core record from West Antarctica. A search in the near-surface firn in West Antarctica found that melt features are extremely rare, and consist of horizontal, laterally continuous, one to a few millimeter thick, ice layers with few air bubbles. The melt layers found date from the 1992/93 and 1991/92 summers. Field experiments to investigate changes in stratigraphy taking place during melt events reproduced melt features as seen in the natural stratigraphy. Melting conditions of varying intensity were created by passively heating the near-surface air for varying lengths of time inside a clear plastic hotbox. Melt layers formed due entirely to preferential flow and subsequent refreezing of meltwater from the surface into near-surface, fine-grained, crust layers. Continuous melt layers were formed experimentally when positive-degree-day values exceeded 1ºC-day, a value corresponding well with air-temperature records from automatic weather station sites where melt layers formed in the recent past.

This paper illustrates the components, capabilities, and some characteristic applications of the Woods Hole Oceanographic Institution Jetyak — a small autonomous surface vehicle (ASV) designed for the collection of oceanographic data from shallow or dangerous waters. The Jetyak is the result of custom modifications to a Mokai jet-powered kayak, including an A-frame and sea chest for installation of instrumentation, servo-driven controls and an Ardupilot autopilot for autonomous operation, an onboard computer for instrument control and data logging, and radios for wireless operation and communications. With these modifications, the Jetyak's cost of replacement is less than $15,000 (excluding the cost of instrumentation payload). The paper addresses the strengths and weaknesses of the Jetyak relative to piloted small boats and jetskis, autonomous underwater vehicles, and existing ASVs. Preliminary data are included from some shallow-water and dangerous Jetyak field campaigns in order to illustrate applications to which the Jetyak is well or uniquely suited.

Interior Antarctica is among the most remote places on Earth and was thought to be beyond the reach of human impacts when Amundsen and Scott raced to the South Pole in 1911. Here we show detailed measurements from an extensive array of 16 ice cores quantifying substantial toxic heavy metal lead pollution at South Pole and throughout Antarctica by 1889 - beating polar explorers by more than 22 years. Unlike the Arctic where lead pollution peaked in the 1970s, lead pollution in Antarctica was as high in the early 20 th century as at any time since industrialization. The similar timing and magnitude of changes in lead deposition across Antarctica, as well as the characteristic isotopic signature of Broken Hill lead found throughout the continent, suggest that this single emission source in southern Australia was responsible for the introduction of lead pollution into Antarctica at the end of the 19 th century and remains a significant source today. An estimated 660 €...t of industrial lead have been deposited over Antarctica during the past 130 years as a result of mid-latitude industrial emissions, with regional-to-global scale circulation likely modulating aerosol concentrations. Despite abatement efforts, significant lead pollution in Antarctica persists into the 21 st century.

Supraglacial lakes play an important role in establishing hydrological connections that allow lubricating seasonal melt water to reach the base of the Greenland Ice Sheet. Here we use new surface velocity observations to examine the influence of supraglacial lake drainages and surface melt rate on ice flow. We find large, spatially extensive speedups concurrent with times of lake drainage, showing that lakes play a key role in modulating regional ice flow. While surface meltwater is supplied to the bed via a geographically sparse network of moulins, the observed ice-flow enhancement suggests that this meltwater spreads widely over the ice-sheet bed. We also find that the complex spatial pattern of speedup is strongly determined by the combined influence of bed and surface topography on subglacial water flow. Thus, modeling of ice-sheet basal hydrology likely will require knowledge of bed topography resolved at scales (sub-kilometer) far finer than existing data (several km).

We present multi-model global datasets of ni- trogen and sulfate deposition covering time periods from 1850 to 2100, calculated within the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP). The computed deposition fluxes are compared to surface wet de- position and ice core measurements. We use a new dataset of wet deposition for 2000–2002 based on critical assessment of the quality of existing regional network data. We show that for present day (year 2000 ACCMIP time slice), the CCMIP results perform similarly to previously published multi-model assessments. For this time slice, we find a multi- model mean deposition of approximately 50 Tg(N) yr−1 from nitrogen oxide emissions, 60 Tg(N) yr−1 from ammo- nia emissions, and 83 Tg(S) yr−1 from sulfur emissions. The analysis of changes between 1980 and 2000 indicates significant differences between model and measurements over the United States but less so over Europe. This difference points towards a potential misrepresentation of 1980 NH3 missions over North America. Based on ice core records, the 1850 deposition fluxes agree well with Greenland ice cores, but the change between 1850 and 2000 seems to be overestimated in the Northern Hemisphere for both nitrogen and sulfur species. Using the Representative Concentration Pathways (RCPs) to define the projected climate and atmo- spheric chemistry related emissions and concentrations, we find large regional nitrogen deposition increases in 2100 in Latin America, Africa and parts of Asia under some of the scenarios considered. Increases in South Asia are especially large, and are seen in all scenarios, with 2100 values more than double their 2000 counterpart in some scenarios and reaching > 1300 mg(N) m−2 yr−1 averaged over regional to continental-scale regions in RCP 2.6 and 8.5, ∼ 30–50 % larger than the values in any region currently (circa 2000). However, sulfur deposition rates in 2100 are in all regions lower than in 2000 in all the RCPs. The new ACCMIP multi- model deposition dataset provides state-of-the-science, con- sistent and evaluated time slice (spanning 1850–2100) global gridded deposition fields for use in a wide range of climate and ecological studies., JRC.H.4-Monitoring Agricultural Resources

[1] Antarctic ice sheet surface melting can regionally influence ice shelf stability, mass balance, and glacier dynamics, in addition to modulating near-surface physical and chemical properties over wide areas. Here, we investigate variability in surface melting from 1999 to 2009 using radar backscatter time series from the SeaWinds scatterometer aboard the QuikSCAT satellite. These daily, continent-wide observations are explored in concert with in situ meteorological records to validate a threshold-based melt detection method. Radar backscatter decreases during melting are significantly correlated with in situ positive degree-days as well as meltwater production determined from energy balance modeling at Neumayer Station, East Antarctica. These results support the use of scatterometer data as a diagnostic indicator of melt intensity (i.e., the relative liquid water production during melting). Greater spatial and temporal melting detected relative to previous passive microwave-based studies is attributed to a higher sensitivity of the scatterometer instrument. Continental melt intensity variability can be explained in part by the dynamics of the Southern Annular Mode and the Southern Oscillation Index, and extreme melting events across the Ross Ice Shelf region may be associated with El Nino conditions. Furthermore, we find that the Antarctic Peninsula accounts for only 20% of Antarctic melt extent but greater than 50% of the total Antarctic melt intensity. Over most areas, annual melt duration and intensity are proportional. However, regional and localized distinctions exist where the melt intensity metric provides greater insight into melting dynamics than previously obtainable with other remote sensing techniques.