NEEM Site Research Articles

Towards reconstructing the Arctic atmospheric methane history over the 20th century: measurement and modelling results for the North Greenland Ice Core Project firn

Abstract. 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.

Atmospheric iron supply and marine productivity in the glacial North Pacific Ocean

Abstract. Iron (Fe) is a key element in the Earth climate system, as it can enhance marine primary productivity in the high-nutrient low-chlorophyll (HNLC) regions where, despite a high concentration of major nutrients, chlorophyll production is low due to iron limitation. Eolian mineral dust represents one of the main Fe sources to the oceans; thus, quantifying its variability over the last glacial cycle is crucial to evaluate its role in strengthening the biological carbon pump. Polar ice cores, which preserve detailed climate records in their stratigraphy, provide a sensitive and continuous archive for reconstructing past eolian Fe fluxes. Here, we show the Northern Hemisphere Fe record retrieved from the NEEM ice core (Greenland), which offers a unique opportunity to reconstruct the past Fe fluxes in a portion of the Arctic over the last 108 kyr. Holocene Fe fluxes (0.042–11.7 ka, 0.5 mg m−2 yr−1) at the NEEM site were 4 times lower than the average recorded over the last glacial period (11.7–108 ka, 2.0 mg m−2 yr−1), whereas they were greater during the Last Glacial Maximum (LGM; 14.5–26.5 ka, 3.6 mg m−2 yr−1) and Marine Isotope Stage 4 (MIS 4; 60–71 ka, 5.8 mg m−2 yr−1). Comparing the NEEM Fe record with paleoceanographic records retrieved from the HNLC North Pacific, we found that the coldest periods, characterized by the highest Fe fluxes, were distinguished by low marine primary productivity in the subarctic Pacific Ocean, likely due to the greater sea ice extent and the absence of major nutrients upwelling. This supports the hypothesis that Fe fertilization during colder and dustier periods (i.e., LGM and MIS 4) was more effective in other regions, such as the midlatitude North Pacific, where a closer relationship between marine productivity and the NEEM Fe fluxes was observed.

Multi-channel and multi-polarization radar measurements around the NEEM site

Abstract. Ice properties inferred from multi-polarization measurements, such as birefringence and crystal orientation fabric (COF), can provide insight into ice strain, viscosity, and ice flow. In 2008, the Center for Remote Sensing of Ice Sheets (CReSIS) used a ground-based VHF (very high frequency) radar to take multi-channel and multi-polarization measurements around the NEEM (North Greenland Eemian Ice Drilling) site. The system operated with 30 MHz bandwidth at a center frequency of 150 MHz. This paper describes the radar system, antenna configurations, data collection, and processing and analysis of this data set. Within the framework derived from uniaxial ice crystal model, we found that ice birefringence dominates the power variation patterns of co-polarization and cross-polarization measurements in the area of 100 km2 around the ice core site. The phase shift between ordinary and extraordinary waves increases nonlinearly with depth. The ice optic axis lies in planes that are close to the vertical plane and perpendicular or parallel to the ice divide depending on depth. The ice optic axis has an average tilt angle of about 11.6° vertically, and its plane may rotate either clockwise or counterclockwise by about 10° across the 100 km2 area, and at a specific location the plane may rotate slightly counterclockwise as depth increases. Comparisons between the radar observations, simulations, and ice core fabric data are in very good agreement. We calculated the effective colatitude at different depths by using azimuth and colatitude measurements of the c axis of ice crystals. We obtained an average effective c axis tilt angle of 9.6° from the vertical axis, very comparable to the average optic axis tilt angle estimated from radar polarization measurements. The comparisons give us confidence in applying this polarimetric radio echo sounding technique to infer profiles of ice fabric in locations where there are no ice core measurements.

Climate information preserved in seasonal water isotope at NEEM: relationships with temperature, circulation and sea ice

