Firn Cores Research Articles

A 60-year ice-core record of regional climate from Adélie Land, coastal Antarctica

Abstract. A 22.4 m-long shallow firn core was extracted during the 2006/2007 field season from coastal Adélie Land. Annual layer counting based on subannual analyses of δ18O and major chemical components was combined with 5 reference years associated with nuclear tests and non-retreat of summer sea ice to build the initial ice-core chronology (1946–2006), stressing uncertain counting for 8 years. We focus here on the resulting δ18O and accumulation records. With an average value of 21.8 ± 6.9 cm w.e. yr−1, local accumulation shows multi-decadal variations peaking in the 1980s, but no long-term trend. Similar results are obtained for δ18O, also characterised by a remarkably low and variable amplitude of the seasonal cycle. The ice-core records are compared with regional records of temperature, stake area accumulation measurements and variations in sea-ice extent, and outputs from two models nudged to ERA (European Reanalysis) atmospheric reanalyses: the high-resolution atmospheric general circulation model (AGCM), including stable water isotopes ECHAM5-wiso (European Centre Hamburg model), and the regional atmospheric model Modèle Atmosphérique Régional (AR). A significant linear correlation is identified between decadal variations in δ18O and regional temperature. No significant relationship appears with regional sea-ice extent. A weak and significant correlation appears with Dumont d'Urville wind speed, increasing after 1979. The model-data comparison highlights the inadequacy of ECHAM5-wiso simulations prior to 1979, possibly due to the lack of data assimilation to constrain atmospheric reanalyses. Systematic biases are identified in the ECHAM5-wiso simulation, such as an overestimation of the mean accumulation rate and its interannual variability, a strong cold bias and an underestimation of the mean δ18O value and its interannual variability. As a result, relationships between simulated δ18O and temperature are weaker than observed. Such systematic precipitation and temperature biases are not displayed by MAR, suggesting that the model resolution plays a key role along the Antarctic ice sheet coastal topography. Interannual variations in ECHAM5-wiso temperature and precipitation accurately capture signals from meteorological data and stake observations and are used to refine the initial ice-core chronology within 2 years. After this adjustment, remarkable positive (negative) δ18O anomalies are identified in the ice-core record and the ECHAM5-wiso simulation in 1986 and 2002 (1998–1999), respectively. Despite uncertainties associated with post-deposition processes and signal-to-noise issues, in one single coastal ice-core record, we conclude that the S1C1 core can correctly capture major annual anomalies in δ18O as well as multi-decadal variations. These findings highlight the importance of improving the network of coastal high-resolution ice-core records, and stress the skills and limitations of atmospheric models for accumulation and δ18O in coastal Antarctic areas. This is particularly important for the overall East Antarctic ice sheet mass balance.

A Firn Densification Process in the High Accumulation Dome of Southeastern Greenland

ABSTRACTWe examine a firn core from a dome in southeast Greenland that exhibits distinct firn densification. The ice was −20.9 °C at 20 m depth, and the core gives an average accumulation rate of 1.0 m w.e. yr−1 in water equivalent. However, the close-off density of 830 kg m−3 occurs at 83.4–86.8 m depth, which is about 20-m shallower than that obtained from two empirical models. Where the density ρ > 750 kg m−3, the densification appears faster than that from the empirical models. As a result, compared to the empirical coefficient, the actual compactive viscosity coefficient is nonlinear and decreases at ρ > 750 kg m−3, indicating that the firn with a higher density is softer than that from the empirical result. We argue here that the high accumulation rate creates a high overburden pressure in a short time. Thus, the relative softness of the firn may arise from (1) there being not enough time to form bonds between grains as strong as those in a lower accumulation-rate area, and similarly, (2) the dislocation density in the firn being relatively high.

Firn Meltwater Retention on the Greenland Ice Sheet: A Model Comparison

