Next Article in Journal
Ciliated Protist Communities in Soil: Contrasting Patterns in Natural Sites and Arable Lands across Italy
Previous Article in Journal
Effectiveness of Bacillus paramycoides for Improving Zinc Nutrition of Rice Irrigated with Alkali Water
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fast High-Resolution pKa Spectrotitrimetry for Quantification of Surface Functional Groups of Retisols

by
Natal’ya V. Matveeva
1,
Andrei V. Garmash
2,
Mikhail A. Shishkin
1,
Alexey A. Dymov
3,
Olga B. Rogova
2,
Dmitry S. Volkov
2 and
Mikhail A. Proskurnin
2,*
1
V.V. Dokuchaev Soil Science Institute, 7 Bld. 2 Pyzhevskiy Per., Moscow 119017, Russia
2
Chemistry Department, M.V. Lomonosov Moscow State University, Moscow 119991, Russia
3
Institute of Biology of Komi Scientific Centre, Ural Branch of Russian Academy of Sciences, Syktyvkar 167982, Russia
*
Author to whom correspondence should be addressed.
Soil Syst. 2024, 8(2), 63; https://doi.org/10.3390/soilsystems8020063
Submission received: 2 April 2024 / Revised: 26 May 2024 / Accepted: 6 June 2024 / Published: 12 June 2024

Abstract

:
Potentiometric titration in a fast and simultaneously high-resolution modality was proposed for the identification and quantification of protolytic groups of variable strength at the surface of primary soil particles. The method is implemented by titrimetric data processing as multicomponent spectra (pKa spectrotitrimetry). Due to the high resolution, the error of acidity-constant assessment (3–5%) is lower, compared to existing approaches; due to the fast titration, the effect of soil hydrolysis is minimized. The soil profiles for acidic Retisols (podzolic soils)—under a spruce crown and in the intercrown space—were studied. These soils, which have similar bulk properties and genesis but developed under different plant covers, were distinguished by pKa spectral features at 4–5; 5.5–6.5; 6.5–8.5; 7.5–8.5; and 9–10, as well as total group concentrations. Differences in acidic and basic-group distribution (carboxyl groups, amorphous aluminosilicates, carbonate species, amino groups, soluble (poly)phenolic compounds, phospholipids) and Al and Fe complex compounds within the same soil profiles and between two Retisols were found and quantified. The acidity constants and group concentrations found by pKa spectrotitrimetry were compared with conventional soil-composition indicators (total organic carbon, oxalate-soluble Fe and Al, and phosphorus), using principal component analysis. The main correlations are between the concentrations of oxalate-soluble Al and groups with pKa values of 5.0–6.5 and 8.5; oxalate-soluble Fe and pKa values of 9.0–10.0; and P2O5 and pKa values of 4.0–6.0 and 6.5–8.5. The method provides a set of major acidity values without a priori information on a soil sample and can be used for screening and identifying similar soils.

1. Introduction

Acidic soils are widespread, including almost all soils in humid areas in northern Europe and Asia. In Russia, these soils occupy more than 50% of the land [1] and are mainly represented by podzols and podzolic soils. These soils are formed in a broad range of climatic conditions, on sandy and loam soil-forming rocks of different mineralogical and granulometric composition. These soils are of low fertility, mainly due to their characteristic low pH values, low organic matter content, geographical location, predominantly in humid (boreal and subboreal) regions, and variable composition, associated with the diversity of soil-forming glacial-fluvio-sedimentary rocks, topography, and vegetation cover. In many regions, they are used for growing industrial and grain crops and forage grasses.
Acidic soil surface properties are connected with compounds of various classes, which are governed by the acid–base properties of their adsorbates and their charges, i.e., surface functional groups of primary soil particles [2,3]. In turn, the acid–base state determines the cation exchange capacity of soils, the mobility of chemical elements [4,5], the composition of exchangeable cations, and enzymatic activity, affecting redox processes, complex formation, precipitation, and dissolution. Soil surfaces also contain organic matter, which determines the soil structure, its composition and water resistivity [6,7,8]. The physical and chemical properties of the surface of soil particles are needed for studies of acidic soils, in addition to understanding changes in the structure and buffering properties of basic clay minerals under the influence of acid or alkaline deposition [9,10,11]. Moreover, such studies are used to improve the doses in soil liming and the methods for predicting their impact on fertility [12].
At present, at least a dozen of different functional groups is attributed to soil organic protolytic groups: carboxyl; phenolic and alcoholic hydroxyl; carbonyl; quinone; methoxy; ester; amino, amide, and imido; and sulfonic, thiol, and disulfide groups [2,13,14,15]. Inorganic fragments also serve as a source of formation of surface charge and pH, which usually occurs during the splitting of mineral particles and deformation of the crystal lattices of clay minerals and amorphous minerals [16]. However, despite multiple studies devoted to various aspects of the formation of soil acidity [17,18], the contribution of acid groups and vegetation cover remains insufficiently studied.
Potentiometric acid–base titration of soil suspensions and extracts is a common technique [19,20,21,22,23,24]. However, the analysis of titration curves is not straightforward in this case, as soil is a multicomponent and polydisperse system, and its components are mostly multifunctional. For both classification and quantification, information on individual components of soil suspensions and extracts is required [25,26]. Their composition is needed to assess the state of the soil cover, as its composition depends on soil and tillage system types [27]. However, the ionization and dissociation constants of these constituents and functional groups are close to each other [28], and most titration approaches cannot identify individual soil components due to the low resolution of titration curves and imperfect algorithms [29]. Thus, titration provides integral parameters only, such as total acid–base buffer capacity [14,30,31,32,33,34] or buffering intensities [28,35,36,37,38].
Moreover, the limited capabilities of acid–base titration to determine the type and quantity of acid or basic groups at the surface of primary soil particles are caused by the time between the titrant addition and equilibrium pH measurement, which should be short enough. Although acid–base equilibria are established very fast, soil hydrolysis may lead to its deep dispersion and transformation, which distorts the titration curve [39,40]. Thus, existing acid–base titration approaches are either low-resolution or slow, and therefore are not reliable for finding the detailed composition of soil surface groups.
We propose a high-resolution and at the same time fast potentiometric titrimetry technique, which maximally details the titration curve, while shortening the analysis time, thus decreasing the amounts of hydrolysis products. This technique is implemented in the form of pKa spectrometry (spectrotitrimetry), i.e., using the processing of the resulting continuous titration curves as arbitrary multicomponent acidity spectra. The technique was previously used for the concentrations and pKa values of ionogenic groups of mono- and polyfunctional polyelectrolytes: polyacrylic and polymethacrylic acids, gelatin, and humic substances [41,42,43,44]. However, so far, it has not been implemented to whole soil samples due to the lack of instruments capable of acquiring high-resolution titration curves at the required speed.
In this study, this technique is used to determine the protolytic groups at the surface of particles of two podzolic soils [45] (Russian soil classification) or Retisols (WRB classification). These soils were selected due to their wide distribution in the taiga region of northern Europe (Finland and Sweden), Canada, and the European part of Russia. The aim of this study was to apply the technique of fast high-resolution pKa potentiometric titrimetry for identifying and quantifying the organic acidity functional groups in acidic soils.

2. Materials and Methods

2.1. Area Description and Soil Sampling

The studies were carried out in Ust-Kulomsky District of Komi Republic (middle taiga), on Zhezhim Parma hill, in a blueberry native old spruce forest (Figure S1, Supplementary Information). The mean annual air temperature is from +0.2 to +1.3 °C. The amount of annual precipitation is from 670 to 790 mm [46]. The profile of albic Glossic Retisol (Siltic, Cutanic) (Retisol I, RI) was studied at intercrown (open) space on a distance of approximately 4 m from the crown of the spruce. The soil profile of albic Ferric Retisol (Siltic, Cutanic) (Retisol II, RII) was opened under the crown (1 m from the stand) of a 260-year-old spruce (Picea obovate Ledeb.). Soil classification is given according to WRB (2022). The distance between pits is about 10 m. Samples were selected for chemical analysis from each genetic horizon. Samples were collected in 2014. For more information, photographs of the spruce forest, where the profiles were opened and studied, and the soil profiles are presented in Figure 1 and Table 1. Air-dried samples were ground and sieved through a 1 mm sieve. More detailed information on the characteristics of the research objects was previously published [47].

2.2. General Soil Analysis