Abstract. Analyzing seasonally resolved δ18O ice core data can aid the interpretation of the climate information in ice cores, also providing insights into factors governing the δ18O signal that cannot be deciphered by investigating the annual δ18O data only. However, the seasonal isotope signal has not yet been investigated in northern Greenland, e.g., at the NEEM (North Greenland Eemian Ice Drilling) ice core drill site. Here, we analyze seasonally resolved δ18O data from four shallow NEEM ice cores covering the last 150 years. Based on correlation analysis with observed temperature, we attribute about 70 and 30 % of annual accumulation to summer and winter, respectively. The NEEM summer δ18O signal correlates strongly with summer western Greenland coastal temperature and with the first principal component (PC1) of summer δ18O from multiple seasonally resolved ice cores from central/southern Greenland. However, there are no significant correlations between NEEM winter δ18O data and western Greenland coastal winter temperature or southern/central Greenland winter δ18O PC1. The stronger correlation with temperature during summer and the dominance of summer precipitation skew the annual δ18O signal in NEEM. The strong footprint of temperature in NEEM summer δ18O record also suggests that the summer δ18O record rather than the winter δ18O record is a better temperature proxy at the NEEM site. Despite the dominant signal of the North Atlantic Oscillation (NAO) and the Atlantic Multidecadal Oscillation (AMO) in the central–southern ice core data, both NAO and AMO exert weak influences on NEEM seasonal δ18O variations. The NEEM seasonal δ18O is found to be highly correlated with Baffin Bay sea ice concentration (SIC) in the satellite observation period (1979–2004), suggesting a connection of the sea ice extent with δ18O at NEEM. NEEM winter δ18O significantly correlates with SIC even for the period prior to satellite observation (1901–1978). The NEEM winter δ18O may reflect sea ice variations of Baffin Bay rather than temperature itself. This study shows that seasonally resolved δ18O records, especially for sites with a seasonal precipitation bias such as NEEM, provide a better understanding of how changing air temperature and circulation patterns are associated with the variability in the δ18O records.

Improvement and performance testing of melting system for measurement of trace elements in firn core drilled at NEEM site, Greenland

ABSTRACTA melting system was improved for measurement of trace elements in a firn core drilled at the North Greenland Eemian Ice Drilling (NEEM) site (77.45°N, 51.06°W). The melting heads of the melting system were specifically constructed with fluorinated ethylene propylene (FEP)-coated copper (type A) and pure nickel (type B) to reduce percolation of melted samples. Firn core sections with densities in the range of ~0.76–0.86 g cm−3 (depth interval: ~58.10–88.83 m) and ~0.56–0.76 g cm−3 (depth interval: ~17.32–58.10 m) were successfully decontaminated by the melting system with type A and B melting heads, respectively, as verified by the volume ratios of melted samples from inner and outer zones and video recorded during the melting procedures. Systematic performance tests of the type A melting head melter were performed using artificial ice cores. The procedural blanks and detection limits of V, Cr, Co, Rb, Cd, Ba, Tl, Th and U were generally lower than those in previous studies, whereas the levels of Mn, Sr and Pb were similar to those reported elsewhere. The concentration ratios of these elements from inner and outer zones of the melting head indicated that the melting system had good decontamination efficiency. Finally, the concentrations of elements in the NEEM firn core measured following mechanical removal of the outer layer with chisels and the type A melting head melter showed a good agreement with each other.

Changes to the chemical state of the Northern Hemisphere atmosphere during the second half of the twentieth century

Abstract. The NOx (NO and NO2) and HOx (OH and HO2) budgets of the atmosphere exert a major influence on atmospheric composition, controlling removal of primary pollutants and formation of a wide range of secondary products, including ozone, that can influence human health and climate. However, there remain large uncertainties in the changes to these budgets over recent decades. Due to their short atmospheric lifetimes, NOx and HOx are highly variable in space and time, and so the measurements of these species are of limited value for examining long-term, large-scale changes to their budgets. Here, we take an alternative approach by examining long-term atmospheric trends of alkyl nitrates, the production efficiency of which is dependent on the atmospheric [NO] ∕ [HO2] ratio. We derive long-term trends in the alkyl nitrates from measurements in firn air from the NEEM site, Greenland. Their mixing ratios increased by a factor of 3–5 between the 1970s and 1990s. This was followed by a steep decline to the sampling date of 2008. Moreover, we examine how the trends in the alkyl nitrates compare to similarly derived trends in their parent alkanes (i.e. the alkanes which, when oxidised in the presence of NOx, lead to the formation of the alkyl nitrates). The ratios of the alkyl nitrates to their parent alkanes increased from around 1970 to the late 1990s. This is consistent with large changes to the [NO] ∕ [HO2] ratio in the Northern Hemisphere atmosphere during this period. Alternatively, they could represent changes to concentrations of the hydroxyl radical, OH, or to the transport time of the air masses from source regions to the Arctic.

