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 1s 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

Methane is an important greenhouse gas that is emitted from multiple natural and anthropogenic sources. Atmospheric methane concentrations have varied on a number of timescales in the past, but what has caused these variations is not always well understood. The different sources and sinks of methane have specific isotopic signatures, and the isotopic composition of methane can therefore help to identify the environmental drivers of variations in atmospheric methane concentrations. Here we present high-resolution carbon isotope data (δ(13)C content) for methane from two ice cores from Greenland for the past two millennia. We find that the δ(13)C content underwent pronounced centennial-scale variations between 100 BC and AD 1600. With the help of two-box model calculations, we show that the centennial-scale variations in isotope ratios can be attributed to changes in pyrogenic and biogenic sources. We find correlations between these source changes and both natural climate variability--such as the Medieval Climate Anomaly and the Little Ice Age--and changes in human population and land use, such as the decline of the Roman empire and the Han dynasty, and the population expansion during the medieval period.