Semi-analytical solutions of a contaminant transport equation with nonlinear sorption in 1D

A new method to determine semi-analytical solutions of one-dimensional contaminant transport problem with nonlinear sorption is described. This method is based on operator splitting approach where the convective transport is solved exactly and the diffusive transport by finite volume method. The exact solutions for all sorption isotherms of Freundlich and Langmuir type are presented for the case of piecewise constant initial profile and zero diffusion. Very precise numerical results for transport with small diffusion can be obtained even for larger time steps (e.g., when the Courant-Friedrichs-Lewy (CFL) condition failed).     

Assessment of gully erosion rates through interviews and measurements: a case study from northern Ethiopia

Gullying has been widespread in the Ethiopian Highlands during the 20th century. It threatens the soil resource, lowers crop yields in intergully areas through enhanced drainage and desiccation, and aggravates flooding and reservoir siltation. Knowing the age and rates of gully development during the last few decades will help explain the reasons for current land degradation. In the absence of historical written or photographic documentation, the AGERTIM method (Assessment of Gully Erosion Rates Through Interviews and Measurements) has been developed. It comprises measurements of contemporary gully volumes, monitoring of gully evolution over several years and semi‐structured interview techniques. Gully erosion rates in the Dogu'a Tembien District, Tigray, Ethiopia, were estimated in three representative case‐study areas. In Dingilet, gullying started around 1965 after gradual environmental changes (removal of vegetation from cropland in the catchment and eucalyptus plantation in the valley bottom); rill‐like incisions grew into a gully, which increased rapidly in the drier period between 1977 and 1990. The estimated evolution of the total gully volume in the other areas show patterns similar to those of the Dingilet gully. Average gully erosion rate over the last 50 years is 6·2 t ha^?1 a^?1. Since 1995, no new gullies have developed in the study area. Area‐specific short‐term gully erosion rates are now on average 1·1 t ha^?1 a^?1. The successful application of the AGERTIM method requires an understanding of the geomorphology of the study area and an integration of the researchers with the rural society. It reveals that rapid gully development in the study area is some 50 years old and is mainly caused by human‐induced environmental degradation. Under the present‐day conditions of ‘normal’ rain and catchment‐wide soil and water conservation, gully erosion rates are decreasing. Copyright © 2006 John Wiley & Sons, Ltd.     

Thermal neutron capture effects in radioactive and stable nuclide systems

Neutron capture effects in meteorites and lunar surface samples have been successfully used in the past to study exposure histories and shielding conditions. In recent years, however, it turned out that neutron capture effects produce a nuisance for some of the short‐lived radionuclide systems. The most prominent example is the ¹⁸²Hf‐¹⁸²W system in iron meteorites, for which neutron capture effects lower the ¹⁸²W/¹⁸⁴W ratio, thereby producing too old apparent ages. Here, we present a thorough study of neutron capture effects in iron meteorites, ordinary chondrites, and carbonaceous chondrites, whereas the focus is on iron meteorites. We study in detail the effects responsible for neutron production, neutron transport, and neutron slowing down and find that neutron capture in all studied meteorite types is not, as usually expected, exclusively via thermal neutrons. In contrast, most of the neutron capture in iron meteorites is in the epithermal energy range and there is a significant contribution from epithermal neutron capture even in stony meteorites. Using sophisticated particle spectra and evaluated cross section data files for neutron capture reactions we calculate the neutron capture effects for Sm, Gd, Cd, Pd, Pt, and Os isotopes, which all can serve as neutron‐dose proxies, either in stony or in iron meteorites. In addition, we model neutron capture effects in W and Ag isotopes. For W isotopes, the GCR‐induced shifts perfectly correlate with Os and Pt isotope shifts, which therefore can be used as neutron‐dose proxies and permit a reliable correction. We also found that GCR‐induced effects for the ¹⁰⁷Pd‐¹⁰⁷Ag system can be significant and need to be corrected, a result that is in contrast to earlier studies.     

Dust ejection from planetary bodies by temperature gradients: Laboratory experiments

Laboratory experiments show that dusty bodies in a gaseous environment eject dust particles if they are illuminated. We find that even more intense dust eruptions occur when the light source is turned off. We attribute this to a compression of gas by thermal creep in response to the changing temperature gradients in the top dust layers. The effect is studied at a light flux of 13 kW/m² and 1 mbar ambient pressure. The effect is applicable to protoplanetary disks and Mars. In the inner part of protoplanetary disks, planetesimals can be eroded especially at the terminator of a rotating body. This leads to the production of dust which can then be transported towards the disk edge or the outer disk regions. The generated dust might constitute a significant fraction of the warm dust observed in extrasolar protoplanetary disks. We estimate erosion rates of about 1 kg s^?1 for 100 m parent bodies. The dust might also contribute to subsequent planetary growth in different locations or on existing protoplanets which are large enough not to be susceptible to particle loss by light induced ejection. Due to the ejections, planetesimals and smaller bodies will be accelerated or decelerated and drift outward or inward, respectively. The effect might also explain the entrainment of dust in dust devils on Mars, especially at high altitudes where gas drag alone might not be sufficient.     