Evaluating the skills of isotope‐enabled general circulation models against in situ atmospheric water vapor isotope observations

AbstractThe skills of isotope‐enabled general circulation models are evaluated against atmospheric water vapor isotopes. We have combined in situ observations of surface water vapor isotopes spanning multiple field seasons (2010, 2011, and 2012) from the top of the Greenland Ice Sheet (NEEM site: 77.45°N, 51.05°W, 2484 m above sea level) with observations from the marine boundary layer of the North Atlantic and Arctic Ocean (Bermuda Islands 32.26°N, 64.88°W, year: 2012; south coast of Iceland 63.83°N, 21.47°W, year: 2012; South Greenland 61.21°N, 47.17°W, year: 2012; Svalbard 78.92°N, 11.92°E, year: 2014). This allows us to benchmark the ability to simulate the daily water vapor isotope variations from five different simulations using isotope‐enabled general circulation models. Our model‐data comparison documents clear isotope biases both on top of the Greenland Ice Sheet (1–11‰ for δ18O and 4–19‰ for d‐excess depending on model and season) and in the marine boundary layer (maximum differences for the following: Bermuda δ18O = ~1‰, d‐excess = ~3‰; South coast of Iceland δ18O = ~2‰, d‐excess = ~ 5‰; South Greenland δ18O = ~4‰, d‐excess = ~7‰; Svalbard δ18O = ~2‰, d‐excess = ~7‰). We find that the simulated isotope biases are not just explained by simulated biases in temperature and humidity. Instead, we argue that these isotope biases are related to a poor simulation of the spatial structure of the marine boundary layer water vapor isotopic composition. Furthermore, we specifically show that the marine boundary layer water vapor isotopes of the Baffin Bay region show strong influence on the water vapor isotopes at the NEEM deep ice core‐drilling site in northwest Greenland. Our evaluation of the simulations using isotope‐enabled general circulation models also documents wide intermodel spatial variability in the Arctic. This stresses the importance of a coordinated water vapor isotope‐monitoring network in order to discriminate amongst these model behaviors.

Accumulation Rates during 1311–2011 CE in North-Central Greenland Derived from Air-Borne Radar Data

Radar–detected internal layering contains information on past accumulation rates and patterns. In this study, we assume that the radar layers are isochrones, and use the layer stratigraphy in combination with ice-core measurements and numerical methods to retrieve accumulation information for the northern part of central Greenland. Measurements of the dielectric properties of an ice core from the NEEM (North Greenland Eemian Ice Drilling) site, allow for correlation of the radar layers with volcanic horizons to obtain an accurate age of the layers. We obtain accumulation patterns averaged over 100 a for the period 1311–2011. Our results show a clear trend of high accumulation rates west of the ice divide and low accumulation rates east of the ice divide. At the NEEM site the accumulation pattern is persistent during our study period and only small temporal variations occur in the accumulation rate. However, from approximately 200 km south of the NEEM drill site, the accumulation rate shows temporal variations based on our centennial averages. We attribute this variation to shifts in the location of the high–low accumulation boundary that usually is aligned with the ice divide, but appears to have moved across the divide in the past.

Greenland during the last interglacial: the relative importance of insolation and oceanic changes

