Outline
Comptes Rendus

Surface geosciences (Pedology)
3D representation of soil distribution: An approach for understanding pedogenesis
Comptes Rendus. Géoscience, Volume 341 (2009) no. 6, pp. 486-494.

Abstracts

Soils are characterised by a spatial variability in the three dimensions (3D) of space. However, 3D studies remain scarce due to the qualitative nature of many soil horizon characteristics, notably the horizon designation. Indeed, existing 3D tools are mainly developed for quantitative data. To solve this difficulty, we propose a new approach based on the interpolation of the horizon thickness to derive digital elevation models for both the upper and the lower limits of each horizon. This approach was applied to Planosols previously extensively studied with 2D approaches. The pseudo 3D obtained representation evidences soil processes that were missed in 2D approaches. As an example, we evidence the impact of differential weathering, resulting from the mineralogical heterogeneity of the parent material, on the subsequent pedogenesis.

Les sols sont des continuums tridimensionnels (3D), caractérisés par une variabilité spatiale dans l’ensemble des directions de l’espace, mais les études 3D demeurent rares. Une des difficultés majeures des approches 3D est due à la nature qualitative de nombreuses caractéristiques des horizons, dans la mesure où la plupart des outils 3D existants sont principalement développés pour l’interpolation 3D de données quantitatives. En substitution à l’interpolation directe de la variable « nom de l’horizon », nous proposons une approche basée sur l’interpolation de l’épaisseur des horizons pour calculer les modèles numériques d’altitude correspondant aux limites de chaque horizon, permettant ainsi une visualisation pseudo 3D. Cette étude, appliquée à des Planosols bien caractérisés en 2D par ailleurs, montre que la représentation 3D des sols permet de souligner des processus omis lors des approches précédentes, tels que l’impact de l’altération différentielle, résultant de l’hétérogénéité initiale du matériau parental, sur une pédogenèse ultérieure.

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/j.crte.2009.04.004
Keywords: 3D representation, Modelling, Horizon, Soil volume, Pedogenesis, Parent material, Soil spatial distribution
Mot clés : Représentation 3D, Modélisation, Horizon, Volume de sol, Pédogenèse, Matériau parental, Distribution spatiale des sols

Frédéric Delarue 1; Sophie Cornu 2; Joël Daroussin 2; Sébastien Salvador-Blanes 3; Hocine Bourennane 2; Patrick Albéric 1; Aurélie Vennink 1; Ary Bruand 1; Dominique King 2

