Surface Spectroscopy of Oxisols, Entisols and Inceptisol and Relationships with Selected Soil Properties

ABSTRACT:

Traditional method of soil survey is expensive, slow, and must be carried out by experienced researchers. Thus, advances in soil observation technologies and the need to obtain information quickly by modern techniques have intensified the use of proximal sensing. This study characterized surface reflectance spectra (A horizon) and related them to traditional soil classification, based on morphological, physical, and chemical properties of representative pedogenetic profiles, developed in two toposequences of the Distrito Federal, Brazil. In the toposequences, 15 soil profiles were selected for a complete morphological description and sampling for laboratory analyses. Soil-landscape relationships were established, and profiles were classified to the fourth level of the Brazilian Soil Classification System (SiBCS). Classes of similar soils were grouped based on their surface spectra, resulting in spectral curves of 10 representative soils in the studied area. The Morphological Interpretation of Reflectance Spectrum (MIRS) and second derivative of the Kubelka-Munk (KM) function were applied to the soil spectra. The clustered soils were similar, mainly in terms of color, textural class, and organic matter content. Groups based on soil physical and chemical properties and on surface and subsurface colors were similar to those determined by surface reflectance. The identification of soil-landscape relationships was fundamental to understand the genesis and distribution of the soils, which had similar chemical and physical properties to their parent materials. The analysis of clusters based on soil surface reflectance proved efficient in determining groups of soil classes with similar properties. Surface reflectance data were related to the soil surface and subsurface properties determined by traditional soil sample analyses, since the two approaches formed similar groups. The simultaneous interaction of soil properties was assessed by MIRS analysis, while the second derivative of the KM function adequately quantified the mineralogy of the spectra.

Keywords:
toposequences; soil-landscape relationships; soil spectroscopy; spectral morphology; Kubelka-Munk

INTRODUCTION

Agriculture and natural resources are under severe pressure. The main causes are population growth, increased food consumption, and intensified use of agricultural land for bioenergy production (OECD/FAO, 2016). The negative impacts resulting from this evolution are numerous and can be related to the expansion (spreading of crops and pastures to new areas) and intensification (increased productivity and yield) of agriculture (Foley et al., 2011). According to the perspective of the Organization for Economic Co-operation and Development (OECD) and Food and Agriculture Organization of the United Nations (FAO) (2016), the agricultural cropland is expected to increase globally by 42 million ha in the next decade, and nearly half of this increase is estimated to occur in Brazil and Argentina.

In view of these changes, there is an urgent need to enhance the understanding of soil and identify the most appropriate for agriculture, which can only be done by mapping them in detail. However, the traditional method of soil mapping has proven to be expensive, time-consuming, and the task must be performed by experienced researchers. The technique provides detailed information on specific locations, but is limited with regard to number, volume, and spatial coverage of data (Adamchuk et al., 2017).

With the advances in observation technologies of the Earth, innovative equipment, user-friendly software, and open-source algorithms have been developed, opening new horizons for soil identification and mapping (Ben-Dor and Demattê, 2015).

Therefore, the need to obtain information rapidly by modern techniques has boosted the use of proximal sensing. The technique consists of analyzing information obtained from a soil sample without a direct contact. A soil sample is hit at light speed by electromagnetic radiation, consisting of different wavelengths. As each wavelength has a specific frequency and energy, some will have a specific affinity with soil properties. This interaction is followed by reflection, which varies between wavelengths, generating the so-called “spectral curve” or “spectral signature” of the target. Each sample is unique and has intrinsic properties represented in the resulting spectrum (Demattê et al., 2015a; Adamchuk et al., 2017).

The spectral bands commonly used for proximal sensing in soil science include visible (VIS: 350-780 nm), near-infrared (NIR: 780-1.100 nm), and short-wave infrared (SWIR: 1.100-2.500 nm) bands (Demattê et al., 2015a). Over the past 50 years, this technique proved to be quick, cheap, non-destructive, reproducible, repeatable, and environmentally friendly (Soriano-Disla et al., 2014).

Soil reflectance is a cumulative property derived from the intrinsically heterogeneous composition of soils, with regard to their mineralogical constitution, organic materials, particle distribution, and chemical composition (Demattê et al., 2012a). This composition is a result of pedogenetic factors and processes and can be represented by soil spectral curves. These signatures become references that are stored in soil spectral libraries, where they can be used for various soil studies (Nocita et al., 2015).

The determination of soil properties by means of reflectance spectroscopy techniques also detected significant correlations between the soil surface and subsurface (Mouazen et al., 2010; Stenberg et al., 2010; Nocita et al., 2015; Mohamed et al., 2017).

Soil classification can also be inferred from reflectance spectroscopy; therefore, authors such as Viscarra Rossel and Webster (2011) and Vasques et al. (2014) related reflectance spectra of all horizons of a profile with their respective classification. This procedure can rationalize processes, although soil subsurface data are required.

Along this line, some authors reported a close relationship between classified profiles and their surface spectra (Nanni et al., 2011, 2012; Cezar et al., 2013; Chicati et al., 2016; Zeng et al., 2017). In fact, the objective is not to classify, but to group similar soils. Such inferences can be useful for digital soil mapping. Therefore, the acquisition of soil spectra and understanding of their behavior in different regions of the country for the elaboration of so-called soil spectral libraries are of interest for the national and international community in Brazil, as stated by Bellinaso et al. (2010), and recently by Viscarra Rossel et al. (2016), with a view to a global spectral library.

As the spectrum of a soil sample is a signature comprising mineralogical information, it is expected to be useful for the grouping of similar soil profiles. In this setting, this study addresses the characterization of surface reflectance spectra (A horizon) and relates them with the respective soil classification by means of morphological, physical, and chemical properties of representative pedogenetic profiles, developed in two toposequences in the Distrito Federal, Brazil.

MATERIALS AND METHODS

Study site and physical environment

Two toposequences with representative soil distribution were selected for the study in the watersheds of the Estanislau Creek and Jardim Stream, both upstream of the water basin of the Jardim River, Distrito Federal (DF), Brazil (Sirgas 2000/UTM 23S, 216,623 - 226,924 m E / 8,261,507 – 8,244,566 m N), covering an area of 10,435 ha (Figure 1).

Figure 1
Watersheds of the Estanislau Creek and Jardim Stream, in the major Jardim River basin, Distrito Federal (DF) (IBGE, 2011), overlaid by a geological map (a) (Freitas-Silva and Campos, 1998), indicating the location of Toposequence 1 (Estanislau) (b) and Toposequence 2 (Alto Jardim) (c). P: profile, numbered in increasing order in relation to the position.

The Jardim River basin developed on rocks of the Paranoá, Canastra, and Bambuí Groups, which account, respectively, for circa 5, 25, and 70 % of the total area (Freitas-Silva and Campos, 1998). The geomorphological compartmentalization of almost the entire basin length is dominated by intermediate plateaus, with flat to gently undulating relief, low drainage density, less than 12 % slope, and elevation from 950 to 1,050 m (Semarh, 2012).

According to the semi-detailed soil mapping (scale 1:50,000) performed by Reatto et al. (2000), the main soils in the Jardim River basin are, according to Brazilian Soil Classification System (SiBCS) (Santos et al., 2013), Latossolos, Cambissolos, Plintossolos, Gleissolos, Argissolos, and Neossolos, which corresponding to Oxisols, Inceptsols, Petroferric Haplustox, Aquents, Ultissols, and Entisols in the American Soil Taxonomy System (ST) (Soil Survey Staff, 2014). The digital soil map generated by Poppiel (2016) shows FFc (Plintossolos Pétricos Concrecionários típicos/Petroferric Haplustox) as a frequent soil class in this basin. These soils are located under residual hills with concentric distribution, occurring throughout the area. They cover major upstream areas in the studied watersheds, together with relicts of the oldest geomorphological surface of the DF, the High Plateaus (called Chapadas Elevadas) (Pinto, 1994). These residual hills represent geomorphological surfaces resulting from the evolution of the Chapadas Elevadas (Morales and Assine, 2015), and their size decreases towards the downstream region of the Jardim River basin, on the eastern border of the DF. These residual relief forms resisted erosion and pedogenesis during the soil evolution of this geomorphological surface, which is active up to the present. Pelitic metasedimentary materials of the Canastra and/or Bambuí Groups, still preserved under FFc, are often observed. The soils evolved from FFc, due to differentiated conditions of erosion and pedogenesis, under the influence of the lithology of the Canastra and Bambuí Groups and local variations in relief features. Therefore, the occurrence of FFc is extremely important for the geomorphological and pedological development in the Jardim River basin, a fact given little attention in studies dealing with pedomorphogeological relations on the Central Plateau of Brazil, e.g., Motta et al. (2002), Barbosa et al. (2009), Campos et al. (2010), and Lacerda and Barbosa (2012), among others.

Field and laboratory activities

Pedomorphogeological relations were assessed in field observations in the watersheds, and two toposequences were selected (Figure 1a). Both had developed from predominantly pelitic lithologies of the Paranoá, Canastra, and Bambuí Groups, with a low metamorphic degree. The first was called Toposequence Estanislau, located in the watershed of the Estanislau Creek, between the highways DF-130 and DF-455, comprising the Rural Center of Tabatinga (Figure 1b). The second, named Alto Jardim, is located in the watershed of Jardim Stream, between the highways DF-455, DF-120, and DF-355 (Figure 1c).

Along the toposequences, 15 profiles (P) were morphologically described as proposed by Santos et al. (2015). Thirty samples were collected, of which 15 were from the A horizon and 15 from diagnostic subsurface horizons (B or C), for soil classification up to the fourth categorical level (subgroups) of the SiBCS (Santos et al., 2013) and the corresponding classification by the ST (Soil Survey Staff, 2014).

Air-dried fine earth (ADFE) fraction of sieved soil samples (<2 mm) was analyzed in the laboratory. To this end, ferruginous concretions (>2 mm) were eliminated from the surrounding soil. Physical and chemical analyses were performed according to methods described by Clasen (1997). Particle size composition was determined by the densimeter method (Bouyoucos, 1926). The chemical analyses determined: pH in water at a ratio of 1:2.5 v/v; potential acidity (H+Al) extracted by calcium acetate; organic matter (OM) determined by the Walkley-Black (1934) method; Ca2²⁺+Mg²⁺ and Al³⁺ extracted by KCl solution; K⁺ determined by flame spectrophotometry; and P₂O₅ assessed by extraction with Mehlich-1 (Vettori, 1969).

Spectroscopic analyses

The ADFE reflectance of surface A horizon, comprising 15 soil samples, was determined. The readings were performed in the laboratory, in a wavelength range from 350 to 2,500 nm (VIS-NIR-SWIR) using a FieldSpec® Pro sensor (Analytical Spectral Devices Inc., Boulder, CO, USA). Thereafter, the samples were placed on Petri dishes and homogenously distributed on the flat surface to perform the readings. Two halogen lights (50 W) were used as light source. The lights were installed at 90° to each other, at a distance of 0.35 m from the sample, and a zenithal angle of 30°, with a beam (uncollimated to the target). A fiber-optic cable, placed vertically at a distance of 0.08 m from the center of the sample surface, captured the energy reflected from an area of about 2 cm². For each sample, the mean reflectance was calculated from three repetitions in different positions, decreasing the shading effect. Each repetition consisted of 100 sensor readings, to maximize the signal-to-noise ratio. The instrument was calibrated at the beginning and then every 20 min, using a Spectralon white plate (Labsphere, North Sutton, NH, USA) as reference, with over 99 % reflectance.