Abstract. Insolation changes during the Eemian (the last interglacial period, 129 000–116 000 years before present) resulted in warmer than present conditions in the Arctic region. The NEEM ice core record suggests warming of 8 ± 4 K in northwestern Greenland based on stable water isotopes. Here we use general circulation model experiments to investigate the causes of the Eemian warming in Greenland. Simulations of the atmospheric response to combinations of Eemian insolation and preindustrial oceanic conditions and vice versa are used to disentangle the impacts of the insolation change and the related changes in sea surface temperatures and sea ice conditions. The changed oceanic conditions cause warming throughout the year, prolonging the impact of the summertime insolation increase. Consequently, the oceanic conditions cause an annual mean warming of 2 K at the NEEM site, whereas the insolation alone causes an insignificant change. Taking the precipitation changes into account, however, the insolation and oceanic changes cause more comparable increases in the precipitation-weighted temperature, implying that both contributions are important for the ice core record at the NEEM site. The simulated Eemian precipitation-weighted warming of 2.4 K at the NEEM site is low compared to the ice core reconstruction, partially due to missing feedbacks related to ice sheet changes and an extensive sea ice cover. Surface mass balance calculations with an energy balance model further indicate that the combination of temperature and precipitation anomalies leads to potential mass loss in the north and southwestern parts of the ice sheet. The oceanic conditions favor increased accumulation in the southeast, while the insolation appears to be the dominant cause of the expected ice sheet reduction. Consequently, the Eemian is not a suitable analogue for future ice sheet changes.

Mineral dust and major ion concentrations in snowpit samples from the NEEM site, Greenland

Polar ice sheets conserve atmospheric aerosols at the time of snowfall, which can be used to reconstruct past climate and environmental conditions. We investigated mineral dust and major ion records in snowpit samples obtained from the northwestern Greenland ice sheet near the North Greenland Eemian Ice Drilling (NEEM) camp in June 2009. We analyzed the samples for mineral dust concentrations as well as stable water isotopes (δ18O, δD, and deuterium excess) and major ions (Cl−, SO42-, methanesulfonic acid (MSA), Na+, and Ca2+). Seasonal δ18O and δD cycles indicate that the snowpit samples covered a six-year period from spring 2003 to early summer 2009. Concentrations of mineral dust, nss-Ca2+, and nss-SO42- showed seasonal deposition events with maxima in the winter–spring layers. On the other hand, the Cl−/Na+ ratio and the concentrations of MSA exhibited maxima in the summer layers, making them useful indicators for the summer season. Moreover, an anomalous atmospheric mineral dust event was recorded at a depth of 165–170 cm corresponding to late winter 2005 to spring 2006. A back trajectory analysis suggests that a major contributor to the Greenland aerosol was an air mass passing over the Canadian Arctic and North America. Several trajectories point to Asian regions as a dust source. The mineral dust deposited at NEEM was strongly influenced by long-range atmospheric transport and dust input from arid source areas in northern China and Mongolia.

Mathematical modeling of ice flow in the north-western Greenland and interpretation of deep drilling data at the NEEM camp

The new ice core obtained from the NEEM site in the north-western Greenland (77.449˚N, 51.056˚W, 2447 m a.s.l) in 2008–2012 was expected to improve our knowledge about the Last interglacial (also known as the Eemian, ca. 115 to 130 kyr before present). A numerical modeling, the overview of which is given in this paper, aims at assistance in interpretation of the NEEM ice core. For this purpose, an area of 400 × 400 km was delineated in northwestern Greenland as the domain for the fine-scale model at 2.5 km resolution. Modeled present-day surface velocity was in good agreement with satellite and GPS measurements. The nested ice-sheet model was run over the last two glacial cycles to reconstruct the flow history relevant for interpreting the NEEM ice core. A Lagrangian backtracing procedure was applied to determine the particle trajectories of the ice drilled at the NEEM site. This provides the places of origin at the time of deposition of NEEM ice, from which the ice chronology and non-climatic biases of the records are determined. The latter biases arise from elevation changes of the ice sheet, advection of higher upstream ice, and from latitudinal contrasts in isotopical composition of the ice. These need to be separated from the isotope records to retrieve the climatic signal. In spite of initial expectations, the segment in the NEEM ice core between 2206 m and 2435 m turned out to be heavily disturbed. The presence of high isotope values below 2206 m depth however reveals that there is ice from the previous warm interglacial. Our model locates the EMBED Equation.DSMT4 maximum of the warmest Eemian ice at the right depth (2400 m). On this basis we suggest that the measured contrast between the present and the Eemian EMBED Equation.DSMT4 must be increased by about 1.5 ‰. Most of this non-climatic bias results from upstream advection over an estimated distance of ~175 km. Except for the disturbed section, we are still confident to be able to provide accurate model-based estimates of the ice chronology and non-climatic biases for the core down to ~87% of relative depth. This requires revision the palaeo-accumulation parameterization to better match thinner annual layer thicknesses observed in the ice core for the last glacial period.

