DEP-derived density from two snow trenches from Kohnen Station, Dronning Maud Land, Antarctica from the 2012/13 field season

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.

Dissolution of a spatially variable DNAPL (dense non aqueous phase liquid) source

Dissolution of non-aqueous phase liquids (NAPLs) or gases into groundwater is a key process, both for contamination problems originating from organic liquid sources, and for dissolution trapping in geological storage of CO2. Dissolution in natural systems typically will involve both high and low NAPL saturations and a wide range of pore water flow velocities within the same source zone for dissolution to groundwater. To correctly predict dissolution in such complex systems and as the NAPL saturations change over time, models must be capable of predicting dissolution under a range of saturations and flow conditions. To provide data to test and validate such models, an experiment was conducted in a two-dimensional sand tank, where the dissolution of a spatially variable, 5x5 cm**2 DNAPL tetrachloroethene source was carefully measured using x-ray attenuation techniques at a resolution of 0.2x0.2 cm**2. By continuously measuring the NAPL saturations, the temporal evolution of DNAPL mass loss by dissolution to groundwater could be measured at each pixel. Next, a general dissolution and solute transport code was written and several published rate-limited (RL) dissolution models and a local equilibrium (LE) approach were tested against the experimental data. It was found that none of the models could adequately predict the observed dissolution pattern, particularly in the zones of higher NAPL saturation. Combining these models with a model for NAPL pool dissolution produced qualitatively better agreement with experimental data, but the total matching error was not significantly improved. A sensitivity study of commonly used fitting parameters further showed that several combinations of these parameters could produce equally good fits to the experimental observations. The results indicate that common empirical model formulations for RL dissolution may be inadequate in complex, variable saturation NAPL source zones, and that further model developments and testing is desirable.