Grouping soil profiles from surface reflectance

Soils with similar properties were grouped, according to the similarity of surface spectra (A horizons). This was possible because the soil surface patterns are related to soil class-specific variations in the underlying horizons and throughout the soil profile. Thus, it is possible to relate surface spectral curves with soil classes, as demonstrated by Demattê et al. (2007, 2009, 2016) and Nanni et al. (2011, 2012). However, these relationships should be evaluated with caution, because in some cases, the soil surface information alone is insufficient for soil classification.

The grouping was based on the Single Linkage method (Sneath and Sokal, 1973). This grouping technique of hierarchical agglomerative clustering is based on the Euclidean distance between all observations. The similarity between spectral behaviors of samples of the soil classes is inversely related to the Euclidean distance between the groups. The mean within each group was calculated, resulting in ten spectra for qualitative assessment.

Qualitative spectral evaluation

The spectral curves of surface soils were characterized qualitatively/descriptively by the Morphological Interpretation of Reflectance Spectrum (MIRS), as proposed by Demattê et al. (2014). This methodology considers aspects of the curve such as shape, intensity, and absorption features. To complement this technique, the second derivative of Savitzky-Golay was applied to the Kubelka-Munk (KM) function (Torrent and Barrón, 2008), using The Unscrambler^® software. This mathematical procedure amplifies small reflectance (positive values) or absorption features (negative values), allowing quantification of the materials (Sellitto et al., 2009).

RESULTS AND DISCUSSION

Classification of soils and physical environment

Based on the morphological, physical, and chemical properties, the soils from Estanislau Toposequence (T1) were classified as: P1-LVd (Latossolo Vermelho Distrófico típico/Rhodic Acrustox); P2-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Acrustox); P3-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox), developed under unit MNPccf (Chlorite Carbonate Phyllite) of the Canastra Group; P4-LVd (Latossolo Vermelho Distrófico típico/Rhodic Acrustox); P5-LVAdc (Latossolo Vermelho-Amarelo Distrófico petroplíntico/Petroferric Acrustox); P6-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Acrustox), originated from unit NPbT (Argillaceous Metasilite/Arkoses) of the Bambuí Group (Tables 1 and 3, Figure 2a).

Table 1
Description of morphological properties of profiles in Toposequence 1 - Estanislau, Distrito Federal, Brazil

Hor.(1)	Depth	Munsell Color (2)	Structure(3)	Consistency(4)
				Dry	Moist	Wet
	m
P1-LVd (Latossolo Vermelho Distrófico típico/Typic RhodRhodic Acrustox)
A	0.00-0.15	2.5YR 4/4	st. vsm. cb.	S.	VFr.	SPL.	SST.
AB	0.15-0.25	2.5YR 4/6	st. vsm. gr.	S.	VFr.	PL.	ST.
BA	0.25-0.40	2.5YR 5/6	st. vsm. gr.	S.	VFr.	PL.	ST.
Bw	0.40+	2.5YR 5/8	st. vsm. gr.	S.	VFr.	PL.	ST.
P2-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Acrustox)
Ac	0.00-0.20	7.5YR 4/4 var. 2.5YR 5/8	st. vsm. cb.	S.	VFr.	PL.	ST.
AF	0.20-0.35	7.5YR 5/4 var. 2.5YR 5/8	st. vsm. gr.	S.	VFr.	PL.	SST.
FA	0.35-0.50	7.5YR 7/8 var. 2.5YR 5/8	st. vsm. gr.	S.	VFr.	PL.	SST.
F	0.50+	5YR 6/8 var. 2.5YR 5/6	st. vsm. gr.	S.	VFr.	PL.	SST.
P3-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox)
Ac	0.00-0.20	10YR 3/6 var. 10YR 5/6	st. vsm. cb.	S.	VFr.	PL.	ST.
F	0.20-0.60	10YR 6/8	st. vsm. gr.	S.	VFr.	PL.	ST.
C	0.60+	10R 5/8 var. 2.5Y 7/4	wk. sm. to lg. bk.	H.	Fi.	PL.	ST.
P4-LVd (Latossolo Vermelho Distrófico típico/Rhodic Acrustox)
A	0.00-0.20	10R 4/3	st. vsm. cb.	S.	Fr.	PL.	ST.
AB	0.20-0.30	10R 4/6	st. vsm. gr.	S.	Fr.	PL.	ST.
BA	0.30-0.51	10R 4/4	st. vsm. gr.	S.	Fr.	PL.	ST.
Bw	0.51+	10R 5/8	st. vsm. gr.	S.	VFr.	PL.	ST.
P5-LVAdc (Latossolo Vermelho-Amarelo Distrófico petroplíntico/Petroferric Acrustox)
A	0.00-0.15	10YR 3/4	st. vsm. gr.	S.	Fr.	PL.	SST.
AB	0.15-0.25	10YR 4/6	st. vsm. gr.	S.	Fr.	PL.	SST.
BA	0.25-0.35	10YR 5/6	st. vsm. gr.	S.	Fr.	PL.	SST.
Bw	0.35-0.70	10YR 6/6	st. vsm. gr.	S.	VFr.	PL.	SST.
F	0.70+	10R 5/8	st. vsm. gr.	S.	Fr.	PL.	SST.
P6-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Acrustox)
Ac	0.00-0.25	7.5YR 4/6	st. vsm. cb.	S.	Fr.	PL.	ST.
F1	0.25-0.75	7.5YR 5/8	st. vsm. gr.	S.	Fr.	PL.	ST.
F2	0.75-1.75	10R 5/8	st. vsm. gr.	S.	Fr.	PL.	ST.
C	1.75+	10R 6/8 var. 7.5YR 7/4	wk. vsm. bk.	H.	Fr.	PL.	ST.

• Hor.= horizon; Var. = variegated.
• ⁽¹⁾
  Horizon suffixes: w = oxic horizon, c = concretionary horizon; master horizons: A = ochric epipedon, F = petroferric contact.
• ⁽²⁾
  Moist color.
• ⁽³⁾
  Structure: st.= strong, wk.= weak, vsm.= very small, sm.= small, lg.= large, cb.= crumb, gr.= granular, bk.= blocky.
• ⁽⁴⁾
  Consistency [determined in air-dried, sieved soil samples (<2 mm)]: S.= soft, H.= hard, Fr.= friable, Fi.= firm, VFr.= very friable, PL.= plastic, ST.= sticky, SPL.= slightly plastic, SST.= slightly sticky.

Table 2
Description of morphological properties of profiles in Toposequence 2 - Alto Jardim, Distrito Federal, Brazil

Hor.(1)	Depth	Munsell Color(2)	Structure(3)	Consistency(4)
				Dry	Moist	Wet
	m
P1-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox)
Ac	0.00-0.18	10YR 3/4	st. vms. gr.	S.	Fr.	PL.	ST.
AF	0.18-0.45	10YR 4/6	st. vms. gr.	S.	Fr.	PL.	ST.
F	0.45-1.20	10YR 6/8	st. vms. gr.	S.	Fr.	PL.	ST.
FC	1.20-1.60	2.5YR 5/8	st. vms. gr.	S.	Fr.	PL.	ST.
C	1.60+	2.5YR 3/6	wk. vms. bk.	SH.	Fi.	SPL.	ST.
P2-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox)
Ac	0.00-0.30	10YR 4/4	st. vms. cb.	L.	Fr.	PL.	ST.
AF	0.30-0.48	10YR 3/4	st. vms. gr.	S.	Fr.	PL.	ST.
F	0.48-1.05	10YR 4/6	st. vms. gr.	S.	Fr.	PL.	ST.
FC	1.05-1.50	5YR 5/6 var. 2.5Y 8/2	st. vms. gr.	H.	Fr.	PL.	ST.
C	1.50+		wk. vms. bk.	H.	Fr.	PL.	ST.
P3-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox)
Ac	0.00-0.20	7.5YR 3/4	st. vms. cb.	S.	Fr.	PL.	SST.
AF	0.20-0.40	7.5YR 4/6	st. vms. gr.	S.	Fr.	PL.	SST.
F	0.40-0.96	7.5YR 4/4	st. vms. gr.	S.	Fr.	PL.	ST.
FC	0.96+	2.5YR 6/6	wk. vms. bk.	H.	Fi.	PL.	ST.
P4-CXbd (Cambissolo Háplico Tb Distrófico típico/Oxic Dystrustepts)
A	0.00-0.10	7.5YR 4/6	wk. vms. bk.	H.	Fi.	PL.	ST.
Bi	0.10-0.20	7.5YR 5/6	wk. vms. bk.	H.	Fi.	PL.	ST.
BC	0.20-0.25	5YR 5/6 var. 2.5YR 4/8	wk. vms. bk.	H.	Fi.	PL.	ST.
C	0.25+	2.5YR 6/4	wk. vms. bk.	VH.	Fi.	PL.	ST.
P5-RRd (Neossolo Regolítico Distrófico típico/Typic Ustorthent)
A	0.00-0.08	10YR 5/4	wk. sm. bk.	H.	Fi.	PL.	ST.
AC	0.08-0.14	10YR 5/3	wk. sm. bk.	H.	Fi.	PL.	ST.
C	0.14+	10R 3/6 var. 7.5YR 8/0	wk. sm. bk.	VH.	Fi.	PL.	ST.
P6-RRd (Neossolo Regolítico Distrófico típico/Typic Ustorthent)
A	0.00-0.10	10YR 5/4	st. sm. bk.	H.	Fi.	PL.	ST.
AC	0.10-0.22	10YR 5/6	st. sm. bk.	H.	Fi.	PL.	ST.
C	0.22+	2.5YR 3/6 var. 2.5Y 7/2	wk. sm. bk.
P7-RRd (Neossolo Regolítico Distrófico típico/Typic Ustorthent)
A	0.00-0.3	5YR 5/5	mo. sm. bk.	H.	Fi.	PL.	ST.
C	0.3-0.70+	10R 4/8	wk. sm. bk.	H.	Fi.	PL.	ST.
P8-GXbd (Gleissolo Háplico Tb Distrófico típico/Typic Fluvaquents)
A	0.00-0.33	10YR 5/2	wk. sm. bk.	VH.	Fi.	PL.	ST.
ACg	0.33-0.40	10YR 6/2 var. 5YR 6/8	wk. sm. bk.	VH.	Fi.	PL.	ST.
Cg	0.40-0.70+	10YR 7/2 var. 5YR 6/8	ma.	VH.	Fi.	PL.	ST.
P9-LVd (Latossolo Vermelho Distrófico típico/Rhodic Acrustox)
A	0.00-0.20	10R 3/6	st. vms. cb.	SH.	Fr.	PL.	ST.
ABw	0.20-0.35	10R 4/6	st. vms. cb.	SH.	Fr.	PL.	ST.
BAw	0.35-0.55	10R 4/8	st. vms. gr.	SH.	Fr.	PL.	ST.
Bw	0.55+	10R 5/8	st. vms. gr.	S.	VFr.	PL.	ST.

• Hor.= horizon; Var. = variegated.
• ⁽¹⁾
  Horizon suffixes: w = oxic horizon, c = concretionary horizon, I = cambic horizon, g = redoximorphic features; master horizons: A = ochric epipedon, F = petroferric contact.
• ⁽²⁾
  Moist color.
• ⁽³⁾
  Structure: st.= strong, mo.= moderate, wk.= weak, vsm.= very small, sm.= small, cb.= crumb, gr.= granular, bk= blocky, ma.= massive.