Runoff has recently become the main source of mass loss from the Greenland Ice Sheet and is an important contributor to global sea level rise. Linking runoff to surface meltwater production is complex, as meltwater can be retained within the firn by refreezing or perennial liquid water storage. To constrain these uncertainties, the outputs of two offline snow/firn models of different complexity (IMAU-FDM and SNOWPACK) are compared to assess the sensitivity of meltwater retention to the model formulation (e.g., densification, irreducible water content, vertical resolution). Results indicate that model differences are largest in areas where firn aquifers form, i.e., particularly along the south-eastern margin of the ice sheet. The IMAU-FDM simulates higher densification rates for such climatic conditions and prescribes a lower irreducible water content than SNOWPACK. As a result, the model predicts substantially lower amounts of refreezing and liquid water storage. SNOWPACK performs better for this area, confirmed both by density profiles from firn cores and radar-inferred observations. Refreezing integrated over the entire ice sheet and averaged for the period 1960 - 2014 amounts to 216 Gt a-1 (IMAU-FDM) and 242 Gt a-1 (SNOWPACK), which is 41% and 46% of the total liquid water input (snowmelt and rainfall). The mean areal extents of perennial firn aquifers for 2010 - 2014 simulated by the models are 55,700 km2 (IMAU-FDM) and 90,200 km2 (SNOWPACK). Discrepancies between modeled firn profiles and observations emphasize the importance of processes currently not accounted for in most snow/firn models, such as vertical heterogeneous percolation, ponding of water on impermeable layers, lateral (sub-)surface water flow, and the issue of ill-constrained refreezing conditions at the base of firn aquifers.

Surface mass balance and water stable isotopes derived from firn cores on three ice rises, Fimbul Ice Shelf, Antarctica

Abstract. Three shallow firn cores were retrieved in the austral summers of 2011/12 and 2013/14 on the ice rises Kupol Ciolkovskogo (KC), Kupol Moskovskij (KM), and Blåskimen Island (BI), all part of Fimbul Ice Shelf (FIS) in western Dronning Maud Land (DML), Antarctica. The cores were dated back to 1958 (KC), 1995 (KM), and 1996 (BI) by annual layer counting using high-resolution oxygen isotope (δ18O) data, and by identifying volcanic horizons using non-sea-salt sulfate (nssSO42−) data. The water stable isotope records show that the atmospheric signature of the annual snow accumulation cycle is well preserved in the firn column, especially at KM and BI. We are able to determine the annual surface mass balance (SMB), as well as the mean SMB values between identified volcanic horizons. Average SMB at the KM and BI sites (0.68 and 0.70 mw. e. yr−1) was higher than at the KC site (0.24 mw. e. yr−1), and there was greater temporal variability as well. Trends in the SMB and δ18O records from the KC core over the period of 1958–2012 agree well with other previously investigated cores in the area, thus the KC site could be considered the most representative of the climate of the region. Cores from KM and BI appear to be more affected by local meteorological conditions and surface topography. Our results suggest that the ice rises are suitable sites for the retrieval of longer firn and ice cores, but that BI has the best preserved seasonal cycles of the three records and is thus the most optimal site for high-resolution studies of temporal variability of the climate signal. Deuterium excess data suggest a possible effect of seasonal moisture transport changes on the annual isotopic signal. In agreement with previous studies, large-scale atmospheric circulation patterns most likely provide the dominant influence on water stable isotope ratios preserved at the core sites.

Brief communication: Organochlorine pesticides in an archived firn core from Law Dome, East Antarctica

Abstract. Organochlorine pesticides (OCPs) were, for the first time, quantified in archived firn cores from East Antarctica representative of 1945–1957 and 1958–1967 (current era, C.E.). The core sections were melted under high-purity nitrogen atmosphere, and the meltwater was analysed. Methods allowed quantification of hexachlorocyclohexanes, heptachlor, trans-chlordane, dieldrin and endrin. While the core presented evidence of nominal contamination by modern-use chemicals, indicating handling and/or storage contamination, legacy OCP concentrations and deposition rates reported are orders of magnitude lower than those from Arctic regions, lending support for their validity. The study further provides a description of equipment used and suggests methods to overcome logistical challenges associated with trace organic contaminant detection in polar regions.

