Hydrogeochemical transport modeling of the infiltration of tertiary treated wastewater in a dune area, Belgium

Managed artificial recharge (MAR) is a well-established practice for augmentation of depleted groundwater resources or for environmental benefit. At the St-André MAR site in the Belgian dune area, groundwater resources are optimised through re-use of highly treated wastewater by means of infiltration ponds. The very high quality of the infiltration water sets this system apart from other MAR systems. The low total dissolved solid (TDS) content in the infiltration water (less than 50 mg/L) compared to the dune aquifer (500 mg/L) triggers a number of reactions, increasing the TDS through soil-aquifer passage. Multi-component reactive transport modelling was applied to analyse the geochemical processes that occur. Carbonate dissolution is the main process increasing the TDS of the infiltration water. Oxic aquifer conditions prevail between the infiltration ponds and the extraction wells. This is driven by the high flow velocities, leaving no time to consume O₂ between the ponds and extraction wells. Cation exchange is important when infiltration water is replaced by native dune water or when significant changes in infiltration-water quality occur. The seasonal variation of O₂ and temperature in the infiltration water are the main drivers for seasonal changes in the concentration of all major ions.     

Source, attenuation and potential mobility of arsenic at New Britannia Mine, Snow Lake, Manitoba

This study was to investigate the source, mobility and attenuation of As at the New Britannia Mine, Snow Lake, Manitoba. One major source of As contamination was determined to be an arsenopyrite residue stockpile (ARS) containing refractory Au in a waste rock impoundment. It appears that As is still moving through glacial clay at the base of the ARS into a confined aquifer even though the pile was capped in the year 2000. Arsenic is also being mobilized from a deposit of tailings, which formed following spills by previous owners, Nor Acme. Arsenic from the tailings is being mobilized by oxidation of arsenopyrite and reduction of arsenate to the more mobile arsenite by arsenate-reducing bacteria. This contamination is affecting a shallow unconfined aquifer and surface water flowing from the tailings through wetlands towards Snow Lake. Arsenic is being attenuated by adsorption to hydrated ferric oxides (HFO) in the tailings, wetland soils and aquatic plants. Although As in surface water, soils and plants along the flow path from the mine to Snow Lake are above Canadian drinking water guidelines, efficient natural attenuation by HFO in soils and plants of the wetlands have limited the concentration in Snow Lake to below drinking water standards.     

The magnesium and calcium contents of sediments,especially pelites,as a function of age and degree of metamorphism

The silicate and carbonate fraction of 98 non-metamorphic shale samples from the Australian platform and of different geological age were analysed for calcium, magnesium, ferrous iron and carbonate. Cainozoic and Mesozoic shales prove to be essentially calcitic, Cambrian and Proterozoic shales are essentially dolomitic and sideritic. A similar trend of high MgO values can be demonstrated for the silicate fraction of the old shales. Extensive literature study confirms these trends for shales and carbonate rocks from all over the world. Slates, hornfelses and schists are Mg rich and Ca poor, whether young or old.Ronov's model of the evolution of the earth's crust ocean and atmosphere, explaining these trends, is critically reviewed but rejected because of impossible storage problems of calcium in the Proterozoic. The increased magnesium content of the old sediments is explained by calcium carbonate sweating out of the sedimentary column, magnesium introduction from altering volcanic rocks within the sedimentary pile and magnesium introduction from connate brines in sandstones. The increasing calcium content of all kinds of sediments with decreasing age is claimed to be related to preferential weathering of extrusive volcanic rocks and sweating out of calcium carbonate from the sedimentary column.     

Evaluating the impact of soil redistribution on the in situ mineralization of soil organic carbon