• ⁽⁴⁾
  Consistency [determined in air-dried, sieved soil samples (<2 mm)]: L.= loose, S.= soft, H.= hard, SH.= slightly hard, VH.= vary hard, Fr.= friable, Fi.= firm, VFr.= very friable, PL.= plastic, ST.= sticky, SPL.= slightly plastic, SST.= slightly sticky.

Table 3
Physical and chemical properties of profiles of Toposequence 1 - Estanislau, Distrito Federal, Brazil

Hor	Clay	Silt	Sand	Textural class	B/A	Si/Cl	pH(H2O)	SB	t	T	Ta	V	m	OM	OC
	---------- g kg-1 ----------							------------ cmolc dm-3 ------------				-------- % --------		------ g kg-1 ------
P1-LVd (Latossolo Vermelho Distrófico típico/Typic Rhodustox)
A	792	87	121	vcly.		0.1	5.0	1.0	1.4	7.9	9.9	13.0	29.0	27.0	16.0
Bw	785	107	108	vcly.	1.0	0.1	5.3	0.4	0.5	2.4	3.1	17.0	20.0	11.0	7.0
P2-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox)
Ac	605	130	265	vcly.		0.2	4.8	1.9	3.0	11.1	18.3	17.0	37.0	34.0	19.0
F	675	65	260	vcly.	1.1	0.1	5.4	0.7	0.8	4.5	6.6	16.0	12.0	20.0	12.0
P3-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox)
Ac	607	213	180	vcly.		0.4	5.2	1.5	2.5	9.1	14.9	16.0	40.0	29.0	16.0
F	587	246	167	cly.	1.0	0.4	4.7	0.8	1.0	4.6	7.8	17.0	20.0	20.0	12.0
P4-LVd (Latossolo Vermelho Distrófico típico/Typic Rhodustox)
A	755	74	171	vcly.		0.1	5.1	1.4	2.1	9.8	13.0	14.0	33.0	36.0	21.0
Bw	773	91	136	vcly.	1.0	0.1	4.9	0.5	0.6	3.8	4.9	13.0	17.0	20.0	11.0
P5-LVAdc (Latossolo Vermelho-Amarelo Distrófico petroplíntico/Plinthic Haplustox)
A	692	133	175	vcly.		0.2	5.1	1.7	2.1	9.1	13.2	19.0	19.0	31.0	18.0
Bw	725	90	185	vcly.	1.1	0.1	5.4	0.5	0.6	2.6	3.6	19.0	17.0	17.0	10.0
P6-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox)
Ac	623	89	289	vcly.		0.1	5.0	1.3	1.7	8.1	13.0	16.0	24.0	28.0	16.0
F1	702	89	209	vcly.	1.1	0.1	5.1	0.7	0.8	3.2	4.5	22.0	12.0	17.0	10.0

• Hor. = horizon. Clay, silt, and sand determined by densimeter method. Textural soil class according to Santos et al. (2015): cly.= clayey, vcly.= very clayey. B/A = B horizon/A horizon relationship; Si/Cl = silt/clay ratio. pH(H₂O): pH in water-saturated soil paste at a ratio of 1:2.5 v/v; K⁺ (Mehlich 1); Ca²⁺, Mg²⁺, and Al³⁺ (KCl 1 mol L^-1); H+Al (Calcium acetate 0.5 mol L^-1 - pH 7.0). SB = sum of bases (Ca²⁺+Mg²⁺+K⁺); t = effective CEC = SB+Al³⁺; T = Cation exchange capacity (SB+H⁺+Al³⁺); Ta = T clay activity = (T × 1000/clay content); V = base saturation (100 × S/T); m = aluminum saturation (Al³⁺ × 100/t); OC = organic carbon (Walkley-Black, 1934); OM = OC × 1.72.

Figure 2
Topographic profile of toposequences Estanislau-T1 (a) and Alto Jardim-T2 (b), and location of the studied soil profiles (P). I: infiltration; R: runoff. T1: P1-LVd (Latossolo Vermelho Distrófico típico/Rhodic Acrustox); P2-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Acrustox); P3-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox), developed under unit MNPccf (Chlorite Carbonate Phyllite) of the Canastra Group; P4-LVd (Latossolo Vermelho Distrófico típico/Rhodic Acrustox); P5-LVAdc (Latossolo Vermelho-Amarelo Distrófico petroplíntico/Petroferric Acrustox); P6-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Acrustox), originated from unit NPbT (Argillaceous Metasilite/Arkoses) of the Bambuí Group. T2: P1-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox); P2-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox); P3-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox); P4-CXbd (Cambissolo Háplico Tb Distrófico típico/Oxic Dystrustepts); P5-RRd (Neossolo Regolítico Distrófico típico/Typic Ustorthent); P6-RRd (Neossolo Regolítico Distrófico típico/Typic Ustorthent); P7-RRd (Neossolo Regolítico Distrófico típico/Typic Ustorthent); P8-GXbd (Gleissolo Háplico Tb Distrófico típico/Typic Fluvaquents); P9-LVd (Latossolo Vermelho Distrófico típico/Rhodic Acrustox), all originated from unit NPbT (Argillaceous Metasilite/Arkoses) of the Bambuí Group.

The soils of the toposequence Alto Jardim (T2) were classified as: P1-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox); P2-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox); P3-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox); P4-CXbd (Cambissolo Háplico Tb Distrófico típico/Oxic Dystrustepts); P5-RRd (Neossolo Regolítico Distrófico típico/Typic Ustorthent); P6-RRd (Neossolo Regolítico Distrófico típico/Typic Ustorthent); P7-RRd (Neossolo Regolítico Distrófico típico/Typic Ustorthent); P8-GXbd (Gleissolo Háplico Tb Distrófico típico/Typic Fluvaquents); P9-LVd (Latossolo Vermelho Distrófico típico/Rhodic Acrustox), all originated from unit NPbT (Argillaceous Metasilite/Arkoses) of the Bambuí Group (Tables 2 and 4, Figure 2b).

Table 4
Physical and chemical properties of profiles of Toposequence 2 - Alto Jardim, Distrito Federal, Brazil

Hor.	Clay	Silt	Sand	Textura class(1)	B/A	Si/Cl	pH(H2O)	S	t	T	Ta	V	m	OM	OC
	-------- g kg-1 --------							---------- cmolc dm-3 ----------				------ % ------		----- g kg-1 -----
P1-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox)
Ac	500	280	220	cly.		0.6	5.1	3.0	3.2	8.4	16.8	36.0	6.0	32.0	18.0
F	651	263	85	vcly.	-	0.4	5.1	0.9	1.3	4.9	7.5	18.0	31.0	21.0	12.0
P2-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox)
Ac	436	453	111	si-cl.		1.0	4.4	0.7	1.6	7.0	16.1	10.0	56.0	28.0	16.0
F	621	312	67	vcly.	-	0.5	4.9	1.0	1.6	4.1	6.5	24.0	38.0	20.0	11.0
P3-FFc (Plintossolo Pétrico Concrecionário típico/Petroferric Haplustox)
Ac	550	381	69	cly.		0.7	5.0	1.8	2.4	8.4	15.2	21.0	25.0	33.0	19.0
F	583	374	43	cly.	-	0.6	5.3	1.1	2.0	7.2	12.3	15.0	45.0	26.0	15.0
P4-CXbd (Cambissolo Háplico Tb Distrófico típico/Oxic Dystrustepts)
A	502	417	82	si-cl.		0.8	4.9	1.3	2.3	6.7	13.3	19.0	43.0	31.0	18.0
Bi	514	397	89	cly.	1.0	0.8	4.6	0.9	2.2	6.3	12.2	14.0	59.0	22.0	13.0
P5-RRd (Neossolo Regolítico Distrófico típico/Typic Ustorthent)
A	493	358	148	cly.		0.7	5.3	3.9	5.0	9.7	19.7	40.0	22.0	33.0	19.0
C	435	422	144	si-cl.	-	1.0	5.0	2.1	8.6	14.5	33.3	14.0	76.0	15.0	9.0
P6-RRd (Neossolo Regolítico Distrófico típico/Typic Ustorthent)
A	489	412	99	si-cl.		0.8	5.0	2.3	4.0	10.1	20.7	23.0	42.0	31.0	18.0
C	528	435	37	si-cl.	-	0.8	4.9	0.7	3.7	6.5	12.3	11.0	81.0	8.0	4.0
P7-RRd (Neossolo Regolítico Distrófico típico/Typic Ustorthent)
A	422	382	197	cly.		0.9	5.0	3.0	4.1	8.8	20.9	34.0	27.0	28.0	16.0
C	415	365	219	cly.	-	0.9	4.7	1.8	4.7	7.6	18.4	24.0	62.0	11.0	7.0
P8-GXbd (Gleissolo Háplico Tb Distrófico típico/Typic Fluvaquents)
A	402	414	185	si-cl.		1.0	5.4	3.9	4.2	7.4	18.6	51.0	7.0	22.0	13.0
Cg	435	422	143	si-cl.	-	1.0	4.8	1.1	2.9	5.7	13.1	19.0	62.0	18.0	11.0
P9-LVd (Latossolo Vermelho Distrófico típico/Rhodic Acrustox)
A	684	248	67	vcly.		0.4	4.7	1.3	1.8	6.2	9.1	19.0	31.0	28.0	16.0
Bw	702	225	73	vcly.	1.0	0.3	4.8	1.0	1.1	2.3	3.3	43.0	9.0	17.0	10.0

• Hor. = horizon. Clay, silt, and sand determined by densimeter method. Textural soil class according to Santos et al. (2015): cly.= clayey, vcly.= very clayey. B/A = B horizon/A horizon relationship; Si/Cl = silt/clay ratio. pH(H₂O): pH in water-saturated soil paste at a ratio of 1:2.5 v/v; K⁺ (Mehlich 1); Ca²⁺, Mg²⁺, and Al³⁺ (KCl 1 mol L^-1); H+Al (Calcium acetate 0.5 mol L^-1 - pH 7.0). SB = sum of bases (Ca²⁺+Mg²⁺+K⁺); t = effective CEC = SB+Al³⁺; T = Cation exchange capacity (SB+H⁺+Al³⁺); Ta = T clay activity = (T × 1000/clay content); V = base saturation (100 × S/T); m = aluminum saturation (Al³⁺ × 100/t); OC = organic Carbon (Walkley-Black, 1934); OM = OC × 1.72.

In the upper and middle thirds of the T1 slopes (Figure 2a), convex-shaped ramps were identified, evolving into a rectilinear ramp in the lower third, incorporated into the drainage network and without floodplain. Vertical water flow prevails across the toposequence, due to the high hydraulic conductivity of the granular and crumb structure of these soils under native vegetation (Campos et al., 2010). The high permeability of the Bw horizon of LV is due to structures and types of associated pores that prevent rapid water circulation (Cunha et al., 2008). In the Mar de Morros region (state of Minas Gerais), with similar relief features, Ippoliti et al. (2005) reported deep soils without drainage limitations. In P2 and P3 (FFc), accumulation of ferruginous concretions was observed, hampering free water infiltration into the profile, originating a wetter pedoenvironment (Motta et al., 2002). According to Barbosa et al. (2009) and Campos et al. (2010), the formation of a concretionary horizon in LVA of DF, as in P5, is related to its distribution and geomorphological evolution from the Chapadas Elevadas. In the subsurface, this horizon developed where the slope increases slightly, with occurrences of water table oscillations, thus resulting in the stabilization of goethite. Currently, near the edges of the Chapadas Elevadas, this horizon crops out in layers of petroplinthites (iron-cemented concretions), constituting the FFc of the study site.

