Jana Votrubova

MR properties of water in saturated soils and resulting loss of MRI signal in water content detection at 2 tesla

This paper reports a systematic MRI study at 2 tesla of 23 soils, each separately saturated with a known amount of water. The percentage of that water which could be detected using various MR methods was determined by comparison with a liquid reference sample. A pulse-acquire sequence gave a bulk detection of between 47 and 94% of the known water content of saturated soil. Also for bulk measurements, the inversion-recovery sequence used for determining T1 values detected a range of 0.7–75% of the existing soil water. The CPMG sequence with an echo time (TE) of 1 ms used for determining the bulk T2 values gave lower values, in the range of 0.4–66% overall. A spin-echo MRI sequence with a TE of 2.9 ms gave an even lower bulk detection, ranging from 0.2 to 57%. These low values for the detectable water content of bulk saturated soil water reflect the loss of water magnetization which occurs even during short echo time MR sequence at 2 tesla field strength. The source of the above findings was investigated by measurements of the longitudinal (T1) and transverse (T2) relaxation times and spectral linewidths of the soil-water protons, and by conventional analysis of soil properties. The MR parameters of critical importance to water quantification are T2 and T2∗, shorter values of which lead to a progressively greater loss of signal intensity for all MR protocols. Those parameters are affected by the following soil chemical and physical features: soil magnetic susceptibility, and the content of free iron oxides, clay, sand, exchangeable cations (K, Na and Ca), and organic matter. The implication of this work is that the only soil water which can be detected quantitatively at 2 tesla using a conventional spin-echo MRI protocol with echo times of 3 ms or longer is that located in the relatively large soil pores. Using the protocols investigated in this work, water in smaller pores will only be detected accurately for soils which have relatively low paramagnetic-metal impurities and/or have low clay content. Future MR studies of soil water should consider the use of other MRI protocols (e.g. for solid state), and measurement at low magnetic fields.

CT derived porosity distribution and flow domains

Computed X-ray tomography images of dry undisturbed soil samples of the coarse sandy loam of Korkusova Hut (KH) and the fine sand of Hupselse Beek (HB) are analysed concerning the image resolution, the direction of scanning and the evaluation of boundaries between the fast and slow flow domains. The image resolution appears to be a less important factor than the direction of scanning. Additional data sets are recommended for partitioning the fast and slow flow domains. Comparison is made with the tracer breakthrough pictures, magnetic resonance (MR) images and outflow volumes obtained from infiltration outflow experiments. The two soils under study appeared to have essentially different flow dynamics. In the KH sandy loam soil, the fast flow occurs in regions of high local porosity. In case of the HB fine sand, the sample appeared to consist of a dense part, where the flow process took place, and a more porous part, which stayed dry. Flow instability was detected on MR images for the KH soil.

The relationships between MR parameters and the content of water in packed samples of two soils

Abstract The relationship between the known water content and that determined using Magnetic Resonance (MR) methods was studied for samples of two soils packed to various dry bulk densities; four basic MR protocols were employed (single 90° pulse, inversion-recovery, Carr–Purcell–Meiboom–Gill (CPMG), and spin-echo (SE) pulse sequences). The soil types were chosen from a previously studied range of soils to represent both an “easy” (fine sand) and a “difficult” (coarse sandy loam) system for MR application. Whereas for the “easy” soil, the MR-detected water relates linearly to the soil water content as expected, for “difficult” soil the observed relationships are non-linear and depend on MR properties of soil water and MR protocols used.

Rainfall interception and spatial variability of throughfall in spruce stand

Abstract The interception was recognized as an important part of the catchment water balance in temperate climate. The mountainous forest ecosystem at experimental headwater catchment Liz has been subject of long-term monitoring. Unique dataset in terms of time resolution serves to determine canopy storage capacity and free throughfall. Spatial variability of throughfall was studied using one weighing and five tipping bucket rain gauges. The basic characteristics of forest affecting interception process were determined for the Norway spruce stand at the experimental area - the leaf area index was 5.66 - 6.00 m2 m-2, the basal area was 55.7 m2 ha-1, and the crown closure above individual rain gauges was between 19 and 95%. The total interception loss in both growing seasons analyzed was 34.5%. The mean value of the interception capacity determined was about 2 mm. Throughfall exhibited high variability from place to place and it was strongly affected by character of rainfall. On the other hand, spatial pattern of throughfall in average showed low variability.