Soil erosion, transport and deposition by water drastically affect the distribution of soil organic carbon (SOC) within a landscape. Moreover, soil redistribution may have a large impact on the exchange of carbon (C) between the pedosphere and the atmosphere. One of the large information gaps within this research domain, concerns the fate of SOC after erosion by water. According to different (mainly laboratory) studies, soil redistribution leads to aggregate breakdown, thereby exposing the contained SOC to mineralization. Our study aims to quantify the extent to which such increased mineralization occurs in a real field situation. Carbon dioxide (CO₂)‐efflux was measured in the field after an important erosion event for a continuous period of 112 days. The specific situation on the field ensured that almost none of eroded SOC was exported from the field. Measurements of CO₂‐efflux were done in areas with sediment deposition, as well as in comparable areas without sedimentation. Comparison of these measurements allowed the net effect of soil deposition on CO₂‐efflux to be assessed. Field data were complemented by measurements on incubated, undisturbed soil core samples, in order to disentangle the contribution of environmental factors (moisture, temperature) from any erosional effect on CO₂‐efflux. Results of these measurements on the field showed that CO₂‐efflux was regulated by a complex interplay of different factors (mostly soil porosity, soil moisture and soil temperature). In combination with the incubation measurements, it could be concluded that the processes of erosion and transport indeed led to an increased mineralization of SOC, as a result of aggregate breakdown and exposure of previously encapsulated SOC. This effect was, however, much smaller than observed in previous laboratory studies. Moreover, it was only important in the first weeks, immediately after the erosion event. The calculated net erosional effect on CO₂‐efflux represented a mere 1·6% of total SOC, originally present in the soil. Copyright © 2010 John Wiley & Sons, Ltd.     

A surface structural model for ferrihydrite II: Adsorption of uranyl and carbonate

The adsorption of uranyl (UO₂²⁺) on ferrihydrite has been evaluated with the charge distribution (CD) model for systems covering a very large range of conditions, i.e. pH, ionic strength, CO₂ pressure, U(VI) concentration, and loading. Modeling suggests that uranyl forms bidentate inner sphere complexes at sites that do not react chemically with carbonate ions. Uranyl is bound by singly-coordinated surface groups present at particular edges of Fe-octahedra of ferrihydrite while another set of singly-coordinated surface groups may form double-corner bidentate complexes with carbonate ions. The uranyl surface speciation strongly changes in the presence of carbonate due to the specific adsorption of carbonate ions as well as the formation of ternary uranyl-carbonate surface complexes. Data analysis with the CD model suggests that a uranyl tris-carbonato surface complex, i.e. (UO₂)(CO₃)₃⁴⁻, is formed. This species is most abundant in systems with a high pH and carbonate concentration. This finding differs significantly from previous interpretations made in the literature. At high pH and low carbonate concentrations, as can be prepared in CO₂-closed systems, the model suggests the additional presence of a ternary uranyl-monocarbonato complex. The binding mode (type A or type B complex) is uncertain. At high uranyl concentrations, uranyl polymerizes at the surface of ferrihydrite giving, for instance, tris-uranyl surface complexes with and without carbonate. The similarities and differences between U(VI) adsorption by goethite and ferrihydrite are discussed from a surface structural point of view.     

A synthetic Earth Gravity Model Designed Specifically for Testing Regional Gravimetric Geoid Determination Algorithms

A synthetic [simulated] Earth gravity model (SEGM) of the geoid, gravity and topography has been constructed over Australia specifically for validating regional gravimetric geoid determination theories, techniques and computer software. This regional high-resolution (1-arc-min by 1-arc-min) Australian SEGM (AusSEGM) is a combined source and effect model. The long-wavelength effect part (up to and including spherical harmonic degree and order 360) is taken from an assumed errorless EGM96 global geopotential model. Using forward modelling via numerical Newtonian integration, the short-wavelength source part is computed from a high-resolution (3-arc-sec by 3-arc-sec) synthetic digital elevation model (SDEM), which is a fractal surface based on the GLOBE v1 DEM. All topographic masses are modelled with a constant mass-density of 2,670 kg/m³. Based on these input data, gravity values on the synthetic topography (on a grid and at arbitrarily distributed discrete points) and consistent geoidal heights at regular 1-arc-min geographical grid nodes have been computed. The precision of the synthetic gravity and geoid data (after a first iteration) is estimated to be better than 30 μ Gal and 3 mm, respectively, which reduces to 1 μ Gal and 1 mm after a second iteration. The second iteration accounts for the changes in the geoid due to the superposed synthetic topographic mass distribution. The first iteration of AusSEGM is compared with Australian gravity and GPS-levelling data to verify that it gives a realistic representation of the Earth’s gravity field. As a by-product of this comparison, AusSEGM gives further evidence of the north–south-trending error in the Australian Height Datum. The freely available AusSEGM-derived gravity and SDEM data, included as Electronic Supplementary Material (ESM) with this paper, can be used to compute a geoid model that, if correct, will agree to in 3 mm with the AusSEGM geoidal heights, thus offering independent verification of theories and numerical techniques used for regional geoid modelling.Electronic Supplementary Material Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s00190-005-0002-z     