The relief of the T2 profile (Figure 2b), from the top toward the drainage course (Jardim Stream) has convex-concave features, with a pronounced slope. Therefore, the soils evolved from a convex feature, from the initial top to a distance of 500 m, to a concave feature in the mid-segment of the western hillside. From the drainage bed, convex-rectilinear relief features (extensive ramp) are observed from about 1,800 m to the last toposequence profile. The occurrence of FFc (P1, P2, and P3 of T2) at the convex top, similar to those in P2, P3, and P6 of T1, represent relicts of the subsurface concretionary horizon of the soils from the first geomorphological surface of the Chapadas Elevadas, where petroplinthites of the F horizon (petroferric contact) hinder pedogenesis of the parent material (Motta et al., 2002). In the sloping areas of concave portions, the beginning of pedogenesis of the parent material is observed, with the development of young soils, such as CXbd (P4) and RRd (P5, P6, and P7), where runoff and land dissection are more intense. In the nearby drainage areas with water accumulation however, GXbd (P8) develops. On the rectilinear eastern side of T2, LVd (P9) occurred, forming an extensive ramp with a flat to gently rugged relief that favors water infiltration into the soil.

In studies on pedomorphogeological relations and distribution of pedomorphs in the DF, Lacerda and Barbosa (2012) observed the presence of Cambissolo Háplico (CX) in the steepest portions of the terrain (slope ≥45 %) and Gleissolo Háplico (GX) on the floodplain with slopes between 0 to 2 %. In the region of Santa Maria, Rio Grande do Sul, Araújo Pedron et al. (2011) described Neossolos Regolíticos (RR) under a hilly to strongly hilly relief, with a low soil water infiltration rate, due to the higher density of the C horizon.

Relief shapes play a decisive role in water dynamics in the soil profile, affecting the pedogenetic processes and consequently formation of different soils (Vidal-Torrado et al., 2005). Thus, the participation of relief features in the evolution of soils of the studied site is determinant, for influencing the amount of water infiltrated into the soil, accelerating chemical reactions of weathering, and promoting the transport of solids or materials in solution, producing effects that result in different soil types in the various positions of the toposequences (Pennock and Veldkamp, 2006).

The occurrence of LV, generally humic, highly evolved, and deep, was reported by Chaves (2016) in the top regions of the Serra da Mantiqueira, preserved in convex and mostly steep slopes (<45 %), related to the regional geomorphological and geological evolution, due to the crustal uplift that originated the mountain chains Serra da Mantiqueira and Serra do Mar, preserving the Paleo-Latossolos (Paleo-Oxisols) developed on the geomorphic surfaces Japi and Sul Americana.

As a general rule, deeper and more evolved soils generally occur on convex slopes, while shallower (saprolitics) and less evolved soils tend to occur on concave slopes (Chagas et al., 2013). The shape (curvature) and orientation (solar radiation incidence) of slopes are important in pedogenesis, for directly affecting the air and soil temperature level, evapotranspiration and soil moisture, as well as other important soil properties (Senthilkumar et al., 2009; Camargo et al., 2010; Oliveira et al., 2013).

Morphological aspects

The morphological characterization of soils (Tables 1 and 2) differentiated each soil class, mainly in terms of color, which clearly shows the action of different pedogenetic processes among soils. The diagnostic horizon of the LVd had a reddish color, while those of LVAdc, FFc, CXbd, and RRd were yellowish. According to Schwertmann (1993) and Gomes et al. (2004), reddish colors are due to higher hematite contents (good drainage conditions), while yellowish colors indicate the dominance of goethite over hematite (wetter pedoenvironment). These colors varied with depth, and were darker at the soil surface because of the influence of OM (Schulze et al., 1993) and redder in the C horizon, due to Fe₂O₃-rich pelitic parent material (Freitas-Silva and Campos, 1998). The grayish color of GXbd is mainly due to the presence of Fe in reduced form (Fe²⁺), indicating a hydromorphic environment (Ribeiro et al., 2012).

In the SiBCS (Santos et al., 2013), color (Munsell Color Chart) is used as an important property of soil classification, from which information about the soil can be inferred. Color variations may be a response to alterations in relief, vegetation, depth, weather, aeration, parent material, weathering degree, mineralogy, and OM content (Schulze et al., 1993; Schwertmann, 1993), linked to the origin of these soils (Lacerda et al., 2008).

In all profiles, the surface diagnostic horizon was moderate (ochric), except for T2P7-RRd, a young soil with lowest thickness (0.03 m) of the A horizon (weak). The LVd and LVAdc were classified as Latossolos (highly weathered Oxisols, mostly classified as Ustox in the ST) due to the presence of a latossólico (oxic) diagnostic B horizon and the observed dominant color (hue 2.5YR or redder for LV and 10YR for LVA). In the FFc, ferruginous concretions were observed in horizons Ac, ABc, and BAc, indicating a concretionary horizon (petroferric contact) in the diagnostic position (F) with variegated color (from 10YR to 10R). They were classified as typical, since no other differentiating property was observed.

The RRd soils contained an A horizon overlying the C horizon and had a maximum depth of 0.22 m, without no other diagnostic horizon. The CXbd had an incipient (cambic) diagnostic B horizon with a depth of 0.25 m and hue 7.5YR. In turn, GXbd had a thick surface A horizon (0.33 m) and transitional gleyed horizon (ACg), overlying the greyish gleyed Cg in diagnostic position. The redoximorphic environment of this horizon is influenced by the high water table or slowly permeable layer in the profile (Santos et al., 2013).

In the LVd, LVAdc, and FFc, strong, very small granular and crumb aggregates prevail in the surface and subsurface horizons. A granular aggregate structure is favored by the presence of Al and Fe oxides, typical of these soil classes in the Cerrado region (tropical savanna on the Central Brazilian Plateau) (Gomes et al., 2004). In turn, the occurrence of LVAdc and FFc is related to a wetter pedoenvironment, due to the proximity of petroferric material from the plateau edges or to groundwater level fluctuations. These conditions are more conducive to the formation of goethite than of hematite, and of plinthite in the subsurface, precursor of ferruginous concretions (Motta et al., 2002). In young soils, such as RRd, CXbd, and GXbd, the structure is variable, with predominance of weak small aggregates in subangular blocks.

Physical and chemical properties

Clay content ranged between 402 and 792 g kg^-1 (Tables 3 and 4) because of the parent material, pre-weathered during genesis (Freitas-Silva and Campos, 1998). The soil textural class ranged from very clayey to clayey and silty clay (Santos et al., 2015). No texture gradients that characterized a textural (argillic) B horizon in soils were observed.

The highest silt content of P3 was recorded in T1 (Table 3), possibly due to the lower thickness of the F horizon (0.20-0.60 m), influenced by the C horizon, which still preserves features of the Canastra Group (Freitas-Silva and Campos, 1998). The shallow depth of this soil, along with the slope of the local relief, contributed to reduce water infiltration into the profile. Therefore, the volume of available soil and the water level mitigated soil weathering by slowing down pedogenetic processes (Kämpf and Curi, 2012).

The silt/clay ratio was very low in LVd, LVAdc, and FFc, due to the high weathering degree of these soils (Santos et al., 2013). In CXbd, RRd, and GXbd, however, the silt/clay ratios exceeded 0.6, indicating a lower evolution degree of these soils. Similar physical properties were described by Barbosa et al. (2009) and Pereira et al. (2010) for soils evolved from pelitic lithology on the Brazilian Central Plateau.

Regarding the chemical composition (Tables 3 and 4), the soils were considered acidic to highly acidic (Santos et al., 2013). The sum of bases (S), cation exchange capacity (CEC), and base saturation (V%) were very low, featuring dystrophic soils. The contents of OM (organic matter) and OC (organic carbon) were high in most soils, as they occurred in areas of native vegetation, confirming observations of Barbosa et al. (2009), Campos et al. (2010), and Lacerda and Barbosa (2012).

Grouping soils based on surface spectra

Cluster analysis, based on the Euclidean distance (Figure 3a) and linkage at each step of the algorithm (Figure 3b), applied to the reflectance data VIS-NIR-SWIR (Figure 3c), allowed the distribution of 15 profiles in 10 groups, representing the soils of the study site. The clustering criterion using linkage distances was based on the identification of a plateau in the vertical direction, relatively large compared to the others (Figure 3b).

Figure 3
Dendrogram of cluster analysis (a) and linkage distance chart indicating the distance of dendrogram cut-offs for the definition of groups (b), based on soil surface spectra (c).

The dendrogram obtained from surface spectral data (Figure 3a) shows the formation of 10 groups constituted by similarities between soils of various classes (Table 5), e.g., of wet color (Munsell), textural class, and OM content. Spectral data from 350 to 2,500 nm provide a reliable measure of soil properties, which can be used in qualitative and quantitative soil analyses (Viscarra Rossel et al., 2016).

Table 5
Main soil properties grouped by surface reflectance

Soil Group(1)	Munsell Color (moist soil)	Textural Class(2)	OM
			g kg-1
LV (Latossolos Vermelhos/Rhodic Acrustox)
LV-I	10R 3/6	Very clayey	28
LV-II	10R 4/3	Very clayey	31.5
LVA (Latossolos Vermelho-Amarelos/Petroferric Acrustox)
LVA	10YR 3/4	Very clayey	31
FF (Plintossolos Pétricos/Petroferric Ustox)
FF-I	10YR 3/6	Clayey	29
FF-II	10YR 3/4	Clayey	31
FF-III	7.5YR 4/4	Very clayey	34
RR (NeossoloS Regolíticos/Typic Ustorthent)
RR-I	10YR 5/4	Silty clay	32
RR-II	5YR 5/5	Clayey	28
CX (Cambissolo Háplico/Oxic Dystrustepts)
CX	7.5YR 4/6	Clayey	31
GX (Gleissolo Háplico/Typic Fluvaquents)
GX	10YR 5/2	Silty clayey	22

• ⁽¹⁾
  Clustered soils with similar properties based on surface reflectance.
• ⁽²⁾
  Textural soil class according to Santos et al. (2015). OM = Organic Matter (determined by Walkley-Black method).

The groups obtained from physical and chemical data of the soil surface and subsurface (Figure 4a) differed slightly from those based on surface reflectance data (Figure 3a). This discrepancy was possibly due to information related to soil mineralogy, which was not included in this analysis. Therefore, the dendrogram grouped basically profiles with similar texture and OM contents.

Figure 4
Dendrograms of the cluster analysis based on surface and subsurface data: physical and chemical properties (a) and soil color (b). T: toposequence; P: profile. ^* Groups obtained from surface reflectance data (Figure 3a).

In the clustering based on surface and subsurface soil color (Figure 3b), the groups were very similar to those based on surface reflectance. The reason is that color is a result of the visible light interaction (350-780 nm) with the different soil properties, carrying this information (Schulze et al., 1993; Schwertmann, 1993). This agrees with the spectroscopic data, which include wavelengths of NIR and SWIR, in addition to VIS, and are strongly related with various soil data, including soil mineralogy (Demattê et al., 2004; Brown et al., 2006; Viscarra Rossel et al., 2009; Soriano-Disla et al., 2014). It is noteworthy that the surface reflectance data were related to the soil surface and subsurface properties determined by traditional analyses of soil samples, since the two approaches formed similar groups.