Representativeness and seasonality of major ion records derived from NEEM firn cores

Abstract. The seasonal and annual representativeness of ionic aerosol proxies (among others, calcium, sodium, ammonium and nitrate) in various firn cores in the vicinity of the NEEM drill site in northwest Greenland have been assessed. Seasonal representativeness is very high as one core explains more than 60% of the variability within the area. The inter-annual representativeness, however, can be substantially lower (depending on the species) making replicate coring indispensable to derive the atmospheric variability of aerosol species. A single core at the NEEM site records only 30% of the inter-annual atmospheric variability in some species, while five replicate cores are already needed to cover approximately 70% of the inter-annual atmospheric variability in all species. The spatial representativeness is very high within 60 cm, rapidly decorrelates within 10 m but does not diminish further within 3 km. We attribute this to wind reworking of the snow pack leading to sastrugi formation. Due to the high resolution and seasonal representativeness of the records we can derive accurate seasonalities of the measured species for modern (AD 1990–2010) times as well as for pre-industrial (AD 1623–1750) times. Sodium and calcium show similar seasonality (peaking in February and March respectively) for modern and pre-industrial times, whereas ammonium and nitrate are influenced by anthropogenic activities. Nitrate and ammonium both peak in May during modern times, whereas during pre-industrial times ammonium peaked during July–August and nitrate during June–July.

10Be climate fingerprints during the Eemian in the NEEM ice core, Greenland

Several deep Greenland ice cores have been retrieved, however, capturing the Eemian period has been problematic due to stratigraphic disturbances in the ice. The new Greenland deep ice core from the NEEM site (77.45°N, 51.06°W, 2450 m.a.s.l) recovered a relatively complete Eemian record. Here we discuss the cosmogenic 10Be isotope record from this core. The results show Eemian average 10Be concentrations about 0.7 times lower than in the Holocene which suggests a warmer climate and approximately 65–90% higher precipitation in Northern Greenland compared to today. Effects of shorter solar variations on 10Be concentration are smoothed out due to coarse time resolution, but occurrence of a solar maximum at 115.26–115.36 kyr BP is proposed. Relatively high 10Be concentrations are found in the basal ice sections of the core which may originate from the glacial-interglacial transition and relate to a geomagnetic excursion about 200 kyr BP.

Archives of total mercury reconstructed with ice and snow from Greenland and the Canadian High Arctic

This study reports total Hg concentration and atmospheric flux data from ice cores and snow/ice shallow pits from two Canadian Arctic and one Greenland glaciers, with the aim of reconstructing a high resolution record of THg deposition extending back into the pre-industrial period. An 88-m ice core (653 samples) from the NEEM glacier site in Northwest Greenland was retrieved in August 2010. The bottom sample was dated to 1748, resulting in a 262year archive. Snow and ice samples (143 samples) were recovered from a 10.3-m pit dug on the Mt. Oxford Icefield, Nunavut, in May 2008, covering 30years. Another 15.5-m short core drilled on the Agassiz Ice Cap, Nunavut, in April 2009 yielded 191 samples covering 74years. Net rates of atmospheric THg deposition (FTHg) were calculated based on THg concentrations and snow accumulation rates. Results from NEEM site show that THg and FTHg range from sub-pg g−1 to 120.6pgg−1 (mean=1.5pgg−1, n=653) and from 0.06 to 1.42μgm−2 year−1 (mean=0.25μgm−2 year−1, n=218) respectively, much lower than those found in other natural media such as sediments, peat bogs and wet precipitation. The discrepancy of FTHg found in glaciers from other natural media could mainly be due to the more severe photo-reduction and reemission of deposited oxidized Hg. This study also demonstrates that reproducible THg archives can be reconstructed with glacier ice and snow samples from Greenland and the Canadian High Arctic. The THg archive reconstructed with the short core from NEEM site is so far the longest with the highest resolution in Greenland and the Canadian High Arctic.

A first chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core