A dynamical solution of the triple asteroid system (45) Eugenia

We present the first dynamical solution of the triple asteroid system (45) Eugenia and its two moons Petit–Prince (diameter ∼ 7 km) and S/2004 (45) 1 (diameter ∼ 5 km). The two moons orbit at 1165 and 610 km from the primary, describing an almost-circular orbit (e ∼ 6 × 10⁻³ and e ∼ 7 × 10⁻² respectively). The system is quite different from the other known triple systems in the main belt since the inclinations of the moon orbits are sizeable (9° and 18° with respect to the equator of the primary respectively). No resonances, neither secular nor due to Lidov–Kozai mechanism, were detected in our dynamical solution, suggesting that these inclinations are not due to excitation modes between the primary and the moons. A 10-year evolution study shows that the orbits are slightly affected by perturbations from the Sun, and to a lesser extent by mutual interactions between the moons. The estimated J₂ of the primary is three times lower than the theoretical one, calculated assuming the shape of the primary and an homogeneous interior, possibly suggesting the importance of other gravitational harmonics.     

Topography and landslides in weathered granitic rock areas—Hai Van landslide in central Vietnam

Many landslides occur every year during heavy rains at the Hai Van Pass and surrounding area in central Vietnam, where granitic rocks are distributed. As is common in granite areas, these landslides often occur as small-scale flow-type and slump-type landslides. However, several horseshoe-shaped loose slopes of widths and lengths of 500 to 800 m, which incorporate these landslides, are observed on slopes across the area. These topographies resemble those formed by past and present large-scale landslides. The presence of such a topography and the repeated occurrences of landslides within this topography are rare in granite areas, where shallow flow-type landslides are generally frequent. To understand the mechanism causing the landslides in the Hai Van region, and as a support for future risk assessment, the factors and processes leading to the formation of such a topography and their relationship with these landslides must be identified and assessed. This study investigated the history of past landslide movement in the Hai Van Pass and surrounding area through observations of drill cores and outcrops, and analysis of the direction of remanent magnetism in the granitic rocks. Mineral compositions, cracks, degrees of weathering, and topographic shapes of the granitic rocks and their relationship to the landslides occurring today were also investigated. The results of the study reveal no variation in the direction of remanent magnetism in the granitic rocks in the region that would indicate disturbance of the ground due to a past large-scale landslide. No evidence of such an event could also be found both in the drill cores and the rock outcrops. Further, results of the analysis of cracks and weathering pattern confirm that the topography of the region is affected by the weathering of the granitic rocks that progresses in concentric circles of various sizes. Thus, it can be concluded that these topographies were not formed by a singular large-scale landslide of the past, but rather by a composite of relatively shallow landslides occurring on the slope of dome structures unique to granite areas, which are formed by differential weathering and denudation regulated by cracks.     

Review paper: Virtual sources and their responses,Part II: data‐driven single‐sided focusing

In Part I of this paper, we defined a focusing wave field as the time reversal of an observed point‐source response. We showed that emitting a time‐reversed field from a closed boundary yields a focal spot that acts as an isotropic virtual source. However, when emitting the field from an open boundary, the virtual source is highly directional and significant artefacts occur related to multiple scattering. The aim of this paper is to discuss a focusing wave field, which, when emitted into the medium from an open boundary, yields an isotropic virtual source and does not give rise to artefacts. We start the discussion from a horizontally layered medium and introduce the single‐sided focusing wave field in an intuitive way as an inverse filter. Next, we discuss single‐sided focusing in two‐dimensional and three‐dimensional inhomogeneous media and support the discussion with mathematical derivations. The focusing functions needed for single‐sided focusing can be retrieved from the single‐sided reflection response and an estimate of the direct arrivals between the focal point and the accessible boundary. The focal spot, obtained with this single‐sided data‐driven focusing method, acts as an isotropic virtual source, similar to that obtained by emitting a time‐reversed point‐source response from a closed boundary.