According to Zeng et al. (2016), the use of spectroscopic data of the soil surface horizons can achieve an accuracy of more than 75 % in individual soil classification. This agrees with Xie et al. (2015), in which five soil classes were discriminated by surface reflectance with 92.2 % accuracy. According to Nanni et al. (2012), using spectral information of the soil surface as indicators, soil classes were discriminated with an overall accuracy between 52 and 70 %.

Interpretation of surface spectra

The spectra obtained from the cluster analysis (Figures 5 to 10) presented spectral features of Fe oxides such as goethite (480 and 950 nm) and hematite (530 and 850 nm) (Sherman and Waite, 1985), kaolinite (1,400 and 2,205 nm), gibbsite (2,265 nm) (Madeira et al., 1995), 2:1 clay minerals, and water adsorbed on the particle surface (1,400 and 1,900 nm) (Ben-Dor et al., 2008; Stenberg et al., 2010).

Figure 5
Soil spectra and diagram of the second derivative of Kubelka-Munk (KM) of the studied Latossolo Vermelho- LV of watersheds of the Estanislau Creek and Jardim Stream, Distrito Federal (DF), Brazil. Soil clusters: LV-I; LV-II. A and B are derived curve regions.

Figure 6
Soil spectrum and diagram of the second derivative of Kubelka-Munk of the studied Latossolo Vermelho-Amarelo (LVA), representative of the LVA of watersheds of the Estanislau Creek and Jardim Stream, Distrito Federal, Brazil. A and B are derived curve regions.

Figure 7
Soil spectra and diagrams of the second derivative of Kubelka-Munk of the Plintossolo Pétrico Concrecionário (FFc) studied in the watersheds of the Estanislau Creek and Jardim Stream, Distrito Federal, Brazil. Soil clusters: FF-I; FF-II; FF-III. A and B are derived curve regions.

Figure 8
Soil spectra and diagrams of the second derivative from Kubelka-Munk (KM) of Neossolos Regolíticos (RR) in the watersheds of Estanislau Creek and Jardim Stream, Distrito Federal, Brazil. Soil clusters: RR-I; RR-II. A and B are derived curve regions.

Figure 9
Soil spectrum and diagram of the second derivative of Kubelka-Munk of the Cambissolo Háplico (CX) studied, representative of CX in the watersheds of the Estanislau Creek and Jardim Stream, Distrito Federal, Brazil. A and B are derived curve regions.

Figure 10
Soil spectrum and diagram of the second derivative of Kubelka-Munk (KM) of the Gleissolo Háplico (GX) studied, representative of the GXs of the watersheds of the Estanislau Creek and Jardim Stream, Distrito Federal (DF), Brazil. A and B are derived curve regions.

In the visible (VIS) range, the second derivative of the KM function (Figures 5 to 10) showed reflectance peaks (positive values) for goethite (450 nm) and hematite (580 nm), as well as absorption feature (negative values) at 480 nm (goethite) and 530 nm (hematite) (Sellitto et al., 2009). In mixtures of goethite and hematite, the bands overlap, resulting in a peak located between the two minerals, at 510 nm (Fernandes et al., 2004). The shift depends on the proportion of these minerals in the spectrum and intensity of absorption bands. In the short-wavelength infrared (SWIR) region, kaolinite appeared as an absorption feature centered at 2,205 nm and gibbsite at 2,265 nm, separated by a maximum reflectance peak (R_max) at ± 2,230 nm (Madeira et al., 1995).

The maximum reflectance factor (RF) of the Latossolo Vermelho (LV) spectra ranged from 0.3 to 0.6, with an almost flat shape throughout the range from 350 to 2,500 nm (Figure 5). Highest values in the LV-I spectrum were recorded where the derived curve indicated a higher proportion of goethite (450 nm) than of hematite (580 nm) (Figure 5a), confirmed by the red color with a lighter tone, with hue 10R, value 3, and Chroma 6 (Table 5). The LV-II spectrum, with lower RF, had a more balanced goethite/hematite ratio in the derived spectrum (Figure 5a), represented by dark-red colors, with hue 10R, value 4, and Chroma 3. This difference in RF can be influenced by the different contents, proportions, and crystallinity degrees of Fe oxides, since the OM texture and content varied little in LVs (Table 5). In an analysis of iron oxides, Fernandes et al. (2004) found that in goethitic soils, the reduction of hematite content increased RF in the spectra. Even in soils with high proportion of goethite, Correa et al. (2008) reported red soil color. According to Schwertmann (1993), goethite is stable under various environmental conditions, tending to be dominant even in reddish soils, but always masked by the pigmenting effect of hematite.

Although OM varied little between the LV, its presence caused a reduction of RF within the range of 350 to 2,500 nm and mitigation of the absorption features from 350 to 1,350 nm (Viscarra Rossel et al., 2016) in the spectrum of LV-II (Figure 4) with 31.5 g kg^-1 of OM (Table 5).

The LV spectra presented absorption features of goethite (between 480 and 950 nm) and hematite (between 530 and 850 nm) (Figure 5). The first hematite feature is wider and the second narrower than the goethite feature. The presence of these Fe minerals in the LV can be associated to the parent material, containing ferruginous cement with different Fe₂O₃ levels (Freitas-Silva and Campos, 1998).

Features of kaolinite (1,400 and 2,205 nm) and gibbsite (2,265 nm) are clearly observable (Figure 5), of which the first has an adsorption “shoulder” or “step” on the left (Demattê et al., 2015b). The derived curve shows that these minerals have balanced proportions in both spectra, LV-I and LV-II (Figure 5). In soil genesis, kaolinite and gibbsite are related with hot and humid climate, high rainfall rates, and free drainage. Therefore, kaolinite is formed when desilication is moderate, while gibbsite occurs under conditions of total desilication. The latter is more abundant and common in the highly evolved Latossolos (Oxisols), developed by latolização (intense loss of basis, silica, and residual concentration of aluminum and iron oxides) (Kämpf et al., 2012). Features at 1,400 and 1,900 nm (Figure 5) are derived from vibrations of H-O-H groups of water adsorbed on the surface of minerals. At 1,400 nm, simultaneous vibrations of the O-H structural groups of clay occur (Stenberg et al., 2010). Thus, the morphology of the reflectance spectra of the studied LV (Figure 5) is similar to that of LV described by Bellinaso et al. (2010) and Demattê et al. (2014).

For the Latossolo Vermelho-Amarelo (LVA), the maximum reflectance factor (RF) was 0.40, with a flattened spectrum shape throughout the 350 to 2,500 nm range (Figure 6). The derived curve showed a higher goethite (450 nm) than hematite content, which had a lower peak amplitude (530 nm) (Figure 6a). Considering the effect of lower hematite contents in goethitic soils on increased reflectance (Fernandes et al., 2004), the smaller RF for LVA than for LV-I (RF 0.6) may be attributed to the higher OM content (Table 5), which caused an intensity decline throughout the entire spectrum (Viscarra Rossel et al., 2016). This effect demonstrates the complexity of the spectral analysis due to the simultaneous interaction of soil properties. Results of clay quantification based on reflectance spectroscopy can be substantially improved if OM is removed from the soil samples (Demattê et al., 2015b). However, from a practical standpoint, this procedure is infeasible.

Goethite was easily identified in the LVA spectrum, represented by an absorption feature of a range between 400 and 500 nm and an enlarged feature, centered at 950 nm (Figure 6). The derived spectrum shows a reflectance peak (450 nm) of large amplitude and an absorption feature at 480 nm (Figure 6a). Hematite was identified by a low-amplitude reflectance peak (580 nm) (Figure 6a). The dominant presence of goethite in the LVA is due to the greater stability of this mineral under deficient internal drainage conditions in the profiles of these soils, subjected to wetting and drying cycles during pedogenesis, which originated the concretionary or litoplíntico horizons in LVA (Motta et al., 2002; Barbosa et al., 2009).

The range of the kaolinite features (1,400 and 2,205 nm) was masked by OM, separated from the gibbsite feature (2,265 nm) by a R_max peak (2,230 nm). Hydroxyl groups were observed at 1,400 and 1,900 nm. The presence of gibbsite in the soil indicates a high weathering degree (Kämpf et al., 2012). Again, the spectral behavior of LVA described above is similar to that of LVA characterized by Bellinaso et al. (2010).

In the Plintossolo Pétrico (FF) soil spectra, the maximum reflectance factor (RF) ranged from 0.2 to 0.5 (Figure 7), influenced mainly by the OM effect. The intensity was highest in FF-I, with the lowest OM content (29 g kg^-1). For having the highest OM content (34 g kg^-1), FF-III had the lowest RF (Table 5). High OM contents in surface horizons caused a reduction in RF and smoothing of the features throughout the spectral range from 350 to 2,500 nm (Dalmolin et al., 2005), as evidenced in the spectrum of the FF-III (Figure 7). The high OM content of FFc was related to higher moisture levels in the profiles, because of the impediment to free water percolation through petroplinthites and the formation of concretionary horizons in these soils, as stated by Barbosa et al. (2009) and Campos et al. (2010).

The range of absorption features corresponding to goethite (400 to 500 nm) in the FF spectra (Figure 7) was softened by the OM effect in the FF-III spectrum. The derived curves (Figure 7a) indicated higher goethite contents by a reflectance peak of major amplitude (450 nm) than of hematite (580 nm). The broadest goethite feature was observed in the FF-I spectrum, centered at 950 nm, attenuated by the OM effect in the other two spectra. According to Moreira and Oliveira (2008), segregation of Fe oxides (hematite and goethite) in the form of nodules or concretions was observed in the FFc, according to Barbosa et al. (2009), resulting from a moister pedoenvironment due to the concretionary horizon. Consequently, hematite hydration occurs and its transformation into goethite, as described by Motta et al. (2002). The predominance of goethite associated with high OM contents induces dark brown soil colors, with hues ranging from 7.5YR to 10 YR, in agreement with studies of Fernandez and Schulze (1992) and Madeira et al. (1997) on soil color.

The FF also presented absorption features of kaolinite (1,400 and 2,205 nm), of hydroxyl groups (1,400 and 1,900 nm) and 2:1 clay minerals (1,900 nm), with low expression of gibbsite at 2,265 nm. Barbosa et al. (2009) verified the presence of illite (2:1 clay mineral) in the clay fraction of concretionary horizons of LVA from DF, along with kaolinite (1:1 clay mineral) and traces of gibbsite, demonstrating the lower degree of pedogenetic evolution than in LV. Despite the tropical environment, alternating conditions of heavy rainfalls and prolonged drought periods in the DF favor the preservation of clay minerals with greater activity in soils with internal drainage restrictions, as in FF (Anjos et al., 2007). In general, the spectral curves of FF are quite similar to those described by Demattê et al. (2012b).

For Neossolos Regolíticos (RR), the RF ranged from 0.35 to 0.50; RR-I was represented by a spectrum with greater intensity of RF, indicating silty-clay texture, and RR-II by a spectrum with less intensity, indicating a clayey texture (Figure 7).

The RR are soils with low evolution degrees and little action of pedogenetic processes, without significant modifications of the parent material (Santos et al., 2013). For this reason, the two RR spectra had lower Fe oxide contents than more evolved soils, such as Latossolos. In the derived spectrum (Figure 8a), a lower proportion of hematite (small amplitude peak at 580 nm) was observed in relation to goethite (450 nm), due to the initial pedogenetic process. The magnitude of this feature indicates a higher goethite content in the RR-II than the RR-I spectrum. The effect of the low hematite content in the presence of goethite induced yellowish colors in these soils, with hues of 10YR to 5YR (Table 5) (Fernandes et al., 2004).