1 Université d’Orléans, université François-Rabelais – Tours, CNRS/INSU, institut des sciences de la Terre d’Orléans–UMR 6113, campus géosciences, 1A, rue de la Férollerie, 45071 Orléans cedex 2, France
2 INRA, unité de science du sol, centre de recherches d’Orléans, avenue de la Pomme-de-Pin, BP 20619, 45166 Olivet cedex, France
3 Université François-Rabelais–Tours, université d’Orléans, CNRS/INSU, institut des sciences de la Terre d’Orléans–UMR 6113, faculté des sciences et des techniques, parc Grandmont, 37200 Tours, France
@article{CRGEOS_2009__341_6_486_0,
     author = {Fr\'ed\'eric Delarue and Sophie Cornu and Jo\"el Daroussin and S\'ebastien Salvador-Blanes and Hocine Bourennane and Patrick Alb\'eric and Aur\'elie Vennink and Ary Bruand and Dominique King},
     title = {3D representation of soil distribution: {An} approach for understanding pedogenesis},
     journal = {Comptes Rendus. G\'eoscience},
     pages = {486--494},
     publisher = {Elsevier},
     volume = {341},
     number = {6},
     year = {2009},
     doi = {10.1016/j.crte.2009.04.004},
     language = {en},
}
TY  - JOUR
AU  - Frédéric Delarue
AU  - Sophie Cornu
AU  - Joël Daroussin
AU  - Sébastien Salvador-Blanes
AU  - Hocine Bourennane
AU  - Patrick Albéric
AU  - Aurélie Vennink
AU  - Ary Bruand
AU  - Dominique King
TI  - 3D representation of soil distribution: An approach for understanding pedogenesis
JO  - Comptes Rendus. Géoscience
PY  - 2009
SP  - 486
EP  - 494
VL  - 341
IS  - 6
PB  - Elsevier
DO  - 10.1016/j.crte.2009.04.004
LA  - en
ID  - CRGEOS_2009__341_6_486_0
ER  - 
%0 Journal Article
%A Frédéric Delarue
%A Sophie Cornu
%A Joël Daroussin
%A Sébastien Salvador-Blanes
%A Hocine Bourennane
%A Patrick Albéric
%A Aurélie Vennink
%A Ary Bruand
%A Dominique King
%T 3D representation of soil distribution: An approach for understanding pedogenesis
%J Comptes Rendus. Géoscience
%D 2009
%P 486-494
%V 341
%N 6
%I Elsevier
%R 10.1016/j.crte.2009.04.004
%G en
%F CRGEOS_2009__341_6_486_0
Frédéric Delarue; Sophie Cornu; Joël Daroussin; Sébastien Salvador-Blanes; Hocine Bourennane; Patrick Albéric; Aurélie Vennink; Ary Bruand; Dominique King. 3D representation of soil distribution: An approach for understanding pedogenesis. Comptes Rendus. Géoscience, Volume 341 (2009) no. 6, pp. 486-494. doi : 10.1016/j.crte.2009.04.004. https://comptes-rendus.academie-sciences.fr/geoscience/articles/10.1016/j.crte.2009.04.004/

Version originale du texte intégral

1 Introduction

Soil is a three-dimensional (3D) continuum characterised by a spatial variability in all the dimensions of space. Its thickness (generally decimetres to metres) is relatively small compared to its lateral extension (from hectometres to kilometres). This represents a substantial difficulty for 3D approaches. Thus, despite their obvious 3D organisation, soils are generally studied in two dimensions by pedological approaches in sequences or by geostatistics.

Pedological approaches in sequences were first theorised by Jenny [14] who studied the impact of a single factor on the evolution of the soil through time and space. For instance, he discussed the impact of the slope by studying soil distribution along toposequences, keeping all other factors as constant as possible. In this case, topography is only represented by one terrain attribute. This approach has been widely used in many pedological studies and remains a valuable tool to study soil evolution as recently discussed by Huggett [13]. However, soil sequences are generally designed to emphasise the impact of one factor of variation – often a single attribute, the slope –, while the others are generally neglected. In order to overtake this limitation, an approach called “structural analysis” was developed [3,4,9,16]. In these approaches, 3D blocks diagrams of soil distribution were derived from the analysis of several toposequences varying in intensity and orientation of the slope. As these block diagrams are based on the interpretation of the toposequences, they do not allow the observation of variations that would have been missed in the 2D approach.

Approaches based on geostatistics were developed in order to acquire an exhaustive quantification of the soil variability [17]. They aim classically at obtaining the spatial estimate of a soil variable. Nevertheless, these approaches deal mainly only with one horizon and more rarely with two that are compared [5].

Salvador-Blanes and coworkers [23,24] demonstrated that maps are useful for understanding pedogenesis by coupling toposequence and geostatistical approaches. We deduce from this work that a 3D representation of soils, as 2D maps, should improve our understanding of pedogenesis.

Despite many attempts in recent years [1,2,12,20,21], 3D numerical modelling of soils can still be considered as being in its infancy and is generally based on soil properties [22,26] less frequently on the organisation of soil volumes in the 3D space.

In this paper, we demonstrate how a 3D representation of soils evidences soil processes that were missed in the two dimensional approaches. For the sake of our demonstration, we will discuss the distribution of Planosols since their pedogenesis is driven by lateral processes [8]. These soils represent thus an ideal case for such a study.

2 Study site and sampling method

2.1 The study site

The study site is located on the Aigurande plateau in the northern part of the French Massif Central (Fig. 1). It is underlain by amphibolite and gneiss. It is a footslope where Planosols developed on B-horizons and alterites of gneiss [23,24]. The upper horizons of the Planosols are allochthonous and derived from colluvial materials of gneissic origin, while the B-horizon (below 50 cm depth) and the alterite are developed from gneiss. Colluvial materials differentiate into four horizons: a clay-loam organic-rich A-horizon (25 to 30 cm thick) and three types of E-horizons [6], noted E1, E2 and E3. E1 contains no nodules unlike E2 and E3, that both are nodule-rich. E3 also contains gravels.

Fig. 1

Location of studied sector after Salvador-Blanes et al. [25].

Fig. 1. Localisation du secteur étudié d’après Salvador-Blanes et al. [25].

2.2 Sampling and available data

In the studied zone, 105 sampling points were selected randomly according to two square cell grids (with 25 m and 40 m cell size respectively) during two sampling campaigns. These points were sampled with a 1.2 m auger. In the case of the sampling points located on the 25 m cell grid, soils were dug to the bedrock or to the end of the auger where the bedrock was deeper than 1.2 m. In the case of the 40 m square grid, soils were only augered to the top of B-horizon to allow surface water sampling. Each sampling point is characterised by its coordinates (X, Y and Z), the type and the thickness of each horizon, and its descriptive characteristics: texture, colour, visual abundance of gravels and nodules. Abundances in both gravels and nodules were classified, based on the soil descriptions, into the following classes: none, weakly abundant, moderately abundant, abundant and strongly abundant. B-horizons and alterites could not be differentiated by augering. They will thus be named in the remainder of this paper as “B-horizons” and characterised by their texture or redox status.

A detailed Digital Elevation Model (DEM) [25] of the area was also obtained with 4600 elevation points measured in the field using a laser tacheometer (Leica TC500).

Manganese was analysed in soil samples from horizon E3 by ICP-AES after fluoridric acid (HF) dissolution and in the free water sampled from the upper horizons (A to E3) by ICP-AES.

3 3D modelling

Horizon designation is a qualitative data. To our knowledge, available 3D interpolation methods are mostly developed for quantitative data [19]. For that reason, only few references of 3D modelling of soil horizons exist in the literature [10,11,27].

Grundwald and Barak [10,11] presented a 3D soil-landscape using virtual reality modeling language (VRML) and interpolating data in the horizontal and vertical dimensions separately.

Santos et al. [27] developed a 3D soil horizon model in a Swiss floodplain using a quadratic finite-element method. Our approach is partially derived from their work [27]. Fig. 2 summarises the various steps performed to obtain “pseudo” 3D visualisations of Planosols. Both ArcInfo® and ArcGIS® softwares were used. The pedological data are converted to a triangular irregular network (TIN) structure and then spatialised using 2D interpolation to cells raster datasets. The resulting raster represents the elevation of the soil surface, the thickness of each horizon and the other characteristics of each horizon. Rasters for the elevation of the soil surface and the horizon thicknesses can be combined to obtain the horizons upper limit elevations. These are then used to obtain “pseudo” 3D representation by horizon thickness extrusion.

Fig. 2

Schematic representation of the 3D approach.

Fig. 2. Représentation schématique de l’approche tridimensionnelle.

4 Impact of the spatial variability of the poorly permeable horizons on the development of the upper horizons

We first studied the 3D variability of the B-horizons, which act as a poorly permeable layer over which the lateral runoff occurs. These horizons are thus responsible for the development of the Planosol [6]. The aim was to identify whether the variability in thickness or characteristics of these horizons could influence the subsequent development of the upper Planosolic horizons and their variability described by Cornu et al. [7]. For that reason, two main characteristics of the B-horizons were considered: their oxidation status considering three classes (oxidized, gley and gley that was reoxydized) and their texture.

The 3D representation evidences the existence of a spatial organisation of the various B-horizons that was not demonstrated in the previous studies [23]. Two main features were evidenced: the presence of a heavy clay hillock (Fig. 3a) that is north-south oriented and located in the middle of the studied area, and the presence of gleyic horizons at depth.

Fig. 3

a: 3D representation of the various B-horizons and subsurface topography; b: spatial variability of E3-horizon thickness; c: abundance in gravels in the E3-horizon and subsurface drainage; d: abundance in nodules in the E3-horizon.

Fig. 3. a : représentation tridimensionnelle de la variabilité spatiale des horizons B et topographie de subsurface ; b : variabilité spatiale de l’épaisseur de l’horizon E3 ; c : abondance en éléments grossiers dans l’horizon E3 et écoulement de subsurface ; d : abondance en concrétions dans l’horizon E3.

Two textural gradients depart from the clayey hillock as the texture becomes coarser from the upper B-horizon to the deepest one and from the hillock to the edges of the studied area. This is interpreted as a weathering gradient – the finer-textured horizons being the most weathered ones – resulting from the original heterogeneity of the gneiss observed by Salvador-Blanes [23].

This subsurface topography results in two distinct subsurface drainage networks roughly north-south oriented (Fig. 4) while the surface runoff is mainly north-west to south-east oriented (Fig. 4d).

Fig. 4

a: manganese concentrations in E3-horizon; b: manganese concentrations in the perched temporary water table; c: thickness of the upper B-horizon; d: surface runoff.

Fig. 4. a : concentrations en manganèse dans l’horizon E3 ; b : concentrations en manganèse dans la nappe perchée temporaire ; c : épaisseur de l’horizon B superficiel ; d : écoulement de surface.

As most of the transfers occurring in this soil are related to the presence of the clayey B-horizon, this spatial distribution of the subsurface drainage network should have consequences on the morphological differentiation of the upper soil horizons of the Planosol. This variability was thus studied in link with that of the underlying B-horizons. Cornu et al. [7] extensively studied the variability of the E-horizons along a toposequence. They evidence three main parameters of variability of these horizons: their thickness, their amount of nodules and of gravels. Gravels are more concentrated at the west and east sides of the heavy clay hillock where the subsurface topography is concave (Fig. 3c) and the E3-horizon is thicker (Fig. 3b); no clear correlation was however derived. The abundance of Fe-Mn nodules, in horizon E3, do not show a clear spatial distribution pattern with the exception of a lower abundance on the east side and of a higher abundance (Fig. 3d) on the west side of the hillock. The thickening of the E3-horizon and the larger abundance in gravels can be explained by filling and venting processes that occurred during the formation of the upper allochtonous horizons in link with the morphology of the B-horizon. This does not however explain the dissymmetry in abundance of nodules on both sides of the hillock.

5 A possible connection of the groundwaters

Horizons E2 and E3 are characterised by the alternation of reducing and oxidizing conditions. These are evidenced by the presence of Fe-Mn nodules [6]. This is due to the presence of a temporary perched water table in these horizons observed in the field in winter conditions (from January to April mainly). The formation of soil nodules by alternation of reducing and oxidizing conditions is well described elsewhere [15,18,29].

These horizons exhibit a thickening along the two sides of the previously described hillock, with a higher content in nodules on the west side resulting in higher concentrations in Mn (Fig. 4a). This enrichment of E3-horizons in Mn also corresponds to a higher richness in the perched temporary water table (Fig. 4b). The question is to determine whether the composition of the E3-horizon is responsible for the composition of the water, that is to say the nodules are undergoing dissolution, or the water composition is responsible for the present formation of the nodules into the E3-horizon with an external source of Mn.

Two potential external sources of dissolved Mn can be hypothesised, either a surface origin or a deep one. Indeed, the north-south drainage axis may drain waters from the neighbouring plot in which soils are developed from amphibolite. These soils are rich in Mn [23,24] and may thus produce soil water rich in Mn. Alternatively, gleyic B-horizons were observed under the clayey hillock, indicating the presence of a permanent and reducing groundwater table. Fig. 4c shows that the thickness of the upper B-horizon decreases in the western part of the studied sector, being only 0 to 10 cm thick, while it is thicker on the remainder of the sector. Thus, a possible connexion between the temporary surface water circulation and the deep groundwater located in the gleyic B-horizon could occur in this area. Indeed, the deep groundwater is reductive, it may be rich in dissolved Mn. Other samplings are needed to confirm that hypothesis as only few of them were dug in this area.

These two hypotheses may explain the spatial distribution of Mn in the water of the perched water table. A detailed characterisation of the spatial variability of these waters is needed to definitely conclude on the validity of these hypotheses.

6 Conclusion: inputs of the 3D approach to the understanding of the Planosols formation

The studied Planosols are developed on B-horizons and alterites of gneiss as already demonstrated elsewhere [7,23]. The parent gneiss material was shown to be mineralogically highly heterogeneous (Fig. 5a; [23]).

Fig. 5

Schematic representation of the Planosols formation at Bois Maillot. a: parent gneiss mineralogically highly heterogeneous; b: differential weathering leading to a contrasted texture B-horizons; c: differential erosion phase; d: deposits of quartz gravels and fine colluvium; e: appearance of subsurface lateral flow at the upper limit of the B-horizons and subsequent formation of nodules; f: current distribution of the Planosols.

Fig. 5. Représentation schématique de la formation de la couverture Planosolique du Bois Maillot. a : gneiss parental très hétérogène minéralogiquement ; b : phase d’altération différentielle à l’origine des horizons B à textures contrastées ; c : phase d’érosion différentielle ; d : dépôt de graviers de quartz et des colluvions fines ; e : apparition de flux latéraux à la surface des horizons B et formation des concrétions ; f : distribution actuelle des Planosols.

On the basis of the B-horizons spatial variability, we suggest the following soil evolution sketch schematically represented in Fig. 5b. The weathering of the heterogeneous parent gneiss led to the development of a vertical and a lateral textural gradient along the studied area. The formation of the subsurface hillock is supposed to be due to a differential erosion phase (Fig. 5c). It could have occurred under a colder climate than that leading to the weathering of the parent gneiss material. Such succession of processes has been previously described in the literature in other pedogenetic contexts [28]. In addition, clayey material are known to be more difficult to erode than sandy ones, thus leading to the development of the present subsurface topography, including two channels on both sides of the hillock.

During the periglacial period, a new erosion phase allowed the deposits of coarse quartz (Fig. 5d). The channels along the hillock may have concentrated the deposit of gravels due to the differential surface topography at that time. Fine colluvium filled the little dimple of subsurface topography, leading to a thicker E3-horizon especially in the western part of the study area.

The textural difference between the gravels and the B-horizon induced preferential lateral flow at the top of the B-horizon and the presence of a temporary perched water table, leading to the formation of the nodules within the E2 and E3-horizons (Fig. 5e). The subsurface topography also induced two preferential water drains at the surface of the B-horizons roughly north-south oriented. However, a discrepancy in the abundance of nodules along both sides of the hillock is observed. Various hypotheses are invoked to explain this phenomenon: a higher flux concentration along the west of the hillock or the existence of another source of water, rich in Mn in that part of the area.

This pseudo 3D approach brings new elements for understanding the spatial distribution of these Planosols (Fig. 5f). Previous studies showed that they were influenced by their position along the slope [6] and by anthropogenic features [25]. We here demonstrate that subsurface topography, resulting from differential weathering that expresses the original mineralogical heterogeneity of the gneiss parent material, also influences their spatial differentiation. As a conclusion, this approach allows evidencing the consequences of the initial weathering phases on a later pedogenesis.

We thus demonstrated that 3D modelling is a useful tool for understanding pedogenesis in addition to traditional 2D approaches. Nevertheless, it should be improved on various aspects. Indeed, the technique of spatialisation performed in this study is very simple. It could be improved by the use of geostatistics that would allow assessing the spatial structure of the studied variables and estimating errors of prediction. Finally, the use of exhaustive data – obtained by geophysical survey notably – could increase our knowledge of the spatial distribution of some of the key horizons. This data could then be re-injected as input data in the 3D model.

However, this study is a first attempt at soil 3D modelling. In the absence of 3D interpolation tool for soil qualitative data, this pseudo 3D approach represents an efficient compromise, quick and simple, which highlights processes that were missed by 2D approaches.

Acknowledgment

The authors are grateful to the Region Centre that funded this work as part of the Metaloe-Metaspace project.


References

[1] M. Ameskamp, Three-dimensional rule-based continuous soil modelling, PhD Thesis, Christian-Albrechts-Universität, Kiel, Germany, 1997.

[2] R. Benz Essai d’approche tridimensionnelle de la couverture pédologique : application à trois stations d’étude dans la région du col du Marchairuz Haut-Jura Vaudois, Université de Neuchâtel, Institut de Botanique, Laboratoire d’Ecologie Végétale, Neuchâtel, 1995

[3] R. Boulet; A. Chauvel; F.X. Humbel; Y. Lucas Analyse structurale et cartographie en pédologie II, Cahiers de l’Orstom, Série Pédologie., Volume XXII (1982), pp. 323-329

[4] R. Boulet; A. Chauvel; F.X. Humbel; Y. Lucas Analyse structurale et cartographie en pédologie II, cahiers de l’Orstom, Série Pédologie, Volume XXII (1982), pp. 309-321

[5] H. Bourennane; S. Salvador-Blanes; S. Cornu; D. King Scale of spatial dependence between chemical properties of topsoil and subsoil over a geologically contrasted area (Massif central, France), Geoderma, Volume 112 (2003), pp. 235-251

[6] S. Cornu; V. Deschatrettes; S. Salvador-Blanes; B. Clozel; M. Hardy; S. Branchut; L. Le Forestier Trace element accumulation in Mn-Fe-oxide nodules of a Planosolic horizon, Geoderma, Volume 125 (2005), pp. 11-24

[7] S. Cornu; I. Cousin; V. Deschatrettes; N. Saby; S. Salvador-Blanes; C. Crouzet; A. Guerin; B. Clozel Consequences of aggregation for the trace element distribution in the subsoil of a Planosol naturally rich in trace metal, Geoderma, Volume 136 (2006), pp. 160-173

[8] FAO, World Reference Base for Soil Resources, Report N°. 84, Food and Agriculture Organization of the United Nations, Rome, 1998.

[9] E. Fritsch; G. Bocquier; R. Boulet; F.X. Humbel Les systèmes transformants d’une couverture ferrallitique de Guyane française: analyse structurale d’une formation supergène et mode de représentation, Cahiers ORSTOM, Série Pédologie., Volume XXII (1986), pp. 361-395

[10] S. Grunwald; P. Barak The use of VRML for virtual soil landscape modelling, Systems Analysis Modelling Simulation, Volume 41 (2001), pp. 755-776

[11] S. Grunwald; P. Barak 3D geographic reconstruction and visualization techniques applied to land resource management, Transactions in GIS, Volume 7 (2003), pp. 231-241

[12] A.W.J. Heijs; C.J. Ritsema; L.W. Dekker Three-dimensional visualization of preferential flow patterns in two soils, Geoderma, Volume 70 (1996), pp. 101-116

[13] R.J. Huggett Soil chronosequences, soil development, and soil evolution: a critical review, Catena, Volume 32 (1998), pp. 155-172

[14] H. Jenny Factors of Soil Formation, A system of quantitative pedology, McGraw-Hill, New York, USA, 1941 (109 p)

[15] F.A. Khan; T.E. Fenton Saturated zones and soil morphology in a Mollisol catena of Central Iowa, Soil Sci. Soc. Am. J., Volume 58 (1994), pp. 1457-1464

[16] Y. Lucas, Systèmes pédologiques en Amazonie brésilienne : équilibres, déséquilibres et transformations, PhD thesis, université de Poitiers, Poitiers, France, 1989.

[17] A.B. McBratney; M.L.M. Santos; B. Minasny On digital soil mapping, Geoderma, Volume 117 (2003), pp. 3-52

[18] R.M. McKenzie, Manganese oxides and hydroxides, in: J.B. and S.B Weed (Eds.), Minerals in soil environments. Eds Dixon, SSSA Book Series: 1, Madison, Wisconsin, USA., pp. 439–465.

[19] M. Neteler, Volume modelling of soils using GRASS GIS 3D tools, in: M. Brovelli, (Ed.), The Geomatics Workbook N. 2, Politecnico di Milano, Italy, 2001.

[20] S.J. Park; P.L.G. Vlek Environmental correlation of three-dimensional soil spatial variability: a comparison of three adaptive techniques, Geoderma, Volume 109 (2002), pp. 117-140

[21] V. Pereira; E.A. Fitzpatrick Three-dimensional representation of tubular horizons in sandy soils, Geoderma, Volume 81 (1998), pp. 259-303

[22] A.M. Petrovic; J.E. Siebert; P.E. Rieke Soil bulk-density analysis in three dimensions by computed tomographic scanning, Soil Sci. Soc. Am. J., Volume 46 (1982), pp. 445-450

[23] S. Salvador-Blanes, Déterminisme de la distribution spatiale des éléments majeurs et traces dans les sols en contexte métamorphique (Plateau d’Aigurande, Massif central, (France), PhD thesis, université François Rabelais, Tours, France, 2002.

[24] S. Salvador-Blanes; S. Cornu; D. King Distribution of chromium in soils developed into metamorphic rocks: variability at the hillslope scale, C. R. Geoscience, Volume 334 (2002), pp. 51-58

[25] S. Salvador-Blanes; S. Cornu; H. Bourennane; D. King Controls of the spatial variability of Cr concentration in topsoils of a central French landscape, Geoderma, Volume 132 (2006), pp. 143-157

[26] A. Samouelian; G. Richard; I. Cousin; R. Guérin; A. Bruand; A. Tabbagh Three-dimensional crack monitoring by electrical resistivity measurement, Eur. J. Soil Sci., Volume 55 (2004), pp. 751-762

[27] M.L.M. Santos; C. Guenat; M. Bouzelboudjen; F. Golay Three-dimensional GIS cartography applied to the study of the spatial variation of soil horizons in a Swiss floodplain, Geoderma, Volume 97 (2000), pp. 351-366

[28] P.M. Vasconcelos, T. Becker, P. Renne, G. Brimhall, K-Ar and Ar/ Ar dating of supergene K-bearing sulphates and Mn oxides; paleoclimatic implications from the time-stratigraphic evolution of deep weathering profiles, in: Anonymous Ed. 29th International Geological Congress, abstracts, 1992, p. 627.

[29] M. Zhang; A.D. Karathanasis Characterization of iron-manganese nodules in Kentucky Alfisols with perched water tables, Clays Clay Miner., Volume 45 (1997), pp. 428-439


Comments - Policy