Seven Asteroids Studied from Modra Observatory in the Course of Binary Asteroid Photometric Campaign

We report lightcurves for five asteroids from the inner main belt—(1703) Barry, (2590) Mourão, (4022) Nonna, (5171) Augustesen, (23031) 1999 XX₇—and for two near-Earth asteroids—(100004) 1983 VA, (144922) 2005 CK₃₈—obtained at Modra Observatory as a part of binary asteroid photometric campaign. Estimated synodic rotational periods and amplitudes of their composite lightcurves were as follows: 107.1 ± 0.5 h, 0.5 mag for (1703) Barry; 15.59 ± 0.01 h, 0.49 mag for (2590) Mourão; 2.5877 ± 0.0005 h, 0.08 mag for (4022) Nonna; 474 ± 10 h, 0.8 mag for (5171) Augustesen; 3.075 ± 0.001 h, 0.43 mag for (23031) 1999 XX₇; 3.1643 ± 0.0009 h, 0.11 mag for (100004) 1983 VA; 4.7894 ± 0.0005 h, 0.27 mag for (144922) 2005 CK₃₈. The slow rotator (5171) Augustesen and possibly also (1703) Barry appear to be tumblers, though their precession periods can not be estimated from the available data.     

Use of unit response functions for management of regional multilayered aquifers: application to the North Aquitaine Tertiary system (France)

Although there are many numerical applications used in aquifer management, few methodologies seem to be sufficiently adapted to provide solutions for regional-scale problems. Their applications for multilayered aquifers are also too limited. A new groundwater management tool adapted to such large physical systems has been developed, using the unit response function (URF) approach. This tool is based on the link between a numerical groundwater flow model and a mathematical optimisation tool. Pre- and post processing are based on a geographical information system (GIS). The developed tool has been applied to a real management case—the regional Aquitaine multilayered aquifer (France). A representative hydrodynamic model was built using MODFLOW 2000 code. First tests show that the linearity condition needed for the application of URF methodology is respected despite the high complexity of the system. Two management scenarios were tested. Optimal solutions thus obtained allow for viable exploitation alternatives to be proposed, which integrate complex environmental constraints.

---

Résumé Malgré un nombre important d’applications existant dans le domaine de la gestion des aquifères, peu de méthodologies semblent adaptées pour proposer des solutions à des problèmes réels de dimension régionale. Leurs application pour des systèmes multicouches sont souvent trop limitées. Un outil de gestion adapté à ces systèmes de large taille a été développé, basé sur la théorie des Fonctions de Réponses Unitaires (FRU). Cet outil est basé sur le couplage d’un modèle numérique d’écoulement et d’un code d’optimisation mathématique. Les entrées et sorties du modèle sont réalisées à l’aide d’un Système d’Information Géographique (SIG). L’outil développé est appliqué à un cas de gestion réel—le système aquifère régional multicouche aquitain (France). Un modèle hydrodynamique représentatif a été construit, basé sur le code MODFLOW 2000. Les premiers tests montrent que les conditions de linéarité nécessaires à l’application de la méthode FRU sont respectées, malgré l’importante complexité du système. Deux scénarii de gestion sont testés. Les solutions optimales ainsi obtenues permettent de proposer des solutions d’exploitation alternatives viables, intégrant des contraintes environnementales complexes.

---

Resumen Aunque hay muchas aplicaciones numéricas usadas en la gestión de acuíferos, unas pocas metodologías parecen estar adaptadas suficientemente para dar las soluciones a los problemas de escala regional. También son limitadas sus aplicaciones para los acuíferos multicapa. Se ha desarrollado una nueva herramienta para el manejo del agua subterránea, adaptada a tales sistemas físicos grandes, usando el avance de la Función de Respuesta Unitaria (URF). Esta herramienta se basa en la combinación entre un modelo numérico de flujo de agua subterránea y una herramienta de optimización matemática. El pre-proceso y el post-proceso se basan en un Sistema de Información Geográfico (GIS). La herramienta desarrollada se ha aplicado a un caso real de manejo en el acuífero regional multicapa de Aquitania (Francia). Se construyó un modelo hidrodinámico representativo, que usó el código MODFLOW 2000. Las primeras pruebas muestran que la condición de linearidad, necesaria para la aplicación de la metodología URF, se respeta a pesar de la complejidad alta del sistema. Se probaron dos escenarios de manejo. Las soluciones óptimas así obtenidas permiten la proyección de alternativas viables de explotación, las cuales integran restricciones medioambientales complejas.

     

