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).     

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.     

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.     

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.     

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.     

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.     

The Indus-Tsangpo suture zone in Ladakh, northwest Himalaya: Further palaeomagnetic data and implications

A total of 81 samples (244 specimens) from Upper Cretaceous Indus Molasse and Middle to Upper Cretaceous Dras Flyschoids of the Indus-Tsangpo suture zone in Ladakh (northwest Himalaya) has been studied by thermal demagnetization methods.Both formations showed a characteristic magnetization component indicative for equatorial to low northern palaeolatitudes of acquisition. Similar palaeolatitudes have been obtained before from secondary magnetization components of Early Tertiary age in the Ladakh Intrusives and in the Tibetan Sedimentary Series of central Nepal. The present characteristic components are interpreted likewise as secondary magnetizations which stabilized between 50 and 60 m.y. ago, during Greater India's collision with Asia's southern margin.The Dras Flyschoids show another magnetic component which, in case of primary origin, indicates acquisition at a low southern palaeolatitude. If correct, this interpretation supports recent suggestions for Late Cretaceous obduction of an island arc on Greater India's northern margin.     

Impact of soil and water conservation measures on catchment hydrological response—a case in north Ethiopia

Impact studies of catchment management in the developing world rarely include detailed hydrological components. Here, changes in the hydrological response of a 200‐ha catchment in north Ethiopia are investigated. The management included various soil and water conservation measures such as the construction of dry masonry stone bunds and check dams, the abandonment of post‐harvest grazing, and the establishment of woody vegetation. Measurements at the catchment outlet indicated a runoff depth of 5 mm or a runoff coefficient (RC) of 1·6% in the rainy season of 2006. Combined with runoff measurements at plot scale, this allowed calculating the runoff curve number (CN) for various land uses and land management techniques. The pre‐implementation runoff depth was then predicted using the CN values and a ponding adjustment factor, representing the abstraction of runoff induced by the 242 check dams in gullies. Using the 2006 rainfall depths, the runoff depth for the 2000 land management situation was predicted to be 26·5 mm (RC = 8%), in line with current RCs of nearby catchments. Monitoring of the ground water level indicated a rise after catchment management. The yearly rise in water table after the onset of the rains (ΔT) relative to the water surplus (WS) over the same period increased between 2002–2003 (ΔT/WS = 3·4) and 2006 (ΔT/WS >11·1). Emerging wells and irrigation are other indicators for improved water supply in the managed catchment. Cropped fields in the gullies indicate that farmers are less frightened for the destructive effects of flash floods. Due to increased soil water content, the crop growing period is prolonged. It can be concluded that this catchment management has resulted in a higher infiltration rate and a reduction of direct runoff volume by 81% which has had a positive influence on the catchment water balance. Copyright © 2010 John Wiley & Sons, Ltd.     

On the origin of rock fragment mulches on Vertisols: A case study from the Ethiopian highlands

Many Vertisols in Tigray, Ethiopia, typically carry a discontinuous rock fragment (RF, size 0.5–> 40 · 10^− 2 m) cover with 10 to 100 RFs m^− 2. Such RF mulches are of agricultural and environmental significance because they influence the water balance in the underlying soils and the crop yield. Natural RF concentrations are mostly considered as eolian or hydraulic lag deposits, or as the result of lateral transport over the soil surface from a rock outcrop, upslope. In cultivated areas RF mulches can develop by tillage.This paper presents the case of a natural RF mulch whose lithology indicates that the RFs are up-squeezed by the local Vertisol. The study site is located in the pass of Enda Maryam, Tigray, Northern Ethiopia (39°8′ E and 13°36′ N). A circular area of 10 m diameter, about 200 m away from the water divide in the valley has been cleared annually between 01/1999 and 05/2003. During this period, 625 RFs, 17 being > 7.5 · 10^− 2 m in size, totalling a mass of nearly 62 kg, have been collected. After correction for measurement procedures, the rate of RF up-warping by the Vertisol at Enda Maryam is assessed at 5 RFs m^− 2 in 3 years. At this rate of appearance, the formation of current RF concentrations on top of active valley Vertisols is only a matter of 10^1–2 years, provided the availability of RFs below the soil surface.Although important underground displacements were measured in the Vertisol between 01/1999 and 05/2002, the supposed link between up-squeezing of RFs and plastic deformations of ‘chimney’, ‘diapir’ or ‘intrusion’-like type in the Vertisol could not be evidenced. Instead, RFs are clearly concentrated on the soil surface as well as in depth, along the existing vertical desiccation cracks, often > 1 m deep which display polygonal configurations at the soil surface. Further, bundles of slickensides containing some RFs, have been mapped at the base of the Vertisol. The slickenside configuration suggests that the RF-bearing substrate is being scraped off.While the underground displacement of RFs along active slickensides seems normal, the process of RFs ascending in ‘upright’ position in the edge of desiccation cracks needs explanation. The closure of a desiccation crack is a peristaltic-like movement, following ascent or descent of the capillary fringe. It is hypothesized that this movement gradually pushes the RF to the surface or to another place or level in the soil profile where the crack closes in last instance.The apparent young age of the valley Vertisol mulches in Ethiopia might indicate the very recent formation of yearly recurrent desiccation cracks of Vertisols in the area. Available information confirms that most valleys in the study area used to be perennially marshy. Under these conditions no movements of RFs in the soil profile are expected to occur. Gullying, leading to pronounced seasonal desiccation of the Vertisols, started in several cases not more than 50 years ago.