ALBERT

Description: Correct estimate of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice cores studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: output of a firn densification model and measurements of δ15N of N2 in air trapped in ice core. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available δ15N measurements performed from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rate and temperature conditions. While firn densification simulations are able to correctly represent most of the δ15N trends over the last deglaciation measured in the EDC, BI, TALDICE and EDML ice cores, they systematically fail to capture BI and EDML δ15N glacial levels, a mismatch previously seen for Central East Antarctic ice cores. Using empirical constraints of the EDML gas-ice depth offset during the Laschamp event (~ 41 ka), we can rule out the existence of a large convective zone as the explanation of the glacial firn model-δ15N data mismatch for this site. The good match between modelled and measured δ15N at TALDICE as well as the lack of any clear correlation between insoluble dust concentration in snow and δ15N records in the different ice cores suggest that past changes in loads of impurities are not the only main driver of glacial-interglacial changes in firn lock-in depth. We conclude that firn densification dynamics may instead be driven mostly by accumulation rate changes. The mismatch between modelled and measured δ15N may be due to inaccurate reconstruction of past accumulation rate or underestimated influence of accumulation rate in firnification models.