Samples were air-dried, and living roots and all particles with a diameter over 1 mm were removed by dry sieving. Organic samples (Oi, Oe, Oa) were ground with a sample grinder WCG75 (Pro Prep, Torrington, CT, USA). Chemical analysis of the soils was performed in the Ecoanalit laboratory and the Soil Science Department of the Institute of Biology, Federal Research Center, Komi Research Center of the Ural Branch of the Russian Academy of Sciences (IB FRC Komi SC UrB RAS). The total contents of carbon (Ctot) and nitrogen (Ntot) were determined by dry combustion on an EA-1110 CHNS-O analyzer (Carlo Erba, Cornareda, Italy). Exchangeable cations were determined by displacing with 1 M KC1, followed by atomic-absorption determination with a Hitachi 180-60 instrument (Tokyo, Japan). The pH values were determined on an Anion-4100 ionomer (Moscow, Russia) at a soil-to-water ratio of 1 (w):2.5 (v) for mineral horizons and 1 (w):25 (v) for organic horizons, respectively [48]. The particle size distribution was determined by the pipetting method, according to Kachinsky, with dispersion and boiling in the presence of NaOH, according to the Russian State Standard (GOST 12536-2014, Soils. Methods of laboratory granulometric (grain-size) and microaggregate distribution). Oxalate-soluble compounds of iron and aluminum were extracted by the Tamm method [49], and dithionite-soluble iron compounds by the method of Mehra and Jackson [50]. The total elemental analysis was determined with a ReSpect energy-dispersive X-ray fluorescence spectrometer (LLC Technoanalytdevice, Moscow, Russia). Phosphorus was determined by inductively coupled plasma atomic emission spectroscopy (ICP SpectroCiros CCD, Spectro Al, Kleve, Germany).
All weighing was carried out on an Ohaus Adventurer Pro AV313C analytical balance (Parsippany, NJ, USA). Glass and plastic chemical workware were used. The autosampler of the AT-710M titrator had 7 positions for beakers: position 1 was used for washing and cleaning the working part of the titrator after each titration; the remaining six positions were used for test samples (volume of vessels, 150 mL, selected for the optimum soil weighed portions of 3 g).

2.3. Reagents and Solvents

Distilled water was used to prepare the reagents and wash the workware. Working solutions of HCl and NaOH were prepared from ampouled reference samples every week; the solutions were stored in closed polypropylene containers connected to the titrator instrument. The titration reagent was collected automatically by the instrument software command. Preparation of saline solutions of acids and titrants was carried out similarly: 0.1M NaOH was prepared from ampouled reference samples and dry NaCl (5.85 g NaCl for a volume of 1 L).
Solutions of 1 M KC1 (puriss.) and 1 M CH3COONa (puriss.) were prepared from weighed portions of 75.5 and 82.0 g, respectively, in 1.000 L of water. All solutions were prepared in volumetric flasks with a volume of 1 L. The solution of 0.1 M H3BO3 was prepared from 6.18 g of a dry sample: a 0.1 M solution of H3PO4 from 5.2 mL of conc. acid (density, 1.88 g/cm3) in 1 L of water or 1 L of 0.1 M NaCl. A 0.1 M solution of borax was prepared from 38.1 g of Na2B4O7·10H2O in 1 L of water. The solution of NaCl was prepared from a 5.8 g weighed portion per 1 L of water.
Kaolinite, Al4Si4O10(OH), was taken from the Prosianovskoye field (Ukraine). It contains 39.5% Al2O3, 46.5% SiO2, and 14.0% H2O. Air-dried kaolinite was ground, sieved through a sieve (diameter, 1 mm), and stored in a stoppered polyethylene container.

2.4. Titration

All potentiometric titrations were made with a KEM AT-710M (Kyoto Electronic Manufacturing, Kyoto, Japan) automatic potentiometric titrator with a pH electrode (KEM Glass electrode C-171) and a CHA-700 autosampler (both, Kyoto Electronic Manufacturing, Kyoto, Japan). This instrument tracks the real-time response of the system to the addition of a titrant in various modes of titrant supply and equilibration time, with continuous output of results to screen and files. The instrument was operated in the incremental mode, when the titrant is added in small equal portions, with a given time to establish the equilibrium after each portion. Titration was carried out with a titrant aliquot of 0.02 mL, every 3 s, with stirring at a rate of 7 (within the instrument stirrer rate range from 1 to 10), with continuous pH measurement. The volume of 10 mL (25 mL for the Oe horizon of R II soil) of titrant was used, thus obtaining 502 pH points for a titration curve. The difference between adjacent points on the titration curve corresponds to hundredths of pH units. To immerse the working part of the instrument into the test sample, at least 80 mL of liquid is required. A single titration took 45 min.

2.4.1. Model Solutions

Verification of the mathematical program described previously [41] was carried out by the titration curves of model mixtures of polybasic acids (S1, 0.0125 M HCl; S2, 0.00625 M H3PO4 and 0.00625 M H3BO3; S3 0.00375 M H3PO4 and 0.00875 M H3BO3; and S4, 0.00875 M H3PO4 and 0.00375 M H3BO3). All solutions were prepared in 0.1 M NaCl. Titration of acid model mixtures was made according to the following scheme: to 10 mL of acid or a mixture of acids, 70 mL of 0.1 M NaCl was added, and the mixture was titrated with 0.1 M NaOH, prepared with 0.1 M NaCl.

2.4.2. Kaolinite

A weighed portion of kaolinite of 3 g was mixed with 80 mL of 0.1 M NaCl. Next, the electrode, the mechanical stirrer, and the feeding tube of the titrant were immersed into the solution, and the sample was measured for 5 s without titration to establish the equilibrium and to measure the initial pH value of the soil suspension with enough precision. The suspension was then titrated with 0.1 M NaOH with 0.1 M NaCl.

2.4.3. Soils

A weighed portion of the soil of 3 g was mixed with 80 mL of 0.1 M NaCl. Titration was made with 0.1 M NaOH prepared with 0.1 M NaCl, similarly to kaolinite titration. To determine hydrolytic acidity, 1 M CH3COONa extracts with pH 8.2 were prepared for a soil-to-solution ratio of 2:5, shaken for 1 h, filtered through a dry folded filter, and titrated with 0.1 M NaOH. Material from poorly mineralized forest litter was not analyzed.

2.5. Data Treatment

2.5.1. Titration-Curve Treatment

The algorithm for computation of the amounts of the acid groups in the soil samples was as follows [41]. Considering a sample as a mixture of monoprotic acids HAi, with acidity constants Ki, the charge neutrality of the solution titrated with NaOH is
i [ A ] = [ H + ] [ O H ] + [ N a + ] .
Equation (1) can be written as follows:
V 0 V 0 + V t i c i K i K i + H + = H + K w H + + c t V t V 0 + V t
Here, V 0 is the initial volume of the solution, V t is the titrant volume, K w is autoprotolysis constant of water, c t is the titrant concentration, and c i is the concentration of the ith acid. Specifying a set of Ki (i = 1, … m) and taking given values of c0, V0, c t and experimental values of V t and [H+] for n > m points of the titration curve, we get an overdetermined system of n equations of (2) linear with respect to c i solved by the algorithm of linear-regression analysis with restrictions on non-negativity of the solutions [51].
The data obtained in the course the titration was transferred from the AT-710 standalone unit to the computer. The computation algorithm, Equation (2), was implemented as Microsoft Excel scripts. The user imports the following data: the titration curve represented as two columns, the initial system volume, titrant concentration, and the value of the activity coefficient of a mono-charged species. The operator also enters the discrete pKa values, which will be included in the calculated pKa spectrum.
We used a pKa range of 2–13 for the model systems and 4–10 for real soil samples. The step was 1 or 0.5 [41,42,43,44]. Contents of groups with pKa equal to 5.0; 6.5; 8.5; 9.0; and 10.0; or ranges 4.0–5.0; 4.0–6.0; 5.5–6.5; 6.5–8.5; 7.5–8.5; and 9.0–10.0 were determined by potentiometric titration of all the horizons.

2.5.2. Principal Component Analysis

Principal component analysis was made with STATISTICA 8.0 software (TIBCO, Santa Clara, CA, USA). All the test variables were continuous. The original data matrix is centered and normalized to standard deviations. Prior to analysis, a Kaiser–Meyer–Olkin (KMO) test was used; the value was KMO = 0.6, which shows that the data are suited for factor analysis [52]. The decomposition was performed using NIPALS as the most standard algorithm providing the minimum numbers of ambiguous alarms and missed alarms [53,54], with the maximum number of iterations of 50 and a convergence criterion of 0.0001. The number of components describing the source set was selected using the cross-validation by the v-fold method; v was equal to 7, which provided the minimum error rate of the following PCA procedure and algorithm efficiency [55]. For most treatments, the first two principal components are selected, which describe 65% of the initial data.

3. Results and Discussion

3.1. Morphological and General Soils Properties

The studied Retisols have typical morphological and chemical properties [56]. The organic horizon was typical for boreal forests and was divided into three subhorizons with different organic matter decomposition levels: Oi, Oe, and Oa. An eluvial albic horizon (E) is located under the organic horizon. The upper part of E was infiltrated by total carbon (Eh). The transitional subeluvial horizon (BEL), with albeluvic glossae, was next. Clearly differentiated texture horizons (argic) occurred under the eluvial part of the soil profile. These had a typical multiordinal structure and numerous cutans at the faces of individual particles. R II soil is characterized by a significant increase in nodules in the middle part of the mineral horizons. This increase is probably associated with stem waters enriched in iron, aluminum, and organic matter, which entered the mineral horizons in significant quantities over several centuries [57]. The soil carbon concentration at organic horizons varied from 35.8 to 46.2% (Table 2). The upper mineral horizons contain significantly less carbon. The upper part of the eluvial horizon is characterized by a carbon content of 2.0 to 3.2%; in the lower horizons, the carbon concentration decreases to 0.1–0.9%. The lowest pH values were observed in the organic horizons, where the pH values varied from 4.0 to 5.2. The most acidic were the lower part of the organic horizon (Oa) and the upper part of the eluvial horizons (E). The studied soils have similar soil texture and a weak profile differentiation of texture fraction contents.

