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

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.     

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.     

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.     

No Higher Criticism of the Bianchi-corrected Wilkinson Microwave Anisotropy Probe data

Higher Criticism (HC) has been proposed by Donoho & Jin as an effective statistic to detect deviations from Gaussianity. Motivated by the success of the Bianchi VII_h model in addressing many of the anomalies observed in the Wilkinson Microwave Anisotropy Probe ( WMAP ) data (Jaffe et al.), we present calculations in real and in wavelet space of the HC statistic of the Bianchi-corrected WMAP first-year data. At the wavelet scale of 5°, the HC of the WMAP map drops from a value above the 99 per cent confidence level (c.l.) to a value below the 68 per cent CL when corrected by the Bianchi template. An important property of the HC statistic is its ability to locate the pixels that account for the deviation from Gaussianity. The analysis of the uncorrected WMAP data pointed to a cold spot in the Southern hemisphere, centred at   l ∼ 209°, b ∼−57°  . The HC of the Bianchi-corrected map indicates that this spot remains prominent, albeit at a level completely consistent with Gaussian statistics. Consequently, it is debatable how much emphasis should be placed on this residual feature, but we consider the effect of modestly increasing the scaling of the template. A factor of only 1.2 renders the spot indistinguishable from the background level, with no noticeable impact on the results published in Jaffe et al. for the low- l anomalies, large-scale power asymmetry or wavelet kurtosis. A trivial interpretation would be that the Bianchi template may require a small enhancement of power on scales corresponding to the wavelet scale of 5°.     

On the character of distribution of chemical elements in nature

Distribution of chemical elements in nature is controlled by the general physical-chemical state of natural processes, a result of the effort to reach thermodynamic equilibrium. Natural processes, forming geologic objects (e.g., rocks, ores) usually do not reach total chemical equilibrium. This is reflected in distributions of chemical elements in nature. We recognize symmetric as well as asymmetric distributions. The question arises, what distribution in nature is most commonly followed? Evidently, this should be a distribution of chemical elements in thermodynamic equilibrium throughout the system, thus with maximum entropy. This paper notes that Gauss's normal distribution may be derived from statistic explanation of entropy as a function of contnt of chemical elements in a system with application of Einstein's formula.