The high OM content in the surface horizon of RR smoothed the spectral features of these soils within the range 350 to 1,350 nm (Figure 8). Kaolinite was detected by the spectral feature at 2,205 nm. Absence of gibbsite was identified, by the missing absorption feature at 2,265 nm, confirming the lower weathering degree of these soils (Figure 8b) (Santos et al., 2013). Adsorption features of hydroxyl groups (1,400 and 1,900 nm) and 2:1 clay minerals (1,900 nm) were observed (Figure 8), due to the occurrence of sparse illite in these soils, as described by Barbosa et al. (2009), Campos et al. (2010), and Lacerda and Barbosa (2012). Some of these 2:1 clay minerals have additional absorption features centered at approximately 2,450 nm (Grove et al., 1992). Similar features of these clay minerals at 2,355 and 2,448 nm were also observed by Demattê and Terra (2014). Moreover, Grove et al. (1992) showed that the presence of 2:1 minerals contribute to a greater intensity of the RF of soil. Therefore, both the silt content and presence of 2:1 clay minerals may have caused a difference in the albedo of the spectra (Figure 8). The presence of weatherable minerals in the silt and sand fraction probably influenced the soil spectral reflectance. In the RR-I spectrum, higher quartz contents in the sand fraction may have influenced the increase in albedo (Hunt and Salisbury, 1970) (Figure 8).

The Cambissolo Háplico (CX) spectrum showed a rectilinear shape along almost the entire spectral range (350 to 2,500 nm), with a maximum RF of 0.50 (Figure 9). This RF can be attributed to the high silt content (417 g kg^-1) (Table 4), with predominance of weatherable minerals in this fraction. The derived spectral curve (Figure 9a) showed low iron oxide contents in CX, with a higher proportion of goethite (450 nm) than hematite (580 nm). A centered feature was observed at 850 nm for hematite, while that of goethite, centered at 950 nm, was masked by OM. Cambissolos (Inceptisols) are young soils with high silt contents, indicating a low evolution degree (Santos et al., 2013). With the ongoing rock weathering, surface layers of CX are developed, where Fe oxides are still being formed (Demattê et al., 2003), for this reason, the adsorption features of goethite and hematite in the spectrum were little defined. The larger proportion of goethite associated with a high OM content gave this soil a brownish color (7.5YR 4/6).

Absorption features of hydroxyl groups (1,400 and 1,900 nm) and 2:1 clay minerals (1,900 and 2,450 nm) were identified. Kaolinite showed absorption feature at 2,205 nm, while gibbsite was absent (2,265 nm) (Figure 9b), confirming the lower evolution degree of CX, as mentioned by Demattê et al. (2014). These descriptions agree with Moller and Kitagawa (1982), who analyzed the mineralogy of Cambissolos from the south-eastern Brazilian Amazon and reported the occurrence of Fe oxides, kaolinite, micas, quartz, and 2:1 clay minerals, but the absence of gibbsite in all soil samples. Studies of Campos (2004) and Barbosa et al. (2009) showed that in soils derived from metasedimentary rocks of pelitic origin, as in studied watersheds, the micaceous clay minerals in the silt fraction, e.g., illite, accounted for the highest proportion.

The linearization of the spectral curve of CX is related to the high OM content in the surface horizon (Table 5), which provides additional absorption of incident energy and attenuates the effect of minerals, making the features less sharp. This agrees with observations of Formaggio et al. (1996) and Demattê et al. (2012a, 2014) on the masking effect of OM on soil mineralogy, causing a flattening of the curve in the 350-2,500 nm range. However, the flat linearity from 1,000 nm onward is an effect of a low quartz content. Therefore, the spectral curve as a whole is the result of multiple effects. This reinforces the importance of spectral morphology analysis by MIRS (Demattê et al., 2014) of the adsorption features and the spectra in general.

The spectroscopic features of CX (Figure 8) were similar to those of other Cambissolos described by several authors (Formaggio et al., 1996; Bellinaso et al., 2010; Demattê and Terra, 2014).

The maximum RFs of the Gleissolo Háplico (GX) spectrum were 0.5, with a positive trend (Figure 10). The derivative spectral curve (Figure 10a) shows goethite (450 nm) and hematite (580 nm). No typical adsorption feature of goethite was detected in the 400 - 500 nm range, as similarly observed in all other spectra (Figures 5 to 9). This mineral had an adsorption feature centered at 950 nm (Figure 10). The presence of goethite and hematite in this soil may have been caused by the beginning of the pedogenetic process of laterization, associated to GX when the water regime is modified, causing a shift from a reducing to an oxidizing system (Santos et al., 2013).

The spectrum pattern of GX can be a consequence of lower Fe oxide and OM contents in the A horizon, expressing the exposure of primary minerals in the silt fraction more intensely, with a quartz content of 414 g kg^-1 (Table 4). The variation in particle size and mineralogical composition of these soils is wide, due to the greater or lesser intensity of flooding events on floodplains and their respective sediment contributions (Reatto et al., 2000; Lacerda and Barbosa, 2012). Within the 1,400-2,205 nm range, the presence of kaolinite was observed with a typical feature of hydroxyl groups (1,400 and 1,900 nm), gibbsite (2,265 nm), and 2:1 clay minerals (1,900 nm) (Figure 9).

During flooding in the rainy periods, the floodplains of the studied watersheds as well as the Jardim River basin, receive sedimentation from local lithological material, since the DF is characterized by the presence of water springs, without major rivers. Thus, the material in hydromorphic soils deposited in the pedogenetic process in the DF consists of sediments originated from the lithology of the Paranoá, Canastra, and Bambuí groups (Freitas-Silva and Campos, 1998). This explains the presence of gibbsite and kaolinite in GX, since the sedimentary parent material has a meta-sedimentary lithology. These rocks, pre-weathered earlier in the genetic sedimentary process, were preserved since the metamorphic degree to which they were subjected was low (Freitas-Silva and Campos, 1998). The morphology of the GX spectrum coincides with descriptions of Demattê et al. (2012a).

CONCLUSIONS

Identification of soil-landscape relationships was fundamental for understanding pedogenetic processes along the toposequences, since the role of relief in the soil evolution of the studied area is extremely important.

The soils in the watersheds of the Estanislau Creek and Jardim Stream had similar physical and chemical properties, with little variations related to the parent materials.

Clustering based on soil surface reflectance was efficient in determining soil groups with similar properties.

The surface reflectance data were related to the soil surface and subsurface properties determined by traditional analyses of soil samples, since the two approaches formed similar groups

Simultaneous interactions of soil properties were assessed by analysis of spectral morphology by MIRS, while the second derivative of the KM function adequately quantified the mineralogy of the spectra.

ACKNOWLEDGMENTS

The authors are indebted to Capes for the Master's scholarship of the first author; to the Faculty of Agronomy and Veterinary Medicine of UnB; the Department of Soil Science of Esalq/USP; the research group GeoCis (http://esalqgeocis.wixsite.com/geocis/equipe); to Finatec for the financial support of this work; to Fapesp (grant number: 2014-22262-0) for providing the equipment used in the analysis of spectroscopy of reflectance (VIS-NIR-SWIR); and are particularly grateful to researcher and Professor Igo Fernando Lepsch for his critical review of the manuscript.

REFERENCES

• Adamchuk V, Allred B, Doolittle J, Grote K, Viscarra Rossel RA. Tools for proximal soil sensing. In: Ditzler C, Scheffe K, Monger HC, editors. Soil survey manual: soil science division staff. Washington, DC: United States Department of Agriculture; 2017. (Handbook, 18). p. 355-56
• Anjos LHC, Pereira MG, Pérez DV, Ramos DP. Caracterização e classificação de Plintossolos no município de Pinheiro-MA. Rev Bras Cienc Solo. 2007;31:1035-44. https://doi.org/10.1590/s0100-06832007000500020
  » https://doi.org/10.1590/s0100-06832007000500020
• Barbosa IO, Lacerda MPC, Bilich MR. Relações pedomorfogeológicas nas chapadas elevadas do Distrito Federal. Rev Bras Cienc Solo. 2009;33:1373-83. https://doi.org/10.1590/s0100-06832009000500029
  » https://doi.org/10.1590/s0100-06832009000500029
• Bellinaso H, Demattê JAM, Romeiro SA. Soil spectral library and its use in soil classification. Rev Bras Cienc Solo. 2010;34:861-70. https://doi.org/10.1590/s0100-06832010000300027
  » https://doi.org/10.1590/s0100-06832010000300027
• Ben-Dor E, Demattê JAM. Remote sensing of soil in the optical domains. In: Thenkabail OS, editor. Land resources monitoring, modeling, and mapping with remote sensing - Remote sensing handbook. Boca Ratom: CRC Press; 2015. v.2. p. 733-87.
• Ben-Dor E, Heller D, Chudnovsky A. A novel method of classifying soil profiles in the field using optical means. Soil Sci Soc Am J. 2008;72:1113-23. https://doi.org/10.2136/sssaj2006.0059
  » https://doi.org/10.2136/sssaj2006.0059
• Bouyoucos GJ. Estimation of the colloidal material in soils. Science. 1926;64:362. https://doi.org/10.1126/science.64.1658.362
  » https://doi.org/10.1126/science.64.1658.362
• Brown DJ, Shepherd KD, Walsh MG, Mays MD, Reinsch TG. Global soil characterization with VNIR diffuse reflectance spectroscopy. Geoderma. 2006;132:273-90. https://doi.org/10.1016/j.geoderma.2005.04.025
  » https://doi.org/10.1016/j.geoderma.2005.04.025
• Camargo LA, Marques Júnior J, Pereira GT. Spatial variability of physical properties of an Alfisol under different hillslope curvatures. Rev Bras Cienc Solo. 2010;34:617-30. https://doi.org/10.1590/s0100-06832010000300003
  » https://doi.org/10.1590/s0100-06832010000300003
• Campos JEG. Hidrogeologia do Distrito Federal: bases para a gestão dos recursos hídricos subterrâneos. Revista Brasileira de Geociências. 2004;34:41-8.
• Campos PM, Lacerda MPC, Silva CL, Sá MAC, Sousa DMG. Drenagem interna como fator de diferenciação de Latossolos do Distrito Federal. Pesq Agropec Bras. 2010;45:306-14. https://doi.org/10.1590/s0100-204x2010000300011
  » https://doi.org/10.1590/s0100-204x2010000300011
• Cezar E, Nanni MR, Chicati ML, Oliveira RB, Demattê JAM. Discriminação entre solos formados em região transicional por meio de sua resposta espectral. Biosci J. 2013;29:644-54.
• Chagas CS, Fernandes Filho EI, Bhering SB. Relação entre atributos do terreno, material de origem e solos em uma área no noroeste do estado do Rio de Janeiro. Soc Nat. 2013;25:147-62. https://doi.org/10.1590/s1982-45132013000100012
  » https://doi.org/10.1590/s1982-45132013000100012
• Chaves AAA. Distribuição dos solos associados ao cultivo de cafés de qualidade especial na Serra da Mantiqueira Mineira [tese]. Brasília, DF: Universidade de Brasília; 2016.
• Chicati ML, Nanni MR, Cézar E, Oliveira RB, Chicati MS. Diffuse reflectance spectroscopy to estimate soil properties of Brazilian wetlands. Int J Environ Agric Res. 2016;2:62-6.
• Claessen MEC, organizador. Manual de métodos de análise de solo. 2. ed. Rio de Janeiro: Embrapa Solos; 1997.
• Correa MM, Ker JC, Barrón V, Fontes MPF, Torrent J, Curi N. Caracterização de óxidos de ferro de solos do ambiente tabuleiros costeiros. Rev Bras Cienc Solo. 2008;32:1017-31. https://doi.org/10.1590/s0100-06832008000300011
  » https://doi.org/10.1590/s0100-06832008000300011