3.2. Titration-Curve Treatment Selection

When selecting a method of direct potentiometric titration for Retisol samples, the continuous variant [34,38,58,59] was selected over the equilibrium method [22,60], as the former is based on the sequential addition of the titrant to the test soil sample and, due to the shorter equilibration time, it makes it possible to investigate the composition of protolytic groups of the surface of primary soil particles [61].
The following methods of analyzing the potentiometric titration curves of soils are used: (i) the calculation of the area under the titration curve to estimate the buffering intensity and calculate the cation-exchange capacity [34,62]; (ii) simulation of the experimental data with titration curves, calculated for a set of possible individual soil components with known composition and acid–base properties [63]; and (iii) processing titration curves to find acidity constants of individual functional groups and to estimate their quantities. The latter two methods form complementary—direct and inverse —computational problems and usually serve as generalized models implemented in so-called acid–base calculators. Such calculators allow both simulating the titration curve for a given system and experimental technique and calculating the number of components different in acidity constant (strength) from the actual titration curve [64,65]. However, they reliably calculate the quantities and acidity constants only for a set of not more than four components, which is not suitable for soils. Volk and Jackson [66] proposed a model approach to identify four acid–base buffering ranges, each of which was attributed to specific soil components appearing in the solution in the corresponding pH range: H3O+ ions and monomeric Al3+ species, polymeric hydroxy-aluminum ions, and other components. The number of components used in the model approach limits the ability to interpret the composition of the system under study. Apart from this study, there are no specialized programs developed directly for the analysis of soil acidity.
A more general approach, which was implemented in this study, describes titration curves as arbitrary acid–base systems without a need of imposing restrictions on the number of components [41,42,43]. The result of the calculations is a pKa spectrum—the dependence of the number of groups on their pKa values specified with a certain ΔpKa step. Garmash and Vorob’eva [43] showed that the adequacy of the calculation results is achieved with a pKa step of at least 0.5.

3.3. Initial Tests of Model Systems

For the initial tests of the titration in the selected mode, the verification of the capabilities of the instrument, as well as the parameters of correctness and reproducibility, simple titration models were selected. They are acids and mixtures of acids of different strengths. We have selected HCl, H3BO3, and H3PO4 as acids with stepwise constants over a wide range. The choice of these model acids was also due to the presence of a wide acidity range of functional groups in soils [19,20,21,22,23,24]. The titration curves of mixtures of weak acids make it possible to check the correctness of the calculation algorithm [41], which was applied subsequently to the treatment curves titration of the soil samples. The titration of the soil samples was carried out using a saline solution with an ionic strength of 0.1 M because this ionic strength is close to that of soil solutions. The curves obtained are presented in Figure S2 (Supplementary information). For 0.0125 M HCl, the concordant concentration value of 0.01245 M was calculated. Table 3 summarizes all non-zero values of pKa calculated for H3BO3/H3PO4 mixtures with a step of 0.5. The second line for each mixture in Table 3 shows (in bold) the total concentration of groups with adjacent pKa.
The value for pKa of 2.5–3.0 is concordant with the concentration of phosphoric acid in the first dissociation stage (pKa 2.15). The amounts for pKa of 7.0 and 7.5 are in agreement with the concentration of phosphoric acid in the second dissociation stage (pKa 7.2). Finally, a sum of concentrations for pKa of 9.0 and 9.5 describes boric acid by its first dissociation stage (pKa 9.15). Thus, in all the cases, the calculations agree well with the composition of the model mixtures (Table 3).

3.4. Kaolinite

Kaolinite, as a mineral of the layered silicate subclass, a widely distributed secondary mineral present in the clay fraction of many soils [67,68], previously applied to model soils [69], was a transitional step from analyzing aqueous acid mixtures to soil analysis. For a ΔpKa step of 1, the pKa range values are 6.0–7.0 and 9.0–10.0 (Figure 2a), which agrees with existing the acidity constants of kaolinite of 6.9 [70] and 8.0 [71]. However, despite the fact that the titration curves of kaolinite (Figure S3, Supplementary Information) differ insignificantly visually, the calculations at this ΔpKa step result in a significant difference between replicate curves, which is inappropriate. As ΔpKa affects the accuracy of processing, a smaller division was tested. Garmash, Ustimova, Kudryavtsev, Vorob’eva, and Polenova [41] recommended a minimum ΔpKa step of 0.5. At a ΔpKa of 0.5 applied for the same experimental curves (Figure 2b), three types of groups with pKa ranges of 6.5–7.0; 8.5–9.5; and 10.0 can be identified clearly and reproducibly.
At the same time, the proposed method of titration-curve processing provided a larger set of acidity constants for this mineral. However, the ΔpKa step cannot be reduced indefinitely to obtain a continuous distribution of acidic or basic components, as it significantly increases the calculation error, and the variation of results becomes no longer significant [41]. In our case, at a ΔpKa of 0.25, noise in the data made it impossible to identify specific groups on the surface (Figure 2c). As soils are more complex compared to kaolinite, a ΔpKa of 0.5 was used for this point on.

3.5. Soil Titration

For model solutions and even kaolinite, we mostly used individual pKa values. However, for soils, it is more expedient to consider pKa ranges because in real soils no group can be described by a single pKa value due to intergroup interactions [2]. According to the state-of-the-art concepts, characteristic pH ranges can be distinguished (Table 4). Previously, pKa spectroscopy has been successfully used to study the acid–base properties of humic [72] and fulvic acids [73], γ-Al2O3 suspensions [74], benthal sediments [75,76], and aqueous soil extracts [29]. A study was conducted on the determination of the pKa spectra of mineral horizons of acid soils of the Republic of Komi [77], which identified four to five groups involved in protolytic equilibria, forming acid–base centers at the solid-phase surface of soils. In this study, the method of equilibrium potentiometric titration was used, and the correlation coefficients between physical and chemical parameters (exchange and non-exchange acidity, Corg) and pKa were found.

3.5.1. Retisol I, Intercrown (Open) Space

To obtain the titration curves for the R I soil samples (Figure S4, Supporting Information), 10 to 25 mL of titrant was taken, depending on the horizon. Figure 3a shows the total distribution of protolytic groups, which may indicate a low humification degree of organic matter. This can be indirectly determined by the C/N ratio: in the upper horizons of this soil, C/N reaches 30, while in the lower horizons it decreases to 5–6. Low values of the indicator indicate the prevailing processes of mineralization of newly arriving organic matter, while high values indicate a slowdown in the decomposition processes in the conditions of boreal forests. At a C/N below 6, humic acids prevail; otherwise, fulvic and labile organic matter predominate. Previously, a high diversity of hydroxy acids was noted in the composition of the litter of podzolic soils in this region [86], which indicates a large number of diverse functional groups, associated mainly with aromatic acids with a wide range of pKa values, from 5 to 10 in the upper part of the profile and from 8 to 10 in the lower part. In general, the distribution of groups also corresponds to the amounts of oxalate-soluble forms of Fe and Al.
The eluvial horizon E1 shows an expected decrease in the number of protolytic groups. Thus, the concentrations of the groups with pKa values of 5.0, 6.5, and 8.0 sharply decrease due to the removal of the active material in the lower horizons; this is implemented with appearance of groups with pKa 5.5, 6.0, 8.0, and 9.0. The next horizon, Eh, is characterized by humus–ferrous films on the surface of mineral grains (i.e., non-silicate forms of sesquioxides and humic substances) and high contents of groups with pKa values of 5.0, 5.5, 6.5, 7.0, 8.0, 8.5, 9.0, and 9.5 (Figure 4a).
The distribution of protolytic functional groups along the soil profile reflects the characteristics of the substance of its genetic horizons. Thus, a high content of components with a pKa value of 9.5 in Eh was noted in the composition of the humified ferruginous eluvial–illuvial horizon. The content of more mobile components with a pKa value of 5.5 in the middle part of the profile is higher than that in the topsoil. This is associated with the migration of easily mobile oxalate-soluble aluminum and iron (Table 1 and Figure 5a). The crystallized species of iron oxides and phosphate compounds with a pKa value of 9.0 (Figure 3a) are distributed by the eluvial–illuvial type, with their accumulation in the illuvial iron-rich horizon of the microprofile and in the lower part of the whole profile.

3.5.2. Retisol II (Under the Crown of Spruce)