Seasonal and spatial variations of 17Oexcess and dexcess in Antarctic precipitation: Insights from an intermediate complexity isotope model

AbstractAn intermediate complexity model (ICM) is used to investigate the sensitivity of water isotope ratios in precipitation, including 17Oexcess, to climate variations in the Southern Hemisphere. The ICM is forced with boundary conditions from seasonal National Centers for Environmental Prediction/Department of Energy II reanalysis data. Perturbations to the surface temperature and humidity fields are used to investigate the isotopic sensitivity. The response of 17Oexcess to a uniform temperature change is insignificant over the ocean, while there is a large magnitude response over the ice sheet, particularly in East Antarctica. A decrease of ocean surface relative humidity produces increased 17Oexcess and dexcess, with a coherent response over both the ocean and Antarctica. For interior East Antarctica, the model simulates a seasonal cycle in 17Oexcess that is positively correlated with δ18O and of large magnitude (~50 per meg), consistent with the observations from Vostok. The seasonal cycle in 17Oexcess for interior West Antarctica is predicted to be considerably smaller in magnitude (12 per meg), and is negatively correlated with δ18O, consistent with new data from a firn core near the West Antarctic Ice Sheet Divide site. Over the ocean, the ICM predicts much smaller seasonal cycles in 17Oexcess. Oceanic source changes (i.e., humidity) are insufficient to explain the amplitude of the simulated seasonal cycle over the Antarctic continent. Spatial differences in the seasonal response of 17Oexcess to local temperature reflect the balance of equilibrium and kinetic fractionation during snow formation.

Layering of surface snow and firn at Kohnen Station, Antarctica: Noise or seasonal signal?

The density of firn is an important property for monitoring and modeling the ice sheet as well as to model the pore close-off and thus to interpret ice core-based greenhouse gas records. One feature, which is still in debate, is the potential existence of an annual cycle of firn density in low-accumulation regions. Several studies describe or assume seasonally successive density layers, horizontally evenly distributed, as seen in radar data. On the other hand, high-resolution density measurements on firn cores in Antarctica and Greenland showed no clear seasonal cycle in the top few meters. A major caveat of most existing snow-pit and firn-core based studies is that they represent one vertical profile from a laterally heterogeneous density field. To overcome this, we created an extensive dataset of horizontal and vertical density data at Kohnen Station, Dronning Maud Land on the East Antarctic Plateau. We drilled and analyzed three 90 m long firn cores as well as 160 one meter long vertical profiles from two elongated snow trenches to obtain a two dimensional view of the density variations. The analysis of the 45 m wide and 1 m deep density fields reveals a seasonal cycle in density. However, the seasonality is overprinted by strong stratigraphic noise, making it invisible when analyzing single firn cores. Our density dataset extends the view from the local ice-core perspective to a hundred meter scale and thus supports linking spatially integrating methods such as radar and seismic studies to ice and firn cores.

An Optical Dye Method for Continuous Determination of Acidity in Ice Cores.

The pH of polar ice is important for the stability and mobility of impurities in ice cores and can be strongly influenced by volcanic eruptions or anthropogenic emissions. We present a simple optical method for continuous determination of acidity in ice cores based on spectroscopically determined color changes of two common pH-indicator dyes, bromophenol blue, and chlorophenol red. The sealed-system method described here is not equilibrated with CO2, making it simpler than existing methods for pH determination in ice cores and offering a 10-90% peak response time of 45 s and a combined uncertainty of 9%. The method is applied to Holocene ice core sections from Greenland and Antarctica and compared to standard techniques such as electrical conductivity measurement (ECM) conducted on the solid ice, and electrolytic meltwater conductivity, EMWC. Acidity measured in the Greenland NGRIP ice core shows good agreement with acidity calculated from ion chromatography. Conductivity and dye-based acidity Hdye+ are found to be highly correlated in the Greenland NEGIS firn core (75.38° N, 35.56° W), with all signals greater than 3σ variability coinciding with either volcanic eruptions or possible wild fire activity. In contrast, the Antarctic Roosevelt Island ice core (79.36° S, 161.71° W) features an anticorrelation between conductivity and Hdye+, likely due to strong influence of marine salts.