Abstract. A stratigraphy-based chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core has been derived by transferring the annual layer counted Greenland Ice Core Chronology 2005 (GICC05) and its model extension (GICC05modelext) from the NGRIP core to the NEEM core using 787 match points of mainly volcanic origin identified in the electrical conductivity measurement (ECM) and dielectrical profiling (DEP) records. Tephra horizons found in both the NEEM and NGRIP ice cores are used to test the matching based on ECM and DEP and provide five additional horizons used for the timescale transfer. A thinning function reflecting the accumulated strain along the core has been determined using a Dansgaard–Johnsen flow model and an isotope-dependent accumulation rate parameterization. Flow parameters are determined from Monte Carlo analysis constrained by the observed depth-age horizons. In order to construct a chronology for the gas phase, the ice age–gas age difference (Δage) has been reconstructed using a coupled firn densification-heat diffusion model. Temperature and accumulation inputs to the Δage model, initially derived from the water isotope proxies, have been adjusted to optimize the fit to timing constraints from δ15N of nitrogen and high-resolution methane data during the abrupt onset of Greenland interstadials. The ice and gas chronologies and the corresponding thinning function represent the first chronology for the NEEM core, named GICC05modelext-NEEM-1. Based on both the flow and firn modelling results, the accumulation history for the NEEM site has been reconstructed. Together, the timescale and accumulation reconstruction provide the necessary basis for further analysis of the records from NEEM.

Drilling a deep ice core at the NEEM site in Greenland

Drilling a deep ice core at the NEEM site in Greenland

Eemian interglacial reconstructed from a Greenland folded ice core

Efforts to extract a Greenland ice core with a complete record of the Eemian interglacial (130,000 to 115,000 years ago) have until now been unsuccessful. The response of the Greenland ice sheet to the warmer-than-present climate of the Eemian has thus remained unclear. Here we present the new North Greenland Eemian Ice Drilling ('NEEM') ice core and show only a modest ice-sheet response to the strong warming in the early Eemian. We reconstructed the Eemian record from folded ice using globally homogeneous parameters known from dated Greenland and Antarctic ice-core records. On the basis of water stable isotopes, NEEM surface temperatures after the onset of the Eemian (126,000 years ago) peaked at 8 ± 4 degrees Celsius above the mean of the past millennium, followed by a gradual cooling that was probably driven by the decreasing summer insolation. Between 128,000 and 122,000 years ago, the thickness of the northwest Greenland ice sheet decreased by 400 ± 250 metres, reaching surface elevations 122,000 years ago of 130 ± 300 metres lower than the present. Extensive surface melt occurred at the NEEM site during the Eemian, a phenomenon witnessed when melt layers formed again at NEEM during the exceptional heat of July 2012. With additional warming, surface melt might become more common in the future.

Gas transport in firn: multiple-tracer characterisation and model intercomparison for NEEM, Northern Greenland

Abstract. Air was sampled from the porous firn layer at the NEEM site in Northern Greenland. We use an ensemble of ten reference tracers of known atmospheric history to characterise the transport properties of the site. By analysing uncertainties in both data and the reference gas atmospheric histories, we can objectively assign weights to each of the gases used for the depth-diffusivity reconstruction. We define an objective root mean square criterion that is minimised in the model tuning procedure. Each tracer constrains the firn profile differently through its unique atmospheric history and free air diffusivity, making our multiple-tracer characterisation method a clear improvement over the commonly used single-tracer tuning. Six firn air transport models are tuned to the NEEM site; all models successfully reproduce the data within a 1σ Gaussian distribution. A comparison between two replicate boreholes drilled 64 m apart shows differences in measured mixing ratio profiles that exceed the experimental error. We find evidence that diffusivity does not vanish completely in the lock-in zone, as is commonly assumed. The ice age- gas age difference (Δage) at the firn-ice transition is calculated to be 182+3−9 yr. We further present the first intercomparison study of firn air models, where we introduce diagnostic scenarios designed to probe specific aspects of the model physics. Our results show that there are major differences in the way the models handle advective transport. Furthermore, diffusive fractionation of isotopes in the firn is poorly constrained by the models, which has consequences for attempts to reconstruct the isotopic composition of trace gases back in time using firn air and ice core records.