Simulation of DNAPL infiltration and spreading behaviour in the saturated zone at varying flow velocities and alternating subsurface geometries

The influence of varying groundwater flow velocities on DNAPL infiltration and spreading behaviour was investigated by multiphase modelling using TMVOC and PetraSim. The multiphase models were calibrated by results of previously conducted laboratory experiments for the complete spatio-temporal range of the experiments. The small scale 2D scenario modelling was applied to qualify and quantify changes in position, architecture, geometry and dissolution of a TCE body in a fully saturated homogeneous sandy medium. The applied flow velocities ranging from 0.05 up to 40.00 m/day exhibited that the DNAPL TCE is affected even at the lowest flow velocity in its position, its size and its architecture. Additionally, several impermeable lenses with simple geometry were assumed in the model, to investigate the influence of stratified subsoil. In the experimental set-ups, the DNAPL body reacts more sensitive to the applied groundwater flow velocities than to the geometrical set-up of the scenarios. A possible consequence can be the transportation and displacement of a DNAPL pool due to natural or anthropogenic induced high groundwater flow velocities, as by Pump and Treat facilities, complicating site investigation process and planning of remediation activities.     

On the meridional circulation and balance of momentum in the Southern Ocean of POP

 The circulation of the Southern Ocean is studied in the eddy-resolving model POP (Parallel Ocean Program) by an analysis of zonally integrated balances. The TEM formalism (Transformed Eulerian Mean) is extended to include topography and continental boundaries, thus deviations from a zonally integrated state involve transient and standing eddies. The meridional circulation is presented in terms of the Eulerian, eddy-induced, and residual streamfunctions. It is shown that the splitting of the meridional circulation into Ekman and geostrophic transports and the component induced by subgrid and Reynolds stresses is identical to a particular form of the zonally integrated balance of zonal momentum. In this balance, the eddy-induced streamfunctions represent the interfacial form stresses by transient and standing eddies and the residual streamfunction represents the acceleration of the zonal current by density fluxes in a zonally integrated frame. The latter acceleration term is directly related to the surface flux of density and interior fluxes due to the resolved and unresolved eddies. The eddy-induced circulation is extremely vigorous in POP. In the upper ocean a shallow circulation, reversed in comparison to the Deacon cell and mainly due to standing eddies, appears to the north of Drake Passage latitudes, and in the Drake Passage belt of latitudes a deep-reaching cell is induced by transient eddies. In the resulting residual circulation the Deacon cell is largely cancelled and the residual advection of the zonal mean potential density is balanced by diapycnal eddy and subgrid fluxes which are strong in the upper few hundred meters but small in the ocean interior. The balance of zonal momentum is consistent with other eddy-resolving models; a new aspect is the clear identification of density effects in the zonally integrated balance. We show that the wind stress and the stress induced by the residual circulation drive the eastward current, whereas both eddy species result in a braking. Finally, we extend the Johnson–Bryden model of zonal transport to incorporate all relevant terms from the zonal momentum balance. It is shown that wind stress and induction by the residual circulation carry an eastward transport while bottom form stress and the stress induced by standing eddies yield westward components of transport. Received: 26 June 2001 / Accepted: 2 November 2001     

The Summer Surface Energy Balance of the High Antarctic Plateau

The summertime surface energy balance (SEB) at Kohnen station, situated on the high Antarctic plateau (75°00′ S, 0°04′ E, 2892m above sea level) is presented for the period of 8 January to 9 February 2002. Shortwave and longwave radiation fluxes were measured directly; the former was corrected for problems associated with the cosine response of the instrument. Sensible and latent heat fluxes were calculated using the bulk method, and eddy-correlation measurements and the modified Bowen ratio method were used to verify these calculated fluxes. The calculated sub-surface heat flux was checked by comparing calculated to measured snow temperatures. Uncertainties in the measurements and energy-balance calculations are discussed. The general meteorological conditions were not extraordinary during the period of the experiment, with a mean 2-m air temperature of −27.5°C, specific humidity of 0.52×10⁻³kg kg⁻¹ and wind speed of 4.1ms⁻¹. The experiment covered the transition period from Antarctic summer (positive net radiation) to winter (negative net radiation), and as a result the period mean net radiation, sensible heat, latent heat and sub-surface heat fluxes were small with values of −1.1, 0.0, −1.0 and 0.7 Wm⁻², respectively. Daily mean net radiation peaked on cloudy days (16 Wm⁻²) and was negative on clear-sky days (minimum of −19 W m⁻²). Daily mean sensible heat flux ranged from −8 to +10 Wm⁻², latent heat flux from −4 to 0 Wm⁻² and sub-surface heat flux from −8 to +7 Wm⁻².     

The Deformational Response of a Viscoelastic Solid Earth Model Coupled to a Thermomechanical Ice Sheet Model