Seasonal accumulation of persistent organic pollutants on a high altitude glacier in the Eastern Alps

The seasonal accumulations of perfluorinated substances (PFAS), polybrominated diphenyl ethers (PBDE) and polycyclic aromatic hydrocarbons (PAH) were measured in a 10 m shallow firn core from a high altitude glacier at Mt. Ortles (Italy, 3830 m above sea level) in South Tyrol in the Italian Eastern Alps. The most abundant persistent organic pollutants of each group were perfluorobutanoic acid (PFBA), perfluorooctanoic acid (PFOA), perfluorononanoic acid (PFNA) (for PFASs); BDE 47, BDE 99, BDE 209 (for PBDEs) and phenanthrene (PHE), fluoranthene (FLA) and pyrene (PYR) (for PAHs). All compounds show different extents of seasonality, with higher accumulation during summer time compared to winter. This seasonal difference mainly reflects meteorological conditions with a low and stable atmospheric boundary layer in winter and strong convective activity in summer, transformation processes during the transport of chemicals and/or post-depositional alterations. Change in the composition of the water-soluble PFCAs demonstrates the influence of meltwater percolation through the firn layers.

Variabilidade do conteúdo iônico da neve e do firn ao longo de um transecto antártico

Com o objetivo de interpretar a variabilidade no conteúdo iônico da neve e do firn entre Patriot Hills (80º18’S, 81º21’W) e o Polo Sul Geográfico, foram determinadas as concentrações dos íons majoritários e as razões isotópicas das 200 primeiras frações de cinco testemunhos de neve e firn (com profundidades de até 46 m), durante a travessia chileno-brasileira ocorrida no verão de 2004–2005. Após as amostras serem limpas (em câmara fria) e derretidas (em sala limpa classe 100), os íons Na+ , K+ , Mg2+, Ca2+, MS- (CH3 SO3 - ), Cl- , NO3 - e SO4 2- foram analisados por cromatografia iônica, com limites de detecção entre 0,3 e 1,8 µg L-1, dependendo do íon. A razão isotópica deutério/hidrogênio (δD) foi determinada com precisão de 0,5‰. Dataram-se as amostras, com precisão anual e exatidão de ±2 anos, a partir da contagem anual de camadas de neve e de firn dos perfis de Na+ , nssSO4 2- (sulfato não proveniente do sal marinho, do inglês non-sea-salt sulfate) e δD. A partir da datação obteve-se as taxas médias de acumulação líquida de neve para os sítios de cada testemunho de neve e firn, as quais apresentaram uma correlação negativa (r > 0,4, em módulo) à medida que se aumentavam a elevação e a distância da costa, respondendo, ainda, às feições superficiais locais e à ocorrência de ventos catabáticos. Constataram-se aerossóis do sal marinho com alteração na razão Na+ /Cl- no registro dos cinco testemunhos de neve e firn, mas a presença de aerossóis produzidos na formação do gelo marinho só foi identificada nos sítios mais próximos à costa. Os registros de NO3 - , assim como os de MS- , indicam a provável ocorrência de processos pós- -deposicionais, ligados tanto a eventos vulcânicos, como a reações fotoquímicas acentuadas por superfícies com formação de esmalte de gelo.

Spatial and temporal variability of snow chemical composition and accumulation rate at Talos Dome site (East Antarctica)