Interpretation of ponded infiltration data using numerical experiments

Abstract Ponded infiltration experiment is a simple test used for in-situ determination of soil hydraulic properties, particularly saturated hydraulic conductivity and sorptivity. It is known that infiltration process in natural soils is strongly affected by presence of macropores, soil layering, initial and experimental conditions etc. As a result, infiltration record encompasses a complex of mutually compensating effects that are difficult to separate from each other. Determination of sorptivity and saturated hydraulic conductivity from such infiltration data is complicated. In the present study we use numerical simulation to examine the impact of selected experimental conditions and soil profile properties on the ponded infiltration experiment results, specifically in terms of the hydraulic conductivity and sorptivity evaluation. The effect of following factors was considered: depth of ponding, ring insertion depth, initial soil water content, presence of preferential pathways, hydraulic conductivity anisotropy, soil layering, surface layer retention capacity and hydraulic conductivity, and presence of soil pipes or stones under the infiltration ring. Results were compared with a large database of infiltration curves measured at the experimental site Liz (Bohemian Forest, Czech Republic). Reasonably good agreement between simulated and observed infiltration curves was achieved by combining several of factors tested. Moreover, the ring insertion effect was recognized as one of the major causes of uncertainty in the determination of soil hydraulic parameters.

Ponded infiltration in a grid of permanent single-ring infiltrometers: Spatial versus temporal variability

Abstract Temporal variability of the soil hydraulic properties is still an open issue. The present study deals with results of ponded infiltration experiments performed annually in a grid of permanent measurement points (18 spatial and 14 temporal replicates). Single ring infiltrometers were installed in 2003 at a meadow site in the Bohemian Forest highlands, the Czech Republic. The soil at the plot is coarse sandy loam classified as oligotrophic Eutric Cambisol. Soil water flow below infiltration rings has distinctly preferential character. The results are marked with substantial interannual changes of observed infiltration rates. Considering just the results from the initial four years of the study, the temporal variability did not exceed the spatial variability detected in individual years. In later years, a shift to extremely high infiltration rates was observed. We hypothesize that it is related to structural changes of the soil profile possibly related to combined effect of soil biota activity, climatic conditions and experimental procedure. Interestingly, the temporal changes can partly be described as fluctuations between seemingly stable infiltration modes. This phenomenon was detected in the majority of rings and was found independent of the initial soil moisture conditions.

Episodic runoff generation at Central European headwater catchments studied using water isotope concentration signals

Abstract Hydrological monitoring in small headwater catchments provides the basis for examining complex interrelating hydraulic processes that govern the runoff generation. Contributions of different subsurface runoff mechanisms to the catchment discharge formation at two small forested headwater catchments are studied with the help of their natural isotopic signatures. The Uhlirska catchment (Jizera Mts., Czech Republic) is situated in headwater area of the Lusatian Neisse River. The catchment includes wetlands at the valley bottom developed over deluviofluvial granitic sediments surrounded by gentle hillslopes with shallow soils underlain by weathered granite. The Liz catchment (Bohemian Forest, Czech Republic) is situated in headwater area of the Otava River. It belongs to hillslope-type catchments with narrow riparian zones. The soil at Liz is developed on biotite paragneiss bedrock. The basic comparison of hydrological time series reveals that the event-related stream discharge variations at the Uhlirska catchment are bigger and significantly more frequent than at Liz. The analysis of isotope concentration data revealed different behavior of the two catchments during the major rainfall-runoff events. At Uhlirska, the percentage of the direct runoff formed by the event water reaches its maximum on the falling limb of the hydrograph. At Liz, the event water related fraction of the direct outflow is maximal on the rising limb of the hydrograph and then lowers. The hydraulic functioning of the Uhlirska catchment is determined by communication between hillslope and riparian zone compartments.