The titration of suspensions of this soil type required 10 mL of the titrant for all horizons (Figure S5, Supplementary Information). The distributions of Fe, Al, and P in the profile are shown in Figure 5b. The physicochemical properties of R II are presented in Table 2. The distribution of acidic components is eluvial–illuvial, with accumulation in the middle part of the profile in the iron-rich transitional horizon (Figure 3b). This is due to the predominance of weak groups with pKa values of 8.5–10.0 (OH groups of dissolved phenols, polyphenols, and their hydrolysis products, hydroxy–organic complexes of Al and Fe, Table 4). Their contents decrease down the profile and are practically aligned in illuvial, texturally differentiated horizons (Figure 4b).
The distribution of stronger carboxyl groups (pKa 4.0–4.5) is accumulative, with the most accumulation in the upper eluvial horizon. The illuvial accumulation of groups with pKa values of 5.0–5.5 down the profile corresponds to the vertical distribution of oxalate-soluble aluminum. In contrast to the R I soil, the distribution of acidity groups along the profile is more uniform in the R II soil, which we associate with a heavier particle size distribution, causing the formation of a pronounced soil particle structure, and less intense vertical migration of the mobile organic, organomineral, and amorphous mineral compounds. Hydrolytic (pH-dependent) acidity is associated with soil components that exhibit the properties of weak acids in general, i.e., not only exchange ions, but also components associated with varying surface charges.
Retisol II is distinguished by a high content of phenolic groups, which is typical for spruce secretions [87] and Fe and Al (columns with pKa 10.0, 8.5 and 5.0). This is due to the chemical composition of the aqueous solutions reaching the soil surface [88]. Thus, Robakidze, Torlopova, and Bobkova [57] showed that in the waters in tree crowns, the concentration of dissolved carbon, iron, and aluminum is significantly higher.

3.6. Principal Component Analysis and Classification

Each of the studied samples was described using more than 35 factors in total (Table 1 and Table 2 and Section 2.3), among which 11 are related to the content of groups with pKa values from the titration curves: 5.0; 6.5; 8.5; 9.0; and 10.0; or ranges 4.0–6.0; 6.5–8.5; 9.0–10.0; 4.0–5.0; 5.5–6.5; and 7.5–8.5. As the number of variables is large, the principal component method was used [89,90,91]. The second reason for PCA application is that it provides a complex view on the properties of soil horizons. In fact, the contents of groups with given pKa values are already the representative parameters, which can be used to study soils (as discussed in the Introduction), but PCA makes it possible to take into account each of these parameters along with others, as well as with conventional soil parameters such as the chemical composition.
Each horizon of two tested soils is described by two principal components: PC1 and PC2 (Table S1, Supplementary Information). It is possible to identify separate regions corresponding to the horizons of R I and R II soils (Figure 6). Therefore, the PCA sample is representative, and the data obtained are enough for an adequate description of the horizons. The Eh horizon of R I falls into the region corresponding to R II, as the properties of this horizon are close to the properties of the E and Eh,c R II horizons. The loading values are summed up in Table S3 (Supplementary Information). The loading diagram is shown in Figure 7. It reflects the known relationships between the properties of soil horizons: strong positive correlations between the total C and N, Al2O3, and Fe2O3, and negative correlation with exchangeable Ca and Mg [89,92,93,94]. There is also a positive correlation between the number of main groups and hydrolytic acidity (pKa 4.0–5.0, 7.5–8.5, and 10.0). It is also clear that the contents of exchangeable calcium and magnesium are associated with the degree of saturation with bases and the pH values of both aqueous and saline solutions (to a lesser extent).
The contents of carbon and nitrogen are associated with a pKa value of 6.5, and the content of sulfur oxide as sulfur is mainly of organic origin in these soils. The C/N ratio, which is an indicator of the humification of soil organic matter, is associated with pKa 9.0–10.0 and, to a lesser extent, with pKa values of 6.5–8.5 and 4.0–6.0. A pKa value of 8.5 is virtually identical to oxalate-soluble alumina, while oxalate-soluble iron oxide is associated with a wide pKa range from 4.0 to 8.5. The correct reflection of these known dependences confirms the correctness of the model. Thus, it can be used to search for new dependences, related to the content of protolytic groups.
It should be noted that positive or negative correlations expected from the loadings plot might be strong or weak as well, or there might be no correlation at all. Thus, we examined all the pairs of the variables for which correlation is expected.

3.7. Correlation of pKa and Soil Indicators

As expected, the single values of pKa did not provide much data, and more information can be mined from pKa ranges (Table 5). This made it possible to indicate which of the conventional soil parameters are responsible for the formation of pKa spectra for soils of the same type but developed under different vegetation.
The two studied soils differ in the degree of correlation between the carbon and nitrogen content and pKa. In the R I soil, formed under the influence of diverse herbaceous vegetation in the intercrown space, correlations with the content of C and N were found, as well as with the degree of humification of C/N for almost all individual pKa values (except for a pKa value of 10, where there is no correlation, and for pKa values of 5.5–6.5, where the correlation is medium and weak). Despite the fact that soil organic matter accumulates in the composition of plant litter in the upper thin layer, the diversity of its composition and the presence of a significant amount of aliphatic compounds provides its mobility for migration along the profile [95,96]. In the R II soil, however, we found correlations between soil carbon and nitrogen with pKa values of 6.5 and 10 only. This soil is developed in the canopy space of spruce; in fact, it was formed under the influence of litter of a constant and homogeneous composition, in the absence of species diversity, mainly under the influence of polyphenolic compounds and their functional groups.
More correlations are found for oxalate-soluble Al2O3 and Fe2O3. The first has a negative correlation with pKa values of 5.0–6.5, 5.5–6.5, and 8.5. This confirms the existing data: the hydrolysis of exchangeable aluminum and polymeric hydroxy complexes of aluminum occur at pH 6.5 [79,80] and 8.5 [14], respectively. The content of oxalate-soluble Fe2O3 has correlations with pKa ranges of 4.0–6.0; 6.5–8.5; and 7.5–8.5, but the strongest one takes place with pKa 9.0–10.0. This result fits well with the hydrolysis of Fe complexes with organic ligands at pH 9.5–10.0 [85].
In addition to sulfur, which strongly correlates with the content and composition of organic matter (carbon, nitrogen, and their ratio), only phosphorus compounds reveal acidic properties in the studied soils. The P2O5 content has positive correlations with pKa values of 4.0–6.0 and 6.5–8.5, which should be attributed to the first two acidic constants of phosphoric acid [97]. In addition, no correlations of phosphorus with organic species should be expected because the acidity of such type of soils does not depend on the content of organic substances, as mentioned above. The acidity ranges 4.0–5.0 and 5.5–6.5 provided no useful information because the amounts of these groups were low.

4. Conclusions

In conclusion, potentiometric spectrotitrimetry at a fast high-resolution mode was applied to find the detailed quantitative distribution of surface protolytic groups in primary soil particles, and the results are directly connected with the soil chemical composition. Due to the high resolution, the error of assessing acidity constants from titration curves is decreased to 3–5% compared to commonly used titrimetry approaches. It is valuable that the approach provides the complete set of major acidity values (pKa spectra) without a priori information on a soil sample. The obtained pKa spectra can be used to distinguish the soils with matching properties and, thus, can be used to identify soils that coincide in properties and genesis but develop under the influence of different types of vegetation.
The shortcoming of this approach is that one has to choose the range of calculated pKa spectra and the step of pKa values, but this does not impose serious difficulties for soil titration because the soil acidic components do not require the precise values of the dissociation constants, but rather acidity range descriptors. Due to the fact that the organic matter entering the soil is always characterized by a high complexity of the composition of compounds, we recommend the analysis of the following pKa values: 4.0–5.0; 4.0–6.0; 5.5–6.5; 6.5–8.5; 7.5–8.5; and 9.0–10.0.
The information about the amounts of the protolytic groups in Retisol soil horizons was treated in two ways. Firstly, the technique allowed us to find that the pKa spectra do depend on the soil chemical composition for both soil types and enabled us to find the differences in the quantitative and qualitative distribution of protolytic groups of different strength both within each soil profile and between the two Retisols under different plant covers. Secondly, the results of the correlation analyses agree with the known soil chemistry and with the distribution of the protolytic groups through the horizons. This feature may be of value for screening and identification applications. We believe that fast high-resolution pKa spectrotitrimetry, guided by chemometric approaches, can be used for various agrochemical, environmental, and genetic studies of acidic soil and soil horizons.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/soilsystems8020063/s1, Figure S1: The location of the soil profiles: 61°45′43.4″ N 54°17′50.3″ E. Figure S2: Titration of acid mixtures with 0.1 NaOH. Explanations are in the text. Figure S3: Titration curves of kaolinite. Figure S4: Titration curves of Retisol I. Figure S5: Titration curves of Retisol II. Figure S6: Correlation coefficients between indicators of acid properties, general physicochemical characteristics of soils and pKa. Table S1: Principal component values for the PCA analysis of soil horizons. Table S2: Bulk content of main elements in soils. Table S3: Loading values for PCA analysis of soil.

Author Contributions