Five snow pits and five firn cores were sampled during the 2003–2004 Italian Antarctic Campaign at Talos Dome (East Antarctica), where a deep ice core (TALDICE, TALos Dome Ice CorE, 1650m depth) was drilled in 2005–2008 and analyzed for ionic content.Particular attention is spent in applying decontamination procedures to the firn cores, as core sections were stored for approximately 10years before analysis. By considering the snow pit samples to be unperturbed, the comparison with firn core samples from the same location shows that ammonium, nitrate and MSA are affected by storage post-depositional losses. All the other measured ions are confirmed to be irreversibly deposited in the snow layer. The removal of the most external layers (few centimeters) from the firn core sections is proved to be an effective decontamination procedure.High-resolution profiles of seasonal markers (nitrate, sulfate and MSA) allow a reliable stratigraphic dating and a seasonal characterization of the samples. The calculated mean accumulation-rate values range from 70 to 85mmw.e.year−1, in the period 2003–1973 with small differences between two sectors: 70–74mmw.e.year−1 in the NNE sector (spanning 2003–1996years) and 81–92mmw.e.year−1 in the SSW sector (spanning 2003–1980years). This evidence is interpreted as a coupled effect of wind-driven redistribution processes in accumulation/ablation areas.Statistical treatment applied to the concentration values of the snow pits and firn cores samples collected in different points reveals a larger temporal variability than spatial one both in terms of concentration of chemical markers and annual accumulation.The low spatial variability of the accumulation rate and chemical composition measured in the five sites demonstrates that the TALDICE ice core paleo-environmental and paleo-climatic stratigraphies can be considered as reliably representative for the Talos Dome area.

Densification of layered firn in the ice sheet at Dome Fuji, Antarctica

ABSTRACTIn order to better understand the densification of polar firn, firn cores from the three sites within ~10 km of Dome Fuji, Antarctica, were investigated using surrogates of density: dielectric permittivities εv and εh at microwave frequencies with electrical fields in the vertical and horizontal planes respectively. Dielectric anisotropy Δε (=εv − εh) was then examined as a surrogate of the anisotropic geometry of firn. We find that layered densification is explained as a result of complex effects of two phenomena that commonly occur at the three sites. Basically, layers with initially smaller density and smaller geometrical anisotropy deform preferentially throughout the densification process due to textural effects. Second, layers having a higher concentration of Cl− ions deform preferentially during a limited period from the near surface depths until smoothing out of layered Cl− ions by diffusion. We hypothesize that Cl− ions dissociated from sea salts soften firn due to modulation of dislocation movement. Moreover, firn differs markedly across the three sites in terms of strength of geometrical anisotropy, mean rate of densification and density fluctuation. We hypothesize that these differences are caused by textural effects resulting from differences in depositional conditions within various spatial scales.

Surface Snow, Firn and Ice Core Composition in Polar Areas in Relation to Atmospheric Aerosol and Gas Concentrations: Critical Aspects

The paper addresses some of the problems surrounding the relation between ice core chemical signals and atmospheric chemical composition in polar areas. The topic is important as the reconstruction of past climate and past atmospheric chemical composition is based on the assumption that chemical concentrations in the air, snow, firn and ice core are correlated. Ice core interpretation of aerosol is more straightforward than that of reactive gases. The transfer functions of gaseous species strongly interacting with ice are complex and additional field and laboratory experiments are required. Ice core chemical signals depend on the chemical composition of precipitations, which are related to the physics of precipitation formation, the chemical composition of the atmosphere, and post-depositional processes. Published papers reporting data on the chemical composition of snow seldom consider the fact that crystal formation and growth in cloud (rimed or unrimed) or near the ground (clear-sky precipitations), hoar-frost formation and surface riming determine different chemical concentrations, even assuming constant background concentration in the atmosphere. This paper discusses the physical and chemical processes affecting the formation of precipitations in polar areas, and the process of scavenging gases from non-growing and growing crystals. Attention is mainly focused on the processes involving nitrate anion in snow, hoar frost and firn. Knowledge of the chemical relationship between surface snow and atmospheric chemical concentration could be enhanced by considering specific events, such as snow falling from cloud, clear sky precipitation, and surface hoar or riming surface, with simultaneous air sampling. In conclusion, field and laboratory experiments are still required to study the scavenging processes during crystal formation.