• Cunha JE, Nóbrega MT, Castro SS. Infiltração da água no solo no sistema pedológico Campus do Arenito, cidade gaúcha, noroeste do estado do Paraná. Rev Bras Cienc Solo. 2008;32:1837-48. https://doi.org/10.1590/s0100-06832008000500005
  » https://doi.org/10.1590/s0100-06832008000500005
• Dalmolin RSD, Gonçalves CN, Klamt E, Dick DP. Relação entre os constituintes do solo e seu comportamento espectral. Cienc Rural. 2005;35:481-9. https://doi.org/10.1590/s0103-84782005000200042
  » https://doi.org/10.1590/s0103-84782005000200042
• Demattê JAM, Araújo SR, Fiorio PR, Fongaro CT, Nanni MR. Espectroscopia VIS-NIR-SWIR na avaliação de solos ao longo de uma topossequência em Piracicaba (SP). Rev Cienc Agron. 2015b;46:679-88. https://doi.org/10.5935/1806-6690.20150054
  » https://doi.org/10.5935/1806-6690.20150054
• Demattê JAM, Bellinaso H, Romero DJ, Fongaro CT. Morphological interpretation of reflectance spectrum (MIRS) using libraries looking towards soil classification. Sci Agric. 2014;71:509-20. https://doi.org/10.1590/0103-9016-2013-0365
  » https://doi.org/10.1590/0103-9016-2013-0365
• Demattê JAM, Campos RC, Alves MC, Fiorio PR, Nanni MR. Visible–NIR reflectance: a new approach on soil evaluation. Geoderma. 2004;121:95-112. https://doi.org/10.1016/j.geoderma.2003.09.012
  » https://doi.org/10.1016/j.geoderma.2003.09.012
• Demattê JAM, Fiorio PR, Clemente C, Nanni MR. Sensoriamento remoto na avaliação da evolução intempérica de solos desenvolvidos sobre granodiorito e gnaisse porfirítico. Magistra. 2003;15:33-43.
• Demattê JAM, Galdos MV, Guimarães RV, Genú AM, Nanni MR, Zullo Junior J. Quantification of tropical soil properties from ETM+/LANDSAT-7 data. Int J Remote Sens. 2007;28:3813-29. https://doi.org/10.1080/01431160601121469
  » https://doi.org/10.1080/01431160601121469
• Demattê JAM, Huete AR, Ferreira Junior LG, Nanni MR, Alves MC, Fiorio PR. Methodology for bare soil detection and discrimination by landsat Tm image. Open Remote Sens J. 2009;2:24-35. https://doi.org/10.2174/1875413901002010024
  » https://doi.org/10.2174/1875413901002010024
• Demattê JAM, Morgan CLS, Chabrillat S, Rizzo R, Franceschini MHD, Terra FS, Vasques GM, Wetterlind J. Spectral sensing from ground to space in soil science: state of the art, applications, potential, and perspectives. In: Thenkabail OS, editor. Land resources monitoring, modeling, and mapping with remote sensing - Remote sensing handbook. Boca Ratom: CRC Press; 2015a. v.2. p. 661-732.
• Demattê JAM, Terra FS. Spectral pedology: a new perspective on evaluation of soils along pedogenetic alterations. Geoderma. 2014;217-218:190-200. https://doi.org/10.1016/j.geoderma.2013.11.012
  » https://doi.org/10.1016/j.geoderma.2013.11.012
• Demattê JAM, Terra FS, Quartaroli CF. Spectral behavior of some modal soil profiles from São Paulo State, Brazil. Bragantia. 2012a;71:413-23. https://doi.org/10.1590/s0006-87052012005000038
  » https://doi.org/10.1590/s0006-87052012005000038
• Demattê JAM, Vasques GM, Corrêa EA, Arruda GP. Fotopedologia, espectroscopia e sistema de informação geográfica na caracterização de solos desenvolvidos do Grupo Barreiras no Amapá. Bragantia. 2012b;71:438-46.
• Demattê JAM, Ramirez-Lopez L, Rizzo R, Nanni MR, Fiorio PR, Fongaro CT, Medeiros Neto LG, Safanelli JL, Barros PPS. Remote sensing from ground to space platforms associated with terrain properties as a hybrid strategy on the development of a pedological map. Remote Sens. 2016;8:826. https://doi.org/10.3390/rs8100826
  » https://doi.org/10.3390/rs8100826
• Fernandes RBA, Barrón V, Torrent J, Fontes MPF. Quantificação de óxidos de ferro de Latossolos brasileiros por espectroscopia de refletância difusa. Rev Bras Cienc Solo. 2004;28:245-57. https://doi.org/10.1590/s0100-06832004000200003
  » https://doi.org/10.1590/s0100-06832004000200003
• Fernandez RN, Schulze DG. Munsell colors of soils simulated by mixtures of goethite and hematite with kaolinite. Z Pflanzenernaehr Bodenk. 1992;155:473-8. https://doi.org/10.1002/jpln.19921550520
  » https://doi.org/10.1002/jpln.19921550520
• Foley JA, Ramankutty N, Brauman KA, Cassidy ES, Gerber JS, Johnston M, Mueller ND, O'Connell C, Ray DK, West PC, Balzer C, Bennett EM, Carpenter SR, Hill J, Monfreda C, Polasky S, Rockström J, Sheehan J, Siebert S, Tilman D, Zaks DPM. Solutions for a cultivated planet. Nature. 2011;478:337-42. https://doi.org/10.1038/nature10452
  » https://doi.org/10.1038/nature10452
• Formaggio AR, Epiphanio JCN, Valeriano MM, Oliveira JB. Comportamento espectral (450-2.450 nm) de solos tropicais de São Paulo. Rev Bras Cienc Solo. 1996;20:467-74.
• Freitas-Silva FH, Campos JEG. Geologia do Distrito Federal. In: Campos JEG, Freitas-Silva FH, organizadores. Inventário hidrogeológico e dos recursos hídricos superficiais do Distrito Federal. Brasília, DF: IEMA/SEMATEC/UnB; 1998. v.1.
• Gomes JBV, Curi N, Schulze DG, Marques JJGSM, Ker JC, Motta PEF. Mineralogia, morfologia e análise microscópica de solos do bioma cerrado. Rev Bras Cienc Solo. 2004;28:679-94. https://doi.org/10.1590/s0100-06832004000400010
  » https://doi.org/10.1590/s0100-06832004000400010
• Grove CI, Hook SJ, Paylor II ED. Laboratory reflectance spectra of 160 minerals, 0.4 to 2.5 micrometers. Pasadena: Jet Propulsion Laboratory; 1992.
• Hunt GR, Salisbury JW. Visible and near infrared spectra of minerals and rocks: I. Silicate minerals. Mod Geol. 1970;1:283-300.
• Instituto Brasileiro de Geografia e Estatística - IBGE. Mapa Índice Digital – MID: mapeamento geral do Brasil [CD-ROM]. 4. ed. Rio de Janeiro: IBGE/DSG; 2011.
• Ippoliti GAR, Costa LM, Schaefer CEGR, Fernandes Filho EI, Gaggero MR. Análise digital do terreno: ferramenta na identificação de pedoformas em microbacia na região de “Mar de Morros” (MG). Rev Bras Cienc Solo. 2005;29:269-76. https://doi.org/10.1590/s0100-06832005000200012
  » https://doi.org/10.1590/s0100-06832005000200012
• Kämpf N, Curi N. Formação e evolução do solo (pedogênese). In: Ker JC, Curi N, Schaefer CEGR, Vidal-Torrado P, editores. Pedologia: fundamentos. Viçosa, MG: Sociedade Brasileira de Ciência do Solo; 2012. p. 207-302.
• Kämpf N, Marques JJ, Curi N. Mineralogia de solos brasileiros. In: Ker JC, Curi N, Schaefer CEGR, Vidal-Torrado P, editores. Pedologia: fundamentos. Viçosa, MG: Sociedade Brasileira de Ciência do Solo; 2012. p. 81-146.
• Lacerda MPC, Barbosa IO. Relações pedomorfogeológicas e distribuição de pedoformas na estação ecológica de águas emendadas, Distrito Federal. Rev Bras Cienc Solo. 2012;36:709-22. https://doi.org/10.1590/s0100-06832012000300003
  » https://doi.org/10.1590/s0100-06832012000300003
• Lacerda MPC, Queménéur JJG, Andrade H, Alves HMR, Vieira TGC. Estudo da relação pedomorfogeológica na distribuição de solos com horizontes B textural e B nítico na paisagem de Lavras (MG). Rev Bras Cienc Solo. 2008;32:271-84. https://doi.org/10.1590/s0100-06832008000100026
  » https://doi.org/10.1590/s0100-06832008000100026
• Madeira J, Bedidi A, Cervelle B, Pouget M, Flay N. Visible spectrometric indices of hematite (Hm) and goethite (Gt) content in lateritic soils: the application of a Thematic Mapper (TM) image for soil-mapping in Brasilia, Brazil. Int J Remote Sens. 1997;18:2835-52. https://doi.org/10.1080/014311697217369
  » https://doi.org/10.1080/014311697217369
• Madeira J, Bédidi A, Pouget M, Cervelle B, Flay N. Spectral (MIR) determination of kaolinite and gibbsite contents in lateritic soils. CR Acad Sci II A. 1995;321:119-28.
• Mohamed ES, Saleh AM, Belal AB, Gad A-A. Application of near-infrared reflectance for quantitative assessment of soil properties. Egypt J Remote Sensing Space Sci. In press 2017. https://doi.org/10.1016/j.ejrs.2017.02.001
  » https://doi.org/10.1016/j.ejrs.2017.02.001
• Moller MRF, Kitagawa Y. Mineralogia de argilas em Cambissolos do sudoeste da Amazônia brasileira. Belém: Embrapa-CPATU; 1982. (Boletim de pesquisa, 34).
• Morales N, Assine ML. Chapada do Araripe: a highland oasis incrusted into the semi-arid region of northeastern Brazil. In: Vieira BC, Salgado AAR, Santos LJC, editores. Landscapes and landforms of Brazil. Netherlands: Springer; 2015. p. 231-42.
• Moreira HL, Oliveira VA. Evolução e gênese de um Plintossolo Pétrico concrecionário êutrico argissólico no município de Ouro Verde de Goiás. Rev Bras Cienc Solo. 2008;32:1683-90. https://doi.org/10.1590/s0100-06832008000400033
  » https://doi.org/10.1590/s0100-06832008000400033
• Motta PEF, Carvalho Filho A, Ker JC, Pereira NR, Carvalho Junior W, Blancaneaux P. Relações solo-superfície geomórfica e evolução da paisagem em uma área do Planalto Central Brasileiro. Pesq Agropec Bras. 2002;37:869-78. https://doi.org/10.1590/s0100-204x2002000600017
  » https://doi.org/10.1590/s0100-204x2002000600017
• Mouazen AM, Kuang B, Baerdemaeker J, Ramon H. Comparison among principal component, partial least squares and back propagation neural network analyses for accuracy of measurement of selected soil properties with visible and near infrared spectroscopy. Geoderma. 2010;158:23-31. https://doi.org/10.1016/j.geoderma.2010.03.001
  » https://doi.org/10.1016/j.geoderma.2010.03.001
