Search

Systematic measurements of atmospheric methane (CH4) mole fractions at the northern high latitudes only began in the early 1980s. Although CH4 measurements from Greenland ice cores consistently covered the period before ∼ 1900, no reliable observational record is available for the intermediate period. We newly report a data set of trace gases from the air trapped in firn (an intermediate stage between snow and glacial ice formation) collected at the NGRIP (North Greenland Ice Core Project) site in 2001. We also use a set of published firn air data at the NEEM (North Greenland Eemian ice Drilling) site. The two Arctic firn air data sets are analysed with a firn air transport model, which translates historical variations to depth profiles of trace gases in firn. We examine a variety of possible firn diffusivity profiles, using a suite of measured trace gases, and reconstruct the CH4 mole fraction by an iterative dating method. Although the reconstructions of the Arctic CH4 mole fraction before the mid-1970s still has large uncertainties (> 30 ppb – parts per billion), we find a relatively narrow range of atmospheric CH4 history that is consistent with both depth profiles of NGRIP and NEEM. The atmospheric CH4 history inferred by this study is more consistent with the atmospheric CH4 scenario prepared for the NEEM firn modelling than that for the CMIP6 (Climate Model Intercomparison Project Phase 6) experiments. Our study shows that the atmospheric CH4 scenario used for the NEEM firn modelling is considered to be the current best choice for the Arctic CH4 history, but it should not be used to tune firn air transport models until being verified by further measurements from sources such as the Arctic ice cores. Given the current difficulty in reconstructing the CH4 history with low uncertainty from the firn air data sets from Greenland, future sampling and measurements of ice cores at a high-accumulation site may be the only way to accurately reconstruct the atmospheric CH4 trend over the 20th century.

The variations of δO2/N2 and δAr/N2 in the Dome Fuji ice core were measured from 112 m (bubbly ice) to 2001 m (clathrate hydrate ice). Our method, combined with the low storage temperature of the samples (−50 ∘C), successfully excludes post-coring gas-loss fractionation signals from our data. From the bubbly ice to the middle of the bubble–clathrate transition zone (BCTZ) (112–800 m) and below the BCTZ (>1200 m), the δO2/N2 and δAr/N2 data exhibit orbital-scale variations similar to local summer insolation. The data in the lower BCTZ (800–1200 m) have large scatter, which may be caused by millimeter-scale inhomogeneity of air composition combined with finite sample lengths. The insolation signal originally recorded at the bubble close-off remains through the BCTZ, and the insolation signal may be reconstructed by analyzing long ice samples (more than 50 cm for the Dome Fuji core). In the clathrate hydrate zone, the scatter around the orbital-scale variability decreases with depth, indicating diffusive smoothing of δO2/N2 and δAr/N2. A simple gas diffusion model was used to reproduce the smoothing and thus constrain their permeation coefficients. The relationship between δAr/N2 and δO2/N2 is markedly different for the datasets representing bubble close-off (slope ∼ 0.5), bubble–clathrate hydrate transformation (∼1), and post-coring gas loss (∼0.2), suggesting that the contributions of the mass-independent and mass-dependent fractionation processes are different for those cases. The method and data presented here may be useful for improving the orbital dating of deep ice cores over the multiple glacial cycles and further studying non-insolation-driven signals (e.g., atmospheric composition) of these gases.

The 13th Symposium on Polar Science/Ordinary sessions [OM] Polar meteorology and glaciology, Wed. 16 Nov.

International audience; Precise ice-core chronologies are essential for identifying the timing and duration of polar climatic changes as well as their phasing with the changes in other parts of the globe. However, existing ice-core chronologies beyond the last 60 kyr show relatively large disagreements with each other and with U-Th chronologies of speleothems. Here, we constructed new ice and gas age scales for the Dome Fuji (DF) core (DF2021) over the last 207 kyr by combining a Bayesian dating model and firn densification model, constrained by various types of chronological and glaciological information including new δO2/N2 age markers, precise synchronization to other high-quality chronologies (volcanic, cosmogenic, and CH4 signals), and high-resolution δ15N of N2 (reflecting past firn thickness). The new chronology is tightly constrained by synchronization to other well-dated records for the last 60 kyr, whereas it is independent from other chronologies for the older period. For the last 60 kyr, the DF2021 chronology agrees with the layer-counted ice core chronologies (GICC05 and a part of WD2014) and U-Th chronologies of speleothems within ∼200 years. For the period 60-130 kyr BP, the timing of all Dansgaard-Oeschger warming events on DF2021 agree with those of corresponding events in the U-Th dated Chinese or European speleothems mostly within 1 kyr (well within 2σ uncertainty of DF2021). The excellent agreement suggests high accuracy of our chronology, and supports the assumption of negligible phasing between the past local summer solstice insolation and δO2/N2 fractionation at bubble close-off (the basis for constructing the δO2/N2 age markers). Between 130 and 207 kyr BP (penultimate glacial period), there is a lower degree of similarity between the variations in atmospheric CH4 and speleothem calcite δ18O than in the last glacial period, making the age comparison challenging. The comparison of DF2021 with 9 U-Th dates at 7 abrupt events shows the mean difference of -0.2 ± 0.8 kyr, which is within the DF2021 uncertainty (on the order of 1.5 kyr). The DF2021 chronology agrees with the AICC2012 chronology within 2 kyr except between ∼103 and 128 kyr BP where AICC2012 is likely too young by up to ∼4 kyr. By analyzing the lag between δ18O of O2 (δ18Oatm) on DF2021 and 65°N summer solstice insolation, the relatively large error in AICC2012 is found to originate in three δ18Oatm age markers assuming a constant lag, which are off by more than 3 kyr. The phasing between δ18Οatm and orbital forcing identified in this study may be useful for future ice core dating, especially where other age constraints are weak or lacking.