Taxonomic characterization and the bio-potential of bacteria isolated from glacier ice cores in the High Arctic.

Glacier ice and firn cores have ecological and biotechnological importance. The present study is aimed at characterizing bacteria in crustal ice cores from Svalbard, the Arctic. Counts of viable isolates ranged from 10 to 7000 CFU/ml (mean 803 CFU/ml) while the total bacterial numbers ranged from 7.20 × 10(4) to 2.59 × 10(7) cells ml(-1) (mean 3.12 × 10(6) cells ml(-1) ). Based on 16S rDNA sequence data, the identified species belonged to seven species, namely Bacillus barbaricus, Pseudomonas orientalis, Pseudomonas oryzihabitans, Pseudomonas fluorescens, Pseudomonas syncyanea, Sphingomonas dokdonensis, and Sphingomonas phyllosphaerae, with a sequence similarity ranging between 93.5 and 99.9% with taxa present in the database. The isolates exhibited unique phenotypic properties, and three isolates (MLB-2, MLB-5, and MLB-9) are novel species, yet to be described. To the best of our knowledge, this is the first report on characterization of cultured bacterial communities from Svalbard ice cores. We conclude that high lipase, protease, cellulase, amylase, and urease activities expressed by most of the isolates provide a clue to the potential industrial applications of these organisms. These microbes, producing cold-adapted enzymes may provide an opportunity for biotechnological research.

Elevation change of the Greenland Ice Sheet due to surface mass balance and firn processes, 1960–2014

Abstract. Observed changes in the surface elevation of the Greenland Ice Sheet are caused by ice dynamics, basal elevation change, basal melt, surface mass balance (SMB) variability, and by compaction of the overlying firn. The last 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–2014. 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. In areas with strong surface melt, the firn model overestimates density. 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 2014 is estimated at −3295 ± 1030 km3 due to firn and SMB changes, corresponding to an ice-sheet average thinning of 1.96 ± 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.

Low-background neutron activation analysis: a powerful tool for atmospheric mineral dust analysis in ice cores

The application of instrumental neutron activation analysis (INAA) for multi-elemental analysis of samples of extremely reduced mass such as dust samples extracted from ice cores requires specific efforts towards the development of a “low level counting” analytical technique. An analytical protocol specifically designed for this kind of samples, based on low-background INAA (LBNAA) is here presented. A first application of the method was successfully performed on samples from the new alpine firn core NextData-LYS12. Sub-ng detection limits were reached for many elements. According to this point the technique is also potentially suitable to be applied to polar ice core samples.

Constraints on the δ<sup>2</sup>H diffusion rate in firn from field measurements at Summit, Greenland

Abstract. We performed detailed 2H isotope diffusion measurements in the upper 3 m of firn at Summit, Greenland. Using a small snow gun, a thin snow layer was formed from 2H-enriched water over a 6 × 6 m2 area. We followed the diffusion process, quantified as the increase of the δ2H diffusion length, over a 4-year period, by retrieving the layer once per year by drilling a firn core and slicing it into 1 cm layers and measuring the δ2H signal of these layers. We compared our experimental findings to calculations based on the model by Johnsen et al. (2000) and found substantial differences. The diffusion length in our experiments increased much less over the years than in the model. We discuss the possible causes for this discrepancy and conclude that several aspects of the diffusion process in firn are still poorly constrained, in particular the tortuosity.

Unlocking annual firn layer water equivalents from ground-penetrating radar data on an Alpine glacier

