Continuum percolation for Gibbs point processes

We consider percolation properties of the Boolean model generated by a Gibbs point process and balls with deterministic radius. We show that for a large class of Gibbs point processes there exists a critical activity, such that percolation occurs a.s. above criticality. For locally stable Gibbs point processes we show a converse result, i.e. they do not percolate a.s. at low activity.

Transport of iodide in structured clay–loam soil under maize during irrigation experiments analyzed using HYDRUS model

Transport of radioactive iodide 131I− in a structured clay loam soil under maize in a final growing phase was monitored during five consecutive irrigation experiments under ponding. Each time, 27 mm of water were applied. The water of the second experiment was spiked with 200 MBq of 131I− tracer. Its activity was monitored as functions of depth and time with Geiger-Müller (G-M) detectors in 11 vertically installed access tubes. The aim of the study was to widen our current knowledge of water and solute transport in unsaturated soil under different agriculturally cultivated settings. It was supposed that the change in 131I− activity (or counting rate) is proportional to the change in soil water content. Rapid increase followed by a gradual decrease in 131I− activity occurred at all depths and was attributed to preferential flow. The iodide transport through structured soil profile was simulated by the HYDRUS 1D model. The model predicted relatively deep percolation of iodide within a short time, in a good agreement with the observed vertical iodide distribution in soil. We found that the top 30 cm of the soil profile is the most vulnerable layer in terms of water and solute movement, which is the same depth where the root structure of maize can extend.

Tracking water pathways in steep hillslopes by d18O depth profiles of soil water

Assessing temporal variations in soil water flow is important, especially at the hillslope scale, to identify mechanisms of runoff and flood generation and pathways for nutrients and pollutants in soils. While surface processes are well considered and parameterized, the assessment of subsurface processes at the hillslope scale is still challenging since measurement of hydrological pathways is connected to high efforts in time, money and personnel work. The latter might not even be possible in alpine environments with harsh winter processes. Soil water stable isotope profiles may offer a time-integrating fingerprint of subsurface water pathways. In this study, we investigated the suitability of soil water stable isotope (d18O) depth profiles to identify water flow paths along two transects of steep subalpine hillslopes in the Swiss Alps. We applied a one-dimensional advection–dispersion model using d18O values of precipitation (ranging from _24.7 to _2.9‰) as input data to simulate the d18O profiles of soil water. The variability of d18O values with depth within each soil profile and a comparison of the simulated and measured d18O profiles were used to infer information about subsurface hydrological pathways. The temporal pattern of d18O in precipitation was found in several profiles, ranging from _14.5 to _4.0‰. This suggests that vertical percolation plays an important role even at slope angles of up to 46_. Lateral subsurface flow and/or mixing of soil water at lower slope angles might occur in deeper soil layers and at sites near a small stream. The difference between several observed and simulated d18O profiles revealed spatially highly variable infiltration patterns during the snowmelt periods: The d18O value of snow (_17.7 ± 1.9‰) was absent in several measured d18O profiles but present in the respective simulated d18O profiles. This indicated overland flow and/or preferential flow through the soil profile during the melt period. The applied methods proved to be a fast and promising tool to obtain time-integrated information on soil water flow paths at the hillslope scale in steep subalpine slopes.

Redox front geochemistry and weathering: theory with application to the Osamu Utsumi uranium mine, Poços de Caldas, Brazil

Numerical calculations describing weathering of the Poços de Caldas alkaline complex (Minas Gerais, Brazil) by infiltrating groundwater are carried out for time spans up to two million years in the absence of pyrite, and up to 500,000 years with pyrite present. Deposition of uranium resulting from infiltration of oxygenated, uranium bearing groundwater through the hydrothermally altered phonolitic host rock at the Osamu Utsumi uranium mine is also included in the latter calculation. The calculations are based on the quasi-stationary state approximation to mass conservation equations for pure advective transport. This approximation enables the prediction of solute concentrations, mineral abundances and porosity as functions of time and distance over geologic time spans. Mineral reactions are described by kinetic rate laws for both precipitation and dissolution. Homogeneous equilibrium is assumed to be maintained within the aqueous phase. No other constraints are imposed on the calculations other than the initial composition of the unaltered host rock and the composition of the inlet fluid, taken as rainwater modified by percolation through a soil zone. The results are in qualitative agreement with field observations at the Osamu Utsumi uranium mine. They predict a lateritic cover followed by a highly porous saprolitic zone, a zone of oxidized rock with pyrite replaced by iron-hydroxide, a sharp redox front at which uranium is deposited, and the reduced unweathered host rock. Uranium is deposited in a narrow zone located on the reduced side of the redox front in association with pyrite, in agreement with field observations. The calculations predict the formation of a broad dissolution front of primary kaolinite that penetrates deep into the host rock accompanied by the precipitation of secondary illite. Secondary kaolinite occurs in a saprolitic zone near the surface and in the vicinity of the redox front. Gibbsite forms a bi-modal distribution consisting of a maximum near the surface followed by a thin tongue extending downward into the weathered profile in agreement with field observations. The results are found to be insensitive to the kinetic rate constants used to describe mineral reactions.