We apply a coupled thermomechanical ice sheet—self-gravitating viscoelastic solid Earth model (SGVEM), allowing for the dynamic exchange of ice thickness and bedrock deformation, in order to investigate the effect of viscoelastic deformation on ice dynamics and vice versa. In a synthetic glaciation scenario, we investigate the interaction between the ice sheet and the solid Earth deformation, the glacial-isostatic adjustment (GIA), accounting for an atmospheric forcing depending on the ice sheet surface altitude. We compare the results from the coupled model to runs with the common elastic lithosphere/relaxing asthenosphere (ELRA) model, where the lithosphere is represented by a thin plate and the mantle relaxes with one characteristic relaxation time, as well as to a rigid Earth without any deformation. We find that the deformational behaviour of the SGVEM on ice dynamics (i.e. stored ice volume, ice thickness and velocity field) is comparable to the ELRA for an optimal choice of the parameters in steady state, but exhibits differences in the transient behaviour. Beyond the ice sheet, in the region of peripheral forebulge, the differences in the transient surface deformation between ELRA and SGVEM are substantial, demonstrating the inadequacy of the ELRA model for interpreting constraints on GIA in the periphery of the ice sheet, such as sea-level indicators and GPS uplift rates.     

Resistance of Lophelia pertusa to coverage by sediment and petroleum drill cuttings

In laboratory experiments, the cold-water coral Lophelia pertusa was exposed to settling particles. The effects of reef sediment, petroleum drill cuttings and a mix of both, on the development of anoxia at the coral surface were studied using O₂, pH and H₂S microsensors and by assessing coral polyp mortality. Due to the branching morphology of L. pertusa and the release of coral mucus, accumulation rates of settling material on coral branches were low. Microsensors detected H₂S production in only a few samples, and sulfate reduction rates of natural reef sediment slurries were low (<0.3 nmol S cm⁻³ d⁻¹). While the exposure to sediment clearly reduced the coral’s accessibility to oxygen, L. pertusa tolerated both partial low-oxygen and anoxic conditions without any visible detrimental short-term effect, such as tissue damage or death. However, complete burial of coral branches for >24 h in reef sediment resulted in suffocation.     

Zircon M257 ‐ a Homogeneous Natural Reference Material for the Ion Microprobe U‐Pb Analysis of Zircon

We introduce and propose zircon M257 as a future reference material for the determination of zircon U‐Pb ages by means of secondary ion mass spectrometry. This light brownish, flawless, cut gemstone specimen from Sri Lanka weighed 5.14 g (25.7 carats). Zircon M257 has TIMS‐determined, mean isotopic ratios (2s uncertainties) of 0.09100 ± 0.00003 for ²⁰⁶pb/²³⁸U and 0.7392 ± 0.0003 for ²⁰⁷pb/²³⁵U. Its ²⁰⁶pb/²³⁸U age is 561.3 ± 0.3 Ma (unweighted mean, uncertainty quoted at the 95% confidence level); the U‐Pb system is concordant within uncertainty of decay constants. Zircon M257 contains ～ 840 μg g^?1 U (Th/U ～ 0.27). The material exhibits remarkably low heterogeneity, with a virtual absence of any internal textures even in cathodoluminescence images. The uniform, moderate degree of radiation damage (estimated from the expansion of unit‐cell parameters, broadening of Raman spectral parameters and density) corresponds well, within the “Sri Lankan trends”, with actinide concentrations, U‐Pb age, and the calculated alpha fluence of 1.66 × 10¹⁸ g^?1. This, and a (U+Th)/He age of 419 ± 9 Ma (2s), enables us to exclude any unusual thermal history or heat treatment, which could potentially have affected the retention of radiogenic Pb. The oxygen isotope ratio of this zircon is 13.9%o VSMOW suggesting a metamorphic genesis in a marble or calc‐silicate skarn.     

Colloidal arsenic composition from abandoned gold mine tailing leachates in Nova Scotia,Canada

Asymmetric-Flow Field-Flow Fractionation (AsFlFFF) coupled to an inductively coupled plasma-mass spectrometer (ICP-MS) was used to determine whether colloidal As exists in mine tailings from abandoned Au mine sites in Nova Scotia. Using this technique, the chemical composition and the size of the colloidal material was determined. Tailings samples were collected from the sites and leachates were analyzed. The resulting fractograms showed the presence of colloidal As. Arsenic co-elutes with Fe and Al suggesting that the As is associated with colloids containing these elements. The colloidal material present in the samples has a calculated median size of 7 nm. The leachates were also analyzed for totalAs, Fe and Al using ICP-MS. The colloidal fraction was determined to account for approximately 20% of the total dissolved As present in these samples. Dissolved As concentrations at the selected mine sites are very high and the arbitrary cutoff for dissolved still incorporates colloidal material. Therefore. it is important to distinguish between truly dissolved and colloidal As as the speciation will affect the toxicity and mobility of As at these locations.