• Nanni MR, Demattê JAM, Chicati ML, Fiorio PR, Cézar E, Oliveira RB. Soil surface spectral data from Landsat imagery for soil class discrimination. Acta Sci-Agron. 2012;34:103-12. https://doi.org/10.4025/actasciagron.v34i1.12204
  » https://doi.org/10.4025/actasciagron.v34i1.12204
• Nanni MR, Demattê JAM, Chicati ML, Oliveira RB, Cézar E. Spectroradiometric data as support to soil classification. Int Res J Agric Sci Soil Sci. 2011;1:109-17.
• Nocita M, Stevens A, Wesemael B, Aitkenhead M, Bachmann M, Barthès B, Ben Dor E, Brown DJ, Clairotte M, Csorba A, Dardenne P, Demattê JAM, Genot V, Guerrero C, Knadel M, Montanarella L, Noon C, Ramirez-Lopez L, Robertson J, Sakai H, Soriano-Disla JM, Shepherd KD, Stenberg B, Towett EK, Vargas R, Wetterlind J. Soil spectroscopy: an alternative to wet chemistry for soil monitoring. In: Sparks DL, editor. Advances in Agronomy. San Diego: Elsevier; 2015. p. 139-59. https://doi.org/10.1016/bs.agron.2015.02.002
  » https://doi.org/10.1016/bs.agron.2015.02.002
• Oliveira DP, Ferreira TO, Romero RE, Farias PRS, Costa MCG. Microrrelevo e a distribuição de frações granulométricas em Cambissolos de origem calcária. Rev Cienc Agron. 2013;44:676-84. https://doi.org/10.1590/s1806-66902013000400003
  » https://doi.org/10.1590/s1806-66902013000400003
• Organisation for Economic Co-operation and Development / Food and Agriculture Organization of the United Nations - OECD/FAO. OECD-FAO Agricultural Outlook 2016-2025. Paris: OECD Publishing; 2016.
• Pedron FA, Fink JR, Rodrigues MF, Azevedo AC. Condutividade e retenção de água em Neossolos e saprolitos derivados de arenito. Rev Bras Cienc Solo. 2011;35:1253-62. https://doi.org/10.1590/s0100-06832011000400018
  » https://doi.org/10.1590/s0100-06832011000400018
• Pennock DJ, Veldkamp A. Advances in landscape-scale soil research. Geoderma. 2006;133:1-5. https://doi.org/10.1016/j.geoderma.2006.03.032
  » https://doi.org/10.1016/j.geoderma.2006.03.032
• Pereira TTC, Ker JC, Schaefer CEGR, Barros NF, Neves JCL, Almeida CC. Gênese de Latossolos e Cambissolos desenvolvidos de rochas pelíticas do grupo Bambuí - Minas Gerais. Rev Bras Cienc Solo. 2010;34:1283-95. https://doi.org/10.1590/s0100-06832010000400026
  » https://doi.org/10.1590/s0100-06832010000400026
• Pinto MN. Caracterização geomorfológica do Distrito Federal. In: Pinto MN, organizador. Cerrado: caracterização, ocupação e perspectivas. 2. ed. Brasília, DF: Editora Universidade de Brasília; 1994. p. 285-344.
• Poppiel RR. Mapeamento pedológico por meio de série histórica Landsat-5 TM e biblioteca espectral na Bacia do Rio Jardim (DF) [dissertação]. Brasília, DF: Universidade de Brasília; 2016.
• Reatto A, Correia JR, Spera ST, Chagas CS, Martins ES, Andahur JP, Godoy MJS, Assad MLCL. Levantamento semidetalhado dos solos da bacia do rio Jardim-DF, escala 1:50.000. Planaltina: Embrapa Cerrados; 2000. (Boletim de pesquisa, 18).
• Ribeiro MR, Oliveira LB, Araújo Filho JC. Caracterização morfológica do solo. In: Ker JC, Curi N, Schaefer CEGR, Vidal-Torrado P, editores. Pedologia: fundamentos. Viçosa, MG: Sociedade Brasileira de Ciência do Solo; 2012. p. 47-80.
• Santos HG, Jacomine PKT, Anjos LHC, Oliveira VA, Oliveira JB, Coelho MR, Lumbreras JF, Cunha TJF. Sistema brasileiro de classificação de solos. 3. ed. rev. ampl. Rio de Janeiro: Embrapa Solos; 2013.
• Santos RD, Santos HG, Ker JC, Anjos LHC, Shimizu SH. Manual de descrição e coleta de solo no campo. 7. ed rev ampl. Viçosa, MG: Sociedade Brasileira de Ciência do Solo; 2015.
• Schulze DG, Nagel JL, Scoyoc GE, Henderson TL, Baumgardner MF, Stott DE. Significance of organic matter in determining soil colors. In: Bigham JM, Ciolkosz EJ, editors. Soil Color. Madison: Soil Science Society of America; 1993. p. 71-90.
• Schwertmann U. Relations between iron oxides, soil color, and soil formation. In: Bigham JM, Ciolkosz EJ, editors. Soil Color. Madison: Soil Science Society of America; 1993. p. 51-70.
• Sellitto VM, Fernandes RBA, Barrón V, Colombo C. Comparing two different spectroscopic techniques for the characterization of soil iron oxides: diffuse versus bi-directional reflectance. Geoderma. 2009;149:2-9. https://doi.org/10.1016/j.geoderma.2008.11.020
  » https://doi.org/10.1016/j.geoderma.2008.11.020
• Secretaria de Estado de Meio Ambiente e Recursos Hídricos - Semarh. Zoneamento ecológicoeconômico do Distrito Federal (ZEE-DF): Subproduto 3.1. Brasília, DF: Semarh; 2012. [Last accessed: October 20, 2016]. Available at: www.zee.df.gov.br/historico/arquivos/doc_download/23-volume-ii-meio-fisico.html
  » www.zee.df.gov.br/historico/arquivos/doc_download/23-volume-ii-meio-fisico.html
• Senthilkumar S, Kravchenko AN, Robertson GP. Topography influences management system effects on total soil carbon and nitrogen. Soil Sci Soc Am J. 2009;73:2059-67. https://doi.org/10.2136/sssaj2008.0392
  » https://doi.org/10.2136/sssaj2008.0392
• Sherman DM, Waite TD. Electronic spectra of Fe³⁺ oxides and oxide hydroxides in the near IR to near UV. Am Mineral. 1985;70:1262-9.
• Sneath PHA, Sokal RR. Numerical taxonomy: the principles and practice of numerical classification. San Francisco: WH Freeman and Company; 1973.
• Soil Survey Staff. Keys to soil taxonomy. 12th ed. Washington, DC: United States Department of Agriculture, Natural Resources Conservation Service; 2014.
• Soriano-Disla JM, Janik LJ, Viscarra Rossel RA, Macdonald LM, McLaughlin MJ. The performance of visible, near-, and mid-infrared reflectance spectroscopy for prediction of soil physical, chemical, and biological properties. Appl Spectrosc Rev. 2014;49:139-86. https://doi.org/10.1080/05704928.2013.811081
  » https://doi.org/10.1080/05704928.2013.811081
• Stenberg B, Viscarra Rossel RA, Mouazen AM, Wetterlind J. Visible and near infrared spectroscopy in soil science. Adv Agron. 2010;107:163-215. https://doi.org/10.1016/s0065-2113(10)07005-7
  » https://doi.org/10.1016/s0065-2113(10)07005-7
• Torrent J, Barrón V. Diffuse reflectance spectroscopy. In: Ulery AL, Drees LR, editors. Methods of soil analysis. Mineralogical methods. Madison: Soil Sci Soc Am; 2008. Pt 5. p. 367-85.
• Vasques GM, Demattê JAM, Viscarra Rossel RA, Ramírez-López L, Terra FS. Soil classification using visible/near-infrared diffuse reflectance spectra from multiple depths. Geoderma. 2014;223-225:73-8. https://doi.org/10.1016/j.geoderma.2014.01.019
  » https://doi.org/10.1016/j.geoderma.2014.01.019
• Vettori L. Métodos de análise de solo. Rio de Janeiro: Equipe de Pedologia e Fertilidade do Solo; 1969. (Boletim Técnico, 7).
• Vidal-Torrado P, Lepsch IF, Castro SS. Conceitos e aplicação das relações pedologia-geomorfologia em regiões tropicais úmidas. In: Vidal-Torrado P, Lepsch IF, Castro SS, organizadores. Tópicos em Ciência do Solo. Viçosa: Sociedade Brasileira de Ciência do Solo; 2005. v.4. p. 145-92.
• Viscarra Rossel RA, Behrens T, Ben-Dor E, Brown DJ, Demattê JAM, Shepherd KD, Shi Z, Stenberg B, Stevens A, Adamchuk V, Aïchi H, Barthès BG, Bartholomeus HM, Bayer AD, Bernoux M, Böttcher K, Brodský L, Du CW, Chappell A, Fouad Y, Genot V, Gomez C, Grunwald S, Gubler A, Guerrero C, Hedley CB, Knadel M, Morrás HJM, Nocita M, Ramirez-Lopez L, Roudier P, Campos EMR, Sanborn P, Sellitto VM, Sudduth KA, Rawlins BG, Walter C, Winowiecki LA, Hong SY, Ji W. A global spectral library to characterize the world's soil. Earth-Sci Rev. 2016;155:198-230. https://doi.org/10.1016/j.earscirev.2016.01.012
  » https://doi.org/10.1016/j.earscirev.2016.01.012
• Viscarra Rossel RA, Cattle SR, Ortega A, Fouad Y. In situ measurements of soil colour, mineral composition and clay content by vis-NIR spectroscopy. Geoderma. 2009;150:253-66. https://doi.org/10.1016/j.geoderma.2009.01.025
  » https://doi.org/10.1016/j.geoderma.2009.01.025
• Viscarra Rossel RA, Webster R. Discrimination of Australian soil horizons and classes from their visible-near infrared spectra. Eur J Soil Sci. 2011;62:637-47. https://doi.org/10.1111/j.1365-2389.2011.01356.x
  » https://doi.org/10.1111/j.1365-2389.2011.01356.x
• Walkley A, Black IA. An examination of the Degtjareff method for determining soil organic matter and a proposed modification of the chromic acid titration method. Soil Science. 1934;37:29-37. http://dx.doi.org/10.1097/00010694-193401000-00003
  » http://dx.doi.org/10.1097/00010694-193401000-00003
• Xie H, Zhao J, Wang Q, Sui Y, Wang J, Yang X, Zhang X, Liang C. Soil type recognition as improved by genetic algorithm-based variable selection using near infrared spectroscopy and partial least squares discriminant analysis. Sci Rep. 2015;5:10930. https://doi.org/10.1038/srep10930
  » https://doi.org/10.1038/srep10930
• Zeng R, Rossiter DG, Yang F, Li D-C, Zhao Y-G, Zhang G-L. How accurately can soil classes be allocated based on spectrally predicted physio-chemical properties? Geoderma. 2017;303:78-84. https://doi.org/10.1016/j.geoderma.2017.05.011
  » https://doi.org/10.1016/j.geoderma.2017.05.011
• Zeng R, Zhang G-L, Li D-C, Rossiter DG, Zhao Y-G. How well can VNIR spectroscopy distinguish soil classes? Biosyst Eng. 2016;152:117-25. https://doi.org/10.1016/j.biosystemseng.2016.04.019
  » https://doi.org/10.1016/j.biosystemseng.2016.04.019

Publication Dates

• Publication in this collection
  2018

History

• Received
  25 Nov 2016
• Accepted
  11 Aug 2017