Abstract. The spatial representation of accumulation measurements is a major limitation for current glacier mass balance monitoring approaches. Here, we present a method for estimating annual accumulation rates on a temperate Alpine glacier based on the interpretation of internal reflection horizons (IRHs) in helicopter-borne ground-penetrating radar (GPR) data. For each individual GPR measurement, the signal travel time is combined with a simple model for firn densification and refreezing of meltwater. The model is calibrated at locations where GPR repeat measurements are available in two subsequent years and the densification can be tracked over time. Two 10.5 m long firn cores provide a reference for the density and chronology of firn layers. Thereby, IRHs correspond to density maxima, but not exclusively to former summer glacier surfaces. Along GPR profile sections from across the accumulation area we obtain the water equivalent (w.e.) of several annual firn layers. Because deeper IRHs could be tracked over shorter distances, the total length of analysed profile sections varies from 7.3 km for the uppermost accumulation layer (2011) to 0.1 km for the deepest (i.e. oldest) layer (2006). According to model results, refreezing accounts for 10% of the density increase over time and depth, and for 2% of the water equivalent. The strongest limitation to our method is the dependence on layer chronology assumptions. We show that GPR can be used not only to complement existing mass balance monitoring programmes on temperate glaciers but also to retrospectively extend newly initiated time series.

Climatic signals from 76 shallow firn cores in Dronning Maud Land, East Antarctica

Abstract. The spatial and temporal distribution of surface mass balance (SMB) and δ18O were investigated in the first comprehensive study of a set of 76 firn cores retrieved by various expeditions during the past 3 decades in Dronning Maud Land, East Antarctica. The large number of cores was used to calculate stacked records of SMB and δ18O, which considerably increased the signal-to-noise ratio compared to earlier studies and facilitated the detection of climatic signals. Considerable differences between cores from the interior plateau and the coastal cores were found. The δ18O of both the plateau and the ice shelf cores exhibit a slight positive trend over the second half of the 20th century. In the corresponding period, the SMB has a negative trend in the ice shelf cores, but increases on the plateau. Comparison with meteorological data from Neumayer Station revealed that for the ice shelf regions, atmospheric dynamic effects are more important than thermodynamics while on the plateau; the temporal variations of SMB and δ18O occur mostly in parallel, and thus can be explained by thermodynamic effects. The Southern Annular Mode (SAM) has exhibited a positive trend since the mid-1960s, which is assumed to lead to a cooling of East Antarctica. This is not confirmed by the firn core data in our data set. Changes in the atmospheric circulation that result in a changed seasonal distribution of precipitation/accumulation could partly explain the observed features in the ice shelf cores.

Observing and modeling the influence of layering on bubble trapping in polar firn

AbstractInterpretation of ice core trace gas records depends on an accurate understanding of the processes that smooth the atmospheric signal in the firn. Much work has been done to understand the processes affecting air transport in the open pores of the firn, but a paucity of data from air trapped in bubbles in the firn‐ice transition region has limited the ability to constrain the effect of bubble closure processes. Here we present high‐resolution measurements of firn density, methane concentrations, nitrogen isotopes, and total air content that show layering in the firn‐ice transition region at the West Antarctic Ice Sheet (WAIS) Divide ice core site. Using the notion that bubble trapping is a stochastic process, we derive a new parameterization for closed porosity that incorporates the effects of layering in a steady state firn modeling approach. We include the process of bubble trapping into an open‐porosity firn air transport model and obtain a good fit to the firn core data. We find that layering broadens the depth range over which bubbles are trapped, widens the modeled gas age distribution of air in closed bubbles, reduces the mean gas age of air in closed bubbles, and introduces stratigraphic irregularities in the gas age scale that have a peak‐to‐peak variability of ~10 years at WAIS Divide. For a more complete understanding of gas occlusion and its impact on ice core records, we suggest that this experiment be repeated at sites climatically different from WAIS Divide, for example, on the East Antarctic plateau.