Conceptualization, D.S.V., A.V.G. and O.B.R.; methodology, A.V.G.; software, A.V.G.; validation, A.V.G. and M.A.P.; formal analysis, M.A.S.; investigation, N.V.M., O.B.R. and A.A.D.; resources, A.A.D.; data curation, N.V.M.; writing—original draft preparation, M.A.S. and A.V.G.; writing—review and editing, O.B.R. and M.A.P.; visualization, M.A.S.; supervision, M.A.P.; project administration, D.S.V.; funding acquisition, M.A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Higher Education of the Russian Federation M.V. Lomonosov Moscow State University No. AAAA-A21-121011590089-1 “Development of highly informative and high-tech methods of chemical analysis for the protection of ecosystems, the creation of new materials and advanced production technologies, the transition to environmentally friendly and resource-saving energy, the development of nature-like technologies, high-tech healthcare and rational use of natural resources”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Shoba, S.A. National Soil Atlas of Russia = Russian Soil Atlas = Soil Atlas of Russian Federation; Astrel: Moscow, Russia, 2011. [Google Scholar]
  2. Sokolova, T.A.; Trofimov, S.Y. Sorption Properties of Soils. Fdsorption. Cation Exchange; University Book: Moscow, Russia, 2009. [Google Scholar]
  3. Thompson, A.; Goyne, K.W. Introduction to the sorption of chemical constituents in soils. Nat. Educ. Knowl. 2012, 4, 7. [Google Scholar]
  4. Huang, P.M.; Wang, S.-L.; Tzou, Y.-M.; Huang, Y.; Weng, B.; Zhuang, S.; Wang, M.K. Physicochemical and biological interfacial interactions: Impacts on soil ecosystem and biodiversity. Environ. Earth Sci. 2012, 68, 2199–2209. [Google Scholar] [CrossRef]
  5. Audette, Y.; Congreves, K.A.; Schneider, K.; Zaro, G.C.; Nunes, A.L.P.; Zhang, H.; Voroney, R.P. The effect of agroecosystem management on the distribution of C functional groups in soil organic matter: A review. Biol. Fertil. Soils 2021, 57, 881–894. [Google Scholar] [CrossRef] [PubMed]
  6. Demydenko, O.; Bulygin, S.; Velychko, V.; Kaminsky, V.; Tkachenko, M. Soil moisture potential of agrocenoses in the Forest-Steppe of Ukraine. Agric. Sci. Pract. 2021, 8, 49–61. [Google Scholar] [CrossRef]
  7. Schindler, M. The structural and compositional complexity of surfaces on minerals and organic matter in soils. Abstr. Pap. Am. Chem. Soc. 2018, 256, 2. [Google Scholar]
  8. Kleber, M.; Sollins, P.; Sutton, R. A conceptual model of organo-mineral interactions in soils: Self-assembly of organic molecular fragments into zonal structures on mineral surfaces. Biogeochemistry 2007, 85, 9–24. [Google Scholar] [CrossRef]
  9. Jeon, I.; Nam, K. Change in the site density and surface acidity of clay minerals by acid or alkali spills and its effect on pH buffering capacity. Sci. Rep. 2019, 9, 9878. [Google Scholar] [CrossRef] [PubMed]
  10. Gang, Y.; Li, Q.; Li, H. Measurement of Surface Charges and Mechanism of Interfacial Processes for Soil Clay Minerals. Eurasian Soil Sci. 2021, 54, 1546–1563. [Google Scholar] [CrossRef]
  11. Hao, W.; Flynn, S.L.; Kashiwabara, T.; Alam, M.S.; Bandara, S.; Swaren, L.; Robbins, L.J.; Alessi, D.S.; Konhauser, K.O. The impact of ionic strength on the proton reactivity of clay minerals. Chem. Geol. 2019, 529, 119294. [Google Scholar] [CrossRef]
  12. Yang, Y.; Wang, Y.; Peng, Y.; Cheng, P.; Li, F.; Liu, T. Acid-base buffering characteristics of non-calcareous soils: Correlation with physicochemical properties and surface complexation constants. Geoderma 2020, 360, 114005. [Google Scholar] [CrossRef]
  13. Sposito, G. The Operational Definition of the Zero Point of Charge in Soils. Soil Sci. Soc. Am. J. 1981, 45, 292–297. [Google Scholar] [CrossRef]
  14. Wang, X.; Tang, C.; Mahony, S.; Baldock, J.A.; Butterly, C.R. Factors affecting the measurement of soil pH buffer capacity: Approaches to optimize the methods. Eur. J. Soil Sci. 2014, 66, 53–64. [Google Scholar] [CrossRef]
  15. Shrestha, B.M.; Certini, G.; Forte, C.; Singh, B.R. Soil Organic Matter Quality under Different Land Uses in a Mountain Watershed of Nepal. Soil Sci. Soc. Am. J. 2008, 72, 1563–1569. [Google Scholar] [CrossRef]
  16. Pinsky, D.L. Ion-Exchange Processes in Soils; Russian Academy of Sciences: Pushchino, Russia, 1997; p. 166. [Google Scholar]
  17. Ulrich, B.; Sumner, M.E. Soil Acidity; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  18. Kelley, A.P. Soil Acidity, an Ecological Factor; University of Pennsylvania: Philadelphia, PA, USA, 1923. [Google Scholar]
  19. Raveh, A.; Avnimelech, Y. Potentiometric Determination of Soil Organic Matter. Soil Sci. Soc. Am. J. 1972, 36, 967–968. [Google Scholar] [CrossRef]
  20. Magdoff, F.R.; Bartlett, R.J. Soil pH Buffering Revisited. Soil Sci. Soc. Am. J. 1985, 49, 145–148. [Google Scholar] [CrossRef]
  21. Bradfield, R. The Nature of the Acidity of the Colloidal Clay of Acid Soils. J. Am. Chem. Soc. 2002, 45, 2669–2678. [Google Scholar] [CrossRef]
  22. Puri, A.N.; Asghar, A.G. Titration Curve and Dissociation Constants of Soil Acidoids. Soil Sci. 1938, 45, 359–368. [Google Scholar] [CrossRef]
  23. Simon, B.; Tolner, L.; Rékási, M.; Michéli, E. Soil Acidity Investigation by Potentiometric Titrations. Cereal Res. Commun. 2006, 34, 283–286. [Google Scholar] [CrossRef]
  24. Xiaocui, W.E.N.; Jiuyu, L.I.; Jie, S.; Liang, T.A.O. Research Progress on the Acid-base Properties of Variable Charge Soils Using Potentiometric Titration. Acta Pedol. Sin. 2022, 59, 910–923. [Google Scholar] [CrossRef]
  25. Food and Agriculture Organization of the United Nations. World Reference Base for Soil Resources; Food & Agriculture Org.: Rome, Italy, 1998; Volume 3. [Google Scholar]
  26. Shishov, L.; Tonkonogov, V.; Lebedeva, I.; Gerasimova, M. Classification and Diagnostics of Soils of Russia; Oikumena: Smolensk, Russia, 2004. (In Russian) [Google Scholar]
  27. Nachtergaele, F. The “soils” to be classified in the World Reference Base for Soil Resources. Eurasian Soil Sci. 2005, 38, S13–S19. [Google Scholar]
  28. Sokolova, T.A.; Tolpeshta, I.I.; Rusakova, E.S. Contributions of separate reactions to the acid–base buffering of soils in brook floodplains (Central Forest State Reserve). Eurasian Soil Sci. 2016, 49, 399–411. [Google Scholar] [CrossRef]
  29. Shamrikova, E.V.; Ryazanov, M.A.; Vanchikova, E.V. Acid-base properties of water-soluble organic matter of forest soils, studied by the pK-spectroscopy method. Chemosphere 2006, 65, 1426–1431. [Google Scholar] [CrossRef] [PubMed]
  30. Cheng, P.; Wang, Y.; Cheng, K.; Li, F.; Qin, H.; Liu, T. The Acid-Base Buffer Capacity of Red Soil Variable Charge Minerals and Its Surface Complexation Model. Acta Chim. Sin. 2017, 75, 637–644. [Google Scholar] [CrossRef]
  31. Nelson, P.N.; Su, N. Soil pH buffering capacity: A descriptive function and its application to some acidic tropical soils. Soil Res. 2010, 48, 201–207. [Google Scholar] [CrossRef]
  32. Xu, R.-K.; Zhao, A.-Z.; Yuan, J.-H.; Jiang, J. pH buffering capacity of acid soils from tropical and subtropical regions of China as influenced by incorporation of crop straw biochars. J. Soils Sediments 2012, 12, 494–502. [Google Scholar] [CrossRef]
  33. Kissel, D.E.; Sonon, L.S.; Cabrera, M.L. Rapid Measurement of Soil pH Buffering Capacity. Soil Sci. Soc. Am. J. 2012, 76, 694–699. [Google Scholar] [CrossRef]
  34. Maksimova, Y.G.; Maryakhina, N.N.; Tolpeshta, I.I.; Sokolova, T.A. The acid-base buffer capacity of podzolic soils and its changes under the impact of treatment with the Mehra-Jackson and Tamm reagents. Eurasian Soil Sci. 2010, 43, 1120–1131. [Google Scholar] [CrossRef]
  35. Shamrikova, E.V.; Sokolova, T.A.; Zaboeva, I.V. Identification of buffer reactions occurring in the course of acid-base titration of water suspensions from virgin and plowed podzolic soils. Eurasian Soil Sci. 2002, 35, 363–373. [Google Scholar]
  36. Hartikainen, H. Acid- and base titration behaviour of finnish mineral soils. Z. Pflanzen. Bodenk. 2007, 149, 522–532. [Google Scholar] [CrossRef]
  37. Kuznetsov, N.B.; Alekseeva, S.A.; Shashkova, G.V.; Dronova, T.Y.; Sokolova, T.A. Buffer capacities of podzolic and peat gleyic podzolic soils to sulfuric and nitric acids. Eurasian Soil Sci. 2007, 40, 353–366. [Google Scholar] [CrossRef]
  38. Rusakova, E.S.; Ishkova, I.V.; Tolpeshta, I.I.; Sokolova, T.A. Acid-base buffering of soils in transitional and transitional-accumulative positions of undisturbed southern-taiga landscapes. Eurasian Soil Sci. 2012, 45, 503–513. [Google Scholar] [CrossRef]
  39. Aringhieri, R.; Pardini, G. Protonation of Soil Particle Surfaces: Kinetics. Can. J. Soil Sci. 1985, 65, 677–686. [Google Scholar] [CrossRef]
  40. Sakurai, K.; Ohdate, Y.; Kyuma, K. Comparison of salt titration and potentiometric titration methods for the determination of zero point of charge (ZPC). Soil Sci. Plant Nutr. 1988, 34, 171–182. [Google Scholar] [CrossRef]
  41. Garmash, A.V.; Ustimova, I.V.; Kudryavtsev, A.V.; Vorob’eva, O.N.; Polenova, T.V. Potentiometric analysis of complex protolytic systems by pK spectroscopy using linear regression. J. Anal. Chem. 1998, 53, 212–218. [Google Scholar]
  42. Garmash, A.V.; Vorob’eva, O.N.; Kudryavtsev, A.V.; Danchenko, N.N. Potentiometric analysis of polyelectrolytes by pK spectroscopy using linear regression analysis. J. Anal. Chem. 1998, 53, 361–367. [Google Scholar]
  43. Garmash, A.V.; Vorob’eva, O.N. Verification of the accuracy of results in the potentiometric analysis of polyelectrolytes by pK spectroscopy. J. Anal. Chem. 1998, 53, 228–233. [Google Scholar]
  44. Lodygin, E.; Shamrikova, E. Use of the pK Spectroscopy Method in the Study of Protolytic Properties of Humic Substances and Other Soil Polyelectrolytes. Agronomy 2021, 11, 1051. [Google Scholar] [CrossRef]
  45. Field Guide to Russian Soils; V.V. Dokuchaev Soil Institute: Moscow, Russia, 2008. (In Russian)
  46. Taskaev, A.I. (Ed.) Atlas Pochv Respubliki Komi; Komi Respublikanskaya Tipografiya: Syktyvkar, Russia, 1997; p. 356. [Google Scholar]
  47. Dymov, A.A. Soils of Native Forest Ecosystems. Eurasian Soil Sci. 2023, 56, S36–S45. [Google Scholar] [CrossRef]
  48. van Reeuwijk, L.P. Technical Paper 09: Procedures for Soil Analysis (6th Edition). In Technical Paper/International Soil Reference and Information Centre; van Reeuwijk, L.P., Ed.; ISRIC: Wageningen, The Netherlands, 2002; Volume 9. [Google Scholar]
  49. Tamm, O. Method of estimation of the inorganic compounds of the gel complex in soils. Medd. Statens. Skoysgforsksanstant. 1922, 19, 387–404. [Google Scholar]
  50. Mehra, O.P.; Jackson, M.L. Iron oxid removal from soils and clays by ditionite-citrate system buffered with sodium bicarbonate. Clays Clay Miner. 1960, 7, 317–327. [Google Scholar]
  51. Lawson, C.L.; Hanson, R.J. Solving Least Squares Problems; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1995. [Google Scholar]
  52. Kaiser, H.F.; Rice, J. Little Jiffy, Mark Iv. Educ. Psychol. Meas. 1974, 34, 111–117. [Google Scholar] [CrossRef]
  53. Li, G.; Qin, S.J.; Zhou, D. Geometric properties of partial least squares for process monitoring. Automatica 2010, 46, 204–210. [Google Scholar] [CrossRef]
  54. Ochoa-Muñoz, A.F.; Contreras-Reyes, J.E. Multiple Factor Analysis Based on NIPALS Algorithm to Solve Missing Data Problems. Algorithms 2023, 16, 457. [Google Scholar] [CrossRef]
  55. Arlot, S.; Celisse, A. A survey of cross-validation procedures for model selection. Stat. Surv. 2010, 4, 40–79. [Google Scholar] [CrossRef]
  56. Startsev, V.V.; Yakovleva, E.V.; Kutyavin, I.N.; Dymov, A.A. Fire Impact on Carbon Pools and Basic Properties of Retisols in Native Spruce Forests of the European North and Central Siberia of Russia. Forests 2022, 13, 1135. [Google Scholar] [CrossRef]
  57. Robakidze, E.A.; Torlopova, N.V.; Bobkova, K.S. Chemical composition of wet precipitation in old-growth middle-taiga spruce stands. Geochem. Int. 2013, 51, 65–75. [Google Scholar] [CrossRef]
  58. López, R.; Gondar, D.; Antelo, J.; Fiol, S.; Arce, F. Study of the acid-base properties of a peat soil and its humin and humic acid fractions. Eur. J. Soil Sci. 2012, 63, 487–494. [Google Scholar] [CrossRef]
  59. Fernández, J.M.; Plaza, C.; Senesi, N.; Polo, A. Corrigendum to: “Acid–base properties of humic substances from composted and thermally-dried sewage sludges and amended soils as determined by potentiometric titration and the NICA-Donnan model” [Chemosphere 69(4) (2007) 630–635]. Chemosphere 2012, 87, 586. [Google Scholar] [CrossRef]
  60. Tipping, E.; Woof, C. Humic substances in acid organic soils: Modelling their release to the soil solution in terms of humic charge. J. Soil Sci. 2006, 41, 573–586. [Google Scholar] [CrossRef]
  61. Skic, K.; Boguta, P.; Sokołowska, Z. Analysis of the sorption properties of different soils using water vapour adsorption and potentiometric titration methods. Int. Agrophys. 2016, 30, 369–374. [Google Scholar] [CrossRef]
  62. Sokolova, T.; Pakhomov, A.; Terekhin, V. Continuous potentiometric titration measurement of acid-base buffering of podzolic soils. Eurasian Soil Sci. 1994, 26, 91–105. [Google Scholar]
  63. Reis, C.D.G.; Fonseca, R.A.D.; Reis, E.L.; Reis, C.; Brighenti, C.R.G.; Matias, A.A. Programa de ajuste multiparamétrico de curvas de titulação potenciométricas de ácidos húmicos. Rev. Bras. Ciência do Solo 2010, 34, 569–573. [Google Scholar] [CrossRef]
  64. Gutz, I.G.R. CurTiPot—pH and Acid-Base Titration Curves: Analisis and Simulation Software, Version 4.3.1. Available online: http://www.iq.usp.br/gutz/Curtipot_.html (accessed on 23 July 2023).
  65. Benjamin, F.T. ProtoFit—A Program for Optimizing Surface Protonation Models. Available online: http://protofit.sourceforge.net/ (accessed on 26 February 2019).
  66. Volk, V.V.; Jackson, M.L. Inorganic pH Dependent Cation Exchange Charge of Soils. Clays Clay Miner. 1963, 12, 281–295. [Google Scholar] [CrossRef]
  67. Bergaya, F.; Theng, B.K.G.; Lagaly, G. Handbook of Clay Science; Elsevier Science Bv: Amsterdam, The Netherlands, 2006; pp. 1–1224. [Google Scholar]
  68. Young, R.A.; Hewat, A.W. Verification of the triclinic crystal structure of kaolinite. Clays Clay Miner. 1998, 36, 225–232. [Google Scholar] [CrossRef]
  69. Kriaa, A.; Hamdi, N.; Srasra, E. Determination of Point of Zero Charge of Tunisian Kaolinites by Potentiometric and Mass Titration Methods. J. Chin. Chem. Soc. 2008, 55, 53–61. [Google Scholar] [CrossRef]
  70. Liu, X.; Lu, X.; Sprik, M.; Cheng, J.; Meijer, E.J.; Wang, R. Acidity of edge surface sites of montmorillonite and kaolinite. Geochim. Cosmochim. Acta 2013, 117, 180–190. [Google Scholar] [CrossRef]
  71. Furda, L.V.; Lebedeva, O.E. The influence of acid-modified montmorillonite containing clays on thermal destruction of polyethylene. Nauk. Pr. DonNTU 2014, 2, 178–184. [Google Scholar]
  72. Danchenko, N.N.; Perminova, I.V.; Garmash, A.V.; Kudryavtsev, A.V. Determination of carboxyl acidity of humic substances by titrimetric methods. Vestn. Mosk. Univ. Seriya 2 Khimiya 1998, 39, 127–131. [Google Scholar]
  73. Ryazanov, M.; Zlobin, D.; Lodygin, E.; Beznosikov, V. Evaluation of the acid-base properties of fulvic acids using pK spectroscopy. Eurasian Soil Sci. 2001, 34, 830–836. [Google Scholar]
  74. Ryazanov, M.A.; Dudkin, B.N. Acid–Base Properties of γ-Al2O3 Suspension Studied by pK Spectroscopy. Colloid J. 2003, 65, 761–766. [Google Scholar] [CrossRef]
  75. Gulmini, M.; Zelano, V.; Daniele, P.G.; Prenesti, E.; Ostacoli, G. Acid-base properties of a river sediment: Applicability of potentiometric titrations. Anal. Chim. Acta 1996, 329, 33–39. [Google Scholar] [CrossRef]
  76. Kremleva, T.A.; Lebedeva, N.N.; Allayarov, D.A.; Novopashin, V.V.; Volchkov, A.V. Determination of the Specific Surface Charge of the Benthal Deposits by the Potentiometric Method. Vestn. Tyumen State Univ. Seriya Khimiya 2014, 5, 63–72. [Google Scholar]
  77. Shamrikova, E.V.; Vanchikova, E.V.; Ryazanov, M.A. Study of the acid-base properties of mineral soil horizons using pK spectroscopy. Eurasian Soil Sci. 2007, 40, 1169–1174. [Google Scholar] [CrossRef]
  78. Rich, C.I. Conductometric and Potentiometric Titration of Exchangeable Aluminum. Soil Sci. Soc. Am. J. 1970, 34, 31–38. [Google Scholar] [CrossRef]
  79. Bergelin, A.; van Hees, P.A.W.; Wahlberg, O.; Lundström, U.S. The acid–base properties of high and low molecular weight organic acids in soil solutions of podzolic soils. Geoderma 2000, 94, 223–235. [Google Scholar] [CrossRef]
  80. Kapoor, B.S.; Goswami, S.C. Acid Character of Dioctahedral Vermiculite: Permanent and pH-Dependent Charge Component of Cation Exchange Capacity. J. Ind. Soc. Soil Sci. 1978, 24, 332–338. [Google Scholar]
  81. Öhman, L.-O.; Sjöberg, S.; Melander, L.; Wahlgren, U.; Watt, C.I.F. Equilibrium and Structural Studies of Silicon(IV) and Aluminium(III) in Aqueous Solution. 8. A Potentiometric Study of Aluminium(III) Salicylates and Aluminium(III) Hydroxo Salicylates in 0.6 M Na(Cl). Acta Chem. Scand. 1983, 37a, 875–880. [Google Scholar] [CrossRef]
  82. Violante, P.; Violante, A. Characterization of H-bentonites and interpretation of the 3rd buffer range appearing on pH titration curves. Agrochimica 1980, 24, 47–58. [Google Scholar]
  83. Orlov, D.S. Humic Acids and General Tejry of Gumification; MSU: East Lansing, MI, USA, 1990; p. 325. [Google Scholar]
  84. Makarov, M.I. Soil Organic Phosphorus; GEOS: Moscow, Russia, 2009; p. 397. [Google Scholar]
  85. Evans, L.J. Some aspects of the chemistry of aluminium in podzolic soils. Commun. Soil Sci. Plant Anal. 2008, 19, 793–803. [Google Scholar] [CrossRef]
  86. Shamrikova, E.V.; Gruzdev, I.V.; Punegov, V.V.; Khabibullina, F.M.; Kubik, O.S. Water-soluble low-molecular-weight organic acids in automorphic loamy soils of the tundra and taiga zones. Eurasian Soil Sci. 2013, 46, 654–659. [Google Scholar] [CrossRef]
  87. Kerienė, I.; Šaulienė, I.; Šukienė, L.; Judžentienė, A.; Ligor, M.; Valiuškevičius, G.; Grendaitė, D.; Buszewski, B. Enrichment of Water Bodies with Phenolic Compounds Released from Betula and Pinus Pollen in Surface Water. Plants 2024, 13, 99. [Google Scholar] [CrossRef] [PubMed]
  88. Ryazanov, M.A.; Dudkin, B.N. The pK spectroscopy study of the acid–base properties of hydrated aluminum oxide sols. Colloid J. 2004, 66, 726–728. [Google Scholar] [CrossRef]
  89. Borůvka, L.; Mládková, L.; Penížek, V.; Drábek, O.; Vašát, R. Forest soil acidification assessment using principal component analysis and geostatistics. Geoderma 2007, 140, 374–382. [Google Scholar] [CrossRef]
  90. Salehi, A.; Zahedi Amiri, G. Study of Physical and Chemical Soil Properties Variations Using Principal Component Analysis Method in the Forest, North of Iran. Caspian J. Env. Sci. 2005, 3, 131–137. [Google Scholar]
  91. Visconti, F.; de Paz, J.M.; Rubio, J.L. Principal component analysis of chemical properties of soil saturation extracts from an irrigated Mediterranean area: Implications for calcite equilibrium in soil solutions. Geoderma 2009, 151, 407–416. [Google Scholar] [CrossRef]
  92. Bolan, N.S.; Loganathan, P.; Saggar, S. Calcium and Magnesium in Soils. In Encyclopedia of Soils in the Environment; Elsevier Ltd.: Amsterdam, The Netherlands, 2005; pp. 149–154. [Google Scholar]
  93. Ge, S.; Xu, H.; Ji, M.; Jiang, Y. Characteristics of Soil Organic Carbon, Total Nitrogen, and C/N Ratio in Chinese Apple Orchards. Open J. Soil Sci. 2013, 03, 213–217. [Google Scholar] [CrossRef]
  94. Paterson, E.; Goodman, B.A.; Farmer, V.C. The Chemistry of Aluminium, Iron and Manganese Oxides in Acid Soils. In Soil Acidity; Ulrich, B., Sumner, M.E., Eds.; Springer: Berlin/Heidelberg, Germany, 1991; pp. 97–124. [Google Scholar]
  95. Shamrikova, E.V.; Punegov, V.V.; Gruzdev, I.V.; Vanchikova, E.V.; Vetoshkina, A.A. Individual organic compounds in water extracts from podzolic soils of the Komi Republic. Eurasian Soil Sci. 2012, 45, 939–946. [Google Scholar] [CrossRef]
  96. Shamrikova, E.V.; Gruzdev, I.V.; Punegov, V.V.; Vanchikova, E.V. Influence of Biota on Low Molecular Weight Organic Acids in Soil Solutions of Taiga and Tundra Soils in the East-European Russia; Springer: Dordrecht, The Netherlands, 2013; pp. 107–111. [Google Scholar]
  97. Pierzynski, G.M.; McDowell, R.W.; Thomas Sims, J. Chemistry, Cycling, and Potential Movement of Inorganic Phosphorus in Soils. In Phosphorus: Agriculture and the Environment; Agronomy Monographs; American Society of Agronomy, Inc.: Madison, WI, USA, 2015; pp. 51–86. [Google Scholar]