Spatial Modeling of Geometallurgical Properties: Techniques and a Case Study

High-resolution spatial numerical models of metallurgical properties constrained by geological controls and more extensively by measured grade and geomechanical properties constitute an important part of geometallurgy. Geostatistical and other numerical techniques are adapted and developed to construct these high-resolution models accounting for all available data. Important issues that must be addressed include unequal sampling of the metallurgical properties versus grade assays, measurements at different scale, and complex nonlinear averaging of many metallurgical parameters. This paper establishes techniques to address each of these issues with the required implementation details and also demonstrates geometallurgical mineral deposit characterization for a copper–molybdenum deposit in South America. High-resolution models of grades and comminution indices are constructed, checked, and are rigorously validated. The workflow demonstrated in this case study is applicable to many other deposit types.     

Stable Cu and Zn isotope ratios as tracers of sources and transport of Cu and Zn in contaminated soil

Copper and Zn metals are produced in large quantities for different applications. During Cu production, large amounts of Cu and Zn can be released to the environment. Therefore, the surroundings of Cu smelters are frequently metal-polluted. We determined Cu and Zn concentrations and Cu and Zn stable isotope ratios (δ⁶⁵Cu, δ⁶⁶Zn) in three soils at distances of 1.1, 3.8, and 5.3 km from a Slovak Cu smelter and in smelter wastes (slag, sludge, ash) to trace sources and transport of Cu and Zn in soils. Stable isotope ratios were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICP-MS) in total digests. Soils were heavily contaminated with concentrations up to 8087 μg g⁻¹ Cu and 2084 μg g⁻¹ Zn in the organic horizons. The δ⁶⁵Cu values varied little (−0.12‰ to 0.36‰) in soils and most wastes and therefore no source identification was possible. In soils, Cu became isotopically lighter with increasing depth down to 0.4 m, likely because of equilibrium reactions between dissolved and adsorbed Cu species during transport of smelter-derived Cu through the soil. The δ⁶⁶Zn_IRMM values were isotopically lighter in ash (−0.41‰) and organic horizons (−0.85‰ to −0.47‰) than in bedrock (−0.28‰) and slag (0.18‰) likely mainly because of kinetic fractionation during evaporation and thus allowed for separation of smelter-Zn from native Zn in soil. In particular in the organic horizons large variations in δ⁶⁶Zn values occur, probably caused by biogeochemical fractionation in the soil-plant system. In the mineral horizons, Zn isotopes showed only minor shifts to heavier δ⁶⁶Zn values with depth mainly because of the mixing of smelter-derived Zn and native Zn in the soils. In contrast to Cu, Zn isotope fractionation between dissolved and adsorbed species was probably only a minor driver in producing the observed variations in δ⁶⁶Zn values. Our results demonstrate that metal stable isotope ratios may serve as tracer of sources, vertical dislocation, and biogeochemical behavior in contaminated soil.     

Impact of draining hilly lands on runoff and on‐site erosion: a case study from humid Ethiopia

The use of drainage ditches on farmland has an impact on erosion processes both on‐site and off‐site, though their environmental impacts are not unequivocal. Here we study the runoff response and related rill erosion after installing drainage ditches and assess the effects of stone bunds in north Ethiopia. Three different land management systems were studied in 10 cropland catchments around Wanzaye during the rainy season of 2013: (1) the exclusive use of drainage ditches (locally called feses), (2) the exclusive use of stone bunds, and (3) a mixture of both systems. Stone bunds are an effective soil and water conservation technique, making the land more resistant against on‐site erosion, and allowing feses to be installed at a larger angle with the contour. The mean rill volumes for the 10 studied cropland catchments during the rainy season of 2013 was 3.73 ± 4.20 m³ ha^?1 corresponding to a soil loss of 5.72 ± 6.30 ton ha^?1. The establishment of feses causes larger rill volumes (R = 0.59, N = 10), although feses are perceived as the best way to avoid soil erosion when no stone bunds are present. The use of feses increases event‐based runoff coefficients (RCs) on cropland from c. 5% to values up to 39%. Also, a combination of low stone bund density and high feses density results in a higher RC, whereas catchments with a high stone bund density and low feses density have a lower RC. Peak runoff discharges decrease when stone bund density increases, whereas feses density is positively related to the peak runoff discharge. A multiple linear relation in which both feses and stone bund densities are used as explanatory variable, performs best in explaining runoff hydrograph peakedness (R² = 83%). Copyright © 2015 John Wiley & Sons, Ltd.