Figure 1. Photographs of vegetation and profiles of the studied soils: (a) albic Glossic Retisol (Siltic, Cutanic) (RI); (b) albic Ferric Retisol (Siltic, Cutanic) II (RII).
Figure 1. Photographs of vegetation and profiles of the studied soils: (a) albic Glossic Retisol (Siltic, Cutanic) (RI); (b) albic Ferric Retisol (Siltic, Cutanic) II (RII).
Soilsystems 08 00063 g001
Figure 2. Molar distribution of protolytic functional groups (FG) for kaolinite with pKa steps (a) 1, (b) 0.5, and (c) 0.25; for (a), errors are shown, for (b,c), replicates are shown individually.
Figure 2. Molar distribution of protolytic functional groups (FG) for kaolinite with pKa steps (a) 1, (b) 0.5, and (c) 0.25; for (a), errors are shown, for (b,c), replicates are shown individually.
Soilsystems 08 00063 g002
Figure 3. Molar distribution of total protolytic functional groups (FG) in profiles of (a) Retisol I and (b) Retisol II. Standard deviation was calculated by 4 replicates.
Figure 3. Molar distribution of total protolytic functional groups (FG) in profiles of (a) Retisol I and (b) Retisol II. Standard deviation was calculated by 4 replicates.
Soilsystems 08 00063 g003aSoilsystems 08 00063 g003b
Figure 4. Molar distribution of protolytic functional groups (FG) of different strengths in mineral horizons of profiles of (a) Retisol I (RI) and (b) Retisol II (RII).
Figure 4. Molar distribution of protolytic functional groups (FG) of different strengths in mineral horizons of profiles of (a) Retisol I (RI) and (b) Retisol II (RII).
Soilsystems 08 00063 g004aSoilsystems 08 00063 g004b
Figure 5. Distribution depth profiles of elements in (a) RI and (b) RII; ox, oxalate-soluble.
Figure 5. Distribution depth profiles of elements in (a) RI and (b) RII; ox, oxalate-soluble.
Soilsystems 08 00063 g005
Figure 6. Principal component values: each point represents a horizon; red points are for Retisol I, and blue for Retisol II.
Figure 6. Principal component values: each point represents a horizon; red points are for Retisol I, and blue for Retisol II.
Soilsystems 08 00063 g006
Figure 7. Loading values; each point represents one variable; ex, exchangeable; ox, oxalate-soluble; and dit, dithionite-soluble. Each point corresponds to a variable, and its coordinates are loadings on PC1 and PC2.
Figure 7. Loading values; each point represents one variable; ex, exchangeable; ox, oxalate-soluble; and dit, dithionite-soluble. Each point corresponds to a variable, and its coordinates are loadings on PC1 and PC2.
Soilsystems 08 00063 g007
Table 1. Description of the research objects.
Table 1. Description of the research objects.
SoilsCoordinatesLocationVegetation of Ground CoverSoil Profile
R IN 61°45′43.6″
E 54°17′50.1″
h = 260
The intercrown spaceVaccinium myrtillus, Gymnocarpium dryopteris, Equisetum sylvaticum, Carex globularis.Oi–Oe–Oa–E/Eh–BEL–Btw–C
R IISpruce crown (260 years old)Hylocomium splendens, Equisetum sylvaticum, Vaccinium myrillus, Pleurozium schereberiOi–Oe–Oa–E/Eh–BELc, i–Btw–C
Table 2. Physicochemical properties of studied soils.
Table 2. Physicochemical properties of studied soils.
HorizonDepth, cmClay Fraction, %pHHydrolytic Acidity, mmol/100 gExchangeable Cations, mmol/kgCtot,
%
Ntot,
%
C:NOxalate Soluble, %
WaterSaltCaMgFe2O3Al2O3
Retisol I (intercrown space)
Oi0–34.63.958.821.47.844.21.6332
Oe3–64.33.356.315.72.9843.21.927
Oa6–84.02.973.711.43.0335.81.5826
E16–11113.92.914.90.600.282.200.148170.180.28
Eh11–14134.13.35.930.710.323.200.22171.270.54
E220–3084.93.88.830.570.2100.410.046100.660.44
30–4085.03.88.280.500.2430.210.03080.350.29
BEL*
BEL**
40–50105.33.86.691.170.770.170.02480.280.24
50–60135.33.86.532.441.680.120.02170.280.26
60–70245.53.85.488.33.400.120.02460.30.27
70–80245.63.84.929.44.30.110.02450.290.31
Btw180–90265.73.94.2310.45.60.140.03050.270.27
90–110285.84.03.7111.55.90.130.02170.280.26
Btw2110–125306.04.13.4012.06.10.140.02570.230.24
Cg130–140265.74.13.0512.16.20.130.02560.220.22
140–150286.14.22.7411.66.00.110.02550.210.21
Retisol II (under the crown of spruce)
Oi0–25.24.550.332.17.945.71.830
Oe2–44.84.060.225.05.146.21.928
Oa4–64.13.188.210.82.5041.81.729
E6–10134.03.114.90.630.322.00.125190.30.34
Eh,c10–20114.33.516.90.760.322.00.128181.280.54
BELc,i20–30174.53.616.90.730.312.10.138181.220.7
30–45224.83.815.50.990.500.910.074140.920.75
Btw145–60265.13.713.92.281.540.290.033100.520.6
Btw260–80275.33.810.16.53.400.290.03890.440.49
Table 3. The pKa values of real titration curves of model acid samples. The total concentrations of groups with adjacent pKa are in bold.
Table 3. The pKa values of real titration curves of model acid samples. The total concentrations of groups with adjacent pKa are in bold.
Acid MixturepKa
2.53.06.57.07.59.09.510.0
0.00625 M H3PO4
and 0.00625 M H3BO3
0.003880.0026900.004390.002000.000000.006190
0.006570.006400.00619
0.00375 M H3PO4
and 0.00875 M H3BO3
0.002480.001350.000330.002990.000000.003840.00449
0.003830.003320.00833
0.00875 M H3PO4
and 0.00375 M H3BO3
0.006340.002780.000000.006020.002860.000000.003140.00067
0.009120.008880.00381
Table 4. Main reactions in soils at different pHs according to the existing data.
Table 4. Main reactions in soils at different pHs according to the existing data.
pH RangeMain Determining Factors
4.0–5.0Reactions of carboxyl groups and H-acids of humic substances [78]
5.0–6.5Hydrolysis of exchangeable Al [79,80]
6.0–7.0HCO3, CO2 from air/water; deprotonation of carboxyl groups of water-soluble organic matter; amorphous aluminosilicates [81,82]
8.0–8.5Amino groups [83] and polymeric hydroxy compounds of Al [14]
8.5–9.5Weak OH groups of soluble phenolic, polyphenolic compounds and their hydrolysis products [83], phospholipids and diesters [84]
9.5–10.0Al and Fe complex compounds with organic ligands and phenolic compounds [85]
Table 5. Correlation between pKa spectra and other soil parameters.
Table 5. Correlation between pKa spectra and other soil parameters.
pKa or pKa RangeR I CorrelationR II CorrelationR I + R II
5.0Strong positive correlations with the amount of C, N, C/N, oxalate-soluble Al2O3 and Fe2O3, dithionite-soluble Fe2O3 and SO3Positive correlations with the amount of oxalate-soluble Al2O3 and dithionite-soluble Fe2O3.
Negative correlation with C, N, C/N
Weak correlations with the amount of C, N. Correlations with the amounts of oxalate-soluble Al2O3, Fe2O3
6.5Strong positive correlations with the amount of C, N, C/N, oxalate-soluble Al2O3 and Fe2O3, dithionite-soluble Fe2O3 and SO3Strong positive correlations with the amount of C and N; negative correlation with oxalate-soluble Al2O3 and Fe2O3 and dithionite-soluble Fe2O3Correlations with the amount of C, N. Correlations with the amounts of oxalate-soluble Fe2O3
8.5Positive correlations with the amount of C, N, C/N, oxalate-soluble Al2O3 and Fe2O3, dithionite-soluble Fe2O3 and SO3Positive correlations with the amounts of oxalate-soluble Al2O3, Fe2O3Correlations with the amounts of oxalate-soluble Al2O3
9.0Weak correlations with the amount of C, N, C/N, oxalate-soluble Al2O3 and Fe2O3, dithionite-soluble Fe2O3 No correlations with the amount of C, N; correlation with the amount of CaO and dithionite-soluble Fe2O3
Positive correlations with C, N, C/N, Al2O3 and Fe2O3, P2O5 and hydrolytic acidity
No correlations
10.0Negative correlations oxalate-soluble Al2O3 and Fe2O3 and dithionite-soluble Fe2O3Correlations with the amounts of oxalate-soluble Al2O3, P2O5
4.0–6.0Correlations with the amount of C, N, oxalate-soluble Al2O3 and Fe2O3Negative correlations with CaO, weak correlation with oxalate-soluble Al2O3, P2O5Correlations with the amount of C, N. Correlations with the amounts of oxalate-soluble Al2O3, Fe2O3, P2O5
6.5–8.5Correlations with the amount of C, N, oxalate-soluble Al2O3 and Fe2O3 and SO3Negative correlations with exchangeable CaO and MgO, positive correlation with oxalate-soluble Al2O3 and Fe2O3 and dithionite-soluble Fe2O3, P2O5Correlations with the amount of C, N. Correlations with the amounts of oxalate-soluble Al2O3, Fe2O3, P2O5
4.0–5.0Correlations with the amount of C, N, oxalate-soluble Al2O3 and Fe2O3 and SO3Positive correlation with the hydrolytic acidity and with oxalate-soluble Al2O3Weak correlations with the amount of C, N. Correlations with the amounts of oxalate-soluble Al2O3
5.5–6.5Weak correlations with the amount of C, N, oxalate-soluble Al2O3 and Fe2O3 and SO3Weak positive correlation with C, N, negative correlation with oxalate-soluble Al2O3 and Fe2O3 and dithionite-soluble Fe2O3No correlations
7.5–8.5Correlations with the amount of C, N, oxalate-soluble Al2O3 and Fe2O3 and SO3Correlation with oxalate-soluble Al2O3, P2O5Weak correlations with the amount of C, N. Correlations with the amounts of oxalate-soluble Al2O3, Fe2O3, P2O5
9.0–10.0Correlations with the amount of C, N, oxalate-soluble Al2O3 and Fe2O3 and SO3Positive correlation with C, N, C/N, P2O5, Al2O3 and Fe2O3 and hydrolytic acidityCorrelations with the amount of C, N, oxalate-soluble Al2O3 and Fe2O3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Matveeva, N.V.; Garmash, A.V.; Shishkin, M.A.; Dymov, A.A.; Rogova, O.B.; Volkov, D.S.; Proskurnin, M.A. Fast High-Resolution pKa Spectrotitrimetry for Quantification of Surface Functional Groups of Retisols. Soil Syst. 2024, 8, 63. https://doi.org/10.3390/soilsystems8020063

AMA Style

Matveeva NV, Garmash AV, Shishkin MA, Dymov AA, Rogova OB, Volkov DS, Proskurnin MA. Fast High-Resolution pKa Spectrotitrimetry for Quantification of Surface Functional Groups of Retisols. Soil Systems. 2024; 8(2):63. https://doi.org/10.3390/soilsystems8020063

Chicago/Turabian Style

Matveeva, Natal’ya V., Andrei V. Garmash, Mikhail A. Shishkin, Alexey A. Dymov, Olga B. Rogova, Dmitry S. Volkov, and Mikhail A. Proskurnin. 2024. "Fast High-Resolution pKa Spectrotitrimetry for Quantification of Surface Functional Groups of Retisols" Soil Systems 8, no. 2: 63. https://doi.org/10.3390/soilsystems8020063

APA Style

Matveeva, N. V., Garmash, A. V., Shishkin, M. A., Dymov, A. A., Rogova, O. B., Volkov, D. S., & Proskurnin, M. A. (2024). Fast High-Resolution pKa Spectrotitrimetry for Quantification of Surface Functional Groups of Retisols. Soil Systems, 8(2), 63. https://doi.org/10.3390/soilsystems8020063

Article Metrics

Back to TopTop